Permanent-Magnet SLM Drive System Using AMRRSPNNB Control System with DGWO

: Because permanent-magnet synchronous linear motors (SLM) still exhibit nonlinear friction, ending e ﬀ ects and time-varying dynamic uncertainties, better control performances cannot be achieved by using common linear controllers. We propose a backstepping approach with three adaptive laws and a beating function to control the motion of permanent-magnet SLM drive systems that enhance the robustness of the system. In order to reduce greater vibration in situations with uncertainty actions in the aforementioned control systems, we propose an adaptive modiﬁed recurrent Rogers–Szego polynomials neural network backstepping (AMRRSPNNB) control system with three adaptive laws and reimbursed controller with decorated gray wolf optimization (DGWO), in order to handle external bunched force uncertainty, including nonlinear friction, ending e ﬀ ects and time-varying dynamic uncertainties, as well as to reimburse the minimal rebuild error of the reckoned law. In accordance with the Lyapunov stability, online parameter training method of the modiﬁed recurrent Rogers–Szego polynomials neural network (MRRSPNN) can be derived by utilizing an adaptive law. Furthermore, to help reduce error and better obtain learning fulﬁllment, the DGWO algorithm was used to change the two learning rates in the weights of the MRRSPNN. Finally, the usefulness of the proposed control system is validated by tested results.


Introduction
Linear motors can provide linear motion in the absence of motion translators, such as ball screws, gears and belts that reduce motor limitations and complexity. Most linear motors have lower load capacity than other types of linear actuators. The permanent-magnet synchronous linear motor (SLM), with a direct drive design, has some good performance characteristics, such as higher speed, higher precision, less friction, no backlash, maintenance-free operation and extra-high thrust force [1]. Hence, the permanent-magnet SLM has been broadly applied in machine tools, semiconductor manufacturing systems and industrial robots [1][2][3]. The permanent-magnet SLM has larger thrust density that is superior to the other linear motors.
The backstepping control kill [4][5][6] has been applied in various kinds of linear feedback systems. First, some proper functions with the state variables are chosen as pseudo-control inputs in the subsystems of the entire system. Then, each backstepping level results in a new pseudo-control variable that is denoted by a term of the pseudo-control design from the preceding design levels. The primary design goal of a feedback design is the true control of input. When the last Lyapunov function from The model error and system disturbance can seriously affect various control performance under the influences of unknown parameters and dynamics system interfaces. In light of the influence of these uncertainties, it is not easy to reach better control fulfillment utilized by the linear controller for the permanent-magnet SLM drive system. Hence, as to enhance robustness the backstepping accuracy with three adaptive laws and a beating function is proposed to govern the motion of permanent-magnet SLM drive system for the tracking of periodic reference trajectories. The backstepping accuracy with three adaptive laws and a beating function, the permanent-magnet SLM drive system holds the merits of good transient control fulfillment and robustness to uncertainties. In addition, so as to reduce great vibration situation with uncertainties action for the before-mentioned control system, we propose the adaptive modified recurrent Rogers-Szego polynomials neural network backstepping (AMRRSPNNB) control system with reimbursed controller and DGWO with two adjusted factors to reckon the external bunched force uncertainty and to reimburse the minimal rebuilt error of the reckoned law. Further, the DGWO algorithm-by using two adjusted factors with cosine function-is an innovative method to quicken convergent speed in this study. Regulating the two variable learning rates of the weights in the MRRSPNN is the use of the DGWO algorithm to speed-up the parameter's delay. Finally, the usefulness and robustness of the proposed AMRRSPNNB control system with reimbursed controller and DGWO with two adjusted factors were validated by some tested results.
The main structure of this study is as follows: Section 2 presents the configuration of permanent-magnet SLM drive. Section 3 presents the AMRRSPNNB control system with reimbursed controller and DGWO with two adjusted factors. Section 4 presents the tested results for the permanent-magnet SLM by using three control methods at five cases. Section 5 presents the conclusions.

Models of Permanent-Magnet SLM
The electric models in the synchronous rotating reference frame for the permanent-magnet SLM with d-q axis model [1][2][3] are denoted by: u qs = r s i qs + L qs . i qs + ω es (L ds i ds + λ ps ) (1) u ds = r s i ds + (L ds . i ds + . λ ps ) − ω es L qs i qs (2) ω es = p s ω rs = p s π v rs /τ (3) v es = p s v rs = 2τ f s (4) where u ds , u qs denote the d and q axis voltages, i ds , i qs denote the d and q axis currents, r s denotes the phase winding resistance, L ds , L qs denote the d and q axis inductances, λ ps denotes the permanent magnet flux linkage, ω es and ω rs denote the electric angular velocity and the mover angular velocity, v rs and v es denote the linear velocity and the electric linear velocity,p s and τ denote the number of pole pairs and the pole pitch, f s denotes the electric frequency. The electromagnetic power P es and the electromagnetic force F es [1][2][3] are denoted by P es = F es v es = 3p s [λ ps i qs + (L ds − L qs )i ds i qs ]ω es /2) F es = 3π p s [λ ps i qs + (L ds − L qs )i ds i qs ]/2τ (6) Thus, the mover dynamic equation can be denoted by W s . v rs + Q s v rs = F es − F ls − F ws − F f s − F cs (7) Energies 2020, 13 where W s and Q s are the entire mass with the moving component system and the viscous friction with iron contact part, F es , F ls , F ws , F f s and F cs are the electromagnetic force, the external load force, the wind resistance force, the flux saturation force and the cogging force, respectively.

Formation of Permanent-Magnet SLM Drive
The common control skill of a permanent-magnet SLM drive system is to adopt the FOC [1][2][3]. Hall sensors can determine if the flux is in place with the d-q axis. The rotor flux is produced in the d-axis only when the current vector is produced in the q-axis for the FOC. When i ds is equal to zero and λ ps is fixed value in Equations (5) and (6), then the electromagnetic force F es will be proportional to i qs in the closed-loop control for a permanent-magnet SLM drive. Since the produced force is linearly proportional to the q-axis current under fixed the d-axis rotor flux, the maximum force per ampere can be reached. The electromagnetic force equation in Equation (6) can be represented by where k s = 3π p s λ ps /2τ is the propulsion coefficient; i * qs is the command of propulsion current. The permanent-magnet SLM drive system can be understood to control the system with H t (s) = 1/(sW s + Q s ) that s is the Laplace's operation. The formation of a field-oriented permanent-magnet SLM drive system is displayed in Figure 1, which makes up a permanent-magnet SLM, a speed control loop, a position control loop, a triangular wave comparison current control, a coordinate transformation, cos θ es / sin θ es generator, a linear scale and three Hall sensors. We used the three Hall sensors to find the flux position of the permanent magnet. The mover of permanent-magnet SLM is equipped on various sizes of iron disks to change the entire mass of the moving component system and the variously viscous friction force.   With the realization of FOC, the permanent-magnet SLM drive system can be understood as the block diagram expressed in Figure 2. The better electromagnetic performance is thus enforced by governing the primary current distributions to lie in the q-axis, i.e., i ds = 0. A linear force per amp characteristic can be reached for the actuator.

Design of the AMRRSPNNB Control System with Reimbursed Controller and DGWO
When the actual permanent-magnet SLM drive system operates with parameter variations, external load disturbance, friction force, wind resistance force, flux saturation force and cogging force is acted, then Equation (7) can be rewritten by The FOC was enforced by the TMS320C32 DSP control system. A main personal computer downloads the program running on the DSP. The specifications of permanent-magnet SLM are 220 V, 1.75 A, 0.5 kW, 112 N in this study. The position control and speed control signals are set at 1 V = 2 mm and 1 V = 2 mm/sec for the convenience of the controller design. The mechanical parameters of the system are as follows W s = 2.4 kg = 0.048 Nsec/V, Q s = 89.54 kg/ sec = 0.179 N/V, k s = 0.104 N/A.

Design of the AMRRSPNNB Control System with Reimbursed Controller and DGWO
When the actual permanent-magnet SLM drive system operates with parameter variations, external load disturbance, friction force, wind resistance force, flux saturation force and cogging force is acted, then Equation (7) can be rewritten by . z s = (c s + ∆c s )z s + (d s + ∆d s )u s + h s (F ls + F ws + F f s + F cs ) (9) . l rs = v rs = z s (10) y s = l rs (11) where l rs and z x denote the mover position and mover velocity of the SLM to be uniformly continuous and bounded; ∆c s and ∆d s denote two parameters uncertainties from W s and Q s to be uniformly continuous and bounded; c s = −Q s /W s ; d s = k s /W s > 0; h s = −1/W s are three known constants to be uniformly continuous and bounded; u s = i qs is the control propulsion to the permanent-magnet SLM drive system, i.e., the propulsion current. Reformulate (10), then where f s1 = ∆c s z s and f s2 = ∆d s u s are two parameters variation to be bounded; f s3 = h s (F f s + F cs ) and f s4 = h s (F ls + F ws ) are the internal bunched uncertainty and external bunched force uncertainty to be bounded. Tracking the reference trajectory y m = l rm asymptotically is the control objective. Hence, as to research good position-tracking performance, the backstepping control system with three adaptive laws and a beating function is devised as follows.
Define the tracking error by e g1 = y m − y s = l rm − l rs (13) Energies 2020, 13, 2914 6 of 25 and take differential Equation (13)  l rm − z s (14) The stabilizing function is defined by where c g1 and c g2 are two positive constants; α g = e g1 (v)dv is the integral factor that makes the tracking error to converge to zero. The virtual tracking error is defined by Take differential Equation (16) as Hence, as to design the control system, the external bunched force uncertainty f s4 is presumed to be bounded and meets the state of affairs f s4 ≤ γ 1 . f s1 , f s2 and f s3 are three unknown parameters. Three estimation errors can be defined by where e s1 , e s2 and e s3 are the estimation errors;f s1 ,f s2 andf s3 are the estimation values of f s1 , f s2 and f s3 . Then the Lyapunov function can be written by A s1 = e 2 g1 /2 + e 2 g2 /2 + e 2 s1 /(2η 1 ) + e 2 s2 /(2η 2 ) + e 2 s3 /(2η 3 ) + c g2 α 2 g /2 (21) By taking the differential of A s1 and by utilized (14)(15)(16)(17)(18)(19)(20) and the integral factor α g = e g1 (v)dv, then (21) In accordance with (22), the control propulsion u s of the backstepping control system with three adaptive laws and a beating function can be devised by where c g3 denotes a positive constant; γ 1 denotes a constant; γ 1 sgn(e g2 )/d s denotes the hitting function term. By utilized Equation (23), then Equation (22) can be rewritten by   (25) .f .f By utilized (25), (26), (27) and f s4 ≤ γ 1 , then Equation (24) can be rewritten by Equation (28) shows .
Utilized by the Barbalat's lemma [29,30], it can be displayed that lim t→∞ B s1 (t) = 0. It means that e g1 and e g2 will converge to zero as t → ∞ . Besides lim t→∞ l rs (t) = l rm and lim t→∞ z s = . l rm . Hence, the stability of the backstepping control system with three adaptive laws and a beating function can be guaranteed. Its block diagram is displayed in Figure 3.
Because the upper bound of the unknown external bunched force uncertainty f s4 is difficult to determine and the reckoned valuef s4 of the unknown external bunched force uncertainty f s4 is not precisely to reckon. Hence, the modified recurrent Rogers-Szego polynomials neural network (MRRSPNN) is proposed to suit the real value of the external bunched force uncertainty f s4 . The MRRSPNN with three-layer constitution is displayed in Figure 4, which makes up the first layer (input layer), the second layer (hidden layer 1) and the third layer (output layer). 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 explained by where x 1 1 = l rm − l rs = e g1 and x 1 2 = e g1 (1 − z −1 ) = ∆e g1 are the speed discrepancy and speed discrepancy alteration, respectively. R is the iteration count. v 1 ik is the recurrent weight through the third layer and the first layer. y 3 k is the output of node at the third layer. The symbol Π is a multiply factor. stability of the backstepping control system with three adaptive laws and a beating function can be guaranteed. Its block diagram is displayed in Figure 3. MRRSPNN with three-layer constitution is displayed in Figure 4, which makes up the first layer (input layer), the second layer (hidden layer 1) and the third layer (output layer). 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 explained by At the second layer, input semaphore and output semaphore are explained by where μ is the recurrent gain at the second layer. Rogers-Szego polynomials function [21,31] is adopted as the activation function. ) ; ( q x RS j is the Rogers-Szego polynomials function with are the 0-, 1-and Figure 3. Block diagram of the backstepping control system with three adaptive laws and a beating function.
Energies 2020, x, x FOR PEER REVIEW 9 of 27 2-order Rogers-Szego polynomials functions, respectively. The first kind type of Rogers-Szego polynomials function at the recurrence relation [21,31] is given by The symbol  is a summation factor.
At the third layer, semaphore and output semaphore are explained by where 2 kj v is the connecting weight through the second layer and the third layer. 3 k g is the linear activation function. The output at the third layer of the MRRSPNN can be denoted by are the weight vector at the third layer and the input vector at the third layer, respectively. is defined as where * Ο is an optimum weight vector that reaches the minimal rebuilt error. Hence, as to make up the minimal rebuilt error 4 s e , the reimbursed controller r u with a reckoned law is proposed.
The positive number . The Lyapunov function is chosen by where 5 η is an adaptive gain.   At the second layer, input semaphore and output semaphore are explained by where µ is the recurrent gain at the second layer. Rogers-Szego polynomials function [21,31] is adopted as the activation function. RS j (x; q) is the Rogers-Szego polynomials function with are the 0-, 1-and 2-order Rogers-Szego polynomials functions, respectively. The first kind type of Rogers-Szego polynomials function at the recurrence relation [21,31] is given by RS n+1 (x; q) = (1 + x)RS n (x; q) + x(q n − 1)RS n−1 (x; q). The symbol is a summation factor. At the third layer, semaphore and output semaphore are explained by where v 2 k j is the connecting weight through the second layer and the third layer. g 3 k is the linear activation function. The output y 3 k (R) at the third layer of the MRRSPNN can be denoted by where are the weight vector at the third layer and the input vector at the third layer, respectively.
The minimal rebuilt error e s4 is defined as where O * is an optimum weight vector that reaches the minimal rebuilt error. Hence, as to make up the minimal rebuilt error e s4 , the reimbursed controller u r with a reckoned law is proposed. The positive number σ s4 is greater than absolute value of e s4 , i.e., σ s4 ≥ |e s4 |. The Lyapunov function is chosen by where η 5 is an adaptive gain. e s5 =ê s4 − e s4 is defined as the reckoned error to be bounded.ê s4 is the reckoned value of minimal rebuilt error e s4 . By taking the derivative of A s2 and by utilized (14)(15)(16)(17)(18)(19)(20) and the integral factor α g = e g1 (v)dv, then (36) is obtained by In accordance with (37), the control propulsion u s =û s of the AMRRSPNNB control system with three adaptive laws and reimbursed controller with DGWO can be designed by By utilized (38), then (37) can be rewritten by By utilized Equations (25), (26), (27) and e s5 =σ s4 − σ s4 , then Equation (39) can be rewritten by To reach A s2 ≤ 0, the adaptive law .
O and the reimbursed controller u r with a reckoned lawσ s4 and an adaptive law . σ s4 to reduce uncertainties effects can be devised by . σ s4 = η 5 e g2 (43) By substituting (41-43) into (40) and by utilized σ s4 ≥ |e s4 |, then (40) can be presented by Equation (44) shows that A s2 (t) is a negative semi-definite function, i.e., A s2 (t) ≤ A s2 (0). It means that e g1 and e g2 are bounded. Utilized by the Barbalat's lemma, it can be presented that −B s1 (t) → 0 at t → ∞ from (29), (30), (44), i.e., e g1 and e g2 will converge to zero at t → ∞ . Thence, the stability of the AMRRSPNNB control system with three adaptive laws and reimbursed controller with DGWO can be ensured. Its block diagram is displayed in Figure 5. The fitted learning rule of the connecting weight utilized by the gradient descent skill with the chain rule is represented by where 2 δ denotes the learning rate. To acquire good delay, the DGWO is applied to look for two changeable learning rates in the MRRSPNN. In addition, for improving delay and looking for two optimal learning rates, the DGWO with two adjusted factors is proposed in this study.  Furthermore, the guaranteed delay of tracking error to be zero does not signify delay of the reckoned value of the unknown external bunched force uncertainty to it real values. The persistent excitation condition [29,30] should be satisfied for the reckoned value to converge to its theoretic value.
The fitted learning rule of the connecting weight utilized by the gradient descent skill with the chain rule is represented by It is known that ∂H 2 /∂y 3 k = −e g2 from Equation (45) and Equation (47). The fitted learning rule of recurrent weight v 1 ik utilized by the gradient descent technology with chain rule then is represented by where δ 2 denotes the learning rate. To acquire good delay, the DGWO is applied to look for two changeable learning rates in the MRRSPNN. In addition, for improving delay and looking for two optimal learning rates, the DGWO with two adjusted factors is proposed in this study.
In the DGWO, the optimization is conducted by α, β and ρ. The DGWO algorithm can be denoted by where G(l 1 + 1) = [δ 1 δ 2 ] is a vector that makes up two learning rates, G 1 (l 1 ), G 2 (l 1 ), G 3 (l 1 ) are denoted by where α(l 1 ), β(l 1 ), δ(l 1 ) are the three vectors as three best solutions; F 1 (l 1 ), F 2 (l 1 ), F 3 (l 1 ) and L 1 (l 1 ), L 2 (l 1 ), L 3 (l 1 ) are denoted by where ϕ 1 and ϕ 2 are two random values. The updated values in the two adjusted factors a 1 (l 1 ) and b 1 (l 1 ) can control the tradeoff between exploration and exploitation. Two adjusted factors a 1 (l 1 ) and b 1 (l 1 ) are updated by a 1 (l 1 ) = 2 − 2 cos((π l 1 )/(2I it1 )) (55) where l 1 is the iteration number; I it1 and I it2 are the total numbers of iteration allowed for the optimization. Finally, G(l 1 + 1) = [δ 1 δ 2 ] is the best solution in connection with the learning rates δ n (t), n = 1, 2 of the two weights in the MRRSPNN. Hence, the better numbers could be optimized by utilized DGWO with two adjusted factors yield two changeable learning rates for two weights to look for two optimal values and to speed-up delay of two weights. Figure 1 displays the control block diagram for the DSP control a permanent-magnet SLM drive system. Figure 6 displays a photo of the experimental setup. The DSP control system includes 2-channel encoder connection port and a 4-channel digital-analog converter (DAC). The proposed accuracy in the real-time control was realized by utilizing the DSP control system. Figure 7 displays the flowchart of the core program and the sub-core interrupt service routine (SCISR) utilized by the DSP control system. First, we deal with some parameters and input/output (I/O) ports initialization and the interrupt interval for the SCISR. When the interrupt was enabled, the core program was employed to observe all the data. The SCISR with a 2-ms sampling interval was employed for reading the mover position of the permanent-magnet SLM drive system from encoder and three-phase currents from A/D converter, computing reference model and position error, enforcing sin/cos generation and coordinate transformation, enforcing the AMRRSPNNB control system with three adaptive laws and reimbursed controller with DGWO, and outputting three-phase current commands to the interlocked and isolated circuits for switching three sets of insulated-gate-bipolar-transistor (IGBT) power modules of the pulse-width-modulation (PWM) voltage source inverter. A triangular-wave comparison current-controlled PWM with a switching frequency of 15 kHz was adopted to govern the PWM voltage source inverter. Additionally, about 0.1 kHz was adopted in the measured bandwidth of position control loop, and about 1 kHz was adopted in the measured bandwidth of current control loop for the permanent-magnet SLM drive system under no-load test. Hence, as to enhance precision of the sampling signals from A/D converter, the sampling interval in the control cycle at the tested results was set at 1 ms (1 kHz) so that the DSP control system has enough time to process control algorithm. Because of inherent uncertainty in the permanent-magnet SLM drive system (e.g., nonlinear friction, ending effects and time-varying dynamic uncertainties), current output limitation and voltage output limitation for DC bus power from AC source, the DC bus power only operated under maximum current and maximum voltage for avoiding burning out the IGBT power modules for the permanent-magnet SLM drive system. Furthermore, so as to prevent operation over-load and the PWM voltage source inverter burnout, the permanent-magnet SLM drive system displayed in Figure 1 had three sets of over current protection circuits, three sets of over voltage protection circuits and three sets of under-voltage protection circuits. The control system was realized by DSP control system. The control goal was to govern the mover to move 6 mm periodically. The reference model was set to be unit gain in the sinusoidal reference trajectory.       Figure 7. Flowchart of the enforcing program utilized by the DSP control system.

Test Results
For comparison control performance with the three control systems, five situations were demonstrated in the tested results. The three control systems were the renowned PI controller as the control system K1, the backstepping control system with three adaptive laws and a beating function as the control system K2, the AMRRSPNNB control system with three adaptive laws and reimbursed controller with DGWO as the control system K3. Five tested situations were as follows. Situation Q1 was the nominal situation under periodic step command from 0 mm to 6 mm. Situation Q2 was the parameters variations, the cogging force, the column friction force and the Stribeck effect force with 4 times the nominal value under periodic step command from 0 mm to 6 mm. Situation Q3 was the nominal situation under periodic sinusoidal command from −3 mm to 3 mm. Situation Q4 was the parameters variations, the cogging force, the column friction force and the Stribeck effect force with 4 times the nominal value under periodic sinusoidal command from −3 mm to 3 mm. Situation Q5 was the adding load force with step disturbance force F ls = 2 N at 6 mm situation.
Hence, as to reach good steady-state and transient-state responses, all the gains of the renowned PI controller as the control system K1 are k sp = 4.1, and k si = k sp /T si = 1.5 via certain heuristic knowledge [32][33][34] on the tuning of the PI controller at the nominal situation for the position tracking.
The tested results utilized by the control system K1 for governing the permanent-magnet SLM drive system at Situation Q1 and Situation Q2 are shown in Figure 8. The mover position responses at Situation Q1 and Situation Q2 are shown in Figure 8a,c, respectively. The associated control propulsion responses at Situation Q1 and Situation Q2 are shown in Figure 8b,d, respectively. The mover velocity responses at Situation Q1 and Situation Q2 are shown in Figure 8e,f, respectively. The tested results utilized by the control system K1 for controlling the permanent-magnet SLM drive system at Situation Q3 and Situation Q4 are shown in Figure 9. The mover position responses at Situation Q3 and Situation Q4 are shown in Figure 9a,c, respectively. The associated control propulsion responses at Situation Q3 and Situation Q4 are shown in Figure 9b,d, respectively. The mover velocity responses at Situation Q3 and Situation Q4 are shown in Figure 9e,f, respectively. The mover position responses with well tracking responses can be obtained utilized by the control system K1 at Situation Q1 and Situation Q3 displayed in Figures 8a and 9a, respectively. Furthermore, the mover position responses with worsen tracking responses displayed in Figures 8c and 9c were very evident in light of the great nonlinear commotion. The mover position responses with stagnant tracking responses were clear for governing the permanent-magnet SLM drive system utilized by the control system K1 from these tested results. Because of no adequate gain adjustment, the control system K1 has inferior robustness under the great nonlinear commotion.
The parameters of the backstepping control system with three adaptive laws and a beating function as the control system K2 are given as c g1 = 2.4, c g2 = 2.5, c g3 = 2.3, γ 1 = 8.2 according to heuristic knowledge [6,7,10] in the nominal situation for the position tracking under the periodic step command from 0 mm to 6 mm so as to reach good steady-state and transient-state responses. The tested results utilized by the control system K2 for governing the permanent-magnet SLM drive system at Situation Q1 and Situation Q2 are shown in Figure 10. The tested results utilized by the control system K2 for governing the permanent-magnet SLM drive system at Situation Q3 and Situation Q4 are shown in Figure 11. The mover position responses at Situation Q1 and Situation Q2 are shown in Figure 10a,c, respectively. The associated control propulsion responses at Situation Q1 and Situation Q2 are shown in Figure 10b,d, respectively. The mover velocity responses at Situation Q1 and Situation Q2 are shown in Figure 10e,f, respectively. The mover position responses at Situation Q3 and Situation Q4 are shown in Figure 11a,c, respectively. The associated control propulsion responses at Situation Q3 and Situation Q4 are shown in Figure 11b,d, respectively. The mover velocity responses at Situation Q3 and Situation Q4 are shown in Figure 11e,f, respectively. The mover position responses with well tracking responses can be obtained utilized by the control system K2 at Situation Q1 and Situation Q3 displayed in Figures  10a and 11a. Meanwhile, the mover position responses with fine tracking responses at Situation Q2 and Situation Q4 displayed in Figures 10c and 11c are evident under the great nonlinear commotion. From these tested results, the mover position responses with good tracking responses can be obtained for governing the permanent-magnet SLM drive system utilized by the control system K2 at Situation Q1, Situation Q2, Situation Q3 and Situation Q4. However, great upper bound with beating function causes to very crucial vibration in the control propulsion. Moreover, the vibration control propulsion will penetrate the bearing institution and may stimulate unstable dynamics occurrence.     Hence, as to display the usefulness of the control system with minimal number of neurons, the MRRSPNN adopts 2, 4 and 1 neuron in the first layer, the second layer and the third layer, respectively. The gains and parameters of the proposed AMRRSPNNB control system with three adaptive laws and reimbursed controller with DGWO as the control system K3 were as follows c g1 = 2.4, c g2 = 2.5, c g3 = 2.3, η 5 = 0.18 µ = 0.2 through some heuristic knowledge [6,7,10] in the nominal situation for the position tracking under the periodic step command from 0 mm to 6 mm so as to reach good steady-state and transient-state responses. The parameter's initialization in the MRRSPNN was adopted in [35] to obtain good initialize the parameters. In addition, the parameter-adjustment process keeps continually active in the tested duration.
The tested results utilized by the control system K3 for governing the permanent-magnet SLM drive system at Situation Q1 and Situation Q2 are shown in Figure 12. The tested results utilized by the control system K2 for governing the permanent-magnet SLM drive system at Situation Q3 and Situation Q4 are shown in Figure 13. The mover position responses at Situation Q1 and Situation Q2 are shown in Figure 12a,c, respectively. The associated control propulsion responses at Situation Q1 and Situation Q2 are shown in Figure 12b,d, respectively. The mover velocity responses at Situation Q1 and Situation Q2 are shown in Figure 12e,f, respectively. The mover position responses at Situation Q3 and Situation Q4 are shown in Figure 13a,c, respectively. The associated control propulsion responses at Situation Q3 and Situation Q4 are shown in Figure 13b,d, respectively. The mover velocity responses at Situation Q3 and Situation Q4 are shown in Figure 13e,f, respectively. The mover position responses with graceful tracking responses that can be obtained by utilizing the control system K2 at Situation Q1 and Situation Q3 displayed in Figures 12a and 13a. Moreover, the mover position responses with excellent tracking responses at Situation Q2 and Situation Q4 are shown in Figures 12c and 13c very conspicuous under the great nonlinear interferences. From these tested results, the mover position responses with better tracking responses are obtained for governing the permanent-magnet SLM drive system utilized by the control system K3 at Situation Q1, Situation Q2, Situation Q3 and Situation Q4.
Furthermore, the vibration was much decreased in the control propulsions utilized by the control system K3 at Situation Q1, Situation Q2, Situation Q3 and Situation Q4 as displayed in Figure 12b,d and Figure 13b,d. However, the robustness of the control performance utilized by the control system K3 under the occurrence of parameter variations for the tracking of periodical trajectories are evident in light of the online adaptive tuning of the MRRSPNN. The control performance utilized by the control system K3 has better control performances than the control system K2 from the tested results.
Finally, tested results for the mover position-regulating responses under adding load force with step disturbance force F ls = 2 N at 6 mm as Situation Q5 are shown in Figure 14. Tested results with the mover position-regulating responses at Situation Q5 utilized by the control system K1, the control system K2 and the control system K3 at Situation Q5 are shown in Figure 14a-c, respectively. From these tested results, performance ability on load force regulation utilized by the control system K3 has better transient response than utilized by the control system K1 and the control system K2. However, the robustness of control performance utilized by the control system K3 was wonderful for governing the permanent-magnet SLM drive system under the occurrences of parameter disturbance and load force interference in light of the online adaptive adjustment of the MRRSPNN with DGWO. From these tested results, control performance responses of the control system K3 was superior to the control system K1 and the control system K2.   K3 has better transient response than utilized by the control system K1 and the control system K2. However, the robustness of control performance utilized by the control system K3 was wonderful for governing the permanent-magnet SLM drive system under the occurrences of parameter disturbance and load force interference in light of the online adaptive adjustment of the MRRSPNN with DGWO. From these tested results, control performance responses of the control system K3 was superior to the control system K1 and the control system K2. In addition, some comparative performances utilized by the control systems K1, K2 and K3 are summarized in Table 1 with regarding to the five tested situations. The maximum errors of      In addition, some comparative performances utilized by the control systems K1, K2 and K3 are summarized in Table 1 with regarding to the five tested situations. The maximum errors of e g1 by utilizing the control systems K1, K2 and K3 at Situation Q1 are 0.55 mm, 0.32 mm and 0.18 mm, respectively. The root mean square errors of e g1 utilized by the control systems K1, K2 and K3 at Situation Q1 are 0.31 mm, 0.21 mm and 0.08 mm, respectively. The maximum errors of e g1 utilized by the control systems K1, K2 and K3 at Situation Q2 are 0.73 mm, 0.41 mm and 0.21 mm, respectively. The root mean square errors of e g1 utilized by the control systems K1, K2 and K3 at Situation Q2 are 0.41 mm, 0.31 mm and 0.11 mm, respectively. The maximum errors of e g1 utilized by the control systems K1, K2 and K3 at Situation Q3 are 0.53 mm, 0.31 mm and 0.18 mm, respectively. The root mean square errors of e g1 utilized by the control systems K1, K2 and K3 at Situation Q3 are 0.31 mm, 0.21 mm and 0.08 mm, respectively. The maximum errors of e g1 utilized by the control systems K1, K2 and K3 at Situation Q4 are 0.72 mm, 0.40 mm and 0.20 mm, respectively. The root mean square errors of e g1 utilized by the control systems K1, K2 and K3 at Situation Q4 are 0.40 mm, 0.30 mm and 0.10 mm, respectively. The maximum errors of e g1 utilized by the control systems K1, K2 and K3 at Situation Q5 are 0.79 mm, 0.42 mm and 0.20 mm, respectively. The root mean square errors of e g1 utilized by the control systems K1, K2 and K3 at Situation Q5 are 0.56 mm, 0.28 mm and 0.11 mm, respectively. As displayed in the table, the control systems K3 produces in smaller tracking error in comparison with the control systems K1 and K2. In the light of the tabulated measurements, the control systems K1 indeed yields the beneficial control performance.

performance.
Furthermore, the feature performance comparisons utilized by the control systems K1, K2 and K3 are summarized in Table 2 from tested results. In Table 2, some feature performances with regarding to the vibration in the control propulsion, the dynamic response, the regulation ability of load force disturbance, the convergent speed, the position tracking error and the rejection ability of parameter disturbance utilized by the control systems K3 are superior to utilized by the control systems K1 and K2. Furthermore, the feature performance comparisons utilized by the control systems K1, K2 and K3 are summarized in Table 2 from tested results. In Table 2, some feature performances with regarding to the vibration in the control propulsion, the dynamic response, the regulation ability of load force disturbance, the convergent speed, the position tracking error and the rejection ability of parameter disturbance utilized by the control systems K3 are superior to utilized by the control systems K1 and K2.  Learning rate None None Vary (two optimal learning rates)

Conclusions
A proposed AMRRSPNNB control system with three adaptive laws and reimbursed controller with DGWO was proposed for governing a permanent-magnet SLM drive system under the conditions of parameter disturbance and load-force interference for the tracking of periodical trajectories.
The significant contributions are as follows: (1) The FOC was successfully used for governing

Conclusions
A proposed AMRRSPNNB control system with three adaptive laws and reimbursed controller with DGWO was proposed for governing a permanent-magnet SLM drive system under the conditions of parameter disturbance and load-force interference for the tracking of periodical trajectories.
The significant contributions are as follows: (1) The FOC was successfully used for governing the permanent-magnet SLM drive system; (2) the backstepping control system with three adaptive laws and a beating function was successfully derived in light of the Lyapunov function to decrease influence under the external lumped force uncertainty disturbances and load force interferences; (3) the proposed AMRRSPNNB control system with three adaptive laws and reimbursed controller with DGWO to reckoned the external lumped force uncertainty and reimburse the reckoned error was successfully derived in the light of the Lyapunov function for reducing the external lumped force uncertainty disturbances and the load force interferences; (4) two optimal learning rates of the MRRSPNN were successfully regulated by utilizing the DGWO algorithm to speed up the parameter's contraction. Furthermore, the proposed AMRRSPNNB control system with three adaptive laws and reimbursed controller with DGWO had better control performances than the renowned PI controller and the backstepping control system with three adaptive laws and a beating function in various feature performances as the vibration in the control propulsion, the convergent speed, the position tracking error, the dynamic response, the regulation ability of load force disturbance and the rejection ability of parameter disturbance from Tables 1 and 2.