| Date 08/09/2010 |
|
||||||
|
| Script S5_2_13.m | |||
|
|||
|
%==============================================
%Grating 2 %============================================== % %wavelength lambda1 and lambda2 in micron lambda1=0.5890; lambda2=0.5896; %their difference Dlambda=lambda2-lambda1; %angular difference in sixtieths of a degree minutes=4.4; %and in radians delta=(minutes/60)*(pi/180); %============================================== % %Calculus of d, the distance between the slits %see Problem (Sec.5.2.13) for formulae lambda1q=lambda1^2; lambda2q=lambda2^2; lambdap=2*lambda1*lambda2; num=lambda1q+lambda2q-lambdap*cos(delta); den=sin(delta)^2; A=num/den Aq=sqrt(A); d=4*sqrt(A) % %============================================== % %d is retrieved using approximations %sin(alfa)=alfa and cos(alfa)=1 if alfa is very small %see Problem (Sec.5.2.13) for formulae % %preliminary calculi for teta1 num=lambda1*delta; den=Dlambda; arg=num/den; %teta1 in radians teta1=atan(arg); %teta1 in degrees teta1g=teta1*180/pi %teta2 in radians teta2=teta1+delta; % and in degrees teta2g=teta2*180/pi %check of their difference in degrees difftetag=teta2g-teta1g %check of their difference in minutes difftetap=(teta2g-teta1g)*60 %d, alias d1, using teta1 d1=4*lambda1/sin(teta1) %d, alias d2, using teta2 d2=4*lambda1/sin(teta2) % %============================================== %N, alias N1 and N2, is retrieved from formula %of resolving power % N1=lambda1/(4*Dlambda) %or from the following formula N2=lambda1/(delta*d1*cos(teta1)) % %============================================== %a check assuming known d d=3; m=4; %true teta1g, alias alfa1g alfa1g=(asin(m*lambda1/d))*180/pi %true teta2g, alias alfa2g alfa2g=(asin(m*lambda2/d))*180/pi %their difference in degrees Dalfa=(alfa2g-alfa1g) %and in minutes Dalfap=Dalfa*60 %errors in minutes Er1g=(abs(alfa1g-teta1g))*60 Er2g=abs(alfa2g-teta2g)*60 %============================================== % |
|||
| Top | |||
