matlab做干涉,matlab模拟干涉衍射

%##########################%

%名称：等厚干涉模拟 %

%作者：程序猿__int64Ago %

%时间：2011-12-16 %

%##########################%

clear %清内存

lamda=632.8e-9; %设定入射光波长

n=1; %缝隙折射率

R=5; %柱透镜半径

ym=5e-3; %设定屏幕的范围,根据实际情况调整

xs=ym; %x方向无变化，单一值

N=1000; %条纹柔和度

ys=linspace(-ym,ym,N); %y轴上坐标离散

for i=1:N %依次求干涉光强

hm=R-sqrt(R^2-ys(i)^2); %柱透镜与底面间隙

delta=2*n*hm+lamda/2; %光程差

phi=2*pi*delta/lamda; %相位差

I(i,:)=4*cos(phi/2).^2; %干涉光强

end %求光强结束

M=255; %设定灰度等级

II=(I/4.0)*M; %设定黑白条纹对比度

subplot(1,2,1); %分割figure1×2，并且在第一块上作图

image(xs,ys,II); %画干涉条纹

colormap(gray(M)); %二维光强分布灰度图

title('光强分布灰度图') %列出标题

subplot(1,2,2); %分割figure1×2，并且在第二块上作图

plot(I,ys); %一维光强分布曲线图

title('光强分布曲线图') %列出标题

%##########################%

%名称：多缝夫朗禾费衍射模拟 %

%作者：程序猿__int64Ago %

%时间：2011-12-16 %

%##########################%

clear %清内存

lamda=500e-9; %设定入射光波长

a=10e-6; %缝宽

d=25e-6; %光栏常数

N0=6; %缝数

D=5; %缝到像面的距离(这个随意，几乎不影响)

ym=1*lamda*D/a; %设定屏幕的范围,根据实际情况调整

xs=ym; %x方向无变化，单一值

N=1000; %条纹柔和度(不能太大)

ys=linspace(-ym,ym,N); %y轴上坐标离散

for i=1:N %依次求衍涉光强

sintheta=ys(i)/D; %以下由书上公式

alpha=pi*a*sintheta/lamda;

beta=pi*d*sintheta/lamda;

I(i,:)=(sin(alpha)./alpha).^2.*...

(sin(N0*beta)./sin(beta)).^2; %上一行续行

end %求光强结束

M=255; %设定灰度等级

II=(I/max(I))*M; %设定黑白条纹对比度

subplot(1,2,1); %分割figure1×2，并且在第一块上作图

image(xs,ys,II); %画干涉条纹

colormap(gray(M)); %二维光强分布灰度图

title('光强分布灰度图') %列出标题

subplot(1,2,2); %分割figure1×2，并且在第二块上作图

plot(I,ys); %一维光强分布曲线图

title('光强分布曲线图') %列出标题