数学建模 微分方程模型

2传染病模型 3战争模型 4最优捕鱼问题 1微分方程模型 微分方程模型 1微分方程模型 一 微分方程模型的建模步骤在自然科学以及工程 经济 医学 体育 生物 社会等学科中的许多系统 有时很难找到该系统有关变量之间的直接关系 函数表达式 但却容易找到这些变量和它们的微小增量或变化率之间的关系式 这时往往采用微分关系式来描述该系统 即建立微分方程模型 我们以一个例子来说明建立微分方程模型的基本步骤 例1某人的食量是10467 焦 天 其中5038 焦 天 用于基本的新陈代谢 即自动消耗 在健身训练中 他所消耗的热量大约是69 焦 公斤 天 乘以他的体重 公斤 假设以脂肪形式贮藏的热量100 地有效 而1公斤脂肪含热量41868 焦 试研究此人的体重随时间变化的规律 模型分析在问题中并未出现 变化率 导数 这样的关键词 但要寻找的是体重 记为W 关于时间t的函数 如果我们把体重W看作是时间t的连续可微函数 我们就能找到一个含有的微分方程 模型假设1 以W t 表示t时刻某人的体重 并设一天开始时人的体重为W0 2 体重的变化是一个渐变的过程 因此可认为W t 是关于连续t而且充分光滑的 3 体重的变化等于输入与输出之差 其中输入是指扣除了基本新陈代谢之后的净食量吸收 输出就是进行健身训练时的消耗 模型建立问题中所涉及的时间仅仅是 每天 由此 对于 每天 体重的变化 输入 输出 由于考虑的是体重随时间的变化情况 因此 可得体重的变化 天 输入 天 输出 天 代入具体的数值 得输入 天 10467 焦 天 5038 焦 天 5429 焦 天 输出 天 69 焦 公斤 天 公斤 69 焦 天 体重的变化 天 W t 公斤 天 当 t 0时 它等于dW dt 考虑单位的匹配 利用 公斤 天 焦 每天 41868 焦 公斤 可建立如下微分方程模型 模型求解用变量分离法求解 模型方程等价于积分得 从而求得模型解就描述了此人的体重随时间变化的规律 现在我们再来考虑一下 此人的体重会达到平衡吗 显然由W的表达式 当t 时 体重有稳定值W 81 我们也可以直接由模型方程来回答这个问题 在平衡状态下 W是不发生变化的 所以这就非常直接地给出了W平衡 81 所以 如果我们需要知道的仅仅是这个平衡值 就不必去求解微分方程了 至此 问题已基本上得以解决 一般地 建立微分方程模型 其方法可归纳为 1 根据规律列方程 利用数学 力学 物理 化学等学科中的定理或许多经过实践或实验检验的规律和定律 如牛顿运动定律 物质放射性的规律 曲线的切线性质等建立问题的微分方程模型 3 模拟近似法 在生物 经济等学科的实际问题中 许多现象的规律性不很清楚 即使有所了解也是极其复杂的 常常用模拟近似的方法来建立微分方程模型 建模时在不同的假设下去模拟实际的现象 这个过程是近似的 用模拟近似法所建立的微分方程从数学上去求解或分析解的性质 再去同实际情况对比 看这个微分方程模型能否刻划 模拟 近似某些实际现象 本章将结合例子讨论几个不同领域中微分方程模型的建模方法 2传染病模型 问题 描述传染病的传播过程 分析受感染人数的变化规律 预报传染病高潮到来的时刻 预防传染病蔓延的手段 按照传播过程的一般规律 用机理分析方法建立模型 已感染人数 病人 i t 每个病人每天有效接触 足以使人致病 人数为 模型1 假设 若有效接触的是病人 则不能使病人数增加 建模 模型2 区分已感染者 病人 和未感染者 健康人 假设 1 总人数N不变 病人和健康人的比例分别为 2 每个病人每天有效接触人数为 且使接触的健康人致病 建模 日接触率 SI模型 模型2 tm 传染病高潮到来时刻 日接触率 tm 病人可以治愈 t tm di dt最大 模型3 传染病无免疫性 病人治愈成为健康人 健康人可再次被感染 增加假设 SIS模型 3 病人每天治愈的比例为 日治愈率 建模 日接触率 1 感染期 一个感染期内每个病人的有效接触人数 称为接触数 模型3 接触数 1 阈值 感染期内有效接触感染的健康者人数不超过病人数 模型2 SI模型 如何看作模型3 SIS模型 的特例 模型4 传染病有免疫性 病人治愈后即移出感染系统 称移出者 SIR模型 假设 1 总人数N不变 病人 健康人和移出者的比例分别为 2 病人的日接触率 日治愈率 接触数 建模 需建立的两个方程 模型4 SIR模型 模型4 SIR模型 相轨线的定义域 在D内作相轨线的图形 进行分析 模型4 SIR模型 相轨线及其分析 s t 单调减 相轨线的方向 P1 s0 1 i t 先升后降至0 P2 s0 1 i t 单调降至0 1 阈值 模型4 SIR模型 预防传染病蔓延的手段 日接触率 卫生水平 日治愈率 医疗水平 传染病不蔓延的条件 s0 1 的估计 降低s0 提高r0 提高阈值1 模型4 SIR模型 被传染人数的估计 记被传染人数比例 小 s0 1 提高阈值1 降低被传染人数比例x s0 1 战争分类 正规战争 游击战争 混合战争 只考虑双方兵力多少和战斗力强弱 兵力因战斗及非战斗减员而减少 因增援而增加 战斗力与射击次数及命中率有关 建模思路和方法为用数学模型讨论社会领域的实际问题提供了可借鉴的示例 第一次世界大战Lanchester提出预测战役结局的模型 3战争模型 一般模型 每方战斗减员率取决于双方的兵力和战斗力 每方非战斗减员率与本方兵力成正比 甲乙双方的增援率为u t v t f g取决于战争类型 x t 甲方兵力 y t 乙方兵力 模型假设 模型 正规战争模型 甲方战斗减员率只取决于乙方的兵力和战斗力 双方均以正规部队作战 忽略非战斗减员 假设没有增援 f x y ay a 乙方每个士兵的杀伤率 a rypy ry 射击率 py 命中率 正规战争模型 为判断战争的结局 不求x t y t 而在相平面上讨论x与y的关系 平方律模型 游击战争模型 双方都用游击部队作战 甲方战斗减员率还随着甲方兵力的增加而增加 f x y cxy c 乙方每个士兵的杀伤率 c rypyry 射击率py 命中率 游击战争模型 线性律模型 混合战争模型 甲方为游击部队 乙方为正规部队 乙方必须10倍于甲方的兵力 设x0 100 rx ry 1 2 px 0 1 sx 1 km2 sry 1 m2 再生资源 渔业 林业等 与非再生资源 矿业等 再生资源应适度开发 在持续稳产前提下实现最大产量或最佳效益 问题及分析 在捕捞量稳定的条件下 如何控制捕捞使产量最大或效益最佳 如果使捕捞量等于自然增长量 渔场鱼量将保持不变 则捕捞量稳定 背景 4最优捕鱼问题 产量模型 假设 无捕捞时鱼的自然增长服从Logistic规律 单位时间捕捞量与渔场鱼量成正比 建模 捕捞情况下渔场鱼量满足 不需要求解x t 只需知道x t 稳定的条件 r 固有增长率 N 最大鱼量 h x Ex E 捕捞强度 x t 渔场鱼量 一阶微分方程的平衡点及其稳定性 一阶非线性 自治 方程 F x 0的根x0 微分方程的平衡点 不求x t 判断x0稳定性的方法 直接法 1 的近似线性方程 产量模型 稳定性判断 x0稳定 可得到稳定产量 x1稳定 渔场干枯 E 捕捞强度 r 固有增长率 产量模型 在捕捞量稳定的条件下 控制捕捞强度使产量最大 图解法 P的横坐标x0 平衡点 P的纵坐标h 产量 产量最大 控制渔场鱼量为最大鱼量的一半 效益模型 假设