mathematica吧 关注：19,994贴子：74,342
- 6回复贴，共1页
不感兴趣
开通SVIP免广告
1. 别用For的问题说得够多的了，这里不多谈了。
2. 计算会卡在半途毫不意外。NSolve的帮助的更多信息里已经说了，NSolve 处理的主要是线性和多项式方程。它在一定程度上是可以解超越方程的，但是限制颇多。所以，如果不是已经确定超越方程存在通解，那就不要一骨脑地用Solve/NSolve去扫一大片参数区域。实际上稍加实验即可发现Solve和NSolve在求解这个方程时都遇到了困难。
3. 对于超越方程的数值求根，FindRoot常常是一个更好的选择。（如果根只有一个的话。）
4. 对于这类问题，先画图观察一下是有好处的：
e1 = 20 - 250 x - ((8 g^2 - 30 (0.0047 - 30 x)) Tanh[
13.63969351909815 Sqrt[(0.0047 - 30 x)^2 + 16 g^2 x]])/(
2 Sqrt[(0.0047 - 30 x)^2 + 16 g^2 x]);
ContourPlot[e1 == 0, {g, 0.01, 1.2}, {x, -0.03, 0.03}, PlotPoints -> 200]
在g较大时方程有多组解。如果这些解你全需要，那么你还需要使用RootSearch程序包：http://library.wolfram.com/infocenter/Demos/4482/
2. 计算会卡在半途毫不意外。NSolve的帮助的更多信息里已经说了，NSolve 处理的主要是线性和多项式方程。它在一定程度上是可以解超越方程的，但是限制颇多。所以，如果不是已经确定超越方程存在通解，那就不要一骨脑地用Solve/NSolve去扫一大片参数区域。实际上稍加实验即可发现Solve和NSolve在求解这个方程时都遇到了困难。
3. 对于超越方程的数值求根，FindRoot常常是一个更好的选择。（如果根只有一个的话。）
4. 对于这类问题，先画图观察一下是有好处的：
e1 = 20 - 250 x - ((8 g^2 - 30 (0.0047 - 30 x)) Tanh[
13.63969351909815 Sqrt[(0.0047 - 30 x)^2 + 16 g^2 x]])/(
2 Sqrt[(0.0047 - 30 x)^2 + 16 g^2 x]);
ContourPlot[e1 == 0, {g, 0.01, 1.2}, {x, -0.03, 0.03}, PlotPoints -> 200]
在g较大时方程有多组解。如果这些解你全需要，那么你还需要使用RootSearch程序包：http://library.wolfram.com/infocenter/Demos/4482/
