数学建模课件_第二部分_动态模型.pptVIP

区域内部有平衡点： x10≈9333,x20≈1333 计算偏导数矩阵在平衡点(x10,x20)处的值A。 求矩阵A的特征值，得 这两个特征值均具有负实部，所以该平衡点为稳定的。 上述过程也可以用计算机求解。 syms x1 x2 f1 = 0.1*x1*(1-(1/10000)*x1- 0.5*(1/10000)*x2); f2 = 0.25*x2*(1-(1/6000)*x2- 0.5*(1/6000)*x1); [x1steady,x2steady] = solve(f1,f2); N=length(x1steady); fprintf(The equilibrium points are\n) disp([x1steady x2steady]) M=15; x1min=0; x1max=11000; x2min=0; x2max=8000; ezplot(f1,[x1min x1max x2min x2max]),hold on ezplot(f2,[x1min x1max x2min x2max]),hold on [X1,X2]=meshgrid(x1min:(x1max-x1min)/M:x1max,x2min:(x2max-x2min)/M:x2max); dX1 = 0.1*X1.*(1-(1/10000)*X1-0.5*(1/