菲涅耳衍射matlab,matlab – 两步菲涅耳衍射

我写了一个简短的matlab脚本文件,假设运行Fresnel的传播(衍射),这样给定一个输入字段U0,它会告诉你字段在距离z0之后的样子.我将结果与教科书结果进行了比较,看起来我的程序运行正常.问题是如果我尝试采取两个传播步骤而不是一个.即,不是采用程序的单次迭代来传播距离z0,而是采用程序的两次迭代来传播每个距离z0 / 2.然后我完全胡说八道,我无法弄清楚问题是什么.我们将非常感激地接受任何建议.

这是代码：

function U = fresnel_advance (U0, dx, dy, z, lambda)

% The function receives a field U0 at wavelength lambda

% and returns the field U after distance z, using the Fresnel

% approximation. dx, dy, are spatial resolution.

k=2*pi/lambda;

[ny, nx] = size(U0);

Lx = dx * nx;

Ly = dy * ny;

dfx = 1./Lx;

dfy = 1./Ly;

u = ones(nx,1)*((1:nx)-nx/2)*dfx;

v = ((1:ny)-ny/2)'*ones(1,ny)*dfy;

O = fftshift(fft2(U0));

H = exp(1i*k*z).*exp(-1i*pi*lambda*z*(u.^2+v.^2));

U = ifft2(O.*H);