Summary
The Z-transform of a complex time signal (or the analytic signal of a real signal) is equal to the Z-transform of a prediction error divided by the Z-transform of the prediction error operator. This inverse is decomposed into a sum of partial fractions, which are used to obtain impulse response operators formed by non-causal filters which complex-conjugate symmetric coefficients. The time-components are obtained by convolving the filters with the original signal, and the peak frequencies, corresponding to the poles of the prediction error operator, are used for mapping the time-components into frequency components. For non-stationary signals, this decomposition is done in sliding time windows, and the signal component values, in the middle of each window, are attributed to the peak value of its frequency response which corresponds to the pole of this partial fraction component. The result is an exact, but non-unique, time-frequency representation of the input signal. A sparse signal decomposition can be obtained by summing along the frequency axis in patches with similar characteristics in the time-frequency domain. The peak amplitude frequency of each new time component is obtained by computing a scalar prediction error operator in sliding time windows, resulting in a sparse time-frequency representation. In both cases, the result is a time-frequency matrix where an estimate of the frequency content of the input signal can be obtained by summation over the time variable. The performance of the new method is demonstrated with excellent results on a synthetic time signal, the LIGO gravitational wave signal, and on seismic field data.