Identification of linear systems using polynomial kernels in the frequency domain

In prior work we presented an identification algorithm using polynomials in the time domain. In this article, we extend this algorithm to include polynomials in the frequency domain. A polynomial is used to represent the imaginary part of the Fourier transform of the impulse response. The Hilbert transform relationship is used to compute the real part of the frequency response and hence the complete process model. The polynomial parameters are computed based on the computationally efficient linear least square method. The order of the polynomial is estimated based on residue decrement. Simulated and experimental results show the effectiveness of this method, particularly for short input/output data sequence with high signal to noise ratio. The frequency domain polynomial model complements the time domain methods since it can provide a good estimate of the time to steady state for time domain FIR (finite impulse response) models. Confidence limits in time or frequency domain can be computed using this approach. Noise rejection properties of the algorithm are illustrated using data from both simulated and real processes.