《多缝的夫琅和费衍射matlab仿真》由会员分享，可在线阅读，更多相关《多缝的夫琅和费衍射matlab仿真（5页珍藏版）》请在金锄头文库上搜索。
1、4、多缝的夫琅和费衍射，使用平行光照明，观察衍射图样随点光源位置（光源 上下移动）的变化图4-1图4-2多缝夫琅禾费衍射如图4-1所示。由于相邻单缝在P点产生的夫琅禾费衍射2 冗.的幅值与中心单缝的相同，只是产生一个相位差5= Udsin9，故，经证明，P 人点处的光强为:任sin5sin 2I(P) = Q)2()2，0其中以=：a sin 9， 5=绑 d sin 9。人人因而，程序代码如下：clearLambda=600*(1e-9);a=0.005*(1e-3);0.005mm，0.02mmd=0.02*(1e-3);f=0.01;N=20;ni=1000;x=linspace(-0.
2、005,0.005,ni); for k=1:nisn=x(k)/sqrt(x(k)2 + f八2);alpha=pi*a*sn/Lambda;%清除原有变量%设置波长为600nm%设置衍射屏参数：缝宽为0.005mm，缝距为0.02mm%汇聚透镜焦距设置为1cm%设置缝数为2 0%将衍射屏按照狭缝方向分为ni个微元%算各微元对应的a和0值delta=2*pi*d*sn/Lambda;I(k)=(sin(alpha)/alpha).八2*(sin(N*delta/2)/sin(delta/2).八2;%求出各处的光强endfigure(gcf);% 显示图像NCLevels=250;Br=I*
3、NCLevels;image(0,x,Br);colormap(gray(NCLevels);title(二维强度分布);运行后结果如图4-2所示。将光源上下移动的结果如图4-3所示:F：图4-3图4-4点光源发出的光经过准直透镜后形成倾斜入射的平行光，倾斜角度为i。此 时，P点强度的公式为：I ( P ) = * A/)2，0以 P其中以=#(sin9 -sini)，6 =孚。- sini)。故而程序代码如下，当i。0时，即光源上下移动，改变其的值，即可仿真出 移动不同距离时的衍射图样。程序代码如下：清除原有变量设置%波长为600nm%设置衍射屏参数：缝宽为0.005mm，clearLamb
4、da=600*(1e-9);a=0.005*(1e-3);缝距为0.02mmd=0.02*(1e-3);f=0.01;N=10;ni=500;i=pi/20;x=linspace(-0.005,0.005,ni); for k=1:nisn=x(k)/sqrt(x(k)2 + f八2);alpha=pi*a*(sn-sin(i)/Lambda;delta=2*pi*d*(sn-sin(i)/Lambda;I(k)=(sin(alpha)/alpha).八2*(sin(N*delta/2)/sin(delta/2).八2;%求出各处的光强%汇聚透镜焦距设置为1cm 设置缝数为10设置平行光入射的
5、倾斜角%将衍射屏按照狭缝方向分为ni个微元%算各微元对应的a和0值endfigure(gcf);NCLevels=250;Br=I*NCLevels;image(0,x,Br);colormap(gray(NCLevels);title(二维强度分布);显示圈像运行程序后结果如图4-4。令i = L，结果如图4-5；令，=-，结果如图4-6。1010故光源上下移动，会使干涉图样分别向右边和左边偏移。清除原有变量设置波长为600nm%设置衍射屏参数：缝宽为0.005mm，缝距为0.02mm%汇聚透镜焦距设置为1cm设置缝魏为10将%扩展光源分为微元%将衍射屏按照狭缝方向分为ni个微元%算各微元对
6、应的a和0值2;%求出各处的光强将各光源洸强累加显示图像2*(sin(N*delta/2)/sin(delta/2)5、多缝的夫琅和费衍射，使用扩展光源照明，前后、上下移动光源，观察衍射 图样的变化使用扩展光源照明，即为4中点光源的光照的叠加。程序代码为：clearLambda=600*(1e-9);a=0.005*(1e-3);d=0.02*(1e-3);f=0.01;N=10;ni=500;I=zeros(1,ni);i=linspace(-pi/1000,pi/1000,200);x=linspace(-0.005,0.005,ni);for t=1:200for k=1:nisn=x(
7、k)/sqrt(x(k)2 + f八2);alpha=pi*a*(sn-sin(i(t)/Lambda;delta=2*pi*d*(sn-sin(i(t)/Lambda;In(k)=(sin(alpha)/alpha).endI=I+In;endfigure(gcf);NCLevels=250;Br=I*NCLevels;image(0,x,Br);colormap(gray(NCLevels);title(二维强度分布);运行结果如图5-1，为左右偏移后的干涉条纹的叠加。二维强度分布二维强度分布图5-1图5-2将光源上下移动，可以通过改变i的取值实现，将上面程序中i=linspace(-pi/1000,pi/1000,200);将扩展光源分为微元 改为i=linspace(-pi/2000,3*pi/2000,200);%将扩展光源分为微元运行后可得结果如图5-2,认为干涉条纹的叠加。若扩展光源足够大，令i=linspace(-pi/2000,3*pi/2000,200);微元则结果如图5-3，为全亮。%将扩展光源分为图5-3
