24小时热门版块排行榜

返回列表

pphuaxia

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.5
帖子: 2
在线: 55分钟
虫号: 928075
注册: 2009-12-14
专业: 流体力学

[交流] 【求助】如何循环求解bessel方程

我用的是2010版的，想要循环求解bessel方程，代码如下：
syms b n
M=10;
lamda1 = 1000;
lamda2 = 1200;
lamda = linspace(lamda1,lamda2,M)*1e-9;
deltalamda=(lamda2-lamda1)/M*1e-9;
tic
for k = 1:M
V(k)=2*pi/lamda(k)*a1*sqrt(n1^2-n2^2)
eq(k)=V(k)*sqrt(1-b(k))*besselj(1,V(k)*sqrt(1-b(k)))/besselj(0,V(k)*sqrt(1-b(k)))-V(k)*sqrt(b(k))*besselk(1,V(k)*sqrt(b(k)))/besselk(0,V(k)*sqrt(b(k)))
b(k)=double(solve(eq(k)))
n(k)=(b(k)*n1^2 - b(k)*n2^2 + n2^2)^(1/2)
end
但是算不出结果,提示说：
??? Error using ==> mupadmex
Error in MuPAD command: Index exceeds matrix dimensions.

Error in ==> sym.sym>sym.subsref at 1366
B = mupadmex('mllib::subsref',A.s,inds{:});
求高手指点，谢谢

回复此楼

» 猜你喜欢

1楼 2010-06-09 10:53:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 pphuaxia 的主题更新

返回列表