圖1中p=3;J是第一類貝塞爾函數(shù);圖2中Re表示取實(shí)部,求和是n取-20:20;theta為積分因子預(yù)期生成的圖像如圖3所示,
下面是我自己編的一個(gè)程序,但是運(yùn)行結(jié)果與預(yù)期圖像差距很大,我編的程序生成的圖像如圖4所示,
希望各位大神幫幫忙,小妹現(xiàn)在這里說聲謝謝了。
下面附上我自己編的程序代碼,希望大神們指出問題所在:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%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)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
![[color=Black][font=宋體]matlab先求和嗎,再積分,求和中有積分因子[/font][/color]]()
圖片1.png
![[color=Black][font=宋體]matlab先求和嗎,再積分,求和中有積分因子[/font][/color]-1]()
圖片2.png
![[color=Black][font=宋體]matlab先求和嗎,再積分,求和中有積分因子[/font][/color]-2]()
圖片3.png
![[color=Black][font=宋體]matlab先求和嗎,再積分,求和中有積分因子[/font][/color]-3]()
圖4.jpg |