5有限差分法教程分析.pptVIP

cem@uestc.edu.cn 超松弛迭代法 求解具有稀疏系数矩阵的大型差分方程组，其中最优的就是超松弛迭代法（Successive Over Relaxation，SOR）。为了说明SOR方法，首先介绍雅可比法和高斯-赛德尔法 雅可比法（Jacobi）就是要使迭代值能精确的满足前一次各点的电位值所能表示的差分方程 cem@uestc.edu.cn 超松弛迭代法 高斯-赛德尔法是雅可比法的改进方法，主要针对减少内存消耗，只需存储一组完整的数组。它采取的措施是对每一次迭代尽量采用最新计算的值来替换上一次迭代的旧值。结果收敛速度比雅可比法快一倍。 第5章 有限差分法 工程电磁场数值方法编程实验－有限差分法 电子科技大学 赖生建 cem@uestc.edu.cn cem@uestc.edu.cn 主要内容 差分和差商 有限差分格式 不同媒质分界面上的差分格式及定解问题的差分格式 有限差分法的求解 场强与电、磁积分量的计算 典型算例分析 cem@uestc.edu.cn 介绍 有限差分方法是一种微分方法，自上世纪五十年代以来得到了广泛的应用，该方法概念清晰，方法简单，直观。虽然其与变分法相结合所形成的有限元法更有效，但有限差分还是以其固有特点在数值计算中有其重要地位，是应用最多的一种数值方法。 为求解由偏微分方程定解问题所构造的数学模型，有限差分法是将定解区域（场区）离散化为网格离散节点的集合。并以各离散点上函数的差商来近似该点的偏导数，使待求的偏微分方程定解问题转化为一组相应的差分方程。根据差分方程组解出各离散点处的待求函数值——离散解。 cem@uestc.edu.cn 1、差分与差商 用差分代替微分，是有限差分法的基本出发点。这一点由微分原理保证的，当自变量的差分趋于零时，差分变成微分 cem@uestc.edu.cn 差分与差商 前向差分 后向差分 中心差分 cem@uestc.edu.cn 差分与差商 通过泰勒公式分析上面差分精度，在点上的一阶导数的逼近度可由泰勒公式展开 两式相减 cem@uestc.edu.cn 差分与差商 前向、后向差分截断于 ，具有h的一阶精度，而中心差分法截断于 ，具有h的二阶精度，中心差分的精度比较高。 前向差分 前向差分 cem@uestc.edu.cn 差分与差商 对偏导数，可仿照上述方法，将表示为： cem@uestc.edu.cn 差分格式 二维Possion方程差分格式 有限差分法的网格划分，通常采用完全有规律的分布方式，这样可使每个离散点上得到相同形式的差分方程，有效的提高解题速度。对能填满平面域的三种规则网格（正方形，正三角形和正六边形）的划分方式，经常采用的是正方形网格划分， cem@uestc.edu.cn 差分格式 一阶偏导数差分格式 可采用待定系数的方法，提高差分格式的精度，它的思路： 1、3结点与0结点在x方向的差分用泰勒公式展开，它们各自占有一定的权系数，以截断误差来计算系数 cem@uestc.edu.cn 差分格式 忽略h3以上的高次幂的项，并且令 项的系数为零，这样处理可以保证得到的差分格式误差为h3量级。系数为零的条件 求出二阶精度精度为一阶偏导数差分格式 cem@uestc.edu.cn 差分格式 二阶偏导数的差分格式 令方程右边的一阶偏导数的系数为0，得到系数间的表达式 代入上式得到精度为O(h3)的二阶偏导数的差分格式 cem@uestc.edu.cn 差分格式 当 时，上式可以简化为 Possion方程五点差分格式 cem@uestc.edu.cn 不同媒质分界面上的差分格式 分界面与网格线重合的情况 两式中 和 是假设“虚”电位，可以利用分界面上场量遵循的边界条件，削去它们 cem@uestc.edu.cn 不同媒质分界面上的差分格式 其次，假设在分界面上没有自由电荷 中心差分格式表示 把前面关于 和 式子代入上式 cem@uestc.edu.cn 不同媒质分界面上的差分格式 分界面与网格线呈对角线的情况 两式中 和 是假设“虚”电位，可以利用分界面上场量遵循的边界条件，削去它们 cem@uestc.edu.cn 不同媒质分界面上的差分格式 其次，假设在分界面上没有自由电荷 对M、N结点应用线性插值 cem@uestc.edu.cn