微分方程建模的若干问题讲稿下.pptVIP

x t x = s ( t ) 0 s0 L 这个问题的 适定性 也已获得证明 。 于是，水、冰两相区域上任意点 x 处的温度 u ( x , t ) 随 时间 t 演变的模型为 （3） 细胞体内氧气的扩散与吸收问题 细胞体内氧气的会向周边 扩散 ，在 扩散 的同时，细胞 体也在 吸收 氧气以维持生命 ；如果细胞得不到氧气的 供给将会死亡。建立一个描绘该 扩散 — 吸收 过程的数 数学模型 。 为简单计，以下只考虑一个一维细胞体模型。 建模假设： （1）假定氧气在细胞体中从氧气浓度大的左边 扩散 至 浓度小的右边；在扩散中，扩散流量 q 的大小与 左、右两点的氧气浓度 c 的差成正比；即： （2）假定任何时刻，每单位立方体的细胞体 吸收 氧气 的速度为一常数 D ； （3）某一时刻起，断绝氧气供给；缺乏氧气的细胞体 即行死亡，不再参与氧气扩散过程 。 （ k 为扩散系数 ） 细胞体末端　　 氧气　　 考虑细胞体在位置 x 处、长为 Δx 的一段细胞上扩散 和吸收氧气情况。 在 Δt 时段内，经扩散进入这段细胞内的氧气数量是 : 经扩散流出这段细胞内的氧气数量是 : 这段细胞内氧气的变化量是： 这段细胞氧气的吸收量是： 进入量、流出量、变化量和吸收量之间应有关系： 根据假设（1）， 氧气扩散、吸收方程 0 x t s0 在细胞体左端，在 t = 0 起断绝氧气输入，故有： 在细胞体右末端 x = s 处，始终有条件： 随着氧气的缺乏， 右末端的细胞逐渐 死亡，故有末端的 位置随时间而变动， 形成一条 自由边界 ： x = s ( t ) . 氧气扩散、吸收问题 ： 寻求未知函数对：{ c ( x , t ) ，s ( t ) } ， 使得它们满足： 在初边至充分光滑情况下，这个问题的 适定性 也可证明。 事实上，若该问题的充分 光滑解为 c ( x , t ) ， 令 u ( x , t ) = c t ( x , t ) , 则有 已知 u ( x , t ) 后 ， 即可得到 ： 0 x t s0 四. 用差分法求解微分方程的数值解 （向前偏心差分）； （中心差分）； （向前偏心差分）； （中心差分）。 1. 用有限差分替代 一阶（偏）导数情况 2. 用有限差分替代 二阶（偏）导数情况 令 （ N 等分区间，步长 h 越小越精确）， 有限差分格式： 练习：分别取 h = 0.1 和 h = 0.05, 求解 （ 与精确解 u ( t ) = tg t – t 作一比较 ） I ) 常微分方程求解： 3. 用差分法求解微分方程数值解实例 令 中心差分五点格式： 其中 练习： (1) 已知 u(x,y) 满足以下 Laplace 方程 的边值问题 ： 求解： II) Poisson 方程 取网格 h = τ = 1/3, 求出 之值。 (2) 用正方形网格 h =τ并分别取 h = 1 和 h = 1/2 求解 Poisson方程的边值问题： III) 热传导方程 令 为网比，上式可以化为 “ 显示格式 ”： 鉴于收敛稳定性考虑，也可以在 t 方向采用另一种差商（向后差商）来离散化微分方程： 令 为网比，上式可以化为 “ 隐式格式 ” ： . 练习： 求解 在 t = 1 和 t = 2 上各节点 u 的近似值 。 IV) 弦振动方程 (2) F( r , t ) 为关于 r , t 的连续可微二元函数，记 (4) 任何时刻 t ，出生婴儿的人口分布密度（即单位时间 内出生的婴儿数）为已知函数