[MATLAB毕业设计] matlab作圆孔的菲涅耳衍射

[lwb6677]

最近在编一程序，想用傅立叶变换直接作圆孔的菲涅耳衍射，但总是不行，望高手指点一下。
a = 2; %圆孔直径
N = 1024;
d = 2/N; %间距
k = -N/2:N/2-1;
[x0,y0] = meshgrid(k*d);
ld=1.055e-3; %波长
z=700; %衍射距离
%===画一圆孔===
Ein = ones(N,N);
for m=1:N
for n=1:N
if (x0(m,n)^2+y0(m,n)^2>(a/2)^2)
Ein(m,n) = 0;
end
end
end

E = Ein.*exp(i*pi*(x0.^2+y0.^2)/ld/z); %傅立叶变换核
Eout = fftshift(fft2(E));
figure,imshow(abs(Eout));

f = 1/a; %观察面频率间隔
[Fx,Fy] = meshgrid(k*f);
x = ld*Fx;
y = ld*Fy; %观察面空间坐标
figure,mesh(x,y,abs(Eout))
主要没搞清楚的地方有：衍射面的坐标（x0,y0）和观察面的坐标（x,y）间如何联系起来？傅立叶变换后的频谱坐标和空间坐标的关系？还有就是是否有坐标缩放的问题？希望大牛们不吝赐教。