Nussbaum Function Based Approach for Tracking Control of Robot Manipulators

Abstract

This paper introduces a novel Nussbaum function-based PID control for robotic manipulators. The integration of the Nussbaum function into the PID framework provides a solution with a simple structure that effectively tackles the challenge of unknown control directions. Stability is achieved through a combination of neural network-based estimation and Lyapunov analysis, facilitating automatic gain adjustment without the need for system dynamics. Our approach offers a gain determination with minimum parameter requirements, significantly reducing the complexity and enhancing the efficiency of robotic manipulator control. The paper guarantees that all signals within the closed-loop system remain bounded. Lastly, numerical simulations validate the theoretical framework, confirming the effectiveness of the proposed control strategy in enhancing robotic manipulator control.

I Introduction

The rapid expansion of robotic systems across various industries has driven the development of advanced control mechanisms to enhance their functionality and adaptability. Despite the progress facilitated by rigorous mathematical frameworks, the nonlinearities and uncertainties intrinsic to robotic operations continue to pose significant challenges. This reality highlights the need for advanced control strategies like fuzzy logic [1], Kalman filter [2], iterative learning control [3], and actor–critic learning [4] to effectively counteract these complexities. However, the inherent complexity of these methods often complicates their application in real-world scenarios, underscoring the ongoing research imperative to refine and simplify control solutions for robotic manipulators.

In the field of control for robotic manipulators, PID control is recognized for its intuitive design and simplicity, positioning it as a pillar for both theoretical exploration and practical application across numerous real-world systems [5, 6]. It addresses the widespread issue of complexity in robotic control strategies. However, traditional PID controls face limitations, particularly in weight updating, stability assurance, and the need for extensive parameter tuning in dynamic environments [7]. Attempts to refine PID control, such as employing optimization techniques for gain tuning, aim to enhance its adaptability [8, 9, 10]. However, these enhancements can introduce new complexities, somewhat undermining the original appeal of PID control’s simplicity.

Understanding the direction of control gain is crucial in robotic system control, as incorrect application can destabilize the system instead of guiding it to the desired state. While various methods exist to address the unknown control direction, such as logic-based switching [11], extreme seeking [12], and nonlinear PI control [13], the Nussbaum function approach [14] is the most studied method [15, 16, 17, 18]. Implementing Nussbaum-based strategies in robotic control is effective; however, their integration into the comprehensive control framework can complicate the overall system. This scenario highlights the importance of developing more accessible strategies that utilize the Nussbaum function’s advantages while avoiding excessive complexity in the control framework.

Given the complexities and challenges previously discussed, our research is driven by the need to develop a control strategy that is simple and low in complexity, effectively managing the critical aspect of unknown control gain direction in robotic manipulators. Acknowledging the limitations of traditional PID controls and the complexities of Nussbaum-based methods, we introduce a novel approach that combines the simple structure of PID control with the capability of the Nussbaum function. This method ensures stability and enhances performance, offering a control mechanism that is inexpensive in online computational demands, and effectively bridging the mentioned gap in robotic control strategies.

In this work, we present a control strategy for robotic manipulators, tackling unknown dynamics and control directions through an adaptive Nussbaum-based control. Our approach, rooted in the PID control structure, simplifies yet enhances the control scheme’s effectiveness. Different from conventional PID methods, we ensure the closed-loop system’s stability through direct Lyapunov analysis. Additionally, our strategy features automatic gain adjustment and utilizes linking parameters, significantly reducing the number of tuning requirements, and thereby enhancing the system’s efficiency and response without complicating the control framework. The main contribution can be summarized as follows:

• We integrate the Nussbaum function into the PID structure for the control of robotic manipulators. The controller is simple, yet it guarantees stability without requiring knowledge of robot dynamics or control direction.

• By employing linked PID terms and establishing adaptive laws, the controller automates gain determination with minimal parameter requirements, thereby facilitating its application in real-time scenarios.

The rest of the paper is structured as follows: Section II delves into problem formulation and preliminaries, setting the theoretical groundwork. Section III introduces the Nussbaum Function-based control design and its stability analysis. Numerical simulations demonstrating the control strategy’s effectiveness are presented in Section IV. The paper concludes with Section V, summarizing the key findings and implications of this research.

II Problem Formulation And Priliminaries

II-A Problem Formulation

Consider the robotic system described in the joint space $q(t)\in\mathbb{R}^{n\times 1}$ as [19, 20],

$$M(q(t))\ddot{q}(t)+C(q(t),\dot{q}(t))\dot{q}(t)+G(t)=\tau(t)\quad$$ (1)

where $M(q(t))\in\mathbb{R}^{n\times n}$ denotes a positive definite inertial matrix; $C(q(t),\dot{q}(t))\in\mathbb{R}^{n\times n}$ represents the Coriolis and centrifugal matrix; $G(t)\in\mathbb{R}^{n\times 1}$ denotes the gravitational force vector. $\tau(t)\in\mathbb{R}^{n\times 1}$ represents the actual torque applied to the robotic system. Furthermore, in practical systems, the actual torque is affected by unknown actuator dynamics and can be expressed as [21],

$$\tau(t)=\kappa u(t),\quad$$ (2)

where the nonlinear matrix $\kappa\in\mathbb{R}^{n\times n}$ is an unknown control direction, and $u(t)$ represents the control input.
The following properties of robot manipulators dynamics (1) are required for control analysis [22],

• P 1: $M(q(t))$ is symmetric positive definite.

• P 2: $\frac{1}{2}(C(q(t),\dot{q}(t)))-\dot{M}(q(t)))$ is skew-symmetric.

• P 3: There exist unknown positive constants $m_{-}$, $m^{+}$, $c$ and $g$ such that $m_{-}\leq M(q(t))\leq m^{+}$, $\|C(q(t),\dot{q}(t))\|\leq c\|\dot{q}(t)\|,$ and $\|G(q(t))\|\leq g<\infty.$

The goal of this paper is to design a control law for a robot dynamics (1) such that

• G 1: the closed-loop system is stable and all the signals remain uniformly bounded.

• G 2: the joint position signal $q(t)$ closely tracks a specified desired trajectory $q_{d}(t)$, so that the limit $\lim_{t\to\infty}|q_{i}(t)-q_{di}(t)|=\eta_{i}$, for $i=1,2,...n$, and $\eta_{i}$ being a small positive constant.

To this end, the following assumptions are needed.

Assumption 1

We consider each link of the manipulator as a slender rod with uniform mass distribution, ignoring additional components’ mass and inertia.

Assumption 2

The desired trajectory $q_{d}(t)$, along with its first $\dot{q}_{d}(t)$ and second derivatives $\ddot{q}_{d}(t)$, are assumed to be smooth, known functions of time and are bounded. Additionally, the robot position vector $q$ is assumed to be accessible for the purposes of control design.

Remark 1

The dynamics of the system, as described by equation (1), are fully unknown, which presents a significant challenge in the control design. This paper addresses the inherent difficulty of unknown control direction by employing a Nussbaum-type function. The Nussbaum approach is particularly adept at handling systems with unknown control coefficients, adapting to the control direction without prior knowledge of the system’s dynamics.

II-B Priliminaries

To approximate the system uncertainties including the manipulator’s unknown dynamics, in its continuous movement, we utilize linear-in-parameter approximators with an a priori-defined basis function vector and a yet-to-be-learned unknown parameter vector [23]. To this purpose, a generic function approximator $F(x)$ is denoted by $F(x,\psi)=\psi^{\top}\phi(x),$ where $\psi\in\mathbb{R}^{n_{p}}$ is the unknown parameter vector of dimension $n_{p}$ and $\phi(x)\in\mathbb{R}^{n_{p}}$ is the user-defined known basis function vector with ${n_{p}}$ being the NN input dimension, and $x\in\mathbb{R}^{m}$ is the neural network input vector with $m$ being the neural network input dimension. In this work, we have used the radial basis function (RBF) given by $\tilde{\phi}(x)=e^{-0.5(x-c)^{\top}B^{-1}(x-c)},$ where $c\in\mathbb{R}^{n}$ is the center and $B\in\mathbb{R}^{n\times n}$ the covariance matrix of the RBF.

Define the position error, $e(t)\in\mathbb{R}^{n}$, as $e(t)=q_{d}(t)-q(t)$ and the velocity error, $\dot{e}(t)\in\mathbb{R}^{n}$, as $\dot{e}(t)=\dot{q}_{d}(t)-\dot{q}(t)$. To further our analysis, we define the generalized intermediate variable $\Psi(t)$ as follows:

$$\Psi(t)=2\gamma e(t)+\gamma^{2}\int_{0}^{t}e(\rho)\,d\rho+\frac{d}{dt}e(t),$$ (3)

where $\gamma>0$. This formulation allows us to address the tracking error by stabilizing $\Psi(t)$ using the following lemma.

Lemma 1

[24] Given the intermediate variable $\Psi(t)$ as defined in (3), if $\Psi(t)\to 0$ as $t\to\infty$, then the tracking errors $e(t)$ and $\dot{e}(t)$, and their integrals are bounded and converge to zero over time.

Definition 1 (Nussbaum Function)

A Nussbaum function $N(\zeta)$ is characterized by its capability to handle unknown control directions within a control system. For a continuously differentiable function $N(\zeta):[0,\infty)\to(-\infty,\infty)$, it is defined via its positive and negative truncated forms, $N^{+}(\zeta)=\max\{0,N(\zeta)\}$ and $N^{-}(\zeta)=\max\{0,-N(\zeta)\}$, respectively. The function satisfies the condition that

$$N(\zeta)=N^{+}(\zeta)-N^{-}(\zeta),$$ (4)

with the properties that for any $\zeta$,

$$\lim_{v\to\infty}\sup\frac{1}{v}\left[\int_{0}^{v}N^{+}(\zeta)d\zeta-\int_{0}^{v}N^{-}(\zeta)d\zeta\right]=\infty,$$ (5)

$$\lim_{v\to\infty}\sup\frac{1}{v}\left[\int_{0}^{v}N^{-}(\zeta)d\zeta-\int_{0}^{v}N^{+}(\zeta)d\zeta\right]=\infty.$$ (6)

A function that satisfies these conditions is utilized for managing the uncertainty in control direction, offering a robust approach to control system design.

Lemma 2

[16, 17] Let $V(t)$ and $\zeta(t)$ be smooth functions defined over the interval $[0,t_{f})$ with $V(t)>0$ for all $t\in[0,t_{f})$. Given $N(\zeta(t))$ as a Nussbaum-type function, if for any $t\in[0,t_{f})$, the following condition is met:

$$V(t)<c_{0}+e^{-c_{1}t}\int_{0}^{t}\left(gN(\zeta(\tau))+1\right)\dot{\zeta}(\tau)e^{c_{1}\tau}d\tau,$$ (7)

where $c_{0}>0$ and $c_{1}>0$ are positive constants, and $g$ is a control parameter within closed intervals $L=[l_{-},l_{+}]$ excluding zero ($0\notin L$), then $V(t)$, $\zeta(t)$, and the integral $\int_{0}^{t}gN(\zeta(\tau))\dot{\zeta}(\tau)e^{c_{1}\tau}d\tau$ are guaranteed to be bounded on $[0,t_{f})$.

III Nussbaum Function-based Control Design and Stability Analysis

We propose a Nussbaum function based PID-like control input for robot manipulator control, delineated as:

$$u(t)= \left(k_{\pi}+\kappa_{\pi}(t)\right)K_{N}(\zeta)e(t)+\left(k_{\iota}+\kappa_{\iota}(t)\right)K_{N}(\zeta)\int_{0}^{t}e(\rho)\,d\rho$$ (8)

$$+\left(k_{\Delta}+\kappa_{\Delta}(t)\right)K_{N}(\zeta)\frac{d}{dt}e(t),$$

with $K_{N}(\zeta)=-N(\zeta)$ is the related Nussbaum function gain. In (8) we employed constant gains $k_{\pi}$, $k_{\iota}$, and $k_{\Delta}$ alongside their time-varying analogs $\kappa_{\pi}(t)$, $\kappa_{\iota}(t)$, and $\kappa_{\Delta}(t)$. This strategy enhances traditional PID controls by integrating adaptive gains that adjust in real-time to the system’s state. To further simplify parameter tuning, we establish the relationships $k_{\Delta}=\frac{k_{\pi}}{2\gamma}=\frac{k_{\iota}}{\gamma^{2}}$ and $\kappa_{\Delta}(t)=\frac{\kappa_{\pi}(t)}{2\gamma}=\frac{\kappa_{\iota}(t)}{\gamma^{2}}$, ensuring the quadratic expression $S^{2}+2\gamma S+\gamma^{2}$ is Hurwitz polynomial, where $S$ is the Laplace operator. Thus, the control design primarily depends on selecting $k_{\Delta}$ and $\kappa_{\Delta}(t)$. The PID-like control input is then refined to:

$$u(t)= \left(k_{\Delta}+\kappa_{\Delta}(t)\right)K_{N}(\zeta)\left(2\gamma e(t)+\gamma^{2}\int_{0}^{t}e(\rho)\,d\rho+\frac{d}{dt}e(t)\right)$$ (9)

$$= -\left(k_{\Delta}+\kappa_{\Delta}(t)\right)N(\zeta)\Psi.$$

This approach significantly reduces the complexity of the gain selection process, focusing on only two key parameters. The time-varying gain $\kappa_{\Delta}(t)$ is adaptively updated by

$$\kappa_{\Delta}(t)=-\alpha\hat{\psi}(t)^{T}\phi(x),$$ (10)

and the adaptive law for $\hat{\psi}(t)$ is given by

$$\dot{\hat{\psi}}(t)=-\Gamma(\alpha\|\Psi(t)\|^{2}\phi(x)+\sigma\hat{\psi}(t)).$$ (11)

with $\alpha$ and $\phi$ are positive control constants, and $\Gamma=\Gamma^{T}>0$.

In this work, we choose $N(\zeta)=\zeta^{2}\cos(\zeta)$ with the property of $N(0)=0$ as our Nussbaum function. The updating law for $\zeta$ is given by

$$\dot{\zeta}(t)=\Psi(t)^{T}(k_{\Delta}+\kappa_{\Delta}(t))\Psi(t).$$ (12)

The above control framework, enhanced by adaptive laws based on the low complexity Nussbaum function PID control, provides us with an effective approach to robotic manipulator control, as formalized in the forthcoming theorem.

Theorem 1

Consider the robot manipulator system described by (1) with the Properties 1-3, satisfying to Assumptions 1 and 2. If the Nussbaum function-based PID-like control law (9), with the adaptive updating laws (11), and (12), and and Lemmas 1, and 2 is implemented, then, with the design parameters are properly chosen, the closed-loop system remains stable, and all signals within the system are uniformly bounded. Furthermore, the joint position signal $q(t)$ closely tracks the desired trajectory $q_{d}(t)$, with the tracking error converging to a small neighborhood around zero.

Proof:

Consider the Lyapunov candidate as

$$V(t)=\frac{1}{2}\Psi(t)^{T}M(q(t))\Psi(t)+\frac{1}{2}\tilde{\psi}(t)^{T}\Gamma^{-1}\tilde{\psi}(t),$$ (13)

where $\tilde{\psi}(t)=\hat{\psi}(t)-\psi^{*}$ with $\psi^{*}(t)$ is the ideal constant weight vector of neural network approximation.

The time derivative of the Lyapunov function $V(t)$ can be bounded as

$${\dot{V}(t)}={\Psi(t)^{T}}M(t)\dot{\Psi}(t)+\frac{1}{2}{\Psi(t)^{T}}\dot{M}(t)\Psi(t)+\tilde{\psi}(t)^{T}\Gamma^{-1}\dot{\hat{\psi}}(t).$$ (14)

Considering dynamics (1) and the definition of the generalized error $\Psi(t)$, one can obtain $M(t)\dot{\Psi}(t)=C(t)\dot{q}(t)+G(t)-\tau(t)+M(t)\bar{e}(\cdot)$, where $\bar{e}(\cdot)=\ddot{q}_{d}(t)+2\gamma\dot{e}(t)+\gamma^{2}e(t)$. Then, utilizing the application of Young’s inequality, and considering Property 2, we obtain $\Psi(t)^{T}C(t)\dot{q}(t)\leq\alpha\|\Psi(t)\|^{2}c^{2}\|\dot{q}(t)\|^{4}+\frac{1}{4\alpha}$, $\Psi(t)^{T}G(t)\leq\alpha\|\Psi(t)\|^{2}g^{2}+\frac{1}{4\alpha}$, $\Psi(t)^{T}M(t)\bar{e}(\cdot)\leq\alpha\|\Psi(t)\|^{2}m^{+^{2}}\|\bar{e}(\cdot)\|^{2}+\frac{1}{4\alpha}$, and $\frac{1}{2}\Psi(t)^{T}\dot{M}(t)\Psi(t)=\Psi(t)^{T}C(t)\Psi(t)\leq\alpha\|\Psi(t)\|^{2}c^{2}\|\dot{q}(t)\|^{2}\|\Psi(t)\|^{2}+\frac{1}{4\alpha}$, where $\alpha>0$ is a design parameter. Accordingly, taking into account equation $\eqref{eq_vd}$ and applying the aforementioned inequalities, one has,

$$\dot{V}(t)\leq\alpha\|\Psi(t)\|^{2}\Lambda(t)-\Psi(t)^{T}\kappa u(t)+\frac{1}{\alpha}+\Gamma^{-1}\tilde{\psi}(t)^{T}\dot{\hat{\psi}}(t),$$ (15)

where

$\Lambda(t)=-\left(c^{2}\|\dot{q}(t)\|^{4}+g^{2}+m^{+2}\|\bar{e}(t)\|^{2}+c^{2}\|\dot{q}(t)\|^{2}\|\Psi(t)\|^{2}\right)$.

Subsequently, considering neural networks approximation $\Lambda(t)=\psi^{*}\phi(x)+\varepsilon(x)$, with $\varepsilon\left(x\right)$ being the unknown approximation error which is upper bounded in the sense that $\left\|\varepsilon\left(x\right)\right\|\leq\varepsilon_{m}$, and utilizing (11), (12), and (9) into (15), and considering $-\tilde{\psi}(t)^{T}{\hat{\psi}}(t)\leq-\tilde{\psi}(t)^{T}\tilde{\psi}(t)+|\psi^{*}\|^{2}$ [4], yields

$$\dot{V}(t) \leq\dot{\zeta}(t)-\|\Psi(t)\|^{2}(k_{\Delta}+\kappa_{\Delta}(t))-\alpha\|\Psi(t)\|^{2}(\psi^{*^{\top}}\phi(x)+\varepsilon)$$

$$+\Psi(t)^{T}\kappa N(\zeta(t))(k_{\Delta}+\kappa_{\Delta}(t))\Psi(t)+\frac{1}{\alpha}$$

$$-\alpha\|\Psi(t)\|^{2}\tilde{\psi}(t)^{T}\phi(x)-\frac{1}{2}\sigma\tilde{\psi}(t)^{T}\tilde{\psi}(t)+\frac{1}{2}\sigma\|\psi^{*}\|^{2}.$$

Considering $\|\Psi(t)\|^{2}\kappa_{\Delta}(t)=-\|\Psi(t)\|^{2}\alpha\hat{\psi}(t)^{T}\phi(x)$, the above equation can be rewritten as

$$\dot{V}(t) \leq\dot{\zeta}(t)+\kappa N(\zeta(t))\dot{\zeta}(t)-\|\Psi(t)\|^{2}(k_{\Delta}+\alpha\varepsilon)$$

$$-\frac{1}{2}\sigma\ \|\tilde{\psi}(t)\|^{2}+\frac{1}{\alpha}+\frac{1}{2}\sigma\bar{\psi}^{2},$$

with $\bar{\psi}$ being the upper bound of the optimal weight $\|\psi^{*}\|$.

Finally, the above equation can be formed as

$$\dot{V}(t)<\dot{\zeta}(t)+\kappa N(\zeta(t))\dot{\zeta}(t)-c_{1}V(t)+c_{2},$$ (16)

where, $c_{1}=\min\left(\frac{2(\alpha\varepsilon_{m}+k_{\Delta})}{m^{+}},\frac{\sigma}{\lambda_{\max}(\Gamma^{-1})}\right)$ and $c_{2}=\frac{1}{\alpha}+\frac{1}{2}\sigma\bar{\psi}^{2}$ and both are positive. Multiplying (16) by $e^{c_{1}t}>0$, yields

$$\frac{d}{dt}\left(V(t)e^{c_{1}t}\right)\leq\dot{\zeta}(t)e^{c_{1}t}+\kappa N(\zeta(t))\dot{\zeta}(t)e^{c_{1}t}+c_{2}e^{c_{1}t},$$ (17)

and integrating both sides of (17) over $[0,t]$, leads to

$$V(t)\leq e^{-c_{1}t}\int_{0}^{t}(\kappa N(\zeta(t))+1)\dot{\zeta}(t)e^{c_{1}\tau}d\tau$$ (18)

$$+\left(V(0)-\frac{c_{2}}{c_{1}}\right)e^{-c_{1}t}+\frac{c_{2}}{c_{1}}.$$

Since $0<e^{-c_{1}t}\leq 1$, the inequality (18) can be rewritten as,

$$V(t)\leq e^{-c_{1}t}\int_{0}^{t}(\kappa N(\zeta(t))+1)\dot{\zeta}(t)e^{c_{1}\tau}d\tau+c_{0},$$ (19)

where, $c_{0}=\frac{c_{2}}{c_{1}}+V(0)$ is a positive constant. Utilizing Lemma 2 we can conclude from (19) that $V(t)$, $\zeta(t)$, and $\int_{0}^{t}\left(\kappa N(\zeta(t))+1\right)\dot{\zeta}(t)\tau d\tau$ are bounded on $[0,t_{f})$. Considering the Lyapunov function (13), it holds that $\Psi(t)$ and $\tilde{\psi}(t)$ are bounded. Thus, since ${\psi}^{*}(t)$ is bounded, then, $\hat{\psi}(t)$ is bounded. In addition, using Lemma 1, the boundedness of $\Psi(t)$ ensures that $e(t)$, $\dot{e}(t)$, and $\int_{0}^{t}e(\cdot)d\tau$ are bounded. Then, considering Assumption 2, the boundedness of $e(t)$, and $\dot{e}(t)$, ensures $q(t)$, and $\dot{q}(t)$ are bounded. Furthermore, considering (10), (11), (12), and (9), and considering the boundedness of the basis function vector $\phi(x)$ (see [4] for the reference), and the boundedness of $\Psi(t)$, $\hat{\psi}(t)$, then $\kappa_{\Delta}(t)$, ${\psi}(t)$, $\zeta(t)$, and control $u(t)$ are bounded. Finally, as $u(t)$ is bounded, the boundedness of $\tau(t)$ is ensured, and accordingly, all closed-loop signals are bounded.

To prove that the error $e(t)$ converges to a small neighborhood of zero, considering the boundedness of $V(0)$ and $\int_{0}^{t}\left(\kappa N(\zeta(t))+1\right)\dot{\zeta}(t)\tau d\tau$, along with Properties 1 and 3, it is concluded from (13) and (19) that

$$\lim_{t\to\infty}\frac{1}{2}\|\Psi(t)\|^{2}\leq\frac{c_{2}}{c_{1}},$$

follows by $\|\Psi(t)\|\leq\sqrt{\frac{2c_{2}}{c_{1}}}:=\Psi_{m}$. Then, considering (3), as detailed in [25], it is concluded that the error $e(t)$ converges to close to zero. ∎

IV Numerical Simulation

In this section, numerical simulations are performed to verify the effectiveness of the proposed Nussbaum-based PID control, as established in Theorem 1. We use a two-link robot manipulator situated in the vertical plane for our simulation study. Physical parameters of the robot were selected as follows: the masses of the links $m_{1}=5\ \text{kg},\ m_{2}=2\ \text{kg}$, the lengths of the links $l_{1}=1\ \text{m},\ l_{2}=0.75\ \text{m}$, and the inertias of the links $I_{1}=1.66\ \text{m}^{2},\ I_{2}=0.37\ \text{m}^{2}$. The reference trajectories are chosen as $q_{d}=[\cos(t);\ -\cos(t)]$, with the initial conditions for each joint given by $q(0)=[\pi/2;-\pi/2]$,  and $\dot{q}(0)=[0;0]$. Control parameters are chosen to be $\alpha=100$, $\Gamma=100$, $k_{\Delta}=0.1$, $\gamma=0.5$, and $\sigma=0.1$. Also, a radial basis function neural network with twenty nodes in each hidden layer is selected so that centers are evenly distributed in the span of the input space $[-12.5,12.5]$, and widths 1. The input vector of the neural network is chosen as $x=\left[e^{T},\dot{e}^{T},q^{T},\Psi^{T}\right]$. The initial points of neural network weights were chosen as $\hat{\psi}\left(0\right)=0$. Simulation results are shown in Figures 1-5. In these figures, indices 1 and 2 denote the first and second links of the robot manipulator, respectively.

The tracking performances of the links are depicted in Figures 1- 3. Figures 1, and 2 illustrate that the actual position and velocity signals closely follow their desired trajectories. Figure 3 illustrates the boundedness of the PID generalized error. Figure 4 shows the boundedness of neural network weights. Finally, the input control is depicted in Figure 5. The above figures show the tracking performance using our developed control. It also shows that the input control vectors are bounded, demonstrating the proposed method’s capability to accomplish the control tasks effectively.

Figure 1: Desired and actual trajectories of joint positions

Figure 2: Desired and actual trajectories of joint velocities

Figure 3: Trajectories of the filtered error $\Psi(t)$

Figure 4: Norms of neural networks weignts

Figure 5: Trajectories of control input

V Conclusion

We presented a novel control strategy for robotic manipulators, leveraging a Nussbaum function-based PID approach to address unknown control directions and dynamics. This strategy simplifies the control design by minimizing the number of tuning parameters and utilizing direct Lyapunov analysis for stability assurance. Our contributions include a simple yet effective PID control framework, automatic gain adjustment through neural network-based estimation, and a reduction in parameter tuning complexity, enhancing adaptability and robustness in uncertain environments. The analytical and numerical validations highlight the effectiveness of our approach in enhancing robotic manipulator control. These findings lay a robust groundwork for future investigations in this field, setting the stage for experimental verification and the exploration of further control innovations.