9

我有一个与轨道力学中的太阳同步共振条件有关的方程。我目前正在学习 Python，所以我尝试使用以下代码在 SymPy 中解决它：

from sympy import symbols,solve

[n_,Re_,p_,i_,J2_,Pe_] = symbols(['n_','Re_','p_','i_','J2_','Pe_'])

del_ss = -((3*n_*(Re_**2)*J2_/(4*(p_**2)))*(4-5*(sin(i_)**2)))-((3*n_*(Re_**2)*J2_/(2*(p_**2)))*cos(i_))-((2*pi)/Pe_)

pprint(solve(del_ss,i_))

对于其中的五个变量，表达式可以成功重新排列，但是当i_在命令中使用变量时solve（如上），会产生错误：

Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "C:\Python27\lib\site-packages\spyderlib\widgets\externalshell\sitecustomize.py", line 479, in runfile
    execfile(filename, namespace)
  File "C:\Users\Nathan\Python\sympy_test_1.py", line 22, in <module>
    pprint(solve(del_ss,i_))
  File "C:\Python27\lib\site-packages\sympy\solvers\solvers.py", line 484, in solve
    solution = _solve(f, *symbols, **flags)
  File "C:\Python27\lib\site-packages\sympy\solvers\solvers.py", line 700, in _solve
    soln = tsolve(f_num, symbol)
  File "C:\Python27\lib\site-packages\sympy\solvers\solvers.py", line 1143, in tsolve
    "(tsolve: at least one Function expected at this point")
NotImplementedError: Unable to solve the equation(tsolve: at least one Function expected at this point

但是，当在 Matlab 中输入相同的表达式并solve调用命令时，它会正确重新排列。我意识到该错误提到了一个未实现的功能，并且这两个功能无疑会有所不同，但如果有更合适的 SymPy 功能可以使用，我仍然会很高兴。任何帮助将不胜感激。

4

1 回答 1

9
  1. 使用 Pi 的 sympy 版本。
  2. cos(i_)用新变量替换，用ci_替换，然后求解。sin(i_)**21-ci_**2ci_

这应该这样做：

from sympy import symbols,solve,sin,cos,pi

[n_,Re_,p_,ci_,J2_,Pe_] = symbols(['n_','Re_','p_','ci_','J2_','Pe_'])

del_ss = -((3*n_*(Re_**2)*J2_/(4*(p_**2)))*(4-5*(1-ci_**2)))-((3*n_*(Re_**2)*J2_/(2*(p_**2)))*ci_)-((2*pi)/Pe_)

pprint(solve(del_ss,ci_))

（编辑是因为我在第一次尝试时只写了一半的解决方案......）

于 2012-11-13T12:54:31.783 回答