第10章-偏微分方程数值解-1.pptVIP

* 第十章 偏微分方程数值解 一、典型的偏微分方程介绍 1. 椭圆型方程 Laplace 方程 Poisson方程 2. 抛物型方程 热传导方程 其中a是常数。它表示长度为L的细杆内，物体温度分布的规律 土壤水运动方程： 溶质运移方程： （水流稳态） （瞬态） 3．双曲型方程 波动方程 它表示长度为L的弦振动的规律。 二、定解问题 决定方程唯一解所必须给定的初始条件和边界条件 叫做定解条件 边界条件 初始条件 计算机只能作有限次的加、减、乘、除运算， 它既不能求导数，更不能解偏微分方程。如果想 在计算机上求得微分方程数值解，它的主要做法 是把偏微分方程中所有的偏导数分别用差商代替， 从而得到一个代数方程组——差分方程组，然后 对差分方程求解，并以所求的解作为偏微分方程 数值解。 10.1 差分法简介 对区域进行剖分，用网格点来代替连续区域， 因此差分法亦称“网格法”。 0 x y 把整体分割成若干个单元来处理问题的方法在 数学上称为“离散化方法” 在结点上采用离散化方法（数值微分、数值积分、 泰勒展开等）将微分方程的初边值问题化成关于 离散变量的相应问题，这个相应问题的解就是方程 在点xi上的数值解f(x)，或在点（xi , ti）上的数值解 U( xi , ti)。一般来说，不同的离散化导致不同的方法。 例：取一边长为1的正方形均匀薄板，上下侧面绝热，四周保持恒温，求板内各点的稳定温定分布。 u=0 u=0 u=0 0 x y Laplace 方程第一边值问题 记 u在这些点满足方程 得到u (i, k)的近似ui,k，所满足的线性代数方程组： 其中 用迭代法来解方程组 简单迭代法 高斯—赛德尔迭代法 i=4 i=3 i=2 i=1 i=0 k=0 0 0 0 0 0 k=1 0.707 0 0.354 0 0 k=2 1 0.75 0 0.25 0 k=3 0.707 0 0.354 0 0 k=4 0 0 0 0 0 表10.1 0 0 0 0 0 0.707 0.453 0.354 0.151 0 1 0.75 0.427 0.25 0 0.707 0.453 0.354 0.151 0 0 0 0 0 0 表10.2 i=0 i=1 i=2 i=3 i=4 k=0 k=1 k=2 k=3 k=4 0 0 0 0 0 0.707 0.453 0.258 0.151 0 1 0.583 0.427 0.182 0 0.707 0.453 0.258 0.151 0 0 0 0 0 0 表10.3 i=0 i=1 i=2 i=3 i=4 k=0 k=1 k=2 k=3 k=4 0 0 0 0 0 0.707 0.453 0.258 0.151 0 1 0.573 0.386 0.182 0 0.707 0.381 0.243 0.134 0 0 0 0 0 0 表10.4 i=0 i=1 i=2 i=3 i=4 k=0 k=1 k=2 k=3 k=4 用差分法解偏微分方程需要考虑三个问题： 1．选用网格，将微分方程离散化为差分方程。 2．当网格步长h ?0时差分方程的准确解是否 收敛于微分方程的解？ 3．如何解相应的代数方程组？ 10.2 椭圆型方程的差分解法 椭圆型方程最简单的典型问题就是拉普拉斯方程 泊松方程 考虑泊松方程第一边值问题： (一) 矩形网格 设?为xy平面上一有界区域，??为其边界， 是分段光滑曲线。 0 x y 正则内点 非正则内点 边界点 （二）五点差分格式 现在假设（i，k）为正则内点。沿着x，y轴方向分别用 二阶中心差商代替uxx，uyy，则得 若以uh，fh表示网函数，记 *