| Date 06/09/2010 |
|
||||||
|
| Script S5_2_9.m | |||
|
|||
|
%==============================================
%Young 8 %============================================== % %angles between -pi/6 and pi/6 teta=linspace(-pi/6,pi/6,200); %corresponding values in degrees tetag=teta*180/pi %wavelength in microns lambda=0.7; %refractive index of the film n=1.38; %the thickness of the film s=lambda/(2*(n-1)) %width of the slits h=3*lambda; %distance between them d=5*lambda; %preliminary calculi cost=(pi/lambda); beta=cost*d*sin(teta); alfa=cost*h*sin(teta); alfa2=alfa.*alfa; coef0=sin(alfa).*sin(alfa); %part of relative intensity due to diffraction Ir1=coef0./alfa2; %part of relative intensity due to interference Ir2=sin(beta).*sin(beta); %the effective relative intensity Ir=Ir1.*Ir2 %plot of the three relative intensities subplot(3,1,1) plot(tetag,Ir1,'ro-'),grid on,axis([-20 20 0 1]) title('Relative intensity due to diffraction') subplot(3,1,2) plot(tetag,Ir2,'bd-'),grid on,axis([-20 20 0 1]) title('Relative intensity due to interference') subplot(3,1,3) plot(tetag,Ir,'g*-'),grid on,axis([-20 20 0 1]) title('The effective relative intensity') % %the second of the two maxima %its value val_Ir and its position %ind in the array Ir [val_Ir ind]=max(Ir) %corresponding angle of the second maximum ang=tetag(ind) %============================================== % %============================================== %============================================== %Corrigenda %change the value "5.2°" for the angular positions %for the maxima of Ir with "4.97°" %This script gives the value "4.9749°" %============================================== %============================================== |
|||
| Top | |||
