python 微分方程组
时间: 2023-09-02 13:09:02 浏览: 100
Hopf振荡器的一种python实现
在 Python 中,可以使用 SciPy 模块中的 `odeint` 函数来求解微分方程组。下面是一个示例代码,解决了一个简单的微分方程组:
```python
import numpy as np
from scipy.integrate import odeint
# 定义微分方程组
def f(y, t):
y0, y1 = y
return [y1, -y0]
# 初始条件
y0 = [1.0, 0.0]
# 时间轴
t = np.linspace(0, 10, 101)
# 求解微分方程组
y = odeint(f, y0, t)
# 绘制图像
import matplotlib.pyplot as plt
plt.plot(t, y[:, 0], 'b', label='y0(t)')
plt.plot(t, y[:, 1], 'g', label='y1(t)')
plt.legend(loc='best')
plt.xlabel('t')
plt.grid()
plt.show()
```
此代码定义了一个简单的微分方程组 `y'' + y = 0`,并使用 `odeint` 函数求解该方程组。然后,结果绘制出来。
可以根据实际情况修改微分方程组的定义和初始条件。
阅读全文