首页 > 其他 > 详细

矩阵QR分解

时间:2020-07-13 17:53:35      阅读:54      评论:0      收藏:0      [点我收藏+]

1 orthonormal 向量与 Orthogonal 矩阵

    orthonormal 向量定义为 技术分享图片,任意向量 技术分享图片 相互垂直,且模长为1;

    如果将  orthonormal 向量按列组织成矩阵,矩阵为 Orthogonal 矩阵,满足如下性质:

    技术分享图片

   当 技术分享图片 为方阵时,技术分享图片为其逆矩阵;当 技术分享图片 为长方形矩阵时,技术分享图片为其左逆;

   当矩阵 Q 为正交矩阵时,对向量变换变换前后点积不发生改变,技术分享图片,证明如下:

   技术分享图片,当 x = y 时,有 技术分享图片

   对任意向量 b ,可以分解为一组正交向量的线性组合,技术分享图片,要求解系数x,可先写成矩阵形式:

   技术分享图片技术分享图片技术分享图片

   因此,向量 b 可分解为 技术分享图片

 

2 Gram-Schmidt 与 QR 分解

   对矩阵 技术分享图片,可以将其转换为正交矩阵 技术分享图片,方法如下:

   1)向量 技术分享图片 方向保持不变,将其长度归一化, 技术分享图片

   2)向量 技术分享图片 可分解为向量 技术分享图片 投影分量与垂直于向量  技术分享图片 的两分量,剔除投影分量得 技术分享图片技术分享图片

   3)同理,剔除向量 技术分享图片 在  技术分享图片, 技术分享图片 上投影分量得 技术分享图片技术分享图片

   4)依照如上方法,可以对所有向量完成正交化。

   以上处理可以使用矩阵表示,矩阵 Q 为矩阵 A 的列进行线性变换结果,故可写为 A=QR 。

   1)向量  技术分享图片 与向量 技术分享图片 具有相同方向,故可表示为 技术分享图片

   2)向量 技术分享图片 被分解为 技术分享图片 方向向量,可表示为 技术分享图片

   3)向量 技术分享图片 被分解为 技术分享图片 方向向量,可表示为 技术分享图片

   4)综上表示为矩阵形式 技术分享图片

 

3 求解 Ax=b

   使用 Gram-Schmidt 可将矩阵 A 转换为正交矩阵 Q,正交矩阵 Q 可简化 Ax=b 运算:

   1)最小二乘法求解 技术分享图片

   2)带入 技术分享图片 得 技术分享图片,化简得 技术分享图片

   3)不管长方形矩阵还是方阵,都有 技术分享图片,故上式可化简为 技术分享图片

   4)由于 R 为上三角矩阵,使用回代法即可求解。

 

    参考资料 Linear Algebra And Its Applications   Gilbert Strang

矩阵QR分解

原文:https://www.cnblogs.com/luofeiju/p/13294501.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!