《矩阵论》学习笔记（四）：4.2 矩阵的QR分解

最新推荐文章于 2024-07-14 00:28:22 发布

熊宝宝爱学习

最新推荐文章于 2024-07-14 00:28:22 发布

阅读量3k

点赞数 3

分类专栏：数学文章标签：线性代数矩阵

本文链接：https://blog.csdn.net/sinat_36831051/article/details/105854179

版权

18 篇文章 12 订阅

订阅专栏

-	Given矩阵与Given变换
几何意义	描述：n维欧氏空间 $R^n$ 中的旋转变换，旋转不改变向量模长
Givens矩阵(初等旋转矩阵): $T_{ij}/ T_{ij}(c,s)$	$T_{ij}= \left[\begin{array}{cccc} I& & & \\&c& ...&s& \\ & &I & \\ &-s& ...&c&\\ && & & I\end{array} \right] ，$ 其中： $c^2+ s^2=1$
$T_{ij}$ 的性质	1） $T_{ij}$ 是正交矩阵 $T_{ij}(c,s)]^{-1}=[T_{ij}(c,s)]^T=[T_{ij}(c,-s)]$ ； 2） $det(T_{ij})=1$ ； 3） $\vec x=(a_1,a_2,...,a_n)，\vec y=T_{ij} \vec x=(b_1,b_2,...,b_n)，$ 则： $b_i=ca_i+sa_j；b_j=-sa_i+ca_j；b_k=a_k(k≠i,j)$
Givens变换定理	任意 $\vec x$ 总可以经过有限个Given矩阵的乘积 $T$ ，到某一方向下的投影， $T\vec x=x模长*\vec z$ ， $\vec z$ 为单位列向量。

-Givens矩阵求解过程：p141

根据Givens变换定理：
> 题目：给定待求 $\vec x=(x_1,x_2,...,x_n)$ 和单位方向向量 $\vec z$ ，求Givens变换矩阵T。
> 求解：对 $\vec x$ 重复操作- 将每个维度的分量逐一合并到第一维度上。从而得到一系列Givens变换后的结果 $T\vec x=\sqrt[]{x_1^2+x_2^2+...+x_n^2}$ . 其中， $T=T_{12}T_{13}...*T_{1n}，(i=1,j∈[2,n])$

-	Householder矩阵与Householder变换
几何意义	描述：n维欧氏空间 $R^n$ 中的反射变换
Householder矩阵 (初等反射矩阵)	$H=I-2\vec u\vec u^T$
$H$ 的性质	1） $H$ 是正交矩阵， $H^TH=I、H=H^T$ ； 2） $d e t (H) = - 1$ ； 3） $\vec y=H\vec x=(I-2\vec u\vec u^T) \vec x$
Householder变换定理	任意 $\vec x$ 总可以经过Householder矩阵 $H$ ，得到关于 $\vec z$ 的正交向量对称的向量， $H\vec x=[x]\vec z，\vec z$ 为单位列向量。

-Householder矩阵求解过程：p143

根据Householder变换定理：
> 题目：给定待求 $\vec x=(x_1,x_2,...,x_n)$ 和单位方向向量 $\vec z$ ，求Givens变换矩阵T。
> 求解： 1- 求 $\vec x$ 模长; 2- 求与 $\vec z$ 正交的 $\vec u$ 的方向; 3- 求 $\vec u$ (单位向量:模长为1); 4- 求得 $H=I-2\vec u\vec u^T$ .

-	矩阵QR分解要点汇总
QR分解的前提	可逆方阵
应用到的知识	正交矩阵的性质
QR分解的定义	对实可逆、方阵A，有： $A = Q R$ ，称作矩阵A的QR分解。其中，Q是正交矩阵( $Q^TQ=QQ^T=I$ )，R是可逆上三角矩阵。
- 存在性	若A满足可逆、方阵，则一定存在。
- 唯一性	A的QR分解不唯一，但有紧密联系。各种分解相差一个对角元素绝对值(模)全等于1的对角矩阵因子。

对任意矩阵 $A_{m*n}$ ，其QR分解一定存在，有：

矩阵的QR分解方法	具体步骤
Schmidt正交化方法	1）求矩阵A列向量的标准正交向量; 2）将矩阵A表示成 $A = Q R$ 的形式; 3）求得正交矩阵Q和下三角阵R。
Given变换方法	A左乘有限个初等旋转变换矩阵到可逆上三角阵R。 1） $R = T * A$ [R的第i行为 $T_iA^{(i-1)}$ 的第一行(i≤n-1)，R的第n行最后一个元素= $T_{n-1}A^{(n-2)}$ 的最后一个元素]; 2） $Q=T^{-1}=T^T$ . -若A为不可逆矩阵， $\vec b$ 中出现0，则跳过这一步进行下一步，得到的R也是不可逆的。
Householder变换方法	A左乘有限个初等反射变换矩阵到可逆上三角阵。 1） $R = S * A$ [R的第i行为 $S_iA^{(i-1)}$ 的第一行(i≤n-1)，R的第n行最后一个元素= $S_{n-1}A^{(n-2)}$ 的最后一个元素]; 2） $Q=S^{-1}=S^T$ . -若A为不可逆矩阵， $\vec b$ 中出现0，则跳过这一步进行下一步，得到的R也是不可逆的。

对于可逆方阵A的存在性与唯一性，有：

LU/三角分解	1- 若满足n个顺序主子式≠0 ⇔三角分解存在且唯一 [充要条件]。 2- 若三角分解不存在，即不满足n个顺序主子式≠0⇔存在P使PA的所有顺序主子式≠0。即A虽不存在，PA存在LU分解。
QR分解	1- 存在性：一定存在。 2-唯一性：A的QR分解不唯一，但有紧密联系。各种分解相差一个对角元素绝对值(模)全等于1的对角矩阵因子。