椭圆型方程差分方法.pptVIP

第三章 椭圆型方程的差分方法;(3.1); 第二边值问题，或称Neumann问题;3.1 正方形区域中的Laplace方程Dirichlet边值问题的差分模拟;因此Laplace方程的五点差分格式为;如图3.1所示; 单位正方形中的内部结点上的 个线性方程(3.8)写成矩阵形式为 AU=K （3.9）;3.2 Neumann边值问题的差分模拟;且除了一个任意常数外，解唯一。因为容易看到，如 果u(x,y)是式(3.10)的解，于是， u(x,y)+C(C是一个任 意常数)也是其解。为了唯一性，需要规定u(x,y)在区域 中某一点上的值。;Neumann边值问题的差分模拟;在x=0上的导数边值条件的差分模拟为;同理，在x=1,y=0,y=1时分别有; 由此，正方形 区域的 个结点上差分 方程解 满足线性方程组;方程组(3.18)中的向量U和g由以下给出：;例 3.1 在单位正方形区域Ω上解Laplace方程Neumann 问题;(3.19);3.3 混合边值问题;在求解区域Ω内由逼近Laplace方程的五点差分公式;对于在x=0上的结点(0,mh)，应用边值条件;在原点(0,0)上，两边值条件相遇，则; 这样在整个计算区域及相应边界网格点上建立了差 分方程：;其中矩阵A为 阶对称方阵。;3.4 非矩形区域;由解u在结点(l,m)上的Taylor展开可得;因此Laplace方程的五点差分格式为;3.5 极坐标形式的差分格式;这样可以将所求解区域映射为;即;且设;对; 再用中心差商代替微商，就得出点(r0,θm)的差 分方程;3.6 矩形区域上的Poisson方程的五点差分逼近的敛速分析; 在边界网格点(l,m)上 满足边值条件;则; 定理 3.1 如果是定义在矩形区域;，; 应用定理3.1的结论，则差分方程解与微分方程解;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;在网格结点;其中; 显然，如果适当选取h，使其满足; ;;记;根据①具有严格对角优势的矩阵是非奇异的， ②具有对角优势的不可约矩阵是非奇异的， ③具有对角优势和正的对角元的不可约对称矩阵 的特征值全是正数，即为对称正定阵。;;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;3.7 一般二阶线性椭圆型方程差分逼近及其性质研究;