找回密码
 快速注册
搜索
查看: 35|回复: 1

正交投影矩阵

[复制链接]

3149

主题

8388

回帖

6万

积分

$\style{scale:11;fill:#eff}꩜$

积分
65401
QQ

显示全部楼层

hbghlyj 发表于 2023-3-3 09:48 |阅读模式
Projection #formulas
点到1维线性子空间(过原点的直线)的正交投影矩阵(外积)$P_u=uu^T$, $\{u\}$是正交规范基(即$u$为直线的单位方向向量). 当$\{u\}$不必须是正交规范基(即长度不为1时)公式$P_u=u(u^Tu)^{-1}u^T$也适用, 这里$(u^Tu)^{-1}=\|u\|^{-2}$是把$u$规范化为$u\over\|u\|$的因子.
例子
点$(3,7,1)^T$ 到 直线$x_1-x_2+3x_3 = 2x_1+ x_2 = 0$ 的正交投影.
方向向量 $(1, -1, 3)^T\times(2, 1, 0)^T=(-3, 6, 3)^T$
单位方向向量 $u=\frac1{\sqrt6}(-1, 2, 1)^T$
$P_u(3,7,1)^T=\frac1{6}(-1, 2, 1)^T(-1, 2, 1)(3,7,1)^T=(-2,4,2)^T$

在任意维线性子空间上的正交投影。设 $\{u_1, …, u_k\}$ 是子空间 $U$ 的一组正交规范基,$A$ 是以这些向量作为列的矩阵,那么投影 $P_{A}=AA^{T}$ 是正交的。
当 $\{u_1, …, u_k\}$ 不必须是正交规范基,公式 $P_{A}=A(A^{T}A)^{-1}A^{T}$ 也适用。这里$(A^{T}A)^{-1}$是规范化因子。
例子
点$(3,-1,1,-1)^T$ 到 $V=\left\{\mathbf{x} \in \mathbb{R}^{4}: x_{1}+x_{2}+x_{3}+x_{4}=x_{1}+2 x_{3}+x_{4}=0\right\}$ 的投影
解线性方程组$\cases{ x_{1}+x_{2}+x_{3}+x_{4}=0\\x_{1}+2 x_{3}+x_{4}=0}$得 $V$ 的一个基为$\{(-1, 0, 0, 1)^T, (-2, 1, 1, 0)^T\}$
则 $A=\pmatrix{-1&-2\\0&1\\0&1\\1&0}$
则 $P_A=\frac14\left(
\begin{array}{cccc}
3 & -1 & -1 & -1 \\
-1 & 1 & 1 & -1 \\
-1 & 1 & 1 & -1 \\
-1 & -1 & -1 & 3 \\
\end{array}
\right)$
则 $P_A\pmatrix{3\\-1\\1\\-1}=\frac12\left(
\begin{array}{c}
5 \\
-1 \\
-1 \\
-3 \\
\end{array}
\right)$

3149

主题

8388

回帖

6万

积分

$\style{scale:11;fill:#eff}꩜$

积分
65401
QQ

显示全部楼层

 楼主| hbghlyj 发表于 2023-3-3 10:19
点 $(1,1,1,1)^T$ 到 $\mathbb R^4$ 的子空间 $W=\left<\frac{1}{\sqrt2}(1,-1,0,0)^T,\frac{1}{\sqrt2}(0,0,1,1)^T\right>$ 的距离?
Solution.
设$\left\{v_1,v_2\right\}$为$W$的规范正交基, 则$P=\pmatrix{v_1&v_2}\pmatrix{v_1&v_2}^T=v_1{v_1}^T+v_2{v_2}^T$是到$W$的投影矩阵, 如$A_W=PA$为$A$到$W$的投影.
[证明: 容易验证$P=P^2$, 所以$(P-I)^TP=0\implies(A_W-A)^TA_W=A^T(P-I)^TPA=0$, 即$A_W-A$与$A_W$正交.]
$A$到$W$的距离等于$\|AA_W\|$.
代入数字, $A_W={\begin{bmatrix}
0.5 & -0.5 & 0 & 0 \\
-0.5 & 0.5 & 0 & 0 \\
0 & 0 & 0.5 & 0.5 \\
0 & 0 & 0.5 & 0.5 \\
\end{bmatrix}} {\begin{bmatrix}
1 \\
1 \\
1 \\
1 \\
\end{bmatrix}}={\begin{bmatrix}
0 \\
0 \\
1 \\
1 \\
\end{bmatrix}}$

手机版|悠闲数学娱乐论坛(第3版)

GMT+8, 2025-3-4 16:28

Powered by Discuz!

× 快速回复 返回顶部 返回列表