| Date 06/09/2010 |
|
||||||
|
| Script S3_2_14.m | |||
|
|||
|
%==============================================
%A quarter-wave plate2 %============================================== % %a beam of polarized light whose direction %of oscillation forms an angle 30° with x axis % %its vector tetg=30 tet=tetg*pi/180; cc1=cos(tet) ss1=sin(tet) v1=[cc1;ss1] %axes of ellipse: "a" parallel to the x axis %and "b" parallel to the y axis a=cc1^2 b=ss1^2 %matrix of the quarter-wave whose optic axis %is parallel to the x axis m1=[1 0;0 -i] %polarizzation form of the beam %emerging from the quarter-wave vs=m1*v1 vst=vs' Is=vst*vs %matrix of the rotating polarizer P2 %whose transmission axis changes %its angle with the x axis %from 0° to 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]; %vector of the beam emerging from P2 vf=m2*m1*v1; %the complex conjugate vft=vf'; %intensity(i) when angle is teta(i) Int(i)=vft*vf; end %array of the intensities Int %plot of intensity emerging from P2 function of angle of rotation tetag plot(tetag,Int,'ro-'),grid on title('Intensity of the beam emerging from polarizer P2 function of tetag') axis([0,360,0,1]) %============================================== % |
|||
| Top | |||
