未定义函数或变量 'k'。 出错 propTF (line 8) H = exp(1j*k*z)*exp(-1j*pi*lambda*z*(X.^2+Y.^2)); 出错 Untitled (line 22) E1 = propTF(E0,L,lambda,z);

非常抱歉，我在之前的代码中漏掉了定义波数 `k` 的代码。请在代码开头添加以下代码：

% 定义常数
lambda = 632.8e-9; % 波长
k = 2*pi/lambda; % 波数
z = 10e-2; % 光程差
d = 1e-3; % 孔径大小
L = 0.1; % 平面波前大小
N = 1024; % CCD像素数目
dx = L/N; % 像素大小

将其中的 `lambda` 和 `z` 值替换为实际值即可。再次非常抱歉造成不便。