有限元法基础-7线性代数方程组的解法论述.ppt

7.4 迭代法 有限元法基础 共轭梯度法迭代公式 7.4 迭代法 有限元法基础 共轭梯度法特点 1）在迭代过程中，rk 是相互正交的 2）对n阶线性代数方程组通过n次迭代得到的rn一定是零 3）xn必定是原方程的解，由于有数值截断误差，计算时不一定得到精确解 另一加速算法为PCG（Preconditioned conjugate gradient method）收敛更快。 7.4 迭代法 有限元法基础 3.预条件共轭梯度法 (Preconditioned conjugate gradient method) 该法简称 PCG法，是一种加速共轭梯度法。PCG法通过引入预条件矩阵M，使方程系数矩阵的条件数降低，以达到提高收敛速度的目的。 7.4 迭代法 有限元法基础 线性方程组 引入对称正定矩阵 原方程转换为 改写为 其中 M称为预条件矩阵。当M为A的近似时， 接近单位矩阵，它的条件数近似为1，然后用CG法求解新方程组。 7.4 迭代法 有限元法基础 PCG法迭代公式 7.4 迭代法 有限元法基础 PCG法对M的要求 1）应是原系数矩阵A的近似 2）应容易求逆，并尽可能少占存储 例如 如何选择M成为PCG法的关键问题。 7.2 带状系数矩阵的直接法 有限元法基础 二.二维等带宽的高斯消去法 工作三角形 由于系数矩阵呈带状 每次消元只涉及包括 主元在内的一个三角 形内的元素，称为工 作三角形 。 7.2 带状系数矩阵的直接法 有限元法基础 二维等带宽高斯消去法公式 7.2 带状系数矩阵的直接法 二维等带宽存储(n×ND) 采用按行分解 I=i , J=j-i+1 新的循环界： r=max(j-ND+1, i-1) 有限元法基础 二维等带宽三角分解 7.2 带状系数矩阵的直接法 有限元法基础 三.一维变列高存储高斯消去法 采用按列分解 7.3 利用计算机外存的直接法 有限元法基础 主要解决计算机内存不足的问题，充分利用外存保存分解后的系数矩阵与未分解的系数矩阵，以达到小内存算大问题的目的。 7.3 利用计算机外存的直接法 有限元法基础 一.高斯消去法的特点 1）第m次消元过程中，所涉及的元素仅在三角形工作区内 消元行元素 7.3 利用计算机外存的直接法 有限元法基础 2）在m次消元过程中，前面的元素不再参加消元，后面的元素尚未参加消元。 3）在整个消元过程中，工作区自上向下运动。 为分块解法奠定基础 7.3 利用计算机外存的直接法 有限元法基础 二.分块解法 设允许使用内存为NA 在每一分块，NQ－ND行集成完毕，可进行消元修正，最后的ND行进入到下一块系数矩阵一起集成，消元修正。 7.3 利用计算机外存的直接法 有限元法基础 分块解法的特点 在每一分块中，系数矩阵的元素是先集成后消元修正 从求解的全过程看，系数矩阵的集成和消元修正是交替进行 7.3 利用计算机外存的直接法 有限元法基础 分块解法简单框图 7.3 利用计算机外存的直接法 有限元法基础 三.波前法（Front Met