Reflactancia


UNIVERSIDAD SANTO TOMAS
DIVISIÓN DE INGENIERÍAS
FACULTAD DE INGENIERÍA DE TELECOMUNICACIONES
OPTOELECTRÓNICA



Julian Morales Ortega.
Para el siguiente análisis tendremos en cuenta que un n1 elcual será Es el índice de refracción del medio de origen n2 que será el índice de refracción en la fibra y n3 que será el índice de refracción del medio de destino.
n1=1;
n2=1.43;
n3=1.33;
Se obtienenpor medio de Matlab las siguientes gráficas:

En los casos de las gráficas analizando n2 y n3 se observan algunas irregularidades puesto que se generan algunos números imaginarios al momento dehallar el coeficiente de reflexión, los cuales generarían valores negativos al hallar la reactancia, por lo cual para evitar esta irregularidad se manejan ángulos hasta un valor de 68.75 y se obtienen lossiguientes resultados:
Código utilizado:
n1=1;
n2=1.43;
n3=1.33;
ang=0:0.01:pi/2;
%polarizacion paralela con n1 y n2r=((n1.*cos(ang)-((n2^2)-(n1^2).*(sin(ang)).^2).^0.5)./(n1.*cos(ang)+((n2^2)-(n1^2).*(sin(ang)).^2).^0.5)).^2;
subplot(2,2,1)
plot(ang,r,'c')
axis([0, pi/2, 0, 1])
title('polarizacion paralela con n1 y n2')
xlabel('Angulo(Rad)')
ylabel('Refractancia')
grid on

%polarizacionparalela con n2 y n3
r1=((n2.*cos(ang)-((n3^2)-(n2^2).*(sin(ang)).^2).^0.5)./(n2.*cos(ang)+((n3^2)-(n2^2).*(sin(ang)).^2).^0.5)).^2;
subplot(2,2,2)
plot(ang,r1,'b')
title('polarizacion paralela con n2 yn3')
xlabel('Angulo(Rad)')
ylabel('Refractancia')
grid on

%polarizacion perpendicular con n1 y n2r2=((-((n2)^2).*cos(ang)+n1.*((n2^2)-(n1^2).*(sin(ang)).^2).^0.5)./((n2^2).*cos(ang)+n1.*((n2^2)-(n1^2).*(sin(ang)).^2).^0.5)).^2;
subplot(2,2,3)
plot(ang,r2,'r')
title('polarizacion perpendicular con n1 y n2')
xlabel('Angulo(Rad)')
ylabel('Refractancia')
axis([0, pi/2, 0, 1])
grid on%polarizacion perpendicular con n2 y n3
ro=(-((n3)^2).*cos(ang)+n2.*((n3^2)-(n2^2).*(sin(ang)).^2).^0.5)./((n3^2).*cos(ang)+n2.*((n3^2)-(n2^2).*(sin(ang)).^2).^0.5);
r3=ro.^2
subplot(2,2,4)
plot(ang,r3,'g')...