研究性教学-有限差分法.pptVIP

数值计算 有限差分法 1. 差分格式 有限差分的网格划分 有限差分的网格划分 有限差分的网格划分 有限差分的网格划分 2. 差分方程组的解 * 静态场的边值问题：在给定的三类边界条件下求解标量位或矢量位的泊松方程或拉普拉斯方程的解的问题。 ? 边值问题的分类： 解析法 数值法 分离变量法 镜像法 复变函数法 保角变换法 格林函数法 有限差分法 有限元法 数值积分法 矩量法 （经典法） 优点：由有限个想构成闭合解或无穷级数，通常具有鲜明的物理意义。 缺点：解题范围窄小，对边界形状十分挑剔 优点：解题范围宽广，对边界的形状没有限制 缺点：数据离散，物理意义深藏其中 ? 理论依据：唯一性定理和叠加原理 在由边界L界定的二维区域D内，电位函数φ满足拉普拉斯方程且给定第一边界条件，则 将区域D划分为正方形网格，网格线的交点称为节点，两相邻平行网格线间的距离称为步距 h 拉普拉斯方程离散化，对于任一点0，有一阶偏导数 对于二阶偏导数 同样 拉普拉斯方程的差分格式为 紧邻边界节点的拉普拉斯方程的差分格式为 紧邻边界节点的网格划分 其中p、q为小于1的正数；1、2为边界上的节点，其值为对应边界点处的值，是已知的。 ? 高斯-赛德尔迭代法 网格节点一般按“自然顺序”排列，即从左到右，再从下到上顺序排列。迭代也按自然顺序进行。 网格节点排列方法 网格节点排列方法 计算步骤 首先对节点 取迭代初值 。再按下式反复迭代 其中 上角标（k）表示k次近似值，下脚标i，j表示节点所在位置，即第i行第j列的交点。 网格节点排列方法 计算步骤 注意：在迭代过程中遇到边界点式，需将边界条件  带入。 所有内节点满足以下条件时迭代停止 W是预定的最大允许误差 ? 逐次超松弛法 高斯-赛德尔迭代法德变形。为加速收敛，相应的迭代公式为 α称为“加速收敛因子”，且1《 α 2。 注意：逐次超松弛法收敛的快慢与α有明显关系。如何选择最佳α，是个复杂问题。 迭代求解程序框图 应用有限差分法求解静电场边值问题的近似解 解：取h=5 作正方形网格，得差分方程 例 利用高斯-赛德尔迭代格式 选取迭代初值 ，6次迭代计算后，可得结果 *