scipy.integrate import odeint
#用来正常显示中文标签
plt.rcParams['font.sans-serif']=['SimHei']
def dySIR(y,t,lamda,mu):# SIR...i
return np.array([di_dt,ds_dt])
#设置模型参数
number=1e6#总人数
lamda =0.2#日接触率, 患病者每天有效接触的易感者的平均人数
sigma...i0=1-s0#i0,患病者比例的初值
Y0=(i0,s0)#微分方程组的初值
ySIR=odeint(dySIR,Y0,t,args=(lamda,mu))#SIR模型...plt.plot(ySIR[:,1],ySIR[:,0])
#绘图
plt.title("SIR模型的相轨迹")
plt.axis([0,1,0,1])
plt.plot([0,1],[1,0],'#ffc0cb...模型的相轨迹是研究解的周期性和稳定性。