M微分方程建模实例概述.ppt

*/33 * 4-1 微分方程的基本理论 4-2 微分方程建模实例 4-3 我国的人口预测模型 Differential equation model and its application 一、简单建模实例 1.1 溶液浓度模型 1.2 时间估计模型 三、小结 2.1 最佳销售时机模型 2.2 传染病模型 1.3 放射性元素衰变模型 二、综合建模实例 溶液浓度问题是工农业生产和治理环境污染中经常要碰到的问题。此类问题通常都可以描述为如下的实验室模型： 1.1 溶液浓度模型 如图所示，一个容器有一个入口，一个出口，里面盛满了某种溶液，如果从入口以不变的速率向容器内注入一定浓度的相同溶液（或清水），搅拌均匀后以同样的速率从出口排出，假设搅拌是在瞬间完成的，那么容器内溶液浓度的变化规律是怎样的呢？ 例1 已知容器内盛有1000公斤的清水，若以每分钟5公斤的速率注入浓度为0.2的盐水且不停地搅拌, 并以同样的速率排出搅拌后的盐水, 那么经过多少时间能使容器内的含盐量达到100公斤？ 基本 假设 假设搅拌是在瞬间完成的，即容器中液体浓度任一时刻都是均匀变化的。 t时刻容液浓度为: 符号说明 y(t)—t时刻容器内的含盐量； y(0)—初始时刻容器内的含盐量；t—任一时刻。 模型建立 y(t)/1000即0.001y(t)。 [t,t+dt]时间间隔内, 进盐量为: 出盐量为: 含盐量的微元为 这是一个一阶线性非齐次微分方程，易求得该方程满足初始条件y(0)=0的特解为 这是容器内含盐量y随t的变化规律. 将y=100代入,可求得 即经过约2小时18分37秒可使容器内的含盐量达到100公斤。 模型求解 1.2 时间估计模型 在凌晨1时警察发现一具尸体，测得尸体温度是29℃，当时环境的温度是21℃，一小时后尸体温度下降到27℃，如果人的正常体温是37℃。请帮助警察估计死者的死亡时间。 问题提出 根据牛顿冷却定律：将温度为T的物体放入温度为T0的介质中，则该物体的温度T的变化速率正比于该物体与周围介质的温度差T-T0，即 设该名死者已经死亡t1小时，于是根据题意得到如下的微分方程边值问题： 模型建立 模型求解 方程(1)的通解为 把环境温度T0=21及边界条件(2)代入通解，得 于是可知该死者死于2.409小时前, 比零晨1时早2.409小时，即前一天夜晚10:35时左右。 该边值问题的解为 通过Matlab作出其 曲线分布如下所示。 上述模型的建立是假设环境的温度恒定为T0=21℃不变，若死者的死亡时间较短则上述模型可以很好的推算出死者的死亡时间；但是若死亡时间较长，显然环境的温度应该是有变化的。从而上述模型可以改进为非自治的情形 。 模型改进 其中T0(t)是描述环境温度随时间变化的函数. 方程(3)是一阶线性的, 其通解为 在已知T0(t)的情况下，类似得方法可以求得t1值。T0(t)可以是根据经验获得的函数，也可以是通过环境温度的变化数据来拟合得到。 1.3 放射性元素衰变模型 1972年8月，湖南长沙出土了马王堆一号墓（注：出土时因墓中女尸历经千年而未腐曾经轰动世界）。经测定，出土的木炭标本中14C的平均原子蜕变速度为29.78次/分，而新砍伐烧成的木炭中14C的平均原子蜕变速度为38.37次/分；如果14C 的半衰期取为5568年（注:14C的半衰期在各种资料中说法不一，分别有5568年、5580年和5730年不等），那么，怎样才能根据以上数据确定这座墓葬的大致年代呢？ 问题 提出 放射性元素的衰变规律 衰变：放射性元素的质量随时间的推移而逐渐 减少(负增长)，这种现象称为衰变。 设放射性元素t时刻的质量为m(t),则其衰变速度即 于是可得 是比例常数，可由该元素的半衰期确定， 负号表明放射性元素的质量m是随时刻t递减的。 若m(0)=m0,则可求方程的特解： 指数递减模型 由物理学定律知，放射性元素任一时刻的衰变速度与该时刻放射性元素的质量成正比。据此原理，我们也可以通过微分方程研究放射性元素衰变的规律。 放射性元素任一时刻的衰变速度与该时刻放射性元素的质量成正比，且比例系数为常数λ。 符号说明 t—任意时刻； t0—初始时刻； △t—时间增量； m(t)—t时的元素质量；m0—t0时的元素质量； λ——比例常数。 模型建立 基本 假设 由基本假设有 放射性元素衰变模型 解微分方程，得 元素衰变微分方程 指数递减模型 初始条件为 14C 解之得 为应用, 令半衰期为T, 故有 故指数递减模型又可表示为 故从初始质量m0衰减到m(