最近在编一程序，想用傅立叶变换直接作圆孔的菲涅耳衍射，但总是不行，望高手指点一下。
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）间如何联系起来？傅立叶变换后的频谱坐标和空间坐标的关系？还有就是是否有坐标缩放的问题？希望大牛们不吝赐教。