Electromagnetic Actuator System Using Witty Control System

: Electromagnetic actuator systems composed of an induction servo motor (ISM) drive system and a rice milling machine system have widely been used in agricultural applications. In order to achieve a ﬁner control performance, a witty control system using a revised recurrent Jacobi polynomial neural network (RRJPNN) control and two remunerated controls with an altered bat search algorithm (ABSA) method is proposed to control electromagnetic actuator systems. The witty control system with ﬁner learning capability can fulﬁll the RRJPNN control, which involves an attunement law, two remunerated controls, which have two evaluation laws, and a dominator control. Based on the Lyapunov stability principle, the attunement law in the RRJPNN control and two evaluation laws in the two remunerated controls are derived. Moreover, the ABSA method can acquire the adjustable learning rates to quicken convergence of weights. Finally, the proposed control method exhibits a ﬁner control performance that is conﬁrmed by experimental results.


Introduction
Compared to other three-phase motors, three-phase induction motors (IMs) are widely used in many industrial and commerce applications due to their simple structures and easy maintenance. In order to achieve better control performance, IMs have served as induction servo motors (ISMs) via structural improvement and encoder installation. Therefore, ISMs have been broadly applied to various servo fields such as computer numerical control (CNC) machine tools and milling machines [1][2][3][4]. Li et al. [1] proposed a new intelligent adaptive CNC system design for a milling machine by using the neural network controller to achieve better control characteristics. Huang et al. [2] proposed an approach for cutting the force control of CNC machines. This approach with a state estimator was executed by using the observed variables and cutting force to achieve robust control. Recently, the developed approach was applied to a milling machine center. Gomes and Sousa [3] proposed the adaptive control of milling machine cutting force by using an artificial neural network to obtain good results. Mikolajczyk [4] proposed a system of a numerical control conventional milling machine with electromagnetic clutches by using VB6 special software to control the machine with G-code. However, these control systems took a long time to fulfill nonlinear disturbances so that they resulted in lower calculation efficiency. Thereby, the aim of the proposed witty control system using a revised recurrent Jacobi polynomial neural network (RRJPNN) control and two remunerated controls with altered bat search algorithm (ABSA) and progressive weight pruning approaches for the ISM driving the rice milling machine is to reduce computing time and to quicken convergence of weights. Meanwhile, the proposed witty control system can increase machining efficiency and control characteristics.
Due to their good learning capability, many neural networks (NNs) have been applied in many linear and nonlinear systems, such as optimal hysteresis modelling methods of magnetic controlled shape memory alloys (MSMAs) [5], control of ionic electroactive polymer actuators [6], predicting the driving force for a multicyclic driving experiment of silicone/ethanol soft composite material actuators [7], and design of a 2-DOF ankle exoskeleton with a bidirectional tendon-driving actuator and a polycentric structure [8]. These NNs take a long time to conduct training and learning of a system. In order to reduce calculating time, some polynomial function NNs have been proposed for applications in function approximation [9], dynamic control for continuously variable transmission [10], approximation and estimation of nonlinear functions [11] and backstepping control of synchronous reluctance motor (SynRM) drive systems [12]. As heuristic comprehension methods were applied to adjust learning rates of weights, these NNs that were applied in control systems resulted in slower convergence. Thereby, the proposed NN that is combined with Jacobi polynomials [13] is a novel NN structure that has not yet been applied to estimate, predict and control nonlinear systems. The feedforward Jacobi polynomial neural network (FJPNN) [14] may not effectively approximate nonlinear dynamics because there is no recurrent path. Therefore, many recurrent NNs have been applied in nonlinear system identification [15], model predictive control for continuous pharmaceutical manufacturing [16], estimation in the effective connectivity of electroencephalography [17] and visual field prediction [18] because of higher accuracy and better identification. Due to it having more advantages than the feedforward FJPNN, the RRJPNN control by adopting the ABSA method with progressive weight pruning approach has not yet been used to control the ISM driving the rice milling machine system in order to cut down computation complexity.
The customary optimization algorithm is applied when solving the classic problem of smaller dimensions that are not easily applicable in reality. The swarm intelligence optimization algorithms (SIOAs) were discovered by adopting simulated natural biological systems. Recently, many researchers have proposed the SIOAs, such as bat algorithm (BA) [19], particle swarm optimization (PSO) [20], social learning optimization algorithm (SLOA) [21], chicken swarm optimization algorithm (CSOA) [22] and ant colony algorithm (ACO) [23]. Yang [24] first proposed the bat algorithm (BA), which is based on the swarm intelligence heuristic search algorithm. Due to it having fewer parameters, powerful robustness, and a simple and easy implementation, BA has attracted more and more attention in the search for a global optimal solution. The BA has been applied in multiobjective function optimizations with a neural network model [19], economic operations [25] and numerical optimizations [26]. However, the local search method by adopting BA has the shortcoming of precocious convergence and does not perform well in the early search stage. Therefore, the proposed altered bat search algorithm (ABSA) method is a novel method to avoid precocious convergence. Moreover, a progressive weight pruning approach based on the ABSA method is proposed to deal with nonconvex optimization problems. The ABSA method with modified loudness and modified pulse rate is used to adjust two optimal learning rates. Thereby, a novel fast-convergence algorithm applied to RRJPNN with two learning rates by adopting the ABSA method with the progressive weight pruning approach that is used to adjust two optimal learning rates and to quicken convergence of weights is proposed in this paper. At last, some tested results show that the fine control performances are confirmed by the proposed control method.
The organization of this research is as follows. Section 2 introduces some explanations of the ISM driving the rice milling machine system. Section 3 advocates the witty control system by adopting the RRJPNN control and two remunerated controls with ABSA. Section 4 promotes the tests and experimental results for the ISM driving the rice milling machine system. Section 5 presents the discussions and explanations of some experimental results. Section 6 describes some conclusions.

Materials
The complete system consists of three subparts that are: the ISM and drive system, the digital signal processor (DSP) control system and the rice milling machine system, as below.

System Description
The ISM and drive system consist of a mounted encoder three-phase ISM with a rotor that adopts a high moment of inertia and low frictional coefficient, a rectifier, an LC filter, an isolated circuit, a lockout-time circuit, a current sensing circuit, an analog-digital conversion, and the voltage fed converter with 3-leg 6-switch transistor power devices. The DSP control system is composed of a speed control, a proportional-integral (PI) current control and an indirect field-oriented control (IFOC) [27][28][29]. The IFOC consists of a spacevector pulse-width-modulation control, sin θ t / cos θ t creation, a coordinate transformation and an inverse coordinate transformation. The rice milling machine system is composed of a feed chopper, milling room, main idler, idler 1, idler 2, chaff suction fan, jet fan, chaff outlet, inlet chopper mouth, rice outlet, thickness adjustment lever, belt 1 and belt 2. The arrangement of the ISM driving the rice milling machine system is shown in Figure 1.

System Model
The voltage equations in the three-phase ISM by using the simplified coordinate frame transforms, from a-b-c to u-v, and the Clarke and Park transformations by using the IFOC [27][28][29] can be represented by where w u and w v are the u − v axis stator voltages. q u and q v are the u − v axis stator currents. q ur and q vr are the u − v axis rotor currents. L u , L v and L uv are the u − v axis self-inductances and mutual inductance, respectively. r a and r b are the stator and equalized rotor resistances. µ v and µ u are the mechanical and electrical angular speeds in the ISM. µ t is the electrical angular speed of synchronous flux [27]. The electromagnetic torque F a [28,29] in the ISM can be described as where ϕ u and ϕ v are u − v axis flux linkages. P t is the number of poles. The ISM driving the rice milling machine system led to a more sluggish performance of the system owing to a nonlinear uncertainty effect. The response of speed control for the ISM driving the rice milling machine system resulted in poor performance. The adumbration view of the ISM and the rice milling machine system is illustrated in Figure 2. Considering that the power loss and sliding loss were insignificant, two dynamic equations with simplified kinematics of the rice milling machine system are described as [30][31][32]: where F v , F d and F b1 are the electromagnetic torque of the ISM, the output torque of idler 2 and the output torque of the main idler, respectively. f a , f c , f d and f b are the four moments of inertia in the ISM, in idler 2, in the main idler and in idler 1, respectively. g v , g c , g d and g b are the four viscid frictional coefficients in the ISM, in idler 2, in the main idler and in idler 1, respectively. φ b is the transposition ratios regarding idler 2 and the main idler for the rice milling machine system. F l b (R vb , F vb , F vl , g b ) is the nonlinear coalescence disturbances function including rolling force R vb , wind force F vb , and braking force F vl . µ v and µ b are the speed in idler 2 and the speed in the main idler. Then, the torque equation can be transformed from the main idler to idler 2 by use of a transformed ratio. The modeling of the rice milling machine can be simplified by omitted sliding losses of two belts; thus, the dynamic equation in the ISM driving the rice milling machine system including the coalescence torque from Equations (6) and (7) can be expressed as where g r = g a + g c and f r = f a + f c are the coalescence viscid friction coefficient and the coalescence moment of inertia including the main idler and the ISM. ∆F a + F b1 is the huge comprehensive coalescence disturbances and parameter variations. F t (F 1c , F 2c , F 3c ) is the coalescence torque [27] including coulomb friction torque F 3c , Stribeck effect torque F 2c and adding load torque F 1c . F b1 represents the comprehensive coalescence disturbances. ∆F a represents the comprehensive parameter variations. The DSP control system with current control and IFOC can fulfill a speed control, an IFOC and a proportional-integral (PI) current control. The IFOC consists of a space-vector pulse-width-modulation control, sin θ t / cos θ t , creation, a coordinate transformation and an inverse coordinate transformation. The control gains of the PI current control are the proportional gain of 19.2 and the integral gain of 8.3 by using the heuristic method [33][34][35] to obtain a finer dynamic response. The drive system was operated under comprehensive coalescence disturbances and comprehensive parameter variations by adopting the DSP control system in this research.

Methods
In order to design the control structure, the dynamic equation of Equation (8) is modified as where ∆F a + F b1 + F e represents the comprehensive coalescence disturbances.
is a friendly constant concerning the coalescence moment of inertia and |h w F o | ≤ R b is assumed to be bounded. h x = f −1 r is a friendly constant concerning the coalescence moment of inertia and h x ≤ R c is assumed to be bounded. R b and R c are two friendly values. l v = F v is the electromagnetic torque of the ISM. The speed difference s a is as follows.
If the comprehensive coalescence disturbances and the comprehensive parameter variations are favorable and affectionate, the excellent control law can be rewritten by where c v is a positive control gain. Equation (12) and l * v = l v are substituted into Equation (11), so the error equation can be rewritten by The system will track the wished state value at t → ∞ and s a (t) → 0 . Nevertheless, the control system will exhibit a sluggish tracking response under the occurrence of uncertainty. Thereby, the proposed witty control system using an RRJPNN control and two remunerated controls with ABSA shown in Figure 3 were developed to control the ISM driving the rice milling machine system in order to enhance the speed of the tracking response.
We can use the differential equation of (11), then we can substitute Equations (8) and (12) into this equation. The error equation can be rewritten by where l x is the RRJPNN control that acts as the main tracking controller to impersonate the excellent control rule. l y is the dominator control that will act in the appropriate region. l z and l w are two remunerated controls that acts as two remunerated controllers to repay the difference between the excellent control and the RRJPNN control. Then, the three-layer RRJPNN, which is shown in Figure 4, consists of the forehead, center and readward layers.
All the informations of all layers are as follows. The input and the output informations in the forehead layer are shown below.
The input and the output informations in the center layer are shown below.
The input and the output informations in the readward layer are shown below.
where Π and Σ are the multiplication and summation symbols. a 1 2 = s a (1 − z −1 ) = ∆s a and a 1 1 = µ * −µ v = s a are the speed difference alteration and the speed difference. m, χ and K are the number of nodes of the center layer, the recurrent gain of the center layer and the iteration number, respectively. ρ 1 rt (K) and ρ 2 ts (K) are the recurrent weight between the readward layer and the forehead layer, and the conjoined weight between the center layer and the readward layer. y 1 r , y 2 s and y 3 t are the three linear activation functions in the forehead, center and readward layers, respectively. d 1 r (K), d 2 s (K) and d 3 t (K) are the information of three outputs of nodes in the forehead, center and readward layers, respectively. P (α,β) s (b) is the Jacobi polynomial function [13,14] with −1 < x < 1 adopted as the activation function in the center layer-i.e., y 2 2 are the 0-, 1-and 2-order Jacobi polynomial functions, respectively. The Jacobi polynomial function with the recurrence relation [13,14] is as below (2s . The output information in the readward layer can be rewritten as d 3 where T are the input information and weight vectors in the readward layer. The dominator control l y can be represented by The dominator control will act in the appropriate region if the RRJPNN control cannot be guaranteed.
To fulfill the remunerated mechanism, a minimum difference λ can be described as δ are greater than zero. By using Equation (22), l x = d 3 t (K) = C T D and l * x = (C * ) T D, then Equation (16) can be described as To obtain two remunerated controls, the attunement law and the two evaluation laws, the Lyapunov function is described as where υ 1 is the learning rate of the conjoined weight; τ 1 and τ 2 are the two positive evaluation rates;γ − γ andδ − δ are the evaluation differences. By using Equations (22) and (23), then the differential equation of (24) can be described as The attunement law dC dt , two remunerated controls, l z and l w , and two evaluation laws, dγ/dt and dδ/dt to fulfill dV x dt ≤ 0 can be described as By using Equations (26)- (30) and (17), then Equation (25) can be rewritten by By using |λ| < γ < δ and δ + ( s a and (C * − C) are represented as bounded when dV x dt ≤ 0, which is a negative semidefinite. Additionally, the uniformly continuous function x a (t) can be described by The integral of x a (t) can be rewritten by The differential of Equation (33) can be described as The limitation of Equation (34) when V x (0) and V x (t) are bounded can be described as As all variables on the right side of Equation (23) are bounded, which implies ds a dt is also bounded and it can be shown that lim t→∞ t 0 x a (ε)dε = 0, thus s a (t) → 0 as t → ∞ by using Barbalat's lemma [36,37]. Therefore, the proposed witty control system is gradually stable from proof. Moreover, the tracking error s a (t) of the system will converge to zero. Therefore, to describe the online training process of the RRJPNN, an objective function can be defined by The conjoined weight by using the backpropagation technology and the gradient descent technology from the attunement law dC dt can be described as By using the two above technologies, the conjoined weight from Equation (38) can be expressed by In comparison with Equations (38) and (39), this can be obtained as The updated law of the conjoined weight can be denoted by [38] By using the two above technologies, the recurrent weight of the attunement law can be expressed by where υ 2 is the learning rate of the recurrent weight. The updated law of the recurrent weight can be denoted by [38] Moreover, the ABSA method was applied to search for two optimal adjustable learning rates and to improve convergent speed of the weights in the RRJPNN in this research. The existing algorithm with the excellent behavior of miniature bat echolocation is the important development of the BA method. By adopting a random technology, this algorithm produces a set of solutions. Then, the optimal solution is found by using the loop search. The local solution is generated by random flight and generates a global optimal solution. For all bats, the position of the bat i is z i (n − 1), the flight velocity is f i (n − 1) and the current global optimal position is z * when their foraging space is part of the d-dimension at n − 1 time.
The flight velocity f i (n) and position z i (n) of bat i at n time can be calculated by where N max is the maximum number of iterations; k max and k min are the maximum and minimum frequencies of the soundwaves produced by the bat. In the initial process, the frequency of the bat's soundwaves is uniformly distributed between k min = 0 and k max = 1. The concerning frequency is obtained by adopting Equation (44). By adopting Equations (45) and (46), the local search is realized. The bat randomly goes along the optimal solution, and the new solution is updated by where σ is a random number at [−1, 1]. z old (n − 1) is the solution selected from the current optimal solution by adopting a random skill. d i (n) is the average loudness from the bat generation at n time. Additionally, it achieves a global search by controlling modified loudness d i (n + 1) and modified pulse rate e i (n + 1). The modified pulse rate e i (n + 1) and modified loudness d i (n + 1) of the bat launch pulse can be updated by where e i (0) is an initial rate and d i (0) is an initial loudness. ς is a constant between 0 and 1 and ξ is a positive constant. When the bat is conscious of the presence of the prey, it will increase its pulse emission rate and reduce the response of its pulsed emission. At last, z i (n), i = 1, 2 is the best solution concerning the learning rates υ i (n), i = 1, 2 of the two weights in the RRJPNN. Thereby, the two adjustable values may be optimized by adopting ABSA method to find the two learning rates of weights. Moreover, a progressive weight pruning approach is based on the ABSA method to quicken convergence of weight.

Tests and Results
The arrangement of the ISM driving the rice milling machine system by adopting a DSP control system is illustrated in Figure 1. In Figure 5, an experimental photo of the ISM and the rice milling machine system is illustrated. The conversion ratio for the rice milling machine system is 2.2. The profile formats are the belt 1 length is 42.2 mm, the belt 2 length is 52.2 mm, the main idler diameter is 92.6 mm, the diameter of idler 1 is 45.2 mm, and the diameter of idler 2 is 64.2 mm. The specification of the ISM is three-phase two-pole 220 V, 60 Hz, 3 kW, 3582 rpm. The position and speed conversion ratios are 1 V = 50 rad and 1 V = 50 rad/s. The internal parameters of the ISM are r a = 1.08 Ω, r b = 1.02 Ω, L u = 8.65 mH, L v = 10.68 mH, L uv = 8.89 mH, f a = 16.22 × 10 −3 Nms 2 = 0.811 Nmsrad/V, and g a = 1.12 × 10 −3 Nms/rad = 0.056 Nm/V. Figure 6 illustrates the control flowchart of executive program by adopting the DSP control system. The program in the experimental tests consists of the basic program (BP) and the auxiliary interrupt routine (AIR). The BP conducts all initializations for the adopted parameters and all settings for the input/output interfaces. The AIR achieves the interrupt execution within 2 ms. The executed processes by AIR are as follows: three-phase currents read from analog-digital conversions, rotor position read from encoder interface circuit, rotor speed computation, speed difference computation, lookup table generation, coordinate transformations realization, PI current control realization, the proposed control system realization, and three-phase space-vector pulse-width-modulation outputs for switching the voltage fed converter. Three discerners A1_x, A1, A and 2 are set as 3, 0 and 0, respectively. The DSP control system with the IFOC applied the discerner A2 to act as the executing number of the proposed control method. If the IFOC is executed less than three times, i.e., A1 < A1_x, the IFOC needs to be enforced repeatedly. The proposed control method is executed one time and the IFOC is executed three times. Then, the AIR will back to the BP.
The experimental results with three test examples are shown to show some of the control performances. Firstly, test JA is the case with huge comprehensive coalescence disturbances and parameter variations ∆F a + F b1 at 1.2 s, starting with a mandate speed of 1600 rpm (167.47 rad/s). Secondly, test JB is the case with double huge comprehensive coalescence disturbances and parameter variations ∆F a + F b1 at 1.2 s, starting with a mandate speed of 3300 rpm (345.40 rad/s). Thirdly, test JC has a mandate speed of 2000 rpm (209.33 rad/s) starting at 2 s and a mandate speed of 3000 rpm (314.00 rad/s) at 10 s, with acceleration and added external load torque disturbance and parameter variations 8 Nm(F 1c ) + F b1 at 14 s, with a mandate speed of 3000 rpm (314.00 rad/s). The PI controller as the TA controller and the proposed witty control system as the TB controller are the two adopted controllers that are compared with control performances. Firstly, two gains of the PI control as the TA controller are the proportional gain of 24.1 and the integral gain of 10.2 by adopting the heuristic method [30][31][32] to obtain finer dynamic response under the requirement of stability consideration. The control gains by using the proposed witty control system as the TB controller are given as c v = 5.41, χ = 0.092, τ 1 = 0.12 and τ 2 = 0.13 to better measure transient performance under the demand of stability planning. In addition, the number of nodes of the RRJPNN, by adopting the progressive weight pruning approach based on ABSA method to quicken convergence of conjoined weight, are 2, 4 and 1 in the forehead, center and readward layers to better measure transient-state and steady-state control properties.
Firstly, Figure 7a,b display speed responses for measured speed µ v , mandate speed µ c and reference model speed µ * via experimental results of test JA by adopting the TA and TB controllers for the ISM driving the rice milling machine system.    Figure 9a,b display the responses of three-phase currents via experimental results at test JA by adopting the TA and TB controllers for the ISM driving the rice milling machine system. Figure 7a,b show that by adopting the TA and TB controllers in test JA, a better speed tracking performance was achieved because of smaller disturbance. However, the increase in tracking error when adopting the TA controller shown in Figure 8a is larger than the increase in tracking error when adopting the TB controller shown in Figure 8b. The response of three-phase currents when adopting the TA controller shown in Figure 9a generates a larger harmonic wave than when adopting the TB controller shown in Figure 9b.
Secondly, Figure 10a,b display speed responses for measured speed µ v , mandate speed µ c and reference model speed µ * via experimental results of test JB by adopting the TA and TB controllers for the ISM driving the rice milling machine system.     Figure 12a,b display responses of three-phase currents via experimental results of test JB by using the controllers TA, and TB for the ISM driving the rice milling machine system. Figure 10a by using the TA controller in test JB appeared to show a dilatory speed response because of no good gain adjustment in the TA controller. In Figure 10b, by adopting the TB controller at test JB, a good speed response is demonstrated, owing to online adjustable method of RRJPNN control and two remunerated controls. However, the increase in tracking error when adopting the TA controller shown in Figure 11a is larger than the increase in tracking error when adopting the TB controller shown in Figure 11b. The response of three-phase currents when adopting the TA controller shown in Figure 12a results in larger harmonics than when adopting the TB controller shown in Figure 12b.
Thirdly, Figure 13a,b display two various speed-regulated responses when adding load torque via experimental results of test JC by using the TA and TB controllers for the ISM driving the rice milling machine system.    The two speed-regulated responses when adding load torque via experimental results by adopting the TA controller, shown in Figure 13a, are worse than the two speed-regulated responses when adding load torque via experimental results by using the TB controller, also shown in Figure 13a. Responses of three-phase currents with adding load torque via experimental results by adopting the TA controller shown in Figure 14a has a greater harmonic than responses of three-phase currents when adding load torque via experimental results by adopting the TB controller shown in Figure 14b.
Moreover, responses of the two learning rates curves in test JB shown in Figure 15a,b using calculated learning rates according to the proposed ABSA method are compared to two learning rates of conjoined weight and recurrent weight by utilizing the PSO method [20] and ACO method [23] to demonstrate the usefulness of this novel technique. This study shows that convergence to optimal values can be achieved by using the proposed ABSA method. The proposed method also achieves faster convergence and less computational complexity.
Additionally, responses of two weights at test JB by using the PSO method [20] and the ACO method [23] and the proposed ABSA method are shown in Figure 16a,b. The convergences of conjoined weight and recurrent weight by using the proposed ABSA method with progressive weight pruning approach are superior to the ACO method and the PSO method. Thereby, the proposed method with progressive weight pruning approach can quicken convergence of weights.  Furthermore, response of the numbers of conjoined weight at test JB by using the PSO method [20] and the ACO method [23] and the proposed ABSA method by adopting the progressive weight pruning approach is shown in Figure 17. The convergence of numbers of conjoined weight by using the proposed ABSA method and the progressive weight pruning approach is superior to the ACO method and the PSO method. Thereby, the proposed method can achieve faster convergence in conjoined weight.

Analyses and Discussion
Dynamic responses for the PI controller as the TA controller and the proposed witty control system as the TB controller at three tested examples via experimental results that are listed in Table 1 are explained as below. For test JA, the maximum differences for the TA and TB controllers are 82 (8.58 rad/s) and 30 rpm (3.14 rad/s), and the quadratic mean differences for the TA and TB controllers are 48 (5.02 rad/s) and 17 rpm (1.78 rad/s). For test JB, the maximum differences for the TA and TB controllers are 128 (13.40 rad/s) and 35 rpm (3.66 rad/s), and the quadratic mean differences for the TA and TB controllers are 53 (5.55 rad/s) and 19 rpm (1.99 rad/s). For test JC, the maximum differences for the TA and TB controllers are 489 (51.18 rad/s) and 192 rpm (20.10 rad/s), and the quadratic mean differences for the TA and TB controllers are 188 (19.68 rad/s) and 46 rpm (4.81 rad/s). The TB controller has better dynamic responses than the TA controller according to the experimental results in tests JA, JB and JC.  (13.40 rad/s) in maximum difference and 35 rpm (3.66 rad/s) in maximum difference. The above performances concerning the harmonic values in the three-phase currents, the dynamic responses, the regulation capabilities for adding load torque, the speed tracking differences and the denial potentialities of parameter disturbance in the TB controller are better than the TA controller. Thereby, the TB controller has better peculiarity performance than the TA controller from experimental results of tests JB and JC.

Conclusions
The proposed witty control system has been applied to control the the ISM driving the rice milling machine system with better robustness. The proposed witty control system that can realize the RRJPNN control, which involves an attunement law, two remunerated controls, which have two evaluation laws, and a dominator control were proposed to obtain a fine control performance.
The contributions of this research are as below. (a) The dynamic models of the ISM driving the rice milling machine system have been developed. (b) The ISM driving the rice milling machine system under huge comprehensive nonlinear synthesized disturbances and parameter variations affect has been controlled by using the proposed witty control method. (c) On the basis of the Lyapunov stability principle, the attunement law in the RRJPNN control and the two evaluation laws in the two remunerated controls have been developed. (d) The ABSA method was utilized to find the learning rates of conjoined and recurrent weights in the RRJPNN to obtain optimal values and to quicken convergence of weights. (e) The proposed witty control system has better sinusoidal shapes than the PI control in terms of the harmonics values of three-phase currents.
Finally, all results show that the proposed witty control system is better than the PI control for the ISM driving the rice milling machine system from all experimental results and control behaviors.

Conflicts of Interest:
The authors declare no conflict of interest.
Abbreviations w u , w v u − v axis stator voltages q u , q v u − v axis stator currents q ur , q vr u − v axis rotor currents L u , L v , L uv u − v axis self inductances, mutual inductance r a , r b stator and equalized rotor resistances.
µ v , µ u , µ t mechanical and electrical angular speeds, electrical angular speed of synchronous flux in the ISM F a electromagnetic torque ϕ u , ϕ v u − v axis flux linkages P t number of pole F v , F d , F b1 electromagnetic torque of the ISM, the output torque of idler 2 and the output torque of the main idler f a , f c , f d , f b four moments of inertia in the ISM, in idler 2, in the main idler and in idler 1 g v , g c , g d , g b four viscid frictional coefficients in the ISM, in idler 2, in the main idler and in idler 1 φ b transposition ratios regarding idler 2 and the main idler for the rice milling machine system F l b (R vb , F vb , F vl , g b ) nonlinear coalescence disturbances function R vb , F vb , F vl rolling force, wind force, braking force µ v , µ b speed in idler 2 and the speed in the main idler. g r = g a + g c coalescence viscid friction coefficient including the main idler and the ISM f r = f a + f c coalescence moment of inertia including the main idler and the ISM ∆F a + F b1 huge comprehensive coalescence disturbances and parameter variations F t (F 1c , F 2c , F 3c ) coalescence torque F 3c , F 2c , F 1c coulomb friction torque, Stribeck effect torque, adding load torque F b1 comprehensive coalescence disturbances ∆F a comprehensive parameter variations ∆F a + F b1 + F e comprehensive coalescence disturbances bounded with functional-bounded value h w = − f −1 r friendly constant concerning the coalescence moment of inertia h x = f −1 r friendly constant concerning the coalescence moment of inertia R b , R c two friendly values with bound l v = F v electromagnetic torque of the ISM s a speed difference c v positive control gain l x , l y , l z , l w RRJPNN control, dominator control, two remunerated controls a 1 2 = s a (1 − z −1 ) = ∆s a , a 1 1 = µ * −µ v = s a speed difference alteration, speed difference m, χ and K node number of the center layer, the recurrent gain of the center layer and the iteration number ρ 1 rt (K), ρ 2 ts (K) recurrent weight, conjoined weight y 1 r , y 2 s , y 3 solution selected from the current optimal solution at n − 1 time d i (n) average loudness from the bat generation at n time d i (n + 1), e i (n + 1) modified loudness at n + 1 time, modified pulse rate at n + 1 time e i (0), d i (0) initial rate, initial loudness ς, ξ constant between 0 and 1, positive constant