[求助] [color=Black][font=宋体]matlab先求和吗，再积分，求和中有积分因子[/font][/color]

图1中p=3；J是第一类贝塞尔函数；图2中Re表示取实部，求和是n取-20:20；theta为积分因子预期生成的图像如图3所示， 下面是我自己编的一个程序，但是运行结果与预期图像差距很大，我编的程序生成的图像如图4所示， 希望各位大神帮帮忙，小妹现在这里说声谢谢了。 下面附上我自己编的程序代码，希望大神们指出问题所在： %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear,clc; rho=3;% fractional vortex x=0; y=0; for n=-20:20 a=1; for S=-3:0.2:3 b=1;     for alpha=-3:0.2:3         psi(a,b)=(n/eps(n-alpha))*sqrt(pi/8)*(-i).^(abs(n)/2).*rho.*exp(i*rho.^2/4).*...     (besselj((abs(n)-1)/2,rho.^2/4)-i.*besselj((abs(n)+1)/2,rho.^2/4)).*...     quad(@(theta) exp(i.*n.*theta),0,2*pi);     x=psi(a,b)+x;     psi(a,b)=x;     psi2(a,b)=sqrt(pi/8)*(-i).^(abs(n)/2).*rho.*exp(i*rho.^2/4).*...      (besselj((abs(n)-1)/2,rho.^2/4)-i.*besselj((abs(n)+1)/2,rho.^2/4)).*...      quad(@(theta) exp(i.*n.*theta),0,2*pi)./eps(n-alpha);    y=psi2(a,b)+y;    psi2(a,b)=y;    b=b+1;     end     a=a+1; end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%psi3=x/y; psi4=real(psi3)/(2*pi); plot(psi4) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 图片1.png 图片2.png 图片3.png 图4.jpg

回复此楼