| Date 08/09/2010 |
|
||||||
|
| Script S3_2_12.m | |||
|
|||
|
%==============================================
%Half-wave plate %============================================== % %natural light incident on the first polarizer P1 %whose transmission axis subtends an angle of 45° %with x axis v1=(1/sqrt(2))*[1;1] %matrix of the half-wave retarder whose optic %axis is parallel to the x axis m1=[1 0;0 -1] %intensity of the beam emerging from retarder vs=m1*v1 vst=vs'; Is=vst*vs %matrix for P2 whose transmission axis rotates %around the z axis. 50 angular positions %between 0° and 360° are considered max=50 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]; %form of polarization of %the beam emerging from POL2 vf=m2*m1*v1; %the complex conjugate vft=vf'; %intensity Int(i)=vft*vf; end %array of the intensities Int %plot of intensity emerging from half-wave plate varying tetag plot(tetag,Int,'ro-'),grid on title('Intensity emerging from half-wave plate varying tetag from 0° to 360°') axis([0,360,0,1]) %============================================== % |
|||
| Top | |||
