| 查看: 421 | 回复: 0 | |||
[交流]
【求助】如何循环求解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{:}); 求高手指点，谢谢 |
回复此楼» 猜你喜欢
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有13人回复
-
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有11人回复
-
表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了
已经有9人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有10人回复
-
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有11人回复