Permanent-Magnet Synchronous Motor Drive System Using Backstepping Control with Three Adaptive Rules and Revised Recurring Sieved Pollaczek Polynomials Neural Network with Reformed Grey Wolf Optimization and Recouped Controller

: Owing to some nonlinear characteristics in the permanent-magnet synchronous motor (SM), such as nonlinear friction, cogging torque, wind stray torque, external load torque, and unmodeled systems, ﬁne control performances cannot be accomplished by utilizing the general linear controllers. Thereby, the backstepping approach adopting three adaptive rules and a swapping function is brought forward for controlling the rotor motion in the permanent-magnet SM drive system to reduce nonlinear uncertainties e ﬀ ects. To improve the chattering phenomenon, the backstepping control with three adaptive rules using a revised recurring sieved Pollaczek polynomials neural network (RRSPPNN) with reformed grey wolf optimization (RGWO) and a recouped controller is proposed to estimate the internal collection and external collection torque uncertainties, and to recoup the smallest fabricated error of the appraised rule. In the light of the Lyapunov stability, the on-line parametric training method of the RRSPPNN can be derived through an adaptive rule. Furthermore, to obtain a beneﬁcial learning rate and improve the convergence of the weights, the RGWO algorithm adopting two exponential-functional adjustable factors is applied to adjust the two learning rates of the weights. Then, the e ﬃ ciency of the used controller is validated by test results.


Introduction
The permanent-magnet synchronous motors (SMs) [1] with many merits are superior to the switched reluctance motors (SRMs) and induction motors (IMs). The permanent-magnet SMs [2] can offer higher efficiency, higher power density, lower power loss, and higher robustness in comparsion with the SRMs and the IMs at the same volume. The permanent-magnet SMs have mostly adopted the field-oriented control technique owing to their easy implementation. Thereby, output torque can result in lower ripple torque in comparison with the SRMs and IMs at the same output torque. On the other hand, the permanent-magnet SMs controlled by field-oriented control [1,2], which can achieve fast four-quadrant operation, are much less sensitive to the parameter variations of the motor. Therefore, they have been widely used in many industrial applications such as robotics [1], computer numerical control (CNC) tools [2], and other mechatronics [3].
The backstepping designs [4] are befitting for a large type of linearizable nonlinear systems. Each backstepping phase can produce a novel fictitious-control design denoted by previous design phases. When the procedure ends, a feedback design can achieve the primitive design aim by utilizing and have been used in certain fields, such as distribution system [26], melanoma detection [27], and feature selection in classification [28], they have poor exploration capability and suffer from local optima stagnation. So, to improve the explorative abilities of GWO, a reformed grey wolf optimization (RGWO) algorithm adopting two exponential-functional adjustable factors is put forward as the novel method in this paper. This newly proposed algorithm makes up two revisions: Firstly, it can explore new regions in the look for space because of diverse locations assigned to the leaders. This can increase the exploration and avoid local perfect stagnation problem. Secondly, an opposition-based learning method has been used in the initial half of iterations to provide diversity among the search agents. To speed up the convergence of weights in the RRSPPNN, the RGWO with two exponential-functional adjustable factors, that is the novel method, is used to adjust the two learning rates of the weights. This novel method can prevent premature convergence and to acquire optimal learning rates with fast convergence.
The better control performance of the permanent-magnet SM drive system cannot be reached by utilizing the linear controller due to the influences of these uncertainties. To heighten robustness, the backstepping approach with three adaptive rules and a swapping function is proposed to control the permanent-magnet SM drive system to trace different periodical references. With the backstepping approach with three adaptive rules and a swapping function, the rotor position of the permanent-magnet SM drive system preserves the merits of fine transient control performance and robustness to uncertainties for the tracedifferent periodical references. Moreover, to improve the large chattering influence under uncertainties, the backstepping control with three adaptive rules by utilizing RRSPPNN with RGWO is proposed to estimate the internal bunched uncertainty and external bunched force uncertainty and the recouped controller to recoup the smallest fabricated error of the appraised rule.
Furthermore, the RGWO algorithm by using two exponential-functional adjustable factors that is applied for regulating two learning rates of the weights in the RRSPPNN is a novel method to speed up the convergence of weights in this paper. Finally, the efficiency of the backstepping control with three adaptive rules using RRSPPNN with RGWO and recouped controller is validated by some test results.
The important issue in this paper is described below. Section 2 presents the models and conformation of the permanent-magnet SM drive. Section 3 describes the backstepping control with three adaptive rules using RRSPPNN with RGWO and the recouped controller. Section 4 is the examination consequences for the permanent-magnet SM utilizing three control methods at five tested events. Section 5 is the conclusion.

Models of Permanent-Magnet SM
For simplicity, in the three-phase as − bs − cs axis coordinate frames via the Clarke and Park transformations, the voltage equations in the coordinate frames transformation from the three-phase as − bs − cs axis to the qs − ds axis in the permanent-magnet SM [1] are typified by u qs = r s i qs + L qs di qs /dt + ω es (L ds i ds + λ ps ) (1) u ds = r s i ds + (L ds di ds /dt + λ ps ) − ω es L qs i qs (2) ω es = p s ω r1 (3) The electromagnetic power P es in the air gap as well as the electromagnetic torque d e1 [1] can be typified by: P es = d e1 ω es = 3p s [λ ps i qs + (L ds − L qs )i ds i qs ]ω es /2 (5) The dynamic equation can be typified by: where u ds − u qs , i ds − i qs and L ds − L qs stand for the qs − ds axis voltages, the qs − ds axis currents, and the qs − ds axis inductances; r s and p s stand for the phase winding resistance and the number of pole pairs; λ ps , ω es , and ω r1 stand for the permanent-magnet flux linkage, the electric angular speed, and the angular speed of the rotor; θ es and θ r1 stand for the electrical angular position and the mechanical angular position of the rotor; d e1 , d l1 , d f 1 , d c1 , and d w1 stand for the electromagnetic torque, the external load torque, the flux saturation torque, the cogging torque, and the wind stray torque; H s and U s stand for total moment of inertia of permanent-magnet SM and the total viscous frictional coefficient of permanent-magnet SM.

Conformation of Permanent-Magnet SM Drive
The decoupled control technology of the permanent-magnet SM drive system is general adopting the field-oriented control (FOC) [1,2] to raise dynamic performance. The electromagnetic torque is produced by the qs− axis current based on the FOC and the rotor flux is generated by the ds− axis current only. When i ds is equal to zero and λ ps is equal to a constant in Equations (5) and (6), then the electromagnetic torque d e1 will be direct ratio to i qs for a permanent-magnet SM drive in the closed-loop control. When the generated torque is linearly direct ratio to the qs− axis current, and the ds− axis rotor flux is a fixed value, the larger torque per ampere can be reached. The electromagnetic torque Equation (6) can be typified by where k s = 3p s λ ps /2 stands for the propulsion coefficient and i * qs stands for the mandate of control current. The permanent-magnet SM drive system can be predigested as W t (s) = 1/(sH s + U s ). Figure 1 is the conformation of FOC permanent-magnet SM drive system, which makes up an encoder and three Hall sensors, a permanent-magnet SM, a sinusoidal pulse-width-modulation (SPWM) current control, a coordinate transformation including inverse Park and Clarke coordinate's transformations, cos/sin generator, a speed control loop and a position control loop. The control technologies in the real-time realization are realized by utilizing the digital signal processor (DSP) controller. Rotor of permanent-magnet SM is equipped on magnet force brake that is mounted with different sizes of iron disks to change the total moment of inertia and the total viscous frictional coefficient, and to add load torque.
The FOC was realized by a digital signal processor (DSP) controller. The nominal values of used permanent-magnet SM are given as 3-phase, 2-poles, 60 Hz, 220 V, 1 kW, 2.8 A, and 3600 r/min. For the convenience of controller design, the position and speed signals in the control loop are set at 1 V = 50 rad and 1 V = 50 rad/s, respectively. The mechanical and electrical parameters of the permanent-magnet SM are given as H s = 2.142 × 10 −3 Nm sec 2 = 0.1071 Nm sec rad/V, U s = 5.86 × 10 −3 Nms/rad = 0.293 Nm/V, r s = 2.5 Ω, L ds = L qs = 4.62 mH, k s = 0.947 Nm/A by means of an open circuit test, short circuit test, blocked rotor test, and added load test. With the fulfillment of FOC [1][2][3], the permanent-magnet SM drive system can be predigested as the control block diagram portrayed in Figure 2. The perfect electromagnetic property for the drive system is hence implemented by controlling the torque current distributions to lie in the qs− axis current when the d-axis current is equal to zero. Then, the torque per amp property for the drive system will generate. The FOC was realized by a digital signal processor (DSP) controller. The nominal values of used permanent-magnet SM are given as 3-phase, 2-poles, 60 Hz, 220 V, 1 kW, 2.8 A, and 3600 r/min. For the convenience of controller design, the position and speed signals in the control loop are set at 1 V = 50 rad and 1 V = 50 rad/s, respectively. The mechanical and electrical parameters of the permanent-  [1][2][3], the permanent-magnet SM drive system can be predigested as the control block diagram portrayed in Figure 2. The perfect electromagnetic property for the drive system is hence implemented by controlling the torque current distributions to lie in the − qs axis current when the d-axis current is equal to zero. Then, the torque per amp property for the drive system will generate.

Design of the Controller
When the permanent-magnet SM drive system with the electromagnetic torque, the wind stray torque, the flux saturation torque, the cogging torque, the parametric variations, and the external load torque disturbance is enacted, then Equation (7) is typified by

Design of the Controller
When the permanent-magnet SM drive system with the electromagnetic torque, the wind stray torque, the flux saturation torque, the cogging torque, the parametric variations, and the external load torque disturbance is enacted, then Equation (7) is typified by where θ r1 and v s stand for the rotor position and rotor speed of the SM to be presumed bounded; ∆e s and ∆ f s stand for two parametric uncertainties from H s and U s to be presumed bounded; e s = −U s /H s ; f s = k s /H s > 0; g s = −1/H s stand for three real numbers to be presumed bounded; l s = i qs is the control propulsion of the permanent-magnet SM drive system, i.e., the propulsion current. Equation (9) can be typified by where w 1 = ∆e s v s and w 2 = ∆ f s d s stand for two parametric variation that are to be presumed bounded; w 3 = g s (d f 1 + d c1 ) and w 4 = g s (d l1 + d w1 ) stand for the internal bunched uncertainty and external bunched force uncertainty to be presumed bounded.
The trace reference locus b m = θ m1 is the control goal. The design procedure is as below.
The trace error is typified by Take differential of (13) by The stabilizing function is typified by: where k m1 and k m2 stand for two positive constants; µ m = e m1 (t)dt stands for the integral factor. The fictitious trace error is typified by: Take the differential of (16) by where w 1 , w 2 , and w 3 stand for three unknown parameters. The estimated errors are typified by where e w1 , e w2 , and e w3 are the estimated errors;ŵ 1 ,ŵ 2 , andŵ 3 are the estimated values of w 1 , w 2 and w 3 . The internal bunched uncertainty and external bunched force uncertainty w 4 satisfies the condition |w 4 | ≤ β 1 and is to be presumed bounded.

Design of the Backstepping Control with Three Adaptive Rules and a Swapping Function
So, to design the backstepping control with three adaptive rules and a swapping function, the Lyapunov function can be typified by: By taking the differential of y m1 and by utilizing Equations (14)- (20) and the integral factor µ m = e m1 (t)dt, Equation (21) can be typified by dy m1 /dt = e m1 de m1 /dt + e m2 de m2 /dt + e w1 ς −1 In accordance with Equation (22), the control propulsion l s in the backstepping control with three adaptive rules and a swapping function can be typified by where k m3 stands for a positive constant; β 1 stands for upper bound that is a constant, and f −1 s β 1 sgn(e m2 ) stands for a swapping function. By utilizing Equation (23), Equation (22) can be typified by: For reaching dy m1 (t)/dt ≤ 0, three adaptive rules dŵ 1 /dt, dŵ 2 /dt, and dŵ 3 /dt can be typified by By utilizing Equations (25)- (27), and |w 4 | ≤ β 1 , then Equation (24) can be typified by Equation (28) shows dy m1 (t)/dt to be negative semi-definite (i.e., y m1 (t) ≤ y m1 (0)), meaning that e m1 and e m2 are bounded. The following term is typified by The integration of Equation (29) is typified by: Because q m1 (e m1 (0), e m2 (0)) is bounded and q m2 (e m1 (t), e m2 (t)) is nonincreasing and presumed bounded, then lim t→∞ t 0 q m1 (v)dτ < ∞. Moreover, dh m1 (t)/dt is presumed bounded, hence h m1 (t) is a uniformly continuous function. By utilizing the Barbalat's lemma [29,30], it can be portrayed that lim t→∞ h m1 (t) = 0. That is, e m1 and e m2 will converge to zero when t → ∞ . Furthermore, lim t→∞ θ r1 (t) = θ m1 and lim t→∞ v s = dθ m1 /dt. The stability of the backstepping control with three adaptive rules and a swapping function can be guaranteed, and consequently, the control block diagram is portrayed in Figure 3.

Design of the Backstepping Control with Three Adaptive Rules Using RRSPPNN with RGWO and Recouped Controller
Because the internal bunched uncertainty and external bunched force uncertainty 4 w is unknown, and its upper bound is troublesome to be decided. The appraised value 4 w of the internal bunched uncertainty and external bunched force uncertainty 4 w is not easy to be estimated, and consequently, the revised recurring sieved Pollaczek polynomials neural network (RRSPPNN) is proposed to adapt the real value of the internal bunched uncertainty and external bunched force uncertainty 4 w .

Constitution of RRSPPNN
The RRSPPNN has a three-layer constitution, with the first layer (input layer), the second layer (hidden layer 1), and the third layer (output layer) portrayed in Figure 4. The semaphore intentions in each node for each layer are explained in the following expression. At the first layer, input semaphore and output semaphore are typified by

Design of the Backstepping Control with Three Adaptive Rules Using RRSPPNN with RGWO and Recouped Controller
Because the internal bunched uncertainty and external bunched force uncertainty w 4 is unknown, and its upper bound is troublesome to be decided. The appraised valueŵ 4 of the internal bunched uncertainty and external bunched force uncertainty w 4 is not easy to be estimated, and consequently, the revised recurring sieved Pollaczek polynomials neural network (RRSPPNN) is proposed to adapt the real value of the internal bunched uncertainty and external bunched force uncertainty w 4 .

Constitution of RRSPPNN
The RRSPPNN has a three-layer constitution, with the first layer (input layer), the second layer (hidden layer 1), and the third layer (output layer) portrayed in Figure 4. The semaphore intentions in each node for each layer are explained in the following expression. In accordance with Equation (37), the control propulsion s s l l= in the backstepping control with three adaptive rules by using RRSPPNN with RGWO and recouped controller can be typified by By utilizing Equation (38), then Equation (37) can be typified by By utilizing Equations (25)  to reduce uncertainties influences can be typified by:  At the first layer, input semaphore and output semaphore are typified by where q 1 1 = θ m1 − θ r1 = e m1 and q 1 2 = e m1 (1 − z −1 ) = ∆e m1 stand for the speed discrepancy and speed discrepancy alteration, respectively. B is the iteration count. b 1 ac stands for the recurring weight through the third layer and the first layer. o 3 c stands for the output of node at the third layer. The symbol stands for a multiply factor.
At the second layer, input semaphore and output semaphore are typified by where δ stands for the recurring gain at the second layer. Sieved-Pollaczek polynomials function [21,31] is adopted as the activation function. SP λ j (q; t) stands for the sieved Pollaczek polynomials in the interval [−1, 1]. SP λ 0 (q; t) = 1, SP λ 1 (q; t) = 2q and SP λ 2 (q; t) = 4q 2 − 1 stand for the zero-order, first order and second order sieved Pollaczek polynomials, respectively. The sieved Pollaczek polynomials may be generated by the recurrence relation [21,31] The symbol stands for a summation factor.
At the third layer, semaphore and output semaphore are typified by where b 2 cb stands for the connecting weight through the second layer and the third layer. g 3 c stands for the linear activation function. The output o 3 c (B) at the third layer of the RRSPPNN can be typified by: where stands for the weight vector at the third layer and the input vector at the third layer, respectively. The smallest fabricated error ϕ w4 is typified by where M * stands for an ideal weight vector that reaches the smallest fabricated error. So as to make up the smallest fabricated error ϕ w , the recouped controller x c with an appraised rule is proposed. It is presumed that the small positive number e w4 stands for greater than absolute value of ϕ w , i.e., e w4 ≥ ϕ w . The Lyapunov function is typified by where σ 1 stands for an adaptive gain. e w4 =ê w4 − e w4 stands for the appraised error to be presumed bounded. By taking the derivative of y m2 utilizing Equations (14)- (20) and the integral factor µ m = e m1 (t)dt, then Equation (36) is typified by In accordance with Equation (37), the control propulsion l s =l s in the backstepping control with three adaptive rules by using RRSPPNN with RGWO and recouped controller can be typified by By utilizing Equation (38), then Equation (37) can be typified by By utilizing Equations (25)- (27) and e w4 =ê w4 − e w4 , then Equation (39) can be typified by

Recouped Controller with an Adaptive Rule
To reach dy m2 /dt ≤ 0, the adaptive rule d M/dt, the recouped controller x c , and the appraised rule dê w4 /dt to reduce uncertainties influences can be typified by: x c =ê w4 sgn(e m2 ) By substituting Equations (41)-(43) into Equation (40) and by utilizing e w4 ≥ ϕ w , then Equation (40) can be typified by: Equation (44) portrays dy m2 (t)/dt to be negative semi-definite, i.e., y m2 (t) ≤ y m2 (0), meaning that e m1 and e m2 are bounded. By utilizing the Barbalat's lemma [29,30], it can be represented that y m2 (t) → 0 at t → ∞ by way of Equations (29), (30) and (44), i.e., e m1 and e m2 will converge to zero at t → ∞. The stability of the backstepping control with three adaptive rules by using RRSPPNN with RGWO and recouped controller can be ensured and consequently the control block diagram is portrayed in Figure 5.

Training of the RRSPPNN
By utilizing the Lyapunov stability and the gradient descent skill with the chain rule, a training skillfulness of parameters in the RRSPPNN can be derived. The RGWO with two adjusted factors is applied to look for two better learning rates in the RRSPPNN to acquire faster convergence. The connecting weight parametric presented in Equation (41) can be typified by: A goal function that explains the online training procedure of the RRSPPNN is typified by

Training of the RRSPPNN
By utilizing the Lyapunov stability and the gradient descent skill with the chain rule, a training skillfulness of parameters in the RRSPPNN can be derived. The RGWO with two adjusted factors is applied to look for two better learning rates in the RRSPPNN to acquire faster convergence. The connecting weight parametric presented in Equation (41) can be typified by: A goal function that explains the online training procedure of the RRSPPNN is typified by The adaptive learning rule of the connecting weight is typified by: It is well-known that ∂L 2 /∂o 3 c = −e m2 by way of Equations (45) and (47). Hence, the adaptive learning rule of recurring weight b 1 ac is typified by where γ 2 stands for the learning rate. To acquire better convergence, the RGWO is applied to look for two changeable learning rates in the RRSPPNN. Additionally, for improving convergence and looking for two perfect learning rates, the RGWO with two adjusted factors is proposed in this study.

Algorithm of Reformed Grey Wolf Optimization (RGWO)
In the RGWO, the optimization is conducted by α, β, and ρ. The RGWO algorithm can be typified by: where H(a 1 + 1) = [γ 1 γ 2 ] is a vector two learning rates, H 1 (a 1 ), H 2 (a 1 ), H 3 (a 1 ) are typified by: where α(a 1 ), β(a 1 ), ρ(a 1 ) stand for the three vectors as the three best solutions. R 1 (a 1 ), R 2 (a 1 ), R 3 (a 1 ) and M 1 (a 1 ), M 2 (a 1 ), M 3 (a 1 ) can be typified by: where ϕ 1 and ϕ 2 stand for two random vectors. The updated numbers of two adjusted factors c 1 (a 1 ) and d 1 (a 1 ) control the tradeoff between exploration and exploitation. Two exponential-functional adjustable factors c 1 (a 1 ) and d 1 (a 1 ) stand for the updated values at iteration according to the following presentation by: where a 1 stands for the iteration number; a t1 and a t2 stand for the total numbers of iteration allowed for the optimization. At last, H(a 1 + 1) = [γ 1 γ 2 ] stands for the best solution in connection with the learning rates γ i (t), i = 1, 2 of the two weights in the RRSPPNN. Hence, the better numbers could be optimized by utilizing RGWO with two adjusted factors that yield two changeable learning rates for two weights to look for two perfect values and to speed-up the convergence of two weights.

Test Results
A block diagram of the FOC permanent-magnet SM drive system utilizing the DSP controller is portrayed in Figure 1. A photo of the examination structure is portrayed in Figure 6. The sampling time of the control program in the examination is set as 2 ms.

Test Results
A block diagram of the FOC permanent-magnet SM drive system utilizing the DSP controller is portrayed in Figure 1. A photo of the examination structure is portrayed in Figure 6. The sampling time of the control program in the examination is set as 2 ms. A DSP controller involves 18 channels of input/output (I/O) ports with 6 channels of pulsewidth-modulation (PWM) ports, 6 channels of analog-digital (A/D) converters, and 2 channel encoder connective ports. The coordinate transformation in the field-oriented control (FOC) is realized by DSP controller. The used control technologies in the real-time realization by utilizing the DSP controller are composed of the core program and the sub-core interrupt service routine (SCISR) in the DSP controller as portrayed in Figure 7. In the core program, parameters and I/O initialization are processed. The interrupt time for the SCISR is set. After permitting the interruption, the core program is used to monitor control data. The SCISR with 2 msec sampling time is used for reading the rotor position of the permanent-magnet SM drive system from encoder and three-phase currents by way of A/D converter circuit, calculating reference model and position error, executing lookup table and coordinate transformation, executing the backstepping control with three adaptive rules using RRSPPNN with RGWO and recouped controller, and outputting three-phase current mandates to swap sinusoidal pulse-width-modulation (SPWM) voltage source inverter with three-sets of insulated-gate bipolar transistor (IGBT) power modules by way of the lockout-time and isolated circuits. The SPWM voltage source inverter with three-sets of IGBT power modules is carried out by SPWM control with a switching frequency of 15 kHz. Additionally, the tested bandwidth of the position control loop and the tested bandwidth of the current control loop are about 90 and 900 Hz for the permanent-magnet SM drive system under the nominal event. The proposed controllers are realized by the DSP controller. The coordinate transformation in the FOC is realized by the DSP controller. The control goal is to control the rotor to rotate 6.28 rad cyclically. Then, when the mandate is a sinusoidal reference locus, the reference is set one. A DSP controller involves 18 channels of input/output (I/O) ports with 6 channels of pulse-width-modulation (PWM) ports, 6 channels of analog-digital (A/D) converters, and 2 channel encoder connective ports. The coordinate transformation in the field-oriented control (FOC) is realized by DSP controller. The used control technologies in the real-time realization by utilizing the DSP controller are composed of the core program and the sub-core interrupt service routine (SCISR) in the DSP controller as portrayed in Figure 7. In the core program, parameters and I/O initialization are processed. The interrupt time for the SCISR is set. After permitting the interruption, the core program is used to monitor control data. The SCISR with 2 msec sampling time is used for reading the rotor position of the permanent-magnet SM drive system from encoder and three-phase currents by way of A/D converter circuit, calculating reference model and position error, executing lookup table and coordinate transformation, executing the backstepping control with three adaptive rules using RRSPPNN with RGWO and recouped controller, and outputting three-phase current mandates to swap sinusoidal pulse-width-modulation (SPWM) voltage source inverter with three-sets of insulated-gate bipolar transistor (IGBT) power modules by way of the lockout-time and isolated circuits. The SPWM voltage source inverter with three-sets of IGBT power modules is carried out by SPWM control with a switching frequency of 15 kHz. Additionally, the tested bandwidth of the position control loop and the tested bandwidth of the current control loop are about 90 and 900 Hz for the permanent-magnet SM drive system under the nominal event. The proposed controllers are realized by the DSP controller. The coordinate transformation in the FOC is realized by the DSP controller. The control goal is to control the rotor to rotate 6.28 rad cyclically. Then, when the mandate is a sinusoidal reference locus, the reference is set one. For a comparison of control performance with the four controllers, five events are provided in the experiment. The four controllers are the popular PI controller as the controller FC1, the backstepping control with three adaptive rules and a swapping function as the controller FC2, the modified recurrent Rogers-Szego polynomials neural network with DGWO [10] as the controller FC3, and the backstepping control with three adaptive rules using RRSPPNN with RGWO and recouped controller as the controller FC4. Five tested events are as follows. Event CQ1 is the nominal event at periodic step command from 0 rad to 6.28 rad. Event CQ2 is the cogging torque, the column friction torque, and the Stribeck effect torque, and the parameters variations event with 4 times the nominal value at periodic step command from 0 rad to 6.28 rad. Event CQ3 is the nominal event due to periodic sinusoidal command from −6.28 rad to 6.28 rad. Event CQ4 is the cogging torque, the column friction torque and the Stribeck effect torque and the parameters variations event with 4 times the For a comparison of control performance with the four controllers, five events are provided in the experiment. The four controllers are the popular PI controller as the controller FC1, the backstepping control with three adaptive rules and a swapping function as the controller FC2, the modified recurrent Rogers-Szego polynomials neural network with DGWO [10] as the controller FC3, and the backstepping control with three adaptive rules using RRSPPNN with RGWO and recouped controller as the controller FC4. Five tested events are as follows. Event CQ1 is the nominal event at periodic step command from 0 rad to 6.28 rad. Event CQ2 is the cogging torque, the column friction torque, and the Stribeck effect torque, and the parameters variations event with 4 times the nominal value at periodic step command from 0 rad to 6.28 rad. Event CQ3 is the nominal event due to periodic sinusoidal command from −6.28 rad to 6.28 rad. Event CQ4 is the cogging torque, the column friction torque and the Stribeck effect torque and the parameters variations event with 4 times the nominal value due to periodic sinusoidal command from −6.28 rad to 6.28 rad. Event CQ5 is the adding load torque disturbance τ lr = 2.5 Nm.
Two control gains of the popular PI controller as the controller FC1 are k pp = 4.5, and k ip = k pp /T ip = 1.8 by using the Kronecker method to construct a stability boundary in the k pp and k ip plane [32][33][34] on the tuning of the PI controller at Event CQ1 in the position trace so as to reach fine steady-state and transient-state control responses.
Some control gains of the modified recurrent Rogers-Szego polynomials neural network with DGWO [10] as the controller FC3 are given as c x1 = 2.2, c x2 = 2.6, c x3 = 2.1, η 1 = 0.1, η 2 = 0.1, η 3 = 0.1, γ = 0.2, τ = 0.1 according to heuristic knowledge [4] at Event CQ1 in the position trace so as to reach fine steady-state and transient-state control responses. Moreover, numbers of neurons in the input layer, the hidden layer, and the output layer of the modified recurrent Rogers-Szego polynomials neural network are 2 neurons, 4 neurons and 1 neuron, respectively, so as to demonstrate the effectiveness of the controller adopting small neuron numbers. The method proposed by Lewis et al. [35] is used to initialize some parameters of the modified recurrent Rogers-Szego polynomials neural network. The adjustment process of these parameters involve a continuous reaction for the duration of the examination.
Some control gains of the backstepping control with three adaptive rules using RRSPPNN with RGWO and recouped controller as the controller FC4 are k m1 = 2.2, k m2 = 2.6, k m3 = 2.1, ς 1 = 0.1, ς 2 = 0.1, ς 3 = 0.1, σ 1 = 0.2, δ = 0.1 according to heuristic knowledge [4] at Event CQ1 in the position trace so as to reach fine steady-state and transient-state control responses. Moreover, the number of neurons in the first layer, the second layer, and the third layer of the RRSPPNN are 2 neurons, 4 neurons, and 1 neuron, respectively, so as to demonstrate the effectiveness of the controller by adopting small neuron numbers. The method proposed by Lewis et al. [35] is used to initialize some parameters of the RRSPPNN. The adjustment process of these parameters is keeping continusly reaction in the duration of the experiments.
All of the experiments obtained by utilizing the four controllers for controlling the permanent-magnet SM drive system at five events are as follows.                        Figures 8a and 14a obtained by utilizing the controller FC1 for controlling the permanent-magnet SM drive system at Event CQ1 and Event CQ3 are displayed as fine trace responses of the rotor positions. Figures 11a and 17a obtained by utilizing the controller FC1 for controlling the permanent-magnet SM drive system at Event CQ2 and Event CQ4 are displayed in sluggish trace responses of the rotor position owing to bigger nonlinear disturbance. Because of inappropriate tuning gains or the degenerate nonlinear effect, the linear controller has weak robustness under bigger nonlinear disturbance. Figures 8b, 11b, 14b and 17b obtained by utilizing the controller FC2 for controlling the permanent-magnet SM drive system at Events CQ1, CQ2, CQ3, and CQ4 are displayed as good trace responses of the rotor positions. However, Figures 10b, 13b, 16b and 19b are displayed as serious vibration in the control propulsions by utilizing the swapping function with large upper bound at Events CQ1, CQ2, CQ3, and CQ4. It is a well-known fact that the control propulsions with serious vibration will wear the bearing mechanism and might excite unstable system dynamics. Figures 8c, 11c, 14c and 17c obtained by utilizing the controller FC3 for controlling the permanent-magnet SM drive system at Events CQ1, CQ2, CQ3, and CQ4 are displayed as better trace responses of the rotor positions due to adaptive mechanism action. Figures 10c, 13c, 16c and 19c displayed a small vibration in the control propulsions at Events CQ1, CQ2, CQ3, and CQ4. Due to the on-line adaptive adjustment of the modified recurrent Rogers-Szego polynomials neural network [10], the magnitudes of vibration in the control propulsions at Events CQ1, CQ2, CQ3, and CQ4 displayed in Figures 10c, 13c, 16c and 19c have been slightly improved. Figures 8d, 11d, 14d and 17d obtained by utilizing the controller FC4 for controlling the permanent-magnet SM drive system at Events CQ1, CQ2, CQ3, and CQ4 are displayed as best trace responses of the rotor positions due to on-line adaptive mechanism action. Figures 10d, 13d, 16d and 19d are displayed as smaller vibrations in the control propulsions at Events CQ1, CQ2, CQ3 and CQ4 due to on-line adaptive mechanism action of the RRSPPNN. Due to on-line adaptive adjustment of the RRSPPNN under bigger nonlinear disturbance the magnitudes of vibration in the control propulsions at Events CQ1, CQ2, CQ3, and CQ4 displayed in Figures 10d, 13d, 16d and 19d have been obviously improved. Figure 20d obtained by utilizing the controller FC4 for controlling the permanent-magnet SM drive system at Event CQ5 under load regulation is better than the controller FC1, FC2, and FC3 displayed in Figure 20a-c.

Discussion and Analysis
Additionally, the control performances displayed in comparsion results by using the controllers FC1, FC2, FC3, and FC4 are listed in Table 1 in connection with five events with some test results. The 0.21, 0.19, 0.15, and 0.10 are the maximum errors of e m1 (rad) by utilizing the controllers FC1, FC2, FC3, and FC4 at Event CQ1, respectively. The 0.11, 0.09, 0.07, and 0.05 are the root-mean-square errors of e m1 (rad) by utilizing the controllers FC1, FC2, FC3, and FC4 at Event CQ1, respectively. The 0.56, 0.36, 0.28, and 0.19 are the maximum errors of e m1 (rad) by utilizing the controllers FC1, FC2, FC3, and FC4 at Event CQ2, respectively. The 0.27, 0.18, 0.13, and 0.09 are the root-mean-square errors of e m1 (rad) by utilizing the controllers FC1, FC2, FC3, and FC4 at Event CQ2, respectively. The 0.21, 0.18, 0.14, and 0.10 are the maximum errors of e m1 (rad) by utilizing the controllers FC1, FC2, FC3, and FC4 at Event CQ3, respectively. The 0.10, 0.09, 0.07, and 0.05 are the root-mean-square errors of errors of e m1 (rad) by utilizing the controllers FC1, FC2, FC3, and FC4 at Event CQ3, respectively. The 0.52, 0.37, 0.27, and 0.18 are the maximum errors of e m1 (rad) by utilizing the controllers FC1, FC2, FC3, and FC4 at Event CQ4, respectively. The 0.25, 0.18, 0.13, and 0.09 are the root-mean-square errors of e m1 (rad) by utilizing the controllers FC1, FC2, FC3, and FC4 at Event CQ4, respectively. The 3.01, 1.57, 1.13, and 0.75 are the maximum errors of e m1 (rad) by utilizing the controllers FC1, FC2, FC3, and FC4 at Event CQ5, respectively. The 1.50, 0.78, 0.56, and 0.32 are the root-mean-square errors of e m1 (rad) by utilizing the controllers FC1, FC2, FC3, and FC4 at Event CQ5, respectively. The controller FC4 has smaller trace error in comparison with the controllers FC1, FC2, and FC3. The controllers FC4 indeed yields the exellent control performance from Table 1.  Table 2 for test results. In Table 2, various performances with regard to the control propulsion with vibration, the dynamic response, the ability of load regulation, the convergence speed, the position trace error, and the rejection ability of parameter disturbance in the controllers FC4 are superior to the controllers FC1, FC2, and FC3. Finally, the robust control performance of the controller FC4 demonstrates outstanding performance for controlling the permanent-magnet SM drive system in the trace of the periodic step and sinusoidal commands under the occurrence of parameter disturbance and load regulation due in large part to the on-line adaptive adjustment of the RRSPPNN. Good (maximum error as 1.57 (rad)with adding load at 6.28 (rad)) Better (maximum error as 1.13 (rad) with adding load at 6.28 (rad)) Best (maximum error as 0.09 (rad) with adding load at 6.28 (rad))

Conclusions
The backstepping control with three adaptive rules and RRSPPNN with RGWO and recouped controller is used to determine the best values for parameters in neural network learning rules, and thereby robustness in learning control can be improved.
The main contribution of this study is as follows. Firstly, the field-oriented control has been smoothly applied to control the permanent-Magnet SM drive system to speed up the control response. Moreover, the backstepping controller with three adaptive rules and a swapping function has been smoothly derived to overcome influences under the external lumped force uncertainty disturbances in light of the Lyapunov function. Further, the backstepping control with three adaptive rules and RRSPPNN with RGWO and recouped controller to estimate the external lumped force uncertainty has been smoothly derived in the light of the Lyapunov function for diminishing the external lumped force uncertainty effect and improving the chattering phenomenon. The error recouped controller to recoup the smallest fabricated error of the error estimation law has been smoothly derived in light of the Lyapunov function. Two optimal learning rates of the RRSPPNN with two exponential-functional adjustable factors have been smoothly calculated by utilizing the RGWO algorithm to speed up the parameter's convergence.
Finally, some control performances regarding the chattering of control propulsion, the position response, the ability of load force adjustment, the position tracing error, and the refusal ability of unknown parameters disturbance by using the backstepping control with three adaptive rules and RRSPPNN with RGWO and recouped controller are more exceptional than the popular PI controller, the backstepping control with three adaptive rules and a swapping function, and the modified recurrent Rogers-Szego polynomials neural network with DGWO [10] from Tables 1 and 2.

Conflicts of Interest:
The author declares no conflict of interest.