[数学]数学建模作业2.doc

数学建模作业-----第二章 ①微分方程稳定性分析 由解得平衡点, , 特征值 ∵, , , ∴是不稳定点 Matlab作图程序: x=[-5:0.1:5]; fora=-5:0.1:5; y=a*x; plot(x,y);hold on; end xlabel(x);ylabel(y); 得到图形： （2）由解得平衡点, , 特征值 ∵, ∴是不稳定点 Matlab作图程序： y=[0:0.1:10]; for a=-5:0.5:5; x=a./(y.^(0.5)); plot(y,x);hold on; end xlabel(y);ylabel(x); 得到如下图形： （3）由解得平衡点, , 特征值 ∵, ∴是不稳定点 用matlab作图程序如下： x=[0:0.01:5]; for a=-10:0.1:10; y=((-2)*x.^2+a).^0.5; plot(x,y);hold on; end xlabel(x);ylabel(y); 得到如下图形： （4）由解得平衡点, , 特征值 ∵,,, ∴是稳定点 用matlab程序作图如下： x=[-5:0.01:5]; for a=-5:0.5:5; y=a*(x.^2); plot(x,y);hold on; end xlabel(x);ylabel(y); 得到如下图形： ②种族增长模型： 由题意，得到菌落的细菌数目N(t)变化的模型为： dN(t)dt=r1N(t)-r2N(t) 设g(t)=r1N(t)-r2N(t), 则令g(t)=0,得到 N=0或N=r22r12 所以此模型的平衡点为N=0或r22r12，并且可以得到： 当0Nr22r12时，N’(t)0; 当Nr22r12时，N’(t)0。 设g’(t)=d2Ndt2=(r1-12Nr2)(r1N-r2N) 令g’(t)=0，得： N=0或r22r12或r224r12. 分析得到： 当0Nr224r12时，N’’(t)0;当r224r12Nr22r12时，N”(t)0; 当Nr22r12时，N”(t)0; 画出趋势图形为： N t 分析：是稳定的，是不稳定的，即表示当种群数量小于r22r12时，病菌会随时间变化逐渐减少0，趋于稳定；当初时病菌数量恰好为r22r12时，数量并不会随时间改变；当初时病菌数量大于r22r12时，则会随时间逐渐增加； ③有限资源竞争模型 （1） 由得到： x1=x2=0或x1=0,x2=b2λ2h2或x1=b1λ1h1，x2=0 三个平衡点： p0(0,0)、p1（0，b2λ2h2）、p2（b1λ1h1，0） 对于，系数矩阵， ，所以不稳定。 2．对于，系数矩阵，，。 ∵ ∴，是稳定的 ∴当时，当时， 3.对于，系数矩阵，， ∵ ∴，不稳定 所以，当时，物种二最终要灭亡。 对于线性方程组在平面上表示两条直线和。 （2） 对于线性方程组在平面上表示两条直线和。 当和，直线和将第一象限分成三个区域Ⅰ，Ⅱ，Ⅲ。 在区域I中，x2’(t)0,随着时间t物种2数量增加，到达直线l2时，有dx2dt=0,即此刻x2趋常数，说明曲线平行于横轴过度到区域II。 在区域II中，物种2数量减少，物种1增加。当经过直线l1时，dx1dt=0，即此刻x 1为一常数，说明曲线是垂直于横轴向下过度到区域III的。 在区域III中，dx2dt0说明物种2随时间变化数量减少，最终消失。综上分析：物种二最终要灭亡。 图如下： ④蝴蝶效应和混沌解 (1)建立描述微分方程组的外部函数（函数名：lorenz.m） functionxdot=lorenz(t,x) xdot=[-8/3*x(1)+x(2)*x(3); -10*x(2)+10*x(3); -x(1)*x(2)+28*x(2)-x(3)]; 输入 [t,x]=ode45(@lorenz,[0 100],[0 0 1e-10]); plot(x(:,2), x(:,1)); xlabel(x2); ylabel(x1); 再输入[t,x]=ode45(@lorenz,[0 100],[0 0 1e-10]); plot(x(:,2), x(:,3)); xlabel(x2); ylabel(x3); 再输入 [t,x]=ode45(@lorenz,[0 100],[0 0 1e-10]); plot(x(:,3), x(:,1)); xlabel(x3); ylabel(x1); 再输入 [t,x]=ode45(@lorenz,[0 100],[0 0 1e-10]); plot(x(:,3), x(:,1)); xlabel(x3); ylabe