16.投影矩阵,最小二乘

server/2024/12/23 1:08:51/

文章目录

1. 投影矩阵

1.1 投影矩阵P

根据上节知识,我们知道当我们在解 A X = b AX=b AX=b的时候,发现当向量b不在矩阵A的列空间的时候,我们希望的是通过投影,将向量b投影到矩阵A的列空间中,这样,我们可以求得一个近似的解,得到如下公式
A T A X ^ = A T b (1) A^TA\hat{X} = A^Tb\tag{1} ATAX^=ATb(1)

  • 我们假设 A T A 可逆, A^TA可逆, ATA可逆,可得到解为:
    X ^ = ( A T A ) − 1 A T b (2) \hat{X}=(A^TA)^{-1}A^Tb\tag{2} X^=(ATA)1ATb(2)
  • 那么可以得到向量b在矩阵A的列空间向量p表示如下:
    p = A ( A T A ) − 1 A T b (3) p=A(A^TA)^{-1}A^Tb\tag{3} p=A(ATA)1ATb(3)
  • 由上可以看出,我们将矩阵 P = A ( A T A ) − 1 A T P=A(A^TA)^{-1}A^T P=A(ATA)1AT代入可得:
    p = P b (4) p=Pb\tag{4} p=Pb(4)
  • 我们发现,向量b为不在矩阵A的列空间中的向量,p为向量b通过投影矩阵P转换后的向量。并且向量p是在矩阵A的列空间中。

1.2 投影向量

对于任意向量b来说,我们可以通过正交分解,将向量b分解到两个垂直的向量空间中,我们考虑两个极端的情况下

  • 假设向量b在矩阵A的列空间中,那么向量b通过投影矩阵P的转换,还是得到其本身
    P b = b (5) Pb=b\tag{5} Pb=b(5)
  • 假设向量b在垂直于矩阵A的列空间中,那么向量b通过投影矩阵P的转换,得到的将是零向量
    P b = 0 (6) Pb=0\tag{6} Pb=0(6)
    那么我们思考下,什么向量空间是垂直于矩阵A的列空间的呢?我们之前学过矩阵A的四个子空间,分别是
  1. Row(A) —> 矩阵A的行空间;2.Colum(A) —> 矩阵A的列空间
  2. N(A) —> 矩阵A的零解空间;4. N ( A T ) N(A^T) N(AT) —> 矩阵 A T A^T AT的零解空间
    我们可以将 A T A^T AT按列向量拆解得到如下
    A T = [ a 1 T a 2 T ⋮ a n T ] ; [ a 1 T a 2 T ⋮ a n T ] [ y 1 y 2 … y n ] = 0 (7) A^T=\begin{bmatrix}a_1^T\\\\a_2^T\\\\\vdots\\\\a_n^T \end{bmatrix};\begin{bmatrix}a_1^T\\\\a_2^T\\\\\vdots\\\\a_n^T \end{bmatrix}\begin{bmatrix}y_1&y_2&\dots&y_n\end{bmatrix}=0\tag{7} AT= a1Ta2TanT ; a1Ta2TanT [y1y2yn]=0(7)
  • 由上述可以看出, A T A^T AT的零解空间是垂直于矩阵A的列空间的,所以我们可以将任意向量b 通过正交分解为一部分投影在列空间的向量p,另一部分投影在 A T A^T AT的零解空间中的e
    在这里插入图片描述
    p = P b e = ( I − P ) b (8) p=Pb\\\\e=(I-P)b\tag{8} p=Pbe=(IP)b(8)

http://www.ppmy.cn/server/42170.html

相关文章

机器人SCI期刊,中科院3区,收稿范围广泛!

一、期刊名称 Journal of Intelligent & Robotic Systems 二、期刊简介概况 期刊类型:SCI 学科领域:计算机科学 影响因子:3.3 中科院分区:3区 出版方式:开放出版 版面费:$2990 三、期刊征稿范围…

java 的类型转换

在 Java 中,类型转换可以分为两类:隐式类型转换(自动类型转换)和显式类型转换(强制类型转换)。 一、隐式类型转换 隐式类型转换是指在编译时自动进行的类型转换。这种转换通常发生在基本数据类型之间&#…

DHCP动态主机配置协议

DHCP概述 DHCP是什么 DHCP:Dynamic Host Configuration Protocol:动态主机配置协议DHCP是一种集中对用户IP地址进行动态管理和配置的技术 DHCP作用: 作用:实现IP地址的动态分配和集中管理优势:避免手工配置IP地址&…

网络学习(二)|网络标识符URI与URL的区别与联系

文章目录 URIURL区别与联系示例结论 在工作中,我们经常会听到URI和URL这两个术语,它们是网络中不可或缺的概念。尽管它们常常被混淆使用,但它们代表着不同的概念和角色。我们将深入探讨URI和URL,帮助彻底理解它们之间的区别和联系…

【二叉树算法题记录】从中序与后序遍历序列构造二叉树,从前序与中序遍历序列构造二叉树

目录 从中序与后序遍历序列构造二叉树题目描述题目分析cpp代码 从前序与中序遍历序列构造二叉树题目描述题目分析cpp代码 从中序与后序遍历序列构造二叉树 题目链接 题目描述 给定两个整数数组 inorder 和 postorder ,其中 inorder 是二叉树的中序遍历&#xff0…

linux编译gdb

下载 我下载了8.3 Index of /gnu/gdb 编译 make cleanmake diskclean./configure \--without-x \--disable-werrormake make install

【SQL】SQL常见面试题总结(2)

目录 1、增删改操作1.1、插入记录(一)1.2、插入记录(二)1.3、插入记录(三)1.4、更新记录(一)1.5、更新记录(二)1.6、删除记录(一)1.7、…

Python学习之路 | Python基础语法(一)

数据类型 Python3 中常见的数据类型有: Number(数字)String(字符串)bool(布尔类型)List(列表)Tuple(元组)Set(集合)Dict…