微分方程本科蒋鲁敏.ppt

* 类似问题 问题: 如图所示,一条猎犬发现它的正西方1里处有一只野兔正朝北边1里处的兔巢逃奔,猎犬随即朝兔子奔走方向追赶.已知兔子奔跑速度为0.42里/分钟,猎犬速度为兔子速度的2倍,试问猎犬能否在野兔逃回兔巢前追上野兔? * 类似问题 问题(1988AMCM_A) 确定毒品走私位置： 走私船走向随机 发现范围是一个区域 寻找最大可能发现毒品船对策略 * 例三：人口问题和生态问题 马尔萨斯人口模型 有资源限制的人口模型 勒斯里模型 * 前两种模型 * 勒斯里模型 人口与年龄、时间有关：两个自变量 归结为含偏导数的微分方程-偏微分方程(参考文献) 离散化后得到差分方程-勒斯里模型 简化了的勒斯里模型 假设： 三种鱼:一龄鱼、二龄鱼、三龄鱼 只有二龄鱼能生育 * 简化了的勒斯里模型 N1的说明 对t积分离散化，按1年周期计算 * 生态模型与 ode45的用法 >> u0=[1;1] u0 = 1 1 >> [t,u]=ode45('aa',[0,6],u0); >> plot(t,u) function rw=aa(t,u) rw=[(3-2*u(2))*u(1);(-2.5+u(1))*u(2)]; * * 相轨线 * 历届数模竞赛中与微分方程有关的问题 注意：不一定纯粹为微分方程问题，和统计、计算方法等知识结合起来 1996A 最优捕鱼策略 2003A SARS的传播 2006B 艾滋病疗法的评价及疗效的预测 2007A 中国人口增长预测 * * * * * * * * * * * i0 i 0 t ? <1 di/dt < 0 发展趋势 * 用matlab计算例子 l=0.4 s=0.8<1 1-1/s=-0.25 i0=0.2(上图);0.1(下图) >>ezplot('1/(-4+9*exp(1/10*t))',0,100)(上图) >>ezplot(‘1/(6+4*exp(-1/15*t))’,0,100)(下图) * 模型的意义: 传染病机制解释 传染病终趋消失的条件——? < 1: 感染因素弱于治疗因素 传染病无法消亡的情况——? > 1: 感染因素强于治疗因素 * 传染病模型IV： 增加假设:病人治愈后不会再被感染 三种人：健康人、病人和病愈者 总人数 N 不变，记病人、健康人和病愈者的比例分别为 记病人的日传染? , 日治愈率?, 传染数 ? = ? / ? 需建立含 的两个方程-常微分方程组 SIR模型 * 建模 ? = ? / ? 变化率 * 解模 无法求出 的解析解 在相平面 上 研究解的性质 (s(t),i(t))可以看做(s,i)平面上的曲线(相轨线)的参数式 * 研究模型：相轨线方程 * 1 1 s i 0 D 在D内作相轨线 (s(t),i(t))可以看做(s,i)平面上的曲线(相轨线)的参数式 相轨线 * 相轨线走势 * 最大病人数条件 * 相轨线分析 s i 1 0 1 D s(t)单调减?相轨线向左变化 P1 s0 im P3 P4 P2 S0 相轨线起点在斜线附近 * 相轨线分析结论 s i 1 0 1 D 传染病不蔓延 P1 s0 im P2: s0<1/s ? i(t)单调 降至0 P3 P4 P2 S0 * 相轨线分析结论 s i 1 0 1 D 传染病有蔓延过程 P1 s0 im P1: s0>1/s ? i(t)先升后降至 0 1/s ~阈值 P3 P4 P2 S0 * 模型的意义: 预防传染病蔓延 ? (日传染率)? ? 隔离 ? (日治愈率)? ? 提高医疗水平 传染病不蔓延的条件——s0<1/? ? 的估计 降低 s0 提高 r0 提高阈值 1/? 降低 ?(=?/?) ? ?, ? ? 打预防针-群体免疫 更多传染病模型(SEIR) * * SEIRD 模型 * 微分方程建模步骤 根据建模目的和问题分析作简化假设 根据函数及其变化率之间的关系建立 方程 定解条件 按照正确方法解微分方程 讨论解的性质 * 方程建模应注意 物理量纲： 在等式两边中应注意物理单位相同，比例系数有时也是有量纲单位的 列出定解条件： 定解条件是系统在某一特定时刻或(和)边界上的信息。独立于方程，用以确定有关通解中的常数。为了完整充分地给出问题的数学陈述，应将这些给定的条件和方程一起列出。 * 微分方程的进一步分析 对初值、参数依赖性的灵敏度分析 可能会引起解的较大变化，也可能不大影响解的值 对求解过程中的误差控制有一定的要求 方法可以写出解对初值、参数的导数公式并作讨论；或用不同初值、参数求解，并作比较