第十五周 微分方程应用建模.pptVIP

微分方程的几个简单实例 Malthus模型与Logistic模型 例5 赝品的鉴定 例6 新产品的推广 §3.3 为什么要用三级火箭来发射人造卫星 §3.4 药物在体内的分布 hy§3.5 传染病模型 §3.6 糖尿病的诊断 §3.7 稳定性问题 §3.8 捕食系统的Volterra方程 §3.9 较一般的双种群生态系统 §3.10 分布参数法建模 例8 人口问题的偏微分方程模型 例9 交通流问题 例7 本章第2节中的Logistic模型 共有两个平衡点：N=0和N=K，分别对应微分方程的两两个特殊解。前者为No=0时的解而后者为No=K时的解。 当NoK时，积分曲线N=N(t)位于N=K的下方；当NoK时，则位于N=K的上方。从图3-17中不难看出，若No0，积分曲线在N轴上的投影曲线（称为轨线）将趋于K。这说明，平衡点N=0和N=K有着极大的区别。 图3-17 定义1 自治系统 的相空间是指以（x1,…,xn）为坐标 的空间Rn。 特别，当n=2时，称相空间为相平面。 空间Rn的点集{(x1,…,xn)}|xi=xi(t)满足(3.28)，i=1,…,n}称为系统的轨线，所有轨线在相空间的分布图称为相图。 定义2 设x0是（3.28）的平衡点，称： （1）x0是稳定的，如果对于任意的ε0，存在一个δ0，只要|x(0)- x0|δ，就有|x(t)- x0|ε对所有的t都成立。 （2）x0是渐近稳定的，如果它是稳定的且 。 微分方程平衡点的稳定性除了几何方法，还可以通过解析方法来讨论，所用工具为以下一些定理。 （3）x0是不稳定的，如果（1）不成立。 根据这一定义，Logistic方程的平衡点N=K是稳定的且为渐近稳定的，而平衡点N=0则是不稳定的。 解析方法 定理1 设xo是微分方程 的平衡点： 若 ,则xo是渐近稳定的 若 ,则xo是渐近不稳定的 证 由泰勒公式，当x与xo充分接近时，有： 由于xo是平衡点，故f(xo)=0。若 ，则当xxo时必有f(x)0，从而x单增；当xxo时，又有f(x)0，从而x单减。无论在哪种情况下都有x→xo，故xo是渐进稳定的。 的情况可类似加以讨论。 高阶微分方程与高阶微分方程组平衡点的稳定性讨论较为复杂，大家有兴趣可参阅微分方程定性理论。为了下两节的需要，我们简单介绍一下两阶微分方程组平衡点的稳定性判别方法。 考察两阶微分方程组： （3.29） 令 ，作一坐标平移，不妨仍用x记x’，则平衡点xo的稳定性讨论转化为原点的稳定性讨论了。将f(x1,x2)、g(x1,x2)在原点展开，（3.29）又可写成: 考察（3.29）的线性近似方程组: （3.30） 其中： 记 λ1、λ2为A的特征值则λ1、λ2是方程： det（A-λI）=λ2- (a+b) λ+ (ad – bc )=0的根 令p=a+d, q=ad-bc=|A|，则 ，记 。 讨论特征值与零点稳定的关系 （1）若△0，可能出现以下情形： ① 若q0，λ1λ20。 当p0时，零点不稳定； 当p0时，零点稳定 若q0，λ1λ20 当c1=0时，零点稳定 当c1≠0时,零点为不稳定的鞍点 ③ q=0，此时λ1=p，λ2=0，零点不稳定。 （2） △=0，则λ1=λ2: λ有两个线性无关的特征向量 当p0时，零点不 稳定 当p0时，零点稳定 ② 如果λ只有一个特征向量 当p≥0时，零点不 稳定 当p0时，零点稳定 （2） △0，此时 若a0，零点稳定 若a=0，有零点为中心的周期解 综上所述：仅当p0且q0时， （3.30）零点才是渐近稳定的；当p=0且q0时（3.30）有周期解，零点是稳定的中心（非渐近稳定）；在其他情况下，零点均为不稳定的。 非线性方程组（3.29）平衡点稳定性讨论可以证明有下面定理成立: 定理2 若（3.30）的零点是渐近稳定的，则（3.29）的平衡点 也是渐近稳定的；若（3.30）的零点是不稳定的，则（3.29） 的平衡点也是不稳定的。 问题背景： 意大利生物学家D’Ancona曾致力于鱼类种群相互制约关系的研究，在研究过