椭圆型方程差分方法.ppt

第三章 椭圆型方程的差分方法 3.1 正方形区域中的Laplace方程Dirichlet边值问题的差分模拟 3.2 Neumann边值问题的差分模拟 3.3 混合边值条件 3.4 非矩形区域 3.5 极坐标形式的差分格式 3.6 矩形区域上的Poisson方程的五点差分逼近的敛速分析 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 3.1 正方形区域中的Laplace方程Dirichlet边值问题的差分模拟 3.3 混合边值问题 3.4 非矩形区域 3.5 极坐标形式的差分格式 3.6 矩形区域上的Poisson方程的五点差分逼近的敛速分析 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 例 3.2 考虑Laplace方程第一边值问题 这里 。采用步长为1/4的正方形网格，差分公式为 结点编号如图3.6所示，则按前法所得方程组 AU=K的系数矩阵为 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 这是一个对称、不可约对角占优矩阵，对角元为正，因此它非奇异，且为对称正定阵。 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 例 3.3 考虑椭圆型方程第一边值问题 这时 采用步长为h=1/3的正方形网格，因此有四个网格内点(如图3.7所示)。 差分格式为 上式中 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 则得四阶线性方程组 AU=K 其中 系数矩阵 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 经具体计算，得 A为严格对角优势，非对称矩阵。 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 例 3.4 自伴线性椭圆型方程第一边值问题 用中心差商近似导数，则 差分方程为 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 即 令 格式可写为 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 网格区域如图3.7所示，则系数矩阵A为 它是对角优势、不可约对称矩阵。 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 例 3.5 考虑区域Ω中椭圆型方程 (3.41) 其中 。为了建立差分格式，在Ω中覆盖一正方形网格区域，步长为h，在区域内点(l,m)上，有 而混合偏导数 一开始用下式 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 代替。这里 。 在(l,m)点Taylor展开，则 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 由此，(3.42)为 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 如果 （3.43） 则 项消失，因此，选择 满足式(3.43)，于是逼近微分方程(3.41)的差分方程能具有截断误差阶 ，这时差分方程为 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 令 格式为 因此(1) 若 ，则选择 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 这时 差分方程为 （3.44.1） 显然，若 ，则式(3.44.1)中 都为正。 (2) 若 ，可令 ，则 3.7 一般二阶线性椭圆型方程差分逼近及其性质研究 差分方程为 (3.44.2) 如果 ，则格式(3.44.2)中 都为正。 上面就不同情况，在个结点上列出了差分方程，设结点按自然次序排列，联立它们得到线性代数方程组 AU=K A为不可约对角优势矩阵。 再用中心差商代替微商，就得出点(r0,θm)的差 分方程 或者 注： 此处没有给出边界条件，因此也没有边界处点 值的确定方法。 考虑泊松方程的第一边值问题 （3.31） 其中 中的网格点集合为 ，而 边界上网格点集合为 ，而 网格点 (记为 ) 上的五点差分逼近是 在边界网格点(l,m)上 满足边值条件 (3.33) 为了讨论差分方程解 与微分方程解 之间的逼 近程度，令 于是 写成 (3.