| Date 08/09/2010 |
|
||||||
|
| Script S3_2_5.m | |||
|
|||
|
%==============================================
%A calcite plate %============================================== % %refractive indices of calcite %for lambda = 0.5893 micron no=1.6584; ns=1.4864; % %the angle of incidence in degrees alfa=45; arg=sin(alfa*pi/180); argo=arg/no; args=arg/ns; %angles of refraction alfa2o=(180/pi)*asin(argo) alfa2e=(180/pi)*asin(args) %distance t between the points where %the ordinary and extraordinary ray %are incident on the lower surface of the plate t=10*(tan(alfa2e*pi/180)-tan(alfa2o*pi/180)) %distance s between the two rays emerging %from the lower surface of the plate s=t*arg % %relative intensity for ordinary (blue line)and extraordinary (red) ray %varying the angle tetag between the optic axis of the plate %and transmission axis of the polarizer from 0° to 90° teta=linspace(eps,pi/2,20) Iro=sin(teta).*sin(teta) Ire=cos(teta).*cos(teta) sum=Ire+Iro tetag=teta*180/pi plot(tetag,Ire,'ro-',tetag,Iro,'bd-',tetag,sum,'g*-'),grid on title('Iro (blue) and Ire (red) varying tetag from 0° to 90°') axis([0,90,0,1.1]) %============================================== % %============================================== %============================================== %Corrigenda %in Fig. 3.19a the angle of incidence is not %alfa1 but alfa %On the third row (pag. 132) change Whel with When %============================================== %============================================== |
|||
| Top | |||
