1

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

问题是v不接受矩阵参数（不可向量化函数，因为它包含拉盖尔相关多项式）

所以如果我为 x,y,l 的值创建一个网格网格

[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 & 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 ，因此我可以使用：

等值面(x,y,l,v) ??

或者有人知道如何通过另一种方式获得提到的等值面v ？

我真的需要我能得到的所有帮助:)

问候，

巴勃罗

4

0 回答 0