椭圆方程数值解.pdfVIP

j. 椭圆方程数值解法 本章考虑椭圆微分方程数值解法。首先以二维二阶椭圆方程为例，给出矩 形网和三角网上的差分法。 然后以一维二阶椭圆方程为例， 简要描述有限元法的 基本思想。 J.1 矩形网上差分方程 考虑二维区域（区域 =连通的开集） G 上的二阶椭圆型偏微分方程第一边 值问题 uxx uyy Cux Du y Eu F x, y G （j.1 ） u x , y x , y 0 其中 C , D , E 是常数； E 0 ； F F x, y C G ； (x, y ) 是给定的光滑函数； 是 G 的边界； G G 。假设（ J.1 ）存在光滑的唯一解。 考虑一种简单情形，即求解区域 G 是矩形区域，并且其四个边与相应坐标 轴平行。令 h 和 h 分别为 x 和 y 方向的步长，用平行于坐标轴的直线段分割区域 1 2 G ，构造矩形网格： G h 为网格内点节点集合， h 为网格边界节点集合， G h G h h 。 对于内点 x , y G ，用如下的差分方程逼近微分方程（ J.1 ）： i j h (J.2) u 2u u u 2u u u u u u i 1, j ij i 1, j i , j 1 ij i , j 1 i 1, j i 1,j i , j 1 i , j 1 C D Eu F 2 2 2 ij ij h h 2h 2h 1 2 1 2 ( , ) 其中 Fij F xi y j 。（J.2 ）通常称为五点差分格式 。 方程（ J.2 ）可以整理改写为 （J.3 ） ai 1, j ui 1, j + ai 1, j ui 1, j + ai , j 1 ui , j 1 + ai , j 1 ui , j 1 +ai , j ui , j Fij 对每一内点 x , y 都可以列出这样一个方程。 方程中遇到边界点时， 注意到边界 i j 点上函数值 u 已知，将相应的项挪到右端去。最后得到以 u 的内点近似值为未知 数的线性方程组。这个方程组是稀疏的，并且当 h 和 h 足够小时是对角占优的。