Site Feedback

Perfect Reconstruction

By Stephen Kruzick, Roy Ha, Justin Romberg


If certain additional assumptions about the original signal and sampling rate hold, then the original signal can be recovered exactly from its samples using a particularly important type of filter. More specifically, it will be shown that if a bandlimited signal is sampled at a rate greater than twice its bandlimit, the Whittaker-Shannon reconstruction formula perfectly reconstructs the original signal. This formula makes use of the ideal lowpass filter, which is related to the sinc function. This is extremely useful, as sampled versions of continuous time signals can be filtered using discrete time signal processing, often in a computer. The results may then be reconstructed to produce the same continuous time output as some desired continuous time system.

Perfect Reconstruction

In order to understand the conditions for perfect reconstruction and the filter it employs, consider the following. As a beginning, a sufficient condition under which perfect reconstruction is possible will be discussed. Subsequently, the filter and process used for perfect reconstruction will be detailed.

Recall that the sampled version xsxs of a continuous time signal xx with sampling period TsTs has a spectrum given by

X s ( ω ) = 1 T s k = - X ω - 2 π k T s . X s ( ω ) = 1 T s k = - X ω - 2 π k T s .

As before, note that if xx is bandlimited to (-π/Ts,π/Ts)(-π/Ts,π/Ts), meaning that XX is only nonzero on (-π/Ts,π/Ts)(-π/Ts,π/Ts), then each period of XsXs has the same form as XX. Thus, we can identify the original spectrum XX from the spectrum of the samples XsXs and, by extension, the original signal xx from its samples xsxs at rate TsTs if xx is bandlimited to (-π/Ts,π/Ts)(-π/Ts,π/Ts).

If a signal xx is bandlimited to (-B,B)(-B,B), then it is also bandlimited to (-π/Ts,π/Ts)(-π/Ts,π/Ts) provided that Ts<π/BTs<π/B. Thus, if we ensure that xx is sampled to xsxs with sufficiently high sampling angular frequency ωs=2π/Ts>2Bωs=2π/Ts>2B and have a way of identifying the unique (-π/Ts,π/Ts)(-π/Ts,π/Ts) bandlimited signal corresponding to a discrete time signal at sampling period TsTs, then xsxs can be used to reconstruct x˜=xx˜=x exactly. The frequency 2B2B is known as the angular Nyquist rate. Therefore, the condition that the sampling rate ωs=2π/Ts>2Bωs=2π/Ts>2B be greater than the Nyquist rate is a sufficient condition for perfect reconstruction to be possible.

The correct filter must also be known in order to perform perfect reconstruction. The ideal lowpass filter defined by G(ω)=Tsuω+π/Ts-uω-π/TsG(ω)=Tsuω+π/Ts-uω-π/Ts, which is shown in Figure 1, removes all signal content not in the frequency range (-π/Ts,π/Ts)(-π/Ts,π/Ts). Therefore, application of this filter to the impulse train n=-xs(n)δ(t-nTs)n=-xs(n)δ(t-nTs) results in an output bandlimited to (-π/Ts,π/Ts)(-π/Ts,π/Ts).

We now only need to confirm that the impulse response gg of the filter GG satisfies our sufficient condition to be a reconstruction filter. The inverse Fourier transform of G(ω)G(ω) is

g ( t ) = sinc ( t / T s ) = 1 t = 0 sin ( π t / T s ) π t / T s t 0 , g ( t ) = sinc ( t / T s ) = 1 t = 0 sin ( π t / T s ) π t / T s t 0 ,

which is shown in Figure 1. Hence,

g ( n T s ) = sinc ( n ) = 1 n = 0 sin ( π n ) π n n 0 = 1 n = 0 0 n 0 = δ ( n ) . g ( n T s ) = sinc ( n ) = 1 n = 0 sin ( π n ) π n n 0 = 1 n = 0 0 n 0 = δ ( n ) .

Therefore, the ideal lowpass filter GG is a valid reconstruction filter. Since it is a valid reconstruction filter and always produces an output that is bandlimited to (-π/Ts,π/Ts)(-π/Ts,π/Ts), this filter always produces the unique (-π/Ts,π/Ts)(-π/Ts,π/Ts) bandlimited signal that samples to a given discrete time sequence at sampling period TsTs when the impulse train n=-xs(n)δ(t-nTs)n=-xs(n)δ(t-nTs) is input.

Therefore, we can always reconstruct any (-π/Ts,π/Ts)(-π/Ts,π/Ts) bandlimited signal from its samples at sampling period TsTs by the formula

x ( t ) = n = - x s ( n ) sinc ( t / T s - n ) . x ( t ) = n = - x s ( n ) sinc ( t / T s - n ) .

This perfect reconstruction formula is known as the Whittaker-Shannon interpolation formula and is sometimes also called the cardinal series. In fact, the sinc function is the infinite order cardinal basis spline ηη. Consequently, the set {sinc(t/Ts-n)|nZ}{sinc(t/Ts-n)|nZ} forms a basis for the vector space of (-π/Ts,π/Ts)(-π/Ts,π/Ts) bandlimited signals where the signal samples provide the corresponding coefficients. It is a simple exercise to show that this basis is, in fact, an orthogonal basis.

Figure 1: The above plots show the ideal lowpass filter and its inverse Fourier transform, the sinc function.
Figure 2: The plots show an example discrete time signal and its Whittaker-Shannon sinc reconstruction.

Perfect Reconstruction Summary

This module has shown that bandlimited continuous time signals can be reconstructed exactly from their samples provided that the sampling rate exceeds the Nyquist rate, which is twice the bandlimit. The Whittaker-Shannon reconstruction formula computes this perfect reconstruction using an ideal lowpass filter, with the resulting signal being a sum of shifted sinc functions that are scaled by the sample values. Sampling below the Nyquist rate can lead to aliasing which makes the original signal irrecoverable as is described in the subsequent module. The ability to perfectly reconstruct bandlimited signals has important practical implications for the processing of continuous time signals using the tools of discrete time signal processing.