椭圆型方程课件.pptVIP

2021/3/23 * 3.2 边值条件的处理 1 讨论第一边值条件： 当 时，取 当 是非正则内点，如同正则内点，建立不等距的差分格式，如在节点 0 处有 其截断误差阶为 ? 2021/3/23 * 如上处理边值条件缺点是破坏了对称正定性,而这一性质是五点差分格式所固有的．为了保持对称正定性,可用下式代替 此时,其截断误差阶为 (3.7) ? 2021/3/23 * 3.3五点差分格式（积分插值法） A B C D 3.3.1 内点 2021/3/23 * A B C D 2021/3/23 * 2021/3/23 * 仅讨论一种特殊情况，假设 中的节点是两族网线的交点．p0(i,j)是边界网点，p1(i+1,j)、p2(i,j-1)是与之相邻的内点，在△ABC上对 (3. 8) 积分，利用Green公式得 3.3.2边界 讨论第二、三边值条件 (3.8) 而 T B A C 2021/3/23 * 于是 (3.9) 就是逼近 (3.8) 的差分方程. (3.9) 于是有 2021/3/23 * 考虑 Poisson 方程 作业 边界点用增设虚拟网点来近似法向导数,试构造5点差分格式 2021/3/23 * 实例 边界点用增设虚拟网点来近似法向导数,试构造5点差分格式,并用对称矩阵与向量表示。 2021/3/23 * 实例 试构造5点差分格式,并用对称矩阵与向量表示。 2021/3/23 * 3.4 三角网差分格式 (3.1) G 的边界 T 为分段光滑曲线. (第一边值条件) (3.2) (第二边值条件) (3.3) (第三边值条件) (3.4) 仍然考虑 Poisson 方程 在 T 上 u 满足下列边值条件之一: 2021/3/23 * 三角剖分: 1）在边界上取一系列的点,以其为顶点连成闭折线 2）将闭折线所围区域分割成有限个三角形之和,满足: ①任意一个三角形的顶点与其它三角形或者不相交,或者仅仅与其他三角形的顶点相交; ②三角形的每个内角不大于。 2021/3/23 * 引入如下术语: 节点—三角形的顶点; 单元—每个三角形; 相邻节点—同一条边上的两个节点互为相邻节点; 相邻单元—有一公共边的两个三角形互为相邻单元; 对偶单元—对于任一节点,考虑所有以它为顶点的三角形单元和以它为顶点的三角形边,过每一条边作中垂线,交于外心,得到围绕该节点的小多边形,称为对偶单元。 对偶剖分—全体对偶单元构成区域G的一个新的网格剖分,称为对偶剖分。 2021/3/23 * 2021/3/23 * 2021/3/23 * §4 极值定理 敛速估计 Poisson 方程 一般二阶线性椭圆型微分方程 2021/3/23 * 记U=(u1,1,u2,1,...uN-1,1,u1,2,u2,2,...uN-1,2,...u1,N-1,u2,N-1,...uN-1,N-1) 上述方程组可以表示为 AU=F,A具有下列性质 1）每行最多5个数非零,稀疏矩阵 2）对角元是正的,非对角元是非正的,且满足 记 2021/3/23 * 定理4.1 设 （负的极小），除非 是 上任一网函数，若 对任意 则 不可能在内点取正的极大 证明 用反证法， 设 中某点达到正 的极大值 M . 由于 联通，必有某一内点 使 且至少有一个相邻网点，比如 使 于是有 矛盾． 注　本定理说明在定理中所给条件下,若有正的极大（负的极小）,则一定在边界上取得． 2021/3/23 * 推论4.1 差分方程有唯一解． 证明:差分方程组可以表示为 AU=F,要证明差分方程有唯一解,只需证明AU=F有唯一解,即证明相应的齐次方程只有平凡解（唯一零解）,与齐次方程对应的差分方程满足:边值和右端都恒为零（齐次问题）。 uij不可能在内点取正的极大或者负的极小,也就是说,正的极大或者负的极小只能在边界处取得,但边界处uij=0,所以uij≡0 2021/3/23 * 推论4.2 若网函数 满足 则 证明： 所以uij不可能在内点取负的极小，也就是说，负的极小只能在边界处取得，所以uij≥0 2021/3/23 * 证明　由条件可知 定理4.2 设 是两个网函数，满足 则 由推论4.2 知,结论成立． 返回 2021/3/23 * 由定理4.2知： 推论4.3 的解 满足 证明　设 是如下问题的解 若 则 若 由定理4.1知， 的最大值只能在 达到，因此 从而结论成立． 2021/3/23 * 4.2 五点