计算物理学案3矩阵特征值运算概述.ppt

西安邮电学院应用物理系 (1) QR基本算法 QR基本算法的迭代过程如下: ① 给定循环步数M和迭代初始阵 A; ② 利用以下迭代公式计算: Ak=QkRk Ak+1=RkQk ③ 当A的特征值满足|? 1||? 2||? 3|≥|? 4|≥…≥|? n|时,经 过 M 次迭代,该算法产生的矩阵序列基本收敛到上三角阵, 对角元素收敛到 A 的特征值。 QR基本算法的求解程序如下: function e=qrtz(A, M) for k=1:M [Q R]=qr(A); A=R*Q; end e=diag(A); 西安邮电学院应用物理系 (2) 海森伯格阵的QR算法 海森伯格阵QR算法的迭代过程如下: ① 将矩阵A经相似变换变为上海森伯格阵; ② 实施QR基本算法, 即 给定循环步数M和迭代初始阵 A1=hess(A); 利用以下迭代公式计算: Ak=QkRk Ak+1=RkQk 海森伯格阵QR算法求解程序如下: function e=hessqrtz(A, M) A=hess(A); for k