Nonnegative or positive band-limited interpolation Given samples of an everywhere non-negative or positive-valued continuous-time signal band-limited to half the sampling frequency, is there some practically applicable way to interpolate it so that there is no danger of producing negative values? The approach should be scalable to arbitrary accuracy so basic linear interpolation won't do. By arbitrary accuracy I mean an arbitrary small difference between the band-limited continuous-time signal and the interpolation.

Perfect sinc interpolation would work because negative contributions would get cancelled by positive contributions. Ideal sinc interpolation is not possible in practice so other filters are used. This could give negative values because usually the filter impulse response has negative values. For example cubic Hermite spline piece-wise polynomial interpolation gives this impulse response:

Figure 1. Cubic Hermite spline impulse response has negative values.

Lagrange polynomials used as splines would also have negative values in the impulse response. B-splines have everywhere non-negative impulse responses, but their pass-band frequency responses deviate increasingly from the brick-wall frequency response of sinc interpolation with increasing B-spline order. The B-spline impulse response tends to a Gaussian function as the order is increased. A Gaussian function in time domain is a Gaussian function in frequency domain, so not a perfect brick-wall filter.

The obvious way to fix the problem would be to 1) take the absolute value of the interpolation output or 2) set all negative output values to zero, but I'm wondering if there is some better (inherently non-negative) method, because both 1) and 2) would create discontinuities in the derivative of the signal. This also creates spectral images.

I don't know how bad the problem stated in this question really is, because if the band-limited signal is non-negative, this also rules out some contradictory discrete signals like $[\dots, 0, 0, 0, 1, 0, 0, 0, \dots].$ The corresponding non-negative continuous signal would then have to be sinc, which is partially negative, which is the contradiction.



Shay Maymon has investigated "bandlimited square-roots" in his 2011 thesis Sampling and Quantization for Optimal Reconstruction. The idea there is to find (not easy) another band-limited signal (complex and at half the bandwidth) the squared magnitude of which equals the wanted band-limited signal. Interpolating that other signal and multiplying the result by its complex conjugate is never negative. That reminds me of this question about interpolation of the magnitude of discrete Fourier transform (DFT).
If we require the impulse response of the interpolation filter to be non-negative then we can obviously guarantee that the interpolated signal is also non-negative. However, it can be shown that the frequency response of a low pass filter with non-negative impulse response is always far from ideal. More specifically, it can be shown that if $\epsilon$ is the maximum negative pass band deviation of the squared magnitude $|H(\omega)|^2$ of the filter's transfer function from its value at $\omega=0$

$$|H(\omega)|^2\ge |H(0)|^2-\epsilon,\quad |\omega|&lt;\omega_c \tag{1}$$

where $\omega_c$ is the cut-off frequency, then the following inequality holds:

$$|H(\omega)|^2\ge |H(0)|^2-4\epsilon,\quad |\omega|&lt;2\omega_c\tag{2}$$

This inequality imposes a severe restriction on the minimum stop band attenuation if the ripple size $\epsilon$ is chosen sufficiently small (as should be the case for a good low pass filter).

As an example, if $\epsilon$ is chosen such that the pass band ripple is $1$ dB, the attenuation at $\omega=2\omega_c$ cannot be better than $7.5$ dB.

Of course, if $\epsilon\ge |H(0)|^2/4$ then $(2)$ represents no restriction, but in that case the maximum pass band deviation of the magnitude of the transfer function is greater than usually required for a good low pass filter.

Since the inequality $(2)$ is independent of the filter order, an interpolation filter with arbitrary accuracy cannot be realized by a filter with non-negative impulse response.

The proof of $(2)$ can be found in Attenuation Limits for Filters with Monotonic Step Response by A. Papoulis, IRE Transactions on Circuit Theory ( Volume: 9, Issue: 1, March 1962).