数学建模,第三章-微分方程模型.ppt

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) 黑龙江科技学院 数 学 建 模 理学院 设待求解的微分方程问题为：  （1） 其中f适当光滑，对y满足Lipschitz（利普希茨）条 件，即存在L使 ，则方程的 解存在且唯一。 常微分方程初值问题的数值解法是近似计算中很重要的部分。重点讨论一阶常微分方程数值解。 黑龙江科技学院 数 学 建 模 理学院 3.6 微分方程数值解 一般情况下，求解微分方程的解析解非常困难甚至有时根本无法求出。对于微分方程而言，更重要的是求数值解，也就是在一系列离散点 上求解析解y=y(x)在点 处的值 的近似值 。通常情况下，为计算方便取等步长h，即 黑龙江科技学院 数 学 建 模 理学院 欧拉方法（折线法）它的基本思想是在小区间 上用差商 代替导数 ，而方程右端函数 中的x在小区间 的端点上取值，得到方程的近似表达式，称为欧拉公式。 向前欧拉公式 黑龙江科技学院 数 学 建 模 理学院 向后欧拉公式 梯形公式 方法：把迭代过程分为两步: (1)由向前欧拉公式算出 的预测值 ；再把它代入梯形公式右端，作为校正，即 梯形公式的预报校正格式 或者 黑龙江科技学院 数 学 建 模 理学院 Runge—kutta（龙格—库塔）法 二阶Runge—kutta公式 基本思想:在 内多预报几个点的斜率值，然后把他们取加权平均作为平均斜率，则可构造出具有高精度的计算格式。 其中 黑龙江科技学院 数 学 建 模 理学院  其中待定系数λi，αi， βi共13个，由于计算过 于复杂所以我们直接给 出一组λi，αi，βi的 值得到 四阶经典Runge—kutta公式 黑龙江科技学院 数 学 建 模 理学院 3.7 建模举例（减肥模型） 试建立一个减肥的数学模型，探讨如下几个问题： 1、建立体重随时间的变化规律，进而提出减肥的一些措施； 2、给出衡量减肥效果的指标； 3、计算出减肥的临界状态（即体重不能再减小的状态，若低于 该指标，就影响身体健康，甚至危及生命）。 想一想影响体重变化的因素有哪些？ 因此我们假