1、利用MATLAB進行夫瑯和費衍射我已經發過相關的帖子，是我以前做過的課程論文。近來看見有很多人回帖說需要程序，故而總結一下方法和共享程序：通過MATLAB軟件編程實現夫瑯和費衍射的方法：(1 )用衍射積分(2)傅立葉變換一、衍射積分相關程序如下：1. 單縫衍射clearlamba=500e-9;% 波長a=1e-3;D=1;ym=3*lamba*D/a;%屏幕上y的范圍n=51;%屏幕上的點數ys=li nspace(-ym ,ym,n);n=51;%屏幕上的點數yp=li nspace(O,a ,n);for i=1: nsin phi=ys(i)/D;alpha=pi*yp*s in ph

2、i/lamba;sumcos=sum(cos(alpha);sums in=sum(s in( alpha);B(i,:)=(sumcosA2+sumsi 門人2)/門人2;endN=256;%確定灰度的等級Br=(B/max(B)*N;subplot(1,2,1)image(ym,ys,Br);colormap(gray(N);% 色調處理subplot(1,2,2)plot(B,ys,'k');陽件*單齢射例K if-0.4-0.20020.41-1 H2. 多縫衍射clearlamda=500e-9; % 波長N=2; %縫數，可以隨意更改變換a=2e-4;D=5;d=5

3、*a; ym=2*lamda*D/a;xs=ym; n=1001;ys=li nspace(-ym ,ym,n);for i=1: nsin phi=ys(i)/D;alpha=pi*a*s in phi/lamda;beta=pi*d*s in phi/lamda;B(i,:)=(s in (alpha)./alpha).A2.*(s in (N*beta)./sin (beta).A2;B仁 B/max(B);endNC=256; %確定灰度的等級Br=(B/max(B)*NC; subplot(1,2,1) image(xs,ys,Br);colormap(gray(NC); %色調處理

4、subplot(1,2,2)plot(B1,ys,'k');-0.0250.025雙縫仙.jpg0.02-0.0160.015.01001-0 0050 0050.01-0.010.00500.G15-0.0150.02-0.020 025 -0.4 -0 200 20 43. 矩孔衍射 clear 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:msin th1=xs(i)/sqrt(xs(i)A2+fA2);X 103x 10si

5、nth2=ys./sqrt(ys.A2+fA2);an gleA=pi*a*si nthl/lamda;an gleB=pi*b*s in th2./lamda;B(:,i)=(si n(a ngleA)42.*si n(a ngleB)42.*5000./ (an gleA.A2.*a ngleB.A2);end subplot(1,2,1) image(xs,ys,B) colormap(gra y(n) subplot(1,2,2) plot(B(m/2,:),ys,'k') 屮4. 正弦光柵clear all xm=10*pi;ys=xm;xs=li nspace(-x

6、m,xm,500);B=cos(xs)+1;N=255;Br=B/2*N;image(xs,ys,Br);colormap(gray(N);-0.025-002-0.0150.010.D050.010.015002-0,0°025-0 400.20.4二、傅里葉變換(1) 基本思想：在傅立葉變換光學中夫瑯和費衍射場的強度分布就 等于屏函數的功率譜。因此我們可以直接將光屏進行傅 立葉變換，再處理得到衍射圖樣。(2) 具體操作：首先，物體圖像的生成：(1)直接由矩陣運算生成(2)也可 利用Windows下的畫圖工具，生成一幅黑白圖像。其次， 并調用命令函數imread ()輸入圖像，輸入

7、的圖像是一 個巨大的二維矩陣，利用MATLAB函數庫中的fft2 ()命令對該矩陣進行二位離散傅立葉變換，得到圖像的頻譜，該頻譜是一個復數矩陣，然后用取模函數abs ()對該復數矩陣取模，得到振幅譜矩陣，利用函數fftshift ( )對取模后的矩陣進行頻譜位移，這是因為變換后的二 維矩陣的直流分量位于圖像的周邊角，該函數交換矩陣 的1、3象限和2、4象限，使直流分量移到頻譜中心，從 而使FFT頻譜可視效果與實際圖像相吻合。最后利用 imshow ()函數將圖像顯示出來。*具體圖像參考：基本的程序：cleara=imread('*');%你的圖像的位置或直接編寫孔生成grid onfigure(1)imshow(a,)afft=fft2(a);aabs=abs(afft);aabss=fftshift(aabs);figure(2)i