11QR分解

1.定义

矩阵A可以化为正交（酉）矩阵Q和上三角矩阵R的乘积，即$A=QR$，则称上式为矩阵的$QR$分解

2.定理

若$A$是一个$n$阶非奇异矩阵，则存在正交（酉）矩阵$Q$和上三角矩阵$R$使得$A=QR$，且除去相差一个对角元素绝对值全为1的对角因子外，上述分解唯一。
回顾一下施密特正交化的过程

可以得到施密特正交化的矩阵为

证明：
设$A=[a_1,a_2...a_n]$,且$a_1,a_2,...a_n$线性无关。应用施密特正交化的办法可以得到$$\left[\begin{array}{cccc}{a}_1& a_2& \cdots & a_n\end{array}\right]=\left[\begin{array}{cccc}{b}_1& b_2& \cdots & b_n\end{array}\right]\left[\begin{array}{cccc}1& b& \cdots & k_{n1}\\ & 1& \cdots & k_{n2}\\ & & \ddots & ⋮\\ & & & k_{nn-1}\end{array}\right]$$对其单位化可得$$\left[\begin{array}{cccc}{a}_1& a_2& \cdots & a_n\end{array}\right]=\left[\begin{array}{cccc}{q}_1& q_2& \cdots & q_n\end{array}\right]\left[\begin{array}{cccc}|b_1|& & & \\ & |b_2|& & \\ & & \ddots & \\ & & & |b_n|\end{array}\right]\left[\begin{array}{cccc}1& b& \cdots & k_{n1}\\ & 1& \cdots & k_{n2}\\ & & \ddots & ⋮\\ & & & k_{nn-1}\end{array}\right]=QR$$接下来我们证明其唯一性：
假设有两个QR分解$$A=QR=Q_1{R}_1$$ $$Q=Q_1{R}_1{R}^{-1}=Q_{1}D$$因为$R_1,R$都为上三角矩阵，而上三角矩阵的逆为上三角矩阵，上三角矩阵成上三角矩阵还是上三角矩阵因此$D$是上三角矩阵。又$$I=Q^{H}Q=D^H{Q}_1^H{Q}_{1}D=D^{H}D$$所以$D$是酉矩阵。所以，$D$只能是对角元素绝对值全为1的对角阵

3.求QR分解的办法

3.1Gives变换（利用Gives变换的$Tx=|x|e_1$）

例题：利用$Givens变换$求$A=\left[\begin{array}{ccc}2& 2& 1\\ 0& 2& 2\\ 2& 1& 2\end{array}\right]$的$QR分解$（正交矩阵的转置就是其逆）
解：首先对第一列进行Givens变换，有$$T_{13}=\left[\begin{array}{ccc}c& 0& s\\ 0& 1& 0\\ -s& 0& c\end{array}\right],c=\frac{a_{11}}{\sqrt{a_{11}^2+a_{31}^2}}=\frac{2}{2\sqrt{2}}=\frac{1}{\sqrt{2}},s=\frac{a_{31}}{\sqrt{a_{11}^2+a_{31}^2}}=\frac{1}{\sqrt{2}}$$所以可得$$T_1=\left[\begin{array}{ccc}\frac{1}{\sqrt{2}}& 0& \frac{1}{\sqrt{2}}\\ 0& 1& 0\\ -\frac{1}{\sqrt{2}}& 0& \frac{1}{\sqrt{2}}\end{array}\right],T_{1}A=\left[\begin{array}{ccc}2\sqrt{2}& \frac{3}{\sqrt{2}}& \frac{3}{\sqrt{2}}\\ 0& 2& 2\\ 0& -\frac{1}{\sqrt{2}}& \frac{1}{\sqrt{2}}\end{array}\right],A^{(1)}=\left[\begin{array}{cc}2& 2\\ -\frac{1}{\sqrt{2}}& \frac{1}{\sqrt{2}}\end{array}\right]$$对$A^{(1)}$再进行Givens变换，可以得到$$T_{12}=\left[\begin{array}{cc}c& s\\ -s& c\end{array}\right],c=\frac{2\sqrt{2}}{3},s=-\frac{1}{3}$$所以可得$$T2=\left[\begin{array}{cc}\frac{2\sqrt{2}}{3}& -\frac{1}{3}\\ \frac{1}{3}& \frac{2\sqrt{2}}{3}\end{array}\right],T_2{A}^{(1)}=\left[\begin{array}{cc}\frac{2\sqrt{2}}{3}& \frac{7\sqrt{2}}{6}\\ 0& \frac{4}{3}\end{array}\right],$$综上可得，$$R=\left[\begin{array}{ccc}2\sqrt{2}& \frac{3}{\sqrt{2}}& \frac{3}{\sqrt{2}}\\ 0& \frac{2\sqrt{2}}{3}& \frac{7\sqrt{2}}{6}\\ 0& 0& \frac{4}{3}\end{array}\right],Q={\left[\left[\begin{array}{cc}1& 0\\ 0& T2\end{array}\right]T_1\right]}^T=\left[\begin{array}{ccc}\frac{1}{\sqrt{2}}& \frac{1}{3\sqrt{2}}& -\frac{2}{3}\\ 0& \frac{2\sqrt{2}}{3}& \frac{7\sqrt{2}}{6}\\ \frac{1}{\sqrt{2}}& -\frac{1}{3\sqrt{2}}& \frac{2}{3}\end{array}\right]$$

3.2Householder变换（利用$Hx=|x|z$）

将$z$矩阵设为e矩阵，u=$\frac{x-|x|z}{|x-|x|z|}$
利用Householder求$A=\left[\begin{array}{ccc}0& 4& 1\\ 1& 1& 1\\ 0& 3& 2\end{array}\right]$的$QR$分解
解：首先对第一列使用Householder变换，可以得到$$H=I-2u{u}^T,|x-|x|e_1|=\sqrt{2},x-|x|e_1=\left[\begin{array}{c}-1\\ 1\\ 0\end{array}\right],u=\frac{1}{\sqrt{2}}\left[\begin{array}{c}-1\\ 1\\ 0\end{array}\right]$$ $$2u{u}^T=\frac{1}{\sqrt{2}}\left[\begin{array}{c}-1\\ 1\\ 0\end{array}\right]\frac{1}{\sqrt{2}}\left[\begin{array}{ccc}-1& 1& 0\end{array}\right]=\frac{1}{2}\left[\begin{array}{ccc}1& -1& 0\\ -1& 1& 0\\ 0& 0& 0\end{array}\right]$$ $$I-2u{u}^T=\left[\begin{array}{ccc}0& 1& 0\\ 1& 0& 0\\ 0& 0& 1\end{array}\right]=H,HA=\left[\begin{array}{ccc}1& 1& 1\\ 0& 4& 1\\ 0& 3& 2\end{array}\right]$$继续做下去，可以得到