| Date 06/09/2010 |
|
||||||
|
| Script S3_2_9A.m | |||
|
|||
|
%==============================================
%Two calcite prisms 2A %============================================== % %refractive indices of calcite %for lambda=0.5893 micron no=1.6584; ne=1.4864; %angle of incidence on the second prism alfa1g=55; alfa1=alfa1g*pi/180; %corresponding angles of refraction %on the surface AC %for the ordinary and extraordinary rays alfa2o=asin((ne/no)*sin(alfa1)); alfa2e=asin((no/ne)*sin(alfa1)); alfa2og=alfa2o*180/pi alfa2eg=alfa2e*180/pi %angles of incidence on the surface CD %for the ordinary and extraordinary rays alfap1o=abs(alfa1-alfa2o); alfap1e=abs(alfa1-alfa2e); alfap1og=alfap1o*180/pi alfap1eg=alfap1e*180/pi %angles of refraction on the surface DC %for the ordinary and extraordinary rays alfap2o=asin(no*sin(alfap1o)); alfap2e=asin(ne*sin(alfap1e)); alfap2og=alfap2o*180/pi alfap2eg=alfap2e*180/pi %the critical angle for the ordinary ray in the first prism alfa1c=asin(ne/no); alfa1cg= alfa1c*180/pi %a minimun angle (one tenth of a degree) alfamin=0.1*pi/180; %a maximum angle obtained subtracting alfamin from the critical angle alfamax=alfa1c-alfamin; %intervall (alfamin,alfamax) is divided in 20 parts alfa1v=linspace(alfamin,alfamax,20); alfa1vg=alfa1v*180/pi; %previous steps are repeated for every angle alfa1v for i=1:20 %angle on incidence on the surface AC alfa1s=alfa1v(i); %angle of the refracted ray on surface AC alfa2e(i)=asin((no/ne)*sin(alfa1s)); alfa2es=alfa2e(i); %angle of incidence on the surface DC alfap1e(i)=abs(alfa1s-alfa2es); alfap1es=alfap1e(i); %angle of refraction on the surface DC alfap2e(i)=asin(ne*sin(alfap1es)); alfap2es=alfap2e(i); alfap2eg(i)=alfap2es*180/pi; end %the angles of refraction on the surface DC varying alfa1vg alfap2eg plot(alfa1vg,alfap2eg,'ro-'),grid on title('Angle emerging from DC varying angle of incidence on AC') %============================================== % |
|||
| Top | |||
| Script S3_2_9B.m | |||
|
|||
|
%==============================================
%Two quartz prisms 2 %============================================== % %refractive indices of quartz for lambda = 0.5893 micron no=1.5443; ne=1.5534; %angle of incidence on the surface AC alfa1g=55; alfa1=alfa1g*pi/180; %angles of refraction on the surface AC alfa2o=asin((ne/no)*sin(alfa1)); alfa2e=asin((no/ne)*sin(alfa1)); alfa2og=alfa2o*180/pi; alfa2eg=alfa2e*180/pi; %angles of incidence on the surface DC alfap1o=abs(alfa1-alfa2o); alfap1e=abs(alfa1-alfa2e); alfap1og=alfap1o*180/pi; alfap1eg=alfap1e*180/pi; %angles of refraction on the surface DC alfap2o=asin(no*sin(alfap1o)); alfap2e=asin(ne*sin(alfap1e)); alfap2og=alfap2o*180/pi alfap2eg=alfap2e*180/pi %============================================== % %============================================== %============================================== %Corrigenda %the angles alfap2og and alfap2eg %are a bit different from those on the first %row of pag. 140 %change 0.74° in 7.492° and 0.70° in 0.740° %The right results are obtained running this %script %============================================== %============================================== |
|||
| Top | |||
