Adaptive Neural Tracking Control of Robotic Manipulators with Guaranteed NN Weight Convergence

Although adaptive control for robotic manipulators has been widely studied, most of them require the acceleration signals of the joints, which are usually difficult to measure directly. Although neural networks (NNs) have been used to approximate the unknown nonlinear dynamics in the robotic systems, the conventional adaptive laws for updating the NN weights cannot guarantee that the obtained NN weights converge to their ideal values, which could degrade the tracking control response. To address these two issues, a new adaptive algorithm with the extracted NN weights error is incorporated into adaptive control, where a novel leakage term is superimposed on the gradient method. By using the Lyapunov approach, the convergence of both the tracking error and the estimation error can be guaranteed simultaneously. In addition, two auxiliary functions are introduced to reformulate the robotic model for designing the adaptive law, and a filter operation is used to avoid measuring the acceleration signals. Comparisons to other well-recognized adaptive laws are given, and extensive simulations based on a 2-DOF SCARA robotic system are given to verify the effectiveness of the proposed control strategy.


Introduction
Robotic manipulators have been widely used to operate some special or repetitive tasks, and accurate modeling and control of robotic manipulators could promote their practical applications.However, nonlinear time-varying dynamics and coupling properties in the robotic systems usually make them a nonlinear multi-input-multi-output (MIMO) system, which is also taken as a preferable experimental platform for researchers to verify different modeling and control algorithms [1].Hence, those facts stimulate significant research interests from both industrial and academic communities to study the modeling and control of robotic manipulators.However, with the rapid and diverse applications of robotic technology, the working environments of robots are becoming more complex and even harsh such that the requirements for robust and flexible control strategies are becoming increasingly demanding, especially in the case with a changing environment, high-velocity motion, and varying load [2].
In the past decades, many studies on the control design for robotic manipulators have been reported in the literature, such as [1,[3][4][5][6][7].According to these studies, the tracking control system for robotic manipulators is usually composed of the following elements: a desired path specifies the ideal response of the controlled system; a numerical model describes the robot dynamics and its interactions with surroundings; a well-designed controller is used to generate appropriate control signals to drive the actuators to create the required torques.Then, in this way, the robotic manipulator can follow the desired trajectory.In order to use adaptive control techniques, the linear parameterization of a robot with nonlinear dynamics was derived [3] since the 1980s.Subsequently, computed torque-based adaptive control [1,3] was also proposed to guarantee the global error convergence.However, the modeling uncertainties existing in the robot manipulators are not considered.In [8][9][10], adaptive tracking control was proposed for trajectory tracking of robotic systems with uncertain kinematics and dynamics.However, these results are usually based on the assumption that uncertainties can be expressed in a linear parameterization form.In this sense, the aforementioned methods are mainly model-based control.Although they can obtain a satisfactory control performance in theory, the requirements for an accurate mathematical model limit their applications [11].Moreover, in most of the conventional adaptive control schemes, the joint acceleration signals are assumed to be available to design the controller [1].However, the measurement of joint accelerations is not practically sensible since it is generally sensitive to the external disturbance and noise [12].
To relax the requirements of the robot model knowledge, neural networks (NNs) [2,[13][14][15] have been used as an effective tool to approximate the nonlinear uncertainties.As it is shown in [16][17][18], any continuous functions can be uniformly approximated by a NN with an arbitrary degree of accuracy.In this framework, an adaptive law should be designed to update the unknown NN weights and to retain the system stability.There have been many adaptation methods for updating the NN weights, such as the gradient algorithm, least square (LS), σ-modification, and projection algorithm [19].Although these methods can achieve fair tracking performance, only a few of them focus on the learning performance and convergence of the estimated NN weights.The authors of [20] pointed out that with inaccurate and slow convergence of the estimated parameters, the effects induced by the unknown robot dynamics cannot be suppressed in transient time, which may degrade the control performance.To guarantee the parameter convergence, the authors of [21,22] proposed a composite learning for parameter estimation.Recently, the authors of [6,23] proposed a novel adaptive law, in which both the tracking control and parameter estimation convergence can be guaranteed simultaneously, while the linearly parameterized form should be imposed again.
In this paper, an adaptive neural control based on a radial basis function neural network (RBFNN) will be proposed for robotic manipulators to achieve guaranteed tracking control and estimation.Firstly, since the measurement of joint accelerations is sensitive to the external noise, we aim to avoid using the acceleration signals directly by reformulating the robotic model.Hence, two auxiliary variables are first designed to reconstruct the robotic model.Then, a low-pass filter is applied to the reformulated model such that the joint acceleration can be avoided in the design of adaptive laws.Moreover, to relax the requirement of system knowledge, a RBFNN is employed to approximate the lumped unknown nonlinear dynamics.To retain the control and estimation convergence, the adaptive parameter estimation proposed in [6] is further tailored and incorporated into an adaptive neural control design.In this new framework, both the tracking error and the NN weight estimation error are used in the adaptive law, and thus, the convergence of the tracking error and weight estimation error can be guaranteed simultaneously.In addition, comparisons to other well-recognized adaptive laws are presented in this paper, which show that better estimation performance with the proposed algorithm can be obtained.Numerical simulations based on a 2-DOF SCARA robot are provided to show the efficiency of the proposed methods.
This paper is organized as follows: The mathematical model of robotic manipulators and a brief introduction of the RBFNN are given in Section 2. An adaptive neural control design based on RBFNN and the stability and robustness analysis are presented in Section 3. Section 4 provides comparisons to other adaptive algorithms.Section 5 shows simulation results based on a 2-DOF SCARA robot.Conclusions of this paper are given in Section 6.

Problem Formulation and Preliminaries
2.1.Robot Manipulator Dynamics.In this paper, we consider robotic manipulators where a set of n rigid bodies are connected in series with the final arm being fixed to the ground.The model of such an n-degrees of freedom (DOF) nonlinear robot manipulator can be expressed as M q q + C q, q q + G q = τ, 1 where q, q, and q ∈ ℝ n are the robot joint position, velocity, and acceleration, respectively; n is the number of DOF; τ ∈ ℝ n is the control torque vector; M q ∈ ℝ n×n is the inertia matrix; C q, q ∈ ℝ n×n denotes the lumped Coriolis/ centripetal torque, viscous and nonlinear damping; and G q ∈ R n represents the gravitational effect.The problem to be addressed in this paper is to design a proper adaptive control such that the output q of robotic manipulator (1) can track a given reference q d , while the convergence of the adaptive laws can be achieved.
According to [1], the following two properties used in the subsequent control design should be presented: Property 1.The matrix M q ∈ R n×n is symmetric and positive definite.
Property 2. The matrix M q − 2C q, q is skew-symmetric; that is, x T M q − 2C q, q , x = 0, ∀x, q, q ∈ ℝ n 2 2.2.RBF Neural Network.Neural networks (NNs) have been widely used as a function approximator for unknown nonlinearities due to their elegant approximation abilities [24].
In this paper, a linearly parameterized radial basis function NN (RBFNN) is employed to approximate a continuous function f X : ℝ l ⟶ ℝ over a compact set Ω z as where z = z 1 , z 2 , ⋯, z l T ∈ Ω z ∈ ℝ l is the input vector of NN; N is the number of the NN nodes; ε ∈ ℝ N is the approximation error; and W * = w * 1 , w * 2 , ⋯, w * N T ∈ ℝ N is the ideal NN weights.As shown in [25], the ideal NN weights W * are defined as where ψ z = ψ 1 z , ψ 2 z , ⋯, ψ N z T is a nonlinear function vector, the components of which can be presented as where ς i = ς i1 , ς i2 , ⋯, ς il T ∈ ℝ l denotes the center of the ith basis function and ρ i is the ith variance.
The following assumptions and definition will be used: Assumption 1.The ideal NN weights W * are bounded such that W * ≤ W holds for a positive constant W, and the approximation error is bounded by ε ≤ ε * for a positive constant ε * .
Assumption 2. The reference signal q d and their derivative q d are smooth and bounded.

Adaptive Control Design Based on RBFNN
In this section, a constructive control design will be presented for a robotic manipulator (1), which has two salient features in comparison to the existing results.On the one hand, a new filter operation and proper mathematical developments are introduced to avoid measuring the joint accelerations directly, which is preferable for applications.On the other hand, a novel adaptive law is developed to update the NN weights online, which can guarantee the exponential convergence of the estimated weights.
3.1.Adaptive Controller Design.As shown in (1), the studied robotic manipulator is modeled using Lagrange-Euler functions, which usually contain the acceleration signal q.This implies that the acceleration signal needs to be available to design the controller.However, the direct measurement or estimation of acceleration signals is generally difficult due to the induced noise.To avoid using the acceleration signal directly, we first define an error variable r q, q as r q, q = e + Λe, 6 where e = q − q d and q d is the desired reference of the robot joint; Λ ∈ ℝ n×n is a positive diagonal matrix.
According to Assumption 2 and the defined variable r in (6), we have e = r − Λe = q − q d and e = r − Λe = q − q d .Then, the velocity and acceleration can be expressed as Substituting ( 6) and ( 7) into the controlled system (1), we have However, the design of control law τ could not be achieved based on (8) because of the nonlinear dynamics induced by M q Λe − q d + C q, q Λe − q d − G q .To accommodate this issue, an RBFNN is employed to approximate these lumped nonlinearities.Thus, system (8) can be described as where W * is the ideal NN weights; ψ z is the basis function with the input vector; z = q T , q T , q T d , q T d , q T d T ; and ε is the approximation error, which are defined in ( 3), (4), and ( 5).
In order to achieve the tracking control, an adaptive neural controller is designed as where K > 0 is a feedback gain matrix and Ŵ is the estimated NN weights of the unknown ideal weights W * , which will be online updated based on the designed adaptive law.Substituting ( 10) into ( 8), the closed-loop control error dynamics with RBFNN can be derived as where W = W * − Ŵ is the NN weight error.The problem to be solved is to design an adaptive law, which can obtain the estimated RBFNN weights Ŵ with guaranteed convergence for implementing the proposed control (10).
Remark 1.Most of existing adaptive laws for updating Ŵ are designed based on the gradient methods on the tracking error dynamics (11), e.g., [1,18,26,27] and the references therein.However, the use of robust leakage terms (σ-modification and projection algorithm) in this framework makes it difficult to retain the convergence of the NN weights.Hence, we will develop a new adaptive algorithm, which guarantees the convergence of the estimated NN weights Ŵ, while retaining the robustness.
From the above analysis, the acceleration q is embedded in the error signal r.To make adaptive law independent of the joint acceleration, two auxiliaries are further defined as Hence, (8) can be further written as It should be noted that the introduced term F 1 z contains the information of the joint acceleration q; hence, it cannot be used directly to design the adaptive law.In this case, inspired by [6], a filter operation is applied on system (9) to eliminate the acceleration.We define the filtered variables F 1f , F 2f , ψ f , and τ f with respect to F 1 , F 2 , ψ, and τ as where k > 0 is a filtered parameter.
Then, applying a filter 1/ ks + 1 on system (9), we have where s is the Laplace operator.From the first equation of ( 14), we can obtain F 1f = F 1 − F 1f /k such that (15) can be rewritten as where ε f is the filtered version of NN error ε by using the same filter operation as (14).Since the NN approximation error ε is bounded, its filtered version ε f is also bounded, i.e., ε f ≤ ε * .Through the above filter operations, only the variable F 1 with q, q and its filtered version F 1f is used in (16), where the acceleration q can be avoided.Thus, we can design an adaptive law based on ( 16) rather than the original tracking error (11).Hence, the rest of this section will present a novel adaptive algorithm based on (16), where a new leakage term containing the estimation error is introduced to achieve a guaranteed convergence of NN weights.By using the above-filtered variables F 1f , F 2f , ψ f , and τ f , two auxiliary matrices P ∈ ℝ N×N and Q ∈ ℝ N×n are defined as where l > 0 is a positive constant serving as the forgetting factor to retain the boundedness of P and Q.
The solution of matrix equation ( 17) is derived as Then, by using the obtained auxiliary matrices P and Q, we further define an auxiliary matrix H as It is shown in (19) that the matrix H can be online calculated based on the auxiliary matrices P and Q by using the filtered variables F 1f , F 2f , ψ f , and τ f and the estimated NN weights Ŵ.Hence, it can be used to design the following adaptive law.Now, an adaptive law with H being a new leakage term is given as where Γ ∈ ℝ N×N is a positive diagonal matrix and κ > 0 is a positive parameter.
To show the merit of the proposed adaptive law (20), we fist present the following fact.Lemma 1.The defined variable in (19) is equivalent to where φ t = − t 0 e −l t−τ ψ f r ε T f dr is a bounded residual error and W = W * − Ŵ is the estimation error.Proof 1.By substituting (16) into the second equation of ( 18), we can verify the fact Q = PW * + φ.Then, by substituting it into (19), the fact shown in ( 21) can be verified.Moreover, since the RBFNN regressor ψ is Lipschitz and the error ε is bounded, the residual error φ is also bounded, i.e., φ ≤ ϖ for a positive constant ϖ > 0.
As shown in (21), one may find that the explicit formulation of the unknown NN weight estimation error W is embedded in the matrix H, which can be online calculated based on the known or measurable robotic dynamics.Specially, when the approximation error is null, i.e., ε = 0, we can verify that H = −PW holds.In this case, H contains the NN weight error to be minimized.Hence, the use of H in the adaptive law can help to guarantee the convergence of the RBFNN weights Ŵ to the ideal unknown value W * .This will be proved in the next subsection.
Remark 2. In the proposed adaptive law (20), the first term ψr T is the conventional gradient term, and the second part κH is the newly introduced leakage term containing the weight error W. With the help of this new leakage term, a 4 Complexity quadratic term of the NN weights error W will be included in the derivative of the Lyapunov function as shown in the next subsection, such that the convergence of tracking error and NN weights error can be guaranteed simultaneously.This is clearly different to other existing adaptive laws.
Remark 3.Although some solutions have been reported to avoid using the acceleration signals q (e.g., [1]) or to use differentiators in the design of adaptive laws, we provide an alternative feasible method by reformulating the robotic model with new variables F 1 and F 2 .In this case, the acceleration signals are not required and the extra differentiator is not required, which leads to a simplified control implementation.

Stability and Convergence
Analysis.Before proving the convergence of the closed-loop system, we first show the positive definiteness of the introduced matrix P. For the sake of simplicity, we define λ min ⋅ and λ max ⋅ as the minimum and maximum eigenvalues, respectively, of the corresponding matrices.Then, we can present the following lemma.
The stability of the closed-loop system and the convergence of the tracking and estimation errors can be given as Theorem 1.Consider robotic manipulator (1) with controller (10) and adaptive law (20), if the RBFNN regressor ψ z is PE (i.e.λ min P t > σ > 0, t > 0); then, all signals in the closedloop system are bounded, and the error variable r, tracking error e, and the NN weights error W all converge to a small compact set around zero.Proof 3. Select a Lyapunov function as Based on Property 2 shown in Section 2, we can calculate its derivative V 1 along ( 11) and (20) as By using Young's inequality [30] for the last two terms in (23), we have where α 1 = min 2λ min K − 1/η 1 /λ max M , 2κσ − κ/2 η 2 /λ max Γ −1 , and β 1 = η 1 /2 ε * 2 + κη 2 /2 ϖ 2 are all positive constants since the constants η 1 and η 2 can be designed to satisfy the condition η 1 ≥ 1/2λ min K and η 2 > 0. Consequently, based on the Lyapunov stability theory, we can obtain that the error variable r and the NN weight error W are bounded.From the (6), we can verify that the tracking error e is also bounded.In addition, since the regressor ψ is bounded, we know that the control signal τ is also bounded based on (10).Furthermore, the solution of inequality ( 24) satisfies V t ≤ V 0 − β 1 /α 1 e −α 1 t + β 1 /α 1 .Hence, the control error r and the NN weights error W will converge to a compact set around zero defined by Ω w ≔ L W, r | W ≤ V 0 e −α 1 t + β 1 /α 1 /λ min P , r ≤ V 0 e −α 1 t + β 1 /α 1 , whose size depends on the RBFNN approximation error ε * , the persistent excitation level σ, and the learning gain Γ.Consequently, based on (6), the tracking error e will also converge to a compact set around zero.Hence, the system output q will converge to a neighborhood around the desired reference q d .This completes the proof.
Remark 4. The online validation of the PE condition through directly checking the regressor vector ψ z is generally difficult in particular for nonlinear robotic systems, because the definition of PE is indeed an integral within a shifting time interval.However, according to Lemma 2, the PE condition of the RBFNN regressor ψ z is equivalent to the positive definiteness of the introduced auxiliary matrix P, such that testing for the positive definiteness of P (e.g., calculate the minimum eigenvalue or rank of matrix P) allows to numerically verify the PE condition online.Thus, this new adaptive algorithm provides a feasible method for online verification of PE condition.Specifically, as shown in [25], the RBFNNbased control system with periodical reference signal allows the PE condition partially fulfilled.

Robust Analysis against Exogenous Disturbance.
In this subsection, we further study the robustness of the proposed adaptive control methods for robotic systems with bounded exogenous disturbance.For this purpose, the robotic manipulator (1) can be represented as M q q + C q, q q + G q = τ + d t , 25 where d t denotes a bounded exogenous disturbance, i.e., d t ≤ ƛ, ƛ > 0.

Complexity
The proposed control (10) is used, while the tracking error (11) is modified as M q r + C q, q r + Kr = W T ψ + ε + d 26 Moreover, to analyze the convergence of the proposed adaptive law (20), ( 9) can be replaced by In applying the same filter operations given in ( 14) on ( 27) and defining the same auxiliary matrices Ρ and Q in (17), then the matrix H in (19) used to drive the adaptive law ( 20) can be given as where is also bounded.Hence, we know that the residual term Ψ is bounded because of the boundedness of ψ f , ε f , and d f ; i.e., Ψ ≤ ϑ holds for a constant ϑ > 0.
Then, we have the following corollary.

Corollary 1.
Considering the robotic manipulator (25) subject to bounded disturbances, the controller (10) and adaptive law (20) are adopted.If the RBFNN regressor ψ z is PE (i.e., λ min P t > σ > 0, t > 0), the closed-loop system is stable.Moreover, the error variable r and the NN weight error W are bounded; converge to a compact set around zero, of which the bound depends on the disturbance bound.
Proof 4. The proof of Corollary 1 can be carried out by selecting the same Lyapunov function V 1 = 1/2 r T Mr + 1/2 tr W T Γ −1 W and calculating its derivative along (26) and ( 20) with (28).The detailed proof is similar to that of Theorem 1 and thus will not be repeated again.The main difference is that the disturbance is involved in the error dynamics (26) and also the residual error Ψ in (28).Hence, the ultimate bound determined by the constant β 1 in the Lyapunov function (24) will be affected by the bound of disturbance ƛ, apart from the NN approximation error ε * .
Remark 5.In practical control implementation, there may be also measurement noise and modeling uncertainties, which could be lumped into the NN approximation error ε * .Hence, ( 27) can be further replaced as where Δ = ε + d ∈ ℝ n denotes the lumped disturbance containing the NN approximation error ε and the effect of disturbance and measurement noise d.Then, the robustness of the proposed control and adaptive law can be analyzed similarly.The only difference is that the residual term Ψ in ( 28) is modified as Ψ t = t 0 e −l t−τ ψ f r Δ T f dr, where Δ f is the filtered version of Δ given by kΔ f + Δ f = Δ, Δ f 0 = 0.

Comparison to Other Adaptive Laws
To show the advantages and superior performance of this new adaptation over classical methods, we will compare it with a widely recognized gradient-based adaptive algorithm and the σ-modification method [31].
4.1.Gradient Method.The classical gradient-based adaptive law is solely driven by the control error r, which is given by Considering the fact W = − Ŵ, the estimation error of ( 30) is derived as we have its derivative along ( 11) and ( 31) as

32
where α 2 = λ min K − 1/2η 3 and β 2 = η 3 /2 ε * 2 , with η 3 > 1/2λ min K .Based on the Lyapunov stability theory, the control error r will converge to a compact set around zero.However, the convergence of the NN weights error W cannot be guaranteed, especially in the absence of PE condition.More specifically, when the system is subject to disturbance and uncertainties, the online obtained NN weights Ŵ may even be unbounded; i.e., it could lead to parameter drift and bursting phenomena [32].Hence, the convergence and robustness of the gradient algorithm (30) are questionable.[31].To guarantee the boundedness of the estimated NN weights, a leakage term σ Ŵ is added in the gradient algorithm, leading to the following σmodification method:

σ-Modification
Then, its estimation error dynamics can be given as We select the Lyapunov function as V σ = 1/2 r T Mr + 1/2 tr W T Γ −1 W , and its derivative along ( 11) and (34) can be derived as

35
where the constants 2 are positive by setting the design parameters η 4 > 1/2λ min K and η 5 > 1/2 .Based on the Lyapunov theory, we can obtain that the error variable r and the estimation error W are all bounded; i.e., the σ-modification method has guaranteed robustness.
In fact, compared with (31), it is shown in (34) that a forgetting factor term −κΓW is involved, such that the estimation error dynamics in (34) are bounded-inputbounded-output (BIBO) stable.Hence, the σ-modification can guarantee the boundedness of the control error r and the NN weight error W.However, since a pure damping term κ Ŵ is added to (33), the estimated NN weights Ŵ will stay in a neighborhood of a preselected value only but will not converge to their true values, because the upper bound W of the unknown NN weights W * is involved in the error bound β 3 in (35) and thus β 3 ≠ 0 even when the NN weight error and tracking error are null (i.e., ε * = r = 0).This can also be observed from the estimation error (34) represented as a transfer function as W = 1/ s + κΓ −Γψr T + κΓW * .4.3.Proposed Adaptive Law.With the help of the new leakage term κH, the NN weight error of the proposed adaptive law (20) can be given as As shown in (36), a forgetting factor −κΓPW stemming from the new leakage term κH is involved to retain its BIBO stability, and thus the boundedness of the NN weight error can be guaranteed.In addition, the estimation error dynamics (36) can be further rewritten as a transfer function as W = 1/ s + κΓP −Γψr T + κΓφ .Consequently, the bound of the NN weight error mainly depends on the residual error φ, which is a function of the NN approximation error ε as shown in (21).Consequently, the NN weight error will converge to zero as long as the control error and the NN approximation error are null (i.e., ε * = r = 0).This implies that the proposed adaptive law achieves better estimation performance compared with the other two adaptive laws, while retaining the same robustness as the σ-modification method.Remark 6.Compared to the conventional adaptive laws, the proposed adaptive law (20) uses filter operations ( 14) to extract the information of the NN weight error.In this framework, the filter coefficient k (14) should be set small to retain fast convergence.On the other hand, the constant k defines the bandwidth of the filter 1/ ks + 1 .Hence, the robustness and transient convergence should be managed when we set this constant.The forgetting factor l in (17) aims to guarantee the boundness of the auxiliary matrices P and Q; thus, l cannot be set too large.The constant κ in (20) affects the convergence of the NN weight error W; in general, κ should not be set too large since the residual term κΓφ may be large with a large κ.Finally, the selection of the learning gain Γ in (20), the diagonal matrix Λ in (6), and the feedback gain matrix K in (10) should be considered the tradeoff between the convergence of the tracking error and the robustness of the closed-loop system.

Simulation
In this paper, a 2-DOF SCARA robot designed in our lab is used to show the validity of the proposed control schemes.The diagram of the 2-DOF SCARA robot can be found in Figure 1, and its schematic diagram can be found in Figure 2.
According to [33], the kinetic energy E and the potential energy Z of the 2-DOF SCARA robot can be given as

Complexity
Then, the dynamics of the 2-DOF SCARA robot are described by In the simulations, the initial joint position q 0 and the initial velocity q 0 are set as q 0 = 0 2 0 3 T and q 0 = 0 0 T , respectively.The desired trajectories to be tracked are q d1 = 0 5 sin t and q d2 = cos 0 5t .The number of the NN node is set as N = 16, and the initial RBFNN weights are set as Ŵ 0 = zeros 16, 2 .The center ς i in ( 5) is uniformly distributed in 1, −1 , and the variance ρ 2 i is set as 50.The other parameters used in the control and the adaptive law are Λ = 5I 2 , k = 0 0001, l = 90, κ = 1, K = 5I 2 , and Γ = 100I 16 .
Figures 3-6 show the tracking performance of robot joint position and velocity.As it is shown in these simulation results, the adaptive control (10) with both the proposed adaptive law (20) and σ-modification method can achieve satisfactory tracking performance.However, as it can be seen from Figures 4 and 6, smaller steady-state errors can be achieved by using the proposed adaptive law (20) in comparison with that of the σ-modification method.This result validates the analysis in Theorem 1 and the discussions in Section 3; that is, a quadratic term of the NN weights error W induced by the proposed leakage term H is incorporated into the derivative of the Lyapunov function, such that the convergence of both the tracking error and the NN weights error can be guaranteed simultaneously.It can be also observed from Figure 7 that the NN weights Ŵ updated by the proposed adaptive law (20) can achieve convergence after a short transient stage,    Moreover, to test the robustness of the proposed method, the bounded exogenous disturbance signals d t = 0 5 sin 16t e −0 2t + 0 03 ; cos 8t e −0 4t + 0 02 are added in the measured position signals of joint 1 and joint 2, respectively.Simulation results given in Figures 8 and 9 show that the proposed adaptive neural control is robust against the modeling uncertainties and disturbances as that achieved by using the σ-modification scheme.From these simulation results, one can claim that the proposed adaptive neural control for robot manipulators can achieve better convergence performance and comparative robustness against the external disturbances.

Conclusions
In this paper, we introduce a new adaptive algorithm to achieve the guaranteed estimation convergence and improve the tracking performance of adaptive neural control for robot manipulators.Two auxiliary functions are first introduced to reconstruct the robotic model such that the acceleration signals are not used in the control implementation and online learning.The unknown nonlinear system dynamics can be effectively approximated by using a RBFNN, where a novel leakage term is superimposed on the gradient-based adaptive law, which is used to online update the RBFNN weights.Through rigorous stability analysis based on the Lyapunov theory, it is proved that the proposed adaptive control could guarantee the convergence of both the tracking error and the NN weight error simultaneously.Therefore, improved tracking control performance can be achieved in comparison to other adaptive methods.In addition, it is shown that even in the presence of bounded disturbances, the proposed algorithm can retain the robustness.Simulations based on a 2-DOF SCARA robotic manipulator are given to validate the efficiency of the proposed adaptive neural control.A future work will focus on relaxing the required PE condition for the adaptive control design.

Figure 3 :
Figure 3: Tracking performance of robotic joint position q 1 and q 2 .

Figure 4 :
Figure 4: Tracking error of robotic joint position q 1 and q 2 .

Figure 8 :Figure 5 :
Figure 8: Tracking performance of the proposed adaptive law with bounded exogenous disturbances.

Figure 6 :
Figure 6: Tracking error of robotic joint velocity q 1 and q 2 .

Figure 7 :
Figure 7: The online obtained weight of RBFNN.

Figure 9 :
Figure 9: Tracking error of the proposed adaptive law with bounded exogenous disturbances.