<figure><img src="image_1.png" /></figure>

HAL
open science

# Adaptive Observer-Based Grid-Synchronization and Sequence Extraction Techniques for Renewable Energy Systems: A Comparative Analysis

Hafiz Ahmed, Mohamed Benbouzid

► To cite this version:

Hafiz Ahmed, Mohamed Benbouzid. Adaptive Observer-Based Grid-Synchronization and Sequence Extraction Techniques for Renewable Energy Systems: A Comparative Analysis. Applied Sciences, 2021, 11, 10.3390/app11020653. [hal-05335355](https://hal-univ-brest.fr/hal-05335355v1)

HAL Id: hal-05335355

https://hal.univ-brest.fr/hal-05335355v1

Submitted on 28 Oct 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Adaptive Observer-Based Grid-Synchronization and Sequence Extraction Techniques for Renewable Energy Systems: A Comparative Analysis †

Hafiz Ahmed <fnref n="1" />,* and Mohamed Benbouzid <fnref n="2" />,<fnref n="3" />

**Abstract:** Various faults and disturbances often make the grid unbalanced and distorted. To make the grid-connected converters resilient to unbalanced and distorted grid, real-time estimation of grid voltage sequences are essential. Adaptive-observer based grid-synchronization techniques are very suitable for this purpose. In this paper, a comparative analysis of three recently proposed adaptive observer-based grid-synchronization technique are presented for frequency estimation and sequence extraction. The considered techniques are: global adaptive observer, gain normalized adaptive observer, and second-order generalized integrator type adaptive observer. Following a fair and systematic gain tuning process, the performances are verified through extensive numerical simulation and experiments. Experimental results demonstrate the merits and demerits of the individual techniques. Academic researchers and industrial practitioners will find the results useful for the selection of adaptive observer-based grid-synchronization and sequence extraction techniques.

**Keywords:** grid-synchronization; frequency-locked loop; sequence extraction; grid-connected renewable energy systems

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Ahmed, H.; Benbouzid, M. Adaptive Observer-Based Grid-Synchronization and Sequence Extraction Techniques for Renewable Energy Systems: A Comparative Analysis. *Appl. Sci.* **2021**, *11*, 653. https://doi.org/10.3390/app11020653

Received: 9 November 2020

Accepted: 1 January 2021

Published: 11 January 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Renewable energy systems (RES) are the way forward towards achieving net-zero carbon emission target as mandated by many countries around the globe. This has contributed to a major push for ever increasing integration of RES into the existing electric power grid. This has attracted a lot of research attention into the control of grid-connected converters (GCC) [1–11]. Grid-connected converters can transfer maximum power to the grid only when the operation is in-phase with the grid voltage signal. Moreover, operation of GCC in unbalanced grid also require real-time extraction of grid voltage sequences. This motivated the researchers to devote significant attention to grid-synchronization based sequence extraction techniques [12–15].

Some of the recent and most popular grid synchronization techniques are: Kalman filter [16–18], demodulation [19,20], phase-locked loop (PLL) [21–26], adaptive observer [27–29], second-order generalized integrator (SOGI) [30–32], open-loop techniques [33], etc. Kalman filter [16] and least square-based techniques [34] typically work by considering a linear-parametric model of the grid voltage. Real-time matrix inversion are required for these techniques. This adds computational complexity for real-time low-cost embedded systems. Similar problems are to some extent true for the maximum likelihood estimator [18]. Demodulation [19] has several low-pass filters to tune, requires the use of trigonometric function, and can only be used for frequency estimation and not for sequence extraction.

Phase-locked loop (PLL) [21–26] has gained lot of attention for grid-synchronization application due to simplicity. Classical PLL generates double frequency oscillation in the

<footnote n="1">¹ School of Engineering and the Built Environment, Birmingham City University, Birmingham B4 7XG, UK</footnote>

<footnote n="2">² UMR CNRS 6027 IRDL, University of Brest, 29238 Brest, France; mohamed.benbouzid@univ-brest.fr</footnote>

<footnote n="3">³ Logistics Engineering College, Shanghai Maritime University, Shanghai 201306, China</footnote>

<footnote>* Correspondence: hafiz.h.ahmed@ieee.org</footnote>

<footnote>† This paper is an extended version of our paper published in the 46th Annual Conference of the IEEE Industrial Electronics Society (IECON-2020), Singapore.</footnote>
presence of unbalanced grid. To avoid this issue, several solutions are proposed in the literature that require additional filtering. Moreover, tuning the PLL in the presence of additional filtering stage can be complicated. PLL is a closed-loop system where the stability is dependent on various parameters. Open-loop grid-synchronization technique [33] has a significant advantage in this regard as unconditional stability can be guaranteed. However, this type of technique often require calculating the derivative. Real-time estimation of signal derivatives are often challenging when the base signal is noisy.

Many nonlinear techniques are presented in the literature in recent times, such as adaptive observer [27–29,35], adaptive notch filter (ANF) [4], self-tuning filter (STF) [5], second-order generalized integrator (SOGI) [30–32], etc. Unlike PLL, these techniques perform the task of grid-synchronization through the help of frequency-locked loop. ANF, SOGI, and STF are similar in nature. All these techniques have band-pass filter property. SOGI has a limited complex conjugate closed-loop poles tuning range [36]. Adaptive observers [27–29] do not have any such limitation. This made adaptive observers-based technique a suitable choice for grid-synchronization in the literature.

Shah and Singh, in Reference [28], proposed the application of adaptive observer for grid-synchronization and sequence extraction. Using Lyapunov function and LaSalle's invariance principle, this observer provides global asymptotic stability of the state and parameter estimation error dynamics. Unfortunately, the frequency update law used in this work do not use any gain normalization. This makes the frequency estimation dynamics dependent on the grid voltage amplitude. This can slow down the convergence of the unknown frequency estimation when the grid voltages suffer significant sag. This issue has been resolved in Reference [27] by introducing gain normalized frequency estimation law. Both [27,28] consider the grid voltage and its derivative as the state variables to obtain the dynamic model. SOGI-type adaptive observers as introduced in Reference [29,36] use the same dynamical model as used by SOGI. In this case, quadrature signals are used as the state variables. This type of observer can further enhance the dynamic tuning range of SOGI.

In this paper, the focus is on comparative analysis of these recently proposed adaptive observers. Such a comparative comparison is not yet available in the literature. To this end, the observers are used for the purpose of grid-synchronization based sequence extraction for unbalanced and distorted grid. The comparison is facilitated by a systematic and fair gain tuning method. Extensive numerical simulations with adaptive notch filter-based frequency-locked loop [37–39] and experimental results are provided to highlight the advantages of the selected technique over the literature. Comprehensive and systemic presentation and extensive performance comparison of three recently proposed adaptive observer-based grid-synchronization techniques are the main contribution of this paper.

The rest of this paper is organized as follows: Problem formulation is given in Section 2, overview of the selected techniques are given in Section 3, numerical simulation and experimental performance comparison are given in Section 4, and, finally, Section 5 concludes this paper.

## 2. Problem Formulation

Unbalanced three-phase grid voltages comprise of positive (+), negative (−) and zero (0) sequence components are given by:
$$ U_a = \underbrace{A^+ \sin(\omega t + \phi^+)}_{U_a^+} + \underbrace{A^- \sin(\omega t + \phi^-)}_{U_a^-} + \underbrace{A^0 \sin(\omega t + \phi^0)}_{U_a^0}, \quad (1) $$

$$ U_b = \underbrace{A^+ \sin\left(\omega t - \frac{2\pi}{3} + \phi^+\right)}_{U_b^+} + \underbrace{A^- \sin\left(\omega t + \frac{2\pi}{3} + \phi^-\right)}_{U_b^-} + \underbrace{A^0 \sin\left(\omega t + \phi^0\right)}_{U_b^0}, \quad (2) $$

$$ U_c = \underbrace{A^+ \sin\left(\omega t + \frac{2\pi}{3} + \phi^+\right)}_{U_c^+} + \underbrace{A^- \sin\left(\omega t - \frac{2\pi}{3} + \phi^-\right)}_{U_c^-} + \underbrace{A^0 \sin\left(\omega t + \phi^0\right)}_{U_c^0}, \quad (3) $$

where the amplitude is denoted by A, the initial phase angle is denoted by φ, and the angular frequency is denoted by ω with the nominal value being given by ωn = 100π. In this paper, the problem is to estimate the unknown angular frequency, ω and extract the positive U+ = [ Ua+ Ub+ Uc+ ]T, negative U− = [ Ua− Ub− Uc− ]T, and zero U0 = [ Ua0 Ub0 Uc0 ]T sequence components from the measured grid voltages (1)–(3). For this purpose, let us consider the 90° phase shifted version of the grid voltages (1)–(3):

$$ qU_a = \underbrace{A^+ \cos(\omega t + \phi^+)}_{qU_a^+} + \underbrace{A^- \cos(\omega t + \phi^-)}_{qU_a^-} + \underbrace{A^0 \cos(\omega t + \phi^0)}_{qU_a^0}, \quad (4) $$

$$ qU_b = \underbrace{A^+ \cos\left(\omega t - \frac{2\pi}{3} + \phi^+\right)}_{qU_b^+} + \underbrace{A^- \cos\left(\omega t + \frac{2\pi}{3} + \phi^-\right)}_{qU_b^-} + \underbrace{A^0 \cos\left(\omega t + \phi^0\right)}_{qU_b^0}, \quad (5) $$

$$ qU_c = \underbrace{A^+ \cos\left(\omega t + \frac{2\pi}{3} + \phi^+\right)}_{U_c^+} + \underbrace{A^- \cos\left(\omega t - \frac{2\pi}{3} + \phi^-\right)}_{U_c^-} + \underbrace{A^0 \cos\left(\omega t + \phi^0\right)}_{U_c^0}. \quad (6) $$

From the grid voltages (1)–(3) and their quadrature version (4)–(6), the positive, negative and zero sequence components can be extracted by the following equations [39]:

$$ \begin{bmatrix} U_a^+ \\ U_b^+ \\ U_c^+ \end{bmatrix} = \frac{1}{3} \begin{bmatrix} 1 & -\frac{1}{2} & -\frac{1}{2} \\ -\frac{1}{2} & 1 & -\frac{1}{2} \\ -\frac{1}{2} & -\frac{1}{2} & 1 \end{bmatrix} \begin{bmatrix} U_a \\ U_b \\ U_c \end{bmatrix} + \frac{1}{2\sqrt{3}} \begin{bmatrix} 0 & 1 & -1 \\ -1 & 0 & 1 \\ 1 & -1 & 0 \end{bmatrix} \begin{bmatrix} qU_a \\ qU_b \\ qU_c \end{bmatrix}, \quad (7) $$

$$ \begin{bmatrix} U_a^- \\ U_b^- \\ U_c^- \end{bmatrix} = \frac{1}{3} \begin{bmatrix} 1 & -\frac{1}{2} & -\frac{1}{2} \\ -\frac{1}{2} & 1 & -\frac{1}{2} \\ -\frac{1}{2} & -\frac{1}{2} & 1 \end{bmatrix} \begin{bmatrix} U_a \\ U_b \\ U_c \end{bmatrix} - \frac{1}{2\sqrt{3}} \begin{bmatrix} 0 & 1 & -1 \\ -1 & 0 & 1 \\ 1 & -1 & 0 \end{bmatrix} \begin{bmatrix} qU_a \\ qU_b \\ qU_c \end{bmatrix}, \quad (8) $$

$$ \begin{bmatrix} U_a^0 \\ U_b^0 \\ U_c^0 \end{bmatrix} = \left\{ \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix} - \frac{1}{3} \begin{bmatrix} 1 & -\frac{1}{2} & -\frac{1}{2} \\ -\frac{1}{2} & 1 & -\frac{1}{2} \\ -\frac{1}{2} & -\frac{1}{2} & 1 \end{bmatrix} \right\} \begin{bmatrix} U_a \\ U_b \\ U_c \end{bmatrix}. \quad (9) $$

From Equations (7)–(9), it is clear that to extract the positive, negative, and zero sequence components, quadrature signals are required. Adaptive observer-based grid-synchronization techniques will be used for this purpose. Details of the selected adaptive-observer based techniques are given in the following section.

### 3. Overview of the Selected Techniques

#### 3.1. Global Adaptive Observer (GAO)

This section briefly presents the globally asymptotically stable adaptive observer (GAO) as proposed in Reference [28,40].
### 3.1.1. Development of Phase Voltage Dynamical Model

Let us first consider the phase a. Models of phase b and c are the same as phase a. To develop the state-space model, let us first consider the phase voltage $U_a$ and its derivative as the following state variables:

$$ \xi_1 = U_a = A^+ \sin(\omega t + \phi^+) + A^- \sin(\omega t + \phi^-) + A^0 \sin(\omega t + \phi^0), $$

$$ \xi_2 = \dot{\xi}_1 = \dot{U}_a = \omega \{ A^+ \cos(\omega t + \phi^+) + A^- \cos(\omega t + \phi^-) + A^0 \cos(\omega t + \phi^0) \}. $$

By considering the state vector, $\xi = [ \xi_1 \ \xi_2 ]^T$, the state-space model of the phase voltage dynamics can be written as:

$$ \dot{\xi} = A\xi, \quad (10) $$

$$ U_a = C\xi, \quad (11) $$

where

$$ A = \begin{bmatrix} 0 & 1 \\ -\omega^2 & 0 \end{bmatrix}, C = \begin{bmatrix} 1 & 0 \end{bmatrix}. $$

State matrix $A$ contains the unknown grid frequency $\omega$. To facilitate the frequency adaptation law development for $\omega$, let us rewrite the frequency in terms of the known nominal frequency as, $\omega^2 = \eta\omega_n^2, \eta > 0$. Observability matrix of system (10) and (11) is given by:

$$ O = \begin{bmatrix} C \\ CA \end{bmatrix} = \begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix}. \quad (12) $$

Since the observability matrix $O$ is of rank 2 which is the same as the rank of the state matrix $A$, the system (10) and (11) is observable. Observer design for this system is described in the following section.

### 3.1.2. Observer Design

Let us consider the following transformation for the state vector $\xi$:

$$ \Xi = T\xi, \quad (13) $$

where $\Xi = [ \Xi_1 \ \Xi_2 ]^T$ is the transformed state vector and the non-singular matrix $T$ is given by:

$$ T = \frac{(1 + \eta)^{-1}}{\omega_n^2} \begin{bmatrix} 1 & -\frac{1}{\omega_n} \\ \eta\omega_n & 1 \end{bmatrix}. $$

For the state vector $\Xi$, the transformed state and output matrices are given by:

$$ A_{\Xi} = T A T^{-1} = A, C_{\Xi} = C T^{-1}. $$

Then, the transformed dynamical model for the phase voltage $U_a$ is given by:

$$ \dot{\Xi} = A_{\Xi}\Xi, \quad (14) $$

$$ U_a = C_{\Xi}\Xi. \quad (15) $$

The following Luenberger observer can be designed for the system (14) and (15):

$$ \dot{\hat{\Xi}} = \hat{\mathcal{A}}_{\Xi}\hat{\Xi} + \mathcal{L}(U_a - \hat{U}_a) \quad (16) $$

where $\hat{\mathcal{A}}$ represents estimated value, $\hat{U}_a = C_{\Xi}\hat{\Xi}, \mathcal{L} \in \mathbb{R}^{2 \times 1}$ is the observer gain matrix, and

$$ \hat{\mathcal{A}}_{\Xi} = \begin{bmatrix} 0 & 1 \\ -\hat{\eta}\omega_n^2 & 0 \end{bmatrix}. $$
To study the convergence of the observer, let us consider the state and parameter estimation errors as:

$$ 
\begin{aligned}
\tilde{\Xi} &= \Xi - \hat{\Xi}, \\
\tilde{\eta} &= \eta - \hat{\eta}.
\end{aligned}
 $$

Then, the closed-loop error dynamics for the system (14) and (15) and the observer (16) are given by:

$$ 
\begin{aligned}
\dot{\hat{\Xi}} &= \hat{\Xi} - \hat{\Xi}, \\
&= A_{\Xi\Xi} - \{\hat{A}_{\Xi\Xi}\hat{\Xi} + \mathcal{L}(C_{\Xi} - C_{\Xi}\hat{\Xi})\}, \\
&= A_{\Xi\Xi} - \{A_{\Xi\Xi}\hat{\Xi} - A_{\Xi\Xi}\hat{\Xi} + \hat{A}_{\Xi\Xi}\hat{\Xi} + \mathcal{L}(C_{\Xi\Xi} - C_{\Xi\Xi}\hat{\Xi})\}, \\
&= (A_{\Xi} - \mathcal{L}C_{\Xi})(\Xi - \hat{\Xi}) + (A_{\Xi} - \hat{A}_{\Xi})\hat{\Xi}, \\
&= (A_{\Xi} - \mathcal{L}C_{\Xi})\tilde{\Xi} + \mathcal{R}\omega_n^2\tilde{\mu}\hat{\Xi}_1,
\end{aligned}
 \quad (17) $$

where $\mathcal{R} = [\begin{matrix} 0 & 1 \end{matrix}]^T$. For suitably designed observer gain matrix $\mathcal{L} \in \mathbb{R}^{2 \times 1}$, the matrix, $A_{\Xi} - \mathcal{L}C_{\Xi}$ can be made Hurwitz, i.e., all the eigenvalues have strictly negative real parts. This makes the closed-loop error system strictly positive real (SPR) [41] and according to the SPR lemma, there exists a matrix $\mathcal{M} = \mathcal{M}^T > 0$, such that

$$ (A_{\Xi} - \mathcal{L}C_{\Xi})^T \mathcal{M} + M(A_{\Xi} - \mathcal{L}C_{\Xi}) = -Q, \quad Q = Q^T > 0, \quad (18) $$

$$ \mathcal{M}\mathcal{R} = C_{\Xi}^T. \quad (19) $$

The following Lyapunov function is considered for system (17):

$$ V(\tilde{\Xi}, \tilde{\eta}) = \tilde{\Xi}^T \mathcal{M} \tilde{\Xi} + \gamma^{-1} \tilde{\eta}^2, \quad \gamma > 0. \quad (20) $$

Then, the derivative of the Lyapunov function becomes:

$$ \dot{V} = -\tilde{\Xi}^T \tilde{\Xi} + 2(U_a - \hat{U}_a)\tilde{\eta}\omega_n^2\hat{\Xi}_1 - 2\gamma^{-1}\tilde{\eta}\dot{\hat{\eta}}. \quad (21) $$

If the frequency update law is chosen as:

$$ \dot{\tilde{\eta}} = \gamma\omega_n^2(U_a - \hat{U}_a)\hat{\Xi}_1. \quad (22) $$

Then, it becomes that $\dot{V} \le 0$, which proves the Lyapunov stability of the estimation error dynamics. From the estimated transformed state vector $\hat{\Xi}$, the state vector $\tilde{\xi}$ can be found through the following equation:

$$ \hat{\xi} = \mathcal{T}^{-1}\hat{\Xi}, \quad (23) $$

where

$$ \mathcal{T}^{-1} = \omega_n^2 \begin{bmatrix} 1 & \frac{1}{\omega_n} \\ -\hat{\eta}\omega_n & 1 \end{bmatrix}. $$

The process described above can be applied in the same way to the grid phase voltages $U_b$ and $U_c$. As the frequency is same for all phases, only one frequency identification law is enough. From the signals $U_a$, $U_b$, and $U_c$ and their derivatives (i.e., quadrature signal) obtained through the observer (16), the sequences can be easily estimated through the Equations (7)–(9).

## 3.2. Gain Normalized Adaptive Observer

This section briefly presents the gain normalized adaptive observer (GNAO) as published in Reference [27]. This observer uses the same dynamic model as presented in Section 3.1.1. For further development, let us define the unknown frequency $\omega$ as,
$\omega = \omega_n + \Delta\omega$, where $\Delta\omega$ is the frequency deviation. Let us consider the following state transformation:

$$ \Xi = \mathcal{T}\xi, \qquad (24) $$

where

$$ \mathcal{T} = \frac{1}{2\omega^3} \begin{bmatrix} \omega & -1 \\ \omega^2 & \omega \end{bmatrix}. $$

Then, the transformed state and output matrices are given by:

$$ A_{\Xi} = \mathcal{T}A\mathcal{T}^{-1} = A, C_{\Xi} = C\mathcal{T}^{-1} = \begin{bmatrix} \omega^2 & \omega \end{bmatrix}. $$

Then, the following Luenberger type linear observer is designed:

$$ \hat{\Xi} = \hat{A}_{\Xi}\hat{\Xi} + \mathcal{L}(U_a - \hat{U}_a), \qquad (25) $$

where $\mathcal{L} = \begin{bmatrix} \mathcal{L}_1 & \mathcal{L}_2 \end{bmatrix}^T$ is the gain matrix for the observer and

$$ \hat{A}_{\Xi} = \begin{bmatrix} 0 & 1 \\ -(\omega_n + \Delta\hat{\omega})^2 & 0 \end{bmatrix}. \qquad (26) $$

The frequency deviation estimation law is given by:

$$ \hat{\Delta\hat{\omega}} = - \frac{\gamma(\mathcal{L}_1 + \mathcal{L}_2)\hat{\omega}^3 \hat{\Xi}_1 (U_a - \hat{U}_a)}{\sqrt{(\hat{\Xi}_1 2\hat{\omega}^3)^2 + (\hat{\Xi}_2 2\hat{\omega}^2)^2}}, \gamma > 0. \qquad (27) $$

The denominator term in the right hand side of Equation (27) scales the frequency deviation update law by the estimate of the grid voltage signal amplitude. This makes the adaptation law less sensitive to large voltage sag/swell. This can be very suitable to provide low voltage ride through capability if the observer is used inside the control module of any grid-connected converter.

Lyapunov function-based stability analysis as presented in Section 3.1.2 cannot be applied to the observer (25) due to the presence of gain normalization term in Equation (27). As such, only local stability results are presented in Reference [27]. From the estimated state vector $\hat{\Xi}$, the state vector $\hat{\xi}$ can be calculated through the following equation:

$$ \hat{\Xi} = \mathcal{T}^{-1}\hat{\xi}, \qquad (28) $$

where

$$ \mathcal{T}^{-1} = \begin{bmatrix} \hat{\omega}^2 & \hat{\omega} \\ -\hat{\omega}^3 & \hat{\omega}^2 \end{bmatrix}. $$

The same observer (25) can be applied to grid phase voltages $U_b$ and $U_c$. From the signals $U_a$, $U_b$, and $U_c$ and their derivatives (i.e., quadrature signal) obtained through the observer (25), the sequences can be easily estimated through the Equations (7)–(9).

### 3.3. SOGI-Type Adaptive Observer

This Section briefly presents the SOGI-type adaptive observer (SAO) result proposed in Reference [29]. Unlike the previous two observers, this observer uses a different dynamic model.

#### 3.3.1. Quadrature Signal-Based Dynamic Model Development

For model development purposes, let us consider that the states are:
$$ \zeta_1 = U_a^{\perp} = -A^+ \cos(\omega t + \phi^+) - A^- \cos(\omega t + \phi^-) - A^0 \cos(\omega t + \phi^0), $$

$$ \zeta_2 = U_a = A^+ \sin(\omega t + \phi^+) + A^- \sin(\omega t + \phi^-) + A^0 \sin(\omega t + \phi^0). $$

Then, the following can be obtained:

$$ \dot{\zeta}_1 = \zeta_2 \omega, $$

$$ \dot{\zeta}_2 = -\zeta_1 \omega. $$

By considering the state vector, $\zeta = [ \zeta_1 \ \zeta_2 ]^T$, the dynamic model for the grid phase voltage $U_a$ is given by:

$$ \dot{\zeta} = A\zeta, \qquad (29) $$

$$ U_a = C\zeta, \qquad (30) $$

where

$$ A = \begin{bmatrix} 0 & \omega \\ -\omega & 0 \end{bmatrix}, C = \begin{bmatrix} 0 & 1 \end{bmatrix}. $$

Observability matrix of system (29) and (30) is given by:

$$ O = \begin{bmatrix} C \\ CA \end{bmatrix} = \begin{bmatrix} 0 & 1 \\ -\omega & 0 \end{bmatrix}. \qquad (31) $$

Since the observability matrix $O$ is of rank 2 which is the same as the rank of the state matrix $A$, the system (29) and (30) is observable. Observer design for this system is described in the following section.

### 3.3.2. Observer Design

Let us consider the following transformation:

$$ \Xi = T\zeta, \qquad (32) $$

where

$$ T = \frac{1}{2\omega} \begin{bmatrix} 1 & 1 \\ -1 & 1 \end{bmatrix}. $$

Then, the transformed state and output matrices are given by:

$$ A_{\Xi} = T A T^{-1} = A, C_{\Xi} = C T^{-1} = \begin{bmatrix} \omega & \omega \end{bmatrix}. $$

Then, the transformed dynamic model for the phase voltage $U_a$ is given by:

$$ \dot{\Xi} = A_{\Xi\Xi}, \qquad (33) $$

$$ U_a = C_{\Xi\Xi}. \qquad (34) $$

For model (33) and (34), the following Luenberger type linear observer is designed:

$$ \dot{\hat{\Xi}} = \hat{A}_{\Xi\Xi} + L(U_a - \hat{U}_a), \qquad (35) $$

where $L = [ \mathcal{L}_1 \ \mathcal{L}_2 ]^T$ is the gain matrix for the observer and

$$ \hat{A}_{\Xi} = \begin{bmatrix} 0 & \omega_n + \Delta\hat{\omega} \\ -(\omega_n + \Delta\hat{\omega}) & 0 \end{bmatrix}. \qquad (36) $$
Moreover, the gain normalized frequency deviation estimation law is given by:

$$ \Delta\dot{\omega} = - \frac{\gamma(\mathcal{L}_1 + \mathcal{L}_2)\hat{\omega}\hat{\Xi}_1(U_a - \hat{U}_a)}{\hat{\Xi}^T\hat{\Xi}}, \quad \gamma > 0. \quad (37) $$

Similar to the gain normalized adaptive observer presented in Section 3.2, only local stability results are presented in Reference [29] due to the presence of the gain normalization term in Equation (37). From the estimated state vector $\hat{\Xi}$, the state vector $\hat{\xi}$ can be calculated through the following equation:

$$ \hat{\xi} = \mathcal{T}^{-1}\hat{\Xi}, \quad (38) $$

where

$$ \mathcal{T}^{-1} = \hat{\omega} \begin{bmatrix} 1 & -1 \\ 1 & 1 \end{bmatrix}. $$

The same observer (35) can be applied to grid phase voltages $U_b$ and $U_c$. From the signals $U_a$, $U_b$, and $U_c$ and their derivatives (i.e., quadrature signal) obtained through the observer (35), the sequences can be easily estimated through the Equations (7)–(9).

## 4. Results and Discussions

In this section, comparative performance verification is considered. The selected observers are tuned by considering the closed-loop poles as: $-1.5\omega_n \pm i\omega_n$. This corresponds to $\mathcal{L} = [0.0012 \ 2.6250]^T$ (GAO and GNAO) and $\mathcal{L} = [0.3750 \ 2.6250]^T$. A 2-cycle convergence criterion is used to tune the parameter $\gamma$ through trial and error. This corresponds to the value of $\gamma$ as 1000, 150 and 0.2 for GAO, GNAO, and SAO, respectively. With these gains, the selected techniques have the same dynamic response in the presence of frequency step change in the ideal grid voltage situation. FLL-based adaptive notch filter (ANF-FLL) [37–39] is selected as the comparison technique. The ANF gain is selected as $1/\sqrt{2}$, while the FLL gain is selected as 18,000, as suggested in Reference [39].

### 4.1. Simulation Results

The following test cases are considered in this section:

* Test-I: Balanced to unbalanced grid,

* Test-II: Voltage sag,

* Test-III: Phase step change, and

* Test-IV: DC offset.

Figure 1 shows the comparative simulation for Test-I. In this test, the negative and zero sequence voltage had a sudden step change of +0.1 p.u. and +0.05 p.u., while the positive sequence voltage had a step change of -0.2 p.u. The frequency estimated by the adaptive observer-based techniques converged in 2-cycles, while ANF took 2.5 cycles with significantly higher overshoot. In addition, the positive, negative, and zero sequence amplitudes are estimated within ≈0.5 cycle by the adaptive observer-based techniques, while it is more than ≈1 cycle for ANF.

Voltage sag of -0.5 p.u. in balanced grid is considered in Test-II. Simulation results are given in Figure 2. As shown in Figure 2c, SAO and GNAO converged within 2 cycles, and GAO took 3 cycles, while, for ANF, it is almost 5 cycles. This shows that adaptive observer-based techniques are suitable for fast and accurate estimation of grid frequency in the presence of large voltage sag. In case of positive sequence amplitude estimation, half cycle convergence can be observed for the adaptive observer-based techniques, while it is ≈1 cycle for the ANF-FLL.

A sudden jump of -45° phase angle in a balanced grid is considered in Test-III. Simulation results as shown in Figure 3 demonstrate that all the techniques are heavily affected by the sudden change in phase angle. However, adaptive observer-based techniques converged very fast within roughly 45 milliseconds, while ANF took more than 60 milliseco-nds. In addition, the peak overshoot for ANF was also ≈+1.5 Hz higher than the adaptive observer-based techniques. This demonstrates the advantage of the adaptive
observer-based techniques for challenging scenarios, such as sudden large change in phase angle.

Due to current transformer saturation, signal conditioning problem, etc., DC offset can often corrupt the measured voltages. This is considered in Test-IV. The considered DC offset values are +0.05 p.u., +0.1 p.u., and -0.05 p.u. in phase a, b, and c, respectively. Results in Figure 4b show that ANF is very sensitive to DC offset. The frequency estimation ripple is between 48.2 Hz and 50.8 Hz for ANF. Among the adaptive observer-based techniques, GAO is the most sensitive with estimation ripple between 49 Hz and 50.78 Hz. GNAO and SAO have similar estimation ripple and the amplitude is lower than ANF and GAO.

By considering the test results as shown in Figures 1–4, it can be claimed that adaptive observer-based techniques have performed better than ANF in almost every situations. As such ANF will not be considered for experimental validation test in Section 4.2. Among adaptive observer-based techniques, GAO under performed compared to GNAO and SAO. It is to be noted here that all the adaptive observer-based techniques have similar dynamic performance.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Numerical simulation results for Test-I.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>(a) Simulated three-phase grid voltage signal</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>(b) Extracted positive sequence amplitudes</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>(c) Estimated frequencies</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>(a) Simulated three-phase grid voltage signal</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>(b) Estimated frequencies</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>(a) Simulated three-phase grid voltage signal</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>(b) Estimated frequencies</figcaption></figure>

<caption>**Figure 4.** Numerical simulation results for Test-IV.</caption>

## 4.2. *Experimental Results*

For the experimental validation, dSPACE MicroLabBox is used to implement the adaptive observer-based techniques in real-time. An overview of the experimental setup is given in Figure 5.

<figure><img src="image_13.png" /><figcaption>**Figure 5.** Experimental setup overview.</figcaption></figure>
### 4.2.1. Unbalanced Grid

In this test, initially the grid is balanced, i.e., no negative and zero sequence components. Suddenly, the positive sequence amplitude had a step change of -0.2 p.u. Moreover, negative sequence and positive sequence voltages of 0.2 p.u. and 0.1 p.u. are added into the grid. In addition, the phase angle of the positive, negative, and zero sequence components had a step change of -22.5°, -22.5°, and 22.5°, respectively. Experimental results for this test are given in Figure 6. From the estimated frequencies, it can be seen that the three techniques have roughly similar dynamic performance and converged in ≈2 cycles. This validates the gain tuning. GAO converged faster than the others, however, with the highest peak overshoot. Similarly, SAO had the lowest peak overshoot and the same time the slowest convergence time among the selected techniques. The sequence extraction performance of all the techniques are similar. The techniques extracted the positive, negative, and zero sequence components from the measured three-phase grid voltage in a fast and accurate manner.

### 4.2.2. Unbalanced and Distorted Grid

Grid distortion is often inevitable. So, this test considers unbalance and distortion at the same time. In this test, initially balanced grid is considered with 1.0 p.u. positive sequence component. After the fault, negative sequence of 0.1 p.u. and zero sequence of 0.05 p.u. are added. In addition, 5th and 13th-order positive sequence harmonics of 0.015 p.u. and 0.012 p.u. and 7th-order negative sequence harmonics of 0.017 p.u. are added. So, the negative sequence component has a distortion of 17% which is a highly challenging test scenario. Experimental results are given in Figure 7. Due to the presence of distortion in the voltage, the estimated frequencies also have estimation ripple. Similar to the previous test, GAO has the highest peak overshoot (≈0.5 Hz) and steady-state estimation ripple (≈0.3 Hz). SAO demonstrated the lowest peak overshoot (≈0.3 Hz) and steady-state estimation ripple (≈0.1 Hz).

Based on simulation and experimental results, a comparative summary of the adaptive observer-based techniques are presented in Table 1.

<table><caption>**Table 1.** Comparative analysis of the adaptive observer-based techniques.</caption><thead><tr><td>Characteristics</td><td>GAO</td><td>GNAO</td><td>SAO</td></tr></thead><tbody><tr><td>Dynamic Performance</td><td>Very Fast</td><td>Very Fast</td><td>Very Fast</td></tr><tr><td>Voltage Sag Sensitivity</td><td>High</td><td>Low</td><td>Low</td></tr><tr><td>Harmonic Disturbance</td><td>Good</td><td>Very Good</td><td>Very Good</td></tr><tr><td>DC Offset Rejection</td><td>Good</td><td>Very Good</td><td>Very Good</td></tr><tr><td>Computational Complexity</td><td>Medium</td><td>High</td><td>Medium</td></tr></tbody></table>
<figure><img src="image_14.png" /><figcaption>(a) Experimental three-phase grid voltage signal</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>(c) Extracted positive sequence by GAO</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>(e) Extracted positive sequence by SAO</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>(g) Extracted negative sequence by GNAO</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>(b) Estimated frequencies</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>(d) Extracted positive sequence by GNAO</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>(f) Extracted negative sequence by GAO</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>(h) Extracted negative sequence by SAO</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>(i) Extracted zero sequence component by the selected techniques</figcaption></figure>

<caption>**Figure 6.** Experimental results for balanced to unbalanced grid.</caption>
<figure><img src="image_23.png" /><figcaption>(a) Experimental three-phase grid voltage signal</figcaption></figure>

<figure><img src="image_24.png" /><figcaption>(c) Extracted positive sequence by GAO</figcaption></figure>

<figure><img src="image_25.png" /><figcaption>(e) Extracted positive sequence by SAO</figcaption></figure>

<figure><img src="image_26.png" /><figcaption>(g) Extracted negative sequence by GNAO</figcaption></figure>

<figure><img src="image_27.png" /><figcaption>(b) Estimated frequencies</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>(d) Extracted positive sequence by GNAO</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>(f) Extracted negative sequence by GAO</figcaption></figure>

<figure><img src="image_30.png" /><figcaption>(h) Extracted negative sequence by SAO</figcaption></figure>

<figure><img src="image_31.png" /><figcaption>(i) Extracted zero sequence component by the selected techniques</figcaption></figure>

<caption>**Figure 7.** Experimental results for unbalanced and distorted grid.</caption>
## 5. Conclusions

Three recently proposed adaptive observer-based grid-synchronization techniques are presented in this paper. These techniques are used to extract the positive, negative, and zero sequence components of unbalanced and distorted three-phase grid voltage signal. A systematic gain tuning procedure has been followed to ensure fair control parameters selection for the selected techniques. Using challenging test cases, the performance of the selected techniques are assessed experimentally. Comparative experimental results revealed the merits and demerits of individual techniques. Comparative experimental results as presented in this paper have the potential to give some guidance on the selection of the adaptive observer-based sequence extraction techniques for power and energy applications.

**Author Contributions:** Conceptualization, H.A. and M.B.; methodology, H.A.; software, H.A.; validation, H.A.; formal analysis, H.A. and M.B.; writing—original draft preparation, H.A.; writing—review and editing, M.B.; visualization, H.A. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data of this study are available from the corresponding author upon reasonable request.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Ouchen, S.; Betka, A.; Abdeddaim, S.; Menadi, A. Fuzzy-predictive direct power control implementation of a grid connected photovoltaic system, associated with an active power filter. *Energy Convers. Manag.* **2016**, *122*, 515–525. [CrossRef]

2. Safa, A.; Berkouk, E.M.; Messlem, Y.; Gouichiche, A. A robust control algorithm for a multifunctional grid tied inverter to enhance the power quality of a microgrid under unbalanced conditions. *Int. J. Electr. Power Energy Syst.* **2018**, *100*, 253–264. [CrossRef]

3. Ngo, T.; Biricik, S.; Basu, M. A Self-tuning Grid Synchronization Method for Active Power Filters. *Electr. Power Components Syst.* **2016**, *44*, 1947–1957. [CrossRef]

4. Biricik, S.; Komurcugil, H.; Tuyen, N.D.; Basu, M. Protection of Sensitive Loads Using Sliding Mode Controlled Three-Phase DVR with Adaptive Notch Filter. *IEEE Trans. Ind. Electron.* **2018**, *66*, 5465–5475. [CrossRef]

5. Biricik, S.; Redif, S.; Özerdem, Ö.C.; Khadem, S.K.; Basu, M. Real-time control of shunt active power filter under distorted grid voltage and unbalanced load condition using self-tuning filter. *IET Power Electron.* **2014**, *7*, 1895–1905. [CrossRef]

6. Biricik, S.; Komurcugil, H.; Ahmed, H.; Babaei, E. Super Twisting Sliding Mode Control of DVR with Frequency-Adaptive Brockett Oscillator. *IEEE Trans. Ind. Electron.* **2020**. [CrossRef]

7. Biricik, S.; Komurcugil, H. Optimized sliding mode control to maximize existence region for single-phase dynamic voltage restorers. *IEEE Trans. Ind. Inform.* **2016**, *12*, 1486–1497. [CrossRef]

8. Merai, M.; Naouar, M.W.; Slama-Belkhodja, I.; Monmasson, E. An adaptive PI controller design for DC-link voltage control of single-phase grid-connected converters. *IEEE Trans. Ind. Electron.* **2019**, *66*, 6241–6249. [CrossRef]

9. Kukrer, O.; Bayhan, S.; Komurcugil, H. Model-Based Current Control Strategy with Virtual Time Constant for Improved Dynamic Response of Three-Phase Grid-Connected VSI. *IEEE Trans. Ind. Electron.* **2019**, *66*, 4156–4165. [CrossRef]

10. Altin, N.; Ozdemir, S.; Komurcugil, H.; Sefa, I. Sliding-mode control in natural frame with reduced number of sensors for three-phase grid-tied LCL-interfaced inverters. *IEEE Trans. Ind. Electron.* **2018**, *66*, 2903–2913. [CrossRef]

11. Rahman, M.M.; Hossain, S.; Rezwan, S. Dispatched Power and Unbalanced Operation of a Three-Phase NPC Based Double-Stage Grid Connected PV System. In Proceedings of the 2019 4th International Conference on Electrical Information and Communication Technology (EICT), Khulna, Bangladesh, 20–22 December 2019; IEEE: Piscataway, NJ, USA, 2019; pp. 1–6.

12. Chatterjee, A.; Mohanty, K.B. Current control strategies for single phase grid integrated inverters for photovoltaic applications—A review. *Renew. Sustain. Energy Rev.* **2018**, *92*, 554–569. [CrossRef]

13. Guerrero-Rodríguez, N.; Rey-Boué, A.B.; Bueno, E.; Ortiz, O.; Reyes-Archundia, E. Synchronization algorithms for grid-connected renewable systems: Overview, tests and comparative analysis. *Renew. Sustain. Energy Rev.* **2017**, *75*, 629–643. [CrossRef]

14. Jaalam, N.; Rahim, N.; Bakar, A.; Tan, C.; Haidar, A.M. A comprehensive review of synchronization methods for grid-connected converters of renewable energy source. *Renew. Sustain. Energy Rev.* **2016**, *59*, 1471–1481. [CrossRef]

15. Jain, B.; Jain, S.; Nema, R. Control strategies of grid interfaced wind energy conversion system: An overview. *Renew. Sustain. Energy Rev.* **2015**, *47*, 983–996. [CrossRef]
16. Amirat, Y.; Oubrahim, Z.; Ahmed, H.; Benbouzid, M.; Wang, T. Phasor Estimation for Grid Power Monitoring: Least Square vs. Linear Kalman Filter. *Energies* **2020**, *13*, 2456. [CrossRef]

17. Ahmed, H.; Biricik, S.; Benbouzid, M. Linear Kalman Filter-Based Grid Synchronization Technique: An Alternative Implementation. *IEEE Trans. Ind. Inform.* **2020**. [CrossRef]

18. Choqueuse, V.; Belouchrani, A.; Auger, F.; Benbouzid, M. Frequency and phasor estimations in three-phase systems: Maximum likelihood algorithms and theoretical performance. *IEEE Trans. Smart Grid* **2018**, *10*, 3248–3258. [CrossRef]

19. Mathuria, K.; Hussain, I.; Singh, B.; Kumar, N. A quadrature oscillator-based DT for accurate estimation of fundamental load current for PV system in distribution network. *IEEE Trans. Ind. Inform.* **2018**, *15*, 3324–3333. [CrossRef]

20. Ahmed, H.; Biricik, S.; Benbouzid, M. Enhanced Frequency Adaptive Demodulation Technique For Grid-Connected Converters. *IEEE Trans. Ind. Electron.* **2020**. [CrossRef]

21. Chedjara, Z.; Massoum, A.; Massoum, S.; Wira, P.; Safa, A.; Gouichiche, A. A novel robust PLL algorithm applied to the control of a shunt active power filter using a self tuning filter concept. In *Proceedings of the 2018 IEEE International Conference on Industrial Technology (ICIT)*; IEEE: Piscataway, NJ, USA, 2018; pp. 1124–1131.

22. Bechouche, A.; Sediki, H.; Abdeslam, D.O.; Haddad, S. An adaptive neural PLL for grid synchronization. In *Proceedings of the IECON 2012-38th Annual Conference on IEEE Industrial Electronics Society*; IEEE: Piscataway, NJ, USA, 2012; pp. 4451–4456.

23. Bobrowska-Rafal, M.; Rafal, K.; Jasinski, M.; Kazmierkowski, M. Grid synchronization and symmetrical components extraction with PLL algorithm for grid connected power electronic converters—A review. *Bull. Pol. Acad. Sci. Tech. Sci.* **2011**, *59*, 485–497. [CrossRef]

24. Meral, M.E.; Çelik, D. Benchmarking simulation and theory of various PLLs produce orthogonal signals under abnormal electric grid conditions. *Electr. Eng.* **2018**, *100*, 1805–1817. [CrossRef]

25. Meral, M.E.; Çelík, D. A comprehensive survey on control strategies of distributed generation power systems under normal and abnormal conditions. *Annu. Rev. Control* **2019**, *47*, 112–132. [CrossRef]

26. Meral, M.E.; Celik, D. DSOGI-PLL based power control method to mitigate control errors under disturbances of grid connected hybrid renewable power systems. *Adv. Electr. Electron. Eng.* **2018**, *16*, 81–91. [CrossRef]

27. Ahmed, H.; Pay, M.L.; Benbouzid, M.; Amirat, Y.; Elbouchikhi, E. Gain normalized adaptive observer for three-phase system. *Int. J. Electr. Power Energy Syst.* **2020**, *118*, 105821. [CrossRef]

28. Shah, P.; Singh, B. Adaptive Observer Based Control for Roof-top Solar PV System. *IEEE Trans. Power Electron.* **2020**, *35*, 9402–9415. [CrossRef]

29. Ahmed, H.; Pay, M.L.; Benbouzid, M.; Amirat, Y.; Elbouchikhi, E. Hybrid estimator-based harmonic robust grid synchronization technique. *Electr. Power Syst. Res.* **2019**, *177*, 106013. [CrossRef]

30. Chishti, F.; Murshid, S.; Singh, B. Weak Grid Intertie WEGS With Hybrid Generalized Integrator for Power Quality Improvement. *IEEE Trans. Ind. Electron.* **2020**, *67*, 1113–1123. [CrossRef]

31. Bendib, A.; Chouder, A.; Kara, K.; Kherbachi, A.; Barkat, S.; Issa, W. New modeling approach of secondary control layer for autonomous single-phase microgrids. *J. Frankl. Inst.* **2019**, *356*, 6842–6874. [CrossRef]

32. Kherbachi, A.; Chouder, A.; Bendib, A.; Kara, K.; Barkat, S. Enhanced structure of second-order generalized integrator frequency-locked loop suitable for DC-offset rejection in single-phase systems. *Electr. Power Syst. Res.* **2019**, *170*, 348–357. [CrossRef]

33. Safa, A.; Berkouk, E.M.; Messlem, Y.; Chedjara, Z.; Gouichiche, A. A Pseudo Open Loop Synchronization technique for heavily distorted grid voltage. *Electr. Power Syst. Res.* **2018**, *158*, 136–146. [CrossRef]

34. Terriche, Y.; Guerrero, J.M.; Vasquez, J.C. Performance improvement of shunt active power filter based on non-linear least-square approach. *Electr. Power Syst. Res.* **2018**, *160*, 44–55. [CrossRef]

35. Hackl, C.M.; Landerer, M. A Unified Method for Online Detection of Phase Variables and Symmetrical Components of Unbalanced Three-Phase Systems with Harmonic Distortion. *Energies* **2019**, *12*, 3243. [CrossRef]

36. Ahmed, H.; Benbouzid, M. On the Enhancement of Generalized Integrator-based Adaptive Filter Dynamic Tuning Range. *IEEE Trans. Instrum. Meas.* **2020**, *69*, 7449–7457. [CrossRef]

37. Chilipi, R.S.R.; Al Sayari, N.; Al Hosani, K.H.; Beig, A.R. Adaptive Notch Filter-Based Multipurpose Control Scheme for Grid-Interfaced Three-Phase Four-Wire DG Inverter. *IEEE Trans. Ind. Appl.* **2017**, *53*, 4015–4027. [CrossRef]

38. Arya, S.R.; Maurya, R.; Naidu, T.A. Amplitude adaptive notch filter with optimized PI gains for mitigation of voltage based power quality problems. *CPSS Trans. Power Electron. Appl.* **2018**, *3*, 313–323. [CrossRef]

39. Yazdani, D.; Mojiri, M.; Bakhshai, A.; Joos, G. A fast and accurate synchronization technique for extraction of symmetrical components. *IEEE Trans. Power Electron.* **2009**, *24*, 674–684. [CrossRef]

40. Ahmed, H.; Benbouzid, M.; Ahsan, M.; Albarbar, A.; Shahjalal, M. Frequency Adaptive Parameter Estimation of Unbalanced and Distorted Power Grid. *IEEE Access* **2020**, *8*, 8512–8519. [CrossRef]

41. Khalil, H.K. *Nonlinear Control*; Pearson Education North Asia Ltd.: Hong Kong, 2014.