Observer-Based Suboptimal Controller Design for Permanent Magnet Synchronous Motors: State-Dependent Riccati Equation Controller and Impulsive Observer Approaches

Abstract

Permanent Magnet Synchronous Motors (PMSMs) with high energy efficiency, reliable performance, and a relatively simple structure are widely utilised in various applications. In this paper, a suboptimal controller is proposed for PMSMs without sensors based on the state-dependent Riccati equation (SDRE) technique combined with customised impulsive observers (IOs). Here, the SDRE technique facilitates a pseudo-linearised display of the motor with state-dependent coefficients (SDCs) while preserving all its nonlinear features. Considering the risk of non-available/non-measurable states in the motor due to sensor and instrumentation costs, the SDRE is combined with IOs to estimate the PMSM speed and position states. Customised IOs are proven to be capable of obtaining quality, continuous estimates of the motor states despite the discrete format of the output signals. The simulation results in this work illustrate an accurate state estimation and control mechanism for the speed of the PMSM in the presence of load torque disturbances and reference speed changes. It is clearly shown that the SDRE-IO design is superior compared to the most popular existing regulators in the literature for sensorless speed control.

1. Introduction

Permanent Magnet Synchronous Motors (PMSMs) have advantages such as constant speed, an adjustable power factor, high efficiency, field winding elimination, DC power supply rejection, Joule loss deletion in excitation circuits, omission in brushes and slip rings, a short radial length, and easy maintenance [1,2,3,4,5]. They have simple configurations and easy-to-understand operation mechanisms which make them ideal to be used as actuators in many applications [6,7].

The control of PMSMs has progressed very well in the last few decades. Examples include works such as [8], where Nonlinear Model Predictive Speed Controllers (NMPSCs) are proposed to predict and control the positions and speed at each sampling period through the variable forecast horizon. In [9], to estimate and control the rotor speed and position, an Extended Kalman Filter (EKF) is designed for a sensor-free vector control drive. In [10,11], linear model prediction controllers (MPCs) with the addition of least mean square (LMS) identification strategies and online parameter detection are presented for permanent magnet synchronous generators (PMSGs). In [12], a Robust Adaptive Backstep** Control (RABC) is suggested for High-Speed Permanent Magnet Synchronous Motors (HSPMSM). In [13], a Sliding Mode Observer (SMO) and Sliding Mode Control (SMC) are, respectively, employed to estimate and control the motor speed, torque, and current with third-degree harmonic current. In [14], a model reference adaptive system (MRAS) for estimating and controlling the velocity in internal PMSMs with bandwidth modulation is offered. Furthermore, in [15], an MRAS controller is utilised for sensorless speed control. The authors used the pseudo-linearisation technique and produced state-dependent matrices. Then, the state-dependent matrices were applied to an MRAS mechanism. It should be noted that only the pseudo-linearisation concept was used, and the Riccati equation technique was not investigated in [15]. Moreover, an adaptive, interconnected observer was suggested in order to estimate the speed. In [16], a robust adaptive observer-based approach for the sensorless speed control of PMSMs is studied to estimate the rotor shaft speed and position. In [17], an Adaptive Nonlinear Continuous Observer (ANCO) is used to estimate the self-induction coefficient, stator resistance, rotor velocity/position, and load torque by using the current measurements for the PMSM. In [18], adaptive fuzzy RST digital control (AFRSTD), digital control, and fuzzy digital RST control are generated for a PMSM system. In [19], a Sliding Mode focused on an Extended State Observer (SMESO) and Fast Terminal Sliding Mode Control (FTSMC) is recommended to develop disturbance elimination capacity and increase the dynamic performance of PMSM drive systems.

Reviewing the literature as summarised above reveals that most of the controllers and observers designed for PMSMs have been based on approaches that depend on linear or linearised models. These approaches, although simple, are not ideal for dealing with the inherent nonlinearities of a PMSM. Conventional linearisation methods such as Jacobian approximations do not preserve critical information about the nonlinear mechanisms and are usually time-invariant to reduce the controller design efforts. However, these aspects all reduce the control performance. In the same context, Linear Quadratic Regulators (LQRs) have been shown to perform very well for many linear systems. Therefore, numerous efforts have been made to expand them to nonlinear dynamical systems. The introduction of the state-dependent Riccati equation (SDRE) [20,21,22,23] technique was an effort to fulfil this desire. The SDRE technique is a state-dependent control method based on a pseudo-linearisation display which produces a optimal plan for the cost function of a suboptimal solution [20,21,22,23]. In the SDRE technique, the nonlinear system is transformed into a pseudo-linearised matrix depending on states without any approximation. This pseudo-linearisation representation might not be unique, and there might be countless versions that would make the degrees of freedom of the controller infinite. This is one of the aspects that requires research for its application in different systems.

The SDRE technique has been applied in many applications, such as robotics, navigation systems, Unmanned Aerial Vehicles (UAVs), submarines, Heating Ventilation and Air Conditioning Systems (HVACs), healthcare, and power systems [21,22,23]. In [24], a suboptimal SMC was obtained from a mix of the SDRE and SMC for a class of nonlinear closed-loop systems for the Scout robot. In [25], an SDRE controller was designed for UAVs using the Artificial Potential Field (APF) approach. In [26], the nonlinear SDRE was applied for the fault-tolerant control of a DC micro-grid (MG) in the power system.

On the other hand, the SDREs face challenges regarding access to the system’s states for feedback design. Here, state observers can provide a solution. In [27], the SDRE method is utilised for both control and observation purposes for an Air Handling Unit (AHU) system. The Riccati equation is solved at every sample time, and a suboptimal control law and an estimated state vector are produced based on its solution. It should be noted that the SDRE observer needs continuous data from output sensors at each sample time. In this regard, impulsive systems and observers are well suited for simulating real-world systems in which internal states change instantaneously. Compared to classical continuous and discrete systems, impulsive systems can contain sudden, abrupt, and intermittent changes at variable time intervals [28,29]. In such a mechanism, an impulse could model the sudden change that happens in the system state. The time intervals between these impulses can be fixed or variable with time, known or unknown, and deterministic or stochastic. In [30], an impulsive observer (IO) is introduced based on state-dependent matrices for nonlinear time-delay systems. Furthermore, an adaptive IO (AIO) is suggested for a class of nonlinear time-delay systems with unknown parameters in [31]. Moreover, in [32], an AIO is investigated based on the feasibility of a centralisation perspective for a class of uncertain nonlinear systems.

Based on the literature review summarised and the challenges listed above, in this paper, a state-dependent, centred pseudo-linearisation representation is combined with an impulsive observer to not only estimate motor speed/position in a so-called “sensorless” scheme but also to control these parameters. The proposed observer can provide continuous estimates of the output states of the system, even in the presence of discrete output measurements. Through this approach, by utilising a pseudo-linearisation display instead of linearisation, not only is the design of a linear control scheme facilitated, but also all the nonlinear specifications of the system are preserved along with an accurate estimate of the states.

It should be mentioned that in [27], the authors used an SDRE controller and SDRE observer for an AHU system. The SDRE observer used continuous output signals to estimate the states, and its structure is of the Luenberger type but with state-dependent matrices. Both controller and observer gains are obtained by solving the SDRE at every sample time. In the present paper, the SDRE controller is used with an observer state-dependent impulsive observer for the PMSM when no speed sensors are available. This observer only utilises the output signal at impulse times and does not require continuously estimated states (between the impulses). In [15], the authors used an MRAS to control an IPMSM. The purpose of the mentioned paper was to use the pseudo-linearisation presentation in the structure of the MRAS controller. To estimate the speed, stator resistance, winding inductance, and torque load, an augmented adaptive interconnected observer (AAIO) was used with a similar structure in [33].

In [30], the basic structure of the state-dependent impulsive observer is introduced for nonlinear time-delay systems. The Congo Ebola disease is addressed and presented as a case study without any input control signal criteria. The combination of the SDRE and an impulsive observer was introduced for the first time by the authors of [34] for a manipulator robot case study. The major difference between the present paper and this previous one is that each joint of the manipulator arm was assumed to have a control actuator, and there were the same number of output and control input signals in this system. Thus, the input matrix was an orthogonal one. Unlike the previous case study, the PMSM system in the present paper has only states that are measurable. Therefore, for the speed and position signals, no direct control is available. Also, here, it is assumed that the angles are measured and available only in impulse time stamps, and only angular velocities are estimated by the observer, i.e., there is no access to physical measurements. Addressing this state access problem while there is a nonlinear dependency between the outputs and speed signals is the key challenge addressed in this work. This challenge is addressed through the following:

• Develo** a pseudo-linearised representation of the PMSM system.
• Designing a controller for optimal tracking of the PMSM’s reference speed with high accuracy and a quick speed without a speed sensor.
• Estimating motor speed in a sensorless framework.
• Addressing the challenge of disturbances during the course of the control.
• Maintaining the function of estimation and speed control during all times of sampling and not just at the impulse sample times.
• Quantifying the effects of impulse intervals (the sample rate) and load torque.

The simulated examples in this work demonstrate the performance of the controlled system despite the changes in the reference speed and load torque and the variable sample rates. To show the superiority of the proposed estimation and control approach, the results are compared with the linear control approach of LQRs. It is worth noting that the PMSM here is only considered as a case study, and the proposed work has the potential to be applied to systems with similar challenges and requirements without a loss in generalisability.

Given the introduction and motivation for this work given above, the rest of the paper is organised as follows: Section 2 is dedicated to the state-dependent Riccati equation principles, and Section 3 summarises the fundamentals of the impulsive state estimator. Section 4 includes the main results and the proposed work; the stability theorems and the proofs are given in this section. Section 5 is dedicated to the PMSM control case study, with simulations and visualisations given in Section 6; this section includes a comprehensive discussion of the results, the sensitivity analysis, and practical considerations as well. Comparisons are also included in this section. The last section includes the conclusion and potential future work.

2. State-Dependent Riccati Equation

In the SDRE method, the nonlinear system is first transformed into a pseudo-linearisation form without any approximation, and then the benefits of the linear representation are used to design the controller, observer, or any nonlinear filters. Its ability to be implemented in a wide range of nonlinear systems, the existence of a methodical process for the transformation, and the possibility of obtaining optimal solutions for control purposes are among the advantages of the SDRE technique [21,22,23].

The dynamics of a nonlinear system are considered as follows [21]:

$$\begin{array}{l}\dot{x}\left(t\right)=f\left(x\left(t\right),u\left(t\right)\right)\\ y\left(t\right)=h\left(x\left(t\right)\right),\end{array}$$

(1)

where $x\in R^n$, $u\in R^p$, and $y\in R^q$ are, respectively, the states, control, and output vectors of the system. $f:R^n\times R^p\to R^n$ and $h:R^n\to R^q$ are the nonlinear functions here. The pseudo-linearised display of the system (1) is as follows:

$$\begin{array}{l}\dot{x}\left(t\right)=A\left(x\left(t\right)\right)x\left(t\right)+B\left(x\left(t\right)\right)u\left(t\right)\\ y\left(t\right)=C\left(x\left(t\right)\right)x\left(t\right),\end{array}$$

(2)

$A\left(x\left(t\right)\right)\in R^{n\times n}$, $B\in \left(x\left(t\right)\right)\in R^{n\times p}$, and $C\left(x\left(t\right)\right)\in R^{q\times n}$ are the system’s state-dependent coefficient matrices, inputs, and outputs, respectively. To track the reference path of the system, the cost function with the SDRE’s infinite time horizon is defined as follows:

$$J={\displaystyle {\int }_0^{\infty }\left\{{\left(Cx\left(t\right)-r\left(t\right)\right)}^{T}Q\left(x\right)\left(Cx\left(t\right)-r\left(t\right)\right)+u^T\left(t\right)R\left(x\right)u\left(t\right)\right\}}dt,$$

(3)

where the weight coefficients $Q\left(x\right)$ and $R\left(x\right)$ are state-dependent [24] and can determine state or input preferences. The control input to minimise the cost function (4) is calculated as follows [20]:

$$u\left(t\right)=-K\left(x\right)x\left(t\right)+R^{-1}\left(x\right)B^T\left(x\right)\nu \left(x\right),$$

(4)

such that $K\left(x\right)=R^{-1}\left(x\right)B^T\left(x\right)P\left(x\right)$ and $P\left(x\right)$ are the solutions to the Riccati equation in the form of (5):

$$A^T\left(x\right)P\left(x\right)+P\left(x\right)A\left(x\right)-P\left(x\right)B\left(x\right)R^{-1}\left(x\right)B^T\left(x\right)P\left(x\right)+C^{T}Q\left(x\right)C=0,$$

(5)

where $\nu \left(x\right)$ is the solution to the following algebraic equation:

$$\begin{array}{l}{\left(A\left(x\right)-B\left(x\right)K\left(x\right)\right)}^T\nu \left(x\right)+C^{T}Q\left(x\right)r\left(t\right)=0\\ \to \nu \left(x\right)=-{\left({\left(A\left(x\right)-B\left(x\right)K\left(x\right)\right)}^T\right)}^{-1}\left(C^{T}Q\left(x\right)r\left(t\right)\right).\end{array}$$

(6)

In the process presented above, optimal tracking can be achieved without worrying about the rank reduction in the controllability and observability matrices.

3. Impulsive Observer

In this section, it is described how the state observer is designed to estimate the states of the system. Through this framework, the controller can perform based on the estimated states. Impulsive systems have two continuous and discrete dynamic behaviours. The continuous dynamic behaviours are described by differential equations that are continuous in time, describing the intervals between impulses. However, the discontinuous dynamic behaviour is described by differential equations that denote the impulse moments when system states abruptly change. According to these features, impulsive systems are very suitable for describing real-world variable processes in which system states change instantaneously at certain points in time. Impulsive dynamic systems can be seen as a subset or class of hybrid systems with precipitate jumps. There is a main continuous or discrete dynamic in impulsive systems, in which moments also occur between impulses and change the states of the system. The time interval between these impulses can be fixed or variable with time, known or unknown, and deterministic or stochastic.

In general, for a nonlinear system, the system is denoted as in Equation (7) [28]:

$$\left\{\begin{array}{ll}\dot{x}\left(t\right)=f\left(t,x\left(t\right)\right);& t\ne t_k\\ \Delta x\left(t\right)=f_I\left(x\left(t\right)\right);& t=t_k\end{array}\right.$$

(7)

with $x\in R^n$ as the continuous state vector and $t$ as the time variable. With, $f:R^{+}\times R^n\to R^n$ and $f_I:R^n\to R^n$ as nonlinear functions with appropriate dimensions, the state variable jumps can be expressed as $\Delta x\left(t_k\right)=x\left(t_k^{+}\right)-x\left(t_k\right),$ in discrete intervals of $k=1,2,\dots$$t_k$, where $t_k>t_{k-1}>0$. $\Delta x\left(t_k\right)$ is also known as the Impulsive Differential Equation (IDE). There are two main approaches to studying stability issues in impulsive systems. The classical approach is based on the positive definite Lyapunov function, with a negative derivative for non-impulsive systems. It should be noted that most of the existing articles have used the classical perspective to examine the stability of a system or its observer. Unlike non-impulsive systems, impulsive systems can have a positive derivative at some point in time but have a descending trend overall [28].

A system with impulsive equations is defined as follows:

$$\left\{\begin{array}{l}\dot{w}\left(t\right)=g\left(t,w\left(t\right)\right);\begin{array}{c}\end{array}\begin{array}{c}\end{array}t\ne t_k\\ w\left(t_k^{+}\right)={\psi }_k\left(w\left(t_k\right)\right),\end{array}\right.$$

(8)

where $V\in V_0$ and $g:R^{+}\times R^{+}\to R$ are continuous, and ${\psi }_k:R^{+}\to R^{+}$ is non-descending. Also, the following equations are established:

$$\left\{\begin{array}{l}{D}^{+}V\left(t,x\right)\le g\left(t,V\left(t,x\right)\right);\begin{array}{cc}& \end{array}t\ne t_k\\ V\left(t,x+\Delta x\right)\le {\psi }_k\left(V\left(t,x\right)\right);\begin{array}{c}\end{array}\begin{array}{c}\end{array}t=t_k\end{array}\right.$$

(9)

Hypothesis 1.

It is assumed that $f\left(t,0\right)=0$, $f_I\left(0\right)=0$, and $g\left(t,0\right)=0$ for all $k$ with $t>0$.

Therefore, the obvious response of the main system and the comparison system is in $\left(\left.t_{k-1},t_k\right]\right.$. The theorem of the comparison system is presented as Theorem 1.

Theorem 1.

The conditions are assumed to be as follows:

-

$V:R^{+}\times S_{\rho }\to R^{+}$, $V\in V_0$, $\rho >0$ in $\left(\left.t_{k-1},t_k\right]\right.$, and $D^{+}V\left(t,x\right)\le g\left(t,V\left(t,x\right)\right)$.

-

There is a ${\rho }_0>0$, so $x\in S_{{\rho }_0}$ and $x+\Delta x\in S_{{\rho }_0}$ for all $k$, and in $t=t_k$, $V\left(t,x+\Delta x\right)\le {\psi }_k\left(V\left(t,x\right)\right)$.

-

$b\left(‖x‖\right)\le V\left(t,x\right)\le a\left(‖x‖\right)$ on $R^{+}\times S_{\rho }$ when $a,b\in \kappa$.

Moreover, the stability characteristics of the obvious solution to the comparison system will result in the stability of the main system [28].

Outcome 1.

Suppose $g\left(t,V\left(t,x\right)\right)=\dot{\xi }\left(t\right)V\left(t,x\right)$ when $\xi \in C^1\left[R^{+},R^{+}\right]$ and ${\psi }_k\left(V\left(t,x\right)\right)=d_{k}V\left(t,x\right)$. The origin of the main system will be asymptotically stable if the following conditions are met:

$$\left\{\begin{array}{l}\dot{\xi }\left(t\right)\ge 0\\ d_k\ge 0;\begin{array}{cc}& \end{array}k=1,2,\dots \\ \xi \left(t_k\right)-\xi \left(t_{k-1}\right)+\ln \left(\gamma d_k\right)\le 0;\begin{array}{c}\end{array}\begin{array}{c}\end{array}\gamma >1\end{array}\right.$$

(10)

State-Dependent Impulsive Observer

A state-dependent impulsive observer (SDIO) is suggested as follows:

$$\left\{\begin{array}{ll}\dot{\widehat{x}}=A\left(\widehat{x}\right)\widehat{x}+B\left(\widehat{x}\right)u;& t\ne t_k\\ \widehat{y}=C\widehat{x}\\ \Delta \widehat{x}=F\left(\widehat{x}\right)\left(y-\widehat{y}\right);& t=t_k\end{array}\right.$$

(11)

where $\widehat{x}$ and $\widehat{y}$ are the estimated state and output vectors, respectively. $\Delta \widehat{x}$ is the jump vector of the estimated states. The observer gain is $F$, which is also determined according to the closed-loop stability theorem [29].

4. Main Design: SDRE Controller Based on State-Dependent Impulsive Observer

The impulsive observer is assumed to be state-dependent (11). The control input when using the state-dependent Riccati controller is as follows:

$$\begin{array}{l}u\left(t\right)=-R^{-1}\left(x\right)B^T\left(x\right)P\left(x\right)x\left(t\right)+R^{-1}\left(x\right)B^T\left(x\right)\nu \left(x\right)\\ \nu \left(x\right)=-{\left({\left(A\left(x\right)-B\left(x\right)K\left(x\right)\right)}^T\right)}^{-1}\left(C^{T}Q\left(x\right)r\left(t\right)\right).\end{array}$$

(12)

Hypothesis 2.

The following Lipshitz conditions are assumed to be valid for the system:

$$\begin{array}{l}‖A_c\left(x\right)x-A_c\left(\widehat{x}\right)\widehat{x}‖\le K_A‖e‖\\ ‖B_c\left(x\right)r-B_c\left(\widehat{x}\right)r‖\le K_B‖e‖\end{array}$$

(13)

where $K_A$ and $K_B$ are constant and positive.

Theorem 2.

The proposed state-dependent impulsive observer’s estimation error for the SDRE controller is $e=x-\widehat{x}$, which will converge to zero asymptotically when the following assumptions and conditions are met:

$$\left[\begin{array}{cc}{\sum }_i& \left({\epsilon }_1^{-1}+{\epsilon }_2^{-1}\right)M_i\\ \left({\epsilon }_1^{-1}+{\epsilon }_2^{-1}\right)M_i& -\left({\epsilon }_1^{-1}+{\epsilon }_2^{-1}\right)I\end{array}\right]\le 0;i=1,2$$

(14)

$${\sum }_i=A_c^T\left(\widehat{x}\right)M_i+M_i{A}_c\left(\widehat{x}\right)+\left(M_1-M_2\right)/{\Delta }_k+\left(2{\epsilon }_1{K}_A^2+{\epsilon }_2{K}_B^2\right)I+2{\epsilon }_1{A}_c^T\left(\widehat{x}\right)A_c\left(\widehat{x}\right)-\alpha M_i$$

(15)

$$\left[\begin{array}{cc}-\sigma M_1& M_2-\left(\overline{F}\left(\widehat{x}\right)C^T\right)\\ M_2-\overline{F}\left(\widehat{x}\right)C& -M_2\end{array}\right]\le 0$$

(16)

$$\alpha {\Delta }_k+\ln \left(\gamma \sigma \right)\le 0$$

(17)

where $\alpha \ge 0$, $\sigma \ge 0$, $\gamma \ge 1$, and $\gamma \sigma \le 1$ are the design parameters, and $F\left(\widehat{x}\right)=\overline{F}\left(\widehat{x}\right)M_2^{-1}$.

Proof of Theorem 2.

The SDRE closed-loop system is written as follows:

$$\begin{array}{l}\dot{x}=A\left(x\right)x+B\left(x\right)\left(-k\left(x\right)x+\nu \left(x\right)r\right)\\ =\left(A\left(x\right)-B\left(x\right)k\left(x\right)\right)x+B\left(x\right)\nu \left(x\right)r=A_c\left(x\right)x+B_c\left(x\right)r,\end{array}$$

(18)

where $A_c\left(x\right)=A\left(x\right)+B\left(x\right)k\left(x\right)$ and $B_c\left(x\right)=B\left(x\right)\nu \left(x\right)$. It can also be written for the observer in $t\ne t_k$:

$$\dot{\widehat{x}}=\left(A\left(\widehat{x}\right)-B\left(\widehat{x}\right)k\left(\widehat{x}\right)\right)\widehat{x}+B\left(\widehat{x}\right)\nu \left(\widehat{x}\right)r=A_c\left(\widehat{x}\right)\widehat{x}+B_c\left(\widehat{x}\right)r$$

(19)

Therefore, the dynamics of the estimation error in $t\ne t_k$ are rewritten as follows:

$$\begin{array}{cc}\hfill \dot{e}& =A_c\left(x\right)x+B_c\left(x\right)r-A_c\left(\widehat{x}\right)\widehat{x}+B_c\left(\widehat{x}\right)r\hfill \\ & =A_c\left(\widehat{x}\right)\left(x-\widehat{x}\right)+\left(A_c\left(x\right)-A_c\left(\widehat{x}\right)\right)x+\left(B_c\left(x\right)-B_c\left(\widehat{x}\right)\right)r\hfill \\ & =A_c\left(\widehat{x}\right)e+\tilde{A}x+\tilde{B}r\hfill \end{array}$$

(20)

where $\tilde{A}=A_c\left(x\right)-A_c\left(\widehat{x}\right)$ and $\tilde{B}=B_c\left(x\right)-B_c\left(\widehat{x}\right)$. In $t=t_k$,

$$\Delta e=F\left(\widehat{x}\right)Ce\left(t_k\right)\to e\left(t_k^{+}\right)=\left(I-F\left(\widehat{x}\right)C\right)e\left(t_k\right).$$

(21)

The positive definite Lyapunov function is introduced as follows:

$$\begin{array}{l}V\left(x,t\right)=e^{T}Me;M=\left(1-\rho \right)M_1+\rho M_2,\\ \rho =\frac{t_k-t}{{\Delta }_k};{\Delta }_k=t_k-t_{k-1};t\in \left(t_{k-1},t_k\right],\end{array}$$

(22)

where $M$ is an alternating matrix in $t\in \left(t_{k-1},t_k\right]$, and $M_1$ and $M_2$ are defined as symmetric, positive definite matrices.

$$M\left(t_k\right)=M_1;M\left(t_k^{+}\right)=M_2,$$

(23)

So, the derivative of the Lyapunov function in the time interval $t\in \left(t_{k-1},t_k\right]$ is equal to

$$\dot{V}\left(x,t\right)=e^T\left(A_c^T\left(\widehat{x}\right)M+M{A}_c\left(\widehat{x}\right)+\frac{M_1+M_2}{{\Delta }_k}\right)e+x^T{\tilde{A}}^{T}Me+e^{T}M\tilde{A}x+r^T{\tilde{B}}^{T}Me+e^{T}M\tilde{B}r.$$

(24)

Using the Yang inequality [30], the result is as follows:

$$\begin{array}{l}{x}^T{\tilde{A}}^{T}Me+e^{T}M\tilde{A}x\le {\epsilon }_1{x}^T{\tilde{A}}^T\tilde{A}x+{\epsilon }_1^{-1}{e}^T{M}^{2}e\\ r^T{\tilde{B}}^{T}Me+e^{T}M\tilde{B}r\le {\epsilon }_2{r}^T{\tilde{B}}^T\tilde{B}r+{\epsilon }_2^{-1}{e}^T{M}^{2}e,\end{array}$$

(25)

And, according to Hypothesis 2, the following result is obtained:

$$\begin{array}{l}{x}^T{\tilde{A}}^T\tilde{A}x={‖\tilde{A}x‖}^2\le 2{‖A_c\left(x\right)x-A_c\left(\widehat{x}\right)\widehat{x}‖}^2+2{‖A_c\left(\widehat{x}\right)\left(x-\widehat{x}\right)‖}^2\\ \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\le 2K_A^2{e}^{T}e+2e^T{A}_c^T\left(\widehat{x}\right)A_c\left(\widehat{x}\right)e\end{array}$$

(26)

$$r^T{\tilde{B}}^T\tilde{B}r={‖\tilde{B}r‖}^2\le K_B^2{e}^{T}e.$$

(27)

Hence, the derivative of Lyapunov’s function will be taken as follows:

$$\begin{array}{l}\dot{V}\le e^T(A_c^T\left(\widehat{x}\right)M+M{A}_c\left(\widehat{x}\right)+\left(M_1+M_2\right)/{\Delta }_k+\left(2{\epsilon }_1{K}_A^2+{\epsilon }_2{K}_B^2\right)I\\ +\left({\epsilon }_1^{-1}+{\epsilon }_2^{-1}\right)M^2+2{\epsilon }_1{A}_c^T\left(\widehat{x}\right)A_c\left(\widehat{x}\right))e.\end{array}$$

(28)

By applying $\pm \alpha V\left(x,t\right)$ to (28), the following result is obtained:

$$\dot{V}\left(x,t\right)\le e^T\overline{\sum }e+\alpha V\left(x,t\right),$$

(29)

which can be written as follows:

$$\begin{array}{c}\hfill {\displaystyle \overline{\sum }}={{\displaystyle \overline{\sum }}}_1+\rho {{\displaystyle \overline{\sum }}}_2+{\rho }^2{{\displaystyle \overline{\sum }}}_3,\hfill \\ {{\displaystyle \overline{\sum }}}_1=A_c^T\left(\widehat{x}\right)M_1+M_1{A}_c\left(\widehat{x}\right)+\left(M_1-M_2\right)/{\Delta }_k+\left(2{\epsilon }_1{K}_A^2+{\epsilon }_2{K}_B^2\right)I+\left({\epsilon }_1^{-1}+{\epsilon }_2^{-1}\right)M_1^2+2{\epsilon }_1{A}_c^T\left(\widehat{x}\right)A_c\left(\widehat{x}\right)-\alpha M_1,\hfill \\ {{\displaystyle \overline{\sum }}}_2=A_c^T\left(\widehat{x}\right)\left(M_2-M_1\right)+\left(M_2-M_1\right)A_c\left(\widehat{x}\right)+2\left({\epsilon }_1^{-1}+{\epsilon }_2^{-1}\right)M_1\left(M_2-M_1\right)-\alpha \left(M_2-M_1\right),\hfill \\ {{\displaystyle \overline{\sum }}}_3=\left({\epsilon }_1^{-1}+{\epsilon }_2^{-1}\right){\left(M_2-M_1\right)}^2.\hfill \end{array}$$

(30)

If $0\le \rho \left(t\right)\le 1$, ${{\displaystyle \overline{\sum }}}_3\ge 0$, and the following conditions are satisfied by the LMIs (14), as well as the definition of $\dot{\lambda }\left(t\right)=\alpha \ge 0$, the first condition of the comparison system theorem is established as follows:

$$\begin{array}{c}{{\displaystyle \overline{\sum }}}_1\le 0\\ {{\displaystyle \overline{\sum }}}_1+{{\displaystyle \overline{\sum }}}_2+{{\displaystyle \overline{\sum }}}_3\le 0.\end{array}$$

(31)

In $t=t_k^{+}$, the Lyapunov function is presented as follows:

$$V\left(x,t_k^{+}\right)=e^T\left(t_k^{+}\right)M\left(t_k^{+}\right)e\left(t_k^{+}\right)=e^T\left(t_k\right){\left(I-F\left(\widehat{x}\right)C\right)}^T{M}_2\left(I-F\left(\widehat{x}\right)C\right)e\left(t_k\right).$$

(32)

This is the case if the following condition is met:

$${\left(I-F\left(\widehat{x}\right)C\right)}^T{M}_2\left(I-F\left(\widehat{x}\right)C\right)\le \sigma M_1$$

(33)

which is the same condition for the LMI (16) as $V\left(x,t_k^{+}\right)\le \sigma V\left(x,t_k\right)$ and $d_k=\sigma \ge 0$.

Therefore, the second condition of the comparison system theorem is established. □

According to the third condition, the upper limit of the impact distance is obtained as follows:

$${\Delta }_{k\max }=-\ln \left(\gamma \sigma \right)/\alpha$$

(34)

5. Case Study: Permanent Magnet Synchronous Motor

The dynamic model of the permanent magnet synchronous motor used in this paper on the $d-q$ axis is as follows [35]:

$$\left\{\begin{array}{l}\frac{d{i}_d}{dt}=-\frac{R_s}{L_d}{i}_d+{\omega }_r{i}_q+\frac{V_d}{L_d}\\ \frac{d{i}_q}{dt}=-\frac{R_s}{L_q}{i}_q-{\omega }_r{i}_d-{\lambda }_f{\omega }_r+\frac{V_q}{L_q}\\ \frac{d{\omega }_r}{dt}=\frac{3}{2}\frac{N_p}{J}{\lambda }_f{i}_q-\frac{B}{J}{\omega }_r-\frac{T_L}{J}\\ \frac{d\theta }{dt}={\omega }_r\end{array}\right.$$

(35)

where $i_d$ and $i_q$ are the currents and $V_d$ and $V_q$ are the voltages of the $d$ and $q$ axes, respectively. $R_s$ is the stator resistance. $L_q$ and $L_d$ are, respectively, the inductances of the winding on the $q$ and $d$ axes. ${\lambda }_f$ and ${\omega }_r$ are the magnetic flux of the permanent magnet and the mechanical speed of the motor, respectively. $\theta$ and $J$ are the angular position of the rotor and the moment of inertia. $B$ is the coefficient of friction. $N_p$ is the number of poles, and $T_L$ is the load torque.

Table 1. The PMSM’s parameters [35].

Parameter	Value
Number of poles $\left(N_p\right)$	$4$
Stator resistance $\left(R_s\right)$	0.0875 Ω
Permanent magnetic flux $\left({\lambda }_f\right)$	1 Wb
Inductance $\left(L_d=L_q\right)$	$0.2155\hspace{0.17em}\mathrm{H}$
Moment of inertia $\left(J\right)$	$0.0151\hspace{0.17em}{\mathrm{kgm}}^2$
Friction coefficient $\left(B\right)$	$0.0378\hspace{0.17em}{\mathrm{kgm}}^2{\mathrm{s}}^{-1}$

presents the values of the engine parameters according to [35]. The pseudo-linearised representation of the PMSM is obtained as follows:

$$\begin{array}{l}A\left(x\right)=\left[\begin{array}{ccc}-\frac{R_s}{L_d}& x_3& 0\\ -x_3& -\frac{R_s}{L_q}& -{\lambda }_f\\ 0& \frac{3}{2}\frac{N_p}{J}{\lambda }_f& -\frac{B}{J}\end{array}\right],\hspace{0.17em}\hspace{0.17em}B\left(x\right)=\left[\begin{array}{cc}\frac{1}{L_d}& 0\\ 0& \frac{1}{L_q}\end{array}\right],\hspace{0.17em}C\left(x\right)=\left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\end{array}\right]\\ \end{array}$$

(36)

where $x={\left[\begin{array}{ccc}{i}_d& i_q& {\omega }_r\end{array}\right]}^T$ is the state vector, $u=\left[\begin{array}{cc}{V}_d& V_q\end{array}\right]$ is the input vector, and $y={\left[\begin{array}{cc}{i}_d& i_q\end{array}\right]}^T$ is the output vector. The advantage of this presentation is that it is both controllable and observable.

An overview of the control, speed estimation, and identification of the synchronous motor’s parameters is shown in Figure 1. In this block diagram, the voltage in the d–q plane is applied to the motor, and the currents are captured as the response. These current values are not necessarily available at each sample point, but to reduce the cost of measurements, they are assumed to be available only for a few time stamps. The impulse observer receives the current and voltage values and estimates the motor speed for all of the time stamps. This speed is then used as an input to the SDRE controller for tracking purposes. Here signals with * reference to the setpoint of the system.

6. Simulation Results and Discussions

6.1. The Main Simulation Results

The design parameters are considered as follows:

$$R\left(x\right)=diag\left(1,1\right),\hspace{0.17em}\hspace{0.17em}Q\left(x\right)=diag\left({10}^3+x_1^2,\hspace{0.17em}{10}^3+10x_2^2,\hspace{0.17em}{10}^2+x_3^2\right),\hspace{0.17em}\hspace{0.17em}\alpha =10,\hspace{0.17em}\hspace{0.17em}\gamma =1.001,\hspace{0.17em}\hspace{0.17em}\sigma =0.95.$$

The initial conditions of the states and their estimations are as follows:

$$x\left(0\right)={\left[0,0,1\right]}^T,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\widehat{x}\left(0\right)={\left[10,10,50\right]}^T.$$

Also, the desired values of the speed and current on the d-axis are assumed to be zero. According to the values of the design parameters, the maximum distance between impulses is calculated as ${\Delta }_{k\max }=-\ln \left(1.001\times 0.95\right)/10=0.005$ s. The sampling time $T_s={10}^{-4}$ s, and the interval between impulses is also assumed to be $\Delta =5\times {10}^{-3}$ s, which means that an output is available every five sampling times. The tracking error is defined as $x^{ref}-\widehat{x}$ and the estimation error is called as $x-\widehat{x}$.

The results were simulated by MATLAB 2023b software with an ASUS laptop (Ausu Laptop Model: K556CPU COREi7/7500u, RAM DDR4, Graphic GFORCE 940MX). These results are displayed in Figure 2 for the speed tracking error, the speed estimation error, and the estimated speed jumps (speed impulses). The figures on the top row are the tracking performance and its zoomed-in version. It is very clear that the reference speed was successfully tracked after only the first 0.02 s. The figures in the second row demonstrate the tracking errors, which are very small after convergence happens. According to the results, it is clear that the proposed observer can estimate the motor speed very well. The proposed controller also tracks reference paths with appropriate accuracy. At first, it is seen that, despite the error between the actual and estimated values, these jumps have a greater amplitude, which gradually decreases.

Further results are plotted in Figure 3 to show the d-current estimation, the d-current tracking error, the d-current estimation error, and the estimated d-current jumps (d-current impulses).

Figure 4 illustrates the q-current estimation, the q-current estimation error, and the estimated q-current jumps (q-current impulses). Eventually, as the signal converges to the true value, the jump amplitude also becomes zero.

The control inputs are shown in Figure 5. As can be seen, all inputs and their changes are smooth and slow. The allowed range for the input control signals was between −250 and +250, so it is clear that the signals are within the prescribed range. Therefore, these inputs can be implemented in practice. Generally, in the proposed method, the control signals can be forced to stay within a range by tuning the controller’s parameters. Also, by making a trade-off between the performance and the control input variations, a design that can be implemented in real applications can be achieved.

The observer gains are also plotted in Figure 6, which are obtained at any time by solving the LMIs of Theorem 2. Obviously, the gains remain unchanged after the conversion.

6.2. The Effect of Impulse Intervals

To show the effects of impulse intervals on the estimation error of the states by the proposed observer, the two criteria of the normalised mean square error (NMSE) and the correlation coefficient (CC) were used. These two criteria can be calculated as follows:

$$NMSE=\frac{1}{N}{{\displaystyle \sum _{j=1}^N{\displaystyle \sum _{i=1}^3\left(\frac{e\left(i,j\right)}{\underset{j}{\max \hspace{0.17em}}\left|e\left(i,j\right)\right|}\right)}}}^2$$

(37)

$$c{c}_i=\frac{{\displaystyle {\sum }_{j=1}^{N}x\left(i,j\right)\widehat{x}\left(i,j\right)}}{\sqrt{{\displaystyle {\sum }_{j=1}^N{x}^2\left(i,j\right)}}\sqrt{{\displaystyle {\sum }_{j=1}^N{\widehat{x}}^2\left(i,j\right)}}}$$

(38)

where $e\left(i,j\right)$ is the error of the $i^{\mathrm{th}}$ state. The NMSE metric shows the strength of the estimation and the correlation coefficient of the similarity of the real and estimated states. If the two signals are completely similar, the correlation coefficient is 1, and if they are completely different, the correlation coefficient is −1.

Table 2. Summary of the impulse interval effect.

Metrices	Impulse Intervals (Δ(s))
	${10}^{-4}$	$5\times {10}^{-4}$	${10}^{-3}$	$5\times {10}^{-3}$	$6\times {10}^{-3}$
$NMSE$	45.6892	128.6873	177.1134	301.8423	2521.4897
$C{C}_1$	0.9992	0.9968	0.9925	0.9725	0.6850
$C{C}_2$	1	0.9999	0.9998	0.9994	0.7982
$C{C}_3$	1	1	0.9999	0.9812	0.8241

shows that the accuracy of the estimation decreases with increasing impulse intervals, although the values are still acceptable up to the maximum impulse intervals determined by the proposed Theorem 2. In the last column, where the impulse interval is selected above the maximum value, a high error is witnessed, and the similarity between the actual and estimated signals is reduced.

6.3. The Effect of Load Torque

In this section, to investigate the effect of the load torque as a disturbance to the performance of the closed-loop system, a load torque of 5 Nm at an equilibrium point from $0.25$ to $0.35$ per second is applied to the motor. The load torque impulsive observer gains, speed tracking and estimation with load torque, speed tracking error with load torque, and speed estimation error with load torque are plotted in Figure 7.

According to the results, it is clear that despite the load torque acting on the system as a disturbance, the performance of the proposed controller and observer in tracking reference paths and estimating the states is satisfactory. The tracking error is directly related to increasing or decreasing load torque amplitude. It is worth noting that the load torque was considered such that a change in the speed or current response is visible. Robust approaches should be used if there is more disturbance or higher accuracy is required.

6.4. Comparisons

In this section, a comparison between the SDRE and LQR optimal controllers is provided. In this comparison, the SDRE controller examines the system nonlinearly, and the LQR, which is a linearised version, examines the system linearly. This comparison shows the better performance of the SDRE method. In the SDRE approach, a pseudo-linearisation matrix and SDC are generated, which will preserve the nonlinear properties of the system, whereas in the LQR version, the system is first linearised around the equilibrium point. So, this linearisation cannot include nonlinear features of the system and eliminates them. Hence, the system does not have a nonlinear nature and will have the same system matrices at all times. Meanwhile, the system matrices in the SDRE controller update their state-dependent system matrices at each sampling time. This solves the state-dependent Riccati equation for each time sample. But in LQRs, the Riccati equation is checked only once. The linearisation of the motor equations around the equilibrium point is performed using the Jacobian method as follows:

$$\begin{array}{c}A\left(x\right)=\left(\begin{array}{ccc}\frac{\partial \left({\dot{x}}_1\right)}{\partial u_1}& \dots & \frac{\partial \left({\dot{x}}_1\right)}{\partial u_3}\\ ⋮& \ddots & ⋮\\ \frac{\partial \left({\dot{x}}_3\right)}{\partial u_1}& \cdots & \frac{\partial \left({\dot{x}}_3\right)}{\partial u_3}\end{array}\right)=\left[\begin{array}{ccc}-\frac{R_s}{L_d}& x_{e3}& x_{e2}\\ -x_{e3}& -\frac{R_s}{L_q}& -{\lambda }_f-x_{e1}\\ 0& \frac{3}{2}\frac{N_p}{J}{\lambda }_f& -\frac{B}{J}\end{array}\right],\\ {B=\left(\begin{array}{ccc}\frac{\partial \left({\dot{x}}_1\right)}{\partial u_1}& \dots & \frac{\partial \left({\dot{x}}_1\right)}{\partial u_3}\\ ⋮& \ddots & ⋮\\ \frac{\partial \left({\dot{x}}_3\right)}{\partial u_1}& \cdots & \frac{\partial \left({\dot{x}}_3\right)}{\partial u_3}\end{array}\right)}_{x=x_e}=\left[\begin{array}{cc}\frac{1}{L_d}& 0\\ 0& \frac{1}{L_q}\\ 0& 0\end{array}\right].\end{array}$$

(39)

The weighting coefficients $R$ and $Q$ are considered to be the same in both the LQR and SDRE methods, and the calculation of these coefficients is arbitrary. The response of the algebraic Riccati equation for the matrices (39) and the LQR state feedback gain is obtained as follows:

$$P=\left[\begin{array}{ccc}2.1259& 0.0012& -0.5347\\ 0.0012& 1.5945& 0.0693\\ -0.5347& 0.0693& 0.5487\end{array}\right],\hspace{0.17em}\hspace{0.17em}K=\left[\begin{array}{ccc}9.8650& 0.0055& -2.4810\\ 0.0055& 7.3993& 0.3217\end{array}\right].$$

The following figures show the outputs and states of the SDRE and LQR controllers. Both controllers assume that all states are available, and they track both paths optimally. Figure 8 displays the speed estimation, currents (d-axis and q-axis), and voltages (d-axis and q-axis). It is clear that the SDRE controller performs much better than the LQR, especially when changing the equilibrium point. Because the LQR controller is linearised around the equilibrium point, by changing the equilibrium point, the linearisation accuracy is significantly reduced. Furthermore,

Table 3. Summary of the comparison between the proposed method and the LQR.

Δ(s) = 5 × 10−3	The Proposed Method	LQR
$NMSE$	301.8423	3128.2159
$C{C}_1$	0.9725	0.6014
$C{C}_2$	0.9994	0.7102
$C{C}_3$	0.9812	0.7371

shows a numerical comparison between the proposed method and the LQR. The numerical results verified the efficiency of the proposed framework compared to the LQR technique in varying the reference speed.

As Table 3 shows, for both cases, the distance between the impulses (where the output signals become available) is 0.005 s. In both cases, a similar impulse observer was designed. However, the LQR performs more weakly compared to the SRDE controller. All the performance measures of the NMSE, c c₁, c c₂, and c c₃, which are the normalised mean squared errors and correlation coefficients across both the d and q axes, are lower for the LQR compared to the SDRE. Based on the performance metrics, this superiority is due to the SDRE being able to take care of the system’s intrinsic nonlinearities, especially when the working point changes during tracking.

7. Conclusions

In this article, an observer-based suboptimal controller for the PMSM according to the SDRE controller and impulsive observer approaches was applied, and a pseudo-linearised representation and SDRE strategy procedure were presented. Moreover, a mathematical formulation of the impulsive state observer was presented. The proposed SDRE technique guarantees the tracking of reference paths with acceptable accuracy and speed. Despite the change in the reference speed, the performance of the closed-loop system was maintained, and the reference paths were well tracked. The results illustrate that despite the load torque, the performance of the system decreases to some extent, but it is still acceptable. A comparison was applied between the SDRE method and an LQR, which was designed based on the Jacobian linearisation approach. According to the results obtained, it is clear that when changing the reference speed, the LQR controller cannot guarantee the performance of the system in tracking the reference path. In the future, suggestions can be made for robust and adaptive SDRE controllers to handle disturbances and uncertainty in system models.