python - 在 python 上求解非线性方程

Question

我一直在尝试解决机械系统的非线性方程，但我找不到正确的方法。

import numpy as np
from sympy import *
import scipy.optimize as opt

a = 0.15
b = 0.4
c = 0.5
d = 0.4
e = 0.3
f = 0.35
beta = 0.3
qDot = 30
q = 0


def eqsSympy():
    teta3 = Symbol('teta3',real=True)
    teta4 = Symbol('teta4',real=True)
    teta5 = Symbol('teta5',real=True)
    s = Symbol('s',real=True)

    e1 = Eq(a*cos(q)+b*cos(teta3)-c*cos(teta4)-d)
    e2 = Eq(a*sin(q)+b*sin(teta3)-c*sin(teta4))
    e3 = Eq(e*cos(teta4-beta)+f*cos(teta5))
    e4 = Eq(e*sin(teta4-beta)+f*sin(teta5)-s)

    print(solve([e1,e2,e3,e4],teta3,teta4,teta5,s))
eqsSympy()

所以我希望用 pi/6 每步递增 0 到 2*pi 的 q 值来求解这些方程。但我不断收到错误：

e1 = Eq(a*cos(q)+b*cos(teta3)-c*cos(teta4)-d)
TypeError: unsupported operand type(s) for -: 'Add' and 'function'

如果有更好的方法来解决这些方程，如果你能帮助我，我将不胜感激。

标签： pythonmathsympynonlinear-optimization