Sign in

username:

password:



Not a member?

Search Online Books



Search tips

Free Online Books

Sponsor

NEW! TMS320C6474: 3x the performance. 1/3 the cost. Three 1 GHz cores on 1 chip.

Chapters

Chapter Contents:

Search Mathematics of the DFT

  

Book Index | Global Index


Would you like to be notified by email when Julius Orion Smith III publishes a new entry into his blog?

  

Projection onto Non-Orthogonal Vectors

Consider this example:

\begin{eqnarray*}
\sv_0 &\isdef & [1,1] \\
\sv_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}_{\sv_0}(x) &=& \frac{x_0 + x_1}{2}\sv_0 \\
{\bf P}_{\sv_1}(x) &=& x_1\sv_1.
\end{eqnarray*}

The sum of the projections is

\begin{eqnarray*}
{\bf P}_{\sv_0}(x) + {\bf P}_{\sv_1}(x) &=&
\frac{x_0 + x_1}...
...\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 $ {\bf 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{\bf 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 $ {\bf C}^N$), we have to find at least $ N$ coefficients of projection (which we may think of as coordinates relative to new coordinate vectors $ \sv_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.


Order a Hardcopy of Mathematics of the DFT

Previous: Projection onto Linearly Dependent Vectors
Next: General Conditions

written by Julius Orion Smith III
Julius Smith's background is in electrical engineering (BS Rice 1975, PhD Stanford 1983). He is presently Professor of Music and Associate Professor (by courtesy) of Electrical Engineering at Stanford's Center for Computer Research in Music and Acoustics (CCRMA), teaching courses and pursuing research related to signal processing applied to music and audio systems. See http://ccrma.stanford.edu/~jos/ for details.


Comments


No comments yet for this page


Add a Comment
You need to login before you can post a comment (best way to prevent spam). ( Not a member? )