[第三讲]数模-微分方程-讲座.pptVIP

常微分方程理论在数学建模中的应用 北京理工大学 王宏洲 1、微分方程的主要适用范围 2、微分方程模型的分析方法 4、建立微分方程模型的依据 5、案例－物体的运动、振动、受力形变 6、案例－生物的数量变化或密度变化 案例一、传染病模型 模型目标 模型假设 1、SI模型（只考虑S和I两类人） (2) 人群数量足够大，只考虑传播过程中的平均效应，即函数s(t)和i(t)可以视为连续且可微的。 (3) 每个I类的人每天“有效接触”的人数（包括病人、健康人）为常数λ。这个常数实际上就是传染率，反映本地区的卫生水平。 (4) 不考虑出生与死亡，以及人群的迁入迁出因素。（简化问题） 构造模型 模型求解 SI模型图形分析 SI模型结果分析 2、SIS模型（可治愈但不免疫模型） 1.假设(前面四条都和模型A一样，再添加一条) （5）病人以固定的比率痊愈，再次成为易感人群。每天被治愈的病人数占病人总数的比例为μ。 构造模型 模型求解 模型结果分析 3、SIR模型（免疫模型） 1、假设：这里的假设类似于模型B，只是引入R类人群。分别记s(t)、i(t)、r(t)为病人、易感人群、移出者在总人口中所占的比例。s(t)+ i(t)+ r(t) = 1。另外，日接触率λ，日治愈率μ。 模型求解 将方程组转化成下面的形式： 模型分析 考察随着时间的推移，s(t)、i(t)、r(t)的变化规律。 首先，t →+∞时，分别以s ? , i ? , r?记各自的极限，这些极限都存在。 模型分析 i ? = 0 ?（用反证法） 假设i ? ? 0 ，那么必然有 i ? = ? 0。 根据极限的定义，对于充分大的t，都应该有i(t)ε/2，把这个结论代入方程组。 模型分析 其次，考虑随着t的变化，i-s平面上解的轨线变化情况。大概的走势图为： 模型分析 1/σ是一个边界点，为了让传染病不蔓延，需要调整s0和1/σ。具体的方法：一是降低s0，如接种疫苗，使S类人群直接变成R类； 二是提高1/σ使之大于s0，σ=λ/μ,也就是降低λ而提高μ，强化卫生教育和隔离病人，同时提高医疗水平。 参数估计 对参数σ的估计： 令解两端同时取t→+∞，因为 i ? = 0 ，得到 其他类型的传染病模型 SIES模型——健康－染病－潜伏期－健康不免疫 SIER模型——健康－染病－潜伏期－移出系统 SIRS模型——健康－染病－短时免疫－健康(易感) 考虑抵抗能力 考虑地域传播 考虑传播途径(接触、空气、昆虫、水源等) 数学模型 ? 微分方程稳定性方法建模 北京理工大学 王宏洲 事物发展的稳定与不稳定 微分、差分方程稳定性理论简介 常微分方程稳定性理论 差分方程稳定性理论 平衡点稳定的几何特征 一阶微分方程 2、二阶微分方程 二阶微分方程 二阶微分方程 二阶微分方程 3、一阶线性差分方程 4、二阶线性差分方程 5、一阶非线性差分方程 一、两个生物种群的竞争模型 考虑两个生物种群竞争同一种有限资源的问题。在自然条件下，适应环境能力弱的种群将趋于灭亡，适应能力强的种群将增长到环境允许的最大数量。种群竞争模型现在已经被广泛地应用到描述企业、国家等社会实体之间的竞争研究中。 下面通过建立模型来解释这种现象，并分析出现各种结局的条件。 2. 稳定性分析（竞争的结局） 2.1 求平衡点 令f(x1,x2)=g(x1,x2)=0，得到四个平衡点： P1(N1,0)，P2(0,N2)， P3(0,0)， 一(2)生物互惠共生模型 第一种情形 平衡点稳定性分析 平衡点有三个：P1(N1, 0), 模型结果分析 一(3) 食饵-捕食者模型 平衡点稳定性分析 稳定性分析 模型结果分析 O N1 x1 x2 N2 2、?11，?21，P2稳定 3、 ?11，?21，P3稳定 O N1 N1/ ?2 x1 P3 N2 N2/ ?1 x2 O N1/ ?2 N2 x1 P3 x2 N2 N2/?1 4、?11，?21，方程的解不存在统一的发展趋势。 甲乙两种群的相互依存有三种形式 1) 甲可以独自生存，乙不能；甲乙一起生存时相互提供食物、促进增长。 2) 甲乙均可以独自生存；甲乙一起生存时相互提供食物、促进增长。 3) 甲乙均不能独自生存；甲乙一起生存时相互提供食物、促进增长。 模型假设 甲可以独自生存，数量变化服从Logistic规律; 甲乙一起生存时乙为甲提供食物、促进增长。 乙不能独自生存；甲乙一起生存时甲为乙提供食物、促进增长