Física II Minas. Contenido Premium.

Material para alumnos matriculados

  • fRESNEL

  • COEFRESNELd100

  • COEFRESNEL

  • FIGURAS

  • CARACTER

FRESNEL

% REFRACCIÓN Y REFLEXIÓN DEL CAMPO EM
% CALLS: functions 'COEFRESNEL', 'FIGURAS', 'CARACTER'
clear
close all
clc
c=3*1e8; e0=1/(36*pi*1e9);
disp('DATOS DE ENTRADA:')
n1=input('Índice de refracción del medio 1, por favor: ');
n2=input('Índice de refracción del medio 2, por favor: ');
Qi=input('Ángulo de incidencia (º), por favor (si es Brewster, introduce 100): ');
if Qi==100
Qi=atan(n2/n1);
Qt=pi/2-Qi;
else
Qi=Qi*pi/180;
Qt=asin((n1/n2)*sin(Qi));
end
disp('--------------------------------------------------')
disp('INTRODUCE LA POLARIZACIÓN DE LA LUZ INCIDENTE: ')
disp('Lineal: 1')
disp('Circular: 2')
disp('Elíptica: 3')
POL=input('Elige 1,2 o 3: ');
disp('--------------------------------------------------')

if POL==1 % LINEAL
disp('Elige tipo de datos de entrada:')
disp('Ángulo PHI e intensidad Ii: 1')
disp('Ángulo PHI y amplitud E0i: 2')
disp('Componentes E0i-paralelo y E0i-perpendicular: 3')
disp('--------------------------------------------------')
TIPO=input('Elige 1,2 o 3: ');
if TIPO==1
PHI=input('Ángulo de E0i con el plano de incidencia (º), por favor: ');
Ii=input('Intensidad incidente (W/m^2), plis :');
E0i=sqrt((2*Ii)/(c*e0*n1));
Ei=[E0i*cosd(PHI);E0i*sind(PHI)];
elseif TIPO==2
PHI=input('Ángulo de E0i con el plano de incidencia (º), por favor: ');
E0i=input('Amplitud incidente (V/m), por favor: ');
Ei=[E0i*cosd(PHI);E0i*sind(PHI)];
elseif TIPO==3
E01=input('Componente parlela del campo (V/m), por favor: ');
E02=input('Componente perpendicular del campo (V/m), por favor: ');
Ei=[E01;E02];
else
disp('ERROR ***************************************')
end
elseif POL==2 % CIRCULAR
disp('Elige tipo de datos de entrada:')
disp('Amplitud E0i y helicidad (dextógira o levógira): 1')
disp('Intensidad Ii y helicidad (dextógira o levógira): 2')
disp('--------------------------------------------------')
TIPO=input('Elige 1 o 2: ');
if TIPO==1
E0i=input('Valor de la amplitud del campo (V/m), por favor: ');
elseif TIPO==2
Ii=input('Valor de la intensidad incidente (W/m^2), por favor: ');
E0i=sqrt(Ii/(c*e0*n1));
else
disp('ERROR ***************************************')
end
disp(' ')
HEL=input('Luz dextrógira (1) o levógira (2): ');
if HEL==1
delta=pi/2;
elseif HEL==2
delta=-pi/2;
else
disp('ERROR ***************************************')
end
Ei=[E0i*exp(1i*delta);E0i];
elseif POL==3 % ELÍPTICA
disp('Elige tipo de datos de entrada:')
disp('E0i-paralelo, E0i-perpendicular: 1')
disp('E0i-paralelo, factor x (E0i-perp=x*E0i-paral): 2')
disp('Intensidad Ii y factor x (E0i-perp=x*E0i-paral): 3')
disp('--------------------------------------------------')
TIPO=input('Elige 1,2 o 3: ');
if TIPO==1
E01=input('Componente E0i-paralela (V/m), por favor: ');
E02=input('Componente E0i-perpendicular (V/m), por favor: ');
elseif TIPO==2
E01=input('Componente E0i-paralela (V/m), por favor: ');
x=input('Factor x, por favor: ');
E02=x*E01;
elseif TIPO==3
Ii=input('Intensidad incidente (W/m^2), por favor: ');
x=input('Factor x, por favor: ');
E01=sqrt(2*Ii/(c*e0*(1+x^2)));
E02=x*E01;
else
disp('ERROR ***************************************')
end
disp(' ')
HEL=input('Luz dextrógira (1) o levógira (2): ');
if HEL==1
delta=pi/2;
elseif HEL==2
delta=-pi/2;
else
disp('ERROR ***************************************')
end
Ei=[E01*exp(1i*delta);E02];
else
disp('ERROR ***************************************')
end

[Rpar,Rper,Tpar,Tper]=COEFRESNEL(n1,n2,Qi);
if abs(Rpar)<1e-10
Rpar=0;
end
R=[Rpar,0;0,Rper];
T=[Tpar,0;0,Tper];
Er=R*Ei;
Et=T*Ei;

caracterI=CARACTER(Ei);
caracterR=CARACTER(Er);
caracterT=CARACTER(Et);

Ii=(1/2)*c*e0*n1*(Ei'*Ei);
Ir=(1/2)*c*e0*n1*(Er'*Er);
It=(1/2)*c*e0*n2*(Et'*Et);

disp(' ')
disp('##################################################')
disp('##################################################')
disp(' ')

disp(['Ángulo de incidencia: ',num2str(Qi*180/pi),'º'])
disp(['Ángulo de refracción: ',num2str(Qt*180/pi),'º'])
disp(['Ángulo de Brewster: ',num2str(atand(n2/n1)),'º'])
if n1>n2
QL=asind(n2/n1);
disp(['Ángulo límite: ',num2str(QL),'º'])
if Qi>pi*QL/180
disp('**************************************************')
disp('¡ATENCIÓN! Incidencia por encima del ángulo límite')
disp('**************************************************')
end
end
disp(' ')
disp('Coeficientes de Fresnel:')
disp(['R_paralelo: ',num2str(Rpar)])
disp(['R_perpendicular: ',num2str(Rper)])
disp(['T_paralelo: ',num2str(Tpar)])
disp(['T_perpendicular: ',num2str(Tper)])
disp(' ')
disp('Luz INCIDENTE --------------------------------------')
disp(['Amplitud Ei_paralela: ',num2str(abs(Ei(1))),' V/m'])
disp(['Amplitud Ei_perpendicular: ',num2str(abs(Ei(2))),' V/m'])
if caracterI(1)==char('L')
disp(['Ángulo PHIi: ',num2str(atand(abs(Ei(2))/abs(Ei(1)))),'º'])
end
disp(['Intensidad Ii (W/m^2): ',num2str(Ii),' W/m^2'])
disp(['Helicidad: ',caracterI])
disp(' ')
disp('Luz REFLEJADA --------------------------------------')
disp(['Amplitud Er_paralela: ',num2str(abs(Er(1))),' V/m'])
disp(['Amplitud Er_perpendicular: ',num2str(abs(Er(2))),' V/m'])
if caracterR(1)==char('L')
disp(['Ángulo PHIr: ',num2str(atand(abs(Er(2))/abs(Er(1)))),'º'])
end
disp(['Intensidad Ir (W/m^2): ',num2str(Ir),' W/m^2'])
disp(['Helicidad: ',caracterR])
disp(' ')
disp('Luz REFRACTADA -------------------------------------')
disp(['Amplitud Et_paralela: ',num2str(abs(Et(1))),' V/m'])
disp(['Amplitud Et_perpendicular: ',num2str(abs(Et(2))),' V/m'])
if caracterT(1)==char('L')
disp(['Ángulo PHIt: ',num2str(atand(abs(Et(2))/abs(Et(1)))),'º'])
end
disp(['Intensidad It (W/m^2): ',num2str(It),' W/m^2'])
disp(['Helicidad: ',caracterT])
disp(' ')
disp(['Comprobación intensidad: ',num2str(Ii-Ir-It*cos(Qt)/cos(Qi))])
disp('(Correcto si es próximo a 0)')

FIGURAS(Ei,Er,Et)

¡Hola! Haga clic a continuación y nos pondremos en contacto con usted lo antes posible.

Si quieres podemos ayudarte