Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 712509, 5 pages

http://dx.doi.org/10.1155/2014/712509

Research Article

Parameters Estimation of Interharmonic Based on State Space Model

Yuan Shiji, Pei Bin, Liu Zhihua, Huang Wenjing, and Sun Mingfeng

Luoyang Electronic Equipment Test Center, Jiyuan, Henan 459000, China

Received 21 May 2013; Revised 6 November 2013; Accepted 7 November 2013; Published 13 February 2014

Academic Editor: Grigori M. Sisoev

Copyright © 2014 Yuan Shiji et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The special Hankel matrix is structured from interharmonic sampling, which is described by state space model. A method of parameters estimation based on state space model is proposed, which can achieve interharmonics frequency, amplitude and, phase of the joint estimation. The simulation results show that the method can effectively restrain white Gaussian noise, with superior performance.

1. Introduction

Traditional power system harmonic analysis method based on Fourier transform (FFT) has the frequency spectrum revelation and the stockade effect, so there are limitations in the detection of harmonic [1, 2]. Some methods have been applied to estimate the damped exponential model parameters based on modern signal processing theory, such as Pisarenko harmonic decomposition, AR model, and multisignal classification (MUSIC) [3–5]. They have high estimation accuracy, robustness, and other advantages, but most of these methods that need to search in the frequency domain spectra did not realize the frequency, amplitude, and other parameters of the joint estimation.

The parameter estimation method based on state space model (SSM) can inhibit the channel between any zero-mean noise, and it also improves the resolution of parameter estimation [6–9]. The special Hankel matrix is structured from interharmonic sampling; then a method of parameter estimation based on state space model has achieved the joint estimation. The simulation results show that the method can effectively restrain Gaussian color noise, with superior performance.

2. Interharmonic Signal Model

The interharmonics are defined as follows in IEC-61000-2-2: there is not an integer multiple signal between the frequency and the fundamental frequency in the voltage and current signal harmonic. Power system harmonics and inter-harmonic detection signal can be expressed as [3].

Consider the following:

where is sampling signal; is sampling points; is inter-harmonic number; is the th inter-harmonic component of amplitude; is the th inter-harmonic frequency; is the initial angle; is a zero-mean white Gaussian noise.

3. Parameters Estimation Algorithm Based on SSM

Equation (1) is transformed for complex frequency signal, because the state space model is transformed from the ARMA model [10]. With the Hilbert transform (1) can be expressed as follows:

Equation (2) can be described using SSM [11, 12]:

where stands for the output of the fourth order mixed cumulants of , stands for the state variables, stands for the pulse function, stands for the system state matrix, and , , and stand for constant matrix. The initial conditions of . can be obtained form the state space theory [11, 12]:

Data matrix is constructed as follows:

where the stands for the correlation window length stands for the data length. Another expression of can be obtained with substitution of (4) into (5):

The and are defined as follows:

The and have the following properties:

where the subscripts 1 and 2 represent deletion of the first row of (column) and the last row (column) of corresponding matrix. Singular value decomposition of matrix is

where the subscripts and stand for the signal and noise subspace and ‘‘’’ stands for conjugate transpose. The left standard matrix is composed by the signal matrix and noise matrix , and the right standard matrix is composed by the signal matrix and noise matrix . is the diagonal matrix of singular values corresponding to the signal, and is the diagonal matrix of singular values ​ ​corresponding to the noise. The matrix should be removed from the noise caused by the singular value and its corresponding eigenvector in order to increase the accuracy of estimated state matrix ; that is, . The -order approximate estimation of and can be obtained from (9) with as the model order. The state transfer matrix approximation can be obtained according to (8).

Consider the following:

The ‘‘’’ represents the conjugate transpose. Let eigenvalues ​ ​of matrix be :

Define a new matrix as follows:

Then , . is the state space corresponding to the transfer function of the pole from the spectral estimation theory.

Consider the following:

where . There is , and is eigenvectors of . The matrix is composed of ,

that is, ,

The substitution into (4) yields .

Make the following:

; . The can be described in (2):

The algorithm steps of estimate model pole based on SSM method will be described as follows:

(i)constructing data matrix according to (5);(ii)singular value decomposition with data matrix ;(iii)calculating the model order according to the literature [11] and -order approximation;(iv)estimation of being obtained from (9);(v)calculating approximation of according to (10);(vi)calculating the eigenvalues of the matrix that is the damped exponential model of pole parameters;(vii)calculating inter-harmonic component of amplitude according to (17).

4. Simulation

The experiment sampling signal that is inter-harmonic is described in (1), and the number of inter-harmonics is ; sampling points are . The parameters are set in Table 1. Experiments with additive white Gaussian noise with variance. The SNR is defined as the signal and the noise ratio of average power:

Data correlation length and .

Table 1: Interharmonic parameters.

1
2
3

0.1
0.2
0.3

0.2
0.4
0.8

π/3 (1.0472)
π/4 (0.7854)
π/5 (0.6283)

The results of this experiment are shown in Table 2 for 200 Monte-Carlo simulations at different SNR on damped exponential model parameter estimates. The method can be used to more accurately achieve damped frequency, damped factor of the joint estimation, and the method is less sensitive to the white Gaussian noise. The accuracy of damped frequency is higher than damped factor, because other parameters are estimated by the damped frequency that is firstly estimated. From this table it can also be seen that the damped frequency estimation is superior as the damped frequency is larger, because the signal model is closer to smooth signal with slowly damping.

Table 2: White Gaussian noise, 200 simulations, and interharmonics parameter estimation bias.

SNR/dB

(π/3)
(π/4)
(π/5)

8

10

15

20

25
—
—
—

The comparison of logarithm mean square error (MSE) of estimation frequency at different SNR between these methods with the traditional SSM algorithm is shown in Figure 1. The traditional SSM parameter estimation method is more sensitive to the Gaussian colored noise, while the SSM parameter estimation method based on FOMC can inhibit the Gaussian colored noise from Figure 1.

(a)

(b)

(c)

(d)

Figure 1: Logarithm mean square error at different method and SNR.

Figure 2 shows that the comparison of logarithm mean square error of estimation frequency at different SNR between the methods with FOC-MUSIC algorithm is proposed in the literature [3]. The performance of this method is superior to FOC-MUSIC algorithm proposed in the literature [3] at low SNR with the same superior performance between the two estimates at high SNR from Figure 2. At the same time, this method can achieve damped frequency, damped factor of the joint estimation with less calculation.

(a)

(b)

(c)

(d)

Figure 2: MSE of estimation frequency at different method and SNR.

5. Conclusion

Traditional inter-harmonic detection method needs to search for peaks in the frequency domain that did not realize the frequency, amplitude, and other parameters of the joint estimation. Parameter estimation method based on SSM achieves damped frequency, damped factor of the joint estimation, which avoids the search in the frequency domain spectra with superior performance.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

W. C. Zhao, W. M. Ma, and A. Hu, “FFT algorithm with high accurary for harmonic analysis in the electric machine,” in Proceedings of the CSEE, vol. 21, pp. 83–87, 2001.

C. J. Qi and X. H. Wang, “Inter harmonics estimation based on interpolation FFT algorithm,” Transactions of China Electro Technical Society, vol. 18, no. 1, pp. 92–95, 2003.

J. W. Zhang, N. C. Zhou, F. Yang, and C. X. Cuo, “Research on the inte-harmonics measure using fourth-order cumulant based multiple signal classification method,” Relay, vol. 36, no. 7, pp. 19–23, 2008.

S. Rouquette and M. Najim, “Estimation of frequencies and damping factors by two-dimensional ESPRIT type methods,” IEEE Transactions on Signal Processing, vol. 49, no. 1, pp. 237–245, 2001.

P. Chevalier, A. Ferréol, and L. Albera, “High-resolution direction finding from higher order statistics: the 2q-MUSIC algorithm,” IEEE Transactions on Signal Processing, vol. 54, no. 8, pp. 2986–2997, 2006.

Z.-J. Liu, X.-T. Kang, L.-L. Zhang, and Y.-W. Shi, “Tls-Hankel method based on the state space model of harmonic retrieval,” Journal of Jilin University (Engineering and Technology Edition), vol. 38, no. 1, pp. 224–228, 2008.

J. E. Piou, “System realization using 2-D output measurements,” in Proceedings of the IEEE American Control Conference (AAC '04), pp. 2840–2845, July 2004.

K. Naishadham and J. E. Piou, “A robust state space model for the characterization of extended returns in radar target signatures,” IEEE Transactions on Antennas and Propagation, vol. 56, no. 6, pp. 1742–1751, 2008.

S. J. Yuan, Research on Feature Extraction of Midcourse Ballistic Target Based on Wideband Information, NUDT.HU Nan, Changsha, China, 2009.

K. Naishadham and J. E. Piou, “A robust state space model for the characterization of extended returns in radar target signatures,” IEEE Transactions on Antennas and Propagation, vol. 56, no. 6, pp. 1742–1751, 2008.

D. Z. Zhang, Linear System Theory, Tsinghua University Press, 2002.

S.-J. Yuan, X.-Z. Gao, G.-H. Jin, and H.-Q. Wang, “Parameter estimation of radar target 1-D scattering centers based on SSM,” Systems Engineering and Electronics, vol. 33, no. 4, pp. 779–782, 2011.