POLDEMO Demonstrates polarization examples

Demonstration of an imperfect polarizer
by Chuck DiMarzio
   Northeastern University
   October 2008
!! This file may be copied, used, or modified for educational and
!! research purposes provided that this header information is not
!! removed or altered, and provided that the book is cited in
!! publications, as DiMarzio, Charles A., Optics for Engineers,
!! CRC Press, Boca Raton, FL, 2011.
!! http://www.crcpress.com
!! Other distribution is prohibited without permission.

(See Chapter 6 for variable names)

echo on
inputpower=5e-3; % Watts
insertionloss=0.08;
extinctionratio=10000;
taux=sqrt(1-insertionloss) % tau values are for fields
tauy=taux/sqrt(extinctionratio)
thetad=(0:.1:90);  % degrees
theta=thetad*pi/180;
transmission=(taux*cos(theta)).^2+(tauy*sin(theta)).^2; % Power
angleout=atan2(tauy*sin(theta),taux*cos(theta));
outputpower=inputpower*transmission;
malus=inputpower*(cos(theta)).^2;  %  Test against Malus Law
echo off
fig1=figure;plot(thetad,outputpower*1e3,'-',thetad,malus*1e3,'--');
grid on;
xlabel('\theta_{in}, Polarization Angle, Deg');
ylabel('P_{out}, Output Power, mW');
legend('Output','Malus Law','Location','SouthWest');

test=(thetad>87); % Don't plot too close to 90 degrees
fig2=figure;
plot(thetad(test),(outputpower(test)*1e6),'-',...
     thetad(test),(malus(test)*1e6),'--');
xlabel('\theta_{in}, Polarization Angle, Deg');
ylabel('P_{out}, Output Power, \mu W');
legend('Output','Malus Law','Location','SouthWest');


fig3=figure;plot(thetad,angleout*180/pi);
xlabel('\theta_{in}, Polarization Angle, Deg');
ylabel('\theta_{out}, Polarization Angle, Deg');
inputpower=5e-3; % Watts
insertionloss=0.08;
extinctionratio=10000;
taux=sqrt(1-insertionloss) % tau values are for fields 
taux =
    0.9592
tauy=taux/sqrt(extinctionratio)
tauy =
    0.0096
thetad=(0:.1:90);  % degrees
theta=thetad*pi/180;
transmission=(taux*cos(theta)).^2+(tauy*sin(theta)).^2; % Power
angleout=atan2(tauy*sin(theta),taux*cos(theta));
outputpower=inputpower*transmission;
malus=inputpower*(cos(theta)).^2;  %  Test against Malus Law
echo off