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
where
is the stator voltage space vector on the stationary reference frame
α-
β; in the same way,
and
are the stator current and flux linkage, respectively, on the
α-
β reference frame. The unit space vector
denotes the position of the coil on the space. The subscript
denotes the three-phase component.
The stator electrical equation on the
α-
β reference frame is
where
is the stator resistance.
Similarly, the rotor electrical equation is
where
and
are the rotor current and flux linkage, respectively, on the
α-
β reference frame;
is the rotor resistance;
p is the number of pole pairs; and
is the rotor shaft angular speed.
The flux linkages of the stator and rotor are
where
and
are the stator and rotor self-inductances, respectively, while
is the mutual inductance between stator and rotor.
On the other hand, the mechanical equation at the rotor shaft is
where
J is the inertia;
is the machine developed torque; and
represents the load torque. The machine developed torque is
where the superscript
denotes the complex conjugate and
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
, the rotor flux linkage
, and the rotor shaft speed
. 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
where
a and
b coefficients are
The above coefficients are in terms of already-named machine parameters and additionally by
where
is known as the stator dynamic transient resistance;
as the leakage coefficient; and
as the rotor time constant. The function
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 (
); finally,
Figure 1d illustrates the dynamic behavior of the load torque.
In compact form, the representation of the IM model in state-space is
Equation (16) is given by (10)–(15);
is the state vector and
is the input vector, in this case
. The output (17) is related to the state vector by the matrix
. Since it is a sensorless drive, only the stator current
is measured, hence
Note that the state variables depend implicitly on time t (i.e., ). 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
where
and
;
is the time interval between samples. The input
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
The IM model (16) approximated by the Euler Method is
For the Euler method, the local truncation error (i.e., the error committed at each time step) is of the order .
3.2. Taylor Method
The derivative of a function is approximated using Taylor series expansions
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
. As simply as possible, the discretization of relative degree
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
is given by
For a second-order Taylor series expansion, the local truncation error is .
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 .
3.3.1. Runge–Kutta 2
Runge–Kutta 2 is a second-order approximation. This method requires two evaluations of
at each time step. The IM model (16) is approximated by
where
and
. Note that
estimation requires future input values. Since the system is assumed to be causal and future inputs are unknown,
is used instead. For the IM,
and
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
in one sample period. The approximation of IM model (16) by RK 4 is
where
Since the input is ZOH and no future input values are taken into account,
is used in
,
, and
. For the IM, the estimate points
,
,
, and
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,
μs and a 4 kW IM was used; its nominal parameters are
380 V;
8.6 A;
50 Hz; two pole pairs;
1.32
;
2.63
;
0.1889 H;
0.1972 H;
0.2012 H; and
0.528 kg·m
2.
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
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
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
increases. Thus, using Euler for approximating the IM model is not advisable when the
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
where
w and
v are Gaussian noise with zero mean and covariance matrix
and
, respectively. The term
w is associated with model uncertainties, while
v 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
, which incorporates all the statistical information about the state
that is contained in the available observation
and the initial condition
[
19,
21]. Using the Bayes theorem, the joint pdf is
where the denominator
is a normalization constant. It is necessary to mention that
and
are the time series of the state vector
and of the measurements
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
, 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
.
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
Second: update. Given the measurement
at time
k, the posterior distribution is calculated by solving
where
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
is the state estimate and
is its covariance matrix. The subscript
k denotes the posterior and
the prior.
is the Kalman gain;
is the identity matrix; and
is the Jacobian matrix of
evaluated at
.
Algorithm 1: EKF Algorithm |
|
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
SP are required, with
. The SPs are computed as follows
where
is a scale factor that modifies the dispersion of the SP around the mean. The term
is computed using Cholesky factorization. Also, weight factors
and
are associated to the mean and covariance, defined as
where typically,
for Gaussian distributions;
;
with
.
The steps to implement the UKF are summarized in Algorithm 2.
is the covariance matrix of the measurement;
is the cross-covariance between the state and the measurement; and
is the filter gain.
Algorithm 2: UKF algorithm |
|
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
μs; the noise covariance matrices are
and
. The SP coefficients were selected empirically as in [
40], where
,
, and
. 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 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.