矩阵分析与应用(17)

Butterfffly

已于 2022-07-29 21:16:15 修改

阅读量2.6k

点赞数 2

分类专栏：矩阵分析与应用文章标签：矩阵

于 2022-07-29 21:15:40 首次发布

本文链接：https://blog.csdn.net/qq_40206924/article/details/126060386

版权

矩阵分析与应用专栏收录该内容

22 篇文章 6 订阅

订阅专栏

学习来源：《矩阵分析与应用》张贤达清华大学出版社

QR 分解及其应用

1. Givens 旋转

Givens 旋转就是初等旋转变换。在 $n$ 维向量空间中，通过固定 $n-2$ 维不变，在剩下的两维构成的平面上做旋转变换。基本结构如下：

$R_{ij}^n=\begin{bmatrix} 1 & & & & & & & & & & \\ & \ddots & & & & & & & & & \\ & & 1 & & & & & & & & \\ & & & cos\theta &0 & \cdots& 0& sin\theta & & & \\ & & &0 &1 & \cdots &0 &0 & & & \\ & & & \vdots&\vdots & &\vdots &\vdots & & & \\ & & & 0&0 &\cdots &1 &0 & & & \\ & & &-sin\theta &0 &\cdots &0 &cos\theta & & & \\ & & & & & & & &1 & & \\ & & & & & & & & &\ddots & \\ & & & & & & & & & &1 \end{bmatrix}$ （1）

矩阵 $R_{ij}^n$ 满足：

$\left\{\begin{matrix} r_{ii}=cos\theta \\ r_{ij}=sin\theta \\ r_{ji}=-sin\theta \\ r_{jj}=cos\theta (i<j)\\ r_{pp}=1(p\neq i,j)\\ r_{pq}=0\quad ,other \end{matrix}\right.$

矩阵 $R_{ij}^n$ 可以看成是 $n$ 阶单位矩阵 $I$ 修改而成：将 $I$ 位于 $(i,i),(ij),(ji),(jj)$ 上的元素分别替换为 $cos\theta ,sin\theta,-sin\theta,cos\theta$ ，其它元素保持不变。其中的 $\theta$ 通常称为旋转角，矩阵 $R_{ij}^n$ 通常叫做初等旋转矩阵，它所确定的线性变换叫做初等线性变换，也可以叫做吉文斯（Givens）变换。

矩阵 $R_{ij}^n$ 具有性质：

1） $detR_{ij}^n=1$ ；

2） $R_{ij}^n$ 对应的初等旋转变换是正交变换， $R_{ij}^n$ 是正交矩阵。

2. 采用 Givens 旋转的 QR 分解

Givens 旋转也可以用来计算 QR 分解。以 $4\times 3$ 矩阵为例，说明 Givens QR 分解的思想：

$\begin{bmatrix} \otimes & \times &\times \\ \times & \times& \times\\ \times & \times & \times\\ \otimes & \times& \times \end{bmatrix}\overset{(1,4)}{\rightarrow}\begin{bmatrix} \otimes & \times &\times \\ \times & \times& \times\\ \otimes & \times & \times\\ 0 & \times& \times \end{bmatrix}\overset{(1,3)}{\rightarrow}\begin{bmatrix} \otimes & \times &\times \\ \otimes & \times& \times\\ 0 & \times & \times\\ 0 & \times& \times \end{bmatrix}\overset{(1,2)}{\rightarrow}\begin{bmatrix} \times & \otimes &\times \\ 0 & \times& \times\\ 0 & \times & \times\\ 0 & \otimes& \times \end{bmatrix}\overset{(1,4)}{\rightarrow}$

$\begin{bmatrix} \times & \otimes &\times \\ 0 & \times &\times \\ 0 & \otimes &\times \\ 0 & 0 & \times \end{bmatrix}\overset{(1,3)}{\rightarrow}\begin{bmatrix} \times & \times &\otimes \\ 0 & \times &\times \\ 0 & 0 &\times \\ 0 & 0 & \otimes \end{bmatrix}\overset{(1,4)}{\rightarrow}\begin{bmatrix} \times & \times &\times \\ 0 & \times &\times \\ 0 & 0 &\times \\ 0 & 0 &0 \end{bmatrix}$