matlab实现夫朗和费矩形和圆孔衍射.docxVIP

2、用MATLA仿真平行光束的衍射强度分布图样。（夫朗和费矩形孔 衍射、夫朗和费圆孔衍射、夫朗和费单缝和多缝衍射。） 理论推导部分 2.（ 1）夫朗和费矩形孔衍射 若衍射孔为矩形则在透镜焦平面上得到的衍射图样如图， 衍射图样的 主要特征为衍射亮斑集中分布在两个相互垂直的方向上， 并且x轴上 的亮斑宽度与y轴亮斑宽度之比，恰与矩形孔在两个轴上的宽度相 反。 2o 2 o (sin a \ f sin /7 “心 ^)=儿 | I分丿I 〃 ) a = kax/(2 f) = (na/ A)(x/ f) = (nasin 0)/ A 其中的B为0x,同样的B中的B为0y,利用 仅二x/f, 0y=y/f进行求解 （2）夫朗和费圆形孔衍射 夫朗和费圆孔衍射的讨论方法和矩形孔衍射的讨论方法相同， 只是由 于圆孔的几何对称性，采用极坐标更为方便。 22 r 2 2 r 2 J (kaO)T- 7 「2 / (少)- ) c I 1 =人 p L “」 ①=ka 9 2. (1)夫朗和费矩形孔衍射 clear all; lamda=500e-9; a=1e-3; b=1e-3; f=1; m=500; ym=8000*lamda*f; ys=li nspace(-ym,ym,m) xs=ys; n=255; for i=1:m si nth2二ys./sqrt(ys42+fA2);% sin th1=xs(i)/sqrt(xs(i).A2+fA2);%xs(i 相当于x/f）作用每给一个an gleA二pi*a*si nthl/lamda;%相当于书上的alfa=kax/2fys值，要遍历到所有的x值 相当于x/f ）作用每给一个 an gleA二pi*a*si nthl/lamda;% 相当于书上的alfa=kax/2f an gleB=pi*b*s in th2./lamda; B(:,i)=(si n(an gleA).A2.*si n(an gleB).A2.*5000./(a ngleA.A2.*a ngleB.A2));%光强度公式 end subplot(1,2,1) image(xs,ys,B) colormap(gray (n)) subplot(1,2,2) plot(B(m/2,:),ys) (2)夫朗和费圆孔衍射 clear lam=500e-9 a=1e-3 f=1 m=300; ym=5*0.61*lam*f/a;%取爱里光斑半径的5倍 ys=li nspace(-ym,ym,m); xs=ys; n=200; for i=1:m r=xs(i)A2+ys.A2;% 相当于r的平方 sin th=sqrt(r./(叶fA2));% 角度 fai=2*pi*a*s in th./lam;%fai=k*a*si nth hh=(2*BESSELJ(1,fai))A2./fai.A2;% 贝塞尔函数 b(:,i)=hh.*5000; end subplot(1,2,1) image(xs,ys,b) colormap(gray (n)) subplot(1,2,2) b(:,m/2) plot(ys,b(:,m/2))