第14章SECTION4偏微分方程的数值解法.docVIP

§4 偏微分方程的数值解法 差分法 差分法是常用的一种数值解法.它是在微分方程中用差商代替偏导数，得到相应的差分方程，通过解差分方程得到微分方程解的近似值. 1. 网格与差商 图14.7 在平面 (x,y)上的一以S为边界的有界区域D上考虑定解问题.为了用差分法求解，分别作平行于x轴和y轴的直线族. (i,j=0,SYMBOL 177 \f Symbol1,SYMBOL 177 \f Symbol2,…,SYMBOL 177 \f Symboln) 作成一个正方形网格，这里h为事先指定的正数，称为步长；网格的交点称为节点，简记为(i,j).取一些与边界S接近的网格节点，用它们连成折线Sh，Sh所围成的区域记作Dh.称Dh内的节点为内节点，位于Sh上的节点称为边界节点（图14.7）.下面都在网格Dh + Sh上考虑问题：寻求各个节点上解的近似值.在边界节点上取与它最接近的边界点上的边值作为解的近似值，而在内节点上，用以下的差商代替偏导数： 注意， 1? 式中的差商称为向后差商，而称为向前差商，称为中心差商.也可用向前差商或中心差商代替一阶偏导数. 2? x轴与y轴也可分别采用不同的步长h,l，即用直线族 (i,j=0, ±1, ±2 ) 作一个矩形网格. 2. 椭圆型方程的差分方法 [五点格式] 考虑拉普拉斯方程的第一边值问题 式中SYMBOL 109 \f Symbol(x,y)为定义在D的边界S上的已知函数. 采用正方形网格，记u(xi,yj)=uij ，在节点(i,j)上分别用差商 代替，对应的差分方程为 （1） 或 即任一节点(i,j)上uij的值等于周围相邻节点上解的值的算术平均，这种形式的差分方程称为五点格式，在边界节点上取 （2） 式中(xi*,yj*)是与节点(i,j)最接近的S上的点.于是得到了以所有内节点上的uij值为未知量的若干个线性代数方程，由于每一个节点都可列出一个方程，所以未知量的个数与方程的个数都等于节点的总数，于是，可用通常的方法（如高斯消去法）解此线性代数方程组，但当步长不很大时，用高斯消去法将会遇到很大困难，可用下面介绍的其他方法求解. 若hSYMBOL 174 \f Symbol0时，差分方程的解收敛于微分方程的解，则称差分方程为收敛的. 在计算过程中，由于进行四则运算引起舍入误差，每一步计算的舍入误差都会影响以后的计算结果，如果这种影响所产生的计算偏差可以控制，而不至于随着计算次数的增加而无限增大，则称差分方程是稳定的. [迭代法解差分方程] 在五点格式的差分方程中，任意取一组初值{uij}，只要求它们在边界节点(i,j)上取以已知值SYMBOL 109 \f Symbol(xi*,yj*)，然后用逐次逼近法（也称迭代法）解五点格式： 逐次求出{uij(n)}.当(i+1,j)，(i－1,j)，(i,j－1)，(i,j+1)中有一点是边界节点时，每次迭代时，都要在这一点上取最接近的边界点的值.当nSYMBOL 174 \f Symbol∞时，uij(n)收敛于差分方程的解，因此n充分大时，{uij(n)}可作差分方程的近似解，迭代次数越多，近似解越接近差分方程的解. [用调节余数法求节点上解的近似值] 以差商代替Δu时，用节点(i+1,j)，(i－1,j)，(i,j+1)，(i,j－1)上u的近似值来表示u在节点(i,j)的值将产生的误差，称此误差为余数Rij，即 图14.8 设在(i,j)上给uij以改变量SYMBOL 100 \f Symboluij，从上式可见Rij将减少4SYMBOL 100 \f Symboluij，而其余含有u(xi,yj)的差分方程中的余数将增加SYMBOL 100 \f Symboluij，多次调整SYMBOL 100 \f Symboluij的值就可将余数调整到许可的有效数字的范围内，这样可获得各节点上u(x,y)的近似值.这种方法比较简单，特别在对称区域中计算更简捷. 例 求Δu=0在内节点A,B,C,D上解的近似值.设在边界节点1,2,3,4上分别取值为1,2,3,4（图14.8） 解 记u(A)=uA，点A,B,C,D的余数分别为 －4uA+ uB+ uc +5=RA uA－4 uB + uD+7=RB