[矩阵的QR分解系列二] 吉文斯(Givens)变换

之前介绍的矩阵的三角分解系列介绍了利用矩阵初等变换解决了矩阵三角化问题以及具体的三角分解。但是以初等变换工具的三角分解方法并不能消除病态线性方程组不稳定问题，而且有时候对于可逆矩阵有可能也不存在三角分解。所以后面为了解决这里问题，发展出来了以正交(酉)变换的矩阵的QR(正交三角)分解，矩阵的正交三角分解是一种对任何可逆矩阵均存在理想分解。进行QR分解需要用到施密特(Schmidt)正交规范化，吉文斯(Givens)变换和豪斯霍尔德(Householder)变换等。这里矩阵的QR分解系列教程主要是针对在学习QR分解时候的涉及到的一些细节，包括很多方法的来源和证明等，以及其中用到的一些矩阵操作的基础知识，主要包括：

• [矩阵的QR分解系列一] 施密特(Schmidt)正交规范化
• [矩阵的QR分解系列二] 吉文斯(Givens)变换
• [矩阵的QR分解系列三] 豪斯霍尔德(Householder)变换
• [矩阵的QR分解系列四] QR(正交三角)分解
• [矩阵的QR分解系列五] Eigen中的QR分解

这个系列后面文章会用到前面文章的理论和技术，所以建议按照顺序查看。

简介

初等旋转变换或者说是吉文斯(Givens)变换是一种正交变换，经过多次吉文斯(Givens)变换可以把矩阵转换成上三角形式，是一种常用的$\mathbf{Q}\mathbf{R}$分解方式。

初等旋转变换

定义

在平面${\mathbf{R}}^2$中旋转变换对应的矩阵为
$$\mathbf{A}=\left[\begin{array}{cc}cos\theta & sin\theta \\ -sin\theta & cos\theta \end{array}\right].$$
一般在$n$维欧式空间$\mathbf{V}$中取一组标准正交基$e_1,e_2,\cdots ,e_n$,沿平面$[e_i,e_j]$旋转，它的旋转矩阵为
$${\mathbf{R}}_{ij}=\left[\begin{array}{ccccccccccc}1& & & & & & & & & & \\ & \ddots & & & & & & & & & \\ & & 1& & & & & & & & \\ & & & cos\theta & 0& \cdots & 0& sin\theta & & & \\ & & & 0& 1& \cdots & 0& 0& & & \\ & & & ⋮& ⋮& & ⋮& ⋮& & & \\ & & & 0& 0& \cdots & 1& 0& & & \\ & & & -sin\theta & 0& \cdots & 0& cos\theta & & & \\ & & & & & & & & 1& & \\ & & & & & & & & & \ddots & \\ & & & & & & & & & & 1\end{array}\right].$$
这里矩阵${\mathbf{R}}_{ij}$里面元素满足：$r_{ii}=cos\theta ,r_{ij}=sin\theta ,r_{ji}=-sin\theta ,r_{jj}=cos\theta (i<j),r_{pp}=1(p\ne i,j)$，除此之外其余元素$r_{pq}=0$。可以看成是$n$阶单位矩阵$\mathbf{I}$修改而成：将$\mathbf{I}$位于$(i,i),(i,j),(j,i),(j,j)$上的元素分别替换成$cos\theta ,sin\theta ,-sin\theta ,cos\theta$，其余元素保持不变。式中的$\theta$通常称为旋转角，而$cos\theta ,sin\theta$常分别记为$c,s$，所以满足$c^2+s^2=1$，所以矩阵${\mathbf{R}}_{ij}$可以简记为
$${\mathbf{R}}_{ij}=[\begin{array}{cccc}1& & & \\ & \ddots & & \\ & & 1& \\ & & & c\\ & & & 0\\ & & & ⋮\\ & & & 0\\ & & & -\\ & & & \\ & & & \\ & & & \end{array}$$