《matlab实现夫朗和费矩形和圆孔衍射.doc》由会员分享，可在线阅读，更多相关《matlab实现夫朗和费矩形和圆孔衍射.doc（4页珍藏版）》请在金锄头文库上搜索。
1、2、用MATLAB仿真平行光束的衍射强度分布图样。（夫朗和费矩形孔衍射、夫朗和费圆孔衍射、夫朗和费单缝和多缝衍射。）理论推导部分2. （1）夫朗和费矩形孔衍射若衍射孔为矩形则在透镜焦平面上得到的衍射图样如图，衍射图样的主要特征为衍射亮斑集中分布在两个相互垂直的方向上，并且x轴上的亮斑宽度与y轴亮斑宽度之比，恰与矩形孔在两个轴上的宽度相反。其中的为x,同样的中的为y,利用x=x/f,y=y/f进行求解。（2）夫朗和费圆形孔衍射夫朗和费圆孔衍射的讨论方法和矩形孔衍射的讨论方法相同，只是由于圆孔的几何对称性，采用极坐标更为方便。=ka2. （1）夫朗和费矩形孔衍射clear all;lamda=50
2、0e-9;a=1e-3;b=1e-3;f=1;m=500;ym=8000*lamda*f;ys=linspace(-ym,ym,m)xs=ys;n=255;for i=1:m sinth2=ys./sqrt(ys.2+f2);%相当于x/f sinth1=xs(i)/sqrt(xs(i).2+f2);%xs(i）作用每给一个ys值，要遍历到所有的x值 angleA=pi*a*sinth1/lamda;%相当于书上的alfa=kax/2f k=2*pi/lamda angleB=pi*b*sinth2./lamda;B(:,i)=(sin(angleA).2.*sin(angleB).2.*50
3、00./(angleA.2.*angleB.2);%光强度公式endsubplot(1,2,1)image(xs,ys,B)colormap(gray(n)subplot(1,2,2)plot(B(m/2,:),ys)（2）夫朗和费圆孔衍射clearlam=500e-9a=1e-3f=1m=300;ym=5*0.61*lam*f/a;%取爱里光斑半径的5倍ys=linspace(-ym,ym,m);xs=ys;n=200;for i=1:m r=xs(i)2+ys.2;%相当于r的平方 sinth=sqrt(r./(r+f2);%角度 fai=2*pi*a*sinth./lam;%fai=k*a*sinth hh=(2*BESSELJ(1,fai).2./fai.2;%贝塞尔函数 b(:,i)=hh.*5000;endsubplot(1,2,1)image(xs,ys,b)colormap(gray(n)subplot(1,2,2)b(:,m/2)plot(ys,b(:,m/2)
