Gráficas de reflectancias de Fresnel y su contraste

La linea negra es el contraste entre señales polarizadas. La linea azul
es la R_perpendicular, y la linea roja R_paralela.

Las ecuaciones de Fresnel son fáciles de graficar y presentar en un reporte. Este es un ejemplo en un formato compacto para las ecuaciones de Fresnel para la intensidad de reflexión:




Para hacer las ecuaciones use Latex, en el sitio codecogs. El cual recomiendo mucho.

Este es el código Matlab de estas ecuaciones para hacer las gráficas de este post

%%%%%%%%%%%%%%%%
clc; clear; close all % limpiando memoria y pantallas

n_inc = 1.0; % índice de refracción incidente
n_trans = 1.5; % índice de refracción transmitido

theta_inc = 0:0.001:90; % ángulo en radianes

rad_inc = theta_inc*pi/180; %conversion a radianes
rad_trans = asin((n_inc/n_trans)*sin(rad_inc)); % el ángulo de transmisión

% Reflexion con polarización perpendicular
R_s = ((sin(rad_trans - rad_inc) )./( sin(rad_trans + rad_inc))).^2;

% Reflexion con polarización paralela
R_p = ((tan(rad_trans - rad_inc) )./(tan(rad_trans + rad_inc))).^2;

[valor_R_p, indice_R_p] = min(R_p);

% Contraste entre R_s y  R_p
Con = (R_s-R_p)./(R_s+R_p);
[valor_Con, indice_Con] = max(Con);

% Graficando
hold on
plot(theta_inc, Con, 'k', 'LineWidth', 6);
plot(theta_inc, R_s, 'b', 'LineWidth', 4);
plot(theta_inc, R_p, 'r', 'LineWidth', 2);
plot(theta_inc(indice_R_p).*Con./Con, theta_inc , '--k', 'LineWidth', 1);
%plot(theta_inc(indice_Con).*Con./Con, theta_inc , '--y', 'LineWidth', 2);
axis([0 90.2 0 1.02])
legend1 = legend('Contraste', 'R_s',  'R_p');
set(legend1,'Location','NorthWest');

hold off
% fin de la rutina