Question

我想绘制函数 v 的等值面。

问题是 v 不接受矩阵参数（不是矢量化函数它包含Laguerre相关多项式）

因此，如果我为x，y，l

的值创建一个meshgrid

[r,th,l]=meshgrid(0:0.5:5,0:pi/2:2*pi,0:2);

x=r.*cos(th);
y=r.*sin(th);

然后我想用循环（计数器）评估 v ，因为我无法绕过 x，y＆amp; l作为参数

但我做错了什么，结果我不是评估网格中的 v ：

for l=0:2
k=k+1;
    for r=0:0.5:5
    i=i+1;
        for th=0:pi/2:2*pi
        j=j+1;      


fun1=@(R)4*real(exp(-r.^2-R^2+2*1i*R*l./r).*(r+1i*R).^(2*l).*...
  (mfun('L',n/2-l/2,l,r.^2+R^2)).^2);

v(i,j,k)=integral(fun1,-inf,inf);

        end
    end  
end

是否有人知道如何进行循环以便我可以获得与x，y，l大小匹配的 v ，以便我可以使用：

isosurface（x，y，l，v）??

或者是否有人知道如何通过其他方式获取上述isosurface v ？

我真的需要我能得到的所有帮助：）

的问候，

巴勃罗

评估网格中的函数，matlab

0 个答案: