import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt
def model(y, u ,t):
2dydt=-y+u
return dydt
y0=1
t=np.linspace(0,10)
u=0
y1=odeint(model, y0, t, args=(u, ))
u=0.2
y2=odeint(model, y0, t, args=(u, ))
.
.
.
u=2.0
y11=odeint(model, y0, t, args=(u, ))
plt.plot(t, y1, 'g:', linewidth=5, label='u=0')
plt.plot(t, y2, 'g:', linewidth=5, label='u=0.2')
.
.
.
plt.plot(t, y11, 'g:', linewidth=5, label='u=2.0')
plt. xlabel('time')
plt. ylabel('y(t)')
plt.show()
不好意思，自学程式中的我，只能贴上来求指正
如果有错的地方，希望可以指点我一下
...的部分代表条件递移的相同程式码，直接省略掉了