| Date 06/09/2010 |
|
||||||
|
| Script S5_2_10.m | |||
|
|||
|
%==============================================
%Young 9 %============================================== % %the distance L (in mm) L=1000; %distance d, width h and wavelength lambda in microns d=50; h=10; lambda=0.5; k=2*pi/lambda; %thickness of the films in microns s=10; %refractive indexes of the films n1=1.78; n2=1.42; %orders between -8 and 23 m=-8:1:23 %============================================== %vertical positions of the fringes %============================================== %see the problem Sec.5.2.10 for definitions %of the following formulae % %first term of arcsin diff=s*(n1-n2); %second term of arcsin prod=lambda*m; %argument of arcsin arg=(diff-prod)/d; %corresponding angle in radians ang=asin(arg); %and in degrees ang_g=ang*180/pi %vertical positions of the fringes y=L*tan(ang) % %============================================== %relative intensity due to interference %============================================== % %first term gamma1=diff; %second term gamma2=d*sin(ang); %gamma gamma=k*(gamma1-gamma2); Ir1=(cos(gamma/2)).*(cos(gamma/2)); % %============================================== %relative intensity due to diffraction %============================================== % %alfa alfa=(h*pi/lambda)*sin(ang); num=sin(alfa).^2; den=alfa.^2; Ir2=num./den; % %============================================== %Effective relative intensity %============================================== % Ir=Ir1.*Ir2 plot(m,Ir,'ro-'),grid on,title('Relative intensity varying m') axis([-9 24 0 1]) figure plot(ang_g,Ir,'bd-'),grid on axis([-8 9 0 1]) title('Relative intensity varying ang-g') %============================================== % |
|||
| Top | |||
