24小时热门版块排行榜    

查看: 402  |  回复: 1

xwndf250

银虫 (小有名气)

[求助] 用matlab求解常微分方程,但无法运行,求助下,另追加10金币

原方程:
S+I+R1+Mi+Ml=1
dS/dt=γβI^2+γ(-2α+β)SI+γβ(M_i+M_l )I+k_is M_l
dI/dt=-γβI^2-γ(-2α+β)SI-γβ(M_i+M_l )I-(γ-k_a)I
dR/dt=γI+k_ir M_i
(dM_i)/dt=k_a I-k_il M_i-k_ir M_i
(dM_l)/dt=k_il M_i-k_ls M_l
求S、I、R、Mi、Ml与时间t的变化(先求的S和t的变化);
用ode45解,0《t《100;当t=0时,S=5/6、I=1/6;R、Mi、Ml=0;
拟合参数为:α=0.1、β=0.5、γ=0.1、r=0.1、ka=0.5、kb=0.5、kir=0.01、kil=0.01、kls=0.1。
本人小白,请高手指教

SIRS1.jpg

[ Last edited by xwndf250 on 2013-2-14 at 12:59 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwndf250

银虫 (小有名气)

已解决.
2楼2013-02-14 16:02:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xwndf250 的主题更新
信息提示
请填处理意见