五点差分格式求解第一边值问题.docxVIP

五点差分格式求解第一边值问题一五点差分格式构造的介绍考虑Poisson方程：， (x,y), 1.1G是xy平面上一有界区域，其边界为分段光滑曲线，在满足边值条件：1.2为解决此问题，做矩形网的差分格式，取定x轴和y轴方向的步长h1和h2，沿x 、y方向分别用二阶中心差商代替微商，则得+]=1.3其中表示节点 ()上的网函数由于差分格式1.3只出现u在(i,j)及其四个邻点上的值，称为五点差分格式，如图示：二、模型问题用五点差分格式求解单位正方形区域的Poisson方程第一边值问题：其中为正方形区域的边界 y y j=4u9 j=3 u6 j=2 u3 j=1 j=0 i=0i=1 i=2 i=3 i=4 x图 1 图 2 特别取正方形网格h1=h2=h,利用Yaylor展式，带入上式后，再用uij代替u(xi,yj), 并略去误差项得= 1.4现假定，将1.4式写成矩阵形式即为U+U=, 1