常 微 分 方 模 型 选 讲.ppt

* 生态模型与 ode45的用法 x0=[1;1] x0 = 1 1 [t,x]=ode45(aa,[0,6],x0); plot(t,x) function rw=aa(t,x) rw=[(3-2*x(2))*x(1);(-2.5+x(1))*x(2)]; * * 相轨线 * 历届数模竞赛中与微分方程有关的问题 注意：不一定纯粹为微分方程问题，和统计、计算方法等知识结合起来 1996A 最优捕鱼策略 2003A SARS的传播 2006B 艾滋病疗法的评价及疗效的预测 2007A 中国人口增长预测 * * * * * * * * * * * 用matlab计算例子 l=0.4 s=0.81 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 不变，记病人、健康人和病愈者的比例分别为 记病人的日传染? , 日治愈率?, 传染数 ? = ? / ? 需建立 的两个方程 * 建模 ? = ? / ? * 解模 无法求出 的解析解 在相平面 上 研究解的性质 (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: s01/s ? i(t)单调 降至0 P3 P4 P2 S0 * 相轨线分析结论 s i 1 0 1 D 传染病蔓延 P1 s0 im P1: s01/s ? i(t)先升后降至 0 1/s ~阈值 P3 P4 P2 S0 * 模型的意义: 预防传染病蔓延 ? (日传染率)? ? 隔离 ? (日治愈率)? ? 提高医疗水平 传染病不蔓延的条件——s01/? ? 的估计 降低 s0 提高 r0 提高阈值 1/? 降低 ?(=?/?) ? ?, ? ? 打预防针-群体免疫 * 微分方程建模步骤 根据建模目的和问题分析作简化假设 根据函数及其变化率之间的关系建立 方程 定解条件 按照正确方法解微分方程 讨论解的性质 * 方程建模应注意 物理量纲： 在等式两边中应注意物理单位相同，比例系数有时也是有量纲单位的 列出定解条件： 定解条件是系统在某一特定时刻或(和)边界上的信息。独立于方程，用以确定有关通解中的常数。为了完整充分地给出问题的数学陈述，应将这些给定的条件和方程一起列出。 * 微分方程的进一步分析 对初值、参数依赖性的灵敏度分析 可能会引起解的较大变化，也可能不大影响解的值 对求解过程中的误差控制有一定的要求 方法可以写出解对初值、参数的导数公式并作讨论；或用不同初值、参数求解，并作比较 对有平衡态的情况，作稳定性分析 可参考任何一本比较详细的常微分方程教科书 * 例二、鱼雷击舰问题 问题: 如图所示,一敌舰在某海域内沿正北方向航行时, 我方战舰恰位于敌舰正西方1海里处, 我舰向敌舰发射制导鱼雷. 敌舰速度为0.42海里/分钟,鱼雷速度为敌舰速度的2倍, 鱼雷的运行方向始终指向敌舰. 试问敌舰航行多远时将被击中? * 建 模 目 的 研究 x、y 随 t 变化的关系 假设： 鱼雷的运行方向始终指向敌舰 设敌舰速度v0、鱼雷速度v为常数 记号:(建立坐标系后) 敌舰速度为 v0；t 时刻敌舰的位置为 (1, v0t) 鱼雷速度为 v=2 v0 ；t 时刻鱼雷的位置为 (x, y) 模型假设与记号 * 1. 建立微分方程模型 * 解模 * 2. 模型求解 问题解答: 当 x=1 时, y=2/3, 时间 t=y/v0 =1.59 分=95.24秒时鱼雷 击中敌舰. * Matlab求解 dsolve(D2u=sqrt(1+(Du)*(