資源描述：
《matlab實(shí)現(xiàn)夫朗和費(fèi)矩形和圓孔衍射》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。
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)圓孔衍射的討論方法和矩形孔衍射的討論方法相同，只是由于圓孔的
2、幾何對(duì)稱(chēng)性，采用極坐標(biāo)更為方便。Ф=kaθ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;
3、%相當(dāng)于書(shū)上的alfa=kax/2fk=2*pi/lamdaangleB=pi*b*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
4、=5*0.61*lam*f/a;%取愛(ài)里光斑半徑的5倍ys=linspace(-ym,ym,m);xs=ys;n=200;fori=1:mr=xs(i)^2+ys.^2;%相當(dāng)于r的平方sinth=sqrt(r./(r+f^2));%角度f(wàn)ai=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
5、(n))subplot(1,2,2)b(:,m/2)plot(ys,b(:,m/2))
