Next Article in Journal
Physical Constraints on Global Social-Ecological Energy System
Previous Article in Journal
Fog Computing Approach for Shared Mobility in Smart Cities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Dynamic Modeling of KF Algorithm for Detecting Distorted AC Signal

1
School of Economics and Management, Nanchang University, Nanchang 330031, China
2
I.H. Asper School of Business, University of Manitoba, Winnipeg, MB R3T 5V4, Canada
3
Department of Energy and Electrical Engineering, Nanchang University, Nanchang 330031, China
*
Author to whom correspondence should be addressed.
Energies 2021, 14(23), 8175; https://doi.org/10.3390/en14238175
Submission received: 5 November 2021 / Revised: 21 November 2021 / Accepted: 2 December 2021 / Published: 6 December 2021
(This article belongs to the Topic Application of Innovative Power Electronic Technologies)

Abstract

:
Kalman filter (KF) is often based on two models, which are phase angle vector (PAV) model and orthogonal vector (OV) model, in the application of distorted grid AC signal detection. However, these two models lack rigorous and detailed derivation from the principle of dynamic modeling. This paper presents a phase angle vector dynamic (PAVD) model and an orthogonal vector dynamic (OVD) model, which are combined with Kalman filter for detecting distorted grid AC signal. They reveal that the state noise covariance of the dynamic model−based KF is related to the sampling cycle, and overcome the defect of more detecting error for conventional model−based KF. Experiment and evaluation results show that the proposed KF algorithms are reasonable and effective. Therefore, this paper contributes a guiding significance for the application of KF algorithm in harmonic detection.

1. Introduction

The Kalman filter (KF) plays an increasingly important role in the real−time detection of grid AC signal [1]. It was first applied to the power signal processing field by Dr. A. A. Girgis in 1981 [2,3,4]. It mainly applies to the following fields: fundamental and harmonic components detection of grid voltage/current [5,6,7,8,9,10], phase−locked loop synchronization [11,12,13,14,15], power quality detection and compensation equipment [16,17,18,19,20], power disturbance feature extraction and machine classification [21,22,23,24], and etc. In these applications, two conventional models, phase angle vector (PAV) model and orthogonal vector (OV) model, are mainly used. The phase angles were used as vectors to establish a state-space model which could be called the phase angle vector (PAV) model through cosine expansion [3]. The orthogonal vector was used to establish a state-space model that could be called the orthogonal vector (OV) model [4]. However, the state noise covariance of the two models is only a given value based on empirical statistics of the AC signal. It is selected as a fixed value, such as Q = 0.05 × I (I represents the unit matrix), when AC signal is the grid voltage [10,11]. It is selected as a fixed value, such as Q = 0.01 × I, when the AC signal is the sampling current [20]. It is mainly to avoid big tracking error and enable to detect response without time delay. With the development of computer hardware and sensor technology, the sampling cycle is getting smaller, and the accuracy requirements are getting higher. The detecting accuracy of the fundamental and harmonics components for the conventional detecting model−based KF algorithm cannot meet the requirements of modern engineering applications.
The current issue is that the detecting accuracy of the fundamental and harmonics components for the conventional detecting model–based KF algorithm cannot meet the requirements of modern engineering applications. There is no fixed rule as references to find the fixed value by applying different sampling cycles. As a result, the detection accuracy cannot achieve the expected goal.
Also, the discretization model lacks of theoretical basis. In view of the fact that the previous literature did not detailedly derive the PAV and OV model from the principle of dynamic modeling. This paper proposes a phase angle vector dynamic (PAVD) model−based KF (PAVD−KF) algorithm and the orthogonal vector dynamic (OVD) model−based KF algorithm (OVD−KF) through model derivation according to the stochastic process theory [12]. It reveals that the state noise covariance of the conventional detecting model is related to the sampling cycle.
This paper contributes:
(1)
We derive the discretization model from the continuous differential equation in accordance with stochastic process theory, to find the value law of fixed value of covariance under different sampling cycles.
(2)
We greatly improve the detecting accuracy of fundamental and harmonics components, enabling the detection results to meet the requirements of modern engineering applications of AC current and voltage signal.
Section 1 introduces the Kalman filter mainly applies to the following fields and model−based KF existing problems. In Section 2, the related works about conventional PAV and OV model are introduced. In Section 3, the PAVD and OVD models are rederived according to of stochastic process theory. Then, the PAVD−KF and OVD−KF algorithms are proposed. In Section 4, the performances of PAVD−KF and OVD−KF algorithms are evaluated and compared by the experimental current and voltage sampling data. Section 5 concludes the paper.

2. Related Work

In this section, we proposed the expression of distorted AC signal first, and introduce the conventional PAV−KF and OV−KF algorithms then.

2.1. Distorted AC Signal

The waveform of the distorted AC signals sampled at time t can be expressed as
s ( t ) = i = 1 N A i ( t ) sin [ i ω 0 T + θ i ( t ) ] + i = 1 N w i ( t )
where i = 1 , 2 , 3 , , N , T denotes the sampling cycle. For the fundamental frequency of 50 Hz, ω 0 = 100 π . i = 1 N A i ( t ) sin [ i ω 0 T + θ i ( t ) ] represents the fundamental and harmonic components. i = 1 N w i ( t ) is the state noise of fundamental and harmonic components.

2.2. PAV Model

The vector of fundamental and harmonic components is expressed as X k = [ x 1 y 1 x n y n ] T , where x 1 = A 1 cos θ 1 , y 1 = A 1 sin θ 1 , x n = A n cos θ n , y n = A n sin θ n .
The measurement equation of model is expressed as [3]
z k = H k X k + V k = [ cos ( ω 0 k T ) sin ( ω 0 k T ) cos ( n ω 0 k T sin ( n ω 0 k T ) ] T [ x 1 y 1 x n y n ] k + V k
The state equation of model is expressed as
X k + 1 = Φ k X k + W k = [ 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 ] [ x 1 y 1 x n y n ] k + [ W x 1 W y 1 W x n W y n ] k
The covariance matrix of state noise W k can be expressed as
Q k = E [ W k W k T ] = σ P A V 2 I 2 n × 2 n

2.3. OV Model

Considering a signal A k sin ( ω k t k + θ k ) with amplitude A k , angular frequency is ω k and phase is θ k [4,11]. Let x k = A k sin ( ω k t k + θ k ) , y k = A k cos ( ω k t k + θ k ) , initially, consider A k + 1 / k A k , ω k + 1 / k ω k and θ k + 1 / k θ k . At the time t k + 1 = t k + T , the signal x k can be expressed as
x k + 1 = A k + 1 sin ( ω k t k + ω k T + θ k ) = x k cos ( ω k T ) + y k sin ( ω k T )
where, T is the sampling cycle. Additionally,
y k + 1 = A k + 1 cos ( ω k t k + ω k T + θ k ) = x k sin ( ω k T ) + y k cos ( ω k T )
the state−space is extended for detection of distorted AC signal existed n frequencies harmonic components.
The vectors are expressed as X k = [ x 1 y 1 x n y n ] k T , where, x 1 k = A 1 k sin ( ω k t k + θ 1 k ) , y 1 k = A 1 k cos ( ω k t k + θ 1 k ) , x n k = A n k sin ( n ω k t k + θ n k ) , y n k = A n k cos ( n ω k t k + θ n k ) .
Then, the state and measurement equation is expressed as [4]
{ X k + 1 = [ Φ 1 0 0 Φ n ] X k + W i k z k = [ 1 0 1 0 ] X k + V k
where, the state transition matrixes Φ i are shown below as
Φ i = [ cos ( i ω k T ) sin ( i ω k T ) sin ( i ω k T ) cos ( i ω k T ) ]
The covariance matrix of state noise W k i can be expressed as
Q i k = E [ W i k W i k T ] = σ O V 2 I 2 n × 2 n
in (4) and (9), I 2 n × 2 n is a unit matrix. σ P A V 2 and σ O V 2 are selected as 0.05, when AC signal is sampling voltage [10,11]. σ P A V 2 and σ O V 2 are selected as 0.01, when AC signal is sampling current [20].

2.4. PAV-KF and OV-KF Algorithm

The discrete PAV−KF and OV−KF algorithm for detecting distorted AC signal can be expressed as follow by combining PAV and OV model.
{ X ^ k + 1 / k + 1 = Φ k X ^ k + 1 / k + K k + 1 ( Z k + 1 H k X ^ k + 1 / k ) K k + 1 = Φ k P k + 1 / k H k T ( H k P k + 1 / k H k T + R k ) 1 P k + 1 / k + 1 = Φ k P k + 1 / k Φ k T K k + 1 H k P k + 1 / k Φ k T + Q k
The amplitude and the phase angle for the fundamental and harmonic components can be calculated as follows.
A i k + 1 = x ^ i k + 1 2 + y ^ i k + 1 2
ϕ i k + 1 = arctan ( x ^ i k + 1 y ^ i k + 1 )

3. PAVD-KF and OVD-KF Algorithm

In Section 3, we construct the AC signal continuous differential equation of phase angle vector and orthogonal vector, and derive the discretization PAVD and OVD model in accordance with stochastic process theory. Next, we combine the proposed models with KF algorithm, to obtain PAVD−KF and OVD−KF algorithm.

3.1. PAVFIGURE

D-KF Algorithm

The phase angle vectors of fundamental components in (2)−(3) are expressed as
{ x ( t ) = A cos θ y ( t ) = A sin θ
The first derivative can be expressed as
d d t [ x ( t ) y ( t ) ] = A [ x ( t ) y ( t ) ] + B [ w x ( t ) w y ( t ) ]
where, w x ( t ) and w y ( t ) are the input white noise, w x ( t ) w y ( t ) . w x ( t ) ( 0 , σ w 2 ) , w y ( t ) ( 0 , σ w 2 ) . Where, A = [ 0 0 0 0 ] , B = [ 1 0 0 1 ] .
Equation (11) is discretized, and the state equation can be obtained as follows [25].
[ x y ] k + 1 = Φ k [ x y ] k + W k
where,
Φ k = e A T = L 1 ( [ s I A ] 1 ) = [ 1 0 0 1 ]
W k = k T ( k + 1 ) T e A [ ( k + 1 ) T τ ] B w ( τ ) d τ
The covariance matrix of state noise W k can be expressed as
Q w = E [ W k W k T ] = σ w 2 [ T 0 0 T ] = σ w 2 T [ 1 0 0 1 ]
from (18), an important conclusion can be known that Q w is related to sampling cycle T .
The state−space dimension of the distorted AC signal is extended under the consideration of fundamental and harmonic components. The state-space formulas of PAVD model are the same as the conventional PAV model (2)−(3). The difference is that the state noise covariance in (4). The state noise covariance of the PAVD model is
Q k = σ P A V D 2 I 2 n × 2 n = ρ T I 2 n × 2 n
where, σ P A V D 2 is the state noise covariance of PAVD model. ρ 1 is a coefficient. When high detection precision is required, ρ is selected as a small value. When dynamic fast detecting is required, ρ is selected to a large value. Then, the discrete PAVD−KF algorithm can be obtained by combining model (2), (3), (15) and KF algorithm.

3.2. OVD−KF Algorithm

In the two−dimensional space, the circular movement of the object along the midpoint is projected on the x, y axis as a pure sine and a pure cosine waveform. It is shown in Figure 1.
It can be expressed as
{ x ( t ) = A m sin ( ω 0 t + θ ) y ( t ) = A m cos ( ω 0 t + θ )
where, A m is the amplitude, ω 0 is the angular frequency, and θ is the phase angle.
The first derivative of (10) can be expressed as
{ x ˙ ( t ) = ω 0 A m cos ( ω 0 t + θ ) + w x ( t ) = ω 0 y ( t ) + w x ( t ) y ˙ ( t ) = ω 0 A m sin ( ω 0 t + θ ) + w y ( t ) = ω 0 x ( t ) + w y ( t )
where, w x ( t ) and w y ( t ) are the input white noise, w x ( t ) w y ( t ) . w x ( t ) ( 0 , σ w 2 ) , w y ( t ) ( 0 , σ w 2 ) .
The state equation can be obtained from (21).
d d t [ x ( t ) y ( t ) ] = A [ x ( t ) y ( t ) ] + B [ w x ( t ) w y ( t ) ] = [ 0 ω 0 ω 0 0 ] [ x ( t ) y ( t ) ] + [ 1 0 0 1 ] [ w x ( t ) w y ( t ) ]
Equation (22) is discretized, and the state equation can be obtained as follows [25].
[ x y ] k + 1 = Φ k [ x y ] k + W k
where,
Φ k = e A T = L 1 ( [ s I A ] 1 ) = [ cos ( ω 0 T ) sin ( ω 0 T ) sin ( ω 0 T ) cos ( ω 0 T ) ]
W k = k T ( k + 1 ) T e A [ ( k + 1 ) T τ ] B w ( τ ) d τ
The state noise covariance matrix can be expressed as
Q w = E [ W k W k T ] = σ w 2 [ T 0 0 T ] = σ w 2 T [ 1 0 0 1 ]
from Equation (16), an important conclusion can be drew that Q w is related to the sampling cycle T .
The state-space dimension of the distorted AC signal is expanded under the consideration of fundamental and harmonic components. The state and measurement equations of OVD model are the same as the conventional OV model (4) and (5). The difference is that state noise covariance matrix of the OV model in (6). That is
Q k = σ O V D 2 I 2 n × 2 n = ρ T I 2 n × 2 n
where σ O V D 2 is the covariance of OVD model. ρ 1 is a coefficient. When high detection precision is required, ρ is selected as a small value. When dynamic fast detecting is required, ρ is selected to a large value.
Then, there are the combinations of the model (7), (8), (27), and KF algorithm. The discrete OVD−KF algorithm can be obtained.

4. Experiment and Evaluation

In this section, we conduct the detection of the current and the voltage with the sampling cycle of 10 μs and 40 μs. Then, comparing the detecting results for fundamental and harmonics components in PAVD−KF and PAV−KF algorithm, and in OVD−KF and OV-KF algorithm. At the last, we come to a conclusion of these results.

4.1. Experiment and Evaluation for AC Current Detection

4.1.1. Experiment Settings

As an experimental example, we sample the AC distorted current signal. Experiment platform consists of AC grid (380 V, 50 Hz), oscilloscope (Tektronix GDS−2102) and nonlinear load (the three−phase rectifier, R = 100 Ω and L = 1.5 mH). It is shown in Figure 2. The AC current signals are sampled by the oscilloscope. Sampling cycles are selected as 10 µs and 40 µs. The experiment sampling current data and spectrogram are shown in Figure 3 and Figure 4. The HRI of harmonic components are measured to be approximately 24% of 5th, 10% of 7th, 10% of 11th, 5% of 13th, 5% of 17th, 4% of 19th by spectrum analysis.
HRI is defined as follows
HRI n = I n I 1 × 100 %
where I 1 is the voltage rms value of fundamental component, I n is the voltage rms value of n−th harmonic components.
The measurement noise covariance R of the algorithms is all selected as 0.06 A, that is the measurement error of the sensor. The initial covariance matrices are selected as P 0 / 0 = 1000 × I 14 × 14 . The state noise covariance is selected as Q k = 0.01 × I 14 × 14 in (4). The state noise covariance is selected as Q k = 0 . 01 × T × I 14 × 14 in (19).

4.1.2. Detection Results for PAVD−KF Algorithm

Two algorithms of PAV−KF and PAVD−KF are evaluated and compared. The evaluation results for sampling cycle 10 µs and 40 µs are shown in Figure 5 and Figure 6.

4.1.3. Detection Results for OVD−KF Algorithm

Two algorithms of OV−KF and OVD−KF are evaluated and compared. The evaluation results for sampling cycle 10 µs and 40 µs are shown in Figure 7 and Figure 8.

4.1.4. Summary for AC Current Detecting

It can be seen from Figure 5 and Figure 8, when Q = 0.01 × I, the detecting errors of the conventional PAV−KF and OV−KF algorithms are big. In particular, the detecting errors of the conventional algorithms for sampling cycle 40 μs are rather big. The errors are too big to affect the normal operation.
Corresponding to it, when Q combines with the sampling cycle, Q = 0.01 × T × I, regardless of the fundamental or harmonic components, the proposed detecting accuracies in PAVD−KF and OVD−KF algorithm are much higher than the conventional algorithms.

4.2. Experiment and Evaluation for AC Voltage Detection

4.2.1. Experiment Settings

The AC signal is grid voltage as an experimental example. Experiment platform is shown in Figure 9. The experiment of superimposing harmonics in fundamental voltage is carried out through a programmable power supply device. The AC voltage signals are sampled and the sampling cycles T are selected by the oscilloscope of Tektronix GDS−2102.
The harmonics ratio for voltage (HRU) is defined as follows
HRU n = U n U 1 × 100 %
where U 1 is the voltage rms value of the fundamental component, and U n is the voltage rms value of n−th harmonic components.
Experiment voltage data for sampling cycle 10 µs is in Figure 10. Experiment voltage data for sampling cycle 40 µs is in Figure 11, and the HRU of harmonic components all are set as 5% of 3rd, 10% of 5th, and 8% of 7th.
The R of the algorithms are all selected as 3.11 V, that is the measurement error of the sensor. The initial covariance matrices of algorithms are selected as P 0 / 0 = 1000 × I 10 × 10 . The state noise covariance matrix of algorithms is selected as Q k = 0.05 × I 10 × 10 . The state noise covariance matrix of algorithms is Q k = 0 . 1 × T × I 10 × 10 .

4.2.2. Detection Results for PAVD−KF Algorithm

PAVD−KF and PAV−KF algorithms are evaluated and compared. The detecting results of the sampling cycle 10 µs and 40 µs are shown in Figure 12 and Figure 13.

4.2.3. Detection Results for OVD−KF Algorithm

OVD−KF and OV−KF algorithms are evaluated and compared. The detecting results of the sampling cycle 10 µs and 40 µs are shown in Figure 14 and Figure 15.

4.2.4. Summary for AC Voltage Detecting

It can be seen from Figure 12 and Figure 13 that the amplitudes of the fundamental components of the conventional model−KF algorithm detecting are fluctuating, the stability of its detection shows poorly. Specifically, when we use the sampling cycle 40 µs, the error is too big to be ignored. From Figure 14 and Figure 15, we know that detecting errors of the harmonics components for the conventional OV−KF algorithm are also big, though the error seems smaller compared to PAV−KF. Corresponding to it, the detecting accuracies of the fundamental and harmonic components for the OVD−KF algorithm are higher than that of the OV−KF algorithm.
Generally, when Q = 0.05 × I, the conventional algorithms could basically meet the engineering requirements with some big errors under different sampling cycle conditions. When Q = 0.05 × T × I, all the novel algorithms could fulfill the requirements.

4.3. Summary for Experiment and Evaluation

In general, when using the fixed value, the error of PAV model is the biggest, neither current nor voltage can meet the demand. OV model performs better, when using the sampling cycle 10 µs, the voltage can be used normally, and the current error is too big to be used. But when using the sampling cycle 40 µs, both current and voltage cannot operate normally.
The proposed algorithms in this paper, whether it is PAVD model or OVD model, it is related to the sampling cycle, so the voltage and current detecting results can meet the engineering requirement.
Take PAV model with 40 μs in current and voltage tracking for example, the conventional methods produce big errors when Q = 0.01 × I detecting the current and Q = 0.05 × I detecting the voltage. Only when Q = 0.00001 × I detecting the current, the accuracy of it would be almost close to that of our proposed algorithm, as Figure 16 shows. When Q = 0.0005 × I detecting the voltage, the accuracy would satisfy the requirement, as Figure 17 shows.

5. Conclusions

This paper presents the PAVD−KF and OVD−KF algorithm derived from the stochastic process theory. It enhances the knowledge about the Kalman algorithms applied to the harmonic detection of the distorted AC grid. It gives a rule that the KF algorithms in the application of distorted grid AC signal detection is related to the sampling cycle and their state noise covariance cannot be simply considered as a constant.
This paper applies the actual current and voltage data of the sampling cycle 10 µs and 40 µs to detection and evaluation. The conventional algorithms have big errors in general and unable to meet engineering requirements. When using a fixed value of covariance, if we are going to obtain the ideal results, we need to spend a lot of time debugging (trying to make up) the parameters. Also, there is no theoretical source to set the fixed value of covariance. To solve it, we propose the algorithms to save much time and workload to get the ideal results. Besides, we make the algorithms more feasible and effective based on stochastic process theory. From which the results show that the proposed PAVD−KF and OVD−KF algorithms are effective and improve the dynamic detecting accuracy of grid AC signal significantly. This paper would provide a reference for the application of KF algorithms to harmonic detection, power quality control and grid synchronization in the future.

Author Contributions

Conceptualization, H.N. and X.N.; methodology, H.N. and X.N.; software, H.N.; validation, H.N. and X.N.; formal analysis, H.N.; investigation, H.N.; resources, H.N.; data curation, H.N.; writing—original draft preparation, H.N.; writing—review and editing, H.N.; visualization, H.N.; supervision, X.N.; project administration, X.N.; funding acquisition, X.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundations of China, grant number 51867017 “Research on Power Quality Dynamic Tracking and Real-time Detection Method for Micro-grid Complex Source Load Environment”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Auger, F.; Hilairet, M.; Guerrero, J.M.; Monmasson, E.; Orlowska-Kowalska, T.; Katsura, S. Industrial Applications of the Kalman Filter: A Review. IEEE Trans. Ind. Electron. 2013, 60, 5458–5471. [Google Scholar] [CrossRef] [Green Version]
  2. Girgis, A.A.; Brown, R.G. Application of Kalman Filtering In Computer Relaying. IEEE Trans. Power Appar. Syst. 1981, 7, 3387–3397. [Google Scholar] [CrossRef]
  3. Girgis, A.A.; Hwang, T.L.D. Optimal Estimation of Voltage Phasors and Frequency Deviation. IEEE Trans. Power Appar. Syst. 1984, 10, 2943–2951. [Google Scholar] [CrossRef]
  4. Girgis, A.A.; Sallam, A.A.; El-Din, A.K. Identification and tracking of harmonic sources in a power system using a Kalman filte. IEEE Trans. Power Deliv. 1998, 13, 414–420. [Google Scholar] [CrossRef]
  5. Nie, X. Detection of Grid Voltage Fundamental and Harmonic Components Using Kalman Filter Based on Dynamic Tracking Model. IEEE Trans. Ind. Electron. 2020, 67, 1191–1200. [Google Scholar] [CrossRef]
  6. Naseri, F.; Kazemi, Z.; Farjah, E.; Ghanbari, T. Fast Detection and Compensation of Current Transformer Saturation Using Extended Kalman Filter. IEEE Trans. Power Deliv. 2020, 34, 1087–1097. [Google Scholar] [CrossRef]
  7. Swain, S.; Subudhi, B. Grid Synchronization of a PV System with Power Quality Disturbances Using Unscented Kalman Filtering. IEEE Trans. Sustain. Energy 2019, 10, 1240–1247. [Google Scholar] [CrossRef]
  8. Uz-Logoglu, E.; Salor, O.; Ermis, M. Real-Time Detection of Interharmonics and Harmonics of AC Electric Arc Furnaces on GPU Framework. IEEE Trans. Ind. Appl. 2019, 55, 1373–1378. [Google Scholar]
  9. Xi, Y.; Tang, X.; Li, Z.; Cui, Y.; Zhao, T.; Zeng, X.; Guo, J.; Duan, W. Harmonic estimation in power systems using an optimised adaptive Kalman filter based on PSO-GA. IET Gener. Transm. Distrib. 2019, 13, 3968–3979. [Google Scholar] [CrossRef]
  10. Bagheri, A.; Mardaneh, M.; Rajaei, A.; Rahideh, A. Detection of Grid Voltage Fundamental and Harmonic Components Using Kalman Filter and Generalized Averaging Method. IEEE Trans. Power Electron. 2016, 31, 1064–1073. [Google Scholar] [CrossRef]
  11. Cardoso, R.; Camargo, R.F.; Pinheiro, H.; Grundling, H.A. Kalman filter based synchronisation methods. IET Gener. Transm. Distrib. 2008, 2, 542–555. [Google Scholar] [CrossRef]
  12. Swain, S.; Subudhi, B. Iterated extended Kalman filter-based grid synchronisation control of a PV system. IET Energy Syst. Integr. 2019, 1, 219–228. [Google Scholar] [CrossRef]
  13. Ferrero, R.; Pegoraro, P.A.; Toscani, S. Synchrophasor Estimation for Three-Phase Systems Based on Taylor Extended Kalman Filtering. IEEE Trans. Instrum. Meas. 2020, 69, 6723–6730. [Google Scholar] [CrossRef]
  14. Golestan, S.; Guerrero, J.M.; Vasquez, J.C.; Abusorrah, A.M.; Al-Turki, Y. Single-Phase FLLs Based on Linear Kalman Filter, Limit-Cycle Oscillator, and Complex Bandpass Filter: Analysis and Comparison With a Standard FLL in Grid Applications. IEEE Trans. Power Electron. 2019, 34, 1064–1073. [Google Scholar] [CrossRef]
  15. Ahmed, H.; Biricik, S.; Benbouzid, M. Linear Kalman Filter-Based Grid Synchronization Technique: An Alternative Implementation. IEEE Trans. Ind. Inform. 2021, 17, 3847–3856. [Google Scholar] [CrossRef]
  16. Alam, S.J.; Arya, S.R. Control of UPQC based on steady state linear Kalman filter for compensation of power quality problems. Chin. J. Electr. Eng. 2020, 6, 52–65. [Google Scholar] [CrossRef]
  17. Xi, Y.; Li, Z.; Zeng, X.; Tang, X.; Liu, Q.; Xiao, H. Detection of power quality disturbance using an adaptive process noise covariance Kalman filte. Digit. Signal Process. 2018, 76, 34–49. [Google Scholar] [CrossRef]
  18. Wu, C.; Magana, M.E.; Cotilla-Sanchez, E. Dynamic Frequency and Amplitude Estimation for Three-Phase Unbalanced Power Systems Using the Unscented Kalman Filter. IEEE Trans. Instrum. Meas. 2019, 68, 3387–3395. [Google Scholar] [CrossRef]
  19. Panigrahi, R.; Subudhi, B.; Panda, P.C. A robust LQG servo control strategy of shunt-active power filter for power quality enhancement. IEEE Trans. Power Electron. 2016, 31, 2860–2869. [Google Scholar] [CrossRef]
  20. Panigrahi, R.; Subudhi, B. Performance enhancement of shunt active power filter using a Kalman filter-based H∞ control strategy. IEEE Trans. Power Electron. 2017, 32, 2622–2630. [Google Scholar] [CrossRef]
  21. Abdelsalam, A.A.; Abdelaziz, A.Y.; Kamh, M.Z. A Generalized Approach for Power Quality Disturbances Recognition Based on Kalman Filter. IEEE Access 2021, 9, 93614–93628. [Google Scholar] [CrossRef]
  22. Xi, Y.; Li, Z.; Tang, X.; Zeng, X. Classification of power quality disturbances based on KF-ML-Aided S-Transform and multilayers feedforward neural networks. IET Gener. Trans. Distrib. 2020, 14, 4010–4020. [Google Scholar] [CrossRef]
  23. Abdelsalam, A.A.; Eldesouky, A.A.; Sallam, A.A. Classification of Power System Disturbances using Linear Kalman Filter and Fuzzy-expert System. Electr. Power Energy Syst. 2012, 43, 688–695. [Google Scholar] [CrossRef]
  24. Abdelsalam, A.A.; Eldesouky, A.A.; Sallam, A.A. Characterization of power quality disturbances using hybrid technique of linear Kalman filter and fuzzy-expert system. Electr. Power Syst. Res. 2012, 83, 41–50. [Google Scholar] [CrossRef]
  25. Brown, R.G.; Hwang, P.Y.C. Introduction to Random Signals and Applied Kalman Filtering with MATLAB Exercises, 4th ed.; John Wiley& Sons, Inc.: Hoboken, NJ, USA, 2011. [Google Scholar]
Figure 1. Diagram of the object circular movement being projected.
Figure 1. Diagram of the object circular movement being projected.
Energies 14 08175 g001
Figure 2. Experiment platform. (a) Current sampling experiment platform (b) Principle diagram of experiment platform.
Figure 2. Experiment platform. (a) Current sampling experiment platform (b) Principle diagram of experiment platform.
Energies 14 08175 g002
Figure 3. Experiment current for sampling cycle 10 µs. (a) The sampling current data. (b) Spectrogram of sampling current data.
Figure 3. Experiment current for sampling cycle 10 µs. (a) The sampling current data. (b) Spectrogram of sampling current data.
Energies 14 08175 g003
Figure 4. Experiment current for sampling cycle 40 µs. (a) The sampling current data. (b) Spectrogram of sampling current data.
Figure 4. Experiment current for sampling cycle 40 µs. (a) The sampling current data. (b) Spectrogram of sampling current data.
Energies 14 08175 g004
Figure 5. PAVD−KF current detecting results for sampling cycle 10 µs. (a) Fundamental components detecting results (b) HRI detecting results.
Figure 5. PAVD−KF current detecting results for sampling cycle 10 µs. (a) Fundamental components detecting results (b) HRI detecting results.
Energies 14 08175 g005
Figure 6. PAVD−KF current detecting results for sampling cycle 40 µs. (a) Fundamental components detecting results (b) HRI detecting results.
Figure 6. PAVD−KF current detecting results for sampling cycle 40 µs. (a) Fundamental components detecting results (b) HRI detecting results.
Energies 14 08175 g006
Figure 7. OVD−KF current detecting results for sampling cycle 10 µs. (a) Fundamental components detecting results (b) HRI detecting results.
Figure 7. OVD−KF current detecting results for sampling cycle 10 µs. (a) Fundamental components detecting results (b) HRI detecting results.
Energies 14 08175 g007
Figure 8. OVD−KF current detecting results for sampling cycle 40 µs. (a) Fundamental components detecting results (b) HRI detecting results.
Figure 8. OVD−KF current detecting results for sampling cycle 40 µs. (a) Fundamental components detecting results (b) HRI detecting results.
Energies 14 08175 g008
Figure 9. Experiment platform. (a) Voltage sampling experiment platform. (b) Principle diagram of experiment platform.
Figure 9. Experiment platform. (a) Voltage sampling experiment platform. (b) Principle diagram of experiment platform.
Energies 14 08175 g009
Figure 10. Experiment voltage for sampling cycle 10 µs. (a) The sampling voltage data. (b) Spectrogram of sampling voltage data.
Figure 10. Experiment voltage for sampling cycle 10 µs. (a) The sampling voltage data. (b) Spectrogram of sampling voltage data.
Energies 14 08175 g010
Figure 11. Experiment voltage for sampling cycle 40 µs. (a) The sampling voltage data. (b) Spectrogram of sampling voltage data.
Figure 11. Experiment voltage for sampling cycle 40 µs. (a) The sampling voltage data. (b) Spectrogram of sampling voltage data.
Energies 14 08175 g011
Figure 12. PAVD−KF voltage detecting results for sampling cycle 10 µs. (a) Fundamental components detecting results (b) HRU detecting results.
Figure 12. PAVD−KF voltage detecting results for sampling cycle 10 µs. (a) Fundamental components detecting results (b) HRU detecting results.
Energies 14 08175 g012
Figure 13. PAVD−KF voltage detecting results for sampling cycle 40 µs. (a) Fundamental components detecting results (b) HRU detecting results.
Figure 13. PAVD−KF voltage detecting results for sampling cycle 40 µs. (a) Fundamental components detecting results (b) HRU detecting results.
Energies 14 08175 g013
Figure 14. OVD−KF voltage detecting results for sampling cycle 10 µs. (a) Fundamental components detecting results (b) HRU detecting results.
Figure 14. OVD−KF voltage detecting results for sampling cycle 10 µs. (a) Fundamental components detecting results (b) HRU detecting results.
Energies 14 08175 g014
Figure 15. OVD−KF voltage detecting results for sampling cycle 40 µs. (a) Fundamental components detecting results (b) HRU detecting results.
Figure 15. OVD−KF voltage detecting results for sampling cycle 40 µs. (a) Fundamental components detecting results (b) HRU detecting results.
Energies 14 08175 g015
Figure 16. PAVD−KF current detecting results for sampling cycle 40 µs. (Q = 0.00001 × I).
Figure 16. PAVD−KF current detecting results for sampling cycle 40 µs. (Q = 0.00001 × I).
Energies 14 08175 g016
Figure 17. PAVD−KF voltage detecting results for sampling cycle 40 µs. (Q = 0.0005 × I).
Figure 17. PAVD−KF voltage detecting results for sampling cycle 40 µs. (Q = 0.0005 × I).
Energies 14 08175 g017
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nie, H.; Nie, X. Research on Dynamic Modeling of KF Algorithm for Detecting Distorted AC Signal. Energies 2021, 14, 8175. https://doi.org/10.3390/en14238175

AMA Style

Nie H, Nie X. Research on Dynamic Modeling of KF Algorithm for Detecting Distorted AC Signal. Energies. 2021; 14(23):8175. https://doi.org/10.3390/en14238175

Chicago/Turabian Style

Nie, Haoyao, and Xiaohua Nie. 2021. "Research on Dynamic Modeling of KF Algorithm for Detecting Distorted AC Signal" Energies 14, no. 23: 8175. https://doi.org/10.3390/en14238175

APA Style

Nie, H., & Nie, X. (2021). Research on Dynamic Modeling of KF Algorithm for Detecting Distorted AC Signal. Energies, 14(23), 8175. https://doi.org/10.3390/en14238175

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop