The ProjectionMatrix(S) command constructs the matrix of the orthogonal linear projection onto the subspace spanned by the vectors in S.  If $B$ is a maximal, linearly independent subset of $S$ and $M$ is the Matrix whose columns are the Vectors in $B$, then

Additional arguments are passed as options to the Matrix constructor which builds the result.

$\mathrm{with}\left(\mathrm{Student}\left[\mathrm{LinearAlgebra}\right]\right)\:$

$S≔\left[⟨1\,2\,3\,4⟩\,⟨4\,3\,2\,1⟩\right]$

$S≔\left[\left[\begin{array}{c}1\\ 2\\ 3\\ 4\end{array}\right]\,\left[\begin{array}{c}4\\ 3\\ 2\\ 1\end{array}\right]\right]$

$P≔\mathrm{ProjectionMatrix}\left(S\right)$

$P≔\left[\begin{array}{cccc}\frac{7}{10}& \frac{2}{5}& \frac{1}{10}& -\frac{1}{5}\\ \frac{2}{5}& \frac{3}{10}& \frac{1}{5}& \frac{1}{10}\\ \frac{1}{10}& \frac{1}{5}& \frac{3}{10}& \frac{2}{5}\\ -\frac{1}{5}& \frac{1}{10}& \frac{2}{5}& \frac{7}{10}\end{array}\right]$

$v≔⟨1\,0\,-1\,3⟩$

$v≔\left[\begin{array}{c}1\\ 0\\ \mathrm{-1}\\ 3\end{array}\right]$

$w≔P·v$

$w≔\left[\begin{array}{c}0\\ \frac{1}{2}\\ 1\\ \frac{3}{2}\end{array}\right]$

$\mathrm{Basis}\left(\left[\mathrm{op}\left(S\right)\,w\right]\right)=S$

$\left[\left[\begin{array}{c}1\\ 2\\ 3\\ 4\end{array}\right]\,\left[\begin{array}{c}4\\ 3\\ 2\\ 1\end{array}\right]\right]=\left[\left[\begin{array}{c}1\\ 2\\ 3\\ 4\end{array}\right]\,\left[\begin{array}{c}4\\ 3\\ 2\\ 1\end{array}\right]\right]$

$w·\left(v-w\right)$

$0$

$\mathrm{ProjectionMatrix}\left(\left\{⟨a\,1⟩\right\}\right)$

$\left[\begin{array}{cc}\frac{a^2}{a^2+1}& \frac{a}{a^2+1}\\ \frac{a}{a^2+1}& \frac{1}{a^2+1}\end{array}\right]$

The Student[LinearAlgebra][ProjectionMatrix] command was introduced in Maple 17.

For more information on Maple 17 changes, see Updates in Maple 17.