应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » [color=Black][font=宋体]matlab先求和吗，再积分，求和中有积分因子[/font][/color]
21/1返回列表
查看: 462  |  回复: 1
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

几度

银虫 (初入文坛)

[求助] [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)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
天行健，君子自强不息
1楼2013-11-18 11:33:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

几度

银虫 (初入文坛)
已经运行出来了
回复此楼
天行健，君子自强不息
2楼2013-11-25 09:26:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 几度 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭