Bandlimited Interpolation of Time-Limited Signals

The previous result can be applied to bandlimited interpolation of arbitrary time-limited signals $x\in{\bf C}^N$ (i.e., not just periodic signals) by (1) replacing the rectangular window $\hbox{\sc Chop}_N()$ with a smoother spectral window $H(\omega_k)$, and (2) using extra zero-padding in the time domain to convert the cyclic convolution between $\hbox{\sc Stretch}_L(x)$ and $h$ into an acyclic convolution between them (recall §7.2.4). The smoother spectral window $H$ can be thought of as the frequency response (sampled) of the FIR^7.15 filter $h$ used as the bandlimited interpolation kernel in the time domain. The number of extra zeros appended to $x$ in the time domain is simply length of $h$ minus 1, and the number of zeros appended to $h$ is the length of $\hbox{\sc Stretch}_L(x)$ minus 1. If $N_x$ denotes the nonzero length of $x$, then the nonzero length of $\hbox{\sc Stretch}_L(x)$ is $L(N_x-1)+1$. Thus, we require the DFT length to be $N\geq L(N_x-1)+N_h$, where $N_h$ is the filter length. In operator notation, we can express bandlimited sampling-rate up-conversion by the factor $L$ for time-limited signals $x\in{\bf C}^{N_x}$ by

$$\zbox {\hbox{\sc Interp}_L(x) \approx \hbox{\sc IDFT}(H\cdot(\hbox{\sc DFT}(\hbox{\sc ZeroPad}_{N}(\hbox{\sc Stretch}_L(x)))).} \protect$$ (7.10)

The approximation symbol `$\approx$' approaches equality as the spectral window $H$ approaches $\hbox{\sc Chop}_{N_x}([1,\dots,1])$ (the frequency response of the ideal lowpass filter passing only the original spectrum $X$), while at the same time allowing no time aliasing (convolution is still acyclic in the time domain).

Equation (7.10) can provide the basis for a high-quality sampling-rate conversion algorithm. Arbitrarily long signals can be accommodated by breaking them into segments of length $N_x$, applying the above algorithm to each block, and summing the up-sampled blocks using overlap-add. That is, the lowpass filter $h$ ``rings'' into the next block and possibly beyond (or even into both adjacent time blocks when $h$ is not causal), and this ringing must be summed into all affected adjacent blocks. Finally, the filter $H$ can ``window away'' more than the top $L-1$ copies of $X$ in $Y$, thereby preparing the time-domain signal for downsampling, say by $M\in{\bf Z}$:

$${\small \zbox {\hbox{\sc Interp}_{L/M}(x) \approx \hbox{\sc Downs... ...T}(H\cdot(\hbox{\sc DFT}(\hbox{\sc ZeroPad}_{N}(\hbox{\sc Stretch}_L(x)))))} }$$

where now the lowpass filter frequency response $H$ must be close to zero for all $\vert\omega_k\vert\geq\pi/\max(L,M)$. While such a sampling-rate conversion algorithm can be made more efficient by using an FFT in place of the DFT (see Appendix A), it is not necessarily the most efficient algorithm possible. This is because (1) $M-1$ out of $M$ output samples from the IDFT need not be computed at all, and (2) $\hbox{\sc Stretch}_L(x)$ has many zeros in it which do not need explicit computation. For an introduction to time-domain sampling-rate conversion (bandlimited interpolation) algorithms which take advantage of points (1) and (2) in this paragraph, see, e.g., Appendix D and [70].