Open access peer-reviewed chapter

Kalman Filters for Parameter Estimation of Nonstationary Signals

Written By

Sarita Nanda

Submitted: 22 April 2017 Reviewed: 23 October 2017 Published: 21 February 2018

DOI: 10.5772/intechopen.71874

From the Edited Volume

Kalman Filters - Theory for Advanced Applications

Edited by Ginalber Luiz de Oliveira Serra

Chapter metrics overview

1,756 Chapter Downloads

View Full Metrics


An adaptive Taylor-Kalman filter with PSO tuning for tracking nonstationary signal parameters in a noisy environment with primary focus on time-varying power signals has been presented in this piece of work. In order to deal with the dynamic envelope of the power signal, second-order Taylor expansion has been used such that the Taylor coefficients are updated with the PSO-tuned Taylor-Kalman Filter algorithm. In addition to this, for fast convergence, a self-adaptive particle swarm optimization technique has been used for obtaining the optimal values of model and measurement error covariances of the Kalman filter. The proposed algorithm is linear and therefore has less computational burden, which is easier to be implemented on a hardware platform like DSP processor or FPGA. The proposed PSO-tuned Taylor-Kalman filter exhibits robust tracking capabilities even under changing signal dynamics, immune to critical noise conditions, harmonic contaminations, and also reveals excellent convergence properties.


  • nonstationary signals
  • power signal frequency and phasor estimation
  • hybrid Kalman approach
  • PSO tuning

1. Introduction

Signal parameter estimation, which dates back to the late 19th century, describes the various methods employed to track amplitude, phase, and frequency-like parameters of a signal. Among all the signal parameters, frequency is the primary concern, as it is a nonlinear function in the received data sequence, and once that is measured accurately, tracking of other parameters like phase, amplitude, and damping factor of a signal can be relatively easier. Most real-world signals are nonstationary in nature, i.e., they have a time-varying frequency behavior. Some of the popular sources of nonstationary signals include speech, audio, sounds of mammals, machine vibrations, electrical power networks, and a variety of biomedical signals like electromyogram (EMG), electroencephalogram (EEG), phonocardiogram (PCG), and vibroarthrogram (VAG)). These signals are rich in information and when analyzed properly provides with information that could be used to improve many aspects of our lives. Hence, the information of interest of the signal can be extracted, which includes the estimation of parameters like amplitude, phase, frequency, and damping factor directly from the discrete measurement in the presence of noise both in stationary and nonstationary environments. Precise and smooth operation of the power generation and distribution system is very much required in the present day scenario. With the increasing demand for power, the number and type of load are having deteriorating effects on the power quality. Power quality is defined as the ability of the electrical grid to deliver clean and stable power to the consumer. Between generation and supply, the power being delivered encounters large number of transformers and several lengths of overhead lines and underground cables. Phenomena like lightning strikes, system faults, load switching, and other such intentional or unintentional events are the main cause of electromagnetic disturbances, which results in voltage or current waveform distortions to propagate in the entire power system. Recently, the increase in the number of power electronic loads in the system causes nonlinear loading effect on the power system signal, leading to degradation of power quality.

Recently, harmonic estimation has become a challenging and critical issue for electrical engineers. Estimating harmonics and other faults is important for maintaining power quality. Research works carried out recently sheds light on various techniques for estimating harmonics. FFT [1]-based techniques are the conventional ones, and they suffer from some pitfalls such as aliasing and picket fence effects, which lead to inaccurate estimation results. There are some other methods suffering from these three problems, and this is because of existing high frequency components measured in the signal; however, truncation of the sequence of sampled data, when only a fraction of the sequence of a cycle exists in the analyzed waveform, can boost leakage problem of the DFT method. So, the need of new algorithms that process the data, sample-by-sample and not in a window as in FFT and DFT, is of paramount importance. Another very robust algorithm for the purpose of estimating sinusoidal signals with unknown noise content is the Kalman filter (KF) [2, 3].

However, when cases related to system dynamics, like sudden changes in frequency, amplitude and phase of a signal, arise, KF exhibits serious drawbacks. Study of several literature shows that single methods employed for the purpose of signal estimation are not efficient on their own, so hybrid methods based on the combination of different need to be formulated. The major contribution of this chapter is the accurate tracking of nonstationary power signal parameters, i.e., phasor, frequency, and harmonics. The power signal is modeled using Taylor series, and the coefficients of the Taylor series are updated using the Kalman Filter [4, 5], which are again utilized to estimate the time varying amplitude, phase, and frequency of the test signal. Moreover, a self-adaptive particle swarm optimization approach is deployed to choose the optimum values of the Kalman filter parameters like model and measurement error covariances, which in turn enables the filter to attain convergence in a faster rate.


2. Literature review

Work on harmonic and parameter estimation has been going since the introduction of AC power generation. Over the course of time, several methods have been proposed to fulfill this particular requirement, but so far, the existing methodologies have exhibited significant drawbacks. Here are some of the research works that have been performed in the last 5 years.

For the estimation of harmonics and interharmonics, a technique using simple techniques like least mean square [4] and a two-stage ADALINE network has been studied [5]. The method utilized here provided a better accuracy even when power frequency deviation and interharmonic components are present in the measured signal. As the conventional ADALINE is unable to detect interharmonics, a two-stage ADALINE is used. The architecture is classified in two parts—the front stage that extracts the frequency value and the back stage that computes amplitude and phase. Here, the adaptive algorithm used in the filter is the RLS algorithm. The method yielded more accurate results in protection and monitoring applications.

Sliding window tracking (SWT) [6] accurately tracks the frequency and amplitude of a signal by processing only three (or more) recent data points. It works for a signal with any nonzero moving average and noise. Teager-Kaiser algorithm (TKA) is a well-known four-point method for online tracking of frequency and amplitude. TKA takes into assumption that the signal is purely harmonic, so any moving average in the signal can totally destroy the accuracy of TKA, whereas SWT uses a pair of windowed regular harmonics to estimate the frequency and amplitude thus eliminating the effect of moving average. In order to start the online tracking of frequency, SWT requires TKA to provide the first estimate of the frequency. The accuracies of SWT and TKA are compared using Hilbert-Huang transform, which is used to extract accurate time-varying frequencies and amplitudes by processing the whole data set without assuming the signal to be harmonic. Tracking accuracy increases when window length is equal to or greater than one quarter of the signal period. If the chosen window length is too long, then the estimated frequency is an average over the window length. The method requires constant frequency and amplitude to accurately track the parameters, and this shows that the dynamic response of the method is very poor and the accuracy deteriorates, if there is no change in the parameter values.

A real-time approach for the estimation of power system frequency based on Newton-type algorithm and least squares method has been used in this paper [7]. The adopted optimization technique has been based on a two-stage mathematical model. A Newton type algorithm has been used to model the first stage for estimating the line to neutral voltage-phase angle and its variation. The second stage has been modeled using LS minimization technique that extracts the power system frequency by processing the information in the phase angles estimated using NTA. The method also studies the modulating effect of time-varying frequency on the online estimation of the phase angle.

Taylor series expansion and Fourier algorithm have been used for frequency estimation [8]. To model the changing envelope of a power signal within an observation, a second-order Taylor series has been used, and the parameters of the model have been estimated using Fourier algorithm. Comparing with the traditional Fourier algorithm, this method introduces more computational load.

A modified ADALINE structure has been used in the paper [9] for online tracking of harmonics. Self-synchronized ADALINE network for power system harmonics estimation relies on the Levenberg Gradient Descent method for updating the system parameters. A faster response and better noise immunity are provided by conventional methods. A high computational load is the only drawback that exists in the proposed approach.

Ensemble Kalman Filter has been used in the proposed method [10] for filtering and estimating signal harmonics and interharmonics. To avoid the problem of singularity and for the computational feasibility of state covariance P, the state covariance P is replaced by a sample covariance C for the computation of Kalman gain.

The proposed method [11] is adopted for real-time estimation of phasor and harmonics. The technique reduces the turnaround time on two different off-the-shelf research and development DSP platforms. The proposed method has been found to be superior to that of ADALINE and RDFT techniques under the presence of noise sub-harmonics and frequency variations. The proposed technique has a computational efficiency that is higher than that of ADALINE and RDFT techniques.

The proposed algorithm in [12] is simple, computational efficient and makes the correction of the signal that enables to reach the mean square error. It provides a new kind of step adaptation for LMS algorithm. Two LMS algorithms have been utilized by this method. The first one has a fixed-step size, and the weight coefficient generated from the first algorithm is used to update the step size of the second algorithm, which has initial step size of 0.001.

An adaptive linear network (ADALINE) [13, 14, 15] for harmonic and interharmonic estimation (Martin) allows the computing of root mean square voltage and total harmonic distortion indices. Classification and detection of sags, swells, outages, and harmonics-interharmonics have been done using the indices computed before. Classification of spikes, notching, flicker, and oscillatory transients has been achieved by using a feed forward neural network through pattern recognition using horizontal and vertical histograms of a specific voltage waveform. The method used in [16] uses noneven item interpolation FFT based on triangular self-convolution window. Variances of frequency estimation are proportional to the energy of the adopted window. By choosing suitable values of length of FFT, sampling frequency, and the shape of the adopted window, the variances of frequency estimation have been determined.


3. PSO-tuned Taylor-Kalman filter

To improve the performance of Kalman Filter in this aspect, a hybrid adaptive filter has been proposed in this thesis work that consists of the combination of Taylor series, Kalman Filter, and self-adaptive PSO. Taylor series is used to model the changing envelope of the sinusoidal signal. The sinusoidal signal is expressed in its trigonometric components, which in turn are expanded by using Taylor series. The Taylor coefficients are stored in the state vector that is further used to estimate the signal and its amplitude, frequency, and phase. In each iteration, the state vector is updated in order to get a better estimate than the previous, and the process continues until convergence is reached. There are two parameters on which the performance of the KF depends—the model and measurement error covariances. In the traditional approach, the values for these parameters are chosen by trial and error that makes the algorithm time consuming and prone to errors. Self-adaptive PSO is used here to select the optimal values of the error covariances in order to achieve fast convergence.

3.1. Signal modeling using Taylor expansion

Let the discrete signal be represented as:

v i = A i . cos i + ϕ i + κ i E1

where A i , ω i , and ϕ i are “the amplitude”, “angular frequency,” and “phase” of the sinusoid, respectively. ω i = 2 πf i and f i is the fundamental frequency of the signal, while κ i is an additive white noise with unknown variance σ g 2 Now let us represent θ i = 2 πfidt + ϕ i . The rate of change of phase angle is equal to frequency. So the signal frequency can be represented as [3]:

f = 1 2 π d dt θ i = f 0 + 1 2 π d dt ϕ i E2

Eq. (1) can be expressed according to trigonometric function as:

v i = Q i cos 2 πf i R i sin 2 πf i E3

where Q i = A i cos ϕ i and R i = A i sin ϕ i .

The coefficient functions Q i and R i express the envelope of the time varying sinusoid and can be expanded using Taylor series [17, 18] as shown:

Q i m 0 + m 1 i + m 2 i 2 + and R i n 0 + n 1 i + n 2 i 2 + E4


m 0 = Q 0 , m 1 = dQ i dt at k = 0; m 2 = d 2 Q i dt 2 and m 3 = d 3 Q i dt 3 at k = 0.

n 0 = R 0 , n 1 = dR i dt at k = 0, n 2 = d 2 R i dt 2 and n 3 = d 3 R i dt 3 at k = 0.

Now we can obtain the amplitude and phase angle of the described given sinusoid using Eq. (3) and (4) as follows at k = 0:

a ̂ = m 0 2 + n 0 2 E5


ϕ ̂ = arctan n 0 / m 0 E6

where m 0 = A 0 . cos ϕ 0

n 0 = A 0 . sin ϕ 0 E7

Similarly for estimating the frequency of the given sinusoid, consider Eq.(4) at k = 0, the first derivative will be:

m 1 = d dt A 0 . cos ϕ 0 E8
n 1 = d dt A 0 . sin ϕ 0 E9

By substituting Eq. (7) in Eq. (8) and (9) and by neglecting

d dt A 0 : d dt ϕ 0 = m 0 n 1 n 0 m 1 m 0 2 + n 0 2 E10

Now from Eq. (2) and Eq. (10), we get the formula for computing the frequency:

f ̂ = f 0 . + 1 2 π m 0 n 1 n 0 m 1 m 0 2 + n 0 2 E11

3.2. Updation of Taylor coefficients using the PSO-tuned Kalman filtering algorithm

Let us consider the following discrete signal:

Y i = a . sin T S + ϕ + n k E12

where a , T s , ω , and ϕ are the amplitude, sampling time, angular frequency, and phase of the signal, respectively, and n k represents measurement noise with a covariance R.

We can represent the state space Eq. (10) of the discrete signal as:

χ ^ i = f i χ i + η i E13
χ i 1 = m 0 ; χ i 2 = m 1 ; χ i 3 = m 2 ; χ i 4 = n 0 ; χ i 5 = n 1 ; χ i 6 = m 2 E14

And the state transition matrix is given by:

f i = 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 E15

The stochastic model of the signal is obtained as

χ 1 i χ 2 i χ 3 i χ 4 i χ 5 i χ 6 i = 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 χ 1 i χ 2 i χ 3 i χ 4 i χ 5 i χ 6 i E16

The measurement model of the signal expressed in Eq. (12) can be calculated as:

S i = h i χ i + v i E17

where the observation matrix can be calculated as:

H i = sin 2 π f 0 idt idt sin 2 π f 0 idt cos 2 π f 0 idt idt cos 2 π f 0 idt E18

The error signal can be obtained as

E i = S H i χ ̂ i E19

Using Eq. (19) the updated state estimate can be obtained from the following equation

χ ^ i = χ ^ i 1 + K i S i H χ ^ i E20

where the Kalman gain K(i) is given as:

K i = P ̂ i 1 H i T H P ̂ i 1 i H i T + r 1 E21

where P ̂ i is the covariance matrix given by

P ̂ i = P ̂ i 1 K i H i P ̂ i 1 E22


P ̂ i + 1 = P ̂ i + q E23

where q is the model noise covariance matrix given by

q = q a 0 0 0 0 0 0 q b 0 0 0 0 0 0 q c 0 0 0 0 0 0 q d 0 0 0 0 0 0 q e 0 0 0 0 0 0 q f E24

r is the measurement noise covariance which is fine tuned by using the error between the desired and estimated signals for

r i = λ g r + E i ^ 2 E25

where λ g is the forgetting factor in the range (0.9–1). Finally using the EKF time updated equations the X i matrix is computed which determines the values of the Taylor series coefficients m 0, m 1, n 0, and n 1.

3.3. Particle swarm optimization-based tuning of the Kalman filter

For fast convergence, optimum values of q and r are selected by applying a self-adaptive PSO [19]. For this purpose, a cost function is formulated, which passed in the PSO algorithm to get the optimum value for q and r. Here, the cost function is:

F = 1 L i = 1 L E i 2 E26

Particle swarm optimization is used to minimize the value of Eq. (26). Each particle is characterized by two attributes:

  1. pbest or Personal best: it holds the best value of position with respect to the previous positions of the particular particle.

  2. gbest or Global best: it holds the best value of position in the entire search space.

The PSO algorithm either minimizes or maximizes the value of gbest. Let x ij and V ij be the position and velocity of the i th particle in the j th dimension at k th instance of time. The personal best value can be determined from Eq. (27).

pbest i k + 1 = pbest i k , if F x i k + 1 > F pbest i k x i k + 1 , if F x i k + 1 < F pbest i k E27

where F indicates the cost function. The value of global best is obtained as:

gbest k = min C .. F pbest 0 k C .. F pbest 1 k C .. F pbest 2 k . C .. F pbest s k E28

For each particle, the updated velocity and position at time k + 1 are given by

V i k + 1 = K α v i k + b 1 rnd 1 pbest i k x i k + b 2 rnd 2 gbest i k x i k E29
x i k + 1 = x i k + V i k + 1 E30

where α is the inertia weight factor, b 1 and b 2 are the acceleration constants, rnd 1 and rnd 2 are random numbers in the range [0, 1], K is a constriction factor given by:

K = 2 2 b b 2 4 b , and b = b 1 + b 2 ; b > 4 E31

The performance of the PSO algorithm is significantly affected by the three factors w , c 1 , and c 2 . In this approach, a detection function defined as: φ k = gbest i x i k / pbest i x i k . The values of the three factors are adjusted dynamically using the following equations

α k = α initial α final 1 + e φ k k 1 + ln φ k L max / μ + α final E32
c 1 k = c 1 φ 1 k E33
c 2 k = c 2 φ k E34

where w initial and w final lie in the range 0 < w < 2 , L max is the final evolutionary generation, and k is the current evolutionary generation.

3.4. PSO-based Taylor-Kalman filter structure

The adaptive filter structure with the proposed adaptive algorithm is shown in the Figure 1 . This particular structure is modeled for only the fundamental component of the signal to be estimated. For a signal with Nth order harmonics, the same structure can be extended to meet the requirements. The signal is modeled using Taylor series up to the second order, so the filter has six weights for the six inputs that are used for the purpose of estimation. The performance of the algorithm is judged on the basis of speed of convergence, which is verified from the simulation results in Section 4.

Figure 1.

Filter structure of the proposed approach.


4. Simulation and results

The performance of the proposed algorithm for power system signals has been shown with the help of three computer simulated examples.

4.1. Tracking of a nonstationary signal with simultaneous change in amplitude, phase, and frequency

A nonstationary test signal as shown in Eq. (35) is generated in MATLAB. The simulation is done over 1000 samples of the signal. To make the signal nonstationary, a double step is introduced in the signal by changing the value of amplitude from 500 to 700 samples. This is done to simulate voltage surge occurrences in real time, where the amplitude increases from that of its desired value for some period of time. Similar disturbances also change the values of frequency and phase which is also simulated to test the tracking accuracy of the proposed algorithm. The results in Figures 2 5 reveal that the accuracy of the proposed algorithm is very high and tracking is achieved within one cycle of the signal.

y i = a i sin t dt + ϕ i + n i E35


a i = { 0.8 p . u , i < 500 samples 1 p . u , 500 < i < 700 samples 0.8 p . u , i > 700 samples

is the signal amplitude.

f i = { 50 Hz , i < 500 samples 51 Hz , i > 500 samples

Figure 2.

Estimated amplitude under 30 dB noise.

Figure 3.

Estimated frequency under 30 dB noise.

Figure 4.

Estimated phase under 30 dB noise.

Figure 5.

Estimated signal under 30 dB noise.

The sampling frequency f s = 2 kHz ϕ i = 0.5 rad 0.45 rad

n i is the noise sequence with power level 30 dB. This signal fed into the algorithm and simulated in MATLAB 2013a environment.

Table 1 contains the estimated values of the amplitude, frequency, and phase of the signal under different noise conditions. In this simulation, three different noise conditions have been considered. The simulation is carried out in a dynamic noise range from high noise (20 dB) to low noise (40 dB) conditions to test the performance of the proposed algorithm under noise. The analysis of the performance under noisy conditions shows that the proposed algorithm is able to track the desired signal very closely even under heavy noise conditions.

Parameter Samples Actual value Estimated value
LMS EKF Proposed method
20 dB 30 dB 40 dB 20 dB 30 dB 40 dB 20 dB 30 dB 40 dB
Amplitude 0–500 0.8 0.7 0.772 0.727 0.76 0.76 0.777 0.795 0.798 0.799
500–700 1 0.885 0.946 0.919 0.89 0.88 0.93 0.896 0.942 0.997
700–1000 0.8 0.7 0.772 0.727 0.78 0.79 0.787 0.789 0.794 0.797
Frequency 0–500 50 50.05 50.02 50.02 50.03 50.02 50.03 50.00 50.01 50.00
500–1000 51 50.38 50.55 50.63 50.43 50.76 51.02 50.94 50.99 51
Phase 0–500 0.5 0.53 0.513 0.521 0.53 0.53 0.53 0.511 0.513 0.533
500–1000 0.45 0.465 0.445 0.47 0.43 0.48 0.47 0.44 0.45 0.45

Table 1.

Estimated values of amplitude, frequency, and phase under different noise conditions.

4.2. Performance of the proposed algorithm in harmonic estimation

In this case, the ability of the proposed algorithm is tested with respect to the tracking of harmonics. The number of harmonic components present in the system is not constant, and it can vary from few to a large number. It is not possible for any method to track infinite number of harmonics but can handle a substantial quantity. In the real-time scenario, harmonics occur as odd multiples of the fundamental frequency, so the simulation is carried out with a system generated signal containing harmonics up to the 19th order.

y i = a 1 i sin i dt + ϕ 1 i + 0.8 i sin i 3 ω i dt + 0.4 + 0.6 i sin i 5 ω i dt + 0.3 + 0.5 i sin i 7 ω i dt + 0.25 + 0.4 i sin i 11 ω i dt + 0.2 + 0.2 i sin i 19 ω i dt + 0.1 E36

The signal parameters are taken as:

a i = 0.1 sin 2 πidt + 0.05 sin 10 πidt a 1 i = 1 + a i f i = 50 Hz f s i = 4 kHz ϕ 1 = 0.5 0.2 sin 2 π 5 idtd + 0.3

The amplitude, frequency, and phase are estimated, and results are shown in Figures 6 11 .

Figure 6.

Estimated fundamental amplitude for case 4.2.

Figure 7.

Estimated third amplitude for case 4.2.

Figure 8.

Estimated fifth amplitude for case 4.2.

Figure 9.

Estimated seventh amplitude for case 4.2.

Figure 10.

Estimated frequency of fundamental for case 4.2.

Figure 11.

Estimated signal for case 4.2.

Table 2 shows the comparison of the absolute errors in amplitude, frequency, and phase estimation for different harmonic components for EKF, LMS, RLS, and the proposed method. The values show that the higher order (>5th order) components exhibit higher error values for all the methods, but the comparison shows that among all the methods compared, the proposed method has the least values of error. This comparison sheds light on the superiority of the proposed method over the other methods.

Parameter Component Absolute error
EKF LMS RLS Proposed
Amplitude (harmonic order) A1 0.01 0.03 0.023 0.007
A3 0.013 0.032 0.04 0.003
A5 0.02 0.047 0.056 0.009
A7 0.04 0.058 0.0856 0.017
A11 0.025 0.0623 0.021 0.017
A19 0.03 0.0875 0.075 0.025
Frequency Fundamental 0.065 0.045 0.0201 0.058
Phase Φ1 0.0029 0.0087 0.0047 0.0023
Φ3 0.006 0.04 0.032 0.0005
Φ5 0.024 0.045 0.054 0.007
Φ7 0.005 0.076 0.072 0.002
Φ11 0.03 0.085 0.088 0.0019
Φ19 0.067 0.083 0.0765 0.04

Table 2.

Comparison of absolute errors.

4.3. Estimation of a power signal in the presence of DC component

When a fault occurs it not only distorts the signal by changing the voltage and current waveforms but some DC component that decays over time also gets added to the signal. DC components are nonperiodic in nature and this simulation shows that the proposed algorithm efficiently tracks nonperiodic components in the signal which is clearly evident from Figures 12 15 . A nonstationary test signal with a decaying DC component as shown in Eq.(37) is considered:

y i = a i sin t dt + ϕ i + n i E37


a i = A exp i / 300 p . u . A = 1 p . u .

is the signal amplitude. f i = 50 Hz

Figure 12.

Estimated amplitude for signal with decaying DC component.

Figure 13.

Estimated frequency for signal with decaying DC component.

Figure 14.

Estimated phase for signal with decaying DC component.

Figure 15.

Estimated signal with decaying DC component.

The sampling frequency f s = 2 kHz and ϕ i = 0.52 rad n i is the 30 dB noise.


5. Conclusion

The traditional Kalman filter has been extended to Taylor-Kalman filter which resulted in filters that are able to have flat magnitude and phase responses. These filters exhibit excellent tracking abilities and accurately estimate the amplitude, frequency and phase of a time varying power signal without any distortion. The further combination of the Taylor-Kalman filter with self-adaptive PSO makes the performance of the proposed method superior to the traditional approach. The methods can be individually used for the purpose of signal and parameter estimation, but individually, they suffer from some drawbacks. By combining the three methods into one hybrid method, the pitfalls of each are compensated by the other and hence much better results are obtained.

Further, the hardware implementation of the proposed method can be attempted for real-time applications [20, 21, 22, 23]. The hardware implementation of the proposed method can be embedded within an integrated circuit that will result in a system on chip that can be installed at power distribution centers, from where power gets distributed to the consumers, thus equipping them with a tool for detecting anomalies in power quality before power is dispatched to the utility network. The objective of developing such a technology is to create a compact and versatile tool. It is a small contribution toward the development of smart grid technology.


  1. 1. Lin HC, Lee CS. Enhanced FFT-based parameter algorithm for simultaneous multiple harmonics analysis. Proceedings of the Institute of Electrical and Electronics Engineers: Generation, Transmission and Distribution Proceedings. May 2001;148(3):209-214. DOI: 10.1049/ip-gtd:20010278
  2. 2. Chen CI, Chang GW, Hong RC, Li HM. Extended real model of Kalman filter for time-varying harmonics estimation. IEEE Transactions on Power Delivery. Jan 2010;25(1):17-26. DOI: 10.1109/TPWRD.2009.2035217
  3. 3. Nanda S, Hasan S, Pujari SS, Dash PK. A fast hybrid adaptive filter for parameter estimation of non-stationary sinusoid under noise. IEEE Conference PCTIC. 2015. DOI: 10.1109/PCITC.2015.7438146
  4. 4. Pradhan K, Routray A, Basak A. Power system frequency estimation using least mean square technique. IEEE Transactions on Power Delivery. Jul. 2005;20(3):1812-1816. DOI: 10.1109/TPWRD.2004.843453
  5. 5. Chen C-I, Chang GW. A two stage ADALINE for harmonic and inter harmonic measurement. IEEE Transactions on Industrial Electronics. 2009;56(6):2220-2228. DOI: 10.1109/TIE.2009.2017093
  6. 6. Pai PF. On-line tracking of instantaneous frequency and amplitude of dynamical system response. Mechanical Systems and Signal Processing. May 2010;24(4):1007-1024. DOI: org/10.1016/j.ymssp.2009.07.014
  7. 7. Sadinezhad I, Agelidis VG. Slow sampling online optimization approach to estimate power system frequency. IEEE Transactions on Smart Grid. June 2011;2(2). DOI: 10.1109/TSG.2011.2114374
  8. 8. Ren J, Kezunovic M. An improved fourier method for power system frequency estimation. North American Power Symposium (NAPS), 4–6 Aug, 2011, Boston, MA. pp. 1-6. DOI: 10.1109/NAPS.2011.6024845
  9. 9. Arghya S, Roy Choudhury S, Sengupta S. A self-synchronized ADALINE network for on-line tracking of power system harmonics. Elsevier Measurement. May 2011;44(4):784-790. DOI: 10.1016/j.measurement.2011.01.009
  10. 10. Ray PK, Subudhi B. Ensemble-Kalman-filter-based power system harmonic estimation. IEEE Transactions on Instrumentation and Measurement. December 2012;61(12). DOI: 10.1109/TIM.2012.2205515
  11. 11. Sadinezhad I, Agelidis VG. Real-time power system phasors and harmonics estimation using a new decoupled recursive-least-squares technique for DSP implementation. IEEE Transactions on Industrial Electronics. 2013;60(6):2295-2308. DOI: 10.1109/TIE.2012.2192895
  12. 12. Marchesan G, Cardoso Jr G, Mariotto L, Morais AP, Oliveira AL. An adaptive step-size least mean square algorithm for electric power systems frequency estimation in protective relays. 12th IET International Conference on IEEE Developments in Power System Protection (DPSP 2014). DOI: 10.1049/cp.2014.0118
  13. 13. Valtierra-Rodriguez M, Romero-Troncoso R d J, Osornio-Rios RA, Garcia-Perez A. Detection and classification of single and combined power quality disturbances using neural networks. IEEE Transactions on Industrial Electronics. May 2014;61(5). DOI: 10.1109/TIE.2013.2272276
  14. 14. Nanda S, Biswal M, Dash PK. Estimation of time varying signal parameters using an improved Adaline learning algorithm. AEU International Journal of Electronics and Communications. 2014;68(2):115-129. DOI: 10.1016/j.aeue.2013.07.014
  15. 15. Nanda S, Dash PK, Chakraborty T, Hasan S. A quadratic polynomial signal model and fuzzy adaptive filter for frequency and parameter estimation of nonstationary power signals. Measurement. June 2016;87:274-293. DOI: 10.1016/j.measurement.2016.03.026
  16. 16. He W, Guo S, Teng Z, Li F, Yang Y. Frequency estimation of distorted and noisy signals in power systems by FFT-based approach. IEEE Transactions on Power Systems. March 2014;29(2). DOI: 10.1109/TPWRS.2013.2283273
  17. 17. Lee DG, Kang SH, Nam SR. Modified dynamic Phasor estimation algorithm for the transient signals of distributed generators. IEEE Transactions on Smart Grid. March.2013;4(1):419-424. DOI: 10.1109/TSG.2012.2233772
  18. 18. de la O Serna JA, Maldonado JR. Taylor Kalman Fourier filters for instantaneous oscillating Phasor and harmonic estimates. IEEE Transactions on Instrumentation and Measurement. April 2012;61(4):941-951. DOI: 10.1109/TIM.2011.2178677
  19. 19. Kingshuk B, Nanda S. An adaptive filtering technique with self-adaptive PSO for estimation of non-stationary signals. 2016 International Conference on Communication and Signal Processing (ICCSP). IEEE, 2016. DOI: 10.1109/ICCSP.2016.7754311
  20. 20. Cardenas A, Guzman C, Agbossou K. Development of a FPGA based real-time power analysis and control for distributed generation Interface. IEEE Transactions on Power Apparatus and Systems. August 2012;27(3):1343-1353. DOI: 10.1109/TPWRS.2012.2186468
  21. 21. Jindapetch N, Chewae S, Phukpattaranont P. FPGA implementations of an ADALINE adaptive filter for power-line noise cancellation in surface electromyography signals. Measurement. 2012;45(3):405-414. DOI: 10.1016/j.measurement.2011.11.004
  22. 22. Sarita N, Dash PK. A gauss-Newton ADALINE for dynamic Phasor estimation of power signals and its FPGA implementation. IEEE Transactions on Instrumentation and Measurement, (Early access). 2016. DOI: 10.1109/TIM.2016.2620841
  23. 23. Sarita N, Dash PK. Field programmable gate array implementation of fuzzy variable step size adaptive linear element for adaptive frequency estimation. IET Signal Processing, (Available online. DOI: 10.1049/iet-spr.2016.0574

Written By

Sarita Nanda

Submitted: 22 April 2017 Reviewed: 23 October 2017 Published: 21 February 2018