此MATLAB代码对应的数学公式% 计算相位因子 P =E0.*exp(-1i * psi); % 进行傅里叶变换 Ef = fftshift(fft2(P)); % 计算传输函数 H = exp(-1i * k * z) ./ (1i * lambda * z) .* exp(1i * k / (2 * z) .* (rho.^2)); % 传输 Ef = Ef .* H; % 反傅里叶变换 E = ifft2(ifftshift(Ef));

这段代码对应的数学公式如下：

- 相位因子：$P(x,y) = E_0(x,y) e^{-i \psi(x,y)}$
- 傅里叶变换：$\tilde{E}(f_x,f_y) = \mathcal{F}[P(x,y)]$
- 传输函数：$H(f_x, f_y) = \frac{e^{-i k z}}{i \lambda z} e^{i k \frac{\rho^2}{2 z}}$
- 传输：$\tilde{E}(f_x,f_y) = \tilde{E}(f_x,f_y) H(f_x,f_y)$
- 反傅里叶变换：$E(x,y) = \mathcal{F}^{-1}[\tilde{E}(f_x,f_y)]$