Sinc Interpolation for Signal Reconstruction

This Demonstration illustrates the use of the sinc interpolation formula to reconstruct a continuous signal from some of its samples. The formula provides exact reconstructions for signals that are bandlimited and whose samples were obtained using the required Nyquist sampling frequency, to eliminate aliasing in the reconstruction of the signal.
You can apply the interpolation formula to a number of continuous signals. Increasing the sampling frequency gives a more accurate reconstruction of the continuous function.
The original signal is shown as a blue solid line and the sample locations are shown by red circles. The reconstructed signal is shown using the dotted magenta line and is superimposed on the original signal to make it easier to see the effect of increasing the sampling frequency on the reconstruction of the original signal from its samples.
Because the number of samples is limited, there will be some error in the reconstruction, even at a higher sampling frequency than the Nyquist frequency. The absolute (point by point) error shown by the original signal and the reconstructed signal is plotted to show how the error decreases as the sampling rate increases.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


The function is defined by for , with . The sinc interpolation formula is defined as , where is the sampling period used to determine from the original signal, and is the reconstructed signal. The above formula represents a linear convolution between the sequence and scaled and shifted samples of the function. In this Demonstration, a limited number of samples are generated, and the above sum is carried out for samples, labeled from to . Due to the shifting of the function by integer multiples of , this results in having the exact value of a sample located at a multiple of . This can be seen by observing that the absolute error is always zero at times which are integer multiples of , in other words at the sample locations. In this implementation, the function is sampled at a much higher rate than the sampling frequency used for the original function, in order to produce a smoother plotted result.
A. V. Oppenheim and R. W. Schafer, Digital Signal Processing, Englewood Cliffs, NJ: Prentice Hall, 1975.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2014 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+