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）间如何联系起来？傅立叶变换后的频谱坐标和空间坐标的关系？还有就是是否有坐标缩放的问题？希望大牛们不吝赐教。

[laosam280]

没有坐标转换这种说法啊，如果你要画出衍射的图形，只要作出圆孔的大小，好友衍射的方程，建立初始点和衍射点的关系，然后在衍射成像的平面上画图就可以啦。傅里叶变换只变化了图形的形状，拓扑结构却不变化。
衍射的结果是应该这样的，我画了下。

[lwb6677]

初始点和衍射点的关系该如何建立呢？能将你画的图贴上来吗？

[wangzl]

衍射面的坐标（x0,y0）和观察面的坐标（x,y）间的联系：在你做2D傅里叶变换的时候，两者已经联系起来了，具体的关系式你可以分析Fresnel公式得到
傅立叶变换后的频谱坐标和空间坐标的关系：这一点也可以从Fresnel衍射公式中得到
至于坐标的缩放，没有这种说法，这取决于你要求的观察面的大小
另外，建议你使用命令mfun('FresnelC',x),mfun('FresnelS',x)计算你要求的衍射图样

[xuezhinvwang]

请问您编出来了吗。我现在也需要这个程序，能否共享一下，(*^__^*) 嘻嘻……