可证,在一定条件下,基本QR方法产生的矩阵序 列{A()}“基本”收敛于一个上三角阵(或分块上 三角阵)。即主对角线(或主对角线子块)及其以下 元素均收敛,主对角线(或主对角线子块)以上元素 可以不收敛。特别的,如果A是实对称阵,则{A()} “基本”收敛于对角矩阵。 因为上三角阵的主对角元(或分块上三角阵中, 主对角线子块的特征值)即为该矩阵的特征值,故当k 充分大时,A()的主对角元(或主对角线子块的特征 值)就可以作为A的特征值的近似 (4)-2)=0 基本的QR方法的主要运算是对矩阵QR分解,分 解的方法有多种。介绍一种Schm正交化方法
可证，在一定条件下，基本QR方法产生的矩阵序 列｛A(k)｝ “基本”收敛于一个上三角阵（或分块上 三角阵）。即主对角线（或主对角线子块）及其以下 元素均收敛，主对角线（或主对角线子块）以上元素 可以不收敛。特别的，如果A是实对称阵，则｛A(k)｝ “基本”收敛于对角矩阵。 因为上三角阵的主对角元（或分块上三角阵中， 主对角线子块的特征值）即为该矩阵的特征值，故当k 充分大时， A(k)的主对角元（或主对角线子块的特征 值）就可以作为A的特征值的近似。 基本的QR方法的主要运算是对矩阵QR分解，分 解的方法有多种。介绍一种Schmit正交化方法。 ( ) ( ) 0 k A u − =
设A为mˉ阶非奇异实矩阵,记A=[a12a2…,an,其中 )(j=1,2,…,n) Exb,a,/,. b2=a-a, b 6=an al (a2a4)(a2a)(a12a1 0 a.a ila h⊥b2取b2=b2则=临=1.(,b2)=0 般地取b=4-∑(a2b为,b=b/(k=2,3…n) 则向量组bb2…,b正交,且 1(k=1,2,…n) 即ak=(a,b)b+…+(a4b1-)b-+|1
1 2 1 2 ' 2 1 1 1 1 2 2 2 1 1 2 1 2 1 2 1 2 1 2 1 1 1 ' 2 1 1 2 1 1 1 1 1 1 ' ' ' 1 2 2 2 2 1 2 [ , , , ], ( , , , ) ( 1, 2, , ). , / . , , , , , , 0 , , / , 1 n T j j j nj A n A a a a a a a a j n a a b a a b a a b b a a a a a a a a a a a a a b b a a a a a a b b b b b b b = = = = = − = − = − = − = ⊥ = = = 设 为 阶非奇异实矩阵，记 其中 取 取 则 1 2 1 ' ' ' 1 1 2 ' 1 1 1 1 , , 0. , , / ( 2,3, ) , , , 1 ( 1, 2, ) , , k k k k i i k k k i n k k k k k k k k b b b a a b b b b b k n b b b b k n a a b b a b b b b − = − − = = − = = = = = + + + 一般地取 则向量组 正交，且 即
即ak=(akb)b1+…+(ab-)b2-1+|bb 于是A=[a12a2…,an]=[b1,b2…,bn a1‖(a2b1 Kan,b,) b2‖ b OR b‖ 这就是用 Schmit正交化方法对矩阵进行的QR分解。 基本QR方法每次迭代都需作一次QR分解与矩阵乘法,计算 量大,而且收敛速度慢。因此实际使用的QR方法是先用一系列 相似变换将A化成拟上三角矩阵(称为上 Hessenberg矩阵),然 后对此矩阵用基本QR方法。因为拟上三角矩阵具有较多零元素, 故可减少运算量。化A为相似的拟上三角阵的方法有多种
' 1 1 1 1 1 2 1 2 1 2 1 1 ' 2 2 ' 1 1 ' , , [ , , , ] [ , , , ] , , , , QR k k k k k k k n n n n n n n n a a b b a b b b b A a a a b b b a a b a b b a b QR b a b b Schmit − − − − = + + + = = = 即 于是 这就是用 正交化方法对矩阵进行的 分解。 基本QR方法每次迭代都需作一次QR分解与矩阵乘法，计算 量大，而且收敛速度慢。因此实际使用的QR方法是先用一系列 相似变换将A化成拟上三角矩阵（称为上Hessenberg矩阵），然 后对此矩阵用基本QR方法。因为拟上三角矩阵具有较多零元素， 故可减少运算量。化A为相似的拟上三角阵的方法有多种
2-10 例用 Schmit正交化方法对A=-12 进行QR分解 0 解:a1=(2,-1,0),a2=(-12,-1),a3=(0,-1,2) 因而有b √5√5 0) b2=a2-(a2,b)b=(-12,-1)+层( /3·√ 0) 1) 55 36 70√70 70 246 (a32b〉b-(a2b2)b2 123 √14~√4√14
1 2 3 1 1 1 2 ' 2 2 2 1 1 ' 2 2 ' 2 2 2 1 0 : 1 2 1 . 0 1 2 : (2, 1,0) , ( 1,2, 1) , (0, 1,2) . 2 1 ( , ,0) , 5 5 4 2 1 3 6 , ( 1,2, 1) ( , ,0) ( , , 1) , 5 5 5 5 5 3 6 5 ( , , ) , 70 70 70 T T T T T T T T Schmit A QR a a a a b a b a a b b b b b b − = − − − = − = − − = − = = − = − = − − + − = − = = − 例用 正交化方法对 进行 分解 解 因而有 ' 3 3 3 1 1 3 2 2 ' 3 3 ' 3 2 2 4 6 , , ( , , ) 777 1 2 3 ( , , ) 14 14 14 T T a a b b a b b b b b = − − = = =
所以A=[a12a2…an]=[b2b2…,bn] (a,b) (a,b 2) 1bm-l( a, bm-1, 3/701454/515 J√56/√72√4‖070/5-16/70 s/703√400240/7
1 2 1 2 1 2 1 1 ' 2 2 ' 1 1 ' [ , , , ] [ , , , ] , , , , 2 5 3 70 1 14 5 4 5 1 5 1 5 6 70 2 14 0 70 5 16 70 0 5 70 3 14 0 0 2 40 7 n n n n n n n n A a a a b b b a a b a b b a b b a b b − − = = − = − − − 所以