| Date 08/09/2010 |
|
||||||
|
| Script S3_2_11.m | |||
|
|||
|
%==============================================
%Three polarizers %============================================== % %a polarized beam with E oscillating in the plane xz v1=[1;0] %first polarizer with transmission axis parallel to the x axis m1=[1 0;0 0] %third polarizer with transmission axis parallel to the y axis m3=[0 0;0 1] %the second rotating polarizer %max angular positions are defined between 0° and 360° max=100; tetag=linspace(0,360,max) teta=tetag*pi/180; for i=1:max %elements of the matrix c1=cos(teta(i)); s1=sin(teta(i)); a11=c1^2; a12=c1*s1; a21=c1*s1; a22=s1^2; m2=[a11 a12;a21 a22]; %beam emerging from the second polarizer vf21=m2*m1*v1; %the complex conjugate transpose vf21t=vf21'; %intensity after the second polarizer Int21(i)=vf21t*vf21; %beam emerging fron the third polarizer vf321=m3*m2*m1*v1; %the complex conjugate transpose vf321t=vf321'; %final intensity Int321(i)=vf321t*vf321; end %array of intensities after the second polarizer Int21 %array of intensities after the third polarizer Int321 plot(tetag,Int21,'ro-'),grid on title('Intensity of the beam emerging from the second polarizer') axis([0,360,0,1]) figure plot(tetag,Int321,'bd-'),grid on title('Final intensity emerging from the third polarizer') axis([0,360,0,0.25]) figure %................................................... %final intensity using directly the Malus'law %................................................... IMalus=0.25*sin(2*teta).^2 plot(tetag,IMalus,'g*-'),grid on title('Final intensity using directly the Malus law') axis([0,360,0,0.25]) %============================================== % |
|||
| Top | |||
