数学模型__姜启源讲述.ppt

数 学 模 型 方程（14）无法求出和的解析解，我们先作数值计算。用Matlab编程如下： function y=ill(t,x) a=1;b=0.3; y=[a*x(1)*x(2)-b*x(1),-a*x(1)*x(2)]; ts=0:50; x0=[0.02,0.98]; [t,x]=ode45(@ill,ts,x0);[t,x] plot(t,x(:,1),t,x(:,2)); grid figure plot(x(:,2),x(:,1)),grid, 运行结果如下： 0 0.0200 0.9800 1.0000 0.0390 0.9525 2.0000 0.0732 0.9019 3.0000 0.1285 0.8169 4.0000 0.2033 0.6927 5.0000 0.2795 0.5438 6.0000 0.3312 0.3995 7.0000 0.3444 0.2839 8.0000 0.3247 0.2027 9.0000 0.2863 0.1493 10.0000 0.2418 0.1145 11.0000 0