数学建模实验二详解.doc

数学建模实验二 微分方程实验 基本实验 微分方程稳定性分析 绘出下列自治系统相应的轨线，并标出随t增加的运动方向，确定平衡点，并按稳定的、渐近稳定的、或不稳定的进行分类： 解答：（1）由平衡点的定义可得，fx x 0，fy y 0，因此平衡点为0,0 ，对照稳定性的情况表，可知平衡点 0, 0 是不稳定的。 由平衡点的定义可得，fx -x 0，fy 2y 0，因此平衡点为0,0 ，对照稳定性的情况表，可知平衡点 0, 0 是不稳定的。 由平衡点的定义可得，fx y 0，fy -2x 0，因此平衡点为0,0 ，对照稳定性的情况表，可知平衡点 0, 0 是不稳定的。 由平衡点的定义可得，fx -x 0，fy -2y 0，因此平衡点为0,0 ，对照稳定性的情况表，可知平衡点 0, 0 是稳定的。 营养平衡问题 营养以每单位时间R个分子的常速流入一个细胞，并且以其内营养浓度成比例的速度离开，比例常数为K，设N为t时刻的浓度，则上述营养变化速度的数学描述为： 即N的变化速度等于营养进入细胞的速度减去它们离开的速度，营养的浓度会达到平衡吗？如果能，平衡解是什么？它是稳定的吗？试用这个方程解的图示解释之。 解答：由题意可得N满足的微分方程为：， 令，可求得方程的平衡点，当时，；当时，. 不难计算出，由题意知,故平衡点是稳定的。 由以上分析可做图示分析如下： 种群增长模型 一个片子上的一群病菌趋向于繁殖成一个圆菌落.设病菌的数目为N，单位成员的增长率为，则由Malthus生长律有，但是，处于周界表面的那些病菌由于寒冷而受到损伤，它们死亡的数量与成比例，其比例系数为，求N满足的微分方程.不用求解，图示其解族.方程是否有平衡解，如果有，是否为稳定的? 考虑单种群开发方程 用数学表达式证明：在稳定状态下，最优捕捞率为 解：由本问题的目标出发，渔场中鱼量达到稳定的平衡状态时的情形，不必知道每一时刻的鱼量变化情况，故不需要解出方程，只需要讨论方程的平衡点并分析其稳定性。 令 解得的两个平衡点为： ， 易得 由定理： 若，则是稳定的；若，则不是稳定的。 应用上述近似判别法，所以有 当E r时，是稳定平衡点，不是稳定平衡点； 当E r时，不是稳定平衡点，是稳定平衡点； 所以，当捕捞适度 即：E r 时,可使渔场产量稳定在，从而获得持续产量，而当捕捞过度（即：E r）时，渔场产量将减至，从而是不可持续的。 令 求导可得： 所以，最优捕捞率为。 设渔场鱼量自然增长服从Gompertz模型： 其中r为固有增长率，N为最大种群数量。若单位时间捕捞量为．讨论渔场鱼量的平衡点及其稳定性，求最大持续产量及获得最大产量的捕捞强度和渔场鱼量水平。 ：令 由，，解得该方程的平衡点为，． ， 可得平衡点是稳定的，而平衡点不稳定. 最大持续产量的数学模型为： 由前面的结果可得 ， 令可得最大产量的捕捞强度，从而得到最大持续产量，此时，渔场鱼量水平。 微分方程 是两个物种为了共同的有限资源而竞争的模型，假设c1 a1，c2 a2。试用微分方程稳定性理论分析： 如果，则如果则用图形分析方法来说明上述两种情况。 令 得方程的平衡点为P0（0,0），P1（,0），P2（0, ）. 对平衡点P0（0,0）， 系数矩阵 又c1 a1，c2 a2则p -[ c1-a1 + c2-a2 ] 0，所以该平衡点不稳定。 对平衡点P1（,0），系数矩阵 则p ，q ，若，且c1 a1，c2 a2，则q 0不稳定。 而对于P2（0,），有p 0，且q 0稳定，此时说明物种1最终要灭亡。 2 如果的情况下，方程在P1（,0）稳定，其他点不稳定，此时说明物种2最终会灭亡。 3 对于线性方程组 ， 其中，直线将第一象限分成三个区域。 当时，P2点稳定，通过分析的单调性可得下图： 此时说明物种1最终要灭亡。 当时，P1点稳定，通过分析的单调性可得下图： 此时说明物种2最终要灭亡。 考虑Lorenz模型 其中，且初值为,，ε为一个小常数，假设，且。 （1）用函数ode45求解，并画出x2~x1,x2~x3,x3~x1的平面图； （2）适当地调整参数σ，ρ，β值，和初始值x10 ，x20 0，x30 ，重复一的工作，看有什么现象发生。 解： （1）创建.m文件，在.m文件中编写下面的语句： 在.m文件中，编写下面的语句：f @ t,x [-8/3*x 1 +x 2 *x 3 ;-10*x 2 +10*x 3 ;-x 1 *x 2 +28*x 2 -x 3 ]; t_final 100;x0 [0;0;1e-10];[t,x] ode45 f,[0,t_final],x0 ; subplot 2,2,1 ,plot x :,2 ,x :,1 xlabel x_2 ;yla