| 查看: 429 | 回复: 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{:}); 求高手指点，谢谢 |
回复此楼» 猜你喜欢
一志愿南京航空航天大学 材料与化工329分求调剂
已经有7人回复
-
一志愿C9的化学工程（085602） 340分，感觉校内调剂无望，求调剂
已经有12人回复
-
312求调剂
已经有4人回复
-
材料专硕322
已经有9人回复
-
求调剂
已经有7人回复
-
085600材料与化工301分求调剂院校
已经有8人回复
-
一志愿 江南大学 085602 化工专硕 338分求调剂
已经有14人回复
-
070300化学学硕311分求调剂
已经有13人回复
-
材料调剂
已经有9人回复
-
一志愿河北工业大学材料工程，初试344求专硕调剂
已经有4人回复
