大规模矩阵问题的krylov子空间方法综述-科学工程培训讲义.pptVIP

大规模矩阵问题的Krylov 子空间方法综述 Krylov子空间方法综述 背景介绍 投影方法 Krylov子空间及其标准正交基 Krylov子空间方法求解线性方程组 Krylov子空间方法求解矩阵的特征值 研究热点和尚未解决的问题 背景 大规模线性方程组的求解 很多科学工程计算问题都转化为求解方程组Ax=b.如偏微分方程组的差分格式，有限元方法离散得到刚度矩阵. 大规模矩阵特征值和特征向量的计算 工程计算领域十分常见。如量子物理中的Kohn-Sham方程求解化为哈密顿矩阵某些关键特征值对的计算. 投影方法 线性方程组的投影方法 方程组Ax=b,A是n×n的矩阵. 给定初始x(0),在m维空间K(右子空间）中寻找x的近似解x(1) 满足残向量r=b-Ax(1) 与m维空间L(左子空间）正交，即 b-Ax(1) ⊥L 此条件称为Petrov-Galerkin条件 . 当空间K=L时，称相应的投影法为正交投影法，否则称为斜交投影法. 解方程组的投影法的矩阵表示 　　设n×m阶矩阵V=[v(1), v(2), …v(m)]与W=[w(1), w(2), …w(m)]的列分别构成K与L的一组基 .记z=x(1)-x(0),z=Vy,有 　 当　　　非奇异时（通常情况成立，见定理1.1），有 从而得到迭代公式 投影方法的最优性 １. (误差投影)