This paper presents an example of practical use of Prony's method for monitoring of power waveform fundamental harmonic fluctuations, which is required for the analysis of window synchronizations in frequency analyses in power monitoring systems. The example presented offers verification of the authors' theoretical considerations published earlier in articles about Prony's method and its opportunities for practical use for real life signals. The investigations shown are based on the least squares Prony's method, which, in connection with digital filtrations, enables estimations of fundamental frequency at the rate of even tens of times per one fundamental harmonic period.

The article presents a new modification of the the least squares Prony method. The so-called variable-frequency Prony method can be a useful tool for estimating parameters of sinusoidal components, which, in the analyzed signal, are characterized by time-dependent frequencies. The authors propose use of the presented method for testing the quality of electric energy. It allows observation of phenomena which, when using traditional methods, are averaged in the analysis window. The proposed modification of least squares Prony method is based on introduction and specific selection of a frequency matrix. This matrix represents frequencies of estimated components and their variability in time.

This paper presents a universal approximation of the unit circle by a polygon that can be used in signal processing algorithms. Optimal choice of the values of three parameters of this approximation allows one to obtain a high accuracy of approximation. The approximation described in the paper has a universal character and can be used in many signal processing algorithms, such as DFT, that use the mathematical form of the unit circle. One of the applications of the described approximation is the DFT linear interpolation method (LIDFT). Applying the results of the presented paper to improve the LIDFT method allows one to significantly decrease the errors in estimating the amplitudes and frequencies of multifrequency signal components. The paper presents the derived formulas, an analysis of the approximation accuracy and the region of best values for the approximation parameters.

This overview paper presents and compares different methods traditionally used for estimating damped sinusoid parameters. Firstly, direct nonlinear least squares fitting the signal model in the time and frequency domains are described. Next, possible applications of the Hilbert transform for signal demodulation are presented. Then, a wide range of autoregressive modelling methods, valid for damped sinusoids, are discussed, in which frequency and damping are estimated from calculated signal linear self-prediction coefficients. These methods aim at solving, directly or using least squares, a matrix linear equation in which signal or its autocorrelation function samples are used. The Prony, Steiglitz-McBride, Kumaresan-Tufts, Total Least Squares, Matrix Pencil, Yule-Walker and Pisarenko methods are taken into account. Finally, the interpolated discrete Fourier transform is presented with examples of Bertocco, Yoshida, and Agrež algorithms. The Matlab codes of all the discussed methods are given. The second part of the paper presents simulation results, compared with the Cramér-Rao lower bound and commented. All tested methods are compared with respect to their accuracy (systematic errors), noise robustness, required signal length, and computational complexity.