How to find peak value of an analog signal efficiently after sampling in the digital domain?

Question

I have a bandlimited analog signal for which I want to find the peak value in real time. The signal is sampled and processed digitally at just enough sampling rate. Since the peak of the analog signal may not be sampled exactly, taking the maximum of the digital samples would give a lower than actual value. So how do I find the peak?

Theoretically, the only way I can think of, is to upsample the signal using an ideal sinc interpolation filter to very high sampling rate, so that we obtain a sample close to the actual peak of the analog signal. Then we can just take the maximum of the digital samples. But this seems roundabout and computationally expensive. Is there is a simpler method to do this?

Answer 1

You are unlikely to be able to evaluate correctly an ideal sinc interpolation filter in real-time (or even with a small time lag), since the future samples of your data are unknown. Plus, when you received your last sample $x[n]$, you generally don't know whether you are on a peak. So, it somehow boils down to

• how real-time you want to be;
• how fast you can estimated you are passed a peak.

There are many local possible interpolations. You can also think about fast option, and close to real-time, only require to collect the next sample $x[n+1]$, and if knowing $x[n-1]$ you consider that a peak is located between $x[n-1]$ and $x[n+1]$ (and there values are positive), you can use the standard frequency peak estimators, which only use the above three samples. For instance:

$$n_{\textrm{peak}} = n+\delta$$ with $$\delta = \frac{(x_{n+1}−x_{n-1})}{(4x_{n} − 2x_{n-1} − 2x_{n+1})}$$

A primer is provided by famous people here: Fast, Accurate Frequency Estimators, Eric Jacobsen and Peter Kootsookos, IEEE signal processing magazine, May 2007. The above equation is only one of those this paper provides.

Answer 2

Assuming the signal is bandlimited to less than Nyquist, the peak between two samples can be determined from the history of prior samples. The number of samples needed is based on the actual response time of the filter or channel that the filter has passed through.

I would suggest using polyphase filters as an efficient interpolation approach that can be done in real time without increasing the sampling rate (real-time with a delay as limited by the depth of the filters). There would be N filter banks corresponding to the additional precision desired (as each filter bank results in an interpolated output including and between the original sampling rate outputs). For each of the N outputs that would be given for each original sample, the maximum value can be selected providing a maximum interpolated sample for every clock update. Further, polyphase filter banks can be relatively short, so the overall delay to the solution is also attractively small. (The number of taps needed for each filter in order to obtain an accurate result is equivalent in time span to the significant response of the channel or filter that the waveform has passed through).

Polyphase interpolators are simple to implement and highly efficient in operation. For more details on their implementation and structure see:

How to implement Polyphase filter?

Answer 3

Do you know the exact filter used to bandlimit that signal? Assuming that it is a symmetrically windowed sinc, I am not sure that you could do much better than just upsampling. Perhaps speculatively upsampling only when you find a «large» digital sample?

If the bandlimiting filter is very brickwall-y, it will contain positive and negative sections. I am not sure that a simple fitted Gaussian is going to be a good model of that.

Do you need to identify every single occurence of peaks accurately, at low complexity and low delay?