当前位置: 首页 > news >正文

复矩阵的QR分解

复矩阵的QR分解

定义:QR分解

\(A\) 是一个 \(m \times n\) 复矩阵,且 \(m \geq n\)。如果存在一个 \(m \times r\) 酉矩阵 \(Q\) 和一个 \(r \times r\) 上三角矩阵 \(R\),使得

\[A = QR \]

则称此分解为 \(A\) 的 QR 分解,其中 \(r=\mathrm{rank}(A)\)。(当 \(Q\)\(m \times n\) 矩阵且满足 \(Q^*Q = I_n\) 时,称为经济型 QR 分解。)

定理:列满秩矩阵QR分解的存在性与唯一性

任意列满秩复矩阵 \(A \in \mathbb{C}^{m \times n}\) 都存在一个 \(m \times n\) 酉矩阵 \(Q\) 和一个 \(n \times n\) 上三角矩阵 \(R\),使得

\[A = QR \]

如果进一步要求 \(R\) 的对角元为正实数,则该分解是唯一的。

证明:存在性

我们使用 Gram-Schmidt 正交化过程来构造证明。设 \(A = [a_1, a_2, \dots, a_n]\),其中 \(a_i \in \mathbb{C}^m\)\(A\) 的列向量。

  1. 正交化过程
    定义:

    • \(u_1 = a_1\)
    • \(q_1 = \dfrac{u_1}{\|u_1\|}\)
    • 对于 \(k = 2, 3, \dots, n\)

      \[u_k = a_k - \sum_{j=1}^{k-1} \langle a_k, q_j \rangle q_j \]

      \[q_k = \dfrac{u_k}{\|u_k\|} \quad (\text{如果 } u_k \neq 0) \]

    其中 \(\langle x, y \rangle = y^\dagger x\) 是复向量空间中的内积。

  2. 构造 QR 分解
    \(Q = [q_1, q_2, \dots, q_n]\),则 \(Q^\dagger Q = I_n\)
    定义上三角矩阵 \(R\) 的元素为:

    \[r_{ij} = \begin{cases} \langle a_j, q_i \rangle, & \text{如果 } i \leq j \\ 0, & \text{如果 } i > j \end{cases}\]

    特别地,\(r_{kk} = \|u_k\|\)

  3. 验证分解
    对于每个 \(k = 1, 2, \dots, n\),有:

    \[a_k = \sum_{i=1}^k r_{ik}q_i = \sum_{i=1}^k \langle a_k, q_i \rangle q_i \]

    因此:

    \[A = [a_1, a_2, \dots, a_n] = [q_1, q_2, \dots, q_n]R = QR \]

证明:唯一性

假设有两个不同的分解 \(A = Q_1 R_1 = Q_2 R_2\),那么:

\[Q_2^* Q_1 = R_2 R_1^{-1} \]

  • 左边是酉矩阵的乘积,仍然是酉矩阵;
  • 右边是两个上三角矩阵的乘积,仍然是上三角矩阵;
  • 一个既是酉矩阵又是上三角矩阵的矩阵必须是对角矩阵;
  • 由于 \(R_1, R_2\) 对角线都是正实数,且乘积为单位矩阵,这个对角矩阵只能是单位矩阵。

因此 \(Q_1 = Q_2, R_1 = R_2\)

定理:一般的复矩阵的QR分解

对于任意复矩阵 \(A \in \mathbb{C}^{m \times n}\),存在一个 \(m \times r\) 酉矩阵 \(Q\),一个 \(r \times n\) 上梯形矩阵 \(R\),一个置换矩阵 \(P \in \mathbb{C}^{n \times n}\),使得

\[AP = QR \]

其中,\(r = \mathrm{rank}(A) \leq n\)

证明

我们用置换矩阵 \(P \in \mathbb{C}^{n \times n}\) 作用于 \(A\) 的列向量组,使得 \(AP\) 的前 \(r\) 个向量线性无关。不妨设 \(AP = [a_1, a_2, \dots, a_n]\),则 \(a_1, \dots, a_r\) 线性无关。

\([a_1, \dots, a_r] = A_r \in \mathbb{C}^{m \times r}\) 列满秩,则根据上一个定理,存在一个 \(m \times r\) 酉矩阵 \(Q_r\) 和一个 \(r \times r\) 上三角矩阵 \(R_r\),使得

\[A_r = Q_r R_r \]

对于 \(k = r + 1, \dots, n\),存在 \(x_k \in \mathbb{C}^r\) 使得 \(a_k = A_r x_k = Q_r R_r x_k\),记 \(y_k = R_r x_k \in \mathbb{C}^r\),令 \(R_r^\prime = [y_{r+1}, \dots, y_n] \in \mathbb{C}^{r \times (n-r)}\),则

\[AP = Q_r \begin{bmatrix} R_r & R_r^\prime \end{bmatrix} := QR \]

其中 \(R\) 为上梯形矩阵。

注记

上面的定理可以进一步扩充:\(Q\) 的列向量可以扩充为一组单位正交基得到 \(\overline{Q} = \begin{bmatrix} Q & Q_{\perp} \end{bmatrix}\)\(m \times m\) 酉矩阵,于是相应的 \(R\) 变成了

\[\overline{R} = \begin{bmatrix} R \\ 0 \end{bmatrix} \]

因此,

\[A = \begin{bmatrix} Q & Q_{\perp} \end{bmatrix} \begin{bmatrix} R \\ 0 \end{bmatrix} \]

这成为完全 QR 分解,但是这种分解是不唯一的。

示例:秩亏缺矩阵的QR分解

考虑秩为1的矩阵:

\[A = \begin{bmatrix} 1 & 2 & 3 \\ 1 & 2 & 3 \\ 1 & 2 & 3 \end{bmatrix} \]

一种可能的 QR 分解为:

\[Q = \begin{bmatrix} \frac{1}{\sqrt{3}} & \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{6}} \\ \frac{1}{\sqrt{3}} & -\frac{1}{\sqrt{2}} & \frac{1}{\sqrt{6}} \\ \frac{1}{\sqrt{3}} & 0 & -\frac{2}{\sqrt{6}} \end{bmatrix}, \quad R = \begin{bmatrix} \sqrt{3} & 2\sqrt{3} & 3\sqrt{3} \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{bmatrix} \]

另一种分解为:

\[Q = \begin{bmatrix} \frac{1}{\sqrt{3}} & \frac{1}{\sqrt{2}} & 0 \\ \frac{1}{\sqrt{3}} & 0 & \frac{1}{\sqrt{2}} \\ \frac{1}{\sqrt{3}} & -\frac{1}{\sqrt{2}} & -\frac{1}{\sqrt{2}} \end{bmatrix}, \quad R = \begin{bmatrix} \sqrt{3} & 2\sqrt{3} & 3\sqrt{3} \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{bmatrix} \]

两者都满足 \(A = QR\),说明分解不唯一。

注记

完全 QR 分解不唯一的根本原因是:\(Q_{\perp}\) 不唯一,其列向量可以随意换位置。

http://www.hskmm.com/?act=detail&tid=36682

相关文章:

  • 权威调研榜单:天津全屋定制整体橱柜方案TOP4榜单好评深度解析
  • 别再手动处理琐事了!用Coze搭建AI工作流,我每天白赚2小时
  • 单时段机组组合优化的粒子群算法实现(MATLAB)
  • 谎言 欺骗 鄙夷 如破碎瓦砾铺满地 利用陷害窒息莫名遭受唾骂遗弃
  • Day21-C:\Users\Lenovo\Desktop\note\code\JavaSE\Basic\src\com\Stream-集合框架(stream)
  • 权威调研榜单:湖南张家界旅游团服务TOP3榜单好评深度解析
  • [GXYCTF2019]Ping Ping Ping 1
  • C# 元组 Tuple ValueTuple
  • Java语言的核心特性与大数据应用研究
  • 扣子Coze智能体万字教程:从入门到精通,一文掌握AI工作流搭建 - Leone
  • 轻量服务器Lighthouse + 1Panel + Halo,三步打造你的专属网站
  • (自用)如何使用 mt19937 生成随机数?
  • 第四章 windows实战-向日葵
  • 第四章 windows实战-emlog
  • 第四章 windows实战-wordpress
  • Docling + LangChain + RAG 构建智能文档问答系统
  • 广义串并联图
  • 2025 种植棚/养殖棚/工程/羊肚菌/保温/园林/加厚/绿化/草苫子推荐榜:济宁泽萌草制品 5 星领跑,适配大棚 / 混凝土 / 园艺多场景需求
  • 2025年口碑好的防锈漆厂家最新推荐排行榜
  • 2025 餐饮/电商/品牌策划优选榜:上海物心 5 星领跑,细分领域这些适配型企业值得 pick
  • 2025 灌装/大桶/桶装/纯净/瓶装/水设备精选推荐榜:路得自动化领衔,这些靠谱品牌值得关注
  • 2025 打捆机捆/包装/绕树干/草绳推荐榜:济宁泽萌草制品 5 星领跑,适配农业 / 园艺 / 建筑 / 物流捆扎需求
  • 遇到一例无法保存编辑后的 Excel 文件的问题
  • 告别难找的本地搜图工具!ImageSearch 2.2 携全新 2.0UI 登场
  • 2025 年焊管源头厂家最新推荐榜单:盘点高性能焊管制造企业,助力下游企业精准选品家具/汽车/高精度内刮/镀铝/不锈钢/高强钢焊管厂家推荐
  • 2025 年钢管厂家最新推荐榜:覆盖精密钢管、汽车钢管、高强钢钢管等品类,为下游采购企业提供权威选品参考
  • ubuntu24.04 server 版本安装xfce 使用web novnc 远程桌面
  • 数据结构做题记录
  • 生成函数入门
  • 2025 年 PPT 生成工具品牌最新推荐榜单:AI 备案技术加持 + 千万用户信赖,优质工具全方位测评ppt一键生成/自动生成ppt/ppt制作ai工具推荐