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.
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 of a continuous time signal with sampling period has a spectrum given by
As before, note that if is bandlimited to , meaning that is only nonzero on , then each period of has the same form as . Thus, we can identify the original spectrum from the spectrum of the samples and, by extension, the original signal from its samples at rate if is bandlimited to .
If a signal is bandlimited to , then it is also bandlimited to provided that . Thus, if we ensure that is sampled to with sufficiently high sampling angular frequency and have a way of identifying the unique bandlimited signal corresponding to a discrete time signal at sampling period , then can be used to reconstruct exactly. The frequency is known as the angular Nyquist rate. Therefore, the condition that the sampling rate 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 , which is shown in Figure 1, removes all signal content not in the frequency range . Therefore, application of this filter to the impulse train results in an output bandlimited to .
We now only need to confirm that the impulse response of the filter satisfies our sufficient condition to be a reconstruction filter. The inverse Fourier transform of is
which is shown in Figure 1. Hence,
Therefore, the ideal lowpass filter is a valid reconstruction filter. Since it is a valid reconstruction filter and always produces an output that is bandlimited to , this filter always produces the unique bandlimited signal that samples to a given discrete time sequence at sampling period when the impulse train is input.
Therefore, we can always reconstruct any bandlimited signal from its samples at sampling period by the formula
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 forms a basis for the vector space of 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.
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.