数学建模作业北工大薛毅实验2.docVIP

第二次作业 平衡点0,0)， 列出A={1,0}， {0,1} 所以特征值为λ1=1，λ2=1， 知p=-(λ1+λ2)0，q=λ1λ20，所以P(0，0 定的 （2）得到平衡点P(0, 0)， A={-1,0} {0,2} 所以其特征值λ1=-1，λ2=2； p=-(λ1+λ2)0，q=-20；所以P(0, 0)是不稳定的 （3）平衡点P(0, 0) A={0,1} {-2,0} 所以特征值λ1= i，λ2=- i； p=-(λ1+λ2)=0，q=λ1λ2=；所以平衡点(0, 0)是不稳定的。 （4）平衡点P(0, 0)， A={-1，0} {0，-2}所以其特征值λ1=-1，λ2=-2； p=-(λ1+λ2) 0，q=λ1λ20；易知平衡点(0, 0)是稳定的。 解：(1) 得 (2) 得N=R/K. 方程为 当t无穷大， R/K. 画出N（t）的图形，如下图所示： 解： (1) 所以N1=0, N2=； 由(1)得： (2) 所以N= 画出N（t）的图形，如下图所示： 其中N1=0是不稳定的；N2=是稳定的。 解：设F(x)= = 0 (1得到：， （2） 带入x0，x1 得 （3） 所以当Er时， x0是稳定平衡点，x1不是； 所以当捕捞适度时,可使渔场产量稳定在， 可以持续获得产量此时考虑最大的产量的条件。 W(x)=Ex（4）x)为最大的产量）将带入（4）式 得到x)= 再由基本不等式可以知道当r/2时 W(x)可以取到最大值 综上：最优捕捞率为 。 解：（1）令 得方程的平衡点为P0（0,0），P1（,0），P2（0, ）. 对平衡点P0（0,0），假设c1a1，c2a2 又c1a1，c2a2则p0，所以该平衡点不稳定。 对平衡点P1（,0）： 则p，q0，P1对于P2（0, ），有p0，且q0稳定，此时， 综上：只有P2点稳定，说明物种1最终要灭亡。而如果的情况下，设c1a1，c2a2， 对平衡点P0（0,0）， 系数矩阵 p=-[(c1-a1)+(c2-a2)] 0，0所以该平衡点不稳定。对于P1（,0）， 则p0，q0，P1 对于P2（0, ）p，q0,所以不稳定。 综上:除了P2点其他点不稳定，此时说明物种2最终会灭亡。 （3） 物种1最终要灭亡 物种2最终会灭亡 解： 1 .建立“Lorenz.m”的M文件function dy = Lorenz(t,y) dy=zeros(3,1); dy(1)=-8*y(1)/3+y(2) *y(3); dy(2)= -10*y(2)+10*y(3); dy(3)=-y(1)*y(2) +28*y(2)-y(3); end 2.建立“Lzdis.m” [t,y]=ode45(Lorenz,[0,30],[ 0 0 1e-10]); axis([10 40 -20 20]); plot(x(:,1),x(:,2)); plot(x(:,1),x(:,3));axis([10 40 -20 20]); plot(x(:,3),x(:,2));axis([-20 20 -20 20]); 对各参数赋值并用ode45函数求解，可得数值解： Columns 1 through 9 0 0.1250 0.2500 0.3750 0.5000 0.5352 0.5705 0.6057 0.6409 0 0.0000 0.0000 -0.0000 -0.0000 0.0000 0.0000 0.0000 0.0000 0 0.0000 0.0000 -0.0000 -0.0000 0.0000 0.0000 0.0000 -0.0000 0.0000 -0.0000 -0.0000 0.0000 0.0000 0.0000 -0.0000 0.0000 0.0000 Columns 10 through 18 0.6761 0.7114 0.7466 0.7818 0.8308 0.8797 0.9286 0.9776 1.0105 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000