椭圆方程数值解简析.docVIP

j. 椭圆方程数值解法 本章考虑椭圆微分方程数值解法。首先以二维二阶椭圆方程为例，给出矩形网和三角网上的差分法。然后以一维二阶椭圆方程为例，简要描述有限元法的基本思想。 J.1 矩形网上差分方程 考虑二维区域（区域=连通的开集）上的二阶椭圆型偏微分方程第一边值问题 （j.1） 其中,是常数；；；是给定的光滑函数；是的边界；。假设（J.1）存在光滑的唯一解。 考虑一种简单情形，即求解区域是矩形区域，并且其四个边与相应坐标轴平行。令和分别为和方向的步长，用平行于坐标轴的直线段分割区域，构造矩形网格： 为网格内点节点集合，为网格边界节点集合，。 对于内点，用如下的差分方程逼近微分方程（J.1）： (J.2) 其中。（J.2）通常称为五点差分格式。 方程（J.2）可以整理改写为 （J.3） ++++ 对每一内点都可以列出这样一个方程。方程中遇到边界点时，注意到边界点上函数值已知，将相应的项挪到右端去。最后得到以的内点近似值为未知数的线性方程组。这个方程组是稀疏的，并且当和足够小时是对角占优的。 用（J.1）的真解在网点上的值、等等分别替换（J.2）中的、等等，然后在点处作Tailor展开，便知差分方程（J.2）逼近微分方程（J.1）的截断误差阶为。另外可以证明，五点差分格式的收敛阶为，并且关于右端和初值都是稳定的。 矩形网格差分格式的优点是计算公式简单直观