|
|
| Line 1: |
Line 1: |
| {{Unreferenced|date=August 2007}}In [[signal processing]], '''reconstruction''' usually means the determination of an original continuous signal from a sequence of equally spaced samples.
| | Hi there, I am Alyson Boon even though it is not the name on my beginning certificate. Credit authorising is exactly where my main income comes from. I've always loved residing in Mississippi. Doing ballet is some thing she would never give up.<br><br>my blog: [http://wellness-guide-videos.com/users/PAviles psychic phone] |
| | |
| This article takes a generalized abstract mathematical approach to signal sampling and reconstruction. For a more practical approach based on band-limited signals, see [[Whittaker–Shannon interpolation formula]].
| |
| | |
| == General principle ==
| |
| | |
| Let ''F'' be any sampling method, i.e. a linear map from the [[Hilbert space]] of square-integrable functions <math>L^2</math> to [[complex number|complex]] space <math>\Bbb C^n</math>.
| |
| | |
| In our example, the vector space of sampled signals <math>\Bbb C^n</math> is ''n''-dimensional complex space. Any proposed inverse ''R'' of ''F'' (''reconstruction formula'', in the lingo) would have to map <math>\Bbb C^n</math> to some subset of <math>L^2</math>. We could choose this subset arbitrarily, but if we're going to want a reconstruction formula ''R'' that is also a linear map, then we have to choose an ''n''-dimensional linear subspace of <math>L^2</math>.
| |
| | |
| This fact that the dimensions have to agree is related to the [[Nyquist–Shannon sampling theorem]].
| |
| | |
| The elementary linear algebra approach works here. Let <math>d_k:=(0,...,0,1,0,...,0)</math> (all entries zero, except for the ''k''th entry, which is a one) or some other basis of <math>\Bbb C^n</math>. To define an inverse for ''F'', simply choose, for each ''k'', an <math>e_k \in L^2</math> so that <math>F(e_k)=d_k</math>. This uniquely defines the (pseudo-)inverse of ''F''.
| |
| | |
| Of course, one can choose some reconstruction formula first, then either compute some sampling algorithm from the reconstruction formula, or analyze the behavior of a given sampling algorithm with respect to the given formula.
| |
| | |
| == Popular reconstruction formulae ==
| |
| Perhaps the most widely used reconstruction formula is as follows. Let <math>\{ e_k \}</math> is a basis of <math>L^2</math> in the Hilbert space sense; for instance, one could use the eikonal
| |
| | |
| :<math>e_k(t):=e^{2\pi i k t}\,</math>, | |
| | |
| although other choices are certainly possible. Note that here the index ''k'' can be any integer, even negative.
| |
| | |
| Then we can define a linear map ''R'' by
| |
| | |
| :<math>R(d_k)=e_k\,</math>
| |
| | |
| for each <math>k=\lfloor -n/2 \rfloor,...,\lfloor (n-1)/2 \rfloor</math>, where <math>(d_k)</math> is the basis of <math>\Bbb C^n</math> given by
| |
| | |
| :<math>d_k(j)=e^{2 \pi i j k \over n}</math>
| |
| | |
| (This is the usual discrete Fourier basis.)
| |
| | |
| The choice of range <math>k=\lfloor -n/2 \rfloor,...,\lfloor (n-1)/2 \rfloor</math> is somewhat arbitrary, although it satisfies the dimensionality requirement and reflects the usual notion that the most important information is contained in the low frequencies. In some cases, this is incorrect, so a different reconstruction formula needs to be chosen.
| |
| | |
| A similar approach can be obtained by using [[wavelet]]s instead of Hilbert bases. For many applications, the best approach is still not clear today.
| |
| | |
| == See also ==
| |
| | |
| * [[Nyquist–Shannon sampling theorem]]
| |
| * [[Whittaker–Shannon interpolation formula]]
| |
| * [[Aliasing]]
| |
| | |
| [[Category:Signal processing|Reconstruction]]
| |
Hi there, I am Alyson Boon even though it is not the name on my beginning certificate. Credit authorising is exactly where my main income comes from. I've always loved residing in Mississippi. Doing ballet is some thing she would never give up.
my blog: psychic phone