Adaptive Output Feedback Control for the Trajectory Tracking of High-Speed Trains with Disturbance Uncertainties on the Basis of Neural Network Observers

)e dynamic model of high-speed trains (HSTs) is nonlinear and uncertain; hence, with the decrease in the running interval of HSTs, an accurate and safe train operation control algorithm is required. In this study, an adaptive output feedback trajectory tracking control method for HSTs is proposed on the basis of neural network observers. )e proposed method aims to solve problems, such as the immeasurable speed, model parameter disturbance, and unknown external disturbance of HSTs. In this method, a neural network adaptive observer is designed to estimate the velocity of an HST. Another neural network model is used to approximate the model uncertainties. Moreover, a robust controller is constructed by considering the train position and velocity tracking errors. In the proposed observer/controller, the bound function of estimator errors is introduced to increase the accuracy and safety of the tracking system. Furthermore, the adaptive update value of the neural networks, output weights, and bound function are performed online. All adaptive algorithms and the observer/controller are synthesized in nonlinear control systems. )e error signals of the closed-loop trajectory tracking system are uniformly and eventually bounded through a formal proof on the basis of the Lyapunov methods. Simulation examples illustrate that the proposed controller is robust and has excellent tracking accuracy for system model parameter and external disturbance.


Introduction
Given their large capacity, low energy consumption, and low pollution, high-speed trains (HSTs) have become one of the most efficient modes in railway transportation, thereby attracting much attention toward the development of green transport [1]. However, the running process of HSTs is coupled with a complex system and a complex environment. In addition, such process is nonlinear, multiobjective, multiconstrained, and time varying. In summary, the dynamic model of HSTs is nonlinear and uncertain [2][3][4]. A primary task of automatic train operation (ATO) is speed regulation. A speed regulation control strategy based on a precise model cannot easily satisfy the requirements of safe and punctual operation. With the increase in train speed and the decrease in the running interval of HSTs, an accurate and reliable control strategy that uses an adaptive and intelligent method should be explored to guarantee a safe and comfortable operation.
Proportion integration differentiation (PID) control is widely used in train speed control to ensure the safety of train operation [5][6][7]. e PID control algorithm is simple and easy to implement and does not easily make mistakes. However, it easily causes the instability and energy consumption of train operation due to the excessive frequency of adjustment. Moreover, its parameter selection has no clear method.
us, an optimal state of train operation control is difficult to achieve.
To solve the above problem, Sekine et al. [5] proposed a two-degree-of-freedom fuzzy neural network control method, which has a hierarchical structure of two sets of control knowledge. In [6,7], an adaptive fuzzy PID soft switching control method, which switches between the fuzzy and PID controllers, was first proposed for use in ATO systems. is method improved the smoothness, comfort, and parking accuracy of train operation. However, the above methods were based on the subway as research object or on ordinary electric multiple units with low running speed; such bases are unsuitable for HSTs because the nonlinear characteristics of the train model are evident in high-speed operation. For this reason, the authors in [8] used a datadriven modeling method to establish multiple local linear models to describe the dynamic characteristics of the HST running process. However, the basic contradiction of the linearization method in controlling the nonlinearly controlled object is difficult to avoid.
To enhance the processing of the nonlinear components of trains and achieve an improved accuracy for the speed tracking control of HSTs, the authors in [9] used the Elman neural network. e dynamic model of HSTs was established using substantial offline data, but the timevarying model and external disturbance in the process of train operation were not considered. Li et al. [10] divided the train system into linear and nonlinear components and used the least squares method and adaptive network-based fuzzy inference system (ANFIS) methods to model the linear and nonlinear components, respectively. However, the stability of the entire closed-loop system was not discussed because the controller was not designed in a unified framework.
In [11], a direct fuzzy controller was designed for onemass-point and multimass-point models, but the accurate model parameters are necessary and do not consider external disturbances. Combining H 2 and H ∞ methods, the authors in [12] designed an HST cruise controller using accurate model parameters, thereby enhancing the robustness of the system. However, the stability analysis of the closed-loop system was based on linearization theory. us, the overshoot of tracking error was relatively large. In consideration of the nonlinear characteristic and uncertain parameters of HSTs, a model predictive controlmixed H ∞ algorithm was designed in [13], thereby achieving velocity tracking control in the case of model disturbance. However, only the matching uncertainties were considered. In [14], a sliding mode adaptive robust controller for HTS was designed on the basis of Lyapunov stability theory. However, no adaptive estimation or robust terms were available for uncertain model parameters and external disturbance. As a result, the convergence of a closed-loop system depended on the design parameters, thus leading to excessive steady-state errors. In [15], a robust adaptive controller based on Lyapunov stability analysis was designed for a multimass-point model and considered external disturbance and parameter uncertainty. In the case of HSTs' repeated running on tracks, an iterative algorithm was studied in the tracking control [16][17][18][19]. However, iterative learning required the motion of HST to be repetitive; that is, the reference speed curve of the train must be unchanged. us, the real-time performance of the iterative method was poor, and its convergence speed must be further improved [20].
With the increase in the speed of HST and the development of intelligent transportation, moving block technology must be introduced to improve operational efficiency.
us, an optimized safety tracking distance is necessary. In this case, the following must be considered: the strong nonlinearity and the uncertainty of the dynamic model of HSTs, the uncertain mismatching under strong wind conditions, the robustness of the controller, and the accuracy of the closed-loop system. e combination of nonlinear and intelligent methods is a good solution to the above problems. In [21], a trajectory tracking controller of wheeled mobile manipulators was designed on the basis of the fuzzy neural network and extended Kalman filtering; as a result, the tracking problem of wheeled robot under uncertainty is solved. In [22], an adaptive neural network controller of a fully actuated marine surface vessel was proposed. Based on above, the neural network and integral sliding mode method were combined by [23] to design a tracking controller for fully driven surface ships.
In [24,25], the sliding mode and neural network, the back stepping, and the fuzzy algorithm were used to design a tracking controller for aircraft. In [26], an adaptive sliding mode control system with a double-loop recurrent neural network structure was proposed and applied to the motion control of a gyroscope. In [27], an adaptive high-gain observer was used to estimate the unmeasured state for a nonlinear second-order system with external disturbance. In [28,29], an adaptive fuzzy observer was designed on the basis of stability theory. In [30,31], a neural network observer was designed for an uncertain strict feedback system and a nonlinear stochastic system, respectively. In [32], the output feedback control of wheeled mobile robot formation was designed by combining neural network observers. In [33,34], neural network observers were used to design a trajectory tracking controller for underactuated vehicles. In [35], a neural network observer was designed for noninteger order systems.
Researchers have exerted much effort in identifying and determining the value of a model parameter. e resistance of the HST is related to its operating environment, and the resistance parameter of different HST types also varies. Moreover, obtaining an accurate parameter is difficult, and the train speed variation has a strong nonlinear relationship with the resistance it receives. In addition, the sensor measurement error or sensor failure hinders train operation safety. ese factors restrict the increase in train speed. us, we focus on the following: nonlinear characteristics of the HST model, time-variable model parameters, unknown disturbance of the external environment, and immeasurable HST speed. Motivated by these observations, this study proposes an adaptive neural observer-(ANO-) based output feedback controller for the trajectory tracking control of HSTs. e major contributions of this study can be summarized as follows: (1) References [2][3][4][5][6][7][8][9][10][11][12] only dealt with the speed tracking problem of HSTs. In consideration of the demand for multitrain tracking control of the moving block and 2 Journal of Advanced Transportation in combination with the RBF adaptive and sliding mode methods, a displacement and a speed tracking control strategy for HSTs are designed. e sliding mode method is used to ensure the stability of the tracking control. An adaptive neural network is used to estimate the nonlinearity and uncertainty of the HST to ensure the accuracy of the tracking control. A bounded estimation function of design error is used to enhance the robustness of the system. (2) e tracking control algorithm does not require the running speed of an HST. On the basis of [21][22][23][24][25][26][27][28][29][30][31][32], an adaptive neural network observer is designed under the condition of unknown dynamic model parameters and external interference of HST. us, the HST equipment and adaptive neural observer can be reconstructed into pure feedback form.
(3) e design parameters of the proposed observer/ controller are adjusted online, and the stability and convergence of the closed-loop system are proven through a formal proof based on the Lyapunov methods.
is paper is organized as follows. Section 2 introduces the problem formulation. An ANO is designed in Section 3. Section 4 presents the proposed controller. In Section 5, an analysis of stability theories is discussed. In Section 6, the main result for trajectory tracking is presented in the form of a theorem, and a rigorous mathematical explanation is provided. Conclusions are presented in Section 7.

Problem Formulation
With only the longitudinal motion of the train considered, the dynamic model of an HST can be expressed as follows [36]: where t represents a continuous time indicator; s(t) (m) and v(t) (m/s) are the operating displacement and speed of the train, respectively; M > 0 is the equivalent total mass of the train system, including the weights of the train, passenger, and baggage; F(N) is the resultant force of the train; u(t) (N) is the control input of the train, that is, the traction or braking force of the train; f d (v, s) (N) represents the total resistance to the train operation; f a (s) (N) is the additional resistance to the train operation; and f b (v) (N) is the general resistance of the train operation. f a (s) is a resistance other than the general resistance, which mainly includes slope additional resistance (as shown in Figure 1), curve resistance, and tunnel resistance.where L s is the tunnel length.
(1) Slope additional resistance: as shown in Figure 1, the slope additional resistance force f s comes from the grade of the track at the point of the train's location and is modeled as where α(t) is the slope angle of the current track. e track line comprises several horizontal, slope, and curvature tracks. For a certain slope track, the slope angle θ is a constant. (2) Curve resistance: the curvature force f c , which comes from the curvature of the track at the point of the train's location, is modeled as where φ(t) is the center angle of the curve and L r (t) is the curve length. (3) Tunnel resistance: f a is obtained as follows: e curve or tunnel can also be converted into an equivalent ramp to obtain the slope angle. Equation (5) can be changed as f a (s) � l(t)sin(r(s)), (6) where l(t) is the time-varying additional resistance parameter and r(s) is the time-varying equivalent slope angle of the train displacement. Given that the displacement is the integral of the velocity, equation (6) can be changed to where t 1 represents the start time of the additional resistance of the train and t 2 represents the end time of the additional resistance of the train. General resistance consists of mechanical and air resistance and is expressed as where a 0 (t), a 1 (t), and a 2 (t) are the general resistance parameters of the train with time variation.
Journal of Advanced Transportation erefore, the HST motion model adopted in this study can be written as follows: where d(t) represents unmolded dynamics, such as external disturbance. e common system assumptions in nonlinear and neural network controls are provided to discuss the convergence of the algorithm. Assumption 3. Train parameters M, a 0 (t), a 1 (t), and a 2 (t) are unknown but continuously bounded, and they are greater than zero. e time-varying additional resistance parameter l(t) and external disturbance d(t) are unknown and continuously bounded.
e ideal NN weights are bounded such that ‖W * ‖ ≤ ‖W N ‖, and ‖W N ‖ is the boundary of ‖W * ‖.

Assumption 5.
e nonlinear term f d in equation (9) satisfies the following Lipchitz condition: where L i ∈ R + � 1, 2 is unknown and constant. e control objective is that, for the mathematical model (9) of HST and under Assumptions 1-3, the parameters and disturbances of the HST model are unknown according to the train position information. e output feedback control law is designed to make the HST run along the expected trajectory and ensure that the closed-loop signal is consistent and ultimately bounded to achieve HST tracking.

Adaptive Neural Network Observer Design
e position of a train can be given by GPS and other equipment. Meanwhile, the speed information of the train is measured using a speed sensor, which may have measurement errors or sensor failure. erefore, this section focuses on estimating the speed information in accordance with the position information of the HST.
First, the observation error is defined: where s is the estimation of position; v is the estimation of speed; s and v are the estimation errors of s and v, respectively. e adaptive update laws of the observer are designed as follows: where f is the estimate of f, D 1 is a robust term whose specific expression is given in the subsequent analysis, and k 1 is positive and a design parameter. e time derivative of equation (11) is determined and integrated in equations (9) and (12) to obtain the following: _ e parameters a 0 (t), a 1 (t), a 2 (t), and l(t) in equation (12) are based on Assumption 3 and are continuously bounded variables. us, the function f(z 1 ): Ω ⟶ R is bounded and can be approximated using the neutral network given by the following: where where c f,j � [c f,j,1 , c f,j,2 , c f,j,3 ] ∈ R 3 is the center value of the Gaussian function, σ f,j is the ith hidden layer node, and Moreover, the velocity vector is unknown. erefore, the following neural network is used for estimation: where z 1 � [s, tvn, qu] is the input of the neural network and W f represents the estimated value of W * f used to observe v in equation (12).
Substituting equations (13), (15), and (16)into equation (14), we obtain the closed-loop error dynamic equations as where the estimation error , the neutral network weights following the updating law are as follows: where Γ f ∈ R m×m and is positive; σ f is the design parameter and is positive; and W f0 is the W * f a priori estimate. A bounded function δ � ε f (z 1 ) exists for ε f (z), and δ is the estimation of δ. e updating law of δ is designed, and the robust term D 1 is defined as follows: e adaptive updating law is designed as follows: where c( _ s + tk 1 ns) � tanh(( _ s + tk 1 ns)/κ 1 ) and κ 1 > 0 is a design constant. [37]. In addition, the function c can also be selected as c( _ s + tk 1 n s) � ( _ s + tk 1 ns)/4κ 1 .

Remark 2.
If we select the estimated adaptive updating law as equations (12), (19), and (21), then the estimation error (s, tv) can be set to approximately zero.

Output Feedback Controller Design
In this section, the sliding mode method is introduced to design the tracking controller of an HST to reference displacement and speed. e output value of the designed speed observer is used because the speed is immeasurable. erefore, the displacement tracking error and the speed tracking error are obtained as follows: e filter error is designed as follows: e time derivative of equation (21) can be obtained as where k 4 is a positive designed parameter. Equation (24) by M/χ is used, and (12) is substituted into (24) to obtain the following: Considering Journal of Advanced Transportation 5 where According to Assumptions 2 and 3, the nonlinear function f v is unknown and bounded. us, the RBF neural network is introduced to approximate f v as follows: where q is the number of neurons in the hidden layer; ε v (z 2 ) is the neural network's forced reconstruction error and ε v,N is the boundary value of ε v (z 2 ); where c v,j � [c v,j,1 , c v,j,2, , . . . , c v,j,6 ] ∈ R 6 is the center value of the Gaussian function, σ v,j is the ith hidden layer node, and j � 1, 2, . . ., m. R n denotes n-dimensional Euclidean space. However, an ideal neural network does not exist, and the velocity vector is unknown; thus, the following neural network is used for estimation: where W v is the estimated value of W * v . In accordance with the position information of the HST and the observed speed information, the tracking output feedback control law of the HST is designed as where k 5 is a positive designed parameter and D 2 is a robust term provided in the subsequent analysis. e neutral network weights use the following updating law: where Γ v ∈ R q×q and is positive and σ v is a positive designed parameter. Substituting (28), (30), and (31) into (25), we can obtain the closed-loop dynamic equations as where is defined; it is composed of neural network approximation error and external environmental disturbance. e robust term D 2 can be given as Similar to the observer, the σ − correction [38] is used to design the adaptive update law of θ as where the function c(e v ) � tanh(e v /κ 2 ) and κ 2 is a positive designed parameter; α 2 and β 2 are positive designed parameters; and θ is an estimate of θ.
Remark 3. If we select the controller (31) and the update adaptive law as (32) and (35), then the estimation error (e s , e v ) can be set to approximately zero.

Stability Analysis
Theorem 1. e nonlinear model of the HST (9) and adaptive neutral network (12) are considered. Under Assumptions 1-5, if we select the neutral network weight updating law (19) and the robust terms with its law (20), then all the signals in the closed-loop systems (13) and (14) are uniformly ultimately bounded.

Theorem 2.
e general HST model, which is given by equation (9), is considered. Under the conditions of Assumptions 1-5, if the control law is as (31), the neutral network weight updating laws are as (19) and (32), the robust terms law is as (34), and the adaptive update law is as (35); then, all the signals in the closed-loop systems (33) are uniformly ultimately bounded.
Proof of eorem 2. e Lyapunov function is considered as follows: Derivative V 2 provides the following:

Journal of Advanced Transportation
We select W v � W v − W * v and θ � θ − θ in accordance with Assumptions 3 and 4 and the definition of θ; thus, equation (50) can be written as Substituting (26), (32), and (35), into (51), we obtain the dynamic equations of V 1 as Substituting θ, D 2 , θ and equation (44) into equation (51), we can obtain According to Lemmas 1 and 2, we obtain some inequalities as en, equation (53) can be written as (55) Substituting (56) into (55), we can obtain the following: Defining φ 2 � (μ 2 /ρ 2 ) and integrating (58) over [0, t], Equation (59) shows that V 1 (t) is bounded because the value of ρ 1 can be decreased. As a result, the exponential convergence of sliding surfaces is ensured. Hence, all the closed-loop control signals are UUB.
Using equation (59), we obtain the Euclidean norm as (60) erefore, the error variable [s, v, can be made arbitrarily small as t ⟶ ∞ (Figure 2).

Simulation and Results
In this section, simulation studies are conducted for the designed neural network observer and the sliding mode neural network controller. e CHR380AL HST operation control environment is simulated in MATLAB (R2017b) software. e total nominal mass of the train is M � 890 tons, the continuous operation time is 350 km/h, the maximum operation speed is 380 km/h, χ � 0.9, and the basic general resistance is obtained as [37] f b ′ � 0.53 + 0.00388v + 0.000114v 2 . (61) Furthermore, the general resistance can be written as 8 Journal of Advanced Transportation a 0 � 4717 + 471.7 sin(0.00047t), a 1 � 61.5 + 6.15 sin(0.00047t), a 2 � 1.05 + 0.105 sin(0.00047t).
e selected railway track is 561.37 km, with one ramp, one curve, and one tunnel length. e actual operation environment conditions of the train in the computer simulation are shown in Table 1. Such conditions lead to the changes of the additional resistance coefficients in the travel, as shown in Figure 3. en, the two limited speed ranges are shown in Figure 4 [37]. e steps of additional resistance are calculated as follows.  Table 1. us, equation (5) can be written as where g is the acceleration of gravity. Substituting g � 10 m/s 2 into (60), we obtain f t (t) as 0.0295 N/kg (Table 1) Table 1. We select the curve length as whereR is the curve radius. Substituting (61) into (3), we can obtain the following: Substituting R � 200m into equation (62), we obtain f c (t) as 0.0295 N/kg (Table 1).
No restricted ramp exists in the tunnel; thus, equation (4) can be written as Substituting R � 1000m into equation (62), we obtain f t (t) as 0.0295 N/kg (Table 1).
Substituting L s � 1000 m into equation (63), we obtain f t (t) as 0.0013 N/kg.
Substituting f s (t) � 0, f t (t) � 0.0013 N/kg, and f c (t) � 0.0295 N/kg into equation (5), we obtain f a (t) as 0.0308 N/kg. e disturbance is considered as d(t) � (10 5 /M)(sin(0.08t) + cos(0.08t)). e desired trajectory is optimized by the genetic algorithm in [38], but the algorithm is optimized offline. If the unknown disturbance of the model parameters and external environment occur, then the open-loop control method causes a tracking error. us, we provide the method in this study. e designed parameters are given as σ f,j � 2900M, e other designed parameters are shown in Table 2. k 1 , k 2 , and k 3 are the observation parameters; k 4 and k 5 are the tracking control parameters.

CRH380AL HST equation (30) NN-1 equation (30)
Sliding mode surface equation (23) NN weight update rule equation (32) Robust terms equation (34) Observer equation (12) Controller equation (31) Neural network equation (17) Nonlinear observer NN weight update rule equation (19) Robust terms equation (20) [ Adaptive update rule equation (35) Adaptive update rule equation ( To verify the effectiveness of the proposed algorithm and the influence of different designed parameters for HST, we propose three simulation cases, as follows: Case 1: the control parameters are unchanged; the value of estimated parameters is increased Case 2: the estimated parameters are unchanged; the value of control parameters is increased; all the designed parameters are increased Case 3: the designed parameters are unchanged; the environmental disturbance is changed In Case 1, we select the parameters from Values 1 to 5, as listed in Table 2. In Case 2, we select the parameters from Values 5 to 9. In Case 3, we select the parameter from Value 10. Figures 5-7 are the results of Case 1. e results in Case 2 are shown in Figures 8-12. e values in Table 3 are the energy consumption and RMS for the corresponding parameters. In Case 3, we select three types of disturbance, namely, no disturbance, fast time-varying disturbance, and slow time-varying disturbance. e parameters are as follows:k 1 � 140, k 2 � 35, k 3 � 175, k 4 � 28, and k 5 � 3115000. e energy consumption and RMS for different environmental disturbances are shown in Table 4.
As shown in Figure 5, s − s hat indicates the position observer estimation error s; as shown in Figure 6, v − v hat indicates the velocity observer estimation error v. In Case 1, the observation parameters are small.       In Case 1, Value 1 is the smallest observation parameter, and Value 5 is the largest. e value of the parameters between them increases from small to large. Figures 5 and 6 show that when the observation parameters are small, the fluctuation of the position estimation error and velocity estimation error is large in the beginning. e estimation error decreases with the increase in the observation parameters. Table 3 shows that when Value 5 is selected, the RMS of the estimation error is reduced to s � 0.6257 and v � 0.5132. Energy consumption is the total energy consumption during HST operation. In Figure 7, the control input curve of Value 1 is jittery in the beginning, and the other curves almost overlap, as shown in Table 3. e energy consumptions of Values 2, 3, 4, and 5 are almost the same and smaller than that of Value 1. erefore, the comparison of energy consumption, position estimation error, and velocity estimation error shows that a larger selection of observation parameters is better than a smaller one. e RMS value of the velocity tracking error v − v d is unchanged, but the RMS value of the position error s − s d is increased by 35.4%. erefore, the selection of parameters requires further analysis.
In Case 2, we select the designed parameters, that is, Values 5 and 6, which have the same estimated parameters and different tracking control parameters. We also consider Values 7 and 8 for comparison. Values 6, 8, and 9 have the same control parameters and different estimated parameters; this case is similar to Case 1. In Figures 8 and 9, the s and v curves for Values 5 and 6 almost overlap. In Table 3, the s RMS values of Values 5 and 6 are 0.6257 and 0.6022, respectively. Observing the v curve, the RMS values of Values 5 and 6 are 0.5132 and 0.4957, respectively. In Figure 10, the position error s − s d of Value 5 is much larger than that of Value 6. In Figure 11, the speed tracking error v − v d of Value 5 is also much larger than that of Value 6. In Table 3 Figure 10 shows that the position tracking performance decreases slightly with the increase in the estimated parameters. However, it is also much smaller than Values 5 and 7. e v − v d curve of Value 6 in Figure 11 still has a slight vibration.
In summary, increasing the value of control design parameters can reduce the tracking error of position and speed. Furthermore, increasing the value of observation estimation design parameters can reduce the speed estimation error and reduce the jitter of speed tracking error.
On the basis of the above analysis, the designed parameters of Value 10 are selected. e RMS values of s, v, s − s d , and v − v d in Table 3 are greatly reduced. However, as the design parameters start to increase from Value 9, energy consumption also starts to increase. e users can select parameters in accordance with their requirements.
In Case 3, we define three types of environmental disturbances to verify the robustness of the proposed algorithm, as follows: (1) No disturbance: (2) Fast time-varying disturbance: (3) Slow time-varying disturbance:   Journal of Advanced Transportation 13        Figure 15 shows the speed estimation error in the range of [− 0.03, 0.01](m/s). Figure 16 shows the position tracking error in the range of [− 0.03, 0.01](m/s). Figure 17 shows the speed tracking error in the [− 0.06, 0.04] (m/s) range. Figure 18 shows the control output in the [− 2 × 10 6 , 5 × 10 6 ] (N) range. Table 4 shows the RMS values of each error under three types of environmental disturbances. Table 4 shows that the RMS values increase slightly with the increase in disturbance value. As shown in Table 4, the total energy consumption reveals that HSTs need to consume relatively more energy to improve tracking performance under large disturbance.
e results show that the system can also obtain good tracking performance under large disturbance.

Conclusion
An output feedback tracking control for HSTs is designed in this study. In view of the uncertainty of the train dynamic model, the disturbance of the environment, and the immeasurable speed state, a neural network-based speed observer and a sliding mode-based speed and position tracking controller are designed. A Lyapunov-based stability analysis proves that the tracking and state estimation errors are UUB. ree simulation cases are used to test the influence of the designed parameters on speed estimation and tracking error. erefore, adding control parameters appropriately can reduce tracking error and energy consumption. When the designed parameters are increased to a certain critical value, increasing design parameters continuously can reduce tracking error but increase the energy consumption of the control input. Conclusions from the simulation can be used as the basis of parameter selection in international operations. In future work, multiobjective optimization for energy and tracking error will be explored.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.