文章目录
前言
本文学习过程来源是《矩阵分析与应用-张贤达》一书. 可以通过 z-lib 下载.
一、 Q R \mathrm{QR} QR 分解的性质
定理 1 ( Q R \mathrm{QR} QR 分解 ): 若 A ∈ R m × n A \in R^{m \times n} A∈Rm×n, 且 m ≥ n m \ge n m≥n, 则存在列正交的矩阵 Q ∈ R m × m Q \in R^{m \times m} Q∈Rm×m 和上三角矩阵 R ∈ R m × n R \in R^{m \times n} R∈Rm×n 使得 A = Q R A = QR A=QR
当 m = n m = n m=n 时, Q Q Q 是正交矩阵. 如果 A A A 是非奇异的 n × n n \times n n×n 矩阵, 则 R R R 的所有对角线元素均为正, 并且在这种情况下 Q Q Q 和 R R R 二者是唯一的. 若 A A A 是复矩阵, 则 Q Q Q 和 R R R 取复值.
( 补充一个小知识: Q Q Q 其实是源于 orthogonal matrices 或 orthonormal basis, 为了避免 O 0 O \ 0 O 0 不好区分的问题; R R R 是指代 right triangular matrices )
引理 1 : 若 A A A 和 B B B 是任意两个 m × n m \times n m×n 矩阵, 则
A H A = B H B (1) A^{\mathrm{H}}A = B^{\mathrm{H}}B \tag{1} AHA=BHB(1)
当且仅当存在一个 m × m m \times m m×m 酉矩阵 Q Q Q, 使得
Q A = B (2) QA = B \tag{2} QA=B(2)
二、采用修正 G r a m − S c h m i d t \mathrm{Gram-Schmidt} Gram−Schmidt 法的 Q R \mathrm{QR} QR 分解
矩阵 A A A 的 Q R QR QR 分解可以利用 G r a m − S c h m i d t \mathrm{Gram-Schmidt} Gram−Schmidt 正交化方法实现.
这个方法在之前就已经提到过, 这里就不再赘述了. 我们最后要得到的就是一组标准正交基 q 1 , q 2 , ⋯ , q n q_1,q_2,\cdots,q_n q1,q2,⋯,qn
容易验证, q i q_i qi 是标准正交基, 即满足
q i H q j = δ i j (3) q_i^{\mathrm{H}}q_j = \delta_{ij} \tag{3} qiHqj=δij(3)
其中, δ i j \delta_{ij} δij 为 K r o n e c k e r δ \mathrm{Kronecker} \ \delta Kronecker δ 函数. 简单说一下这个函数, 它和 D i r a c δ \mathrm{Dirac} \ \delta Dirac δ 非常相似.
δ i j = { 0 i f i ≠ j 1 i f i = j (4) \delta_{ij} = \left\{\begin{matrix} 0& \mathrm{if} \ i \neq j\\ 1& \mathrm{if} \ i = j \end{matrix}\right. \tag{4} δij={01if i=jif i=j(4)
如果令 m × n m \times n m×n 矩阵 A A A 的列向量为 a 1 , a 2 , ⋯ , a n a_1,a_2,\cdots,a_n a1,a2,⋯,an, 则以 q 1 , q 2 , ⋯ , q n q_1,q_2,\cdots,q_n q1,q2,⋯,qn 为列向量的矩阵 Q Q Q 与 A A A 之间有下列关系:
A = Q R (5) A=QR \tag{5} A=QR(5)
又由于 q i q_i qi 组成标准正交基, 所以
Q H Q = I n Q^{\mathrm{H}}Q=I_n QHQ=In
这里提出一种修正正交化算法, 之前是使上三角矩阵 R R R 的元素不是按列, 而是按行计算.
和经典正交化算法不同的是, 经典算法是一个一个基计算, 计算到某个基时就减去与之前求得的基相平行的分量. 修正后的算法就是每算出一个基就对后面未正交化的基执行减去平行分量的操作.
修正后的正交化算法减小误差的效果更明显.
三、 H o u s e h o l d e r \mathrm{Householder} Householder Q R \mathrm{QR} QR 分解
H o u s e h o l d e r \mathrm{Householder} Householder 变换可以实现任意 m × n m \times n m×n 矩阵 A A A 的 Q R QR QR 分解, 其原理时使用变维向量的 H o u s e h o l d e r \mathrm{Householder} Householder 变换, 使得该向量除第一个元素外, 其他元素皆变成 0.
使一个 p p p 维向量 x = [ x 1 , x 2 , ⋯ , x p ] T x=[x_1,x_2,\cdots,x_p]^{\mathrm{T}} x=[x1,x2,⋯,xp]T 的第 1 个元素后面的所有元素变为 0, 则 p p p 维的 H o u s e h o l d e r \mathrm{Householder} Householder 向量应取
ω = x − β e 1 β ˉ ( β − x 1 ) (6) \omega = \frac{x - \beta e_1}{\sqrt{\bar{\beta} (\beta - x_1)}} \tag{6} ω=βˉ(β−x1)x−βe1(6)
式中
β ˉ = − ∣ x 1 ∣ ∥ x ∥ , β = − x 1 ∣ x 1 ∣ ∥ x ∥ (7) \bar{\beta} = - |x_1| \lVert x \rVert, \quad \beta = - \frac{x_1}{|x_1|} \lVert x \rVert \tag{7} βˉ=−∣x1∣∥x∥,β=−∣x1∣x1∥x∥(7)
首先令 x = a 1 = [ a 11 , a 21 , ⋯ , a m 1 ] T x=a_1=[a_{11},a_{21},\cdots,a_{m1}]^{\mathrm{T}} x=a1=[a11,a21,⋯,am1]T, 并取 p = m p=m p=m, 按照式子 (6) 和式子 (7), 可以计算得到 u 1 = ω m u_1=\omega_m u1=ωm. 此时,
H 1 = I − u 1 u 1 T → A 1 = H 1 A = [ a 1 ( 1 ) , a 2 ( 1 ) , ⋯ , a n ( 1 ) ] (8) H_1 = I - u_1 u_1^{\mathrm{T}} \to A_1 = H_1A = [a_1^{(1)},a_2^{(1)},\cdots,a_n^{(1)}] \tag{8} H1=I−u1u1T→A1=H1A=[a1(1),a2(1),⋯,an(1)](8)
变换后矩阵 A 1 A_1 A1 的第 1 列 a 1 ( 1 ) a_1^{(1)} a1(1) 的第一个元素等于 ( a 11 2 + a 21 2 + ⋯ + a m 1 2 ) 1 / 2 (a_{11}^2 + a_{21}^2 + \cdots + a_{m1}^2)^{1/2} (a112+a212+⋯+am12)1/2 (像第二范数). 其他元素全为 0.
接下来就是针对矩阵 A 1 A_1 A1 的第 2 列 a 2 ( 1 ) a_2^{(1)} a2(1), 令 p = m − 1 p = m-1 p=m−1 和
x = [ a 22 ( 1 ) , a 32 ( 1 ) , ⋯ , a m 2 ( 1 ) ] T x = [a_{22}^{(1)},a_{32}^{(1)},\cdots,a_{m2}^{(1)}]^{\mathrm{T}} x=[a22(1),a32(1),⋯,am2(1)]T
按照式子 (6) 和式子 (7), 可以计算得到 ( m − 1 ) (m-1) (m−1) 维向量 ω m − 1 \omega_{m-1} ωm−1. 此时, 取 u 2 = [ 0 ω m − 1 ] u_2 = \begin{bmatrix} 0\\ \omega_{m-1}\end{bmatrix} u2=[0ωm−1], 又可得到
H 2 = I − u 2 u 2 T → A 2 = H 2 A 1 = H 2 H 1 A = [ a 1 ( 1 ) , a 2 ( 2 ) , ⋯ , a n ( 2 ) ] (9) H_2 = I - u_2 u_2^{\mathrm{T}} \to A_2 = H_2A_1 = H_2H_1A = [a_1^{(1)},a_2^{(2)},\cdots,a_n^{(2)}] \tag{9} H2=I−u2u2T→A2=H2A1=H2H1A=[a1(1),a2(2),⋯,an(2)](9)
可以发现矩阵 A 2 A_2 A2 的第 1 列与 A 1 A_1 A1 的第 1 列相同, 而第 2 列 a 1 ( 1 ) a_1^{(1)} a1(1) 的第一个元素等于 a 12 ( 1 ) a_{12}^{(1)} a12(1), 第二个元素等于 [ ∣ a 22 ( 1 ) ∣ 2 + ∣ a 32 ( 1 ) ∣ 2 + ⋯ + ∣ a m 2 ( 1 ) ∣ 2 ] 1 / 2 [ |a_{22}^{(1)}|^2 + |a_{32}^{(1)}|^2 + \cdots + |a_{m2}^{(1)}|^2 ]^{1/2} [∣a22(1)∣2+∣a32(1)∣2+⋯+∣am2(1)∣2]1/2 , 而该列的其他元素全部为 0, 以此类推.
假定矩阵 A A A 经过 k − 1 k -1 k−1 次 H o u s e h o l d e r \mathrm{Householder} Householder 变换后, 已变成 A ( k − 1 ) A^{(k-1)} A(k−1) 即有
A ( k − 1 ) = H k − 1 A ( k − 2 ) = H k − 1 ⋯ H 1 A = [ a 1 ( k − 1 ) , a 2 ( k − 1 ) , ⋯ , a n ( k − 1 ) ] , k = 2 , 3 , ⋯ \begin{aligned} A^{(k-1)} &= H_{k-1}A^{(k-2)}=H_{k-1} \cdots H_1A &= [a_1^{(k-1)},a_2^{(k-1)},\cdots,a_n^{(k-1)}], \quad k = 2,3,\cdots \end{aligned} A(k−1)=Hk−1A(k−2)=Hk−1⋯H1A=[a1(k−1),a2(k−1),⋯,an(k−1)],k=2,3,⋯
并且其前 k − 1 k-1 k−1 列具有以下变换结果:
a j ( k − 1 ) = [ a 1 j ( k − 1 ) , ⋯ , a j j ( k − 1 ) , 0 , ⋯ , 0 ] T , j = 1 , 2 , ⋯ , k − 1 a_j^{(k-1)} = [a_{1j}^{(k-1)},\cdots,a_{jj}^{(k-1),0,\cdots,0}]^{\mathrm{T}}, \quad j=1,2,\cdots,k-1 aj(k−1)=[a1j(k−1),⋯,ajj(k−1),0,⋯,0]T,j=1,2,⋯,k−1
因此, 第 k k k 次 H o u s e h o l d e r \mathrm{Householder} Householder 变换的目的就是保证前 k − 1 k-1 k−1 列不变, 实现 A ( k − 1 ) A^{(k-1)} A(k−1) 的第 k k k 列的下述变换:
H k ~ [ a k , k ( k − 1 ) a k + 1 , k ( k − 1 ) ⋮ a m , k ( k − 1 ) ] = [ a k , k ( k ) 0 ⋮ 0 ] \tilde{H_k}\begin{bmatrix} a_{k,k}^{(k-1)}\\ a_{k+1,k}^{(k-1)}\\ \vdots \\ a_{m,k}^{(k-1)} \end{bmatrix} = \begin{bmatrix} a_{k,k}^{(k)}\\ 0\\ \vdots \\ 0 \end{bmatrix} Hk~⎣⎢⎢⎢⎢⎡ak,k(k−1)ak+1,k(k−1)⋮am,k(k−1)⎦⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎡ak,k(k)0⋮0⎦⎥⎥⎥⎤
这相当于对矩阵 A ( k − 1 ) A^{(k-1)} A(k−1) 进行 H o u s e h o l d e r \mathrm{Householder} Householder 变换 H k A ( k − 1 ) H_kA^{(k-1)} HkA(k−1) 时取
H k = [ I k − 1 0 0 H k ~ ] H_k = \begin{bmatrix} I_{k-1} & 0\\ 0 & \tilde{H_k} \end{bmatrix} Hk=[Ik−100Hk~]
取 n n n 次 H o u s e h o l d e r \mathrm{Householder} Householder 变换后, 即可实现 Q R \mathrm{QR} QR 分解.