数学建模讲座-微分方程模型.ppt

一般模型 每方战斗减员率取决于双方的兵力和战斗力 每方非战斗减员率与本方兵力成正比 甲乙双方的增援率为u(t), v(t) f, g 取决于战争类型 x(t) ~甲方兵力，y(t) ~乙方兵力 模型假设 模型 正规战争模型 甲方战斗减员率只取决于乙方的兵力和战斗力 双方均以正规部队作战 忽略非战斗减员 假设没有增援 f(x, y)=?ay, a ~ 乙方每个士兵的杀伤率 a=ry py, ry ~射击率， py ~命中率 0 正规战争模型 为判断战争的结局，不求x(t), y(t)而在相平面上讨论 x 与 y 的关系 平方律 模型 乙方胜 游击战争模型 双方都用游击部队作战 甲方战斗减员率还随着甲方兵力的增加而增加 忽略非战斗减员 假设没有增援 f(x, y)=?cxy, c~ 乙方每个士兵的杀伤率 c = ry py ry~射击率 py ~命中率 py=sry /sx sx ~ 甲方活动面积 sry ~ 乙方射击有效面积 0 游击战争模型 线性律 模型 0 混合战争模型 甲方为游击部队，乙方为正规部队 乙方必须10倍于甲方的兵力 设 x0=100, rx/ry=1/2, px=0.1, sx=1(km2), sry=1(m2) 4 药物在体内的分布与排除 药物进入机体形成血药浓度(单位体积血液的药物量) 血药浓度需保持在一定范围内——给药方案设计 药物在体内吸收、分布和排除过程 ——药物动力学 建立房室模型——药物动力学的基本步骤 房室——机体的一部分，药物在一个房室内均匀分布(血药浓度为常数)，在房室间按一定规律转移 本节讨论二室模型——中心室(心、肺、肾等)和周边室(四肢、肌肉等) 中心室 周边室 给药 排除 模型假设 中心室(1)和周边室(2),容积不变 药物在房室间转移速率及向体外排除速率，与该室血药浓度成正比 药物从体外进入中心室，在二室间相互转移,从中心室排出体外 模型建立 线性常系数非齐次方程 对应齐次方程通解 模型建立 几种常见的给药方式 1.快速静脉注射 t=0 瞬时注射剂量D0的药物进入中心室,血药浓度立即为D0/V1 给药速率 f0(t) 和初始条件 2.恒速静脉滴注 t T, c1(t)和 c2(t)按指数规律趋于零 药物以速率k0进入中心室 0 T t ￡ ￡ 吸收室 中心室 3.口服或肌肉注射 相当于药物( 剂量D0)先进入吸收室，吸收后进入中心室 吸收室药量x0(t) 参数估计 各种给药方式下的 c1(t), c2(t) 取决于参数k12, k21, k13, V1,V2 t=0快速静脉注射D0 ,在ti(i=1,2,?n)测得c1(ti) 由较大的 用最小二乘法定A,? 由较小的 用最小二乘法定B,? 参数估计 进入中心室的药物全部排除 * 解法指令 解题类型 特 点 适合场合 ode45 非刚性 采用4、5阶Runge－Kutta法 大多数场合的首选算法 ode23 非刚性 采用Adams算法 较低精度（10－3）场合 ode113 非刚性 多步法；采用Adams算法；高低精度均可（10－3～10－6） ode45计算时间太长时取代ode45 ode23t 适度刚性 采用梯形法则算法 适度刚性 ode15s 刚性 多步法；采用2阶Rosenbrock算式，精度中等 当ode45失败时使用；或存在质量矩阵时 ode23s 刚性 一步法；采用2阶Rosenbrock算式，低精度 低精度时，比ode15s有效；或存在质量矩阵时 ode23tb 刚性 采用梯形法则－反向数值微分两阶段算法，低精度 低精度时，比ode15s有效；或存在质量矩阵时 各种solver 解算指令的特点 五、 微分方程模型在 模型分析 中的主要问题 之 —— 稳定性分析 用微分方程方法建立的动态模型问题 模型分析 中的一个 重要问题是：当时间充分长后 ，动态过程的 变化趋势 是什么？ 微分方程模型中 , 方程 ( 组 ) + 初始条件 → 解 初始条件的作用在于确定解, 它的微小变化会产生不同的 解，换言之，对解的发展性态变化 , 往往具有影响作用 . 问题是这种对解的发展性态的影响作用是 长期存在 的 , 还是当时间充分大以后 , 影响作用会 “消逝 ” ? （1）微分方程模型的稳定性及其实际意义 有时候 , 初始条件的微小变化会导致解的性态随时间变 大后 , 产生显著的差异 , 这时称 系统是不稳定 的 ; 有时候 , 初始条件变化导致解的性态差异会随时间变大 后而消失 , 这时称该 系统是稳定 的. 在实际问题中, 初始状