数值分析23(免费阅读).ppt

第二章 线性方程组求解 差商和差分的基本概念 拉普拉斯方程和泊松方程的差分形式 求解静态场第一类边值问题 差商和差分的基本概念 有限差分法: 数值分析中最早应用的一种计算方法 简单、直观 应用广泛：常微分方程、偏微分方程、初值问题、边值问题、二阶线性方程和非线性方程等等，转化为代数方程组后用计算机求解 差商和差分的基本概念 设函数 f(x)，其独立变量 x 有一很小的增量△x，则相应函数的增量为： △f (x)是函数f (x) 在 x 处沿 x正方向的改变量， 称为函数 f(x)的一阶向前差分 又因为△f (x)为有限量，故又称为有限差分 差商和差分的基本概念 取向前、向后差分的平均值（中心差分）： 同样，一阶向后差分： 如果 很小，差分和微分之间的差异很小，这是因为 差商和差分的基本概念 用x的增量△x去除差分，便得差商（一阶）： 向前差商： 向后差商： 中心差商： 它们对一阶导数的逼近度可通过泰勒公式的展开式得到 差商和差分的基本概念 泰勒公式： 则 向前差商： 一阶精度 向后差商： 一阶精度 上两式相减： 中心差商： 二阶精度 差商和差分的基本概念 二阶差商：对一阶差商再求差商 拉普拉斯方程和泊松方程的差分形式 二维电磁场泊松方程的边值问题 场域D为正方形： ，等步长 x y 0 二维静态场泊松方程： 拉普拉斯方程和泊松方程的差分形式 推导差分方程的方法： 采用泰勒级数法，将场量展开成泰勒级数，用差商代替偏导，得到差分方程。 不对称星形节点 拉普拉斯方程和泊松方程的差分形式 将任一点的位函数φ(x,y)沿x轴方向展开为O点位函数φ0的泰勒级数： 将节点A和C的坐标代入上式： 拉普拉斯方程和泊松方程的差分形式 为消去式中的一阶偏导 当h很小时，忽略h的三阶以上高次项，得到 同样可得 拉普拉斯方程和泊松方程的差分形式 将上两式代入O点的泊松方程 得到 上述过程，直接用差商逐项逼近微分方程中的微商来推导差分方程，称为逐项逼近法。 拉普拉斯方程和泊松方程的差分形式 另外，由二阶差商的定义 代替 。同样，可得 。 将它们代入微分方程，便可同样得到差分方程。 拉普拉斯方程和泊松方程的差分形式 正方形网格（对称星形节点），即 代入到 无源区域，可得到拉普拉斯方程 或 静态场第一类边值问题求解 第一类边值问题： φ|C = g 0 1 2 3 4 h C 静态场第一类边值问题求解 y x 0 V= V0 V=0 V=0 V=0 二维电磁场泊松方程 求解方槽中的电位分布 只考虑网格和边界重合的情况 静态场第一类边值问题求解 所满足的差分方程 (1) 雅可比迭代: 缺点： 需要两套存储单元，占用内存较大； 收敛速度较慢 静态场第一类边值问题求解 (2) 高斯－赛德尔迭代: 特点： 占用内存较小； 收敛速度较快； 但是当网格数很大时，收敛速度仍然较慢 静态场第一类边值问题求解 (3) 超松弛迭代（SOR迭代） 特点： 适当的松弛因子 将大大加快收敛速度 高斯－赛德尔迭代的值作为一个中间结果 对 和 加权平均，得 静态场第一类边值问题求解 选用何种坐标系； 选用何种网格（包括边界）； 选用何种差分格式（特别是边界）； 问题的精确度要求； 初始值选择； 建立差分方程组； 选用何种迭代法，及参数选择； 检验迭代解的收敛和控制迭代结束； 对问题所得的数值解进行检验。 步骤： 静态场第一类边值问题求解 流程图： 静态场第一类边值问题求解 对正方形场域的第一类边界条件： l、m分别是每边的网格数 对矩形场域的第一类边界条件： l是每边的网格数 超松弛迭代（SOR迭代） 方槽中的电位分布 迭代次数： G-S迭代：8424 SOR迭代：296 静态场第一类边值问题求解 课程设计：“计算金属槽中的电位分布” 分别取网格步长为：0.2a、0.02a、0.002a，用雅克比迭代法、G-S迭代、SOR迭代和共轭梯度法（博士生要求必须做共轭梯度法）进行计算，并比较它们收敛的迭代次数，总结收敛次数与网格数之间的关系。精度要求为10-6. 参考《电子科技大学学报》投稿格式 1. 题目和摘要（中英文） 2. 正文： a 引言; b 原理； c 数值算例； d 结论。 3. 附件：程序 国庆假期结束后上课时交打印版和电子版（发邮件：xsyang@126c.om。 每章约一个设计题目，每人总共做2个设计题目。 格式要求