\(\DeclarePairedDelimiterX{\Set}[2]{\{}{\}}{#1 \nonscript\;\delimsize\vert\nonscript\; #2}\) \( \DeclarePairedDelimiter{\set}{\{}{\}}\) \( \DeclarePairedDelimiter{\parens}{\left(}{\right)}\) \(\DeclarePairedDelimiterX{\innerproduct}[1]{\langle}{\rangle}{#1}\) \(\newcommand{\ip}[1]{\innerproduct{#1}}\) \(\newcommand{\bmat}[1]{\left[\hspace{2.0pt}\begin{matrix}#1\end{matrix}\hspace{2.0pt}\right]}\) \(\newcommand{\barray}[1]{\left[\hspace{2.0pt}\begin{matrix}#1\end{matrix}\hspace{2.0pt}\right]}\) \(\newcommand{\mat}[1]{\begin{matrix}#1\end{matrix}}\) \(\newcommand{\pmat}[1]{\begin{pmatrix}#1\end{pmatrix}}\) \(\newcommand{\mathword}[1]{\mathop{\textup{#1}}}\)
Needs:
Real Affine Set Representations
Real Vector Projections
Projections On Subspaces
Real Matrix Range
Needed by:
Minimum Residual Affine Sets
Minimum Residual Affine Sets
Links:
Sheet PDF
Graph PDF

Projections On Affine Sets

Why

What is the projection of a vector onto an affine set?

Result

Suppose $a \in \R ^n$ and $U \in \R ^{n \times k}$ with $U^\top U = I$. Define the affine set $W(a, U) = a + \range(U)$. Then

\[ \textstyle \proj_{W(a, U)}(x) = UU^\top x + (I - UU^\top )a. \]

The minimizer of $J: \R ^n \to \R $ defined by

\[ J(z) = \norm{a + Uz - x} = \norm{Uz - (x -a)}, \]

is $z^\star = U^\top (x-a)$.
So the projection of $x$ onto $W(a, U)$ is

\[ a + Uz^\star = a + UU^\top (x - a) = UU^\top x + (I - UU^\top )a. \]

Copyright © 2023 The Bourbaki Authors — All rights reserved — Version 13a6779cc About Show the old page view