|X||is the input signal|
|M||mode ('q' performs quadratic interpolation, 'v' finds valleys instead of peaks)|
|W||is the width tolerance; a peak will be eliminated if there is a higher peak within +-W samples|
|K||are the peak locations in X (fractional if M='q')|
|V||are the peak amplitudes: if M='q' the amplitudes will be interpolated whereas if M~='q' then V=X(K).|
Outputs are column vectors regardless of whether X is row or column. If there is a plateau rather than a sharp peak, the routine will place the peak in the centre of the plateau. When the W input argument is specified, the routine will eliminate the lower of any pair of peaks whose separation is <=W; if the peaks have exactly the same height, the second one will be eliminated. All peak locations satisfy 1<K<length(X).
If no output arguments are specified, the results will be plotted.