计算物理学案3 矩阵特征值运算.pptVIP

西安邮电学院应用物理系 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 1:M [Q R] qr A ; A R*Q; end e diag A ; 西安邮电学院应用物理系 例5: QR基本算法和海森伯格阵QR算法求矩阵A的特征值 A [1 5 6; 4 7 0; 8 11 4]; e1 qrtz A,20 e2 qrtz A,23 e3 hessqrtz A,20 输出计算结果分别为: e1 13.4695 -3.8565 2.3869 e2 13.4695 -3.8566 2.3871 e3 13.4695 -3.8566 2.3871 e4 13.4695 -3.8566 2.3871 西安邮电学院应用物理系 3 位移QR算法 位移QR算法的迭代过程如下: ① 给定循环步数M和迭代初始阵 A1 hess A ; ② 利用以下迭代公式计算: Ak -? kI QkRk Ak+1 RkQk + ? kI ③ ? k 的选择一般有以下两种: 瑞利商位移 威尔金森位移: 如果子矩阵 的两个特征值为实数, 则选离 最近的那个作为 。 西安邮电学院应用物理系 瑞利商位移QR算法求解程序如下: function e rqrtz A, M A hess A ; for k 1:M n size A,1 ; u A n,n ; [Q R] qr A-u*eye n ; A R*Q+u*eye n ; end e diag A ; [Q, R] qr A-uI A RQ+uI 开始 k 1, 2, 3, …, M A hess A n size A,1 u Ann e diag A 若 k M 若 k M 西安邮电学院应用物理系 威尔金森位移QR算法求解程序如下: function e wilkqrtz A, M A hess A ; for k 1:M n size A,1 ; A1 A n-1:n,n-1:n ; t eig A1 ; t diag t ; if imag t 1 0 imag t 2 0 if abs t 1 -A n,n abs t 2 -A n,n u t 1 ; else u t 2 ; end else u A n,n ; end [Q R] qr A-u*eye n ; A R*Q+u*eye n ; end e diag A ; [Q, R] qr A-uI A RQ+uI 开始 k 1, 2, 3, …, M A1 A n-1:n,n-1:n 求出A1本征值 t u t1 t ? R A hess A ? 1 |t1-ann| ? 2 |t2-ann| n size A,1 u t2 t ? C u Ann ? 1 ? 2 ? 1 ? 2 e diag A 若 k M 若 k M 西安邮电学院应用物理系 例6: 位移QR算法求矩阵A的特征值 A [1 5 6; 4 7 0; 8 11 4]; e1 rqrtz A,13 e2 wilkqrtz A,20 输出计算结果分别为: e1 13.4695 -3.8566 2.3871 e2 13.4695 -3.8566 2.3871 e3 13.4695 -3.8566 2.3871 西安邮电学院应用物理系 7. Schur分解和奇异值分解 Schur分解可以直接给出矩阵的特征值; 而对于对称矩 阵, 奇异值分解给出矩阵特征值的绝对值。 求解方法如下: [S U