五、数学建模--微分方程模型.pptVIP

1）烟雾浓度C(x,y,z,t)的变化规律 初始条件 2）穿过烟雾光强的变化规律 I(l)~沿l方向的光强 C(l)~沿l方向的烟雾浓度 光强的减少与烟雾浓度成正比 l方向 3）仪器灵敏度与烟雾明暗界线 烟雾浓度连续变化 烟雾中光强连续变化 穿过烟雾进入仪器的光线只有明暗之分，明暗界线由仪器灵敏度确定 不透光区域有扩大、缩小、消失的过程 ~不透光区 域边界线 仪器 4）不透光区域边界的变化规律 不透光区域边界半径 结果分析 观察到不透光区域边界达到最大的时刻为t1,可以预报烟雾消失的时刻为t2 rm t1 t2 t r(t) 5.8万有引力定律的发现 背景 航海业发展 天文观测精确 “地心说”动摇 哥白尼“日心说” 伽里略：落体运动 开普勒：行星运动三定律 变速运动的计算方法 牛顿：一切运动有力学原因 牛顿运动三定律 牛顿：研究变速运动，发明微积分（流数法） 开普勒三定律 牛顿运动第二定律 万有引力定律 《自然科学之数学原理》(1687) 模型假设 极坐标系(r,?) 太阳 (0,0) 行星位置： O (太阳) P (行星) 1.行星轨道 a—长半轴，b—短半轴，e—离心率。 模型建立 O (太阳) ? P (行星) r 向径 的基向量 模型建立 万有引力定律 需证明 4A2/p =kM （与哪一颗行星无关） A~单位时间 扫过面积 O (太阳) ? P (行星) r (习题) l p / / 2 2 = p A 正规战争模型 为求战争的结局，不求x(t), y(t)，而在相平面上讨论x和y之间关系 平方律模型 游击战争模型 双方都用游击部队作战 甲方的战斗减员率还随着甲方的兵力增加而增加 f (x,y) = -cxy，c~乙方士兵杀伤力 忽略非战斗减员 假设没有增援 游击战争模型 线性律模型 混合战争模型 甲方用游击部队，乙方为正规部队 n 0 乙方胜 乙方必须10倍于甲方的兵力 5.4、药物在体内分布与排除 药物进入机体形成血药浓度(单位体积血液药物含量) 血药浓度必须保持在一定范围——给药方案设计 研究药物在体内吸收、分布和排除的过程——药物动力学 建立房室模型——药物动力学的基本步骤 房室——机体一部分，药物在一个房室内均匀分布(血药浓度)为常数，在房室间按一定规律转移 本节讨论二室模型——中心室(心、肺、肾等)和周边室(四肢、肌肉等) 模型假设 中心室(1)和周边室(2)，容积不变。 药物从体外进入中心室，在二室间相互转移，从中心室排除体外。 药物在房室间转移率及向体外排除率于该室的血药浓度成正比。 模型建立 xi(t)~药量ci(t)~浓度Vi~容积 i =1,2 中心室 周边室 给药 排除 模型建立 线性常系数非齐次方程 对应非其次方程组通解 其中 几种常见的给药方式 给药速率f0(t)和初始条件的讨论 1、快速静脉注射 t=0瞬间注入剂量D0的药物进入中心室，血药浓度立即为D0/V1 2、恒速静脉滴注 0≤t ≤ T药物以速率k0进入中心室 3、口服或肌肉注射 相当于药物(剂量D0)先进入吸收室，吸收后进入中心室 吸收室 中心室 参数估计 各种给药方式下，c1(t), c2(t)取决于参数k21, k12, k13, V1, V2等 t=0快速静脉注射D0，在ti (i =1,2,…,n) 时测得c1(ti) 参数估计 进入中心室得药物全部排除 5.5、香烟过滤嘴模型 问题 过滤嘴的作用与材料、长度的关系 建立吸烟过程的数学模型 分析吸入的毒物量与哪些因素有关，定量关系如何？ 模型分析 分析吸烟时毒物进入人体的过程 设想一个机器人在典型环境下吸烟，其方式和外界环境不变 模型假设 1）l1 ~ 烟草长度，l2 ~ 过滤嘴长度， l= l1 + l2毒物 M均匀分布，密度w0=M/l1 2）点燃处毒物随烟雾进入空气和沿香烟穿行的数量比是a’:a，a+a’=1 3）未点燃的烟草和过滤嘴对随烟雾穿行的毒物吸收率分别是b和β 4）烟雾随香烟穿行速度是常数v，香烟燃烧的速度为u，且vu 定性分析 Q ~ 吸一支烟毒物进入人体的总量 模型建立 t=0, x=0, 点燃香烟 q(x,t) ~毒物流量 w(x,t) ~毒物密度 w(x,0) =w0 1)、求q(x,0)=q(x) 1)、求q(x,0)=q(x) 2)、求q(l,t) t时刻，香烟燃至x=ut H(t)=uw(ut,t) 3)、求w(ut,t) 4)、计算Q 结果分析 1)、Q与a, M成正比，aM是毒物集中在x=l处的吸入量 4)、与另一支不带过滤嘴香烟比较。w0, b, a, v, l 均相同，且吸至x=l1处扔掉。 带过滤嘴 不带过滤嘴 5.6、按年龄分布的人口模型 年龄分布对人口预测的重要性 只考虑自然出生与死亡，不考虑迁移 人