Next Article in Journal
A Model of the Control Problem of the Thermal Effect of a Laser Beam on a Two-Layer Biomaterial
Previous Article in Journal
Multi-Target Feature Selection with Adaptive Graph Learning and Target Correlations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Accurate Discrete-Time Dynamic Models of an Induction Machine

by
Ramón Herrera Hernández
1,
Carlos Reusser
1,*,
María Coronel
2 and
Rodrigo Carvajal
1
1
School of Electrical Engineering, Pontificia Universidad Católica de Valparaíso (PUCV), Av. Brasil 2147, Valparaíso 2362804, Chile
2
Department of Electricity, Universidad Tecnológica Metropolitana (UTEM), Av. Jose Pedro Alessandri 1242, Santiago 7800002, Chile
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(3), 373; https://doi.org/10.3390/math12030373
Submission received: 6 December 2023 / Revised: 13 January 2024 / Accepted: 15 January 2024 / Published: 24 January 2024

Abstract

:
Induction machines have become the standard for highly demanding industrial applications. This has led to the utilization of modern discrete-time control techniques (such as model predictive control) that require the estimation of internal variables that are not subject to measurement (such as the rotational velocity in sensorless applications). From this point of view, it is fundamental to have accurate discrete-time models of induction machines, particularly given their nonlinear nature, so that control techniques perform according to design requirements. In spite of the above, the modeling of induction machines has not received much attention in the literature, even though more powerful machines and faster microcontrollers are currently being used. To better understand induction machine models for control, in this paper, we develop and compare various discrete-time models of induction machines based on Euler, Taylor, and Runge–Kutta methods. In addition, we compare the Extended Kalman Filter and Unscented Kalman Filter for state estimation in terms of accuracy and computational burden. The models are derived and compared through extensive Monte Carlo simulations and the state estimation techniques are compared in terms of root mean squared error, execution time, and maximum absolute error. Our simulations show that, in general, the Taylor method yields more accurate models than both the Runge–Kutta and Euler methods. In particular, the Taylor method results in a root mean square error that is one order of magnitude smaller than the Euler method for stator current and rotor flux linkages. For rotor angular speed, the Runge–Kutta methods are more accurate than both the Taylor and Euler methods, resulting in a root mean square error that is two orders of magnitude smaller than the Euler method. On the other hand, the Extended Kalman Filter results in smaller execution time than the Unscented Kalman Filter, up to two orders of magnitude. In terms of root mean squared error and maximum absolute error, both filtering algorithms perform similarly.
MSC:
37M10; 93C55; 93C57; 60G35

1. Introduction

Induction machines (IM) are basically an alternating current (AC) polyphase machine connected to an AC power source, either in the stator or in the rotor. The winding arrangement on the stator produces a rotating magnetic field in the machine airgap, which induces circulating currents that result in an induced rotating magnetic field in the rotor [1].
Induction machines have become the standard for highly demanding industrial applications. This is due to the simplicity in their construction and working principle, robustness, and wide speed range operation [2]. On the other hand, control strategies such as Field Oriented Control or Direct Torque Control are today’s industrial standard control schemes, since they allow for high energy efficiency, i.e., they develop high torque per Ampere ratio [1,2]. Among modern control methods, Model Predictive Control (MPC) has been widely studied in the literature [3,4]. Moreover, technological advances have allowed for the implementation of MPC strategies, see, e.g., [5,6]. It is important to note that power converters play a key role in driving the machine, delivering the necessary input voltage according to the control requirements. In general, one of the practical limitations in machine drives is the switching frequency of the power converter, which in turn limits the maximum sampling rate. This limitation arises from the fact that the higher the voltage the power converter handles, the lower the switching frequency it can withstand [7]. This sampling rate restriction may limit the performance of the controller since a low sampling rate could result in an inaccurate discrete-time model.
In variable-speed drives, speed is conventionally measured using a sensor, which is susceptible to electromagnetic noise; it requires maintenance and it is unsuitable for abrasive environments [8,9]. For these reasons, speed sensors are not used today. Instead, the rotor speed is estimated using an observer. In addition, the benefits of speed estimation encompass diminished hardware complexity and costs, reduced physical space requirements for the drive, and decreased maintenance demands, to name a few [8,9]. On the other hand, modern control techniques typically require knowledge of the rotor magnetic flux position. Given the cost implications and challenges associated with integrating electromagnetic flux sensors directly into the machine, the flux position is commonly estimated. Therefore, only the stator currents are measured whilst rotor speed and magnetic fluxes are estimated. This kind of drive is known as sensorless [8].
To estimate the rotor speed and rotor flux linkages, model-based methods are typically used, such as model adaptive systems [10,11], full-order observers [12], Luenberger observers [13], and sliding mode observers [14,15], among others. The methods previously mentioned are based on a deterministic formulation of the system. However, deterministic models do not take into account parameter offset from the nominal value, friction, lubricant ageing, or the effects of increased temperature (under normal operation) on the windings, all of which can be understood as model uncertainties. These uncertainties are typically modeled using random variables [16,17,18,19]. In this framework, the unknown variables are estimated using Bayesian filtering. In linear systems, the Kalman Filter (KF) is used, whereas for addressing nonlinear state estimation, the Extended Kalman Filter (EKF) and the Unscented Kalman Filter (UKF) are typically used [20,21,22,23,24].
The well-known Kalman Filter has been used in state estimation of induction machines in several works under the assumption of constant speed, see, e.g., [25,26,27]. This assumption results in a linear system, and since the Kalman Filter is the optimal linear filter [28,29], it is the best option in this case. However, in our problem of interest, the speed is not constant, particularly in the numerical examples where we have considered direct start and load application. This non-constant speed results in a nonlinear system, which in turn requires the utilization of nonlinear state estimation techniques. To the best of our knowledge, the Extended Kalman Filter is the most widely used technique for state estimation in induction machines, see, e.g., [30,31,32,33,34]. Additionally, the Euler method is typically used for the attainment of discrete-time induction machine models, and whilst at high sampling rates it can yield adequate models, our work suggests that higher-order methods yield more accurate models that result in more accurate state estimation in terms of root mean square error.
On the other hand, the state estimation using UKF for an induction machine is presented in [35,36,37], where it is compared with EKF. In [35,37], it is determined that UKF has a slightly better performance than EKF in state estimation for the induction machine. Conversely, in [36] it is reported that the EKF yields a slightly better performance than UKF. Moreover, UKF has been used in sensorless control [37,38,39,40]. In [38], in addition to estimating the states, the rotor resistance and the load torque are also estimated. In [40], UKF is compared with EKF, and UKF provides slightly better estimation. However, in terms of execution time, the UKF algorithm takes six times longer than the EKF algorithm.
Currently, due to the digital nature of controllers and processors, it has become necessary to discretize continuous-time models for their implementation. Historically, in machine drives, the Euler method was widely adopted in the literature to obtain a discrete-time model of the machine due to its simplicity, see, e.g., [41,42,43]. In fact, in Bayesian filtering applications on IM drives using EKF [27,32,33,34] and UKF [35,37,38,39,40], Euler discretization methods are typically considered. Higher-order methods like the Taylor or Runge–Kutta, which involve greater computational costs and usually yield more complex systems, are not commonly utilized. The implementation of Taylor discretization, for instance, requires explicit calculations of higher-order derivatives. In contrast, Runge–Kutta employs estimated intermediate points between each sample to construct the model’s evolution [44]. Notably, Runge–Kutta methods are chosen to solve the nonlinear differential equations that define the machine behavior, resulting in “approximate” continuous-time currents, magnetic fluxes, and rotor speed. Moreover, Runge–Kutta methods are extensively utilized in simulation software such as Simulink and Plecs [45,46].
The contribution of this article is twofold: (i) a formulation and comparison of discrete-time models for induction machines using Euler, Taylor, and Runge–Kutta methods are presented, and (ii) a thorough comparison of EKF and UKF applied to the induction machine using the different discrete-time models is presented. Our analysis shows that the Euler method accurately represents the machine’s dynamics, but at very high sampling rates. On the other hand, at lower sampling rates, various discrete-time models based on the Taylor and Runge–Kutta methods can accurately represent the machine model. Finally, for the different discrete-time models, EKF and UKF algorithms are compared through extensive Monte Carlo simulations, in terms of root mean squared error (RMSE), execution time, and maximum absolute error.
The reminder of this article is organized as follows. In Section 2, the continuous-time model of the induction machine is described. In Section 3, the Euler, Taylor, and Runge–Kutta discrete-time models of the induction machine are formulated and compared with the continuous-time model in terms of RMSE. In Section 4, EKF and UKF algorithms are presented based on a stochastic model of the machine. In Section 5, EKF and UKF algorithms using different discrete-time models are compared in terms of execution time, RMSE, and maximum absolute error. Finally, we present our conclusions in Section 6.

2. Continuous-Time Induction Machine Model

The main parts of the induction machine are the stator and the rotor. The stator is made up of three symmetrically distributed coils. On the other hand, the rotor consists of aluminum bars shorted at the ends forming a cylinder. The three-phase coils are represented on a complex space, whose axes are named α-β. The stator space vectors are
v s α β = v s α + j v s β = 2 3 a 0 v a + a v b + a 2 v c
i s α β = i s α + j i s β = 2 3 a 0 i a + a i b + a 2 i c
ψ s α β = ψ s α + j ψ s β = 2 3 a 0 ψ a + a ψ b + a 2 ψ c
where v s α β is the stator voltage space vector on the stationary reference frame α-β; in the same way, i s α β and ψ s α β are the stator current and flux linkage, respectively, on the α-β reference frame. The unit space vector a = e j 2 π 3 denotes the position of the coil on the space. The subscript a , b , c denotes the three-phase component.
The stator electrical equation on the α-β reference frame is
v s α β = R s i s α β + d d t ψ s α β
where R s is the stator resistance.
Similarly, the rotor electrical equation is
0 = R r i r α β + d d t ψ r α β + j p ω r ψ r α β
where i r α β and ψ r α β are the rotor current and flux linkage, respectively, on the α-β reference frame; R r is the rotor resistance; p is the number of pole pairs; and ω r is the rotor shaft angular speed.
The flux linkages of the stator and rotor are
ψ s α β = L s i s α β + L m i r α β
ψ r α β = L r i r α β + L m i s α β
where L s and L r are the stator and rotor self-inductances, respectively, while L m is the mutual inductance between stator and rotor.
On the other hand, the mechanical equation at the rotor shaft is
J d d t ω r = T d T l
where J is the inertia; T d is the machine developed torque; and T l represents the load torque. The machine developed torque is
T d = 3 2 p Re j ψ r α β i s ( α β )
where the superscript denotes the complex conjugate and Re { · } denotes the real part of a complex number.

State-Space Representation

A state-space representation is used to analyze the IM model given in (2)–(9). For convenience, the selected state variables are the stator current i s α β , the rotor flux linkage ψ r α β , and the rotor shaft speed ω r . By manipulating these variables, torque and speed control is achievable [2]. Furthermore, in sensorless drives only the stator current is measured.
Typically, the state vector is expanded and the load torque dynamics is modeled as well [32,39,40]. Including the load torque as a state variable, the IM dynamic equations are
d d t i s α = a 1 i s α + a 2 ψ r α + a 3 ω r ψ r β + b 1 v s α
d d t i s β = a 1 i s β + a 2 ψ r β a 3 ω r ψ r α + b 1 v s β
d d t ψ r α = a 4 i s α a 5 ψ r α a 6 ω r ψ r β
d d t ψ r β = a 4 i s β a 5 ψ r β + a 6 ω r ψ r α
d d t ω r = a 7 ψ r α i s β ψ r β i s α a 8 T l
d d t T l = g ( T l , ω r )
where a and b coefficients are
a 1 = R σ σ L s a 2 = L m σ L s τ r L r a 3 = L m σ L s L r p a 4 = L m τ r a 5 = 1 τ r a 6 = p a 7 = 3 2 1 J p L m L r a 8 = 1 J b 1 = 1 σ L s
The above coefficients are in terms of already-named machine parameters and additionally by
R σ = R s + R r L m 2 L r 2 σ = 1 L m 2 L s L r τ r = L r R r
where R σ is known as the stator dynamic transient resistance; σ as the leakage coefficient; and τ r as the rotor time constant. The function g ( T l , ω r ) represents any load torque function. The load torques that are generally found in electric drives are given in Table 1, where K is a proportional constant.
The dynamic equations are depicted in the block diagrams shown in Figure 1, where the α-β components of the space vectors are jointly represented. The states are obtained by integrating the right-hand side of Equations (10)–(15). Specifically, Figure 1a illustrates the joint behavior of the α-β stator current, while Figure 1b captures the joint dynamic behavior of the α-β rotor flux linkages. Figure 1c represents the dynamic behavior of the rotor shaft angular speed, with ⊗ representing the cross product of the rotor flux linkages and stator currents ( ψ r α β × i s α β ); finally, Figure 1d illustrates the dynamic behavior of the load torque.
In compact form, the representation of the IM model in state-space is
d d t x = f x , u
y = H x
Equation (16) is given by (10)–(15); x is the state vector and u is the input vector, in this case v s α β . The output (17) is related to the state vector by the matrix H . Since it is a sensorless drive, only the stator current i s α β is measured, hence
H = 1 0 0 0 0 0 0 1 0 0 0 0
Note that the state variables depend implicitly on time t (i.e., x = x ( t ) ). Also, note that the IM state Equation (16) exhibits a nonlinear behavior.

3. Discrete-Time Induction Machine Model

Due to the discrete nature of controllers, mostly discrete signals are processed today. In order to achieve a precise and accurate discretization, the discrete-time model must be as close as possible to the continuous-time model; hence, the error between continuous-time and discrete-time models must be close to zero. The discrete-time model is described in the form
x k + 1 = f d x k , u k
y k = H x k
where x k = x ( k T s ) and x k + 1 = x ( k T s + T s ) ; T s is the time interval between samples. The input u k is considered as a zero-order hold (ZOH) sinusoidal wave.
The goal is to approximate numerically the IM dynamic model. There are several methods to numerically approximate differential equations. In this section, three methods are discussed: Euler, Taylor, and Runge–Kutta.

3.1. Euler Method

The Euler method is intuitive and straightforward, and is the basis of the two other methods. It approximates the solution to the differential equation as follows
x k + 1 x k + T s d d t x | x k
The IM model (16) approximated by the Euler Method is
i s , k + 1 α = i s , k α + T s a 1 i s , k α + a 2 ψ r , k α + a 3 ω r , k ψ r , k β + b 1 v s , k α
i s , k + 1 β = i s , k β + T s a 1 i s , k β + a 2 ψ r , k β a 3 ω r , k ψ r , k α + b 1 v s , k β
ψ r , k + 1 α = ψ r , k α + T s a 4 i s , k α a 5 ψ r , k α a 6 ω r , k ψ r , k β
ψ r , k + 1 β = ψ r , k β + T s a 4 i s , k β a 5 ψ r , k β + a 6 ω r , k ψ r , k α
ω r , k + 1 = ω r , k + T s a 7 ψ r , k α i s , k β ψ r , k β i s , k α a 8 T l , k
T l , k + 1 = T l , k + T s g T l , k , ω r , k
For the Euler method, the local truncation error (i.e., the error committed at each time step) is of the order O ( T s 2 ) .

3.2. Taylor Method

The derivative of a function is approximated using Taylor series expansions
x k + 1 x k + T s d d t x | x k + T s 2 2 d 2 d t 2 x | x k + T s 3 3 ! d 3 d t 3 x | x k +
The more terms in the series considered for the approximation, the more accurate the model. According to [48], in order to capture completely the effect of a ZOH-input at the next sampling instant, the input must appear explicitly in the state equations, that is, a discrete-time relative degree 1 . As simply as possible, the discretization of relative degree = 1 is achieved by expanding by second-order Taylor series the rotor flux, speed, and load torque given in (12)–(15). The discrete-time IM model with relative degree = 1 is given by
i s , k + 1 α = i s , k α + T s a 1 i s , k α + a 2 ψ r , k α + a 3 ω r , k ψ r , k β + b 1 v s , k α
i s , k + 1 β = i s , k β + T s a 1 i s , k β + a 2 ψ r , k β a 3 ω r , k ψ r , k α + b 1 v s , k β
ψ r , k + 1 α = ψ r , k α + T s a 4 i s , k α a 5 ψ r , k α a 6 ω r , k ψ r , k β + T s 2 2 a 4 a 1 i s , k α + a 2 ψ r , k α + a 3 ω r , k ψ r , k β + b 1 v s , k α + T s 2 2 a 6 ω r , k a 4 i s , k β a 5 ψ r , k β + a 6 ω r , k ψ r , k α + ψ r , k β a 7 ψ r , k α i s , k β ψ r , k β i s , k α a 8 T l , k T s 2 2 a 5 a 4 i s , k α a 5 ψ r , k α a 6 ω r , k ψ r , k β
ψ r , k + 1 β = ψ r , k β + T s a 4 i s , k β a 5 ψ r , k β + a 6 ω r , k ψ r , k α + T s 2 2 a 4 a 1 i s , k β + a 2 ψ r , k β a 3 ω r , k ψ r , k α + b 1 v s , k β + T s 2 2 a 6 ω r , k a 4 i s , k α a 5 ψ r , k α a 6 ω r , k ψ r , k β + ψ r , k α a 7 ψ r , k α i s , k β ψ r , k β i s , k α a 8 T l , k T s 2 2 a 5 a 4 i s , k β a 5 ψ r , k β + a 6 ω r , k ψ r , k α
ω r , k + 1 = ω r , k + T s a 7 ψ r , k α i , s k β ψ r , k β i s , k α a 8 T l , k + T s 2 2 a 7 i s , k β a 4 i s , k α a 5 ψ r , k α a 6 ω r , k ψ r , k β + T s 2 2 a 7 ψ r , k α a 1 i s , k β + a 2 ψ r , k β a 3 ω r , k ψ r , k α + b 1 v s , k β T s 2 2 a 7 i s , k α a 4 i s , k β a 5 ψ r , k β + a 6 ω r , k ψ r , k α T s 2 2 a 7 ψ r , k β a 1 i s , k α + a 2 ψ r , k α + a 3 ω r , k ψ r , k β + b 1 v s , k α T s 2 2 a 8 g ( T l , k , ω r , k )
T l , k + 1 = T l , k + T s g ( T l , k , ω r , k ) + T s 2 2 g ˙ ( T l , k , ω r , k )
For a second-order Taylor series expansion, the local truncation error is O ( T s 3 ) .

3.3. Runge–Kutta Method

The Runge–Kutta (RK) method approximates the solution of the differential equation based on estimates of the solution at different points. In this subsection, two RK methods are taken into account: RK 2 and RK 4. The RK 2 approximation considers two estimates, whereas RK 4 considers four. The advantage that Runge–Kutta has over Taylor is that it is not necessary to explicitly calculate higher-order derivatives of f ( x , u ) .

3.3.1. Runge–Kutta 2

Runge–Kutta 2 is a second-order approximation. This method requires two evaluations of f ( x , u ) at each time step. The IM model (16) is approximated by
x k + 1 = x k + T s 2 r 1 + r 2
where
r 1 = f ( x k , u k )
r 2 = f x k + T s r 1 , u ( t k + T s )
and t k = k T s . Note that r 2 estimation requires future input values. Since the system is assumed to be causal and future inputs are unknown, u k is used instead. For the IM, r 1 and r 2 are found in Appendix A.1.

3.3.2. Runge–Kutta 4

Runge–Kutta 4 is a fourth-order method that requires four evaluations of f ( x , u ) in one sample period. The approximation of IM model (16) by RK 4 is
x k + 1 = x k + T s 6 r 1 + 2 r 2 + 2 r 3 + r 4
where
r 1 = f ( x k , u k )
r 2 = f x k + T s 2 r 1 , u ( t k + T s 2 )
r 3 = f x k + T s 2 r 2 , u ( t k + T s 2 )
r 4 = f x k + T s r 3 , u ( t k + T s )
Since the input is ZOH and no future input values are taken into account, u k is used in r 2 , r 3 , and r 4 . For the IM, the estimate points r 1 , r 2 , r 3 , and r 4 are found in Appendix A.2.

3.4. Comparison between Models

The previously mentioned discrete-time IM models were contrasted with each other and with a continuous-time model. The simulations were carried out in Matlab. To compare the discrete-time models with the continuous-time one, the latter was obtained using a fifth-order Dormand–Prince (DOPRI 5) method with the same discretization time as the discrete-time models. More information about this method is given in Appendix B. For the simulations, T s = 200 μs and a 4 kW IM was used; its nominal parameters are V = 380 V; I = 8.6 A; f = 50 Hz; two pole pairs; R s = 1.32 Ω ; R r = 2.63 Ω ; L m = 0.1889 H; L s = 0.1972 H; L r = 0.2012 H; and J = 0.528 kg·m2.
To illustrate the differences between discrete-time and continuous-time models, two scenarios are taken into account. First, the IM is directly started from the grid. Second, a 15 Nm step-type load impact occurs suddenly at 4 s. The total simulation time was 6 s.
Table 2 contrasts the discrete-time models with the continuous-time model by presenting the computed root mean squared error (RMSE) between them, calculated as follows
R M S E ( x ) = 1 N i = 1 N x i , c o n t x i , d i s c 2
To begin with, in Table 2, the best approximation for the stator current and rotor flux linkage is obtained by Taylor. The RK 2 and RK 4 models are not as accurate as Taylor, due to the fact that RK 2 and RK 4 do not take into account future input estimation for the approximation. Conversely, for the rotor speed RK 4 and RK 2 provide the best approximation. Finally, Euler results in the worst approximation as expected.
Figure 2 shows the α-β stator currents and rotor flux linkages for the discrete-time and continuous-time models with no load condition at steady state. Figure 2a shows that the Taylor and RK 4 methods yield discrete-time values that are almost identical to the continuous-time values. In the same figure, we can observe that the RK 2 method yields a larger difference since the amplitude of the currents is smaller than the continuous-time currents, whilst the Euler method results in even larger errors since the resulting currents have a much smaller amplitude. In Figure 2b, the rotor flux linkages for each method are shown. The Taylor, RK 2, and RK 4 signals resemble the continuous-time signals. However, the Euler method results in a slightly larger difference from the continuous-time flux than the Taylor and RK methods. On the other hand, Figure 3 illustrates the rotor angular speed and zero load torque for all models described. As shown, both the Taylor and Euler models are less accurate under no-load conditions than the RK methods. In particular, RK 4 provides the most accurate speed approximation, followed by RK 2.
On the other hand, notice that after the sudden step-type load torque impact occurs, stator current increases while rotor speed decreases. Figure 4 shows the stator currents and rotor flux linkages for the discrete-time and continuous-time models after the torque impact occurs. In this case, the Euler method yields a bad approximation, since the resulting currents not only exhibit a smaller amplitude but also a large phase error. In contrast, the Taylor, RK 2, and RK 4 models resemble the continuous-time model. Figure 5 illustrates the rotor angular speed and the load torque impact for all models. We can observe that, once the speed reaches the steady state, RK 4 and Taylor result in a very similar curve to the continuous-time model, followed by RK 2 and finally Euler. In addition, since load torque does not depend on rotor angular speed or any other state variable, the differences between the discrete-time models and the continuous-time model are, in general, small, as shown in Figure 5b.
In general, accurate models are desired since, in modern industrial systems, modern control techniques are applied. If the model is not accurate, model-based observers may fail to represent the actual behavior of the states, providing inaccurate feedback to the controller, hence resulting in poor control performance. Since with the Euler method the model error is reflected not only in the magnitude but also in the phase of the signals of interest, the overall control loop’s performance decreases greatly. This situation is naturally avoided in order to satisfy control system performance and requirements.

Changes in Sample Time

In this section, the behavior of the RMSE between the discrete-time models and the continuous-time model is observed as the discretization time varies. The discretization time is adjusted until the error of the model discretized by Euler is nearly zero. Two scenarios are considered; first, the IM is directly started, and then a sudden load impact occurs at 4 s. This follows the same approach as the previous section.
Figure 6 illustrates how RMSE varies according to T s for α-β stator currents, α-β rotor flux linkages, and rotor shaft speed. Since the error between the α and β components of the stator current and rotor flux linkage components is similar, the average value is used in this figure. As expected, as the sampling time decreases, error decreases as well. In particular, the Taylor discrete-time model is better than RK 2 and RK 4 in terms of stator current and rotor flux linkage RMSE. Conversely, RK 2 and RK 4 are more accurate for speed approximation. Regarding the load torque approximation, RK 2 and RK 4 are better; nevertheless, Euler and Taylor are near zero. On the other hand, Euler RMSE rises dramatically as T s increases. Thus, using Euler for approximating the IM model is not advisable when the T s exceeds 50 μs. To conclude, Taylor is the best approximation, prioritizing stator currents and rotor flux linkage approximations.

4. Filtering

The continuous-time deterministic model of the IM was described in (16) and (17). Even so, those equations do not include model uncertainties or unknown disturbances, such as winding temperature rises, mechanical vibrations, switching noise, etc. All these uncertain effects can be represented by adding Gaussian noise in the state Equation (16) and output Equation (17) [16,17,18,19]. Also, considering the discrete nature of the measurement and control, a discrete-time representation of the model is used. Therefore, the discrete-time stochastic model of the IM is as follows
x k + 1 = f d x k , u k + w k
y k = H x k + v k
where w k and v k are Gaussian noise with zero mean and covariance matrix Q and R , respectively. The term w k is associated with model uncertainties, while v k represents measurement noise.
Under a stochastic point of view, the state estimation is obtained from noisy measurements. In Bayesian inference, the estimation corresponds to calculating the joint probability density function (pdf) of the states given the measurements P ( x 0 : T | y 1 : T ) , which incorporates all the statistical information about the state P ( x 0 : T ) that is contained in the available observation P ( y 1 : T | x 0 : T ) and the initial condition P ( x 0 : T )  [19,21]. Using the Bayes theorem, the joint pdf is
P ( x 0 : T | y 1 : T ) = P ( y 1 : T | x 0 : T ) P ( x 0 : T ) P ( y 1 : T )
where the denominator P ( y 1 : T ) is a normalization constant. It is necessary to mention that x 0 : T and y 1 : T are the time series of the state vector { x 0 , , x T } and of the measurements { y 1 , , y T } vector, respectively. For real-time applications, it is unnecessary and computationally inefficient to compute the full joint distribution. For this reason, the marginal distribution of the current state given all the measurements P ( x k | y 1 : k ) , also known as posterior, is computed instead. Additionally, assuming the IM model to be Markovian, it is possible to compute the posterior recursively [19,21]. The recursion is based on calculating the posterior distribution in two steps, starting from an initial prior distribution P ( x 0 ) .
  • First: predict. The prediction of the distribution of the states at the current time k is computed from the dynamic model by means of the Chapman–Kolmogorov equation
    P ( x k | y 1 : k 1 ) = P ( x k | x k 1 ) P ( x k 1 | y 1 : k 1 ) d x k 1
  • Second: update. Given the measurement y k at time k, the posterior distribution is calculated by solving
    P ( x k | y 1 : k ) = 1 Z k P ( y k | x k ) P ( x k | y 1 : k 1 )
    where Z k is a normalization constant.
Unfortunately, since the IM model is nonlinear, it is not possible to solve the filtering Equations (46) and (47) in closed form. Nevertheless, there exist state estimation methods for nonlinear models; two of the classic methods are based on the Kalman Filter (KF): the Extended Kalman Filter (EKF) and Unscented Kalman Filter (UKF). Both are addressed in this section.

4.1. Extended Kalman Filter

The most well-known nonlinear filtering method is the Extended Kalman Filter, which propagates the mean and covariance through a linearization of the system resulting in a Gaussian posterior pdf. However, EKF does not ensure the convergence of the estimation, since it depends on the linearization of the model around an operating point that is also estimated [20,21,22,23,24]. The implementation steps of EKF for N samples are summarized in Algorithm 1, where x ˆ is the state estimate and P is its covariance matrix. The subscript k denotes the posterior and k | k 1 the prior. K is the Kalman gain; I is the identity matrix; and F k | k 1 is the Jacobian matrix of f d x , u evaluated at k 1 .
Algorithm 1: EKF Algorithm
Mathematics 12 00373 i001

4.2. Unscented Kalman Filter

On the other hand, there is the Unscented Kalman Filter, which by means of the Unscented Transform (UT) propagates representative points of the pdf through the nonlinearity of the system to estimate the state. The advantage that UKF has over EKF is that the representative points are adjustable by the user, attaining in general a better estimation [20,21,22,23,24]. These representative points, named Sigma Points (SPs), ensure the propagation of the mean and the covariance [49]. A total of 2 n + 1 SP are required, with n = dim y . The SPs are computed as follows
X ( 0 ) = x
X ( i ) = x + ( n + λ ) P i i = 1 , . . . , n
X ( i ) = x ( n + λ ) P i i = n + 1 , . . . , n
where λ is a scale factor that modifies the dispersion of the SP around the mean. The term P i is computed using Cholesky factorization. Also, weight factors W i ( m ) and W i ( c ) are associated to the mean and covariance, defined as
W 0 ( m ) = λ n + λ W 0 ( c ) = λ n + λ + ( 1 α 2 + β ) W i 0 ( m ) = λ 2 ( n + λ ) W i 0 ( c ) = λ 2 ( n + λ )
where typically, β = 2 for Gaussian distributions; 0 α 1 ; λ = α 2 ( n κ ) n with κ = 3 n .
The steps to implement the UKF are summarized in Algorithm 2. S k is the covariance matrix of the measurement; C k is the cross-covariance between the state and the measurement; and K k is the filter gain.
Algorithm 2: UKF algorithm
Mathematics 12 00373 i002

5. Numerical Examples

In this section, EKF and UKF were implemented for each IM discrete-time model discussed above. The simulations were carried out in Matlab. For the simulations, all the initial condition are zero; the discretization time is T s = 200 μs; the noise covariance matrices are Q = d i a g 2.12 × 10 2 , 2.12 × 10 2 , 10 6 , 10 6 , 10 3 , 9.64 × 10 4 and R = d i a g 1 / 9 , 1 / 9 . The SP coefficients were selected empirically as in [40], where α = 10 1 , β = 2 , and κ = 3 . The same noise covariance matrices were used for EKF and UKF. To test the IM models in filtering, a thousand Monte Carlo simulations were carried out, which is equivalent to measuring the same event a thousand times, in this case a direct start from a 50 Hz grid. All models in both filtering algorithms were compared in terms of execution time RMSE, where the actual system was simulated using DOPRI 5. Additionally, the maximum absolute error for each model using EKF and UKF was acquired to compare the maximum deviation from the actual value for each state, particularly at startup and after starting.
Table 3 and Table 4 show the average RMSE from the Monte Carlo simulations in filtering using EKF and UKF, respectively. Both tables consider the Euler, Taylor, and Runge–Kutta IM models. From the tables, for the same IM model a considerable similarity between EKF and UKF is seen. In regard to IM models in both algorithms, there was a significant improvement in state estimation with the Taylor and Runge–Kutta models as expected. In currents and flux linkages, Taylor has less error than Runge–Kutta, due to the fact that RK models do not take into account future values of the input since they are unknown. Nevertheless, Runge–Kutta speed estimation is better than Taylor. Finally, load torque estimation is nearly the same in all discrete-time models.
On the other hand, Table 5 shows average execution time per sample of the EKF and UKF algorithms for different IM models. This analysis is of particular interest in practical applications since IMs are usually coupled with modern control techniques that rely on state estimation, among other things. Hence, a low execution time allows for carrying out other computing tasks, such as optimization in a Model Predictive Control setup [6]. From Table 5, we can observe that as the complexity of the model increases, so does the execution time when utilizing both EKF and UKF. In particular, when utilizing UKF, apart from Euler, RK 2 has the smallest execution time followed by Taylor and finally RK 4. As for EKF and UKF, clearly the former requires much less execution time than the latter, no matter which discretization technique is used. In fact, our results show that UKF requires more than half of the sampling interval during the execution of the algorithm. This could leave not enough time for other tasks. However, more research is needed in this area in a practical setup since filtering can be implemented in different computing platforms with different architectures, such as a microcontroller or a field programmable gate array.
Table 6 and Table 7 show the maximum absolute error at the startup of each state for each model using EKF and UKF algorithms, respectively, at startup. The maximum error is not symmetric, especially in the α - β components of stator currents and rotor flux linkages. This asymmetry is primarily due to the high uncertainty of the rotor angular speed. At the start, all initial conditions are zero; therefore, obtaining an accurate approximation is challenging. Notice that the maximum absolute errors are similar for all the discretization methods and both EKF and UKF. By themselves, these values only show what the maximum deviation is but not the overall behavior of each technique. However, they show what to expect in terms of the worst estimation in a single experiment. Moreover, if we also consider the average root mean squared error from Table 3 and Table 4, we can conclude that the Taylor and RK methods exhibit a reduced spread around the mean when compared to the Euler method. On the other hand, Table 8 and Table 9 illustrate the maximum absolute error after the startup of each state for each model using EKF and UKF algorithms, respectively. For the Euler, Taylor, and RK 2 models, EKF and UKF’s maximum absolute error are similar; in particular, UKF RK 4 is slightly less than EKF RK 4. Once the rotor speed reaches the steady state, the estimation shows better performance compared to the startup.

6. Conclusions

In this article, we formulated and compared different discretization techniques, namely, Euler. The results show that the Taylor method exhibits slightly lower RMSE compared to RK methods for stator currents and rotor flux linkages. This is mainly because, in the RK methods, it was assumed that the future values of the input were equal to the input at the current time instant. Since a direct start from the grid was simulated, the retention of the input does not fully reflect the continuous-time sinusoidal input behavior. However, when implementing these methods with an input provided by an inverter, input retention is fulfilled at every time instant. Therefore, a more accurate representation of the system is expected in practical setups.
Overall, Taylor exhibits the smallest α-β stator current and rotor flux linkage RMSE, followed by RK 4 and then RK 2. Conversely, the RK 4 model demonstrates the smallest rotor shaft speed error, followed by the RK 2 and Taylor models. The Taylor method resulted in a similar rotor shaft speed error to the Euler method under no-load conditions; however, when a load torque was applied, the Taylor method exhibited a performance similar to RK models and signals that resembled the continuous-time signals. Regarding the load torque, all models performed similarly. These behaviors persisted as T s varied. Also, it was demonstrated that the Euler method can only accurately represent the machine at very high sampling rates.
We also compared state estimation techniques, namely, EKF and UKF, utilizing the discrete-time models from the Euler, Taylor, and RK methods. EKF and UKF show similar RMSE; however, in terms of execution time, EKF required fewer computations. In most cases, as the complexity of the model increases, the execution time increases as well. Nonetheless, for the UKF algorithm, the RK 2 model is faster than the Taylor model. The execution time required for one iteration of the EKF and UKF utilizing the Taylor and RK methods is less than the sampling period. Therefore, these techniques are suitable for implementation in modern control techniques. In contrast, when utilizing the Euler method at high sampling rates, both filtering techniques result in an execution time per iteration that exceeds the sampling period. Finally, in terms of the maximum error, EKF and UKF are similar; therefore, both algorithms are similar in the maximum deviation from the actual state. Particularly, UKF RK 4 is slightly better than EKF RK 4.

Author Contributions

Each author played an essential role in the development of this research. Conceptualization, R.H.H., C.R. and R.C.; formal analysis, R.H.H., C.R., M.C. and R.C.; funding acquisition, M.C. and R.C.; investigation, R.H.H., C.R., M.C. and R.C.; methodology, R.H.H., C.R., M.C. and R.C.; resources, C.R. and R.C.; software, R.H.H.; validation, R.H.H., C.R., M.C. and R.C.; visualization, R.H.H., C.R., M.C. and R.C.; writing—original draft preparation, R.H.H., C.R., M.C. and R.C.; writing—review and editing, R.H.H., C.R., M.C. and R.C. All authors have read and agreed to the published version of the manuscript.

Funding

DI Regular 039.309/2023, Pontificia Universidad Católica de Valparaíso; Pontificia Universidad Católica de Valparaíso Postgraduate Research Scholarship; ANID-FONDECYT grant number: 3230398, 11201187, and 1211676; Advanced Center for Electrical and Electronic Engineering, AC3E, Basal Project FB0008, ANID, Chile.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
IMInduction Machine
KFKalman Filter
EKFExtended Kalman Filter
UKFUnscented Kalman Filter
UTUnscented Transform
MPCModel Predictive Control
RMSERoot Mean Squared Error
ZOHZero-Order Hold
RKRunge–Kutta
DOPRIDormand–Prince
pdfProbability Density Function
SPSigma Points

Appendix A. Runge–Kutta Estimate Points

In this appendix section, Runge–Kutta estimate points of the evaluation of the IM model are presented. As discussed, future estimation points of the input are not possible; therefore, input values at the actual sample instant k are taken into account instead.

Appendix A.1. RK 2

For Runge–Kutta 2, the two estimate points are
r 1 = a 1 i s , k α + a 2 ψ r , k α + a 3 ω r , k ψ r , k β + b 1 v s , k α a 1 i s , k β + a 2 ψ r , k β a 3 ω r , k ψ r , k α + b 1 v s , k β a 4 i s , k α a 5 ψ r , k α a 6 ω r , k ψ r , k β a 4 i s , k β a 5 ψ r , k β + a 6 ω r , k ψ r , k α a 7 ψ r , k α i s , k β ψ r , k β i s , k α a 8 T l , k g ( T l , k , ω r , k )
r 2 = a 1 i s , k α + T s r 1 ( 1 ) + a 2 ψ r , k α + T s r 1 ( 3 ) + a 3 ω r , k + T s r 1 ( 5 ) ψ r , k β + T s r 1 ( 4 ) + b 1 v s , k α a 1 i s , k β + T s r 1 ( 2 ) + a 2 ψ r , k β + T s r 1 ( 4 ) a 3 ω r , k + T s r 1 ( 5 ) ψ r , k α + T s r 1 ( 3 ) + b 1 v s , k β a 4 i s , k α + T s r 1 ( 1 ) a 5 ψ r , k α + T s r 1 ( 3 ) a 6 ω r , k + T s r 1 ( 5 ) ψ r , k β + T s r 1 ( 4 ) a 4 i s , k β + T s r 1 ( 2 ) a 5 ψ r , k β + T s r 1 ( 4 ) + a 6 ω r , k + T s r 1 ( 5 ) ψ r , k α + T s r 1 ( 3 ) a 7 ψ r , k α + T s r 1 ( 3 ) i s , k β + T s r 1 ( 2 ) ψ r , k β + T s r 1 ( 4 ) i s , k α + T s r 1 ( 1 ) a 8 T l , k g ( T l , k + T s r 1 ( 6 ) , ω r , k + T s r 1 ( 6 ) )
where r ( i ) is the i-th row of the vector r .

Appendix A.2. RK 4

For Runge–Kutta 4, the first estimate point r 1 is equal to (A1); the other three are
r 2 = a 1 i s , k α + T s 2 r 1 ( 1 ) + a 2 ψ r , k α + T s 2 r 1 ( 3 ) + a 3 ω r , k + T s 2 r 1 ( 5 ) ψ r , k β + T s 2 r 1 ( 4 ) + b 1 v s , k α a 1 i s , k β + T s 2 r 1 ( 2 ) + a 2 ψ r , k β + T s 2 r 1 ( 4 ) a 3 ω r , k + T s 2 r 1 ( 5 ) ψ r , k α + T s 2 r 1 ( 3 ) + b 1 v s , k β a 4 i s , k α + T s 2 r 1 ( 1 ) a 5 ψ r , k α + T s 2 r 1 ( 3 ) a 6 ω r , k + T s 2 r 1 ( 5 ) ψ r , k β + T s 2 r 1 ( 4 ) a 4 i s , k β + T s 2 r 1 ( 2 ) a 5 ψ r , k β + T s 2 r 1 ( 4 ) + a 6 ω r , k + T s 2 r 1 ( 5 ) ψ r , k α + T s 2 r 1 ( 3 ) a 7 ψ r , k α + T s 2 r 1 ( 3 ) i s , k β + T s 2 r 1 ( 2 ) ψ r , k β + T s 2 r 1 ( 4 ) i s , k α + T s 2 r 1 ( 1 ) a 8 T l , k g ( T l , k + T s 2 r 1 ( 6 ) , ω r , k + T s 2 r 1 ( 6 ) )
r 3 = a 1 i s , k α + T s 2 r 2 ( 1 ) + a 2 ψ r , k α + T s 2 r 2 ( 3 ) + a 3 ω r , k + T s 2 r 2 ( 5 ) ψ r , k β + T s 2 r 2 ( 4 ) + b 1 v s , k α a 1 i s , k β + T s 2 r 2 ( 2 ) + a 2 ψ r , k β + T s 2 r 2 ( 4 ) a 3 ω r , k + T s 2 r 2 ( 5 ) ψ r , k α + T s 2 r 2 ( 3 ) + b 1 v s , k β a 4 i s , k α + T s 2 r 2 ( 1 ) a 5 ψ r , k α + T s 2 r 2 ( 3 ) a 6 ω r , k + T s 2 r 2 ( 5 ) ψ r , k β + T s 2 r 2 ( 4 ) a 4 i s , k β + T s 2 r 2 ( 2 ) a 5 ψ r , k β + T s 2 r 2 ( 4 ) + a 6 ω r , k + T s 2 r 2 ( 5 ) ψ r , k α + T s 2 r 2 ( 3 ) a 7 ψ r , k α + T s 2 r 2 ( 3 ) i s , k β + T s 2 r 2 ( 2 ) ψ r , k β + T s 2 r 2 ( 4 ) i s , k α + T s 2 r 2 ( 1 ) a 8 T l , k g ( T l , k + T s 2 r 2 ( 6 ) , ω r , k + T s 2 r 2 ( 6 ) )
r 4 = a 1 i s , k α + T s r 3 ( 1 ) + a 2 ψ r , k α + T s r 3 ( 3 ) + a 3 ω r , k + T s r 3 ( 5 ) ψ r , k β + T s r 3 ( 4 ) + b 1 v s , k α a 1 i s , k β + T s r 3 ( 2 ) + a 2 ψ r , k β + T s r 3 ( 4 ) a 3 ω r , k + T s r 3 ( 5 ) ψ r , k α + T s r 3 ( 3 ) + b 1 v s , k β a 4 i s , k α + T s r 3 ( 1 ) a 5 ψ r , k α + T s r 3 ( 3 ) a 6 ω r , k + T s r 3 ( 5 ) ψ r , k β + T s r 3 ( 4 ) a 4 i s , k β + T s r 3 ( 2 ) a 5 ψ r , k β + T s r 3 ( 4 ) + a 6 ω r , k + T s r 3 ( 5 ) ψ r , k α + T s r 3 ( 3 ) a 7 ψ r , k α + T s r 3 ( 3 ) i s , k β + T s r 3 ( 2 ) ψ r , k β + T s r 3 ( 4 ) i s , k α + T s r 3 ( 1 ) a 8 T l , k g ( T l , k + T s r 3 ( 6 ) , ω r , k + T s r 3 ( 6 ) )
where r ( i ) is the i-th row of the vector r .

Appendix B. DOPRI 5

Dormand–Prince (DOPRI) belongs to the family of Runge–Kutta methods; however, it is generally employed as an adaptive method, adjusting the sampling time by comparing the error between a fourth- and fifth-order approximation, with the latter being deemed more accurate. In this particular case, given the fixed sampling time, only the fifth-order approximation is utilized.
DOPRI requires evaluating the function six times, resulting in a longer execution time compared to the methods presented earlier. However, DOPRI aims to emulate the continuous-time model; hence, the execution time is not taken into account.
The DOPRI equations are
x k + 1 = x k + T s 35 384 r 1 + 500 1113 r 3 + 125 192 r 4 2187 6784 r 5 + 11 84 r 6
where
r 1 = f t , x
r 2 = f t + 1 5 T s , x + T s 1 5 r 1
r 3 = f t + 3 10 T s , x + T s ( 3 40 r 1 + 9 40 r 2 )
r 4 = f t + 4 5 T s , x + T s ( 44 45 r 1 56 15 r 2 + 32 9 r 3 )
r 5 = f t + 8 9 T s , x + T s ( 19372 6561 r 1 25360 2187 r 2 + 64448 6561 r 3 212 729 r 4 )
r 6 = f t + T s , x + T s ( 9017 3168 r 1 355 33 r 2 + 46732 5247 r 3 + 49 176 r 4 5103 18656 r 5 )

References

  1. Novotny, D.W.; Lipo, T.A. Vector Control and Dynamics of AC Drives; Monographs in Electrical and Electronic Engineering; Clarendon Press: Oxford, UK, 1996. [Google Scholar]
  2. Sul, S. Control of Electric Machine Drive Systems; Wiley: Hoboken, NJ, USA, 2010. [Google Scholar] [CrossRef]
  3. Rawlings, J.B. Model Predictive Control Theory and Design; Nob Hill Pub.: Madison, WI, USA, 2009; 533p. [Google Scholar]
  4. Wang, L. Model Predictive Control System Design and Implementation Using MATLAB®; Springer: London, UK, 2009. [Google Scholar]
  5. Cortes, P.; Kazmierkowski, M.; Kennel, R.; Quevedo, D.; Rodriguez, J. Predictive Control in Power Electronics and Drives. IEEE Trans. Ind. Electron. 2008, 55, 4312–4324. [Google Scholar] [CrossRef]
  6. Pérez, J.R. Predictive Control of Power Converters and Electrical Drives; Wiley: Hoboken, NJ, USA, 2012. [Google Scholar]
  7. ONSEMI. IGBT Technologies and Applications Overview: How and When to Use an IGBT. White Paper, 2023. Document ID: TND6235/D. Available online: https://www.onsemi.com/design/resources/technical-documentation#ZHQ9V2hpdGUgUGFwZXJzO3A9MTs= (accessed on 12 December 2023).
  8. Holtz, J. Sensorless control of induction motor drives. Proc. IEEE 2002, 90, 1359–1394. [Google Scholar] [CrossRef]
  9. Wang, G.; Yu, Y.; Zhang, G.; Wang, B.; Yang, M.; Xu, D. Sensorless Control of Motor Drives. In Control of Power Electronic Converters and Systems; Elsevier: Amsterdam, The Netherlands, 2018; pp. 331–369. [Google Scholar] [CrossRef]
  10. Smith, A.N.; Gadoue, S.M.; Finch, J.W. Improved Rotor Flux Estimation at Low Speeds for Torque MRAS-Based Sensorless Induction Motor Drives. IEEE Trans. Energy Convers. 2016, 31, 270–282. [Google Scholar] [CrossRef]
  11. Pal, A.; Das, S.; Chattopadhyay, A.K. An Improved Rotor Flux Space Vector Based MRAS for Field-Oriented Control of Induction Motor Drives. IEEE Trans. Power Electron. 2018, 33, 5131–5141. [Google Scholar] [CrossRef]
  12. Qu, Z.; Hinkkanen, M.; Harnefors, L. Gain Scheduling of a Full-Order Observer for Sensorless Induction Motor Drives. IEEE Trans. Ind. Appl. 2014, 50, 3834–3845. [Google Scholar] [CrossRef]
  13. Zorgani, Y.A.; Koubaa, Y.; Boussak, M. MRAS state estimator for speed sensorless ISFOC induction motor drives with Luenberger load torque estimation. ISA Trans. 2016, 61, 308–317. [Google Scholar] [CrossRef] [PubMed]
  14. Zhao, L.; Huang, J.; Liu, H.; Li, B.; Kong, W. Second-Order Sliding-Mode Observer With Online Parameter Identification for Sensorless Induction Motor Drives. IEEE Trans. Ind. Electron. 2014, 61, 5280–5289. [Google Scholar] [CrossRef]
  15. Comanescu, M. Design and Implementation of a Highly Robust Sensorless Sliding Mode Observer for the Flux Magnitude of the Induction Motor. IEEE Trans. Energy Convers. 2016, 31, 649–657. [Google Scholar] [CrossRef]
  16. Söderström, T. Discrete-Time Stochastic Systems Estimation and Control; Springer: London, UK, 2002; 375p. [Google Scholar]
  17. Bertsekas, D.P.; Tsitsiklis, J.N. Introduction to Probability, 2nd ed.; Athena Scientific: Nashua, NH, USA, 2008; 416p. [Google Scholar]
  18. Stark, H.; Woods, J.W. Probability and Random Processes with Applications to Signal Processing; Pearson Education, Limited: London, UK, 2011; 704p. [Google Scholar]
  19. Jazwinski, A.H. Stochastic Processes and Filtering Theory; Elsevier Science and Technology Books: Amsterdam, The Netherlands, 1970; 376p. [Google Scholar]
  20. Haykin, S. Kalman Filtering and Neural Networks; Wiley-Interscience: Hoboken, NJ, USA, 2001; 304p. [Google Scholar]
  21. Särkkä, S. Bayesian Filtering and Smoothing; Cambridge University Press: Cambridge, UK, 2013; 254p. [Google Scholar]
  22. Anderson, B.D.O. Optimal Filtering; Dover Publications: Mineola, NY, USA, 2004. [Google Scholar]
  23. Simon, D. Optimal State Estimation; Wiley-Interscience: Hoboken, NJ, USA, 2006; 552p. [Google Scholar]
  24. Candy, J.V. Bayesian Signal Processing Classical, Modern, and Particle Filtering Methods; John Wiley & Sons: Hoboken, NJ, USA, 2016; 640p. [Google Scholar]
  25. Chen, F.; Dunnigan, M. Comparative study of a sliding-mode observer and Kalman filters for full state estimation in an induction machine. IEE Proc.-Electr. Power Appl. 2002, 149, 53. [Google Scholar] [CrossRef]
  26. Zahraoui, Y.; Akherraz, M. Kalman Filtering Applied to Induction Motor State Estimation. In Dynamic Data Assimilation-Beating the Uncertainties; IntechOpen: London, UK. [CrossRef]
  27. Singh, K.; Singh, M. Design of Kalman Filter for induction motor drive. In Proceedings of the 2013 Students Conference on Engineering and Systems (SCES), Allahabad, India, 12–14 April 2013; pp. 1–6. [Google Scholar] [CrossRef]
  28. Anderson, B.; Moore, J. Optimal Filtering; Dover Books on Electrical Engineering; Dover Publications: New York, NK, USA, 2012. [Google Scholar]
  29. Kailath, T. Lectures on Wiener and Kalman Filtering; CISM International Centre for Mechanical Sciences Series; Springer: Berlin/Heidelberg, Germany, 1981. [Google Scholar]
  30. Kim, Y.R.; Sul, S.K.; Park, M.H. Speed sensorless vector control of induction motor using extended Kalman filter. IEEE Trans. Ind. Appl. 1994, 30, 1225–1233. [Google Scholar] [CrossRef]
  31. Ramachandiran, G.; Veeranan, S.; Reddy, B. Sensorless Control of Induction Motor with Extended Kalman Filter on TMS320F2812 Processor. Full Pap. Int. J. Recent Trends Eng. 2009, 2, 14. [Google Scholar]
  32. Alonge, F.; D’Ippolito, F. Extended Kalman Filter for sensorless control of induction motors. In Proceedings of the 2010 First Symposium on Sensorless Control for Electrical Drives, Padua, Italy, 9–10 July 2010. [Google Scholar] [CrossRef]
  33. Basheer, O.; Obaid, M. Flux and Speed Estimation of Induction Motors using Extended Kalman Filter. Int. J. Comput. Appl. 2018, 181, 27–31. [Google Scholar] [CrossRef]
  34. Et-Taaj, L.; Boulghasoul, Z.; Elkharki, A.; Kandoussi, Z.; Elbacha, A. Extended Kalman Filter for High performances Sensorless Induction Motor drive. In Proceedings of the 2019 4th World Conference on Complex Systems (WCCS), Ouarzazate, Morocco, 22–25 April 2019. [Google Scholar] [CrossRef]
  35. Akin, B.; Orguner, U.; Ersak, A. State estimation of induction motor using unscented Kalman filter. In Proceedings of the 2003 IEEE Conference on Control Applications, Istanbul, Turkey, 25 June 2003; Volume 2, pp. 915–919. [Google Scholar] [CrossRef]
  36. Kumar, S.; Prakash, J.; Kanagasabapathy, P. A critical evaluation and experimental verification of Extended Kalman Filter, Unscented Kalman Filter and Neural State Filter for state estimation of three phase induction motor. Appl. Soft Comput. 2011, 11, 3199–3208. [Google Scholar] [CrossRef]
  37. Jafarzadeh, S.; Lascu, C.; Fadali, M.S. State Estimation of Induction Motor Drives Using the Unscented Kalman Filter. IEEE Trans. Ind. Electron. 2012, 59, 4207–4216. [Google Scholar] [CrossRef]
  38. Dominguez, P.A.; Silva, C.A.; Yuz, J.I. State and Resistance Estimation in Sensorless FOC Induction Motor Drive Using a Reduced Order Unscented Kalman Filter. In Proceedings of the 2012 VI Andean Region International Conference, Cuenca, Ecuador, 7–9 November 2012. [Google Scholar] [CrossRef]
  39. Rigatos, G.; Siano, P. Sensorless nonlinear control of induction motors using Unscented Kalman Filtering. In Proceedings of the IECON 2012-38th Annual Conference on IEEE Industrial Electronics Society, Montreal, QC, Canada, 25–28 October 2012. [Google Scholar] [CrossRef]
  40. Yildiz, R.; Barut, M.; Zerdali, E. A Comprehensive Comparison of Extended and Unscented Kalman Filters for Speed-Sensorless Control Applications of Induction Motors. IEEE Trans. Ind. Inform. 2020, 16, 6423–6432. [Google Scholar] [CrossRef]
  41. Yan, L.; Wang, F.; Dou, M.; Zhang, Z.; Kennel, R.; Rodríguez, J. Active Disturbance-Rejection-Based Speed Control in Model Predictive Control for Induction Machines. IEEE Trans. Ind. Electron. 2020, 67, 2574–2584. [Google Scholar] [CrossRef]
  42. Yan, L.; Wang, F. Observer-Predictor-Based Predictive Torque Control of Induction Machine for Robustness Improvement. IEEE Trans. Power Electron. 2021, 36, 9477–9486. [Google Scholar] [CrossRef]
  43. Liu, S.; Wang, Z.; Chen, Y.; Wei, G. Protocol-Based Unscented Kalman Filtering in the Presence of Stochastic Uncertainties. IEEE Trans. Autom. Control. 2020, 65, 1303–1309. [Google Scholar] [CrossRef]
  44. Kharab, A.; Guenther, R. An Introduction to Numerical Methods, 4th ed.; CRC Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  45. Mathworks. Simulink User’s Guide, 23.2 ed.; The Mathworks Inc.: Natick, MA, USA, 2023; Available online: https://www.mathworks.com/help/pdf_doc/simulink/simulink_ug.pdf (accessed on 12 December 2023).
  46. Plexim GmbH. PLECS—The Simulation Platform for Power Electronic Systems, 4.7 ed.; Plexim GmbH: Zurich, Switzerland, 2023; Available online: https://www.plexim.com/files/plecsmanual.pdf (accessed on 12 December 2023).
  47. El-Sharkawi, M. Fundamentals of Electric Drives (Electrical Engineering); Brooks/Cole: San Francisco, CA, USA, 2000; 336p. [Google Scholar]
  48. Silva, C.A.; Yuz, J.I. On sampled-data models for model predictive control. In Proceedings of the IECON 2010-36th Annual Conference on IEEE Industrial Electronics Society, Glendale, AZ, USA, 7–10 November 2010. [Google Scholar] [CrossRef]
  49. Julier, S.; Uhlmann, J.K.; Durrant-Whyte, H. A new approach for filtering nonlinear systems. In Proceedings of the 1995 American Control Conference-ACC’95, Seattle, WA, USA, 21–23 June 1995; Volume 3, pp. 1628–1632. [Google Scholar]
Figure 1. Block diagrams of IM continuous-time model. (a) α-β Stator currents dynamic. (b) α-β Rotor flux linkages dynamic. (c) Rotor shaft angular speed dynamic. (d) Load torque dynamic.
Figure 1. Block diagrams of IM continuous-time model. (a) α-β Stator currents dynamic. (b) α-β Rotor flux linkages dynamic. (c) Rotor shaft angular speed dynamic. (d) Load torque dynamic.
Mathematics 12 00373 g001
Figure 2. No-load machine. (a) Stator α-β current. (b) Rotor α-β flux linkage.
Figure 2. No-load machine. (a) Stator α-β current. (b) Rotor α-β flux linkage.
Mathematics 12 00373 g002
Figure 3. No-load machine. (a) Rotor shaft speed. (b) Load torque.
Figure 3. No-load machine. (a) Rotor shaft speed. (b) Load torque.
Mathematics 12 00373 g003
Figure 4. Torque impact. (a) Stator α-β current. (b) Rotor α-β flux linkage.
Figure 4. Torque impact. (a) Stator α-β current. (b) Rotor α-β flux linkage.
Mathematics 12 00373 g004
Figure 5. Torque impact. (a) Rotor shaft speed. (b) Load torque.
Figure 5. Torque impact. (a) Rotor shaft speed. (b) Load torque.
Mathematics 12 00373 g005
Figure 6. Variation of RMSE as a function of sampling time. (a) α-β Stator current RMSE. (b) α-β Rotor flux linkage RMSE. (c) Rotor shaft speed RMSE. (d) Load torque RMSE.
Figure 6. Variation of RMSE as a function of sampling time. (a) α-β Stator current RMSE. (b) α-β Rotor flux linkage RMSE. (c) Rotor shaft speed RMSE. (d) Load torque RMSE.
Mathematics 12 00373 g006
Table 1. Typical load torque functions [47].
Table 1. Typical load torque functions [47].
Load TorqueModel
Constant T l = T l 0
Linear T l = K ω r
Quadratic T l = K ω r 2
Inverse T l = K ω r 1
Table 2. Models’ root mean squared error.
Table 2. Models’ root mean squared error.
StateEulerTaylorRK 2RK 4
i s α 2.32880.37430.58300.4188
i s β 2.32860.37230.59850.4177
ψ r α 0.05670.00910.02450.0191
ψ r β 0.05670.00890.02860.0190
ω r 21.691411.31171.99970.1401
T l 0.00460.00467.0356 × 10 5 7.0171 × 10 9
Table 3. EKF average root mean squared error.
Table 3. EKF average root mean squared error.
StateEKF
EulerTaylorRK 2RK 4
i s α 0.36120.19770.20290.2026
i s β 0.35770.19670.20170.2013
ψ r α 0.07770.03770.04330.0433
ψ r β 0.07840.03790.04560.0456
ω r 28.406327.210124.276224.5003
T l 0.10380.10380.10420.1042
Table 4. UKF average root mean squared error.
Table 4. UKF average root mean squared error.
StateUKF
EulerTaylorRK 2RK 4
i s α 0.36110.19780.20290.2026
i s β 0.35750.19660.20160.2012
ψ r α 0.07770.04120.04310.0429
ψ r β 0.07840.04250.04410.0443
ω r 28.798228.030724.699224.8631
T l 0.10380.10380.10420.1042
Table 5. EKF and EKF average execution times per time step.
Table 5. EKF and EKF average execution times per time step.
MethodEKF Execution TimeUKF Execution Time
Euler6.9456 × 10 6 s9.9939 × 10 5 s
Taylor7.4185 × 10 6 s1.3481 × 10 4 s
RK 28.0369 × 10 6 s1.1440 × 10 4 s
RK 49.7486 × 10 6 s1.4439 × 10 4 s
Table 6. Max absolute error EKF at startup.
Table 6. Max absolute error EKF at startup.
StateEKF
EulerTaylorRK 2RK 4
i s α 3.16733.16733.16733.1673
i s β 1.32791.32791.32801.3280
ψ r α 6.32246.37656.67736.6699
ψ r β 5.89925.89995.90025.9041
ω r 103.6104.4497.95198.636
T l 0.38060.38040.38050.3805
Table 7. Max absolute error UKF at startup.
Table 7. Max absolute error UKF at startup.
StateUKF
EulerTaylorRK 2RK 4
i s α 3.16733.16733.16733.1673
i s β 1.32791.32791.32801.3239
ψ r α 6.31666.32116.67586.4512
ψ r β 5.89925.90005.90015.8696
ω r 103.92105.1798.55499.191
T l 0.38060.38040.38050.3805
Table 8. Max absolute error EKF after starting.
Table 8. Max absolute error EKF after starting.
StateEKF
EulerTaylorRK 2RK 4
i s α 1.44411.03301.03951.0388
i s β 1.48681.05031.04831.0443
ψ r α 0.12870.04130.06820.0683
ψ r β 0.12970.03770.06200.0621
ω r 23.81811.50714.46513.498
T l 0.55160.55190.55180.5518
Table 9. Max absolute error UKF after starting.
Table 9. Max absolute error UKF after starting.
StateUKF
EulerTaylorRK 2RK 4
i s α 1.44401.03561.03950.7169
i s β 1.48661.05141.04840.5494
ψ r α 0.12860.04060.06820.0344
ψ r β 0.12960.03570.06200.0049
ω r 23.81211.78414.4593.116
T l 0.55160.55190.55180.0867
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Herrera Hernández, R.; Reusser, C.; Coronel, M.; Carvajal, R. On Accurate Discrete-Time Dynamic Models of an Induction Machine. Mathematics 2024, 12, 373. https://doi.org/10.3390/math12030373

AMA Style

Herrera Hernández R, Reusser C, Coronel M, Carvajal R. On Accurate Discrete-Time Dynamic Models of an Induction Machine. Mathematics. 2024; 12(3):373. https://doi.org/10.3390/math12030373

Chicago/Turabian Style

Herrera Hernández, Ramón, Carlos Reusser, María Coronel, and Rodrigo Carvajal. 2024. "On Accurate Discrete-Time Dynamic Models of an Induction Machine" Mathematics 12, no. 3: 373. https://doi.org/10.3390/math12030373

APA Style

Herrera Hernández, R., Reusser, C., Coronel, M., & Carvajal, R. (2024). On Accurate Discrete-Time Dynamic Models of an Induction Machine. Mathematics, 12(3), 373. https://doi.org/10.3390/math12030373

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

Article Metrics

Back to TopTop