| 查看: 422 | 回复: 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{:}); 求高手指点，谢谢 |
回复此楼» 猜你喜欢
材料与化工一志愿南昌大学327求调剂推荐
已经有8人回复
-
0703 物理化学调剂
已经有5人回复
-
化学工程321分求调剂
已经有7人回复
-
0703化学调剂 290分有科研经历，论文在投
已经有7人回复
-
308求调剂
已经有3人回复
-
本人考085602 化学工程 专硕
已经有10人回复
-
焦虑
已经有9人回复
-
0856求调剂
已经有3人回复
-
283求调剂
已经有11人回复
-
深圳大学硕士招生（2026秋，传感器方向，仅录取第一志愿）
已经有8人回复
