模擬夫瑯禾費(fèi)衍射的matlab源代碼

1、精品源代碼：N=512;disp( ' 衍射孔徑類型1.圓孔2.單縫3.方孔 ')kind=input('please input衍射孔徑類型：' );% 輸入衍射孔徑類型whilekind=1&kind=2&kind=3disp( ' 超出選擇范圍，請(qǐng)重新輸入衍射孔徑類型' );kind=input('please input衍射孔徑類型：');% 輸入衍射孔徑類型endswitch (kind)case 1r=input('please input衍射圓孔半徑(mm) ： ' );% 輸入衍射圓

2、孔的半徑I=zeros(N,N);m,n=meshgrid(linspace(-N/16,N/16-1,N);D=(m.2+n.2).(1/2);I(find(D<=r)=1;subplot(1,2,1),imshow(I);title( '生成的衍射圓孔');case 2a=input('please input衍射縫寬： ' );% 輸入衍射單縫的寬度b=1000;% 單縫的長(zhǎng)度I=zeros(N,N);-可編輯 -精品m,n=meshgrid(linspace(-N/4,N/4,N);I(-a<m&m<a&-b<n&

3、amp;n<b)=1;subplot(1,2,1);imshow(I);title( '生成的衍射單縫');case 3a=input('please input方孔邊長(zhǎng)： ' );% 輸入方孔邊長(zhǎng)I=zeros(N,N);m,n=meshgrid(linspace(-N/4,N/4,N);I(-a/2<m&m<a/2&-a/2<n&n<a/2)=1;subplot(1,2,1),imshow(I);title( '生成的方孔 ');otherwisekind=input('pleas

4、e input衍射孔徑類型：' );% 輸入衍射孔徑類型end% 夫瑯禾費(fèi)衍射的實(shí)現(xiàn)過(guò)程L=500;x,y=meshgrid(linspace(-L/2,L/2,N);lamda_1=input('please input衍射波長(zhǎng) (nm) ： ');% 輸入衍射波長(zhǎng);lamda=lamda_1/1e6k=2*pi/lamda;z=input('please input衍射屏距離衍射孔的距離(mm) ： ' );%衍射屏距離衍射孔的距離h=exp(1j*k*z)*exp(1j*k*(x.2+y.2)/(2*z)/(1j*lamda*z);% 脈沖相應(yīng)-可

5、編輯 -精品H =fftshift(fft2(h);% 傳遞函數(shù)B=fftshift(fft2(I);% 孔頻譜G=fftshift(ifft2(H.*B);subplot(1,2,2),imshow(log(1+abs(G),);title( '衍射后的圖樣');figuremeshz(x,y,abs(G);title( '夫瑯禾費(fèi)衍射強(qiáng)度分布')實(shí)驗(yàn)輸入：衍射孔徑類型1.圓孔2.單縫3.方孔please input衍射孔徑類型：1please input衍射圓孔半徑 (mm) ： 3please input衍射波長(zhǎng) (nm) ： 632lamda =6.3200e-04please input衍射屏距離衍射孔的距離(mm) ： 1000000實(shí)驗(yàn)結(jié)果：-可編輯 -精品程序說(shuō)明：本實(shí)驗(yàn)可以選擇孔徑類型、孔徑半徑、輸入波長(zhǎng)、衍射屏和衍射孔的距離等。當(dāng)衍射屏和衍射孔的距離相對(duì)較小時(shí)，此衍射為菲涅耳衍射，當(dāng)距離相對(duì)較大時(shí)滿足夫瑯禾費(fèi)