from sympy import symbols,diff import sympy as sp from sympy.abc import beta,omega,u import math k1 = 150.0 k2 = 100.0 a = 0.6 b = 0.77 L = 1.37 step_size = 0.01 time = 10 m = 256 vx = 30.0 Iz = 0.1 vx_plot = [] omega_plot = [] vx = vx sigma = 0.5 eq1 = sp.Eq((k1+k2)*beta+omega/u*(a*k1-b*k2)-k1*sigma -m*(0.09+u*omega),0) eq2 = sp.Eq((a*k1-b*k2)*beta+omega/u*(a**2*k1-b**2*k2)-a*k1*sigma -Iz*omega,0) eq3 = sp.Eq((u-vx/math.tan(beta)),0) sol = sp.solve((eq1,eq2,eq3),(beta,omega,u)) print(sol)
时间: 2024-04-06 22:28:48 浏览: 100
sympy-1.5.1-py2.py3-none-any.whl
您的代码看起来是在使用SymPy库对一个非线性方程组进行求解。您定义了三个方程,分别是eq1、eq2和eq3,然后使用solve()函数来求解这三个方程组成的方程组,这是正确的。最终的解存储在sol中,并且可以打印出来。
不过,我注意到您没有定义任何初始值,这可能会影响求解结果。您可以使用symbols()函数来定义变量,使用diff()函数来定义微分方程,使用subs()函数来代入初始值,从而获得更准确的结果。
阅读全文