Self-Organizing Interval Type-2 Fuzzy Neural Network Compensation Control Based on Real-Time Data Information Entropy and Its Application in n-DOF Manipulator

In order to solve the high-precision motion control problem of the n-degree-of-freedom (n-DOF) manipulator driven by large amount of real-time data, a motion control algorithm based on self-organizing interval type-2 fuzzy neural network error compensation (SOT2-FNNEC) is proposed. The proposed control framework can effectively suppress various types of interference such as base jitter, signal interference, time delay, etc., during the movement of the manipulator. The fuzzy neural network structure and self-organization method are used to realize the online self-organization of fuzzy rules based on control data. The stability of the closed-loop control systems are proved by Lyapunov stability theory. Simulations show that the algorithm is superior to a self-organizing fuzzy error compensation network and conventional sliding mode variable structure control methods in control performance.


Introduction
The n-DOF manipulator, as the earliest application of robotic systems, still plays an important role in various fields such as industry, the military, medicine, aerospace, and so on. In the context of the future Industry 4.0 era, we have higher and higher requirements for the control accuracy of n-DOF manipulators [1,2]. The mechanical and electrical systems of the manipulator have strong uncertainty in practical applications (including non-modeled uncertainty and modeling parameter uncertainty [3]), strong coupling, large influence of control signal time delay error [4], large influence of control signal error [5,6], etc. These problems are the difficulties of current control research [7,8].
In response to the above problems, many scholars have conducted a lot of valuable research work. As a reference, a manipulator error compensation method based on deep belief network and error similarity was proposed by [9]. This method requires data reprocessing and uses particle swarm optimization (PSO) to optimize the network structure and parameters, which cannot theoretically guarantee the stability of the control process. An adaptive fuzzy trajectory tracking method was proposed by [10]. When this method is applied, the steady-state error of the step response is good, but the response speed is too slow, and there is a large time delay error when the sinusoidal input response is used. Otherwise, a new type of fuzzy control algorithm was proposed by [11]. In this method, a large peak error occurs when the manipulator moves to the point of commutation, and the jitter is obvious during the rotation process. The overall error is about 0.06 degrees, and the control effect is not very ideal. An adaptive fuzzy control method with asymptotic tracking performance was proposed by [12]. The step response time of this method is about 4 s, which cannot meet the control effect in practical applications. A control method based on the combination of dual encoders was proposed by [13]. Instead of compensating the actuator, it actually reduces the error generated by the encoder. This approach actually improves the parameter accuracy of the object model. A predictive control method based on feedback linearization to compensate the time delay error in the control process of the manipulator system was proposed by [14]. This method can restore the stability of the control system when the time delay of the control signal reaches about 1 s, but the control accuracy is not improved. In addition, during the control process, load, speed, motion direction, and ambient temperature all have an impact on the control accuracy. Using deep learning to achieve high-precision compensation, the time complexity is too high and the robustness is poor.
The sliding-mode control method has advantages in electromechanical robot control due to its own characteristics. The most important advantage is that sliding-mode control itself has a certain robustness [15], so the controller part of the multi-degree-of-freedom manipulator control system is suitable to use the method based on sliding-mode variable structure control plus error compensation. In recent years, sliding-mode control has been increasingly used as the core control method of various complex control systems [16,17]. The main reason is that sliding-mode control has advantages in stability proof, which facilitates stability analysis after adding intelligent control algorithms to it.
In the aspect of manipulator error compensation, the fuzzy control method is a suitable method [18,19]. The main advantage of a fuzzy system is that it can express the experience and knowledge of human understanding by mathematical formulas [20][21][22]. However, with the increasing requirements of control precision, the traditional simple fuzzy control system based on expert experience has been unable to meet the needs [23]. The real-time self-organizing fuzzy error compensation method based on control data can break through the shackles of expert experience to a large extent, make full use of various data in the control process, realize self-organization and self-adaptation of fuzzy rules, and make the system grow into the most suitable rules-governing objects [24,25]. However, the existing self-organizing control algorithms based on real-time data drive still fail to effectively solve the problems of huge amounts of information, increasing entropy, and high computing cost in the control process.
Through the above content, we can easily find out that some problems in trajectory tracking of the n-DOF manipulator systems that need to be addressed are summarized as follows.
(1) Most of the previous studies are aimed at suppressing or compensating for a certain type of interference [26][27][28][29]. However, as a multi-coupling complex system with fast state change, the manipulator system has multiple types of disturbances in the movement process, and the information entropy of real-time data in the control process is high, and most of the disturbances are time-varying random disturbances [30]. (2) While pursuing high precision, previous studies usually inevitably increase the complexity of the control algorithm and increase the time complexity of the control system calculation [31][32][33], thus making the response time to significantly longer. In the n-DOF manipulator control system, the state response speed is very fast and the controller response speed is very high. It is difficult for the existing research to be competent in this kind of control work. (3) In the existing research of the fuzzy self-organizing control framework, the initial rules and structure of fuzzy neural network have a great influence on the control effect [34,35]. The fuzzy neural network usually is self-organized in the control process after setting one or several good initial rules. However, in complex object control, it is usually impossible to set a suitable initial rule in advance, and the network needs to learn by itself from the first rule.
Regarding the aforementioned problems, a new self-organizing type-2 fuzzy feedback neural network compensation control framework is proposed to handle the trajectory tracking problem of mechanical n-DOF manipulator system with many types of uncertain disturbances. The main contributions of this paper are listed as follows: (1) An error compensation algorithm for a n-DOF manipulator is proposed, which compensates various types of error perturbations in various mechanisms by directly comparing input and output errors, and reduces the computational complexity of data-driven fuzzy algorithms when the amount of real-time data and information entropy are large. At the same time, the stability of the closed-loop control systems are proved by Lyapunov stability theory, and the simulation also confirms its effect. (2) A fuzzy network structure self-organization method and parameter self-adaptation method are proposed for the control object with high real-time requirements of electromechanical system. The time complexity and response time of this method can meet the requirements of real-time control. (3) A fuzzy rule self-organizing increase and decrease index function is proposed. This function does not need the overall information of the data set (such as the overall variance, mean, etc., of the data set) nor does it need to store the previous data information. It only needs the current data information at each moment to guide the self-organizing network to achieve reasonable rule increase and decrease. At the same time, this paper explains the reason why the function is suitable for fuzzy self-organizing network from the theoretical structure level of the function.
The remainder of the paper is organized as follows. The first section introduces the research background and research status of the multi-DOF manipulator and its error compensation method. The second section mainly introduces the kinematic modeling and dynamics modeling of the control object of the manipulator used in this paper. The third section mainly introduces the controller design method of the manipulator, including sliding mode variable structure control and self-organizing interval type-2 fuzzy neural network design. The fourth section analyzes the stability of the control framework, which shows that the system is stable and convergent. The fifth section shows the operation results of the control system in the simulation environment. The last section summarizes the research presented in this paper and gives an outlook on future research directions.
To show the mentioned symbols clearly, their corresponding explanations are presented in Nomenclature Section.

Kinematic Modeling of the Manipulator
Kinematics studies the position, velocity, and acceleration of a rigid body in motion without considering the forces that produce the motion. In this paper, the D-H method is used to model the object's n-DOF manipulator. The advantages of this method are that it involves few parameters, the model is consistent, and the coordinate system has 6 degrees of freedom. As shown in Figure 1, there are four parameters of the Denavit-Hartenberg (D-H) method [5,36]: a, α, d, and θ.
The kinematics problem of a multi-degree-of-freedom manipulator studies the mapping between the joint space and the end operation space of the manipulator. The forward kinematics problem refers to solving the motion parameters of the end of the manipulator given the input motion parameters of the joints. The inverse kinematics problem is the opposite, in which the motion parameters of the input joints are solved given the motion parameters of the end of the manipulator. In practical application, the user can only give the expected end motion trajectory and parameters of the manipulator, so the solution of inverse kinematics is more important. The kinematics model established in this paper is based on the right-handed coordinate system [37], as shown in Figure 1. To sum up, the forward and inverse kinematic equations of the n-DOF manipulator designed in this paper are shown in (1) and (3), respectively.

Dynamic Modeling of the Manipulator
The dynamic model of the n-DOF manipulator reflects the relationship between the motion of the manipulator and the force/torque it is subjected to. The three-degree-offreedom manipulator discussed in this paper is a free-moving manipulator; that is, the movement of the manipulator is not affected by environment obstacles. The position of this manipulator can be uniquely determined by the joint variable q, and each q is independent of each other. Therefore, we choose q = [q 1 , q 2 , . . . , q n ] as the generalized coordinate of the robot. In this paper, a three-DOF manipulator is modeled based on the second Lagrange equation.
For a three-DOF manipulator, as shown in Figure 1, the moment of inertia and inertia product in the three-dimensional direction are as follows [38]: Because link i dm = m i , a symmetric constant matrix can be obtained to describe the mass distribution of the rod: Since there are identities for any pair of indicator functions, The kinetic energy and potential energy of the entire robotic arm can be described as From (6) and (7), the Lagrange equation and its partial derivatives can be obtained as It follows from the Lagrange equation that Available: Writing (13) in matrix form, the dynamic model of the ideal n-DOF manipulator can be mathematically described as In the actual environment, there are measurable and immeasurable external disturbances such as friction. Therefore, the dynamic mathematical model of the multi-degree-offreedom manipulator is Among them, H(q) is the inertia matrix, C(q,q) is the centripetal force matrix and the Coriolis force matrix, G(q) is the gravity matrix, F(q,q,q) is the uncertainty item composed of friction, and λ is the uncertainty item caused by parameter uncertainty and external interference. The dynamics model of the three-degree-of-freedom manipulator discussed in this paper is as follows [39]: The three-degree-of-freedom manipulator system shown in (15) has the following characteristics [40]: (1) The inertia matrix H is a symmetric positive definite matrix and bounded, that is (2)Ḣ(q) − 2C(q,q) is an obliquely symmetric matrix, that is

Design of Robotic Arm Controller
In this paper, aiming at the manipulator object to be controlled, the sliding mode control is used as the coarse tracking controller, T2-SOFNNEC is used as the fine tracking controller, and the total control signal is the output control signal of the two controllers. The block diagram of the control system is shown in Figure 2.

Design of Sliding Mode Controller
In this paper, the manipulator shown in Section 2 (15) is used as the control object to design the sliding mode controller. The goal of robot trajectory tracking control is to require the actual joint angular displacement vector q = q 1 , q 2 , q 3 T to track the desired joint angular displacement vector q d = q d1 , q d2 , q d3 T as well as possible.
The tracking error is defined as The sliding surface is designed as Among them, s = s 1 , s 2 , s 3 T , Λ = diag(r 1 , r 2 , r 3 ) is a positive diagonal matrix, and r i > 0, i = 1, 2, 3 is the slope of the sliding mode surface. Since the uncertain items in the model cannot be known, the object model is simplified as The silding mode controller is designed as

Design of Self-Organizing Interval Type-2 Fuzzy Error Compensation Controller (T2-SOFNNEC)
The type-2 fuzzy system is further fuzzified on the basis of the type-1 fuzzy system.The single membership function in Type-1 fuzzy systems is extended to the upper membership function and the lower membership function. Its fuzzy logic system can be described by the following "IF-THEN" statement: where R (j) represents the jth rule; A is the antecedent interval type-2 fuzzy set; A i n represents the ith rule of the jth input; B j represents the consequent type-1 fuzzy set.
The most important difference between the type-2 fuzzy system and the type-1 fuzzy system is the membership function. The type-2 fuzzy is a three-dimensional surface, x is the input, u is the primary degree of membership, and µ A (x, u) is the secondary degree of membership. In the interval type-2 fuzzy set µ A (x, u) ≡ 1, this greatly reduces the amount of computation.
In order to make full use of the information related to the error in the control process, this paper uses the error, the error rate of change, the angle, and the angular velocity as the reference input of the error compensation algorithm, and the output is a compensation control signal. Obviously, the error compensation loop is a multi-input single-output structure. At the same time, the response ability of the proposed control algorithm should be ensure at the structural level. Therefore, this paper adopts the structure design compensation loop algorithm of IT2-FNN. The structure of the IT2-FNN is shown in Figure 3. The structure is mainly divided into five layers: (1) Input layer: The input layer is the transmission of the reference input, and each neuron represents an input dimension. (2) Membership layer: Each neuron in the membership layer represents the antecedent of a fuzzy rule, and each membership function takes the Gaussian function as the basis function, which can be expressed as Among them, c ij represents the center of the jth rule input in the ith dimension, σ ij represents the width of the jth rule input in the ith dimension, and ∆σ represents the half interval width of the membership function. (3) Rule layer: The value of each neuron in the rule layer is normalized by the corresponding membership degrees of all inputs, which can be expressed as (4) Consequent layer: The output of the consequent layer is calculated by multiplying the rule value by the weight: (5) Output layer: The function of the output layer is to downgrade the interval type-2 fuzzy system, and its output value is the product of the consequent value and the proportional coefficient: According to the above description, the biggest difference between the algorithm proposed in this paper and other type-2 fuzzy neural network structures is that although this paper adopts the interval type-2 fuzzy structure with a large amount of calculation, the depth and number of nodes of the fuzzy neural network are small, and there is no crossover between different rule antecedents. In fact, the algorithm proposed in this paper makes use of the stronger ability to deal with uncertain items of interval type-2 fuzzy sets and adopts a concise network structure so that the controller can achieve the best control effect with the fewest rules in the control process, thereby reducing the number of rules.
A large number of crossover operations are used to reduce the time complexity of the control algorithm. In order to realize the real-time self-organization in the control process, a rule increase and decrease function is proposed in this paper, and a rule increase and decrease selforganization algorithm is designed based on this function. The self-organizing algorithm and its nested structure in the IT2-FEC algorithm are shown in Figure 4. At the beginning of each control cycle, if there is no rule, a new rule is directly generated, and then the fuzzy network is calculated. If there is already a rule, we first update the network structure parameters according to the current input data and calculate the cost function and the rule increase or decrease index function. Then, according to the index function, it judges whether the rules need to be added or deleted. In order to ensure that the control algorithm always has a fast real-time response, each control cycle can perform rule addition or deletion at most once. After completing the self-organization of the rules, the control algorithm obtains the control signal through the fuzzy neural network calculation, and so on. The important components of the proposed control algorithm are described below.
(1) Estimation error In order to guide the network structure and parameter learning, a network structure risk assessment model is proposed, which is defined as [41]   where K is the current total number of rules, T is the current moment, and t g is the current cycle number. The estimation error is an evaluation of the risk of the network structure.
(2) Experience error The empirical error of the fuzzy network is defined as [41] Obviously, the empirical error index evaluates the compensation effect of the current network on the error on the one hand, and measures the adaptability of the network to new samples on the other hand.

(3) Rule increase and decrease function
In the self-organizing algorithm, when the controller increases the rules and when the rules are reduced is one of the most important steps and requires designing a reasonable and reliable rule increase and decrease function to help the controller determine whether to increase or reduce the rules. However, in the past, there was no paper to systematically design and evaluate the rules indicator function, and no paper explained the working mechanism of the rules to increase or decrease index functions.Therefore, a rule increase and decrease function is proposed in this paper. Otherwise, the working mechanism of the function during the control process is explained later.
Synthesizing the estimated error index and empirical error index of the network, the rule increase and decrease function is defined as The graph of this function is shown in Figure 5, from which we can see that the function value increases with the number of rules, while the partial derivative with time is different under different number of rules. It can be seen from Figure 6 that when the time is small, if the number of rules is large, the function value is too large, and the system adjusts by deleting rules; if the number of rules is small, the function value will be too small, and the system will adjust by adding rules. When the time is long, the function value tends to be stable, avoiding the system frequently adding and deleting rules and wasting computation.  (4) Rule growth stage If K d < α 1 , the network will add a new membership layer neuron to improve the generalization performance of the network. The parameters of the neurons newly added to the network are c new = x(t) where c new is the center of the new neuron, x(t) is the current input value, θ new is the width of the new neuron, and ω new is the weight of the new neuron. When the network structure is adjusted, its parameters are updated as where C(t), Θ(t), Ω(t) are the regular center matrix, width matrix and weight matrix of the network, respectively. It can be seen from (34) that every time a rule is added to the fuzzy set, it means that the data information of the new rule precedents appears in the discussion domain area with low probability. In this case, the fuzzy system covers the situation of the system information entropy being larger, which is the self-organization principle of the proposed algorithm.

(5) Rule deletion stage
In order to avoid the rule explosion caused by the continuous entropy increase of the system, the rule deletion stage plays an important role in the self-organizing algorithm. If k d > α2, the redundant neuron with the smallest activation function in the network will be deleted. The redundant parameter membership neuron will be set to: In order to improve the generalization performance of the fuzzy neural network, this paper uses the gradient descent method to optimize the structural parameters in each cycle. The optimization method is as follows: In the past, because the designed control object is a slow -changing object, the real-time requirements of the control algorithm are not high. Many various T2-FNN were proposed with a complicated self-organizing algorithm, which causes the algorithm to sacrifice more fastness under the requirements of the pursuit of control accuracy [42]. In fact, this is also a difficult point in the current high-level self-organized fuzzy algorithm. However, The above-mentioned self-organizing algorithm proposed in this article only increases or deletes the rules in a cycle cycle, and there is only one parameter correction, which not only greatly improves the real-time response capacity of the control algorithm but also obtains the requirements in organizing learning. It is obvious that the proposed algorithm can meet the complicated systems such as fastness and control accuracy of the robotic arm movement control system.

Stability Analysis of Sliding-Mode Control Loop
For (25), choose the Lyapunov function as Substituting (21) into (40) givesV Bringing the above formula into (25), we can obtain Then, Bringing (43) into (41) givesV Since the matrix P is positive definite, it can be obtained from (44) thaṫ It can be seen from (39) and (45) that the Lyapunov function V is positive definite and theV is negative define. According to Lyapunov's global stability theorem, the slidingmode control system designed in this paper is globally asymptotically stable at the equilibrium point.

Stability Analysis of Self-Organized Growth Stage
In the growth stage, a new fuzzy rule is generated which changes the structure of the self-organizing interval type-2 fuzzy error compensation network. According to (34), when the network adds a new rule, its prediction error is According to (27)-(30), the above formula can be extended as where f M+1 andf M+1 are both less than unity.
))) 2 According to the above analysis, when a new fuzzy rule is added to the network, its convergence speed is accelerated.

Stability Analysis of Self-Organizing Cut-Out Stage
When the network removes a rule, the prediction error of the network is According to (36), the above formula can be written as Therefore, the convergence rate of the network remains unchanged after pruning a fuzzy rule.
Based on the above analysis and discussion, the three-DOF manipulator trajectory tracking control algorithm based on self-organizing interval type-2 fuzzy error compensation proposed in this paper can ensure its stability. The self-organizing link can ensure the convergence and even speed up the convergence in the parameter update stage, the rule addition stage and the rule deletion stage.

Simulation
In this paper, the control simulation is performed on the MATLAB2021 platform. The parameters of the target manipulator are shown in Table 1. The manipulator is controlled by SMC, T1-SOFNNEC [41] and T2-SOFNNEC.

Simulation with Base Jitter
First of all, the effectiveness of the proposed control algorithm is verified in the case of only simulated base shaking. In the simulation, the pedestal jitter is expressed in the form of modeling uncertainty; that is, a time-varying random matrix of [0, 10] is added to the plant dynamics model. The control effects of SMC, T1-SOFNNEC and T2-SOFNNEC under the same conditions are compared as shown in Figure 7.
It can be seen from Figure 7 that all three control methods can realize the control of the manipulator. However, compared with the sliding-mode control, the self-organized fuzzy compensation control has a smaller error and higher control accuracy when the rotation direction of the manipulator joint changes. In addition, it can be seen from the figure that the self-organization method has a faster learning speed, and the control effect during the self-organization is almost the same as that of the sliding mode control.

Simulation with Base Jitter and Signal Interference
In this part of the simulation, signal interference is added to the base jitter. In the simulation process, the manifestation of the signal interference is that the time-varying random interference signal of [0, 10] is added on the basis of the output control signal of the controller. The simulation results are shown in Figure 8.
It can be seen from Figure 8 that under the action of the self-organized fuzzy error compensation algorithm, the error fluctuation range of the three joints of the manipulator is significantly reduced, and the control accuracy is significantly improved. At the same time, the robustness of the self-organizing fuzzy error compensation algorithm is also improved; that is, the error jitter in Figure 8 is significantly reduced. The RMSE of the three control framework in the presence of both jitter and control signal interference are shown in Table 2. It can be seen that the proposed T2-SOFNNEC has outstanding performance. This is due to the lifelong learning ability of the self-organizing method, which can adjust the controller structure according to the actual situation during the control process. At the same time, the interval type-2 fuzzy system itself has a strong ability to deal with uncertain terms, which has an advantage when multiple types of disturbances exist at the same time. The change in the rule number is shown in Figure 9. It can be seen that both algorithms are in the stage of adding rules at the beginning of the control, but the interval-type two algorithm gradually adjusts to only need one rule as the control time increases. However, the type-1 algorithm needs to be increased to six rules to compensate. In addition, the change of the rule index function in the control process of the proposed T2-SOFEC is shown in Figure 10, from which it can be seen that the proposed T2-SOFEC can regularly increase or decrease the index function during the control period between the thresholds. Therefore, the self-organizing interval type-2 fuzzy error compensation algorithm can achieve almost the same compensation effect with fewer rules than the self-organizing type-1 fuzzy error compensation algorithm, which makes it significantly reduce the time complexity in the control process, enhance the computational efficiency of the controller, and speed up the response time.
From this part of the simulation, it can be seen that the self-organizing fuzzy compensation algorithm has greater advantages than the previous algorithm and can compensate for errors to achieve better control effects in the presence of multiple disturbances. At the same time, under the same interference conditions, T2-SOFNNEC to T1-SOFNNEC uses fewer rules, which further reduces the storage and computation of the algorithm.

Simulation with Base Jitter, Signal Interference, and Time Delay
Based on the simulation in the previous part, this part adds the time delay interference of the control signal to simulate the delay of the control signal communication in the actual control process. This situation generally exists in the control process of bus communication.
In the simulation, the control delay shows that there is no control signal within one second when the control starts, and the control object receives the control signal output by the controller only after one second. When the system control signal delay is one second, the motion system error of the manipulator is shown in Figure 11. It can be seen from Figure 11 that the system is in the uncontrollable divergence phase in the previous second. After one second, the T2-SOFNNEC can quickly stabilize the system and restore to a stable and controllable state. In the presence of time-delay error, the interval type-2 fuzzy system can make use of its more powerful advantages in dealing with uncertainty and make up for the limitation of insufficient number of rules. From another point of view, interval type-2 fuzzy error compensation trades less rules for its stronger advantage in dealing with the structure of uncertain terms.

Coppeliasim-MATLAB Co-Simulation
In order to further verify the feasibility of the method, a dynamic simulation environment such as Figure 12 is established on the Coppeliasim platform, and the control method proposed in this paper is verified by the MATLAB-Coppeliasim co-simulation method. Coppeliasim is a multi-platform robot simulation software that covers robot modeling, programming, and simulation functions. It has the characteristics of close to the real real control environment and supports multiple programming languages.
The environment for co-simulation is the same as in the previous section, and the cosimulation results are shown in Figures 13 and 14. It can be seen from Figure 14 that the initial state of the robotic arm is not at the origin; that is, the initial angle of the three joints is not zero, which is very consistent with the actual control process. In the first 2 s, the robotic arm did not receive a signal from the control system, and the robotic arm was in a disabled state at this time. At a certain moment, the robotic arm receives the time-delayed control signal sent by the controller, and at this time the robotic arm eliminates the error and suppresses the interference. It can be seen from Figure 13 that the control algorithm can eliminate the error near zero in a short time. It can also be seen from Figure 14 that the control method has a practical effect in the co-simulation and has a good control effect. Last but not least, the longest response time of the simulation control signal is 2 ms, which can fully meet the needs of real-time control.

Conclusions
The difficulty of trajectory tracking control of n-DOF manipulators under strong interference lies in how to achieve stable and high-precision control of the system under the condition of various types of interference, various interference, and large amounts of realtime information data. The proposed T2-SOFNNEC in this paper can better compensate for various disturbances in the system and their resulting effects under the condition that the structure is more concise, the calculation cost of self-organization and parameter tuning is lower, and the rules used are fewer so that it satisfies the effect of both real-time control and accuracy. Meanwhile, this paper also proposes a more interpretable self-organizing rule increase and decrease index function. This function does not need to store the previous data information or the overall information of the data set and can realize the self-organization of the number of rules under the condition of low storage cost and computational cost. The proposed algorithm in this paper uses the fuzzy neural network structure to ensure the stability of the control system and uses the structure self-organization and parameter adaptive methods to ensure high-precision control. In summary, the proposed algorithm based on T2-SOFNNEC is suitable for high-precision motion control of n-DOF manipulators with strong interference in practical situations.
Author Contributions: Y.S. and T.Z. contributed to the manuscript conception and design of the study, performed the experiment and the data analysis, and wrote the manuscript; N.L. contributed to the manuscript preparation and important revisions to the paper. All authors have read and agreed to the published version of the manuscript.
Funding: This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
Institutional Review Board Statement: Not applicable.

Data Availability Statement:
The research data supporting this publication are provided within this paper.

Conflicts of Interest:
The authors declare no conflict of interest. Nomenclature c ij The ijth term of the centrifugal force matrix C Centrifugal force matrix c j The center of the jth rule e Control-error signal f j The upper output of the rule layer f j The lower output of the consequent layer f i Force of friction F Static and dynamic friction matrix g i The ith gravitational vector G Gravitational vector matrix h ij The ijth term of the inertial matrix H Inertial matrix k d The index function l i The length of the ith arm q d Expected angle q i The actual angle of the ith joint y The upper output of the consequent layer y The lower output of the consequent layer y The output of the FNN α The parameter of rule adaptive θ j The width of the jth rule λ The uncertain disturbance term caused by modeling error and external interferencē µ ij The value of Upper membership function µ ij The value of Lower membership function ρ 1 Estimation error ρ 2 Experience error τ Control moment