模擬夫瑯禾費衍射的matlab源代碼

1、 源代碼：N=512; disp(衍射孔徑類型 1.圓孔 2.單縫 3.方孔)kind=input(please input 衍射孔徑類型：);% 輸入衍射孔徑類型while kind=1&kind=2&kind=3 disp(超出選擇范圍，請重新輸入衍射孔徑類型); kind=input(please input 衍射孔徑類型：);% 輸入衍射孔徑類型 endswitch(kind)case 1 r=input(please input 衍射圓孔半徑(mm)：);% 輸入衍射圓孔的半徑I=zeros(N,N);m,n=meshgrid(linspace(-N/16,N/16-1,N);D=(

2、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;% 單縫的長度I=zeros(N,N);m,n=meshgrid(linspace(-N/4,N/4,N);I(-am&ma&-bn&nb)=1;subplot(1,2,1);imshow(I);title(生成的衍射單縫); case 3a=input(please input 方孔邊長：);% 輸入方孔邊長I=zeros(N,N);m,n=meshg

3、rid(linspace(-N/4,N/4,N);I(-a/2m&ma/2&-a/2n&na/2)=1;subplot(1,2,1),imshow(I);title(生成的方孔); otherwise kind=input(please input 衍射孔徑類型：);% 輸入衍射孔徑類型end % 夫瑯禾費衍射的實現(xiàn)過程 L=500;x,y=meshgrid(linspace(-L/2,L/2,N);lamda_1=input(please input 衍射波長(nm)：);% 輸入衍射波長;lamda=lamda_1/1e6k=2*pi/lamda;z=input(please input

4、衍射屏距離衍射孔的距離(mm)：);% 衍射屏距離衍射孔的距離h=exp(1j*k*z)*exp(1j*k*(x.2+y.2)/(2*z)/(1j*lamda*z);%脈沖相應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(夫瑯禾費衍射強度分布) 實驗輸入：衍射孔徑類型 1.圓孔 2.單縫 3.方孔please input 衍射孔徑類型：1please input 衍射圓孔半徑(mm)：3please input 衍射波長(nm)：632lamda =6.3200e-04please input 衍射屏距離衍射孔的距離(mm)：1000000實驗結(jié)果：程序說明：本實驗可以選擇孔徑類型、孔徑半徑、輸入波長、衍射屏和衍射孔的距離等。當衍射屏和衍射孔的距離相對較小時，此衍射為菲涅耳衍射，當距離相對較大時滿足夫瑯禾費衍射的條件，