ODEINT 求解常微分方程(3）

最新推荐文章于 2024-03-14 11:36:47 发布

weixin_30433075

最新推荐文章于 2024-03-14 11:36:47 发布

阅读量548

点赞数

CC 4.0 BY-SA版权

文章标签： python

原文链接：http://www.cnblogs.com/conpi/p/11102719.html

本文介绍了一个使用Python中的odeint函数求解微分方程组的例子。通过定义函数model来描述微分方程，初始条件z0设置为[0,0]，并利用numpy生成时间点t。通过调用odeint函数求解微分方程，并使用matplotlib绘制结果。此例展示了如何使用odeint解决实际问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt

# function that returns dz/dt
def model(z,t):
    dxdt = 3.0 * np.exp(-t)
    dydt = -z[1] + 3
    dzdt = [dxdt,dydt]
    return dzdt

# initial condition
z0 = [0,0]

# time points
t = np.linspace(0,5)

# solve ODE
z = odeint(model,z0,t)

# plot results
plt.plot(t,z[:,0],'b-',label=r'$\frac{dx}{dt}=3 \; \exp(-t)$')
plt.plot(t,z[:,1],'r--',label=r'$\frac{dy}{dt}=-y+3$')
plt.ylabel('response')
plt.xlabel('time')
plt.legend(loc='best')
plt.show()