資源描述:
《matlab實(shí)現(xiàn)夫朗和費(fèi)矩形和圓孔衍射》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、2、用MATLAB仿真平行光束的衍射強(qiáng)度分布圖樣。(夫朗和費(fèi)矩形孔衍射、夫朗和費(fèi)圓孔衍射、夫朗和費(fèi)單縫和多縫衍射。)理論推導(dǎo)部分2.(1)夫朗和費(fèi)矩形孔衍射若衍射孔為矩形則在透鏡焦平面上得到的衍射圖樣如圖,衍射圖樣的主要特征為衍射亮斑集中分布在兩個(gè)相互垂直的方向上,并且x軸上的亮斑寬度與y軸亮斑寬度之比,恰與矩形孔在兩個(gè)軸上的寬度相反。其中的θ為θx,同樣的β中的θ為θy,利用θx=x/f,θy=y/f進(jìn)行求解。(2)夫朗和費(fèi)圓形孔衍射夫朗和費(fèi)圓孔衍射的討論方法和矩形孔衍射的討論方法相同,只是由于圓孔的幾何對(duì)稱性,采用極坐標(biāo)更為方便。Ф=kaθ
2、2.(1)夫朗和費(fèi)矩形孔衍射clearall;lamda=500e-9;a=1e-3;b=1e-3;f=1;m=500;ym=8000*lamda*f;ys=linspace(-ym,ym,m)xs=ys;n=255;fori=1:msinth2=ys./sqrt(ys.^2+f^2);%相當(dāng)于x/fsinth1=xs(i)/sqrt(xs(i).^2+f^2);%xs(i)作用每給一個(gè)ys值,要遍歷到所有的x值angleA=pi*a*sinth1/lamda;%相當(dāng)于書上的alfa=kax/2fk=2*pi/lamdaangleB=pi*b*
3、sinth2./lamda;B(:,i)=(sin(angleA).^2.*sin(angleB).^2.*5000./(angleA.^2.*angleB.^2));%光強(qiáng)度公式endsubplot(1,2,1)image(xs,ys,B)colormap(gray(n))subplot(1,2,2)plot(B(m/2,:),ys)(2)夫朗和費(fèi)圓孔衍射clearlam=500e-9a=1e-3f=1m=300;ym=5*0.61*lam*f/a;%取愛里光斑半徑的5倍ys=linspace(-ym,ym,m);xs=ys;n=200;fo
4、ri=1:mr=xs(i)^2+ys.^2;%相當(dāng)于r的平方sinth=sqrt(r./(r+f^2));%角度fai=2*pi*a*sinth./lam;%fai=k*a*sinthhh=(2*BESSELJ(1,fai)).^2./fai.^2;%貝塞爾函數(shù)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))