使用MATLAB的“int”两次:
信谊一个(x1, y1 x y zλ
f = exp (- (x1 ^ 2 + x2 ^ 2) / a ^ 2) * exp(我*π/(λ* z) * ((x-x1) ^ 2 + (y-y1) ^ 2));
E = int (int (f, x1, 3 * 3 * 1), x2, 3 *, 3 * 1)
另外:
信谊一个(x1, y1 x y zλ
f1 = exp (x1 ^ 2 / ^ 2 + 1我*π/(λ* z) * (x-x1) ^ 2);
f2 = exp (y1 ^ 2 / ^ 2 + 1我*π/(λ* z) * (y-y1) ^ 2);
E = int (f1, x1, 3 * 3 *) * int (f2 x2, 3 * 3 * 1)
最好的祝愿
Torsten。