1微分方程建模.pptVIP

微分方程求解与建模 华中农业大学数学建模基地系列课件 内容提要 问题背景和实验目的 Euler 折线法 Runge-Kutta 方法 Matlab 解初值问题 数值求解举例 内容提要 微分方程模型 冷却或加热模型 车间空气的清洁 赝品的鉴定 药物在体内的分布 人口问题的偏微分方程模型 交通流问题 问题背景和实验目的 Euler 折线法 Euler折线法 Euler 折线法 Euler 折线法源程序 Euler折线法举例（续） Runge-Kutta 方法 Runge-Kutta 方法 四阶 R-K 方法源程序 Runge-Kutta 方法 Euler 法与 R-K法误差比较 Matlab 解初值问题 dsolve 求解析解 dsolve 的使用 dsolve 举例 dsolve 举例 dsolve 举例 Matlab函数数值求解 Matlab提供的ODE求解器 参数说明 数值求解举例 数值求解举例 Matlab 求解微分方程小结 微分方程模型 涉及“改变”、“变化”、“增加”、“减少”、“衰变”、“边际”、“速度”、 “运动”、“追赶”、“逃跑”、、、等等词语的确定性连续问题。 1、寻找改变量 一般说来微分方程问题都遵循这样的文字等式 变化率（微商）=单位增加量--单位减少量 等式通常是利用已有的原则或定律。 3、用微元法建立微分方程； 4、确定微分方程的定解条件（初边值条件）； 5、求解或讨论方程（数值解或定性理论）； 6、模型和结果的讨论与分析。 例1 一个较热的物体置于室温为18℃的房间内，该物体最初的温度是60 ℃ ，3分钟以后降到50 ℃ .想知道它的温度降到30 ℃需要多少时间？10分钟以后它的温度是多少？ 例3 赝品的鉴定 例3 赝品的鉴定 例3 赝品的鉴定 例3 赝品的鉴定 例3 赝品的鉴定 4 药物在体内的分布 4 药物在体内的分布 4 药物在体内的分布 5 人口问题的偏微分方程模型 5 人口问题的偏微分方程模型 6 交通流问题 6 交通流问题 c、微分方程建模的基本规则 2、对问题中的特征进行数学刻画 微分方程模型 微分方程模型 例1 冷却或加热模型 牛顿冷却（加热）定律：将温度为T的物体 放入处于常温 m 的介质中时，T的变化速率 正比于T与周围介质的温度差. 分析：假设房间足够大，放入温度较低或较 高的物体时，室内温度基本不受影响，即室温 分布均衡,保持为m，采用牛顿冷却定律是一个 相当好的近似. 建立模型：设物体在冷却过程中的温度为T(t)，t≥0， 例1 冷却或加热模型 “T的变化速率正比于T与周围介质的温度差” 翻译为 数学语言 建立微分方程 其中参数k 0，m=18. 求得一般解为 例1 冷却或加热模型 ln(T－m)=－k t+c, 代入条件，求得c=42 ， , 最后得 T(t)=18+42 , t ≥0. 结果 ：T(10)=18+42 =34.97℃， 该物体温度降至30℃ 需要13.82分钟. 例1 冷却或加热模型 例2 车间空气的清洁 问题：已知一个车间体积为V立方米，其中有一台机器每分钟能产生r立方米的二氧化碳（CO2），为清洁车间里的空气，降低空气中的CO2含量，用一台风量为K立方米/分钟的鼓风机通入含CO2为m%的新鲜空气来降低车间里的空气的CO2含量。假定通入的新鲜空气能与原空气迅速地均匀混合，并以相同的风量排出车间。又设鼓风机开始工作时车间空气中含x0%的CO2.问经过t时刻后，车间空气中含百分之几的CO2?最多能把车间空气中CO2的百分比降到多少？ 设t时刻（单位为分钟）车间每立方米空气含CO2的百分比为x(t),考虑时间区间 并利用质量守恒定律： 内车间空气含CO2量的“增加”等于 时间内进入的新鲜空气中含CO2的量加上机器产生的CO2的量减去排出空气中CO2的量。用数学公式表示出来就是 分析和建模 例2 车间空气的清洁 于是，令 得 其中， 解为 这就是t时刻空气中含CO2的百分比。 通常 否则含CO2的量只会增加。 例2 车间空气的清洁 令 得 这表明车间空气中含CO2的量最多只能降到 讨论：如果设V=10000立方米，r=0.3立方米/分钟，K=1500立方米/分钟，m=0.04%，x0=0.12%。试问： （1）需多少分钟后，车间空气中含CO2的百分比低于0.08%? （2）车间空气中含CO2的百分比最多只能降到多少？ 例2 车间空气的清洁 在第二次世界大战比利时解放以后，荷兰野战军保安机关开