Home

Introduction

Contents

Scripts

FeedBack

 Date 08/09/2010
Script S4_2_6A.m
Download Script S4_2_6A.m
%==============================================
%Michelson 2A
%a minor number m of fringes (1001)
%and distance d (about 295 microns)
%are considered
%==============================================
%
%wavelengths for a sodium-vapor lamp
lambda1=0.58900;
lambda2=0.58958;
%preliminary calculi
lambda=(lambda1+lambda2)/2;
k=2*pi/lambda;
deltak=(2*pi)*((lambda2-lambda1)/(lambda1*lambda2));
%
%condition for the first minimum of visibility d_star*deltak=pi/2
d_star=pi/(2*deltak)
fringes=d_star*2/lambda
%
m=0:1000;
max=length(m);
for i=1:max
%maximum of intensity
dmax=m(i)*lambda/2;
arg1max=deltak*dmax;
Amax=cos(arg1max);
arg2max=2*k*dmax;
Bmax=cos(arg2max);
Imax(i)=4*(1+Amax*Bmax);
%minimum of intensity
dmin=(2*m(i)+1)*lambda/4;
arg1min=deltak*dmin;
Amin=cos(arg1min);
arg2min=2*k*dmin;
Bmin=cos(arg2min);
Imin(i)=4*(1+Amin*Bmin);
%visibility
V(i)=(Imax(i)-Imin(i))/(Imax(i)+Imin(i));
end
dmax
dmin
Imax;
Imin;
V;
plot(m,Imax,'ro-',m,Imin,'bd-'),grid on
title('Imax (red) and Imin (blue) function of the order m')
figure
plot(m,V,'bd-'),grid on
title('visibility (first mode) function of the order m'),figure
dd=1000*lambda/2
d1=0:1.5:dd;
Ad=abs(cos(d1*deltak));
Bd=cos(2*k*d1);
Ir=4*(1+Ad.*Bd);
plot(d1,Ir,'ro-'),grid on
title('Whole intensity function of the distance d (in micron)',figure
plot(d1,Ad,'bd-'),grid on
title('visibility (second mode) function of the distance d (in micron)'
%==============================================
%
Top
Script S4_2_6B.m
Download Script S4_2_6B.m
%==============================================
%Michelson 2B
%a major number m of fringes (5001)
%and distance d (about 1473 microns)
%are considered
%==============================================
%
%wavelengths for a sodium-vapor lamp
lambda1=0.58900;
lambda2=0.58958;
lambda=(lambda1+lambda2)/2;
k=2*pi/lambda;
deltak=(2*pi)*((lambda2-lambda1)/(lambda1*lambda2));
m=0:5000;
max=length(m)
for i=1:max
%maximum of intensity
dmax=m(i)*lambda/2;
arg1max=deltak*dmax;
Amax=cos(arg1max);
arg2max=2*k*dmax;
Bmax=cos(arg2max);
Imax(i)=4*(1+Amax*Bmax);
%minimum of intensity
dmin=(2*m(i)+1)*lambda/4;
arg1min=deltak*dmin;
Amin=cos(arg1min);
arg2min=2*k*dmin;
Bmin=cos(arg2min);
Imin(i)=4*(1+Amin*Bmin);
%visibility
V(i)=abs((Imax(i)-Imin(i)))/(Imax(i)+Imin(i));
end
dmax
dmin
Imax;
Imin;
V;
%
plot(m,Imax,'r-',m,Imin,'b-'),grid on
title('Imax (red) and Imin (blue) function of the order m')
figure
plot(m,V,'b-'),grid on
title('visibility (first mode) function of the order m'),figure
dd=5000*lambda/2
d1=0:1.5:dd;
Ad=cos(d1*deltak);
Bd=cos(2*k*d1);
Ir=4*(1+Ad.*Bd);
plot(d1,Ir,'k-'),grid on
title('Whole intensity function of the distance d (in micron)'),figure
% axis([-10,500,0,8])
plot(d1,Ad,'k-'),grid on
title('visibility (second mode) function of the distance d (in micron)')
% axis([-10,300,0,8])
%==============================================
%
Top