level 1
六指琴魔的琴
楼主
代码如下
f[Subscript[k, i] _,
Subscript[\[Omega],
i] _] = ((\[Omega] - k*I)^2*BesselI[m, Subscript[n, 1]])/(
Subscript[n, 1]*
D[BesselI[m, Subscript[n, 1]], Subscript[n, 1]]) - (
Q*\[Omega]^2*BesselK[m, Subscript[n, 2]])/(
Subscript[n, 2]*D[BesselK[m, Subscript[n, 2]], Subscript[n, 2]]) +
We*(k^2 + m^2 - 1) /. {Subscript[n,
1] -> (k^2 + Subscript[Ma, 1]^2*(\[Omega] - k*I)^2)^0.5,
Subscript[n,
2] -> (k^2 + Subscript[Ma, 2]^2*\[Omega]^2 )^0.5} /. {\[Omega] ->
Subscript[\[Omega], r] + Subscript[\[Omega], i]*I,
k -> Subscript[k, r] + Subscript[k, i]*I}
m = 0;
We = 8*10^-6;
Q = 0.0014;
Subscript[Ma, 1] = 0.2;
Subscript[Ma, 2] = 0.57;
Subscript[\[Omega], r] = 0;
Subscript[k, r] = 150;
FindRoot[{Re[f[Subscript[k, i] _, Subscript[\[Omega], i] _]] == 0,
Im[f[Subscript[k, i] _, Subscript[\[Omega], i] _]] ==
0}, {{Subscript[k, i], 1.5}, {Subscript[\[Omega], i], 160}}]
求解会报错,而且求不出结果,请问这是什么原因?该怎么解决?



2019年02月20日 08点02分
1
f[Subscript[k, i] _,
Subscript[\[Omega],
i] _] = ((\[Omega] - k*I)^2*BesselI[m, Subscript[n, 1]])/(
Subscript[n, 1]*
D[BesselI[m, Subscript[n, 1]], Subscript[n, 1]]) - (
Q*\[Omega]^2*BesselK[m, Subscript[n, 2]])/(
Subscript[n, 2]*D[BesselK[m, Subscript[n, 2]], Subscript[n, 2]]) +
We*(k^2 + m^2 - 1) /. {Subscript[n,
1] -> (k^2 + Subscript[Ma, 1]^2*(\[Omega] - k*I)^2)^0.5,
Subscript[n,
2] -> (k^2 + Subscript[Ma, 2]^2*\[Omega]^2 )^0.5} /. {\[Omega] ->
Subscript[\[Omega], r] + Subscript[\[Omega], i]*I,
k -> Subscript[k, r] + Subscript[k, i]*I}
m = 0;
We = 8*10^-6;
Q = 0.0014;
Subscript[Ma, 1] = 0.2;
Subscript[Ma, 2] = 0.57;
Subscript[\[Omega], r] = 0;
Subscript[k, r] = 150;
FindRoot[{Re[f[Subscript[k, i] _, Subscript[\[Omega], i] _]] == 0,
Im[f[Subscript[k, i] _, Subscript[\[Omega], i] _]] ==
0}, {{Subscript[k, i], 1.5}, {Subscript[\[Omega], i], 160}}]
求解会报错,而且求不出结果,请问这是什么原因?该怎么解决?


