用matlab怎么求衍射后的涡旋光的轨道角动量谱（OAM谱）

类似于图片中的

主要不知道该怎么用数值积分写下列公式,

该怎样用积分函数quad写啊

E是什么表达式？

syms x y z k x0 y0 z0 a
E0 = 1;
sigma = 1;
l = 1;
E = E0^abs(l)*exp(-(x^2+y^2)/sigma^2)...
    *exp(1i*l*atan(y/x));
F = exp(1i*k/(2*z0)*((x0-x)^2-(y0-y)^2));
E = int(int(E*F,x,-a/2,a/2),y,-inf,inf)

T题主你好，我也遇到了这个困难，不知道如可对他进行傅里叶变换，请问您解决了吗

有看到这种类似的程序，如果有需要可联系我

题主解决了吗？我也有这样的疑问