特殊孔径的衍射模拟图像

sunbader

请教各位好汉：
          我从网上找了一个所谓的通用程序,用来画任意特殊孔径的衍射图像的。我试了一下，可是结果老是不对，不知道问题出在哪里。请大家帮忙看看，互相讨论一下。谢谢

A=imread('file','format')     %file为图片名称，format为图片格式
imshow(A)                    %显示图片
B=fft2(A);C=abs(B).^2;
D=fftshift(C);
E=max(D(:));
G=D/E;                      %光强归一化
figure,
image(sqrt(G))          %灰度非线性处理
axis square;  
axis off

下面是孔径图像和衍射结果。孔径的图像是我用windows自带的画图工具随便画的。
K:\matlab\file.bmp
K:\matlab\jieguo.bmp

sunbader

sunbader