[理学]第五章 微分方程模型.ppt

一、简单的微分方程模型 二、差分方程 三、人口模型 四、方程组一 战争模型 五、方程组二、传染病模型 六、食饵—捕食者模型 七、微分方程的稳定性分析 应用 田鼠和猫头鹰 田鼠和其天敌猫头鹰在田间的群种消长规律也服从 Volterra模型. 根据一段时间的观察, 得到下面一组数据对: 184 152 104 88 72 80 96 112 144 136 168 111.1 230.6 268.5 183.3 92.5 48.1 33.6 22.6 10.9 25.7 41.9 192 184 128 104 88 72 96 120 120 152 168 95.0 195.7 272.7 215.8 124.6 57.6 37.6 21.7 20.7 34.0 39.8 192 176 144 104 88 96 88 104 128 猫头鹰 69.6 162.3 269.6 236.0 134.2 69.1 32.5 33.1 29.7 田鼠 田鼠和猫头鹰数据 试建立相应的模型和曲线形式. 模型分析 首先做出散点图: 可以看到, 该曲线符合Volterra模型特征. 如同前面分析, 设 是食饵—田鼠的数量, 而 是其天敌 猫头鹰的数量, 则有方程 问题是如何得到其中的未知参数. 将上面的方程变形, 得到 即有 用数值微分取代微分, 再用矩阵形式的曲线拟合方法, 得 相应的曲线为 相轨线为 微分方程的稳定性理论是微分方程模型中的重要内容之 一, 它对于模型解的讨论有着重要的意义. 1.稳定性的概念. 设微分方程 称使 的点 为方程的平衡点. 若 则称 为方程稳定的平衡点. ⑶ 自动力系统 2.稳定性的判定 因 由一阶泰勒公式: 代入方程有 ⑷ 该方程的通解为 ⑸ 线性近似 所以 是稳定平衡点 二阶微分方程的稳定性 一个二阶微分方程一般可以写成一个一阶微分方程组 的形式, 即有 的形式. ⑹ 该微分方程的平衡点为满足 的点 若有关系 则点 称为稳定平衡点. 稳定性分析 将⑹改写成近似线性关系 为简单起见, 做个平移变换, 上式可以简单写成 ⑺ 此时平衡点为 记 并假设 则 点是方程的惟一平衡点. 点 的稳定性由矩阵 的特征方程 的根 （特征根）决定. 又 此时 相应的特征值为 方程⑺的解具有形式 或 容易看到, 点 是稳定平衡点 又由已知条件 由此得到: 不稳定 中心 不稳定 不稳定焦点 稳定 稳定焦点 不稳定 不稳定退化结点 稳定 稳定退化结点 不稳定 鞍点 不稳定 不稳定结点 稳定 稳定结点 稳定性 平衡点类型 由特征方程决定的平衡点类型和稳定性 由该表可以得到方程稳定性讨论的一般法则. 应用 原微分方程为 此时方程有平衡点 （舍去）及 食饵—捕食者模型的稳定性分析 由方程组 所以 该点的稳定性不能用近似线性方程进行分析 用相轨线对稳定性进行分析 由方程组 消去 变形为 积分后得到 取指数后有 其中的 值由初始条件确定. 则相轨线为 在相平面上讨论 相线的图形 结论 当 时相轨线不存在. 故设 相轨线 中心 相轨线退化为 点 取 再令 即有点 对于点 因 存在 使 从而有 所以曲线为封闭曲线族. 注意到: 相轨线是封闭曲线 是周期函数. 函数在一个周期内的平均值 同理: 轨线中心: 即 注 在问题一中, 按此公式计算结果为 这与数值积分的结 果基本吻合. 传染病传染过程研究有悠久的历史, 其数学模型的建立、 发展和改善和人口模型一样有一个过程. 最早的传染病模 型是Danel Bollinie 1760年的天花模型来评估种痘效果, 以 后众多学者不懈努力进行不断修正. 下面, 我们通过三个 模型看一看这个进化过程. 这些模型都是研究传染病传染 过程的, 都是基于一个基本传染机制, 即假定病人是通过 与他人接触而将病原体传染给他人的. 模型Ⅰ 模型 模型假定 ⑴疫区封闭, 即总人数 常数, 其中病人数 （ill people） 并假设函数有连续导数; ⑵在单位时间内一个病人能接触到的人数为定量, 记作 称为接触率, 并将接触到的人中的健康人传染成病人; ⑶初始时