数学建模实验2.docx

数学建模实验二2微分方程实验2.2 基本实验1.微分方程稳定性分析对于一般的二元微分方程组 (1)称代数方程组 (2)的实根x1=x10,x2=x20为微分方程(1)的平衡点，记作P0(x10,x20)。对于一般方程，可以考虑函数在平衡点(x10,x20)处展开，因此微分方程与线性方程 (3)有相同的平衡点。令 (4)对于(1)平衡点为(0,0)p= -tr(A)=-2;q=det(A)=1;由于p0,则(0,0)是不稳定的。对于(2)平衡点为(0,0)p=-tr(A)=-1;q=det(A)=2;由于p0,则(0,0)是不稳定的。对于(3)平衡点为(0,0)p=-tr(A)=0;q=det(A)=2;由于p=0,则(0,0)是不稳定的。x与y的关系图对于(4)平衡点为(0,0)p=-tr(A)=3;q=det(A)=2;由于p0,q0,则(0,0)是稳定的。2.种群增长模型得到N满足的微分方程为 (5)令得到则方程有平衡解令，则平衡解N1是稳定的。,则平衡解N2是不稳定的。3.单种群开发模型(1)令得到两个平衡点则x0是稳定的，x1是不稳定的；当画出的线族图如下:(2)令， (5)则平衡点是y=g(x)与直线y=Ex交点的横坐标从g(x)的数学表达式中可以看出，其定点坐标为(N/2,hm)。hm=Nr/4,当hhm时，无平衡点，并优于横有dx/dt0,从而x(t)下降，最后导致无鱼可捕，因而最大的捕捞量为hm=Nr/4，令其与(5)式中的h相等。得到稳定状态下最优捕捞率为。4. Compertz模型单种群开发方程如下：令得到两个平衡点则x0是不稳定的，x1是稳定的。令当时，鱼量x必0。则是g(x)的极大值对应的横坐标则最大持续产量此时则捕捞强度其中g(x)曲线如下：5.有限资源竞争模型 (6)知平衡点为1)对于平衡点P0所以P0点是不稳定的2)对于平衡点P1若，则p0,q0，则P1是稳定的。3)对于平衡点P2若，则p0,q0，则P2是稳定的。综上所述若，对于平衡点P2是稳定的，当若，对于平衡点P1是稳定的，当对于线性方程组在平面上表示两条直线l1和l21)当，对应的曲线如下：直线l1和l2将平面分成三个区域：在I区域：在II区域：在III区域：从图中看出当2)当，对应的曲线如下：直线l1和l2将平面分成三个区域：在I区域：在II区域：在III区域：从图中看出当6.蝴蝶效应与混沌解(1)建立微分方程组的外部函数(lorenz.m)function xdot=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)];再调用ode45()求解，输入：[t,x]=ode45(@lorenz,[0,100],[0 0 1e-10]);plot(x(:,1),x(:,2));axis([-1 50 -20 20 ]);xlabel(x1),ylabel(x2)x1与x2的关系图如下：x2与x3的关系图如下：x1与x3的关系图如下：令=10，=10，=8/3，初始点为(0 0 10-10)x1与x2的关系图如下：x2与x3的关系图如下：x1与x3的关系图如下：令=5，=28，=8/3初始点为(0 0 10-10)x1与x2的关系图如下：x2与x3的关系图如下：x1与x3的关系图如下：令=10，=28，=10 初始点为(1 1 10-10)x1与x2的关系图如下：x2与x3的关系图如下：x1与x3的关系图如下：2.3 加分实验餐厅废物的堆肥优化问题1 问题分析堆肥是利用微生物的分解作用将有机废物转化成无害稳定形式的生物化学过程，要提高堆肥产量方法之一就是通过增强细菌的生长繁殖能力以提高分解率。细菌群体的增长一般要经历延滞期，加速生长其、对数生长期、减数生长期、稳定期、加速死亡期和对数死亡期。另外，对数生长期培养基中所有养分都过剩，细菌可以充分繁殖，其倍增速率恒定，取决于底物浓度、温度、水活度、供氧量。对于当前该餐厅来说底物浓度由每天的剩余食物、蔬菜叶和碎纸屑决定。碎纸屑是吸收水分的调理剂，微菌呼吸所需要的溶解氧由蔬菜提供，水活度可以通过测定相对湿度来决定，其关系式是：相对湿度：B=P/P0*100%水活度：aw=P/P0其中P为该容易蒸汽压，P0为纯水蒸汽压在这个堆肥系统中，可供微菌消耗的营养底物和溶解氧都是有限的，他们的消耗会对微菌生长率产生重要影响，一种微菌消耗营养底物的速度和底物浓度之间的关系式为：ds/dt=-KmSX/(Ks+s),式中，ds