Next  |  Prev  |  Up  |  Top  |  Index  |  JOS Index  |  JOS Pubs  |  JOS Home  |  Search

Projection onto Non-Orthogonal Vectors

Consider this example:

\begin{eqnarray*}
\underline{s}_0 &\isdef & [1,1] \\
\underline{s}_1 &\isdef & [0,1]
\end{eqnarray*}

These point in different directions, but they are not orthogonal. What happens now? The projections are

\begin{eqnarray*}
{\bf P}_{\underline{s}_0}(x) &=& \frac{x_0 + x_1}{2}\underline{s}_0 \\
{\bf P}_{\underline{s}_1}(x) &=& x_1\underline{s}_1.
\end{eqnarray*}

The sum of the projections is

\begin{eqnarray*}
{\bf P}_{\underline{s}_0}(x) + {\bf P}_{\underline{s}_1}(x) &=&
\frac{x_0 + x_1}{2}\underline{s}_0 + x_1\underline{s}_1 \\
&\isdef & \frac{x_0 + x_1}{2}(1,1) + x_1 \cdot (0,1) \\
&=& \left(\frac{x_0 + x_1}{2},
\frac{x_0 + 3x_1}{2}\right) \\
&\neq& x.
\end{eqnarray*}

So, even though the vectors are linearly independent, the sum of projections onto them does not reconstruct the original vector. Since the sum of projections worked in the orthogonal case, and since orthogonality implies linear independence, we might conjecture at this point that the sum of projections onto a set of $ N$ vectors will reconstruct the original vector only when the vector set is orthogonal, and this is true, as we will show.

It turns out that one can apply an orthogonalizing process, called Gram-Schmidt orthogonalization to any $ N$ linearly independent vectors in $ \mathbb{C}^N$ so as to form an orthogonal set which will always work. This will be derived in Section 5.10.4.

Obviously, there must be at least $ N$ vectors in the set. Otherwise, there would be too few degrees of freedom to represent an arbitrary $ x\in\mathbb{C}^N$ . That is, given the $ N$ coordinates $ \{u(n)\}_{n=0}^{N-1}$ of $ x$ (which are scale factors relative to the coordinate vectors $ \underline{e}_n$ in $ \mathbb{C}^N$ ), we have to find at least $ N$ coefficients of projection (which we may think of as coordinates relative to new coordinate vectors $ \underline{s}_k$ ). If we compute only $ M<N$ coefficients, then we would be mapping a set of $ N$ complex numbers to $ M<N$ numbers. Such a mapping cannot be invertible in general. It also turns out $ N$ linearly independent vectors is always sufficient. The next section will summarize the general results along these lines.


Next  |  Prev  |  Up  |  Top  |  Index  |  JOS Index  |  JOS Pubs  |  JOS Home  |  Search

[How to cite this work]  [Order a printed hardcopy]  [Comment on this page via email]

``Mathematics of the Discrete Fourier Transform (DFT), with Audio Applications --- Second Edition'', by Julius O. Smith III, W3K Publishing, 2007, ISBN 978-0-9745607-4-8
Copyright © 2024-04-02 by Julius O. Smith III
Center for Computer Research in Music and Acoustics (CCRMA),   Stanford University
CCRMA