1

我有一个方程，我试图NSolve在Mathematica中使用它来求解。这是我的做法：

T == (0.000242895 E^(-(2.09472*10^11/(5.70068*10^8 + 
            7.76206*10^12 T))))/(1 + 
    0.969073 E^(-(4.18945*10^11/(5.70068*10^8 + 7.76206*10^12 T))) - 
    1.96883 Cos[8.77331*10^6/(2.28027*10^9 + 7.76206*10^12 T)])

NSolve[T == (0.000242895 E^(-(2.09472*10^11/(5.70068*10^8 + 
             7.76206*10^12 T))))/(1 + 
     0.969073 E^(-(4.18945*10^11/(5.70068*10^8 + 7.76206*10^12 T))) - 
     1.96883 Cos[8.77331*10^6/(2.28027*10^9 + 7.76206*10^12 T)])
 , {T}, Reals]

问题是需要很长时间（约 1 分钟）才能找到解决方案。有没有比使用更快的方法来获得解决方案NSolve？

4

1 回答 1

2

您的表达式的分母在大约 T=215 处通过零。另一个解是 T=0。如果您想要非零解决方案，只需使用

FindRoot[denom,{T,T0}]

其中 T0 几乎是任何正数。这大约需要一毫秒。

于 2013-05-01T15:08:53.653 回答