第五章数学建模.docVIP

第五章 稳定性方法建模 对于某些实际问题，建立数学模型的目的并不是要寻求动态过程的每个时刻的性态，而往往是要研究经过充分长的时间后动态过程的变化趋势.即是否会越来越接近某个确定的值.若是，则称为稳定，否则称为不稳定.在分析这种稳定性时，常常不用解微分方程,而可利用微分方程的稳定性理论，直接研究. §5.1 预备知识 （一）一阶微分方程的平衡点及稳定性 现在我们来介绍一阶微分方程的平衡点及稳定性的概念． 设 ， (5.1.1) 即右端不显含自变量t，则方程 f(x)=0 (5.1.2) 的实根称为微分方程(5.1.1)的平衡点. 显然是方程(5.1.1)的一个解. 如果当x(0)充分接近x0时，方程(5.1.1)的解x(t)满足 (5.1.3) 则称平衡点是稳定的，否则称是不稳定的. 对于平衡点的稳定性有如下的判别法: 设是(5.1.1)的平衡点.若 ,则是稳定的, 若,则是不稳定的. 假设一个质点的运动可以用方程(5.1.1)描述, 即当它在x处时运动速度为f(x),仅与它的位置有关. 设, 则在附近, f(x)递减, 故若x稍小于, 则有, 就是说, 当质点在x0 的左面时, 它会向右运动. 而当x稍大于时, 则有, 故当它在x0 的右面时, 它会向左运动. 总之, 当质点稍微偏离平衡点时,会自动回到平衡点,故是稳定平衡点. 反之, 若, 则当质点在x0 的左面时, 它会向左运动, 偏离平衡点越来越远, 当质点在x0 的右面时, 它也会偏离平衡点越来越远,故不是稳定平衡点. 根据稳定平衡点的定义,(5.1.1) 的稳定平衡点可以作为该方程的一个近似解.这给一些实际问题的研究带来方便. (二）二阶微分方程的平衡点与稳定性 平衡点的概念 此时方程可用两个一阶方程表为 (4) 右边不显含t，方程组 (5) 的实根，称为方程（4）的平衡点，记为 稳定性的概念 若不论初始条件如何，方程（4）的解都满足 , (6) 则称平衡点是稳定的，否则，称是不稳定的。 平衡点稳定性判别法 记，则当使tr(A)0且det(A)0时，是稳定的； 当使tr(A)0或det(A)0时，是不稳定的。 §5.2 捕鱼业的产量模型 5.2.1问题描述 设渔场的鱼量按一定规律增长,在渔场捕鱼时，如何控制捕捞率，使单位时间捕鱼量达到最大. 5.2.2模型假设 1．设时刻t渔场的鱼量为x(t),在无捕捞条件下，x(t)服从Logistic规律.即 (5.2.1) 其中，r为固有增长率.N为环境允许的最大鱼量.f(x)表示单位时间鱼的增长量（自然增长率）. 2、设单位时间捕鱼量为h(x)，它与渔场鱼量x(t)成正比，比例系数为E，E称为捕捞率.即 h(x) = E x (0≤E≤1) 3、假设是独家捕捞的情形 5.2.3建模构建 记F(x)=f(x)–h(x)，则F(x)表示该渔场鱼量的净增长率. 在有捕捞情况下，该渔场鱼量满足方程 (5.2.2) 5.2.4模型分析 一般情况,我们不一定要解方程(5.2.2)求得x(t)的动态变化过程,而只希望知道x(t)的趋向, 并由此确定最大持续产量. 为此我们求(5.2.2)的平衡点和分析其稳定性. 先求(5.2.2)的平衡点. (x≥0) 令F(x) = 0，得出两个平衡点：x0 = N ( 1 – E/r )，x1 = 0 因要求x0≥0, 故知0≤E≤r. 再分析平衡点稳定性 若E r, F’(x0)0, F’(x1)0, 故x0是稳定平衡点，x1不稳定. 若E=r，则,(5.2.2)化为，立刻解得 , 因此. 于是总有 (5.2.3) 现在考虑在稳定条件下(即当t充分大时)，捕捞率E为多大时, h(x)达到最大？ 由(5.2.3)式可知，在稳定条件下，要使h(x) = E x 达到最大 则E应满足0Er,又t充分大时， 此时 h(x0)=x0= 把其视为自变量E的函数,记为 的驻点为，此时0，故当时，最大，也即最大，此时 ， 说明若把捕捞率E控制在鱼量固有增长率的一半时，则可获得最大的持续产量. 注意：Er时， (5.2.2)的解为 对应的曲线形状： §5.3 捕鱼业的效益模型 5.3.1问题描述 从经济角度分析，不应追求单位时间产量最大，而应追求单位时间获利最大.接上节考虑如下问题： 捕捞率控制在什么水平时单位时间获利最大？ 5.3.2模型假设 1.鱼的销售单价为常数P； 2.捕鱼的成本与捕捞率E成正比，比例系数为常数C； 3.单位时间所获利润为R(E)； 4.独家捕捞. 5.3.3模型构建 在稳定条件下,, 从销