一个游击战问题模型论文2（精选）

ction y=ill(t,x)Рa=5.3*10^(-3);Рb=4.6*10^(-6);Рy=[-a*x(2),-b*x(1)*x(2)]';Р② function y=ill1(t,x)Рa=5.61*10^(-3);Рb=4.34*10^(-6);Рy=[-a*x(2),-b*x(1)*x(2)]';Р③ function y=ill2(t,x)Рa=7.48*10^(-3);Рb=5.13*10^(-6);Рy=[-a*x(2),-b*x(1)*x(2)]';Р④ function y=ill3(t,x)Рa=7.48*10^(-3);Рb=5.37*10^(-6);Рy=[-a*x(2),-b*x(1)*x(2)]';Р⑤ function y=ill4(t,x)Рa=3.74*10^(-3);Рb=4.06*10^(-6);Рy=[-a*x(2),-b*x(1)*x(2)]';Р主程序Рplot(x(:,2),x(:,1),'r') Рhold onР[t,x]=ode45('ill1',ts,x0);Рplot(x(:,2),x(:,1),'b') Рhold onР[t,x]=ode45('ill2',ts,x0);Рplot(x(:,2),x(:,1),'g')Рhold onР[t,x]=ode45('ill3',ts,x0);Рplot(x(:,2),x(:,1),'k') Рhold onР[t,x]=ode45('ill4',ts,x0);Рplot(x(:,2),x(:,1),'b') Рtitle('双方人数随时间递减')Рxlabel('游击队人数变化')Рylabel('日军人数变化')Р参考文献Р[1]姜启源等编。数学模型.高等教育出版社.2011.1:149-154Р[2]张德丰等编。matlab数值分析.机械工业出版社.2012.3