|
|
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
这样的问题 最好还是找下师兄或者老师呀
第一个方程没有唯一解 看数值解情况 大概象是个T的周期函数，
下面是代码 求出T(t) 和k(t) 取了0<T<10的第一个解
另外 epsilon_2 > epsilon_1
画图速度比较慢 建议自己打点输出结果 然后用ListPlot
---------------------------------------------------
\[Epsilon]1 = 2;
\[Epsilon]2 = 1;
\[Mu]0 = 3;
tf[\[Tau]_?NumericQ] :=
t /. NSolve[{\[Epsilon]1*BesselJ[1, t]/(
t*BesselJ[0, t]) + \[Epsilon]2*BesselJ[1, \[Tau]]/(
t*BesselJ[0, \[Tau]]) == 0, 0 < t < 100}, t]
tf[1.0]
Differences[%]
ParametricPlot[{tf[\[Tau]][[1]],
Sqrt[(\[Tau]^2 +
tf[\[Tau]][[1]]^2)/(\[Mu]0*(\[Epsilon]1 - \[Epsilon]2))]}, {\
\[Tau], 0, 1}]
------------------------------------------------- |
|
回复此楼