Learning-Driven Intelligent Passivity Control Using Nonlinear State Observers for Induction Motors

1. Introduction

The control of induction motors (IMs) continues to attract strong research interest due to their nonlinear structure, sensitivity to parameter variations, and widespread use in industrial applications. Traditional control methods often rely on accurate measurements of rotor speed and flux, which may not always be available or reliable in sensorless configurations. In this context, passivity-based control has emerged as a promising framework due to its solid theoretical foundation and robustness to system uncertainties. Early works demonstrated the use of passivity theory to ensure stable torque and flux regulation in IMs. For example, Cecati [1] proposed a passivity-based position controller, while Duarte-Mermoud [2] introduced an adaptive passivity approach for robust speed control. Travieso-Torres [3] extended this framework with adaptive structures for enhanced robustness under parameter variations. Moreover, Mansouri [4] applied passivity theory to sensorless IM control by integrating it with nonlinear observers for improved state reconstruction. In parallel, adaptive control strategies have played a pivotal role in addressing variations in machine parameters, such as rotor resistance and load torque, which frequently degrade control performance. Marino [5] introduced an adaptive input-output linearizing controller for IMs, achieving asymptotic tracking despite parametric uncertainties. Similarly, Chenafa [6] and Espinosa-Pérez [7] developed globally stable adaptive schemes using nonlinear feedback structures. In the sensorless domain, Belkacem Bekhiti [8], [9] proposed adaptive Luenberger and extended Kalman filter observers to enhance rotor speed estimation without mechanical sensors. More recently, Bekhiti [10] extended the adaptation ideas by validating a hyper-stability-based multivariable nonlinear adaptive control framework on induction motors, demonstrating strong experimental performance under non-ideal conditions. The integration of nonlinear state observers has further strengthened sensorless control by enabling full-state estimation based solely on measurable stator currents and voltages. Gauthier [11] and Busawon [12] contributed to the theoretical foundations of nonlinear observers, while Mansouri [13] tailored these concepts to IMs operating under field-oriented control. These observer designs significantly improve robustness and observability, especially at low speeds and under load disturbances.

More recently, artificial intelligence (AI) and neural networks have emerged as powerful tools in intelligent motor control. Neural models can approximate complex nonlinear mappings and adapt online, making them particularly suitable for estimation and control in dynamic environments. Beltran-Carbajal [14] and Moghadasian [15] proposed neural network-based controllers for IMs that adapt to system variations in real time. Sujatha and Vaisakh [16], along with Brandstetter and Kuchar [17], demonstrated the effectiveness of adaptive neuro-fuzzy inference systems (ANFIS) and radial basis function (RBF) networks in sensorless IM control. Acikgoz [18] implemented a real-time speed controller using a type-2 fuzzy neural network, while Hussain [19] introduced a type-II adaptive neuro-fuzzy controller with strong transient response capabilities.

Several studies have specifically integrated AI-based observers into MRAS or passivity-based frameworks. For instance, Medjadji [20] and Govindharaj [21] enhanced speed estimation using fuzzy and neuro-fuzzy observers. Parimalasundar [22] combined artificial neural networks with sliding mode control for power steering applications, and Mekrini [23] applied fuzzy logic in the control of asynchronous machines. Recent contributions by Sun [24] and Mienye [25] explored adaptive learning and recurrent neural network architectures to further strengthen robustness. Wu [26] addressed neural optimization by improving the convergence of online learning through gradient normalization—relevant for real-time motor drive applications.

Despite significant progress in adaptive, passivity-based, and learning-driven control of induction motors, a critical gap remains in unifying these approaches into a robust and practical sensorless control framework. Existing methods often treat passivity and learning separately or rely on fixed observers with limited adaptability. Neural and neuro-fuzzy controllers, while effective in handling nonlinearities, frequently lack formal stability guarantees or disrupt the energy-preserving structure of the system. These limitations become especially pronounced under low-speed operation, load transients, and parameter variations, which are common in industrial applications. Moreover, many nonlinear observers are structurally complex or computationally intensive, making real-time implementation difficult. This motivates the need for a cohesive control architecture that integrates the stability of passivity theory, the adaptability of neural learning, and the practical efficiency of nonlinear observer design in a fully sensorless setting.

This paper addresses the above challenges and presents a unified learning-driven passivity control strategy for induction motors using a nonlinear neural-adaptive observer. The main contributions are summarized as follows:

• A novel neural-based adaptive observer is proposed to estimate rotor flux and speed in real time, enhancing estimation accuracy under partial observability and eliminating the need for mechanical sensors.
• A nonlinear passivity-based control law is developed that integrates the neural learning loop while preserving passivity and ensuring global asymptotic stability.
• Experimental validation is conducted on a 1.1 kW induction motor under dynamic benchmarks (speed reversal, torque steps, and parameter variations), demonstrating up to 88.6% reduction in torque error and 75.3% improvement in flux estimation accuracy compared to classical methods.

This architecture bridges the gap between model-based passivity control and data-driven adaptation, offering a scalable and robust solution for high-performance industrial IM drives. The remainder of this paper is organized as follows: Section 2 presents the nonlinear dynamic model of the IM motor and the passivity framework. Section 3 details the neural adaptive observer design and learning algorithm. Section 4 introduces the passivity-based control strategy and its integration with the observer. Section 5 describes the experimental platform and test scenarios. Section 6 presents and discusses the results. Finally, Section 7 concludes the paper and outlines future directions.

2. Mathematical Framework and Problem Statement

Before deriving the control laws, it is essential to establish a physically consistent mathematical model of the induction motor that captures the key electromechanical dynamics relevant for observer and controller design. The modeling approach adopted here is based on an energy-based formulation that reflects the internal coupling between electrical and mechanical subsystems. This framework serves as the foundation for the passivity-based and adaptive structures developed later in the paper.

2.1. Dynamic Model of an Asynchronous Motor

In the mathematical analyses, it is standard practice to substitute the squirrel-cage motor, (rotor is characterized by uniform distributed conductors), with a theoretical equivalent rotor. This theoretical model matches the stator in the number of phases and features conductors distributed in a sinusoidal pattern. This method operates under the assumption that only the fundamental harmonic of the rotor’s magnetomotive force (MMF) is considered. Practical findings confirm that analyses and controller designs using this reduced model are applicable to the actual machine. However, special attention is required when modeling wound rotors using sinusoidally distributed windings [27,29]. The dynamic model is derived through the direct use of Euler-Lagrange formalism, as outlined in [10] which results to (differential Lagrangian system):

$$•\mathrm{E}\mathrm{l}\mathrm{e}\mathrm{c}\mathrm{t}\mathrm{r}\mathrm{i}\mathrm{c}\mathrm{a}\mathrm{l} \mathrm{P}\mathrm{a}\mathrm{r}\mathrm{t} {\mathit{D}}_e\left(q_m\right){\ddot{\mathit{q}}}_e+{\mathit{W}}_1\left(q_m\right){\dot{q}}_m{\dot{\mathit{q}}}_e+{\mathit{R}}_e{\dot{\mathit{q}}}_e={\mathit{M}}_e\mathbf{u}$$

(1)

$$•\mathrm{M}\mathrm{e}\mathrm{c}\mathrm{h}\mathrm{a}\mathrm{n}\mathrm{i}\mathrm{c}\mathrm{a}\mathrm{l} \mathrm{P}\mathrm{a}\mathrm{r}\mathrm{t}{D}_m{\ddot{q}}_m+R_m{\dot{q}}_m=\mathcal{T}-{\mathcal{T}}_L={\displaystyle \frac{1}{2}}{\dot{\mathit{q}}}_e^{\top }{\mathit{W}}_1\left(q_m\right){\dot{\mathit{q}}}_e-{\mathcal{T}}_L$$

(2)

Or more compactely we write:

$$\left[\begin{array}{cc}{\mathit{D}}_e\left(q_m\right)& 0\\ 0& D_m\end{array}\right]\left[\begin{array}{c}{\ddot{\mathit{q}}}_e\\ {\ddot{q}}_m\end{array}\right]+\left[\begin{array}{c}{\mathit{W}}_1\left(q_m\right){\dot{q}}_m{\dot{\mathit{q}}}_e\\ -{\displaystyle \frac{1}{2}}{\dot{\mathit{q}}}_e^{\top }{\mathit{W}}_1\left(q_m\right){\dot{\mathit{q}}}_e\end{array}\right]+\left[\begin{array}{cc}{\mathit{R}}_e& 0\\ 0& R_m\end{array}\right]\left[\begin{array}{c}{\dot{\mathit{q}}}_e\\ {\dot{q}}_m\end{array}\right]=\left[\begin{array}{c}{\mathit{M}}_e\mathbf{u}\\ -{\mathcal{T}}_L\end{array}\right]$$

(3)

The conventional two-phase αβ model, where the stator axes are stationary and the rotor axes rotate at the electrical angular velocity, is commonly used for a p pole-pair squirrel-cage induction motor with a uniform air-gap. This model defines the electrical parameters that are essential for analysis and control design.

$${\mathit{D}}_e\left(q_m\right)=\left[\begin{array}{c}{L}_s{\mathit{I}}_2\\ L_m{e}^{-\mathcal{J}p{q}_m}\end{array} \begin{array}{c}{L}_m{e}^{\mathcal{J}p{q}_m}\\ L_r{\mathit{I}}_2\end{array}\right]; {\mathit{R}}_e=\left[\begin{array}{c}{R}_s{\mathit{I}}_2\\ 0\end{array} \begin{array}{c}0\\ R_r{\mathit{I}}_2\end{array}\right]; {\mathit{M}}_e=\left[\begin{array}{c}{\mathit{I}}_2\\ 0\end{array}\right]; \mathcal{J}=\left[\begin{array}{c}0\\ 1\end{array} \begin{array}{c}-1\\ 0\end{array}\right]$$

$$e^{\mathcal{J}p{q}_m}=\left[\begin{array}{cc}\mathrm{c}\mathrm{o}\mathrm{s}\left(p{q}_m\right)& -\mathrm{s}\mathrm{i}\mathrm{n}\left(p{q}_m\right)\\ \mathrm{s}\mathrm{i}\mathrm{n}\left(p{q}_m\right)& \mathrm{c}\mathrm{o}\mathrm{s}\left(p{q}_m\right)\end{array}\right]; e^{-\mathcal{J}p{q}_m}={\left(e^{\mathcal{J}p{q}_m}\right)}^{\top }; \mathcal{J}=-{\mathcal{J}}^{\top }$$

$${\mathit{W}}_1\left(q_m\right)≔{\displaystyle \frac{\partial {\mathit{D}}_e\left(q_m\right)}{\partial q_m}}=\left[\begin{array}{cc}0& p{L}_m{e}^{\mathcal{J}p{q}_m}\\ -p{L}_m{e}^{-\mathcal{J}p{q}_m}& 0\end{array}\right]; \mathcal{J}{e}^{-\mathcal{J}p{q}_m}=e^{-\mathcal{J}p{q}_m}\mathcal{J}$$

The current vector is: ${\dot{\mathit{q}}}_e≔{\left[ {\dot{\mathit{q}}}_s^{\top }, {\dot{\mathit{q}}}_r^{\top }\right]}^{\top }={\left[ {\dot{q}}_{s\alpha }, {\dot{q}}_{s\beta }, {\dot{q}}_{r\alpha }, {\dot{q}}_{r\beta }\right]}^{\top }$ with$L_s, L_r, L_m>0$, $R_s, R_r>0$, $D_m>0$, $R_m\ge 0$. The flux vector $\mathit{\lambda }$ is related to the current vector ${\dot{\mathit{q}}}_e$ via $\mathit{\lambda }≔{\mathit{D}}_e\left(q_m\right){\dot{\mathit{q}}}_e={\left[ {\mathit{\lambda }}_s^{\top }, {\mathit{\lambda }}_r^{\top }\right]}^{\top }={\left[ {\lambda }_{s\alpha }, {\lambda }_{s\beta }, {\lambda }_{r\alpha }, {\lambda }_{r\beta }\right]}^{\top }$. To carry out the controller design, a suitable linear factorization of these workless forces into a form

$$\left[\begin{array}{c}{\mathit{W}}_1\left(q_m\right){\dot{q}}_m{\dot{\mathit{q}}}_e\\ -{\displaystyle \frac{1}{2}}{\dot{\mathit{q}}}_e^{\top }{\mathit{W}}_1\left(q_m\right){\dot{\mathit{q}}}_e\end{array}\right]=\mathit{C}\left(\mathit{q},\dot{\mathit{q}}\right)\dot{\mathit{q}}$$

(4)

In particular, the matrix $\mathit{C}\left(\mathit{q},\dot{\mathit{q}}\right)$ must be chosen such that $\dot{\mathit{D}}\left(\mathit{q}\right)-2\mathit{C}\left(\mathit{q},\dot{\mathit{q}}\right)={\mathit{C}}^{\top }\left(\mathit{q},\dot{\mathit{q}}\right)-\mathit{C}\left(\mathit{q},\dot{\mathit{q}}\right)$ is skew-symmetric, and its 3^rd-4^th rows are independent of the quantity${\dot{\mathit{q}}}_e$. These conditions are necessary for the subsequent stability analysis. By utilizing the transpose of $\mathcal{T}$ and${\mathit{W}}_1\left(q_m\right)$, it becomes apparent that the desired objectives can be attained with the following selection

$$\mathit{C}\left(\mathit{q},\dot{\mathit{q}}\right)=\left[\begin{array}{ccc}0& 0& \mathbf{f}\left(\mathit{q},\dot{\mathit{q}}\right)\\ -p{L}_m\mathcal{J}{e}^{-\mathcal{J}p{q}_m}& 0& 0\\ {\mathbf{f}}^{\top }\left(\mathit{q},\dot{\mathit{q}}\right)& 0& 0\end{array}\right];\mathrm{a}\mathrm{n}\mathrm{d}\mathbf{f}\left(\mathit{q},\dot{\mathit{q}}\right)=p{L}_m\mathcal{J}{e}^{\mathcal{J}p{q}_m}{\dot{\mathit{q}}}_r$$

(5)

This decomposition results in the next compact model [13]:

$$\mathit{D}\left(\mathit{q}\right)\ddot{\mathit{q}}+\mathit{C}\left(\mathit{q},\dot{\mathit{q}}\right)\dot{\mathit{q}}+\mathit{R}\dot{\mathit{q}}=\mathit{M}\mathbf{u}+\mathit{\xi }$$

(6)

with $\mathit{D}\left(q\right)=\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left({\mathit{D}}_e\left(q_m\right), D_m\right)$,$\mathbf{R}=\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left({\mathit{R}}_e, R_m\right)$, $\dot{\mathit{q}}={\left[{\dot{\mathit{q}}}_e^{\top }, {\dot{q}}_m\right]}^{\top }$, $\mathit{M}={\left[{\mathit{M}}_e^{\top }, 0\right]}^{\top }$, and $\mathit{\xi }={\left[0, -{\mathcal{T}}_L\right]}^{\top }$

2.2. Problem Formulation

Consider the induction motor model given by Eq(6) with outputs $\mathcal{T}$ and $‖{\mathit{\lambda }}_r‖$ to be controlled. Assume that:

A1 ${\mathcal{T}}_L$ is an unknown constant (switched step).

A2${\dot{q}}_{s\alpha }, {\dot{q}}_{s\beta }$ are accessible.

A3 All parameters are well-known and $R_m>0$.

Let ${\mathcal{T}}_d$ ​ be a differentiable function, with a known bounded, first-order derivative, and let β(t)>0 be the desired norm of the rotor flux, which is a twice differentiable with known, bounded 1^st and 2^nd-order derivatives. Under these assumptions, the goal is to design a nonlinear control law that guarantees the stability (i.e. internally) and the asymptotic tracking of both ${\mathcal{T}}_d\left(t\right)$ and β(t). This control law must guarantee that, for any initial conditions, all signals within the closed-loop system remain uniformly bounded, thereby achieving the desired system performance.

$$\underset{t\to \infty }{\mathrm{l}\mathrm{i}\mathrm{m}}=\left|\mathcal{T}-{\mathcal{T}}_d\right|=0, \underset{t\to \infty }{\mathrm{l}\mathrm{i}\mathrm{m}}=\left|\begin{array}{c} \\ ‖{\mathit{\lambda }}_r‖-\beta \end{array}\right|=0$$

(7)

To address the need for accurate flux and torque estimation using only stator-side measurements, we proceed to design a reduced-order nonlinear observer formulated in the stationary αβ reference frame. The observer is structured to be consistent with the motor’s energy-based dynamics and is later augmented with a recurrent neural network to enhance estimation accuracy under nonlinearities and uncertainties.

3. Nonlinear Observer Design

In this section, we review several flux observer designs from the literature, followed by the development of a tailored observer formulated specifically for the nonlinear framework considered in this work.

3.1. Flux Observer in a Closed Loop

A foundational reference in this area is the observer introduced by Verghese and Sanders [32], which has since inspired several variants, including those proposed by De Luca and Ulivi [33], Garcia in [34], and by Mansouri [4]. The general structure of this classical observer is given by:

$${\displaystyle \frac{d}{dt}}\left[\begin{array}{c}{\widehat{\mathit{i}}}_s\\ {\widehat{\mathit{\lambda }}}_r\end{array}\right]=\left[\begin{array}{c}-\gamma {\mathit{I}}_2\\ \left(L_m/T_r\right){\mathit{I}}_2\end{array}\begin{array}{c}\left(K/T_r\right){\mathit{I}}_2-Kp\mathsf{\Omega }\mathcal{J}\\ \left(-1/T_r\right){\mathit{I}}_2+p\mathsf{\Omega }\mathcal{J}\end{array}\right]\left[\begin{array}{c}{\widehat{\mathit{i}}}_s\\ {\widehat{\mathit{\lambda }}}_r\end{array}\right]+{\displaystyle \frac{1}{\sigma L_s}}\left[\begin{array}{c}{\mathit{I}}_2\\ 0_2\end{array}\right]{\mathbf{u}}_s+\left[\begin{array}{c}{k}_1{\mathit{I}}_2+k_{2}p\mathsf{\Omega }\mathcal{J}\\ k_3{\mathit{I}}_2+k_{4}p\mathsf{\Omega }\mathcal{J}\end{array}\right]{\mathit{e}}_s$$

(8)

where ${\mathit{i}}_s=d{\dot{\mathit{q}}}_s/d$, $\mathsf{\Omega }={\dot{q}}_m$, ${\widehat{\mathit{i}}}_s={[{\widehat{i}}_{s\alpha },{\widehat{i}}_{s\beta }]}^{\top }$, ${\widehat{\mathit{\lambda }}}_r={[{\widehat{\mathsf{\lambda }}}_{r\alpha },{\widehat{\mathsf{\lambda }}}_{r\beta }]}^{\top }$, ${\mathbf{u}}_s={[u_{s\alpha },u_{s\beta }]}^{\top }$, the $k_i’s$ being scalars and ${\mathit{e}}_s=\left({\widehat{\mathit{i}}}_s-{\mathit{i}}_s\right)$. Note that the gains depend on the speed in (8). We show the diagram block of this observer in Figure 1.

The resulting model for the observer error dynamics is then

$${\displaystyle \frac{d\mathit{e}}{dt}}=\left[\begin{array}{cc}\left(k_1-\gamma \right){\mathit{I}}_2+k_{2}p\mathsf{\Omega }\mathcal{J}& \left(K/T_r\right){\mathit{I}}_2-Kp\mathsf{\Omega }\mathcal{J}\\ k_3+\left(L_m/T_r\right){\mathit{I}}_2+k_{4}p\mathsf{\Omega }\mathcal{J}& \left(-1/T_r\right){\mathit{I}}_2+p\mathsf{\Omega }\mathcal{J}\end{array}\right]\mathit{e} \mathrm{w}\mathrm{i}\mathrm{t}\mathrm{h}\mathit{e}=\left[\begin{array}{c}{\widehat{\mathit{i}}}_s-{\mathit{i}}_s\\ {\widehat{\mathit{\lambda }}}_r-{\mathit{\lambda }}_r\end{array}\right]$$

(9)

Note that we can freely determine the scalar coefficients in equation (9). If $k_1$ and $k_3$ are selected such that $k_1-\gamma =-K/T_r$, $k_3+\left(L_m/T_r\right)=-k_4/T_r$ the error dynamics become

$${\displaystyle \frac{d\mathit{e}}{dt}}={\mathit{A}}_e\mathit{Q}\left(\mathsf{\Omega }\right)\mathit{e} \mathrm{w}\mathrm{i}\mathrm{t}\mathrm{h} {\mathit{A}}_e=\left[\begin{array}{cc}{k}_2{\mathit{I}}_2& -K{\mathit{I}}_2\\ k_4{\mathit{I}}_2& {\mathit{I}}_2\end{array}\right]; \mathit{Q}\left(\mathsf{\Omega }\right)=\left[\begin{array}{cc}p\mathsf{\Omega }\mathcal{J}-\left({\mathit{I}}_2/T_r\right)& 0\\ 0& p\mathsf{\Omega }\mathcal{J}-\left({\mathit{I}}_2/T_r\right)\end{array}\right]$$

(10)

We select $k_2$ and $k_4$ to place the eigenvalues of ${\mathit{A}}_e$ in arbitrary positions. Note that the characteristic polynomial of ${\mathit{A}}_e$ is ${\left[s^2-\left(1+k_2\right)s+k_2+k_{4}K\right]}^2$. If the eigenvalues of ${\mathit{A}}_e$ are $s_1$ (twice) and $s_2$ (twice), then the eigenvalues of time varying matrix ${\mathit{A}}_e\mathit{Q}\left(\mathsf{\Omega }\right)$ are $\left[\left(-1/T_r\right)\pm p\mathsf{\Omega }j\right]s_1$ and $\left[\left(-1/T_r\right)\pm p\mathsf{\Omega }j\right]s_2$ using the slow dynamics theory.

3.2. The Proposed Nonlinear Intelligent Observer Design

Passivity-based control strategies are widely recognized for their reliance on rotor flux information, which is often challenging to measure directly. Consequently, considerable research efforts have been directed toward developing methods for estimating these values. Now, we introduce the design of a reduced-order flux observer for an asynchronous motor operating in the αβ reference frame. The proposed observer requires inputting only the stator voltage, current, and rotor speed. Specifically, this new nonlinear-observer incorporates the speed measurements, ensuring that only the electrical quantities are required in the estimation process [13]. Consider the αβ model for induction motor:

$$\left\{ \begin{array}{c}{\displaystyle \frac{d}{dt}}\left[\begin{array}{c}\begin{array}{c}\begin{array}{c}{\dot{q}}_{s\alpha }\\ \end{array}\\ \end{array}\\ \begin{array}{c}\begin{array}{c}\begin{array}{c}{\dot{q}}_{s\beta }\\ \end{array}\\ \end{array}\\ \begin{array}{c}\begin{array}{c}{\lambda }_{r\alpha }\\ \end{array}\\ \end{array}\end{array}\\ {\lambda }_{r\beta }\end{array}\right]=\left[\begin{array}{c}-\gamma {\dot{q}}_{s\alpha }+{\displaystyle \frac{K}{T_r}}{\lambda }_{r\alpha }+Kp{\dot{q}}_m{\lambda }_{r\beta }+{\displaystyle \frac{1}{\sigma L_s}}{u}_{s\alpha }\\ \begin{array}{c}-\gamma {\dot{q}}_{s\beta }-Kp{\dot{q}}_m{\lambda }_{r\alpha }+{\displaystyle \frac{K}{T_r}}{\lambda }_{r\beta }+{\displaystyle \frac{1}{\sigma L_s}}{u}_{s\beta }\\ \begin{array}{c}{\displaystyle \frac{L_m}{T_r}}{\dot{q}}_{s\alpha }-{\displaystyle \frac{1}{T_r}}{\lambda }_{r\alpha }-p{\dot{q}}_m{\lambda }_{r\beta }\\ {\displaystyle \frac{L_m}{T_r}}{\dot{q}}_{s\beta }+p{\dot{q}}_m{\lambda }_{r\alpha }-{\displaystyle \frac{1}{T_r}}{\lambda }_{r\beta }\end{array}\end{array}\end{array}\right]\\ {\left[y_1{y}_2\right]}^{\top }={\left[{\dot{q}}_{s\alpha }{\dot{q}}_{s\beta }\right]}^{\top } \end{array} \right\}$$

(11)

where

$$T_r={\displaystyle \frac{L_r}{R_r}}; \sigma =1-{\displaystyle \frac{L_m^2}{L_s{L}_r}}; K={\displaystyle \frac{L_m}{\sigma L_s{L}_r}}; \gamma ={\displaystyle \frac{R_s}{\sigma L_s}}+{\displaystyle \frac{R_r{L}_m^2}{\sigma L_s{L}_r^2}}$$

The system (11) is the form

$$\dot{\mathit{z}}=\mathbf{F}\left(\mathsf{\Omega }\right)\mathit{z}+\mathbf{G}\left(\mathbf{u},\mathsf{\Omega },\mathit{z}\right) \mathrm{a}\mathrm{n}\mathrm{d}\mathbf{y}=\mathbf{C}\mathit{z}$$

(12)

with

$$\left\{\begin{array}{c}\mathbf{F}\left(\mathsf{\Omega }\right)=\left[\begin{array}{cc}0& {\mathbf{F}}_1\left(\mathsf{\Omega }\right)\\ 0& 0\end{array}\right];\hspace{0.17em}\mathbf{G}\left(\mathbf{u},\mathsf{\Omega },\mathit{z}\right)=\left[\begin{array}{c}{\mathbf{g}}_1(u,\Omega ,z)\\ {\mathbf{g}}_2(u,\Omega ,z)\end{array}\right];\hspace{0.17em}\mathbf{u}=\left[\begin{array}{c}{u}_{s\alpha }\\ u_{s\beta }\end{array}\right]; \mathit{z}=\left[\begin{array}{c}{\mathit{z}}_1\\ {\mathit{z}}_2\end{array}\right];\hspace{0.17em}\hspace{0.17em}\mathbf{y}=\left[\begin{array}{c}{\dot{q}}_{s\alpha }\\ {\dot{q}}_{s\beta }\end{array}\right]; {\mathit{z}}_1=\left[\begin{array}{c}{\dot{q}}_{s\alpha }\\ {\dot{q}}_{s\beta }\end{array}\right];\hspace{0.17em}{\mathit{z}}_2=\left[\begin{array}{c}{\lambda }_{r\alpha }\\ {\lambda }_{r\beta }\end{array}\right];\hspace{0.17em}\mathsf{\Omega }={\dot{q}}_m\\ {\mathbf{F}}_1\left(\mathsf{\Omega }\right)=\left[\begin{array}{cc}{\displaystyle \frac{K}{T_r}}& Kp\mathsf{\Omega }\\ -Kp\mathsf{\Omega }& {\displaystyle \frac{K}{T_r}}\end{array}\right]; {\mathbf{g}}_1\left(\mathbf{u},\mathsf{\Omega },\mathit{z}\right)=\left[\begin{array}{c}-\gamma {\dot{q}}_{s\alpha }+{\displaystyle \frac{1}{\sigma L_s}}{u}_{s\alpha }\\ -\gamma {\dot{q}}_{s\beta }+{\displaystyle \frac{1}{\sigma L_s}}{u}_{s\beta }\end{array}\right]; {\mathbf{g}}_2\left(\mathbf{u},\mathsf{\Omega },\mathit{z}\right)=\left[\begin{array}{c}{\displaystyle \frac{L_m}{T_r}}{\dot{q}}_{s\alpha }-{\displaystyle \frac{1}{T_r}}{\lambda }_{r\alpha }-p\Omega {\lambda }_{r\beta }\\ {\displaystyle \frac{L_m}{T_r}}{i}_{s\beta }+p\Omega {\lambda }_{r\alpha }-{\displaystyle \frac{1}{T_r}}{\lambda }_{r\beta }\end{array}\right]; \mathrm{a}\mathrm{n}\mathrm{d} {\mathbf{C}}^{\top }=\left[ \begin{array}{c}1\\ \begin{array}{c}1\\ 0\end{array}\\ 0\end{array}\right]\end{array}\right\}$$

Here, we make the following assumptions as posed by K. Busawon (see [12])

(A4) There exists a class $\mathcal{U}$ of bounded admissible controls, a compact set $\mathcal{K}\subset {\mathbb{R}}^{2\times 2}$ and positive constants α, β > 0 such that for every $\mathbf{u}\in \mathcal{U}$ and every output $\mathbf{y}\left(t\right)$ associated with $\mathbf{u}$ and with an initial state $\mathit{z}\left(0\right)\in \mathcal{K}$, we have $0<\alpha {\mathit{I}}_2\le {\mathbf{F}}_1^{\top }{\mathbf{F}}_1\le \beta {\mathit{I}}_2$.

(A5) $\mathsf{\Omega }\left(t\right)$ and its time derivative $d\mathsf{\Omega }\left(t\right)/dt$are bounded.

(A6) The matrix ${\mathbf{F}}_1\left(\mathsf{\Omega }\right)$ is of class ${\mathcal{C}}^r$, $r \ge 1$ with respect to their arguments.

(A7) Functions ${\mathbf{g}}_1$, ${\mathbf{g}}_2$ are global Lipschitz with respect to $\mathit{z}$ uniformly in $\mathbf{u}$ and $\mathsf{\Omega }$.

We characterize the observer design for the system (12) in the following theorem

THEOREM: 1 Assume that the system (12) satisfies Assumptions (A4) to (A7). Then there exists $\theta >0$ such that the system

$$\dot{\widehat{z}}=\mathbf{F}\left(\mathsf{\Omega }\right)\widehat{\mathit{z}}+\mathbf{G}\left(\mathbf{u},\mathsf{\Omega },\widehat{\mathit{z}}\right)-{\mathbf{L}}_{\theta }\left(\mathsf{\Omega }\right)\left(\widehat{\mathit{z}}-\mathit{z}\right) \mathrm{w}\mathrm{i}\mathrm{t}\mathrm{h} {\mathbf{L}}_{\theta }\left(\mathsf{\Omega }\right)={\mathsf{\Lambda }}^{-1}\left(\mathsf{\Omega }\right){\mathbf{S}}_{\theta }^{-1}{\mathbf{C}}^{\top }\mathbf{C}$$

(13)

is an exponential observer for the system (12), where ${\mathit{S}}_{\theta }$ is the unique solution of the algebraic Lyapunov equation $\theta {\mathit{S}}_{\theta }+{\mathit{A}}^{\top }{\mathit{S}}_{\theta }+{\mathit{S}}_{\theta }\mathit{A}-{\mathit{C}}^{\top }\mathit{C}=0$, $\mathsf{\Lambda }\left(\mathsf{\Omega }\right)=\mathrm{b}\mathrm{l}\mathrm{k}\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left[{\mathit{I}}_2, {\mathit{F}}_1\left(\mathsf{\Omega }\right)\right]$, ${\mathbf{S}}_{\theta }^{-1}{\mathbf{C}}^{\top }={\left[2\theta {\mathit{I}}_2 {\theta }^2{\mathit{I}}_2\right]}^{\top }$, $\mathit{A}={\left[{\left[0_2 {\mathit{I}}_2\right]}^{\top }, {\left[0_2 0_2 \right]}^{\top }\right]}^{\top }$ and the explicit formula of ${\mathbf{L}}_{\theta }\left(\mathsf{\Omega }\right)$ is

$${\mathbf{L}}_{\theta }\left(\mathsf{\Omega }\right)={\mathsf{\Lambda }}^{-1}\left(\mathsf{\Omega }\right){\mathbf{S}}_{\theta }^{-1}{\mathbf{C}}^{\top }\mathbf{C}={\left[{\left[2\theta {\mathit{I}}_2 0_2\right]}^{\top }, {\left[{\theta }^2{\mathbf{F}}_1^{-1}\left(\mathsf{\Omega }\right) 0_2\right]}^{\top }\right]}^{\top },\mathrm{a}\mathrm{n}\mathrm{d} \theta >0$$

(14)

Moreover, we can make the dynamics of this observer arbitrarily fast. ∎

So, by multiplying the left and right-hand sides of (14) by $\mathsf{\Lambda }\left(\mathsf{\Omega }\right)$ and ${\mathsf{\Lambda }}^{\top }\left(\mathsf{\Omega }\right)$ respectively, the following algebraic equation holds:

$$\theta {\widehat{\mathit{S}}}_{\theta }+{\mathit{F}}^{\top }\left(\mathsf{\Omega }\right){\widehat{\mathit{S}}}_{\theta }+{\widehat{\mathit{S}}}_{\theta }\mathit{F}\left(\mathsf{\Omega }\right)-{\mathit{C}}^{\top }\mathit{C}=0 \mathrm{w}\mathrm{i}\mathrm{t}\mathrm{h}{\widehat{\mathit{S}}}_{\theta }={\mathsf{\Lambda }}^{\top }\left(\mathsf{\Omega }\right){\mathit{S}}_{\theta }\mathsf{\Lambda }\left(\mathsf{\Omega }\right)$$

(15)

The choice of θ permits the pole placement of the observer according to the speed.

In summary the system can be written as ${\dot{\mathbf{x}}}_{\mathrm{I}\mathrm{M}}\left(t\right)={\mathit{A}}_{\mathrm{I}\mathrm{M}}\left(\mathsf{\Omega }\right){\mathbf{x}}_{\mathrm{I}\mathrm{M}}\left(t\right)+{\mathit{B}}_{\mathrm{I}\mathrm{M}}\left(\mathsf{\Omega }\right){\mathbf{u}}_s\left(t\right)$ or

$${\displaystyle \frac{d}{dt}}\left[\begin{array}{c}{\dot{\mathit{q}}}_s\\ {\mathit{\lambda }}_r\end{array}\right]=\left[\begin{array}{c}-\gamma {\mathit{I}}_2\\ \left(L_m/T_r\right){\mathit{I}}_2\end{array}\begin{array}{c}\left(K/T_r\right){\mathit{I}}_2-Kp\mathsf{\Omega }\mathcal{J}\\ \left(-1/T_r\right){\mathit{I}}_2+p\mathsf{\Omega }\mathcal{J}\end{array}\right]\left[\begin{array}{c}{\dot{\mathit{q}}}_s\\ {\mathit{\lambda }}_r\end{array}\right]+{\displaystyle \frac{1}{\sigma L_s}}\left[\begin{array}{c}{\mathit{I}}_2\\ 0_2\end{array}\right]{\mathbf{u}}_s;\mathbf{y}={\dot{\mathit{q}}}_s$$

The observer is ${\dot{\mathbf{x}}}_{\mathrm{O}}\left(t\right)={\mathit{A}}_{\mathrm{I}\mathrm{M}}\left(\mathsf{\Omega }\right){\mathbf{x}}_{\mathrm{O}}\left(t\right)+{\mathit{B}}_{\mathrm{I}\mathrm{M}}\left(\mathsf{\Omega }\right){\mathbf{u}}_s\left(t\right)+{\mathit{L}}_{\theta }\left(\mathsf{\Omega }\right){\mathit{e}}_s\left(t\right)$ or

$${\displaystyle \frac{d}{dt}}\left[\begin{array}{c}{\widehat{\dot{\mathit{q}}}}_s\\ {\widehat{\mathit{\lambda }}}_r\end{array}\right]=\left[\begin{array}{c}-\gamma {\mathit{I}}_2\\ \left(L_m/T_r\right){\mathit{I}}_2\end{array}\begin{array}{c}\left(K/T_r\right){\mathit{I}}_2-Kp\mathsf{\Omega }\mathcal{J}\\ \left(-1/T_r\right){\mathit{I}}_2+p\mathsf{\Omega }\mathcal{J}\end{array}\right]\left[\begin{array}{c}{\widehat{\dot{\mathit{q}}}}_s\\ {\widehat{\mathit{\lambda }}}_r\end{array}\right]+{\displaystyle \frac{1}{\sigma L_s}}\left[\begin{array}{c}{\mathit{I}}_2\\ 0_2\end{array}\right]{\mathbf{u}}_s+\left[\begin{array}{cc}2\theta {\mathit{I}}_2& 0\\ {\theta }^2\mathbf{h}\left(\mathsf{\Omega }\right)& 0\end{array}\right]{\mathit{e}}_s$$

So the error equation is ${\dot{\mathbf{e}}}_{\mathrm{s}}\left(t\right)={\mathit{A}}_{\mathrm{O}}\left(\mathsf{\Omega }\right){\mathbf{e}}_{\mathrm{s}}\left(t\right)$ with ${\mathit{A}}_{\mathrm{O}}\left(\mathsf{\Omega }\right)={\mathit{A}}_{\mathrm{I}\mathrm{M}}\left(\mathsf{\Omega }\right)-{\mathit{L}}_{\theta }\left(\mathsf{\Omega }\right)$ or

$${\displaystyle \frac{d}{dt}}{\mathit{e}}_s=\left[\begin{array}{c}2\theta {\mathit{I}}_2-\gamma {\mathit{I}}_2\\ {\theta }^2\mathbf{h}\left(\mathsf{\Omega }\right)+\left(L_m/T_r\right){\mathit{I}}_2\end{array}\begin{array}{c}\left(K/T_r\right){\mathit{I}}_2-Kp\mathsf{\Omega }\mathcal{J}\\ \left(-1/T_r\right){\mathit{I}}_2+p\mathsf{\Omega }\mathcal{J}\end{array}\right]{\mathit{e}}_s\mathrm{a}\mathrm{n}\mathrm{d}\mathbf{h}\left(\mathsf{\Omega }\right)={\displaystyle \frac{\left(K/T_r\right)\mathit{I}+Kp\mathsf{\Omega }\mathcal{J}}{{\left(K/T_r\right)}^2+{\left(Kp\mathsf{\Omega }\right)}^2}}$$

To incorporate learning-based compensation for model uncertainties and to enhance the robustness of the nonlinear observer, we introduce a radial basis function neural networks (RBFNN) term into the observer dynamics. The network approximates unknown disturbances or residual nonlinearities using a regressor $\mathit{\varphi }\left(t\right)\in {\mathbb{R}}^{n_{\varphi }}$ composed of observable and estimated signals. Let $\mathsf{\Delta }\left(\mathit{\varphi }\right)\in {\mathbb{R}}^{n_z}$ represent an unknown nonlinear function (e.g., modeling error or disturbance). The universal approximation theorem states that, for any continuous function $\mathsf{\Delta }\left(\mathit{\varphi }\right)$, there exists a feedforward neural network (three-layer NN) of the form: $\mathsf{\Delta }\left(\mathit{\xi }\right)\approx \widehat{\mathsf{\Delta }}\left(\mathit{\xi }\right)={⅀}_{i=1}^{n_h}{\mathcal{W}}_i{\psi }_i\left(\mathit{\xi }\right)$ with $\mathit{\xi }=\mathcal{V}\mathit{\varphi }+\mathcal{B}$ or we omit the intermediary variable $\mathit{\xi }$ as $\mathsf{\Delta }\left(\mathit{\varphi }\right)={\mathcal{W}}^{\top }\mathsf{\Psi }\left(\mathit{\varphi }\right)$ where

• $\mathsf{\Psi }\left(\mathit{\varphi }\right)={\left[{\psi }_1\left(\mathit{\varphi }\right),\dots ,{\psi }_{n_h}\left(\mathit{\varphi }\right)\right]}^{\top }\in {\mathbb{R}}^{n_h}$ is the activation vector,
• $\mathcal{W}=\arg \underset{\mathbf{W}\in {\mathbb{R}}^{n_h\times n_z}}{\min }\left\{\underset{\mathit{\varphi }\in {\mathsf{\Omega }}_{\varphi }}{\mathrm{s}\mathrm{u}\mathrm{p}}\left|\mathsf{\Delta }\left(\mathbf{x}\right)-{\mathit{W}}^{\top }\mathsf{\Psi }\left(\mathit{\varphi }\right)\right|\right\}$ is the output weight matrix,
• ${\psi }_i\left(\mathit{\varphi }\right)$ are nonlinear sigmoids activation functions (e.g., tanh, ReLUs, etc.).
• $\mathcal{V}\in {\mathbb{R}}^{n_h\times n}$: stacked weight matrix ($i^{th}$ row is ${\mathbf{v}}_i^{\top }$), $\mathcal{B}\in {\mathbb{R}}^{n_h}$: stacked biases,
• The exponential and division are element-wise.

The coefficients ${\mathbf{v}}_i^{\top }$ and $b_i$ are the hidden layer parameters (${\mathbf{v}}_i^{\top }\in {\mathbb{R}}^n$ Input-to-hidden weights "nonlinear layer" and $b_i\in \mathbb{R}$ are the Bias of neuron $i$). Now, the resulting enhanced observer is given by:

$$\dot{\widehat{z}}=\mathbf{F}\left(\mathsf{\Omega }\right)\widehat{\mathit{z}}+\mathbf{G}\left(\mathbf{u},\mathsf{\Omega },\widehat{\mathit{z}}\right)-{\mathsf{\Lambda }}^{-1}\left(\mathsf{\Omega }\right){\mathbf{S}}_{\theta }^{-1}{\mathbf{C}}^T\left(\mathbf{C}\widehat{\mathit{z}}-\mathbf{y}\right)+{\mathcal{W}}^{\top }\mathsf{\Psi }\left(\mathit{\varphi }\left(t\right)\right)$$

(16)

Here, $\mathcal{W}\in {\mathbb{R}}^{n_h\times n_z}$ is the adaptive weight matrix of the neural network and the vector $\mathsf{\Psi }\left(\mathit{\varphi }\left(t\right)\right)\in {\mathbb{R}}^{n_h}$ is a vector of radial basis activation functions, typically Gaussian (${\mathsf{\Psi }}_i=\exp \left[-‖\mathit{\varphi }\left(t\right)-{\mathit{c}}_i‖/{\sigma }_i^2\right]$ and $\mathit{\varphi }\left(t\right)$ is the regressor vector, e.g., $\mathit{\varphi }\left(t\right)={\left[{\mathit{z}}^{\top },{\mathbf{u}}^{\top },\mathsf{\Omega }\right]}^{\top }$). The network parameters are updated online via adaptation rules derived from Lyapunov stability analysis to ensure convergence and passivity preservation. The RBF neural term ${\mathbf{f}}_{RBF}\left(\mathit{\varphi }\left(t\right)\right)={\mathcal{W}}^{\top }\mathsf{\Psi }\left(\mathit{\varphi }\left(t\right)\right)$ is added as a correction to the standard observer dynamics, to compensate for unknown nonlinearities or disturbances not captured by $\mathbf{F}$ and $\mathbf{G}$. This term acts like a learning system and is updated online (see Figure 2).

So the RBFNN becomes: $\mathsf{\Delta }\left(\mathit{\varphi }\right)={⅀}_{i=1}^p{\mathbf{w}}_i{\psi }_i\left(\mathit{\varphi }\right)={\mathcal{W}}^{\top }\mathsf{\Psi }\left(\mathit{\varphi }\right)$. In order to train the RBFNN we should to minimize the squared error objective function: $\mathcal{J}={‖{\mathbf{e}}_n‖}^2={‖{\widehat{\mathbf{y}}}_n-{\mathbf{y}}_n‖}^2$. To adapt ${\mathit{c}}_j$, perform $\partial \mathcal{J}/\partial {\mathit{c}}_j=-{\mathbf{e}}_n^{\top }{\mathbf{w}}_j.\partial {\psi }_j\left({\mathit{\phi }}_n\right)/{\mathit{c}}_j=-{\mathbf{e}}_n^{\top }{\mathbf{w}}_j.{\psi }_j\left({\mathit{\phi }}_n\right).\left({\mathit{\phi }}_n-{\mathit{c}}_j\right)/{\sigma }_j^2$. Again, for ${\sigma }_j$ we have​: $\partial \mathcal{J}/\partial {\sigma }_j=-{\mathbf{e}}_n^{\top }{\mathbf{w}}_j.\partial {\psi }_j\left({\mathit{\phi }}_n\right)/\partial {\sigma }_j=-{\mathbf{e}}_n^{\top }{\mathbf{w}}_j.{\psi }_j\left({\mathit{\phi }}_n\right).{‖{\mathit{\phi }}_n-{\mathit{c}}_j‖}^2/{\sigma }_j^3$. Here is a formal description of this training algorithm

3.3. The Proposed Fuzzy-MRAS Estimator

Rotor speed $\widehat{\mathsf{\Omega }}$ can be estimated using a model reference adaptive system (MRAS) structure composed of two estimators: a reference model and an adaptive model. Both models independently compute the rotor flux linkage components ${\lambda }_{r\alpha }$ and ${\lambda }_{r\beta }$ in the stationary αβ reference frame. By comparing the flux-linkage estimates, the adaptive model adjusts its internal speed estimate to minimize the discrepancy, thereby converging to the true rotor speed [8,9]. The reference model, derived from the stator voltage equations, does not depend on rotor speed and is therefore suitable for speed-independent flux estimation. Its expression in the stationary frame is given by:

$${\displaystyle \frac{d{\lambda }_{r\alpha }^{\mathrm{r}\mathrm{e}\mathrm{f}}}{dt}}={\displaystyle \frac{L_r}{L_m}}\left[u_{s\alpha }-R_s{i}_{s\alpha }-\sigma L_s{\displaystyle \frac{d{i}_{s\alpha }}{dt}}\right];{\displaystyle \frac{d{\lambda }_{r\beta }^{\mathrm{r}\mathrm{e}\mathrm{f}}}{dt}}={\displaystyle \frac{L_r}{L_m}}\left[u_{s\beta }-R_s{i}_{s\beta }-\sigma L_s{\displaystyle \frac{d{i}_{s\beta }}{dt}}\right]$$

(17)

In contrast, the adaptive model, formulated from the rotor voltage equations, explicitly depends on the rotor flux components and the estimated speed $\widehat{\Omega }$. It is expressed as:

$${\displaystyle \frac{d{\widehat{\lambda }}_{r\alpha }}{dt}}={\displaystyle \frac{L_m}{T_r}}{i}_{s\alpha }-{\displaystyle \frac{1}{T_r}}{\widehat{\lambda }}_{r\alpha }-p\widehat{\mathsf{\Omega }}{\widehat{\lambda }}_{r\beta }; {\displaystyle \frac{d{\widehat{\lambda }}_{r\beta }}{dt}}={\displaystyle \frac{L_m}{T_r}}{i}_{s\beta }-{\displaystyle \frac{1}{T_r}}{\widehat{\lambda }}_{r\beta }+p\widehat{\mathsf{\Omega }}{\widehat{\lambda }}_{r\alpha }$$

(18)

The angular discrepancy between the estimated and reference flux vectors serves as the adaptation signal, which is processed by a linear PI controller to adjust the estimated speed. Applying the Popov $\mathrm{h}\mathrm{y}\mathrm{p}\mathrm{e}\mathrm{r}\mathrm{s}\mathrm{t}\mathrm{a}\mathrm{b}\mathrm{i}\mathrm{l}\mathrm{i}\mathrm{t}\mathrm{y}$ criterion [10,20] to guarantee global asymptotic stability of the closed-loop estimation system, the adaptive law is given by:

$${\displaystyle \frac{d\widehat{\mathsf{\Omega }}}{dt}}=K_p{\displaystyle \frac{d\epsilon \left(t\right)}{dt}}+K_i\epsilon \left(t\right);\mathrm{w}\mathrm{i}\mathrm{t}\mathrm{h} \epsilon \left(t\right)=\left[{\widehat{\lambda }}_{r\alpha }{\lambda }_{r\beta }^{\mathrm{r}\mathrm{e}\mathrm{f}}-{\widehat{\lambda }}_{r\beta }{\lambda }_{r\alpha }^{\mathrm{r}\mathrm{e}\mathrm{f}}\right]$$

(19)

Fuzzy logic is particularly advantageous for systems that are difficult to model analytically, due to limited data availability, incomplete information, or the inherent complexity of the underlying process [16]. These systems can be easily extended by adding new rules, allowing for performance enhancements and the inclusion of additional features. In conventional MRAS-based speed observers, a proportional-integral (PI) controller is typically employed within the adaptation mechanism to estimate rotor speed by minimizing the error between the reference and adaptive models [18]. In the proposed approach, the PI controller is replaced by a fuzzy logic controller, as depicted in Figure 3. The complete block diagram of the fuzzy controller is presented in Figure 4, and its internal structure is illustrated in Figure 5. The fuzzy logic system operates by first processing the input variables through a fuzzification stage, where crisp input values are mapped to fuzzy sets using predefined membership functions and scale transformations. This scale mapping adjusts the range of the input variables to match the domain expected by the controller. The rule base stores the fuzzy inference rules, while the inference mechanism applies these rules to derive fuzzy conclusions. Common inference strategies include the Mamdani method, which typically uses the logical AND operator for rule evaluation, and the Takagi-Sugeno Model (TSM), which generates functional outputs based on the input variables [23,36]. Finally, the defuzzification stage converts the fuzzy output set into a crisp control signal suitable for direct application to the system.

Knowing that gains $K_1$, $K_2$ and $K_3$ are scaling factors, they should be chosen carefully after doing trial and error technique to get the optimal performance of the controller. From the expert estimator, the following memberships of fuzzy sets are used. The rules supporting this system are expressed in the Table 1, where NB: Negative Big, PB: Positive Big, NS: Negative Small, PS: Positive Small and ZE: Zero Equal. The number of inputs to the fuzzy controller can be 2 or 3. In the proposed method to reduce the complexity the two inputs are considered. The linguistics values that can be considered for each input are 3, 5, 7 or 9 as the case may be. Since there are two inputs, the numbers of rules that are possible are $\mathrm{N}\times \mathrm{N}$ or ${\mathrm{N}}^2$ which gives 9, 25, 49 or 81 rules for the above linguistic values respectively. The choice of linguistic value should be made keeping into consideration that the control should be good and computation taken by the fuzzy controller should not be high make the process slow. The linguistic value 3 which gives 9 rules is too small to get the better control. Figure 6 shows the fuzzy membership functions and control surface of the proposed adaptive system.

, change in $\epsilon$ and the control surface of the fuzzy logic estimator

The fuzzy inference rules used in the proposed system are summarized in Table 1.

4. Passivity Based Output-Feedback Controller

4.1. Non-Intelligent Sensorless Adaptive Passivity-Based Controller Design

The fundamental concept of passivity lies in shaping the total energy of the system and introducing an additional damping term to guide its dynamic behavior. Using the Euler–Lagrange modeling (ELM) framework, the total energy expression can be systematically formulated and then shaped toward a desired minimum. If the control input is designed to inject a supplementary dissipative term into the system, the convergence to the equilibrium state can be significantly improved beyond what is achieved through the system’s natural dissipation alone. In essence, passivity is an intrinsic property of certain physical systems and is closely linked to the notion of input–output stability [10,35].

: ${\mathbb{R}}^n\to \mathbb{R}$, called a storage function, with $V\left(0\right)=0$, as for all $\mathbf{u}\in \mathcal{U}$, $\mathbf{x}\left(0\right)={\mathbf{x}}_0\in \mathcal{X}$ and $t \ge 0$, it is satisfied

$$V\left({\mathbf{x}}_t\right)-V\left(0\right)\le {\int }_0^t{\mathbf{y}}^{\top }\left(\xi \right)\mathbf{u}\left(\xi \right)d\xi$$

(20)

where $V$ is a continuous storage function with continuous derivatives of order $r$. Condition (20) can also be expressed as $\dot{V}\left({\mathbf{x}}_t\right)\le {\mathbf{y}}^{\top }\left(t\right)\mathbf{u}\left(t\right)$ [29]. The system is said to be lossless when $\dot{V}\left({\mathbf{x}}_t\right)={\mathbf{y}}^{\top }\left(t\right)\mathbf{u}\left(t\right)$. ∎

In cases where the nature of the system does not involve any motion of coordinates or changes in reference frames (systems that are not explicitly time dependent), the Hamiltonian $\mathcal{H}$ is simply the total energy of the system. The Hamiltonian of IM (i.e. conservative system with zero potential energy) equals the sum of the kinetic energy $\mathcal{K}$ and the potential energy $\mathcal{P}$

$$\mathcal{H}\left({\dot{\mathit{q}}}_e,{\dot{q}}_m, q_m\right)={\displaystyle \frac{1}{2}}{\dot{\mathit{q}}}_e^{\top }{\mathit{D}}_e\left(q_m\right){\dot{\mathit{q}}}_e+{\displaystyle \frac{1}{2}}{D}_m{\dot{q}}_m^2$$

(21)

The time derivative of this Hamiltonian function along the trajectories of the system can be expressed as

$${\displaystyle \frac{d\mathcal{H}}{dt}}={\displaystyle \frac{1}{2}}{\dot{\mathit{q}}}_e^{\top }\left\{{\displaystyle \frac{\partial {\mathit{D}}_e}{\partial q_m}}{\dot{q}}_m\right\}{\dot{\mathit{q}}}_e+{\dot{\mathit{q}}}_e^{\top }{\mathit{D}}_e\left(q_m\right){\ddot{\mathit{q}}}_e+D_m{\ddot{q}}_m{\dot{q}}_m$$

(22)

Using the previous equations ${\mathit{D}}_e{\ddot{\mathit{q}}}_e={\mathbf{u}}_e-{\mathit{W}}_1{\dot{q}}_m{\dot{\mathit{q}}}_e-{\mathit{R}}_e{\dot{\mathit{q}}}_e;$ $D_m{\ddot{q}}_m=\mathcal{T}-{\mathcal{T}}_L-R_m{\dot{q}}_m$ and $\mathcal{T}=\left[{\dot{\mathit{q}}}_e^{\mathrm{\top }}{\mathit{W}}_1{\dot{\mathit{q}}}_e\right]/2$ we get

$$\left\{\begin{array}{c}{\displaystyle \frac{d\mathcal{H}}{dt}}=-{\mathcal{T}}_L{\dot{q}}_m+{\dot{\mathit{q}}}_e^{\top }{\mathbf{u}}_e-\left[{\dot{\mathit{q}}}_e^{\top }{\mathit{R}}_e{\dot{\mathit{q}}}_e+{\dot{q}}_m^{\top }{R}_m{\dot{q}}_m\right]\\ =-{\left[{\dot{\mathit{q}}}_e^{\top } {\dot{q}}_m^{\top }\right]}^{\top }\left[\begin{array}{cc}{\mathit{R}}_e& 0\\ 0& R_m\end{array}\right]\left[\begin{array}{c}{\dot{\mathit{q}}}_e\\ {\dot{q}}_m\end{array}\right]+\left[{\dot{\mathit{q}}}_e^{\top } {\dot{q}}_m^{\top }\right]\left[\begin{array}{c} {\mathbf{u}}_e\\ -{\mathcal{T}}_L\end{array} \right]\end{array}\right\}$$

(23)

So, the derivative of Hamiltonian function of IM is given by:

$$d\mathcal{H}/dt=-{\dot{\mathit{q}}}^{\top }\mathit{R}\dot{\mathit{q}}+{\mathbf{y}}^{\top }\mathbf{u};\mathrm{w}\mathrm{i}\mathrm{t}\mathrm{h}\mathbf{u}=\left[\begin{array}{c} {\mathbf{u}}_e\\ -{\mathcal{T}}_L\end{array} \right];\mathbf{y}=\left[\begin{array}{c}{\dot{\mathit{q}}}_e\\ {\dot{q}}_m\end{array} \right]$$

(24)

$$\begin{array}{c}\mathrm{S}\mathrm{t}\mathrm{o}\mathrm{r}\mathrm{e}\mathrm{d}\mathrm{E}\mathrm{n}\mathrm{e}\mathrm{r}\mathrm{g}\mathrm{y}=\mathrm{S}\mathrm{u}\mathrm{p}\mathrm{p}\mathrm{l}\mathrm{i}\mathrm{e}\mathrm{d}\mathrm{E}\mathrm{n}\mathrm{e}\mathrm{r}\mathrm{g}\mathrm{y}-\mathrm{D}\mathrm{i}\mathrm{s}\mathrm{s}\mathrm{i}\mathrm{p}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{d}\mathrm{E}\mathrm{n}\mathrm{e}\mathrm{r}\mathrm{g}\mathrm{y}\\ \mathcal{H}\left(t\right)-\mathcal{H}\left(0\right)={\int }_0^t{\mathbf{y}}^{\top }\left(\tau \right)\mathbf{u}\left(\tau \right)d\tau -{\int }_0^t{\dot{\mathit{q}}}^{\top }\left(\tau \right)\mathit{R}\dot{\mathit{q}}\left(\tau \right)d\tau \end{array}$$

(25)

Because $\mathcal{H}\left(t\right)\ge 0$ we deduce that

$$\mathcal{H}\left(t\right)=\mathcal{H}\left(0\right)+{\int }_0^t{\mathbf{y}}^{\top }\left(t\right)\mathbf{u}\left(t\right)dt-{\int }_0^t{\dot{\mathit{q}}}^{\top }\left(t\right)\mathit{R}\dot{\mathit{q}}\left(t\right)dt \ge 0$$

(26)

$\mathrm{t}\mathrm{h}\mathrm{e}\mathrm{r}\mathrm{e}\mathrm{f}\mathrm{o}\mathrm{r}\mathrm{e}{\int }_0^t{\mathbf{y}}^{\mathrm{\top }}\left(\tau \right)\mathbf{u}\left(\tau \right)d\tau \ge {\int }_0^t{\dot{\mathit{q}}}^{\mathrm{\top }}\left(\tau \right)\mathit{R}\dot{\mathit{q}}\left(\tau \right)d\tau -\mathcal{H}\left(0\right)$

$${\int }_0^t{\mathbf{y}}^{\top }\left(\tau \right)\mathbf{u}\left(\tau \right)d\tau \ge {\lambda }_{min}\left\{\mathit{R}\right\}{\int }_0^t{‖\dot{\mathit{q}}\left(\tau \right)‖}^{2}d\tau -\mathcal{H}\left(0\right)$$

(27)

This outcome proves that the system (IM) is passive, hence the $\mathrm{p}\mathrm{a}\mathrm{s}\mathrm{s}\mathrm{i}\mathrm{v}\mathrm{i}\mathrm{t}\mathrm{y}$ of the induction motor is proved. Also we note that the vector given by $\mathit{W}=\mathit{C}\left(\mathit{q},\dot{\mathit{q}}\right)\dot{\mathit{q}}$ contains the forces which do not produce work. As the IM system is an interconnection of two subsystems (electrical and mechanical see Figure 7) then it can be shown that each of the two subsystems are passive; so we arrive at a new approach to the structure of passivity based on the fact that the interconnection of two passive subsystems is passive.

We consider the total energy ${\mathcal{H}}_e$ of the electrical subsystem: ${\mathcal{H}}_e=\left[{\dot{\mathit{q}}}_e^{\top }{\mathit{D}}_e\left(q_m\right){\dot{\mathit{q}}}_e\right]/2$, then the time derivative of this ${\mathcal{H}}_e$ along the trajectories of the electrical system can be expressed as

$$\left\{\begin{array}{c}{\displaystyle \frac{d{\mathcal{H}}_e}{dt}}={\displaystyle \frac{1}{2}}{\dot{\mathit{q}}}_e^{\top }\left\{{\mathit{W}}_1{\dot{q}}_m\right\}{\dot{\mathit{q}}}_e+{\dot{\mathit{q}}}_e^{\top }{\mathit{D}}_e\left(q_m\right){\ddot{\mathit{q}}}_e \\ \begin{array}{c}={\displaystyle \frac{1}{2}}{\dot{\mathit{q}}}_e^{\top }\left\{{\mathit{W}}_1{\dot{q}}_m\right\}{\dot{\mathit{q}}}_e+{\dot{\mathit{q}}}_e^{\top }{\mathbf{u}}_e-{\dot{\mathit{q}}}_e^{\top }{\mathit{W}}_1{\dot{q}}_m{\dot{\mathit{q}}}_e-{\dot{\mathit{q}}}_e^{\top }{\mathit{R}}_e{\dot{\mathit{q}}}_e\\ \begin{array}{c}={\dot{\mathit{q}}}_e^{\top }{\mathbf{u}}_e-\mathcal{T}{\dot{q}}_m-{\dot{\mathit{q}}}_e^{\top }{\mathit{R}}_e{\dot{\mathit{q}}}_e\\ =\left[{\dot{\mathit{q}}}_e^{\top } \mathcal{T}\right]\left[\begin{array}{c}{\mathbf{u}}_e\\ -{\dot{q}}_m\end{array}\right]-{\dot{\mathit{q}}}_e^{\top }{\mathit{R}}_e{\dot{\mathit{q}}}_e\end{array} \end{array}\end{array} \right\}$$

(28)

So, the derivative of this Hamiltonian function is given by $d{\mathcal{H}}_e/dt=-{\dot{\mathit{q}}}_e^{\top }{\mathit{R}}_e{\dot{\mathit{q}}}_e+{\mathbf{y}}_1^{\top }{\mathbf{u}}_1$ by integrating this, we obtain

$${\mathcal{H}}_e\left(t\right)-{\mathcal{H}}_e\left(0\right)={\int }_0^t{\mathbf{y}}_1^{\top }\left(\tau \right){\mathbf{u}}_1\left(\tau \right)d\tau -{\int }_0^t{\dot{\mathit{q}}}_e^{\top }\left(\tau \right){\mathit{R}}_e{\dot{\mathit{q}}}_e\left(\tau \right)d\tau$$

(29)

Because ${\mathcal{H}}_e\left(t\right)\ge 0$ we deduce that

$${\mathcal{H}}_e\left(0\right)+{\int }_0^t{\mathbf{y}}_1^{\top }\left(\tau \right){\mathbf{u}}_1\left(\tau \right)d\tau -{\int }_0^t{\dot{\mathit{q}}}_e^{\top }\left(\tau \right){\mathit{R}}_e{\dot{\mathit{q}}}_e\left(\tau \right)d\tau \ge 0$$

(30)

This inequality gives us:

$${\int }_0^t{\mathbf{y}}_1^{\top }\left(\tau \right){\mathbf{u}}_1\left(\tau \right)d\tau \ge {\lambda }_{min}\left\{{\mathit{R}}_e\right\}{\int }_0^t{‖{\dot{\mathit{q}}}_e\left(t\right)‖}^{2}dt-{\mathcal{H}}_e\left(0\right)$$

(31)

which proves the passivity of the electrical subsystem. On the other hand the mechanical subsystem is positive real for all $D_m$ and $R_m$ which is also passive.

PROPOSITION: 1 Consider the asynchronous motor given by equation(6) with outputs $\mathcal{T}$ and $‖{\lambda }_r‖=\beta$ to be controlled, and assumptions A1 to A3. Let the control law be defined as,

$$\mathbf{u}\left(t\right)=L_s{\ddot{\mathit{q}}}_s^{\mathrm{r}\mathrm{e}\mathrm{f}}+L_m{e}^{\mathcal{J}p{q}_m}{\ddot{\mathit{q}}}_r^{\mathrm{r}\mathrm{e}\mathrm{f}}+p{L}_m\mathcal{J}{e}^{\mathcal{J}p{q}_m}{\dot{\widehat{\mathit{q}}}}_r{\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}+R_s{\dot{\mathit{q}}}_s^{\mathrm{r}\mathrm{e}\mathrm{f}} -K_1\left({\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}\right){\dot{\tilde{\mathit{q}}}}_s$$

(32)

The superscript $\mathrm{r}\mathrm{e}\mathrm{f}$ is the desired variable, with

${\dot{\mathit{q}}}_r^{\mathrm{r}\mathrm{e}\mathrm{f}}=-e^{\mathcal{J}{\rho }_d}\left[\begin{array}{c}\begin{array}{c}\dot{\beta }/R_r\\ \end{array}\\ {\mathcal{T}}_d/\left(p\beta \right)\end{array}\right]; {\dot{\mathit{q}}}_s^{\mathrm{r}\mathrm{e}\mathrm{f}}={\displaystyle \frac{1}{L_m}}{e}^{\mathcal{J}(p{q}_m+{\rho }_d)}\left[\begin{array}{c}\begin{array}{c}\beta +\left(L_r/R_r\right)\dot{\beta }\\ \end{array}\\ L_r{\mathcal{T}}_d/\left(p\beta \right)\end{array}\right]; {\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}=\left[\begin{array}{c}{\dot{\mathit{q}}}_e^{\mathrm{r}\mathrm{e}\mathrm{f}}\\ {\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}\end{array}\right]; {\dot{\mathit{q}}}_e^{\mathrm{r}\mathrm{e}\mathrm{f}}=\left[\begin{array}{c}{\dot{\mathit{q}}}_s^{\mathrm{r}\mathrm{e}\mathrm{f}}\\ {\dot{\mathit{q}}}_r^{\mathrm{r}\mathrm{e}\mathrm{f}}\end{array}\right]$and with the following controller dynamics ${\dot{\rho }}_d\left(t\right)=\left(R_r/p\beta \right){\mathcal{T}}_d\left(t\right), {\rho }_d\left(0\right)=0$ and ${\ddot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}\left(t\right)=\left[-p{L}_m{\dot{\widehat{\mathit{q}}}}_r^{\top }\mathcal{J}{e}^{-\mathcal{J}p{q}_m}{\dot{\mathit{q}}}_s^{\mathrm{r}\mathrm{e}\mathrm{f}}-R_m{\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}-{\widehat{\mathcal{T}}}_L+K_2\left({\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right){\dot{\tilde{q}}}_m\right]/D_m, {\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}\left(0\right)={\dot{q}}_m\left(0\right)$.

The gains $K_1\left({\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)$ and $K_2\left({\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)$ are given as

$$K_1\left({\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}\right):={\displaystyle \frac{p^2{L}_m^2}{4{\epsilon }_1}}{\left({\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)}^2+k_1; K_2\left({\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right):={\displaystyle \frac{p^2{L}_m^2}{4{\epsilon }_1}}\left[{‖{\dot{\mathit{q}}}_e^{\mathrm{r}\mathrm{e}\mathrm{f}}‖}^2+{\left({\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)}^2\right]+k_2;$$

(33)

with $0<{\epsilon }_1<R_r/2, k_1\ge 0, k_2>0$. while the state estimator and the load adaptation law are

$${\mathit{D}}_e\left(q_m\right){\ddot{\widehat{\mathit{q}}}}_e+{\mathit{W}}_1\left(q_m\right){\dot{q}}_m{\dot{\widehat{\mathit{q}}}}_e+{\mathit{R}}_e{\dot{\widehat{\mathit{q}}}}_e={\mathit{M}}_e\mathbf{u}-\mathit{L}\left({\dot{q}}_m\right){\dot{\mathit{e}}}_e$$

(34)

$${\dot{\widehat{\mathcal{T}}}}_L=-{\gamma }_L{\dot{\tilde{q}}}_m, {\gamma }_L>0$$

(35)

where ${\dot{\mathit{e}}}_e={\dot{\widehat{\mathit{q}}}}_e-{\dot{\mathit{q}}}_e$, ${\dot{\tilde{q}}}_m={\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}-{\dot{q}}_m$ and

$$\mathit{L}\left({\dot{q}}_m\right)=\left[\begin{array}{cc}2\theta {\mathit{I}}_2& 0_2\\ {\theta }^2\mathbf{h}\left({\dot{q}}_m\right)& 0_2\end{array}\right]; \mathrm{a}\mathrm{n}\mathrm{d}\mathbf{h}\left({\dot{q}}_m\right)={\mathbf{F}}_1^{-1}\left({\dot{q}}_m\right)={\displaystyle \frac{\left\{\frac{K}{T_r}\mathit{I}+Kp{\dot{q}}_m\mathcal{J}\right\}}{{\left(K/T_r\right)}^2+{\left(Kp{\dot{q}}_m\right)}^2}}\in {\mathbb{R}}^{2\times 2}$$

Under these conditions, the closed-loop system guarantees global tracking of torque and flux norm, while ensuring all signals remain uniformly bounded. ∎

Proof: 

Because the control law utilizes the estimated states rather than the actual ones, as opposed to the ideal scenario, the error dynamics in this case are:

$$\mathit{D}\left(\mathit{q}\right)\ddot{\tilde{\mathit{q}}}\left(t\right)+\left[\mathit{C}\left(\mathit{q},\dot{\mathit{q}}\right)+\mathit{R}+\mathit{K}\left({\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)\right]\dot{\tilde{\mathit{q}}}\left(t\right)=\mathit{S}\left(q_m,{\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right){\dot{\mathit{e}}}_e\left(t\right)+\tilde{\mathit{\xi }}\left(t\right)$$

(36)

where $\mathit{K}\left({\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)=\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left\{K_1\left({\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}\right){\mathit{I}}_2,\mathrm{0,0},K_2\left({\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)\right\}$, $\tilde{\mathit{\xi }}:={\left[0_{1\times 4},{\tilde{\mathcal{T}}}_L\right]}^{\top }={\left[\mathrm{0,0},\mathrm{0,0},({\widehat{\mathcal{T}}}_L-{\mathcal{T}}_L)\right]}^{\top }$and

$$\mathit{S}\left(q_m,{\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)=\left[\begin{array}{cc}{0}_2 & p{L}_m\mathcal{J}{e}^{\mathcal{J}p{q}_m}{\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}\\ 0_{1\times 2}& 0_2\\ 0_2 & -p{L}_m{\left({\dot{\mathit{q}}}_s^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)}^{\top }\mathcal{J}{e}^{\mathcal{J}p{q}_m}\end{array}\right]\in {\mathbb{R}}^{5\times 4}; \dot{\tilde{\mathit{q}}}\left(t\right)=\left[\begin{array}{c}\begin{array}{c}{\dot{\mathit{q}}}_s\left(t\right)\\ {\dot{\mathit{q}}}_r\left(t\right)\end{array}\\ {\dot{q}}_m\left(t\right)\end{array}\right]-\left[\begin{array}{c}\begin{array}{c}{\dot{\mathit{q}}}_s^{\mathrm{r}\mathrm{e}\mathrm{f}}\left(t\right)\\ {\dot{\mathit{q}}}_r^{\mathrm{r}\mathrm{e}\mathrm{f}}\left(t\right)\end{array}\\ {\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}\left(t\right)\end{array}\right]\in {\mathbb{R}}^{5\times 1}$$

Alternatively, based on equations (34) and (1), the observation error ${\dot{\mathit{e}}}_e$​ satisfies the following equation: ${\mathit{D}}_e\left(q_m\right){\ddot{\mathit{e}}}_e+\left[{\mathit{W}}_1\left(q_m\right){\dot{q}}_m+\mathit{L}\left({\dot{q}}_m\right)\right]{\dot{\mathit{e}}}_e+{\mathit{R}}_e{\dot{\mathit{e}}}_e=0$. So, consider the following Lyapunov function candidate

$$V={\displaystyle \frac{1}{2}}{\dot{\tilde{\mathit{q}}}}^{\top }\mathit{D}\left(\mathit{q}\right)\dot{\tilde{\mathit{q}}}+{\displaystyle \frac{1}{2}}{\dot{\mathit{e}}}_e^{\top }{\mathit{D}}_e\left(q_m\right){\dot{\mathit{e}}}_e+{\displaystyle \frac{1}{2{\gamma }_L}}{\tilde{\mathcal{T}}}_L^2$$

(37)

Knowing that $\dot{\mathit{D}}\left(\mathit{q}\right)-2\mathit{C}(\mathit{q},\dot{\mathit{q}})$ is skew-symmetric and ${\dot{\mathit{D}}}_e\left(q_m\right)-2\left[{\mathit{W}}_1\left(q_m\right){\dot{q}}_m+\mathit{L}\left({\dot{q}}_m\right)\right]$ which yields the derivative

$$\dot{V}=-{\dot{\tilde{\mathit{q}}}}^{\top }\left[\mathit{R}+\mathit{K}\left({\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)\right]\dot{\tilde{\mathit{q}}}+{\dot{\tilde{\mathit{q}}}}^{\top }\mathit{S}\left(q_m,{\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right){\dot{\mathit{e}}}_e-{\dot{\mathit{e}}}_e^{\top }{\mathit{R}}_e{\dot{\mathit{e}}}_e+{\displaystyle \frac{1}{{\gamma }_L}}{\tilde{\mathcal{T}}}_L{\dot{\tilde{\mathcal{T}}}}_L+{\dot{\tilde{\mathit{q}}}}^{\top }\tilde{\mathit{\xi }}$$

(38)

Use of equation (35) in the equation above and let we define$z:={\left[{\dot{\tilde{\mathit{q}}}}^{\top },{\dot{\mathit{e}}}_e^{\top }\right]}^{\top }$, which results in the next quadratic function $\dot{V}=-{\mathit{z}}^{\top }\mathit{M}\mathit{z}$ with

$$\mathit{M}=\left[\begin{array}{cc}\mathit{R}+\mathit{K}\left({\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)& 0\\ 0& {\mathit{R}}_e\end{array}\right]-{\displaystyle \frac{1}{2}}\left[\begin{array}{cc}0& \mathit{S}\left(q_m,{\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)\\ {\mathit{S}}^{\top }\left(q_m,{\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)& 0\end{array}\right]$$

Verifying that $\mathit{K}\left({\dot{\mathit{q}}}^{\mathrm{r}\mathrm{e}\mathrm{f}}\right)$ guarantee the strict positive definiteness of$\mathit{M}$, i.e. $\mathit{M}\ge \delta {\mathit{I}}_9>0$∎

To enable full sensorless operation prior to incorporating intelligent learning, the control framework is first augmented with a conventional PI-based speed estimator. This allows rotor speed reconstruction from electrical measurements alone, eliminating the need for mechanical sensors while preserving the structure and stability of the original passivity-based control loop.

The estimation flux is used to drive a PI-like update law for ${\widehat{q}}_m$. A standard structure is:

$${\widehat{q}}_m=K_p\epsilon \left(t\right)+K_i\int \epsilon \left(t\right)dt \mathrm{w}\mathrm{i}\mathrm{t}\mathrm{h} \epsilon \left(t\right)={\dot{\mathit{q}}}_s^{\top }\mathcal{J}{\mathit{\lambda }}_r=\left[{\widehat{\lambda }}_{r\alpha }\left(t\right){\dot{q}}_{s\beta }\left(t\right)-{\widehat{\lambda }}_{r\beta }\left(t\right){\dot{q}}_{s\alpha }\left(t\right)\right]$$

(39)

$\epsilon \left(t\right)={\dot{\mathit{q}}}_s^{\top }\mathcal{J}{\mathit{\lambda }}_r$ is the error between the real and estimated rotor fluxes and currents.

The overall architecture of the non-intelligent sensorless adaptive passivity-based controller is illustrated in Figure 8.

Although the previous controller ensures tracking and robustness via passivity control, its performance may degrade under unmodeled dynamics or uncertainties. To enhance adaptability, a recurrent neural network is integrated to learn residual dynamics in real time, yielding an intelligent adaptive control.

4.2. The Proposed Intelligent Adaptive Passivityt Based Controller Design

The proposed intelligent adaptive passivity-based controller integrates neural learning for enhanced sensorless performance. As shown in Figure 9, a recurrent neural network estimates both the load torque ${\widehat{\mathcal{T}}}_L$​ and rotor speed ${\dot{\widehat{q}}}_m$​, enabling full sensorless operation. These estimates drive a passivity-based control loop to generate stator voltage commands, ensuring robustness in both steady-state and dynamic conditions.

The proposed scheme employs a recurrent neural network to estimate the load torque ${\mathcal{T}}_L^{\star }$​ and the rotor speed ${\widehat{\dot{q}}}_m$. By capturing temporal dependencies, the RNN effectively models the motor’s nonlinear dynamics, outperforming conventional feedforward networks. Let $e_{\omega }\left(k\right)={\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}\left(k\right)-{\widehat{\dot{q}}}_m\left(k\right)$ and $e_{\mathcal{T}}\left(k\right)={\widehat{\lambda }}_{r\alpha }\left(k\right){\dot{q}}_{s\beta }\left(k\right)-{\widehat{\lambda }}_{r\beta }\left(k\right){\dot{q}}_{s\alpha }\left(k\right)$ denote the speed and the electromagnetic torque errors, at discrete time step $k$. The internal structure of the RNN is governed by the following equations [25,26]:

1.

Hidden State Update: $\mathbb{h}\left(k\right)=\mathsf{\Psi }({\mathcal{W}}_{\mathrm{i}\mathrm{n}}\mathit{}e\left(k\right)+{\mathcal{W}}_{\mathrm{r}\mathrm{e}\mathrm{c}}\mathbb{h}\left(k-1\right)+{\mathbf{b}}_h\mathit{})$ where:

• $\mathbb{h}\left(k\right)\in {\mathbb{R}}^n$: hidden state vector, ${\mathbf{b}}_h\in {\mathbb{R}}^n$: bias vector,
• ${\mathcal{W}}_{\mathrm{i}\mathrm{n}}\in {\mathbb{R}}^{n\times 1}$: input weight matrix, ${\mathcal{W}}_{\mathrm{r}\mathrm{e}\mathrm{c}}\in {\mathbb{R}}^{n\times n}$: recurrent weight matrix,
• $\mathsf{\Psi }\left(\cdot \right)$: element-wise activation function (e.g., tanh or ReLU).

2.

The Output Estimation: $y_{\mathrm{n}\mathrm{e}\mathrm{t}}^{\star }\left(k\right)={\mathcal{W}}_{\mathrm{o}\mathrm{u}\mathrm{t}}\mathbb{h}\left(k\right)+b_{\mathrm{o}}$ where:

• ${\mathcal{W}}_{\mathrm{o}\mathrm{u}\mathrm{t}}\in {\mathbb{R}}^{1\times n}$: output weight vector, $b_{\mathrm{o}}\mathbb{}\in \mathbb{R}$: output bias.

The complete mapping from input error to the nedded estimate is compactly expressed as one function: $y_{\mathrm{n}\mathrm{e}\mathrm{t}}^{\star }\left(k\right)={\mathbf{f}}_{\mathrm{R}\mathrm{N}\mathrm{N}}\left(e\left(k\right),\mathbb{h}\left(k-1\right);\mathit{\theta }\right)$ where $\mathit{\theta }=\left\{{\mathcal{W}}_{\mathrm{i}\mathrm{n}},{\mathcal{W}}_{\mathrm{r}\mathrm{e}\mathrm{c}},{\mathcal{W}}_{\mathrm{o}\mathrm{u}\mathrm{t}},{\mathbf{b}}_h,b_{\mathrm{o}}\right\}$ is the set of all trainable parameters [30]. This RNN-based estimator enables the controller to adaptively generate commands that account for the system's nonlinearities, time delays, and dynamic uncertainties, thus enhancing overall control performance. The RNN parameters $\mathit{\theta }$ are optimized to minimize the mean squared error loss: $\mathcal{L}=\left[{⅀}_{k=1}^N{\left\{y_{\mathrm{n}\mathrm{e}\mathrm{t}}^{\star }\left(k\right)-y_{\mathrm{n}\mathrm{e}\mathrm{t}}\left(k\right)\right\}}^2\right]/N$. First we compute hidden states and output estimates sequentially: $\mathbb{h}\left(k\right)=\mathsf{\Psi }({\mathcal{W}}_{\mathrm{i}\mathrm{n}}\mathit{}e\left(k\right)+{\mathcal{W}}_{\mathrm{r}\mathrm{e}\mathrm{c}}\mathit{}\mathbb{h}\left(k-1\right)+{\mathbf{b}}_h\mathit{})$ and $y_{\mathrm{n}\mathrm{e}\mathrm{t}}^{\star }\left(k\right)={\mathcal{W}}_{\mathrm{o}\mathrm{u}\mathrm{t}}\mathbb{h}\left(k\right)+b_{\mathrm{o}}$ then compute gradients $\mathsf{\Delta }=\partial \mathcal{L}/\partial \mathit{\theta }$​ by unrolling the network through time and applying the chain rule. Finally, update parameters using gradient descent or an advanced optimizer (e.g., Adam): $\mathit{\theta }=\mathit{\theta }+\eta \mathsf{\Delta }$ where $\eta$ is the learning rate. This iterative process continues until convergence, enabling the RNN to learn the dynamic mapping from error sequences to the output commands effectively [31]. If we define $\delta \left(k\right)=\partial \mathcal{L}/\partial y_{\mathrm{n}\mathrm{e}\mathrm{t}}^{\star }=2\left[y_{\mathrm{n}\mathrm{e}\mathrm{t}}^{\star }\left(k\right)-y_{\mathrm{n}\mathrm{e}\mathrm{t}}\left(k\right)\right]/N$ then the gradients for output layer are given by ${\mathbf{g}}_{\mathrm{o}\mathrm{u}\mathrm{t}}=\partial \mathcal{L}/\partial {\mathcal{W}}_{\mathrm{o}\mathrm{u}\mathrm{t}}={⅀}_{k=1}^N\delta \left(k\right){\mathbb{h}}^{\mathit{\top }}\left(k\right)$ and ${\mathbf{g}}_{b_{\mathrm{o}}}=\partial \mathcal{L}/\partial b_{\mathrm{o}}={⅀}_{k=1}^N\delta \left(k\right)$. Also, if we back-propagate through time (BPTT) into hidden states $\mathit{ϵ}\left(k\right)=\partial \mathcal{L}/\partial \mathbb{h}\left(k\right)\in {\mathbb{R}}^n$ then we compute recursively (i.e. $k=N:1$) as $\mathit{ϵ}\left(k\right)={\mathcal{W}}_{\mathrm{o}\mathrm{u}\mathrm{t}}^{\mathit{\top }}\delta \left(k\right)+{\mathcal{W}}_{\mathrm{r}\mathrm{e}\mathrm{c}}^{⟙}\left[\mathit{ϵ}\left(k+1\right)⨀{\mathsf{\Psi }}^{\prime }\left(\mathit{a}\left(k+1\right)\right)\right]$ where $\mathit{a}\left(k\right)={\mathcal{W}}_{\mathrm{i}\mathrm{n}}e\left(k\right)+{\mathcal{W}}_{\mathrm{r}\mathrm{e}\mathrm{c}}\mathbb{h}\left(k-1\right)+{\mathbf{b}}_h$, ${\mathsf{\Psi }}^{\prime }\left(\mathit{a}\right)$ is the element-wise derivative of activation function and $⨀$ is the element-wise multiplication. Now use $\mathit{ϵ}\left(k\right)$ to compute the following gradients: ${\mathbf{g}}_{\mathrm{r}\mathrm{e}\mathrm{c}}=\partial \mathcal{L}/\partial {\mathcal{W}}_{\mathrm{r}\mathrm{e}\mathrm{c}}={⅀}_{k=1}^N\mathbf{s}\left(k\right).{\mathbb{h}}^{\mathit{\top }}\left(k-1\right) \in {\mathbb{R}}^{n\times n}$, ${\mathbf{g}}_{\mathrm{i}\mathrm{n}}=\partial \mathcal{L}/\partial {\mathcal{W}}_{\mathrm{i}\mathrm{n}}={⅀}_{k=1}^N\mathbf{s}\left(k\right).e^{\mathit{\top }}\left(k\right)\in {\mathbb{R}}^{n\times 1}$ and ${\mathbf{g}}_{{\mathbf{b}}_h}=\partial \mathcal{L}/\partial {\mathbf{b}}_h={⅀}_{k=1}^N\mathbf{s}\left(k\right).\in {\mathbb{R}}^n$ with $\mathbf{s}\left(k\right)=\left[\mathit{ϵ}\left(k\right)⨀{\mathsf{\Psi }}^{\prime }\left(\mathit{a}\left(k\right)\right)\right]$. Finally, if we let $\eta$ be the learning rate, then we can updtae: ${\mathcal{W}}_{\mathrm{o}\mathrm{u}\mathrm{t}}\leftarrow {\mathcal{W}}_{\mathrm{o}\mathrm{u}\mathrm{t}}-\eta {\mathbf{g}}_{\mathrm{o}\mathrm{u}\mathrm{t}}$, ${\mathcal{W}}_{\mathrm{r}\mathrm{e}\mathrm{c}}\leftarrow {\mathcal{W}}_{\mathrm{r}\mathrm{e}\mathrm{c}}-\eta {\mathbf{g}}_{\mathrm{r}\mathrm{e}\mathrm{c}}$ etc. All parameters are updated this way. The structure of this learning procedure is given in the flowchart shown in Figure 10

The RNN parameter update steps using BPTT are described in the following algorithm

5. Experimental Results and Discussion

The neural passive controller has been tested on a test setup (Figure 11) which consists of:

• Three-phase induction motor with rated values:

• Synchronous machine with powder brakefor for loading the IM.

• Electronic power converter (without any control system):

• 3$\varphi$ (three-phase) diode rectifier and

• VSI composed of three IGBT modules

• Electronic card for monitoring the stator phase with

• Instantaneous voltage sensors (LEM LV 25-P) and

• Instantaneous current sensors (LEM LA 55-P)

• Voltage sensor for monitoring the instantaneous value of

• The dc-link voltage. (model LEM CV3-1000)

• Incremental encoder only for comparison measurements.

• Model RS 256-499, 2500 pulses per round, (for validation-comparison, not for control)

• $\mathrm{d}\mathrm{S}\mathrm{P}\mathrm{A}\mathrm{C}\mathrm{E}$ card ($\mathrm{D}\mathrm{S}1104$) with a PowerPC $604\mathrm{e}$ at $400\mathrm{M}\mathrm{H}\mathrm{z}$ and floating-point digital signal processor $\mathrm{D}\mathrm{S}\mathrm{P}\mathrm{⎯}\mathrm{T}\mathrm{M}\mathrm{S}320\mathrm{F}240$. During the real-time operation of the control, the supervision/capturing of the important data can be done by the $\mathrm{C}\mathrm{O}\mathrm{N}\mathrm{T}\mathrm{R}\mathrm{O}\mathrm{L}\mathrm{D}\mathrm{E}\mathrm{S}\mathrm{K}$ software provided with $\mathrm{D}\mathrm{S}\mathrm{P}$ board.

Figure 11. Experimental setup for the proposed neural-PBC daptive control of the IM.

Figure 11. Experimental setup for the proposed neural-PBC daptive control of the IM.

Simulations in MATLAB verified the proposed control scheme using a rated flux reference and $500\mathrm{V}$ DC-link voltage. The feedback signal reconstruction was confirmed in simulation, and experimental results further validated the effectiveness-robustness of the sensorless control strategy. The electrical and mechanical parameters of the induction motor are listed in Table 2.

The “Control/Command” part is based on the $\mathrm{D}\mathrm{S}1104$ R&D Controller Board developed by the German company dSPACE and housed in a computer. This control card is made up of two processors. The master processor manages the application while the slave processor, a DSP (“Digital Signal Processor”) from TEXAS INSTRUMENT (type $\mathrm{T}\mathrm{M}\mathrm{S}320\mathrm{F}240$), generates the $\mathrm{P}\mathrm{W}\mathrm{M}$ (Pulse Width Modulation) control signals in 0/5 V TTL logic. This constitutes the “hardware” part of $\mathrm{d}\mathrm{S}\mathrm{P}\mathrm{A}\mathrm{C}\mathrm{E}$.. The "software" part consists of two software programs. The first, $\mathrm{M}\mathrm{A}\mathrm{T}\mathrm{L}\mathrm{A}\mathrm{B}/\mathrm{S}\mathrm{i}\mathrm{m}\mathrm{u}\mathrm{l}\mathrm{i}\mathrm{n}\mathrm{k}$, allows easy programming of the real-time application under Simulink by using specific blocks (Belonging to the "Real Time Interface ($\mathrm{R}\mathrm{T}\mathrm{I}$) toolbox") allowing the configuration of the inputs/outputs of the $\mathrm{D}\mathrm{S}1104$ card. The second software, $\mathrm{C}\mathrm{o}\mathrm{n}\mathrm{t}\mathrm{r}\mathrm{o}\mathrm{l}\mathrm{D}\mathrm{e}\mathrm{s}\mathrm{k}$, allows the program code to be loaded onto the card (written in graphical form in Simulink, compiled and transformed into C code), to create a complete experimental environment and in particular a graphical interface for controlling the real-time process, to process the data and save it in a format compatible with $\mathrm{M}\mathrm{A}\mathrm{T}\mathrm{L}\mathrm{A}\mathrm{B}$ (for later processing) or even to monitor in real time the evolution of the measured or calculated data using graphical or digital displays. The exchange of information between the two parts described is carried out via an external connection box (Connector Panel $\mathrm{C}\mathrm{P}1104$ from $\mathrm{d}\mathrm{S}\mathrm{P}\mathrm{A}\mathrm{C}\mathrm{E}$) connected to the card via a shielded cable and receiving the analog signals via $\mathrm{B}\mathrm{N}\mathrm{C}$ connectors, an interface for conditioning the $\mathrm{P}\mathrm{W}\mathrm{M}$ control signals and any error signals returned by the Semikron converter and a measurement environment consisting of various sensors. The signal conditioning interface converts the signals from $\mathrm{T}\mathrm{T}\mathrm{L}$ 0/5 V logic to $\mathrm{C}\mathrm{M}\mathrm{O}\mathrm{S}$ 0/15 V logic and vice versa. This modification is essential because the $\mathrm{D}\mathrm{S}1104$ control card works with $\mathrm{T}\mathrm{T}\mathrm{L}$ 0/5 V logic signals, while these must be in $\mathrm{C}\mathrm{M}\mathrm{O}\mathrm{S}$ 0/15 V logic for the voltage inverter. The measurement environment consists of LEM type LA25TP sensors (closed-loop current sensors using the Hall effect) for current measurements, LEM type LV100-500 sensors (closed-loop voltage sensors using the Hall effect) for voltage measurements and an incremental encoder to measure the motor rotation speed. Finally, the measured analog currents are intended for digital processing and must therefore be sampled. Therefore, in order to avoid any spectral aliasing phenomenon, it is necessary to insert a guard filter (with an estimated cut-off frequency of 500 Hz, one order of magnitude above the fundamental frequency of 50 Hz) between each sensor and the analog-to-digital converter.

Analyzing the operation of an IM supplied by a voltage inverter is complex due to its complicated nonlinear nature. Moreover, neural (or neuro-fuzzy) adaptive control of an IM requires solid theoretical foundations in areas such as electrical machines, power electronics, and control systems. Using the $\mathrm{d}\mathrm{S}\mathrm{P}\mathrm{A}\mathrm{C}\mathrm{E}$-instrumented platform, the proposed control strategy is first implemented in Simulink, then deployed in real-time via $\mathrm{d}\mathrm{S}\mathrm{P}\mathrm{A}\mathrm{C}\mathrm{E}$. The $\mathrm{P}\mathrm{W}\mathrm{M}$-driven inverter operates at a switching frequency of $9.5\mathrm{k}\mathrm{H}\mathrm{z}$. A detailed view of the experimental setup is provided in Figure 12.

5.1. Study of the Proposed Nonlinear Intelligent Observer

The performance of the proposed nonlinear observer is first evaluated in open-loop and then integrated into a closed-loop configuration with the nonlinear controller of the induction motor. To align the observer dynamics with those of the plant, a tuning parameter $5\le \theta \le 30$ is used. The trajectories of the motor and observer poles, are experimentally analyzed under this setting. Figure 13 illustrates the rotor flux observation error obtained at various operating speeds (0.0125, 3.4, and 150 rad/s), with $\theta =30, \mathrm{a}\mathrm{n}\mathrm{d} \theta =30$, demonstrating the accuracy and robustness of the observer. The given architecture is tested under different perturbation scenarios, including parameter variations and unmodeled dynamics, to validate the strength and reliability of the proposed intelligent strategy.

The convergence characteristics of the proposed nonlinear observer are significantly affected by the tuning parameter $\theta$, which shapes the placement of its eigenvalues relative to those of the motor. As shown in the upper subplots of Figure 14, increasing $\theta$ shifts the observer poles (blue) to follow the motor poles (red) more closely across the operating range $\mathsf{\Omega }=0$ to 400 rad/s. For $\theta =5$, the observer exhibits a moderately damped response, while $\theta =30$ yields faster convergence with reduced transients. The lower subplots of Figure 14 illustrate the corresponding rotor flux observation errors under three representative speeds (0.0125, 3.4, and 150 rad/s). In both cases, the observer demonstrates strong convergence properties, but higher $\theta$ improves transient rejection and tracking precision, especially at low and nominal speeds.

To evaluate the sensitivity of the proposed nonlinear observer to rotor resistance variations, Figure 14 presents experimental results under three resistance values: $3.047\mathrm{O}\mathrm{h}\mathrm{m}$ (nominal), $6.4571\mathrm{O}\mathrm{h}\mathrm{m}$ (+50%), and $8.094 \mathrm{O}\mathrm{h}\mathrm{m}$ (+100%). The observer is initialized with rotor flux components set to ${\lambda }_{r\alpha 0}={\lambda }_{r\beta 0}=0.5 \mathrm{W}\mathrm{b}$.

Figure 14 demonstrates that, even with $R_r$ increased from $3.047\mathrm{O}\mathrm{h}\mathrm{m}$ (nominal) to $6.457\mathrm{}\mathrm{O}\mathrm{h}\mathrm{m}$ (+50%) and $8.094\mathrm{}\mathrm{O}\mathrm{h}\mathrm{m}$ (+100%), the nonlinear observer ensures rapid flux convergence. At both nominal (150 rad/s) and low speed (60 rad/s), transients settle within 0.6 s, showing strong disturbance rejection. A small steady-state offset remains (about 0.15 Wb at +50% and 0.30 Wb at +100%) indicating sensitivity to parameter mismatch. To address this, an RBFNN-based observer is introduced (Figure 14), incorporating an adaptive term that continuously compensates for rotor resistance errors and effectively eliminates steady-state flux bias under large uncertainties.

Figure 15 highlights the impact of the RBFNN augmentation on the observer’s transient performance under large rotor resistance variations. In the baseline observer, the flux error exhibits pronounced oscillations and peak deviations nearing 0.8 Wb for a +100 % resistance change, with settling oscillations persisting beyond 0.6s. With the RBFNN-enhanced observer, the peak error is reduced by over 50 %, and oscillations are significantly damped—settling in under 0.4 s across all resistance cases. This clear reduction in transient amplitude and faster damping confirms that the online neural adaptation effectively compensates for parameter mismatches, yielding smoother and more responsive flux estimation throughout the speed range.

Mismatch.

5.2. Comprehensive Evaluation of the Proposed Nonlinear Intelligent Passivity-Based Control

This section presents an overall assessment of the proposed nonlinear intelligent passivity-based control, which integrates a neural PBC controller, an RNN-based adaptive observer, and a sensorless speed estimator. The complete system is tested in closed-loop with the induction motor under various conditions, including load disturbances, speed reversals, and parameter variations. The results highlight the dynamic performance of key signals such as torque, rotor flux, speed, and control voltages, confirming the robustness, adaptability, and precision of the proposed scheme. Figure 16 presents these responses under all test scenarios.

The performance of the proposed nonlinear intelligent PBC strategy is illustrated in Figure 16, which presents both control and estimation results under time-varying speed and load torque conditions. The top-left subplot shows that the controlled electromagnetic torque ${\mathcal{T}}_d\left(t\right)$ tracks the reference smoothly with minimal steady-state error. The top-right plot confirms accurate load torque estimation ${\mathcal{T}}_L\left(t\right)$, closely matching the timing and magnitude of applied disturbances. The middle-left subplot shows that both estimated and measured speeds follow the commanded reference ${\dot{q}}_m^{\mathrm{r}\mathrm{e}\mathrm{f}}\left(t\right)$ throughout acceleration, cruising, deceleration, and re-acceleration phases. The speed error, shown bottom-left, remains within ±8 rad/s, demonstrating the precision of the RNN-based speed estimator. The control voltages $V_a$ and $V_b$, in the middle-right plot, reach ±250 V with consistent modulation behavior. The bottom-right plot shows torque error bounded within ±0.5 Nm, confirming high tracking accuracy and effective disturbance rejection. These results validate the proposed architecture’s capability for precise tracking, reliable estimation, and smooth control under sensorless operation.

Figure 17 was performed on a sensorless induction motor using the proposed PBC strategy for torque and flux control. The system faced load variations, noise, and rotor resistance mismatch, with all estimations handled observer-based. Control ran in real time with PWM injection, and key signals were precisely recorded for evaluation. The applied torque includes direction reversals to assess robustness.

The experimental results clearly confirm the effectiveness of the proposed nonlinear intelligent passivity-based control strategy. The motor tracks the torque reference with high accuracy, including during rapid transients and full torque reversal, maintaining a peak error within ±0.2 Nm. The estimated load torque closely matches the real applied load, capturing both magnitude and timing under dynamic conditions. The speed response aligns well with the reference, and the speed estimation remains precise even during acceleration and sharp changes. The flux estimation error is well confined within ±0.1 Wb, indicating strong observer performance. Input voltages exhibit clean modulation without saturation or chattering, demonstrating the quality of PWM-based voltage control. Despite load disturbances, parameter uncertainties, and full sensorless operation, the controller preserves stability, accurate tracking, and robustness, confirming its suitability for advanced industrial applications.

5.3. Comparative Study and Performance Evaluation

To evaluate the performance of the proposed control, we use standard electrical engineering (EE) control metrics based on the response to torque amplitude or flux amplitude (i.e. |final − initial|):

• $t_{2\%}$ is the time at which the response signal has covered 2% of the ramp amplitude.

• $t_{95\%}$ is the time after which the response remains at less than 5% of the target value.

• $D_{\%}$ Peak overshoot relative to target (% of ramp amplitude).

• $E_{\mathrm{s}\mathrm{s}}$ (steady-state error settling) the error once the $t_{95\%}$ has been reached.

• $E_{\mathrm{f}\mathrm{o}\mathrm{l}}$ Error when reference reaches 50% of ramp.

• ${∆\mathcal{T}}_{\mathrm{m}\mathrm{a}\mathrm{x}}$ Max torque deviation during speed ramp.

And in order analyze the performance at global scope, we employ standard machine learning (ML) metrics: mean absolute error (MAE), symmetric mean absolute percentage error (SMAPE), and coefficient of determination ($R^2$). These metrics are defined as: $\mathrm{M}\mathrm{A}\mathrm{E}=\left[{⅀}_{t=1}^T\left|y_t-{\widehat{y}}_t\right|\right]/T$; $\mathrm{S}\mathrm{M}\mathrm{A}\mathrm{P}\mathrm{E}=100.{⅀}_{t=1}^T\left[\left|y_t-{\widehat{y}}_t\right|/\left(\left|y_t\right|+\left|{\widehat{y}}_t\right|\right)\right]/T$ and $R^2=1-{⅀}_{t=1}^T{\left[{\overline{y}}_t-{\widehat{y}}_t\right]}^2/{⅀}_{t=1}^T{\left[{\overline{y}}_t-y_t\right]}^2$ where $y_t$ is the ground truth, ${\widehat{y}}_t$ is the predicted output of the model at time $t$, and $T$ is the total experiment duration. ${\overline{y}}_t$ denotes the mean of $y_t$. These metrics jointly characterize accuracy, relative error, and explained variance, enabling a comprehensive assessment of model fidelity.

To comprehensively assess the effectiveness of the proposed adaptive neuro-PBC sensorless control scheme under realistic and demanding operating conditions, two benchmark trajectories were designed:

• Pseudo-Static Benchmark: The torque reference ${\mathcal{T}}_d$ ramps from 0 to 8 Nm in the first second and remains constant for the rest of the experiment. During the period from 1 s to 3 s and 5 s to 8.5 s, a flux reference $\beta$ step of 8 Wb is applied with superimposed ripples to emulate disturbances. This test evaluates steady-state tracking and disturbance rejection.
• Dynamic Benchmark: Here, the torque ${\mathcal{T}}_d$​ follows a complex trajectory: it ramps from 0 to 8 Nm in the first second, holds steady until 5 s, then reverses sharply to approximately –8 Nm, followed by a final ramp back to 6 Nm by 10 s. Simultaneously, the flux reference $\beta$ features a sequence of four step changes (switch on to 8 Wb in some time intervals), testing the controller’s performance under rapid transients and flux variations.

Throughout the entire profile, small random perturbations are superimposed to simulate external disturbances, ensuring a challenging scenario for evaluating control robustness, responsiveness, and adaptation under dynamic and load-varying conditions. Figure 18 shows the profiles used during the experimental tests.

Now we present a comparative evaluation of the proposed neural adaptive PBC controller against four advanced control strategies. These include:

• CNO-PBC: classical PBC controll with a nonlinear observer,

• CFE-PBC: classical PBC controll with a fuzzy enhanced MRAS estimator

• CRBF-PBC: classical PBC controll with a RBF network enhanced observer, and a

• DRNN-PBC: deep learning adaptive PBC strategy with RNN network observer.

To evaluate steady-state accuracy and low-dynamic performance, all four control strategies were tested under the pseudo-static benchmark, where the torque reference ${\mathcal{T}}_d$ ramps from 0 to 8 Nm within the first second and is then held constant. During the intervals $\left(1\mathrm{s}\to 3\mathrm{s}\right)$ and $\left(5\mathrm{s}\to 8.5\mathrm{s}\right)$, a flux reference step of 8 Wb is applied, with superimposed ripple signals to emulate realistic disturbance conditions. This scenario is designed to assess tracking precision, overshoot suppression, robustness to steady-state perturbations, and the quality of torque/flux estimation in a low-dynamic environment. A summary of comparative performance under this benchmark is presented in Table 3.

The results in Table 3 confirm a progressive improvement in performance from classical observer-based control to advanced neural network-enhanced strategies. The baseline CNO-PBC controller shows the highest errors in both torque and flux estimation, with a torque MAE of 0.68 Nm and SMAPE of 21.40%, along with the highest overshoot (1.92%) and torque deviation (36.83%). Introducing fuzzy enhancement in CFE-PBC improves tracking accuracy significantly, reducing torque MAE to 0.15 Nm and improving flux SMAPE by nearly 8.4 percentage points. Further enhancement with radial basis function networks (CRBF-PBC) yields even lower errors (torque MAE = 0.09 Nm, SMAPE = 15.13%), shorter rise and settling times, and better rejection of steady-state disturbances. The best overall performance is achieved by the proposed DRNN-PBC controller, which achieves the lowest torque and flux estimation errors (MAE = 0.06 Nm and 0.10 Wb), fastest response ($t_{2\%}$ = 13 ms), and minimal steady-state error (0.05 Hz). It also maintains the smallest overshoot (0.43%) and torque deviation (24.81% of nominal), confirming its robustness, adaptability, and estimation fidelity under pseudo-static conditions. All controllers maintain high model fit with R²≈0.998.

To evaluate the controllers under dynamic and nonlinear conditions, a challenging dynamic benchmark was executed. The torque reference ${\mathcal{T}}_d$​ ramps from 0 to 8 Nm within the first second, holds steady until 5 s, then reverses sharply to approximately –8 Nm before ramping back to 6 Nm by 10 s. Simultaneously, the flux reference β undergoes multiple step changes up to 8 Wb, simulating rapid flux variations. Throughout the test, small random perturbations emulate external disturbances, creating a demanding scenario to assess the controllers’ robustness, responsiveness, and adaptation under rapid transients and load variations. This benchmark emphasizes each control strategy’s ability to maintain estimation accuracy and transient response quality under simultaneous nonlinear torque and flux variations. The results of this dynamic benchmark, including the speed reversal and load variation, are presented in Table 4.

Under dynamic torque and flux variation conditions involving rapid torque ramps, reversals, and flux step changes, the performance disparities between control strategies become more pronounced. The classical CNO-PBC controller exhibits relatively high tracking errors and slower transient response, with a torque disturbance index ($D_{\%}$) of 4.02% and steady-state frequency error ($E_{\mathrm{s}\mathrm{s} }$) of 1.10 Hz. Conversely, the adaptive neuro-based controllers significantly improve tracking precision and responsiveness. Notably, the DRNN-PBC achieves the best overall performance, attaining the lowest torque MAE (0.18 Nm) and SMAPE (18.5%), alongside minimized maximum transient torque deviation (${∆\mathcal{T}}_{\mathrm{m}\mathrm{a}\mathrm{x}}$ = 30.85% nominal). Intermediate controllers, such as CRBF-PBC and CFE-PBC, also provide substantial enhancements in flux and torque estimation accuracy. These findings confirm the superior adaptability and robustness of proposed neural adaptive PBC schemes in handling rapid and complex torque and flux dynamics, ensuring stable and precise control during challenging real-world transient conditions.

6. Conclusions

This work introduced a learning-driven intelligent passivity control strategy for induction motors, integrating a nonlinear state observer with an adaptive recurrent neural framework. The proposed observer significantly enhances the estimation accuracy of key motor variables such as torque, speed, and rotor flux, while simplifying the tuning process by requiring only a single gain adjustment over a wide speed range. Unlike traditional control schemes, this approach ensures rapid convergence and robust performance even under low-torque (i.e. low-speed) and highly dynamic operating conditions. Extensive simulation studies validated the global stability of the closed-loop system using Lyapunov-based analysis and demonstrated the observer’s robustness against various torque disturbances and flux variations. The integration of neural adaptive mechanisms into the passivity-based control framework yielded superior transient response, disturbance rejection, and steady-state tracking performance. Future work will focus on real-time experimental validation and extending the methodology to other induction motor topologies, including doubly-fed and coaxial configurations, to establish the broad applicability and scalability of this intelligent control approach in advanced industrial drive systems.