performance - Mathematica:代码在几天前提供了输出,但现在一直处于“运行”状态
问题描述
我求解了一个方程,其中输出为感兴趣变量的 5 个解。花了一段时间,但它昨天跑了。我今天尝试使用ClearSystemCache[]
和ClearGlobal
选项重新运行代码,但它似乎在过去 2 小时内“正在运行”。我中止并再次尝试,重新启动计算机,清除计算机缓存并再次尝试......没有运气。当我使用数值求解时,NDSolve
我可以看到它生成了 5 个解决方案,就像前一个一样。但是,我希望使用解析解,而不是数值解。什么可能导致这种窒息?
Eq1 = (1/(e1^2 ep^2))(-(((a1 ep f f1 fp R - e1 ep fp^2 m1 (h + R) +
ep^2 f1 (h + R) (ab fp - a1 mp R) + Sqrt[
ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)]) (ep^2 f1 (h +
R) (ab fp + a1 mp R) +
ep fp (a1 f f1 R -
e1 m1 (fp - 2 a1 mp) (h +
R)) + \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)))^4 (-ep^2 f1 (h +
R) (ab f f1 fp R - a1 f f1 mp R^2 +
ab e1 fp m1 mp (h + R) +
e1 m1 mp (2 fp - a1 mp) R (h + R)) -
ep fp (a1 f^2 f1^2 R^2 - e1 f f1 m1 (fp - a1 mp) R (h + R) +
e1^2 fp m1^2 mp (h + R)^2) +
e1 h m1 mp \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)) +
f f1 R \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)) +
e1 m1 mp R \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)))^4 (fp + (1/(
2 ep (h + R) (e1 fp m1 + ep f1 R)))(a1 ep f f1 fp R -
e1 ep fp^2 m1 (h + R) +
ep^2 f1 (h + R) (ab fp -
a1 mp R) + \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp +
ep mp (h +
R)))^2)))))/(512 a1^2 ep^9 f1^2 fp^6 mp^4 (h +
R)^9 (e1 fp m1 + ep f1 R)^5 (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R))^4)) -
e1^2 ep m1 (ep mp + (1/(
2 a1 (h + R) (e1 fp m1 + ep f1 R)))(a1 ep f f1 fp R -
e1 ep fp^2 m1 (h + R) + ep^2 f1 (h + R) (ab fp - a1 mp R) +
Sqrt[ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)]) - (1/(
2 a1 (h + R) (e1 fp m1 +
ep f1 R)))(ep^2 f1 (h + R) (ab fp + a1 mp R) +
ep fp (a1 f f1 R - e1 m1 (fp - 2 a1 mp) (h + R)) + Sqrt[
ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)])) (R - (-ep^2 f1 (h +
R) (ab f f1 fp R - a1 f f1 mp R^2 +
ab e1 fp m1 mp (h + R) +
e1 m1 mp (2 fp - a1 mp) R (h + R)) -
ep fp (a1 f^2 f1^2 R^2 - e1 f f1 m1 (fp - a1 mp) R (h + R) +
e1^2 fp m1^2 mp (h + R)^2) +
e1 h m1 mp Sqrt[
ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)] +
f f1 R Sqrt[
ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)] +
e1 m1 mp R Sqrt[
ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)])/(2 ep^2 f1 mp (h +
R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R))) - ((a1 ep f f1 fp R -
e1 ep fp^2 m1 (h + R) + ep^2 f1 (h + R) (ab fp - a1 mp R) +
Sqrt[ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)]) (-ep^2 f1 (h +
R) (ab f f1 fp R - a1 f f1 mp R^2 +
ab e1 fp m1 mp (h + R) +
e1 m1 mp (2 fp - a1 mp) R (h + R)) -
ep fp (a1 f^2 f1^2 R^2 - e1 f f1 m1 (fp - a1 mp) R (h + R) +
e1^2 fp m1^2 mp (h + R)^2) +
e1 h m1 mp \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)) +
f f1 R \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)) +
e1 m1 mp R \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp +
ep mp (h + R)))^2))))/(4 ep^3 f1 fp mp (h +
R)^2 (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)))) +
ab e1 ep (-(((ep^2 f1 (h + R) (ab fp + a1 mp R) +
ep fp (a1 f f1 R -
e1 m1 (fp - 2 a1 mp) (h +
R)) + \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2))) (-ep^2 f1 (h +
R) (ab f f1 fp R - a1 f f1 mp R^2 +
ab e1 fp m1 mp (h + R) +
e1 m1 mp (2 fp - a1 mp) R (h + R)) -
ep fp (a1 f^2 f1^2 R^2 -
e1 f f1 m1 (fp - a1 mp) R (h + R) +
e1^2 fp m1^2 mp (h + R)^2) +
e1 h m1 mp \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)) +
f f1 R \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)) +
e1 m1 mp R \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp +
ep mp (h + R)))^2))))/(4 a1 ep fp mp (h +
R)^2 (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)))) +
e1 m1 ((1/(
2 a1 (h + R) (e1 fp m1 +
ep f1 R)))(ep^2 f1 (h + R) (ab fp + a1 mp R) +
ep fp (a1 f f1 R - e1 m1 (fp - 2 a1 mp) (h + R)) + Sqrt[
ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)]) -
ep (mp + (1/(
2 a1 ep (h + R) (e1 fp m1 + ep f1 R)))(a1 ep f f1 fp R -
e1 ep fp^2 m1 (h + R) +
ep^2 f1 (h + R) (ab fp -
a1 mp R) + \[Sqrt](ep^2 (-4 a1 fp mp (h +
R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp +
ep mp (h + R)))^2)))))) - (e1 (ep^2 f1 (h +
R) (ab fp + a1 mp R) +
ep fp (a1 f f1 R - e1 m1 (fp - 2 a1 mp) (h + R)) + Sqrt[
ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)]) (-ep^2 f1 (h +
R) (ab f f1 fp R - a1 f f1 mp R^2 + ab e1 fp m1 mp (h + R) +
e1 m1 mp (2 fp - a1 mp) R (h + R)) -
ep fp (a1 f^2 f1^2 R^2 - e1 f f1 m1 (fp - a1 mp) R (h + R) +
e1^2 fp m1^2 mp (h + R)^2) +
e1 h m1 mp Sqrt[
ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)] +
f f1 R Sqrt[
ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)] +
e1 m1 mp R Sqrt[
ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h + R) -
e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)]) (ep R (ep mp - (1/(
2 a1 (h + R) (e1 fp m1 +
ep f1 R)))(ep^2 f1 (h + R) (ab fp + a1 mp R) +
ep fp (a1 f f1 R -
e1 m1 (fp - 2 a1 mp) (h +
R)) + \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp +
ep mp (h + R)))^2)))) - ((-ep^2 f1 (h +
R) (ab f f1 fp R - a1 f f1 mp R^2 +
ab e1 fp m1 mp (h + R) +
e1 m1 mp (2 fp - a1 mp) R (h + R)) -
ep fp (a1 f^2 f1^2 R^2 -
e1 f f1 m1 (fp - a1 mp) R (h + R) +
e1^2 fp m1^2 mp (h + R)^2) +
e1 h m1 mp \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)) +
f f1 R \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)) +
e1 m1 mp R \[Sqrt](ep^2 (-4 a1 fp mp (h + R) (e1 fp m1 +
ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp +
ep mp (h + R)))^2))) (-((m1 (a1 ep f f1 fp R -
e1 ep fp^2 m1 (h + R) +
ep^2 f1 (h + R) (ab fp -
a1 mp R) + \[Sqrt](ep^2 (-4 a1 fp mp (h +
R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2))) (ep^2 f1 (h +
R) (ab fp + a1 mp R) +
ep fp (a1 f f1 R -
e1 m1 (fp - 2 a1 mp) (h +
R)) + \[Sqrt](ep^2 (-4 a1 fp mp (h +
R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp +
ep mp (h + R)))^2)))^2 (-ep^2 f1 (h +
R) (ab f f1 fp R - a1 f f1 mp R^2 +
ab e1 fp m1 mp (h + R) +
e1 m1 mp (2 fp - a1 mp) R (h + R)) -
ep fp (a1 f^2 f1^2 R^2 -
e1 f f1 m1 (fp - a1 mp) R (h + R) +
e1^2 fp m1^2 mp (h + R)^2) +
e1 h m1 mp \[Sqrt](ep^2 (-4 a1 fp mp (h +
R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)) +
f f1 R \[Sqrt](ep^2 (-4 a1 fp mp (h +
R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2)) +
e1 m1 mp R \[Sqrt](ep^2 (-4 a1 fp mp (h +
R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2))) (fp + (1/(
2 ep (h + R) (e1 fp m1 +
ep f1 R)))(a1 ep f f1 fp R -
e1 ep fp^2 m1 (h + R) +
ep^2 f1 (h + R) (ab fp -
a1 mp R) + \[Sqrt](ep^2 (-4 a1 fp mp (h +
R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp +
ep mp (h +
R)))^2)))))/(16 a1^2 ep^3 f1 fp^2 mp (h +
R)^4 (e1 fp m1 + ep f1 R)^3 (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)))) +
ep (ep mp + (1/(
2 a1 (h + R) (e1 fp m1 + ep f1 R)))(a1 ep f f1 fp R -
e1 ep fp^2 m1 (h + R) +
ep^2 f1 (h + R) (ab fp -
a1 mp R) + \[Sqrt](ep^2 (-4 a1 fp mp (h +
R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp + ep mp (h + R)))^2))) - (1/(
2 a1 (h + R) (e1 fp m1 +
ep f1 R)))(ep^2 f1 (h + R) (ab fp + a1 mp R) +
ep fp (a1 f f1 R -
e1 m1 (fp - 2 a1 mp) (h +
R)) + \[Sqrt](ep^2 (-4 a1 fp mp (h +
R) (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)) + (ab ep f1 fp (h +
R) - e1 fp m1 (fp - 2 a1 mp) (h + R) +
a1 f1 R (f fp +
ep mp (h + R)))^2))))))/(2 ep^2 f1 mp (h +
R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R)))))/(4 ep^2 fp^2 mp (h +
R)^2 (e1 fp m1 + ep f1 R) (a1 f f1 R -
e1 m1 (fp - a1 mp) (h + R))))
这一直在持续运行。以前我得到了 5 个解决方案的输出。
Eq1solve = Solve[Eq1==0, fp]
当我使用数值求解器时,这一点得到了验证:
Pars = {f -> 1, h -> 1, ab -> 0.1, a1 -> 0.1, f1 -> 0.1, e1 -> 0.1,
m1 -> 0.1, ep -> 0.1, mp -> 0.01,R -> 1000};
NSolve[{F3 == 0} /. Pars, fp]
Plot[F3 /. Pars, {fp, -1, 5}]
#Output
{{fp -> -0.193087}, {fp -> 0.0123643 - 0.17559 I}, {fp ->
0.0123643 + 0.17559 I}, {fp -> 0.167316}, {fp -> 1.1011}}
解决方案
推荐阅读
- python - 如何捆绑 Python 包(案例研究:Microsoft Visual Studio 2017)
- r - 如何对数据框进行分组,然后汇总每个组的总数,然后除以该组中的项目数?
- python - python字典检查键
- python - 拆分单元格,创建列,计数
- python - Web使用Python抓取显示为1970年1月而不是真实值的日期值
- sql - 如何为加入的“列表”表实现 EVERY 和 NO 操作
- rabbitmq - 在 SignalR 应用程序中订阅 RabbitMq 事件的位置
- r - 将每个观察结果相互匹配
- mysql - 给定一个酒吧的定价,获得一个人的消费总和
- ruby-on-rails - 允许更新但阻止更新中允许列?