| Date 06/09/2010 |
|
||||||
|
| Script S3_2_15.m | |||
|
|||
|
%==============================================
%A quarter-wave plate 3 %============================================== % %a polarized light whose direction of oscillation %forms an angle of 0° with the x axis tetg=0 tet=tetg*pi/180; cc1=cos(tet); ss1=sin(tet); v1=[cc1;ss1] %a quarter-wave plate whose optic axis subtends %an angle of 30° with x axis %its matrix fi=30*pi/180 SS=sin(fi)^2 CC=cos(fi)^2 b11=CC-i*SS; b12=(sin(fi)*cos(fi))*(1+i) b21=b12; b22=-i*CC+SS; m1=[b11 b12;b21 b22] %polarization form of the beam %emerging from the quarter-wave vs=m1*v1 % %================================================== %transform of the complex numbers of the vector vs %to polar from rectangular form %first number a=vs(1) aM=abs(a) ateta=angle(a) atetag=ateta*180/pi %second number b=vs(2) bM=abs(b) bteta=angle(b) btetag=bteta*180/pi %ellipticity eta=bM/aM %phase shift fidiff=ateta-bteta fidiffg=fidiff*180/pi %================================================= % %intensity of the beam emerging from quarter wave vst=vs' Is=vst*vs %matrix of the rotating polarizer P2 %whose transmissionn axis varies 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]; %beam emerging from polarizer P2 vf=m2*m1*v1; %the complex conjugate vft=vf'; %intensity Int(i)=vft*vf; end Int %plot intensity emerging from polarizer P2 varying tetag from 0° to 360° plot(tetag,Int,'ro-'),grid on title('Intensity of the beam emerging from polarizer P2 varying tetag') axis([0,360,0,1]) %============================================== % |
|||
| Top | |||
