Imagine we have band limited signal with bandwidth of $B$, so the required Nyquist rate would be $f_{nyq}>2B$ that is oversampled with rates $f_s$ where $f_s = M*f_{nyq}$ and $M$ is random and $M>=1$. By that, I mean the sampling rate is at least at Nyquist rate but most of the time its way more than Nyquist and it is random. I have heard and read about iterative My questions are:

  1. Is there a method to construct the signal as if it was taken in Nyquist more, without degradation of signal quality? After all, we have more samples than needed, so I do not expect any degradation, right?

  2. Is there are some loss, how much is it, what factors it depends on?

  3. Is there a toolbox or something that I can feed my data in and have some experiments of my own? ,

    4.Is the technique used in MATLAB lossy?

I appreciate any kind of contribution here :)

enter image description here


2 Answers 2


If you assume the signal was strictly bandlimited to below some Nyquist frequency, then it can be decomposed into some number (N) of DFT basis vectors over the sampling aperture (although that decomposition will include rectangular windowing artifacts if the signal wasn't integer periodic in the aperture width). If you have enough sample points (M >= N) of that signal, then this becomes a problem of fitting M equations to N unknowns. IIRC, the farther the sample points are from being equally spaced over the aperture, the more sensitive any computed solution might be to noise and numerical issues. Once deconstructed into DFT basis vectors, any other sample points of that strictly bandlimited signal can be interpolated using a summation of the resultant complex exponential coefficients.

  • $\begingroup$ So the solution is not a perfect reconstruction? $\endgroup$
    – MimSaad
    Aug 11, 2017 at 20:25
  • $\begingroup$ The reconstruction will be as perfect as the band-limiting and quality of the solution or estimation of the DFT basis vectors allows. $\endgroup$
    – hotpaw2
    Aug 11, 2017 at 20:29
  • $\begingroup$ so the more bandlimited signal is, the less unkowns and the more samples we have the more equations are added to our overdetermined system of equations, but we always have leakage so even bandlimited signal might occupy the whole spectrum $\endgroup$
    – MimSaad
    Aug 11, 2017 at 20:38
  • $\begingroup$ Windowing artifacts are always an issue with signals that are not integer periodic within a finite aperture, even with equally spaced sampling above the Nyquist rate. $\endgroup$
    – hotpaw2
    Aug 11, 2017 at 21:51
  • $\begingroup$ As far as I understood your solution is Fourier based interpolation, true? If we choose other kernels, does this windowing effect cause problem? $\endgroup$
    – MimSaad
    Aug 12, 2017 at 10:25

For the moment restrict the discussion to 1D signals. Sampling operation can be performed as uniformly or non-uniformly. The uniform sampling is the most obvious, simplest and the preferred form, unless otherwise stated.

Nonuniform sampling can be performed in a number of ways, periodically or aperiodically, sampling instant given by a formula or being completely random.

Depending on the nonuniform sampling form, the reconstruction criterions may change. For example when you have a formula based nonuniform sampling instant $t_n$, you can (if possible) find a kernel $\phi_n(t)$ for the reconstruction integral that provides exact (perfect) recovery of a continuous time signal $x(t)$ from its nonuniform samples $x_n$. In this context the $sinc()$ kernel is a special case of the generalized sampling theorem WKS (Whittaker-Kotelnikov-Shannon) with sampling instants given by $t = n T_s$. This was historically called as the Nyquist sampling, who most successfully emplyed it into the earliest digital pulse communication systems.

The primary requirement that Nyquist (uniform) sampling theorem would provide exact reconstruction of a bandlimited baseband signal $x(t)$ from its uniform samples is that the sampling rate satisfy the Nyquist criterion; $F_s \geq 2 F_c$ where $F_c$ is the signal's bandwidth.

Now for the non-uniform sampling strategies, the most common restriction for perfect signal recovery is that on the average the Nyquist rate associated with the bandlimited signal is maintained. So when there is a deficiency of samples at some short interval $I_1$ then there must be an excess of samples at another closeby short interval $I_2$ so that the average Nyquist rate preserved.

The above condition can be stated as the local Nyquist rate criterion.

Note that unlike uniform sampling case, for nonuniform sampling strategies and exlcusively true for the random sampling strategy, the sampling times $t_n$ in addition to the samples $x_n$ should be known in order to reconstruct the signal perfectly.

You can look for Whittaker's, Yen's, Shannon's papers for classical view on the nonuniform sampling of signals.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.