高等数值算法与应用-wenjianyu.pdfVIP

高等数值算法与应高等数值算法与应用 （八） AdvancedAdvanced NumericalNumerical Alggorithms Apppplications 计算机科学与技术系计算机科学与技术系 喻文健喻文健 内容概要内容概要 线性方程组的迭代解法（三） 对称矩阵的对称矩阵的AArnoldildi过程过程—LLanczos过程过程 简 从从Lanczos过程到共轭梯度法过程到共轭梯度法 单单 第三类Krylov子空间迭代法基础— 介 LanczosLanczos双正交化过程双正交化过程 绍绍 构造预条件矩阵的技术 Wenjian Yu 2 Krylov Methods - Symmetric Lanczos PProceddure andd CGCG Wenjian Yu 3 Arnoldi过程过程: 为为Krylov子空间子空间K (A, v )构造构造一组组 m 1 单位正交基向量 Arnoldi过程的矩阵表示： (H 为上Hessen- m bberg矩阵矩阵)) 00 AV V H m m+1 m TT T VV AVAV HH V H =+w e m m m m m m m 考虑对实对称矩阵A应用Arnoldi过程，以及L=K 的的正交投影方法交投影方法，其系数矩阵为其系数矩阵为H (如果A对称正定， mm 此方法不会恶性中断) T T T H V A V H 对称矩阵！ m m m m + H 为对称三对角矩阵 上Hessenberg矩阵 m Wenjian Yu 4 Arnoldi过程过程(针对实对称阵) 1)1) 选选vv 使使 vv1 11 1 2 2) For j=1, …, m z := Av j For i=1, …, jj -1, jj