1=diff(g,'x1');gx2=diff(g,'x2');Рdisp('');A=[fx1,fx2;gx1,gx2]%显示结果Рp=subs(-(fx1+gx2),{x1,x2},{P(:,1),P(:,2)));%替换Рp=simple(p)%简化符号表达式pРq=subs(det(A),{x1,x2},{P(:,1),P(:,2)});Рq=simple(q);Рdisp('');[Ppq]%显示结果Р%得到教材p225表1的前3列，经测算可得该表的第4列，即稳定条件РРРР2.2(验证、编程)计算与验证p227Р微分方程组Р,、“x1x2、РX1⑴=r1x1(1--f-)РjNiN2Рx1x2РX2(t)=搂2(1-02)РNiN2Р、=0.5,二2二1.6,r1=2.5,r2=1.8,N1=1.6,N2=1Р(1)(验证)当x1(0)=x2(0)=0.1时，求微分方程的数值解，将解的数值分别画Р出x1(t)和x2(t)的曲线，它们同在一个图形窗口中。Р程序：Р%命令文件Рts=0:0.2:8;Рx0=[0.1,0.1];Р[t,x]=ode45cp227',ts,x0);Рplot(t,x(:,1),t,x(:,2));%x(:1)是x1(t)的一组函数值,x(:,2)对应x2(t)Рgridon;Рaxis([0802]);Рtext(2.4,1.55,'x1(t)');Рtext(2.2,0.25,'x2(t)');РР%函数文件名：p227.mРfunctiony=p227(t,x)Рk1=0.5;k2=1.6;r1=2.5;r2=1.8;N1=1.6;N2=1;Рy=[r1*x(1)*(1-x(1)/N1-k1*x(2)/N2),r2*x(2)*(1-k2*x(1)/N1-x(2)/N2)]';Р☆(1)运行程序并给出结果(比较［227］图2(a)):
全文预览
