第1讲 微分方程模型.pptVIP

1.5 烟雾的扩散与消失 现象和 问题 炮弹在空中爆炸，烟雾向四周扩散，形成圆形不透光区域。 不透光区域不断扩大，然后区域边界逐渐明亮，区域缩小，最后烟雾消失。 建立模型描述烟雾扩散和消失过程，分析消失时间与各因素的关系。 问题分析 无穷空间由瞬时点源导致的扩散过程，用二阶偏微分方程描述烟雾浓度的变化。 观察的烟雾消失与烟雾对光线的吸收，以及仪器对明暗的灵敏程度有关。 模型假设 1）烟雾在无穷空间扩散，不受大地和风的影响；扩散服从热传导定律。 2）光线穿过烟雾时光强的减少与烟雾浓度成正比；无烟雾的大气不影响光强。 3）穿过烟雾进入仪器的光线只有明暗之分，明暗界限由仪器灵敏度决定。 模型建立 1）烟雾浓度 的变化规律 热传导定律：单位时间通过单位法向面积的流量与浓度梯度成正比 曲面积分的奥氏公式 1）烟雾浓度 的变化规律 初始条件 Q~炮弹释放的烟雾总量 ? ~单位强度的点源函数 对任意t, C的等值面是球面 x2+y2+z2=R2; R??C? 仅当 t??, 对任意点(x,y,z), C?0 1）烟雾浓度 的变化规律 2）穿过烟雾光强的变化规律 光强的减少与烟雾浓度成正比 3）仪器灵敏度与烟雾明暗界限 烟雾浓度连续变化 烟雾中光强连续变化 仪器 z -? ? 设光源在z=-?, 仪器在z=?,则观测到的明暗界限为 不透光区域有扩大、缩小、消失的过程 穿过烟雾进入仪器的光线只有明暗之分，明暗界限由仪器灵敏度决定。 ~不透光区域边界 4）不透光区域边界的变化规律 对任意t, 不透光区域边界是圆周 不透光区域边界半径 r(t) rm 0 t1 t2 t 结果分析 观测到不透光区域边界达到最大的时刻t1，可以预报烟雾消失的时刻t2 1.6 万有引力定律的发现 背景 航海业发展 天文观测精确 “地心说”动摇 哥白尼：“日心说” 伽里略：落体运动 开普勒：行星运动三定律 变速运动的计算方法 牛顿：一切运动有力学原因 牛顿运动三定律 牛顿：研究变速运动，发明微积分（流数法） 开普勒三定律 牛顿运动第二定律 万有引力定律 《自然科学之数学原理》(1687) 模型假设 极坐标系 (r,?) 太阳 (0,0) 1. 行星轨道 a~长半轴, b~短半轴, e~离心率 3. 行星运行周期 T 行星位置：向径 O (太阳) ? P (行星) r 2. 单位时间 扫过面积为常数 A m ~ 行星质量 ? ~ 绝对常数 4. 行星运行受力 模型建立 O (太阳) ? P (行星) r 向径 的基向量 模型建立 万有引力定律 需证明 4A2/p =kM （与哪一颗行星无关） A~单位时间 扫过面积 O (太阳) ? P (行星) r (习题) l p / / 2 2 = p A 模型的解析求解 Logistic 模型 令 i=1/2 得 结果分析 1/2 tm i i0 1 0 t (1) ? (日接触率)? ? tm? Logistic 模型 病人可以治愈！ ?! 失败原因： 模型3—第二次改进 3）病人治愈成为健康人，病人每天治愈的比例为?。 增加假设 SIS 模型 4）健康人对传染病无免疫性，可再次被感染。 ? ~日治愈率 建模 ? ~ 日接触率 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 3）假设传染病有免疫性——病人治愈后即移出感染系统，称为移出者. 病人、健康人和移出者的比例分别为 SIR模型 改变假设 4）假设总人数N不变： 5）病人的日接触率? , 日治愈率?, 接触数 ? = ? / ? 建模 需建立 的两个方程 模型4--第三次改进 SIR模型 无法求出 的解析解 在相平面 上研究解的性质 模型4 消去dt SIR模型 相轨线 的定义域 相轨线 1 1 s i 0 D 在D内作相轨线 的图形，进行分析 s i 1 0 1 D 模型4 SIR模型 相轨线 及其分析 传染病蔓延 传染病不蔓延 s(t)单调减?相轨线的方向 P1 s0 im P1: i(t)先升后降至0 P2: