网站大量收购独家精品文档,联系QQ:2885784924

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

  1. 1、本文档共77页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数学建模讲座-微分方程模型

一般模型 每方战斗减员率取决于双方的兵力和战斗力 每方非战斗减员率与本方兵力成正比 甲乙双方的增援率为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)微分方程模型的稳定性及其实际意义 有时候 , 初始条件的微小变化会导致解的性态随时间变 大后 , 产生显著的差异 , 这时称 系统是不稳定 的 ; 有时候 , 初始条件变化导致解的性态差异会随时间变大 后而消失 , 这时称该 系统是稳定 的. 在实际问题中, 初始状

文档评论(0)

yaocen + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档