icon

Sealos,在云桌面中运行分布式应用程序,像使用个人电脑一样使用云

icon

去看看👀

icon

扫码加入微信群,和云原生大佬们一起探讨云原生和不可描述的事情!

wechat qr code
 
今日天气
 
舔狗日记

子空间投影
  1. 博客/

子空间投影

·480 字·1 分钟· · ·
数学 向量 矩阵
米开朗基杨
作者
米开朗基杨
云原生搬砖师 & Sealos 开发者布道师 & FastGPT 熟练工
Table of Contents
gptgod
FastGPT
Laf
Contact me

为了弄明白子空间投影是怎么一回事,我们遵循从低维到高维的规律,先从二维开始讲起。

1. 二维空间
#


如下图所示(我随手画的,不要介意),设向量p是向量b在向量a上面的投影,向量e垂直于向量p及a。

图片描述: 1.png

于是我们可以得到这样的一个等式:

$$ a^Te = 0 $$

$$a^T(b-xa) = 0\tag{1}$$

解得:

$$x = \frac{a^Tb}{a^Ta}$$

于是向量p可表示为:

$$p = xa = a\frac{a^Tb}{a^Ta}\tag{2}$$

现在我们设

$$p = Pb\tag{3}$$

我们把这个矩阵P称为投影矩阵

比较式(2)和式(3),立即可以知道:

$${P = \frac{a \cdot a^T}{a^T \cdot a}}\tag{4}$$

2. 三维空间
#


为了让你们能够有一个直观的认识,我仍然用我高超的画艺画了一幅美图:

图片描述: 3.png

假设图中的那个平面由向量$a_1$和$a_2$构成,令

$$A = \begin{bmatrix} a_1 & a_2 \end{bmatrix}$$

由于向量p在平面上,所以p可以表示为:

$$p = \hat{x_1}a_1 + \hat{x_2}a_2\tag{5}$$

$$p = A\hat{x}\tag{6}$$

与二维空间类似,设向量p是向量b在平面上的投影,向量e垂直于那个平面,当然也垂直于向量p,同样也垂直于向量$a_1$和$a_2$,于是可以得到方程组:

$$ \begin{cases} a_1^T(b - A\hat{x}) = 0 \\ a_2^T(b - A\hat{x}) = 0 \\ \end{cases}\tag{7} $$

$$ \begin{bmatrix} a_1^T \\a_2^T \end{bmatrix}(b - A\hat{x}) = \begin{bmatrix} 0 \\0 \end{bmatrix} $$

进一步化简得到:

$$A^T(b - A\hat{x}) = 0\tag{8}$$

解得:

$$\hat{x} = (A^TA)^{-1}(A^Tb)\tag{9}$$

将(9)代入(6)得:

$$p = A\hat{x} = A(A^TA)^{-1}A^Tb\tag{10}$$

与二维空间类似,我们设

$$p = Pb\tag{11}$$

比较式(10)和式(11),立即可以得到:

$${P = A(A^TA)^{-1}A^T}\tag{12}$$

这就是投影矩阵的表达式!

-------他日江湖相逢 再当杯酒言欢-------

相关文章

最小二乘法的本质
·1262 字·3 分钟·
数学 向量 矩阵

公众号二维码