数学建模之微分方程方法.pptVIP

数学建模及其应用微分方程模型 在研究实际问题时，常常会联系到某些变量的变化率或导数， 这样所得到变量之间的关系式就是微分方程模型。微分方程模型反映的是变量之间的间接关系，因此，要得到直接关系，就得求解微分方程。 【高崖估算问题】假如你站在崖顶且身上带着一只具有跑表功能的计算器，你也许会出于好奇心想用扔下一块石头听回声的方法来估计山崖的高度，假定你能准确地测定时间，你又怎样来推算山崖的高度呢?请分析一下这一问题。 一 问题分析与假设 假定空气阻力不计，可以直接利用自由落体运动的公式来计算. 例如，设t=4秒，g=9.81米/秒2，则可求得h≈78.5米。 ● 但实际上空气阻力却是不得不考虑的因素。根据流体力学知识，此时可设空气阻力正比于石块下落的速度，阻力系数k为常数。 二 模型建立与求解 由牛顿第二定律可得方程： 方程的解为 初始条件v(0)=0，代入可得C=-mg/k，于是 求解得 又根据问题可知，初始条件h(0)=0，代入可得 ，于是 内容纲要 一 微分方程的平衡点及稳定性 二 案例分析：传染病模型 三 案例分析：捕鱼业的持续收获模型 方程(3)的一般解具有形式 或 c1, c2为任意常数. 则 二 案例分析：传染病模型 求解微分方程有三种方法： 1）求精确解；2）求数值解（近似解）；3）定性理论方法。 微分方程模型 人口模型 生物种群模型 兰彻斯特（Lanchester） 作战模型等 The End 微分方程的求解方法 Matlab软件介绍 建立微分方程模型的方法 （1）根据规律列方程 利用数学、力学、物理、化学等学科中的定理或经过实验检验的规律等来建立微分方程模型。 （2）微元分析法 利用已知的定理与规律寻找微元之间的关系式，与第一种方法不同的是对微元而不是直接对函数及其导数应用规律。 --------------- （3）模拟近似法 在生物、经济等学科的实际问题中，许多现象的规律性不很清楚，即使有所了解也是极其复杂的，建模时在不同的假设下去模拟实际的现象，建立能近似反映问题的微分方程，然后从数学上求解或分析所建方程及其解的性质，再去同实际情况对比，检验此模型能否刻画、模拟某些实际现象。 --------------- --------------- * 三 模型讨论 一 微分方程的平衡点及稳定性 １．平衡点的概念 * * 问题：如何来断别平衡点的稳定性呢？ * 2. 一阶方程的平衡点及稳定性 * * * * 3.平面方程的平衡点及稳定性 * * 则方程组(3)的特征方程为 * * * * * * * * * * * * 已感染人数 (病人) i(t) 每个病人每天有效接触(足以使人致病)人数为? 模型1 假设 若有效接触的是病人，则不能使病人数增加 必须区分已感染者(病人)和未感染者(健康人) 建模 ? * * 模型2 区分已感染者(病人)和未感染者(健康人) 假设 1）总人数N不变，病人和健康 人的 比例分别为 2）每个病人每天有效接触人数为?, 且使接触的健康人致病 建模 ? ~ 日 接触率 SI 模型 * * 模型2 1/2 tm i i0 1 0 t tm~传染病高潮到来时刻 ? (日接触率)? ? tm? Logistic 模型 (阻滞增长模型) 病人可以治愈！ ? t=tm, di/dt 最大 * * 模型3 传染病无免疫性——病人治愈成为健康人，健康人可再次被感染 增加假设 SIS 模型 3）病人每天治愈的比例为? ? ~日治愈率 建模 ? ~ 日接触率 1/? ~感染期 ? ~ 一个感染期内每个病人的有效接触人数，称为接触数。 * * 模型3 i0 i0 接触数? =1 ~ 阈值 感染期内有效接触感染的健康者人数不超过病人数 1-1/? i0 模型2(SI模型)如何看作模型3(SIS模型)的特例 i di/dt 0 1 ? 1 0 t i ? 1 1-1/? i 0 t ? ?1 di/dt 0 * * 模型4 传染病有免疫性——病人治愈后即移出感染系统，称移出者 SIR模型 假设 1）总人数N不变，病人、健康人和移出者的比例分别为 2）病人的日接触率? , 日治愈率?, 接触数 ? = ? / ? 建模 需建立 的两个方程 * * 模型4 SIR模型