Next Article in Journal
Effective Noise Reduction in NDR Systems: A Simple Yet Powerful Apriori-Based Approach
Next Article in Special Issue
Transportation Mode Detection Using Learning Methods and Self-Contained Sensors: Review
Previous Article in Journal
EHNet: Efficient Hybrid Network with Dual Attention for Image Deblurring
Previous Article in Special Issue
Deep Learning Based Over-the-Air Training of Wireless Communication Systems without Feedback
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Active Disturbance Rejection Control via Neural Networks for a Lower-Limb Exoskeleton

by
Karina I. Espinosa-Espejel
1,
Yukio Rosales-Luengas
1,
Sergio Salazar
1,
Ricardo Lopéz-Gutiérrez
2 and
Rogelio Lozano
1,*
1
Department of Research and Multidisciplinary Studies, Center for Research and Advanced Studies of the National Polytechnic Institute, Mexico City 07360, Mexico
2
CONAHCYT-INAOEP, San Andrés Cholula 72840, Mexico
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(20), 6546; https://doi.org/10.3390/s24206546
Submission received: 30 July 2024 / Revised: 1 October 2024 / Accepted: 9 October 2024 / Published: 11 October 2024

Abstract

:
This article presents the design of a control algorithm based on Artificial Neural Networks (ANNs) applied to a lower-limb exoskeleton, which is aimed to carry out walking trajectories during lower-limb rehabilitation. The interaction between the patient and the exoskeleton leads to model uncertainties and external disturbances that are always present. For this reason, the proposed control considers that the non-linear part of the model is unknown and is perturbed by external disturbances, which are estimated by an active disturbance rejection control via Artificial Neural Networks. To validate the proposed approach, a numerical simulation and an experimental implementation of the ANN-Controller are developed.

1. Introduction

Exoskeletons are used to assist patients with motor disabilities, such as paraplegia or quadriplegia, during their therapy sessions. These robotic devices facilitate the execution of different movements to regain locomotion ability. Lower-limb exoskeletons provide benefits during walking-assisted execution. During therapy, the features of the exoskeleton are adjusted to the patient’s requirements, such as velocity and the number of walking cycles. Human–machine interaction leads to the appearance of uncertain dynamics that are not considered in the dynamical model of an exoskeleton. Moreover, there are external disturbances, such as the vibrations of a treadmill during the walking cycle, that have an impact on the exoskeleton performance. Therefore, a research area of interest is the development of controllers with active disturbance rejection [1].
The standard Active Disturbance Rejection Controller (ADRC) depends on the system model and a virtual variable that contains internal uncertainties, unmodeled dynamics, and external disturbances. The total uncertainty is estimated by an Extended State Observer (ESO), and the disturbances are rejected in real time [2]. The ADRC is a strategy that can be applied even when the precise mathematical description of the model is not available. The three parts of the scheme are the Tracking Differentiator (TD), the ESO, and the ESO-based feedback control. The TD provides the reference system, the ESO estimates the disturbances, and the ESO-based feedback control compensates the disturbances and tracks the desired trajectories. Some of the ADRC approaches developed in recent years use different types of ESO and feedback control schemes, such as linear, nonlinear, and fractional order schemes [3,4].
In [5], an active disturbance rejection controller with a linear Extended State Observer, a tracking differentiator, and a nonlinear state error feedback applied to a lower-limb exoskeleton of 2 Degrees of Freedom (DoFs) is developed, and the controller eliminates external disturbances and parameter variations. They test the controller with a sinusoidal function noise of −0.5° to +0.5° of variance with +20% parameter variation. Furthermore, they test a random disturbance between −1 to 1, a constant disturbance, and an harmonic disturbance.
The authors in [6] developed an active disturbance rejection controller based on a fast terminal sliding mode controller for a lower-limb exoskeleton, which is able to improve tracking performance and to converge to a bounded region. Moreover, a fractional proportional derivative-based active disturbance rejection control of a knee exoskeleton device for rehabilitation is developed in [7].
An active disturbance rejection controller based on the flatness property and a control Lyapunov function is presented in [8], which is implemented in an actuated ankle orthosis for paraplegic patients. The stability analysis is carried out with the Input-to-State Stability (ISS) method and an Extended State Observer (ESO), which estimates unknown disturbances, and this estimation is added to the control input. In [9], a piecewise Extended State Observer to estimate disturbance and uncertainty is implemented in a double-joint dexterous hand, which is driven by pneumatic artificial muscles. A proportional derivative controller with disturbance estimation is applied to the system to improve angle tracking and disturbance rejection precision. To verify the disturbance rejection ability, an object of 50 g is added to the first joint as an abrupt external disturbance.
Furthermore, a Robust Unknown Input Observer (RUIO) that obtains the decoupled estimation of the unknown disturbance, which is implemented in a robotic manipulator, is presented in [10]. In that paper, the authors develop a state-feedback control using the linear parameter-varying method to avoid torque saturation through gain shifting. Then, the authors in [11] use a control strategy based on an ESO and a finite time stable tracking differentiator (FTSTD) in an upper limb exoskeleton. The exoskeleton performs flexion and extension for the shoulder and elbow in the sagittal plane. Additionally, ADRC schemes are used in applications such as a fluid-driven hand rehabilitation device, which implements an ADRC based on an ESO to achieve the control purpose due to the disturbances generated by the unconscious tremors and the difficulty in obtaining the exact model [12].
Artificial Neural Networks (ANNs) are universal estimators of a given function [13]. The ANN estimates a function by updating the weights to reduce the error between the output and the desired function through a learning rule. Due to the ANN estimation ability, they are used in ADRC schemes, and some control applications implement ANN as estimators using the tracking error to update the weights. Different applications have implemented ADRC schemes based on ANN. In the field of biomedicine and rehabilitation, artificial intelligence is gaining ground. In [14], the authors propose a novel online real-time gait terrain detection algorithm from the measurements of a foot-mounted Inertial Measurement Unit (IMU), using a shallow cascaded Convolutional Neural Network and Long Short-Term Memory neural network (CNN-LSTM).
In [15], the authors employ a human-in-the-loop scheme with a combination of a Recurrent Neural Network (RNN) and an Adaptive Non-singular Fast Terminal Sliding Mode Controller (ANFTSMC) strategy to classify the user’s movements to control the trajectories of an exoskeleton. The authors in [16] present a neural PID control with application in an upper limb exoskeleton. They obtain a semi-globally stable neural PID controller and a locally asymptotically stable neural PID controller with a velocity observer and test standard weight training algorithms.
An adaptive command filtering control based on neural networks is designed in [17], which performs trajectory tracking of an arm robot with flexible pneumatic artificial muscles. A barrier term is introduced so that the tracking errors converge to neighborhoods of the origin, and these are also bounded. They obtained a continuous controller that is able to work with unmodeled dynamics, uncertainties, and multiple input constraints. The robustness of the method is validated on a humanoid arm robot, the control requires the outputs to be measurable, and noises caused by compressed air should be suppressed. In [18], the authors present an adaptive neural network based on the integral non-singular terminal sliding mode control method to perform a fixed-time position tracking control in a wheelchair upper-limb exoskeleton under external disturbance.
The authors in [19] use an ANN to optimize the tuning of PID gains in a multi-joint lower-limb exoskeleton for gait rehabilitation, and the different tuning methods are compared with each other. The ANN-based methods improve the results and reduce the steady-state errors. A lower limb of 12 DoFs is proposed in [20], which uses a controller based on the Time-Delay Estimation (TDE), a Computed Torque Control (CTC), and a Robust Adaptive Radial Basis Neural Network (RBNN) control. The TDE is used to estimate unmodeled dynamics and external disturbances. The robust adaptive compensator based on the RBNN approximates and compensates the error of the TDE. The TDE-based CTC shows improved results in comparison with the CTC and the sliding mode based on the CTC.
An adaptive neural network control based on an RBNN for motion tracking of a lower-limb exoskeleton is presented in [21], and the RBNN ensures the accuracy and robustness of the controller, and a disturbance observer eliminates the external disturbances online. In [22], the authors present a data-driven reinforcement learning control to work with different hemiplegic patients and to deal with unpredictable disturbances. A policy interaction algorithm learns the optimal control to achieve online adaptation.
Compared with other works, the ADRC based on an ANN considers a feedback controller linearization and the ANN that is capable of estimating the unknown function by using the tracking error. Moreover, the ADRC-ANN requires less information about the model and has a good learning capacity for nonlinear characteristics and external disturbances. Additionally, the ANN structure can be established by the designer; thus it properly estimates the disturbance.
The present paper proposes a robotic device for walking rehabilitation concerning the lower-limb exoskeleton (LLE). The LLE is used to assist patients with motor disability during the walking cycle, and this promotes the motor recovery. The contributions of this paper are listed as follows:
  • A control scheme based on an ANN is proposed to estimate the external disturbances, the nonlinear exoskeleton model, and the human behavior. Furthermore, a feedback linearization controller is proposed to perform a smooth tracking trajectory to prevent injuries.
  • The stability and robustness of the proposed approach are obtained using Lyapunov stability analysis.
  • A rejection disturbance strategy ensures the feasibility and safety of the robot during walking rehabilitation, dealing with disturbances such as vibrations of the treadmill.
Unlike other works, we added a platform to move the patient from a wheelchair, and the system has elastic joints to smooth the movement and sense the intention of motion.
This paper is organized as follows: the dynamic model is presented in Section 2. Section 3 presents the design of the ANN-based controller, and the numerical simulation and experimental implementation results are reported in Section 4. Finally, the conclusions are presented in Section 5.

2. Exoskeleton Dynamic Model and Properties

2.1. System Description and Modeling

The exoskeleton used in this work was designed by the Research and Advanced Studies Center (https://umi.cinvestav.mx/ accessed on 20 September 2024). The lower-limb exoskeleton, manufactured using 6061-T6 aluminum alloy, has 6 active degrees of freedom (DoFs); 4 of them correspond to the lower limbs (2 for the knees and 2 for the hip), and two of them are part of the standing mechanism. The ankle joint is considered for future tests, and a mechanical support is attached to this joint to sustain the foot, since the first tests are performed on a healthy subject. The actuators for the hip and the knees are of rotational elastic type, which consists of harmonic drive motors with springs. A property of these actuators is the damping during the walking cycle. In order to obtain the angular position of each joint, two AMT203-V, CUI-devices, Oregon UE. absolute rotary magnetic encoders were used, one to read the angular position of the motor and one to read the angular position of the link.
Figure 1 shows the final prototype and the general distribution of the actuators. The prototype is divided in two sections: lifting and lower-limb sections. The lifting section consists of a parallel mechanism and a safety harness attached to a crane that is used to lift up the patient. The lower-limb section consists of an exoskeleton that supports the leg, and the axis of the actuators coincides with the axis of rotation of the joints.
The exoskeleton has an anthropomorphic approach to simplify the adjustment to the patient’s legs, and it executes different exercises such as walking, sitting, and getting up. A schematic diagram of the exoskeleton is presented in Figure 2, in which there are two links corresponding to the thigh (link 1) and the leg (link 2) with lengths l 1 and l 2 , with masses m 1 and m 2 , respectively.
The angular positions q 1 ( t ) of the hip and  q 2 ( t ) of the knee constitute the vector of angular positions q ( t ) = [ q 1 ( t ) q 2 ( t ) ] T . The dynamics of the links are actuated by the angular motion of the motors q m ( t ) = [ q m 1 ( t ) q m 2 ( t ) ] T through flexible joints, and the dynamics of the motors are driven by the motor torques τ . Then, the dynamic equations of the exoskeleton are as follows:
M ( q ) q ¨ + C ( q , q ˙ ) q ˙ + g m ( q ) + K [ q q m ] = 0 J q ¨ m + B q ˙ m K [ q q m ] = τ + d ,
where M ( q ) R 2 × 2 is the inertia matrix, C ( q , q ˙ ) R 2 × 2 is the Coriolis matrix, g m ( q ) R 2 × 1 is the vector of gravitational forces, and τ R 2 × 1 are the torques of the actuators. J is a positive definite matrix with values in its diagonal equal to the product of the inertia moments of the motors and the square of the gear ratio J = d i a g { J 1 r 1 , 1 2 , J 2 r 1 , 2 2 } . K is a diagonal positive definite matrix with the elastic constants on the diagonal defined as K = d i a g { k s 1 , k s 2 } , B is a positive definite diagonal matrix, and its diagonal elements are the viscous friction parameters of each motor B = d i a g { b 1 , b 2 } . Additionally, d R 2 × 1 is an unknown term, which contains non-modeled dynamics, external noise, and human interaction. Then, the vector d is defined as follows:
d = τ h + D + ρ ,
where τ h R 2 × 1 is the vector of the torque generated by the human, D R 2 × 1 is the vector of external disturbances, and ρ is the unknown part of the model.
Furthermore, the elements of M ( q ) , C ( q , q ˙ ) , and g m ( q ) are as follows:
M 11 ( q ) = I 1 + I 2 + l 1 2 m 2 + l c 1 2 m 1 + l c 2 2 m 2 + 2 l 1 l c 2 m 2 c o s ( q 2 ) M 12 ( q ) = I 2 + l c 2 2 m 2 + l c 2 2 m 2 + l 1 l c 2 m 2 c o s ( q 2 ) M 21 ( q ) = I 2 + l c 2 2 m 2 + l c 2 2 m 2 + l 1 l c 2 m 2 c o s ( q 2 ) M 22 ( q ) = I 2 + m 2 l c 2 2 C 11 ( q , q ˙ ) = m 2 l 1 l c 2 s i n ( q 2 ) q ˙ 2 C 12 ( q , q ˙ ) = m 2 l 1 l c 2 s i n ( q 2 ) [ q ˙ 1 + q ˙ 2 ] C 21 ( q , q ˙ ) = m 2 l 1 l c 2 s i n ( q 2 ) q ˙ 1 C 22 ( q , q ˙ ) = 0 g m 1 ( q ) = [ m 1 l c 1 + m 2 l 1 ] g s i n ( q 1 ) + m 2 g l c 2 s i n ( q 1 + q 2 ) g m 2 ( q ) = m 2 g l c 2 s i n ( q 1 + q 2 )
The desired trajectories corresponding to the hip and knee joints are denoted by q d 1 and q d 2 , respectively.
Solving the second derivative of links and motor position, and rewriting the dynamic equations in (1), the following expressions for q ¨ and q ¨ m are obtained:
q ¨ = M 1 [ C q ˙ g m K ( q q m ) ]
q ¨ m = J 1 [ τ + d B q ˙ m + K ( q q m ) ]
By using (4a) and (4b), the following expression is obtained:
q ( 4 ) = M ¨ 1 [ C q ˙ g m K ( q q m ) ] + 2 M ˙ 1 [ C q ¨ C ˙ q ˙ g ˙ m K ( q ˙ q ˙ m ) ] + M 1 [ 2 C ˙ q ¨ C q ( 3 ) C ¨ q ˙ g ¨ m K ( q ¨ J 1 [ τ + d B q ˙ m + K ( q q m ) ] ) ]
Furthermore, the input control τ depends on q, q ˙ , q ¨ , q ( 3 ) , and q ( 4 ) . The state variables are assigned as x 1 = q , x 2 = q ˙ , x 3 = q ¨ , and  x 4 = q ( 3 ) . Hence, the variable state vector is x = [ q q ˙ q ¨ q ( 3 ) ] T R 8 × 1 , and the state space representation is as follows:
x ˙ = x ˙ 1 x ˙ 2 x ˙ 3 x ˙ 4 = q ˙ q ¨ q ( 3 ) q ( 4 ) = x 2 x 3 x 4 f ( x ) + g ( x ) τ ,
where
f ( x ) = M ¨ 1 [ C x 2 g m K ( x 1 q m ) ] + 2 M ˙ 1 [ C x 3 C ˙ x 2 g ˙ m K ( x 2 q ˙ m ) ] + M 1 [ 2 C ˙ x 3 C x 4 C ¨ x 2 g ¨ m K ( x 3 J 1 [ d B q ˙ m + K ( x 1 q m ) ] ) ] g ( x ) = M 1 K J 1
Assumption 1. 
The external disturbance d is bounded.

2.2. Artificial Neural Network Structure for the LLE

An Artificial Neural Network (ANN) is applied to approximate the unknown function of the robot. The ANN is a universal estimator, and the structure of the neural network is as follows:
Δ Q = W T ϕ + ϵ R 2 × 1 ,
where W R N × 2 is the weight matrix and N is the number of nodes in the hidden layer.
According to the approximation property, given a continuous real function y : Ω R in a compact set Ω R n , by choosing ϕ j , j = 1 , , N , there exists an ideal weight matrix that can approximate the function y with a bounded approximation error ϵ > 0 . The matrix of the weights is W = [ w j k ] , where j = { 1 , 2 , , N } and k = { 1 , 2 } ; ϵ is an approximation error ϵ = [ ϵ 1 ϵ 2 ] T .
Additionally, ϕ R N is the activation function vector, and its components are as follows:
ϕ = [ ϕ 1 ϕ 2 ϕ N ] T
In this paper, the activation function is selected to be a sigmoid function, defined by the following:
ϕ j ( v j ) = a j 1 + e b j v j c j ,
where a j , b j , and c j are the parameters of the sigmoid functions.
The implemented ANN has three layers, the input layer, the hidden layer with the activation functions, and the output layer. The inputs of the neural network are denoted by v j . The ANN approximates the unknown function of the uncertainties and the non-linearities of the system. As the weight matrix W is unknown, it is estimated on-line, and it is denoted by W ^ . Moreover, an adaptive law is obtained to update the weight matrix. Therefore, the output of the on-line artificial neuron is as follows:
Δ Q ^ = W ^ T ϕ R 2 × 1
Remark 1. 
We consider that the next parameters are bounded; the approximation error ϵ < ϵ N , the weight matrix W = W M , and the activation function ϕ ( v ) σ . Additionally, we consider that ϵ N , W M , and σ are positive constants.

3. ANN-Based Control Design

In this section, we describe the design of an active disturbance rejection control via artificial neural networks for the lower-limb exoskeleton (LLE). The control objective is to track the desired trajectory such that the closed-loop system is uniformly ultimately bounded and the output of the system q follows the desired trajectory q d .
The diagram in Figure 3 describes the proposed control approach. The control scheme has three parts: the linear feedback control, the ANN estimation of a variable that represents the external disturbances and unmodeled dynamics ( Δ Q ), and the flexible joint robot model (Equation (1)). The linear feedback control includes the terms of desired position, velocity, acceleration, and the fourth derivative (jerk) of the position error.
The desired trajectories for the walking cycle are trigonometric polynomials of the 8th order [23]. Afterwards, the desired trajectory of the hip is q d 1 and the desired trajectory of the knee is q d 2 , as follows:
q d 1 = 17.2 s i n ( 0.08 ω t + 0.24 ) + 69.61 s i n ( 1.00 ω t + 1.72 ) + 6.90 s i n ( 0.17 ω t + 2 ) + 4 s i n ( 2 ω t 1.72 ) + 1.68 s i n ( 2.99 ω t 0.01 ) + 1.63 s i n ( 0.36 ω t + 2.31 ) + 0.70 s i n ( 0.44 ω t + 4.24 ) + 48.17 s i n ( 1.01 ω t + 4.81 ) q d 2 = 44.7 s i n ( 0.08 ω t + 0.36 ) + 21.13 s i n ( 0.99 ω t 2.92 ) + 16.07 s i n ( 2 ω t 0.94 ) + 20.08 s i n ( 0.15 ω t + 2.44 ) + 4.37 s i n ( 2.99 ω t 0.02 ) + 20.92 s i n ( 0.34 ω t + 2.98 ) + 18.45 s i n ( 0.35 ω t + 6 ) + 1.27 s i n ( 3.98 ω t 0.65 ) ,
where ω is the frequency of the walking cycle.

3.1. ANN-Based Control Design

In order to design the ANN-based control, the auxiliary linear dynamics are proposed as follows:
q ^ ( 4 ) = τ
The error between the auxiliary dynamics and the exoskeleton dynamics is denoted as Δ Q = q ( 4 ) q ^ ( 4 ) and is described as follows:
Δ Q = M ¨ 1 [ C q ˙ g m K ( q q m ) ] + 2 M ˙ 1 [ C q ¨ C ˙ q ˙ g ˙ m K ( q ˙ q ˙ m ) ] + M 1 [ 2 C ˙ q ¨ C q ( 3 ) C ¨ q ˙ g ¨ m K ( q ¨ J 1 [ τ + d B q ˙ m + K ( q q m ) ] ) ] τ
Then,
q ( 4 ) = q ^ ( 4 ) + Δ Q q ( 4 ) = τ + Δ Q
A control input is proposed as follows:
τ * = v τ Δ Q
Substituting the above into (15) leads to the following:
q ( 4 ) = v τ
Considering the state variable assignment proposed in Section 2.1: x 1 = q , x 2 = q ˙ , x 3 = q ¨ and x 4 = q ( 3 ) , and Equation (17), it follows:
x ˙ 1 = x 2 x ˙ 2 = x 3 x ˙ 3 = x 4 x ˙ 4 = v τ y = x 1
The input v τ is proposed such that q follows the smooth desired trajectory q d when t . Therefore, v τ is proposed as follows:
v τ = q d ( 4 ) K 1 ( q q d ) K 2 ( q ˙ q ˙ d ) K 3 ( q ¨ q ¨ d ) K 4 ( q ( 3 ) q d ( 3 ) ) ,
where K i = d i a g { k i , k i } with i = 1 , 2 , 3 , 4 are the gain matrices such that each gain matrix is a diagonal positive definite matrix. Substituting v τ from (19) into Equation (17), we have the following:
q ( 4 ) q d ( 4 ) = K 1 ( q q d ) K 2 ( q ˙ q ˙ d ) K 3 ( q ¨ q ¨ d ) K 4 ( q ( 3 ) q d ( 3 ) )
The error variables are defined as q ˜ = q q d , q ˜ ˙ = q ˙ q ˙ d , q ˜ ¨ = q ¨ q ¨ d , and  q ˜ ( 3 ) = q ( 3 ) q d ( 3 ) . Let us define the error vector as follows:
e = [ q ˜ T q ˜ ˙ T q ˜ ¨ T q ˜ ( 3 ) T ] T R 8 × 1 ,
and the derivative of the error vector is as follows:
e ˙ = [ q ˜ ˙ T q ˜ ¨ T q ˜ ( 3 ) T q ˜ ( 4 ) T ] T R 8 × 1
Using the error vector and its derivative in (21) and (22) with (20), the dynamic equation for the tracking errors is as follows:
e ˙ = A e ,
where A R 8 × 8 is a positive definite matrix given by the following:
A = 0 I 2 × 2 0 0 0 0 I 2 × 2 0 0 0 0 I 2 × 2 K 1 K 2 K 3 K 4
The error between the exoskeleton model and the linear system Δ Q is assumed to be unknown. This error is estimated by a neural network, and the estimation Δ Q ^ has the form of Equation (11).
Therefore, using the neural network estimation, the actual control law τ is as follows:
τ = v τ Δ Q ^ ,
where v τ is given in Equation (19) and Δ Q ^ is the estimation given by (11). Substituting the actual control law in Equation (15), we have the following:
q ( 4 ) = v τ + Δ Q Δ Q ^
Since the output of the neural network is Δ Q = W T ϕ ( v ) + ϵ and the estimation is Δ Q ^ = W ^ T ϕ ( v ) , then the estimation error Δ Q ˜ is as follows:
Δ Q ˜ = Δ Q Δ Q ^ = W T ϕ ( v ) + ϵ W ^ T ϕ ( v ) = W ˜ T ϕ ( v ) + ϵ ,
where the estimation error of the weights W ˜ is as follows:
W ˜ = W W ^ R N × 2
Therefore, Equation (26) becomes the following:
q ( 4 ) = v τ + W ˜ T ϕ ( v ) + ϵ
Substituting v τ from Equation (19) in the above equation, the error dynamics are as follows:
q ( 4 ) q d ( 4 ) = K 1 ( q q d ) K 2 ( q ˙ q ˙ d ) K 3 ( q ¨ q ¨ d ) K 4 ( q ( 3 ) q d ( 3 ) ) + W ˜ T ϕ ( v ) + ϵ
Using Δ Q ˜ = W ˜ T ϕ ( v ) + ϵ , the error vector (21), and its derivative (22), the error dynamics in matrix form are as follows:
e ˙ = A e + B Δ Q ˜ ,
where A has the form of (24) and B = [ 0 2 × 2 0 2 × 2 0 2 × 2 I 2 × 2 ] T R 8 × 2 .

3.2. Stability Analysis

The stability analysis of the closed-loop system is developed in the present section.
Theorem 1. 
Consider the closed-loop system (31), the control law (25), and the adapting law as W ^ ˙ = γ ϕ ( x ) e T P B based on [24]. Then, for any bounded initial conditions, the closed-loop system signals remain bounded and the steady-state tracking error converges to a neighborhood of the origin by choosing the parameters adequately as in Section 3.3.
Proof. 
To perform the stability analysis of the closed-loop system, the following Lyapunov function candidate is proposed:
V ( e , W ˜ ) = 1 2 e T P e + 1 2 γ t r ( W ˜ T W ˜ ) ,
where γ is a positive constant, W ˜ denotes the weight estimation error given in (28), and P R 8 × 8 is a positive definite symmetric matrix that satisfies the following Lyapunov equation:
A T P + P A = Q ,
with Q 0 R 8 × 8 and A given by (24). Differentiating (32), we obtain V ˙ ( e , W ˜ ) as follows:
V ˙ ( e , W ˜ ) = 1 2 e T P e ˙ + e ˙ T P e + 1 γ t r ( W ˜ T W ˜ ˙ )
Substituting the tracking error dynamics from Equation (31), we have the following:
V ˙ ( e , W ˜ ) = 1 2 e T P ( A e + B Δ Q ˜ ) + ( A e + B Δ Q ˜ ) T P e + 1 γ t r ( W ˜ T W ˜ ˙ ) = 1 2 e T P A e + e T P B Δ Q ˜ + e T A T P e + Δ Q ˜ B T P e + 1 γ t r ( W ˜ T W ˜ ˙ ) = 1 2 e T ( P A + A T P ) e + e T P B Δ Q ˜ + Δ Q ˜ B T P e + 1 γ t r ( W ˜ T W ˜ ˙ ) = 1 2 e T Q e + e T P B ( W ˜ T ϕ ( v ) + ϵ ) + 1 γ t r ( W ˜ T W ˜ ˙ ) = 1 2 e T Q e + e T P B W ˜ T ϕ ( v ) + e T P b ϵ + 1 γ t r ( W ˜ T W ˜ ˙ )
Using the trace property, we have e T P B W ˜ T ϕ ( v ) = t r ( W ˜ T ϕ ( v ) e T P B ) . Additionally, W ˜ = W W ^ ; then W ˜ ˙ = W ^ ˙ , and V ˙ ( e , W ˜ ) is as follows:
V ˙ ( e , W ˜ ) = 1 2 e T Q e t r ( W ˜ T ϕ ( v ) e T P B ) + e T P B ϵ + 1 γ t r ( W ˜ T W ˜ ˙ ) = 1 2 e T Q e t r ( W ˜ T ϕ ( v ) e T P B + 1 γ W ˜ T W ˜ ˙ ) + e T P B ϵ = 1 2 e T Q e t r W ˜ T ϕ ( v ) e T P B 1 γ W ^ ˙ + e T P B ϵ
Since our aim is that ϕ ( v ) e T P B 1 γ W ^ ˙ = 0 , then the updating law of the weights W ^ ˙ is as follows:
W ^ ˙ = γ ϕ ( v ) e T P B
Substituting (37) into (36), then the following holds:
V ˙ ( e , W ˜ ) = 1 2 e T Q e + e T P B ϵ
We consider λ m i n ( Q ) e 2 e T Q e λ m a x ( Q ) e 2 and the bounds B = I 2 × 2 F = t r ( I T I ) 1 / 2 = 2 with the subindex F that represents the Frobenius norm of a matrix, P λ m a x ( P ) and ϵ < ϵ N ; then we have the following:
V ˙ ( e , W ˜ ) 1 2 λ m i n ( Q ) e 2 + e λ m a x ( P ) 2 ϵ N = 1 2 e λ m i n ( Q ) e 2 2 ϵ N λ m a x ( P )
If λ m i n ( Q ) e 2 2 ϵ N λ m a x ( P ) 0 holds, then V ˙ 0 . Therefore, the ultimate bound for e is [24] as follows:
e = 2 2 ϵ N λ m a x ( P ) λ m i n ( Q )

3.3. Control Gains Tuning

The auxiliary linear dynamics of Equation (13) are proposed, and an unknown function Δ Q is estimated by the neural network to compensate the nonlinear part of the model and external disturbances. Likewise, the control law consists of a linear controller v τ and the compensation of the neural network, as shown in Equation (25). Therefore, the gains of the linear controller in Equation (19) should be tuned. The gains k i with i = 1 , 2 , 3 , 4 are selected such that the roots of the characteristic polynomial π c ( s ) are in the left semi-plane of the complex plane.
π c ( s ) = s 4 + k 4 s 3 + k 3 s 2 + k 2 s + k 1
We assign the characteristic polynomial to the following 4th-order polynomial with the four roots in 1 , as follows:
π d = ( s + 1 ) 4 = s 4 + 4 s 3 + 6 s 2 + 4 s + 1
Then, the gains are k 4 = 4 , k 3 = 6 , k 2 = 4 k 1 = 1 . This dynamic behavior guaranties an over-damping response.
The following pseudocode (Algorithm 1) describes the algorithm of the neural network:
Algorithm 1 Pseudocode of neural network algorithm
  • Initialize the weight matrix of the neural network as W ^ ( 0 ) = o n e s ( N , 2 ) since W ^ R N × 2
  • For j = 1 to N, compute (Equation (10)):
    ϕ j ( v j ) = a j 1 + e b j v j c j
    to construct the vector ϕ R N × 1
  • Compute the output of the neuron as Δ Q ^ = W ^ T ϕ R 2 × 1 (Equation (11))
  • Compute the error vector e = [ q ˜ T q ˜ ˙ T q ˜ ¨ T q ˜ ( 3 ) T ] T R 8 × 1 (Equation (21))
  • Set the gain values for Equation (37) as γ = 2 , P = e y e ( 8 , 8 ) R 8 × 8 B = o n e s ( 8 , 2 ) R 8 × 2 .
  • Compute the derivative of the weight estimation matrix (Equation (37))
    W ^ ˙ = γ ϕ ( v ) e T P B
  • Integrate the derivative of the weight estimation matrix W ^ ˙ to obtain the weight estimation W ^
  • Output W ^ and Δ Q ^

4. Numerical and Experimental Results

In this section, we describe the numerical and experimental results obtained by using walking trajectories and implementing the proposed approach to the exoskeleton prototype.

4.1. Numerical Results

For the simulation tests, two gait cycle trajectories for the hip and knee joints are used. External disturbances are introduced as sinusoidal type signals, which represent different events in walking exercises such as muscle spasms and unmodeled dynamics of the exoskeleton. The vector of initial conditions for simulation is [ x 1 ( 0 ) , x 2 ( 0 ) , x 3 ( 0 ) , x 4 ( 0 ) ] T = [ 0.5 , 0.1 , 0.5 , 0.1 , 0.5 , 0.1 , 0.5 , 0.1 ] T , which represents the state initial conditions. The external disturbance vector is d = [ 0.2 s i n ( t ) 0.5 s i n ( t ) ] T . Moreover, the values of the neural network are a j = 2 , b j = 2 , and c j = 0 . The initial weight of the neural network is W ^ ( 0 ) = 0 , and considering a number of activation functions N = 5 with only one hidden layer, the weight matrix is as follows:
W = w 1 , 1 w 1 , 2 w 5 , 1 w 5 , 2
The parameter N was selected through an experimental test by reducing and increasing the number of activation functions. Then, the least number of activation functions with good performance was chosen, in this case, N = 5 . Figure 4 presents the trajectory tracking for two walking cycles, and the robot trajectory follows the desired trajectory with over-damping response. In Figure 4A, the red line is the desired trajectory of the hip q d 1 and the blue line is the trajectory of the hip robot joint q 1 . In Figure 4B, the red line is the desired trajectory of the knee q d 2 and the blue line is the trajectory of the knee robot joint q 2 generated by the implemented control law.
Figure 5A shows that the tracking error converges close to zero in approximately 2.5 s for both cases; the red line represents the hip joint error, and the blue line represents the knee joint error. Moreover, the means of absolute errors are q ˜ 1 = 0.93 ° and q ˜ 2 = 1.33 ° for the hip and the knee, respectively. Additionally, Figure 5B shows the unknown model parameter estimation by a neural network. Finally, Figure 6A and Figure 6B show the hip torque τ 1 and the knee torque τ 2 , respectively, generated by the control implemented in the exoskeleton.

4.2. Experimental Setup

In this subsection, we describe the experimental setup. The control algorithm was programmed on a myRIO card from National Instruments, using LabView software 2018.
The walking exoskeleton was tested on a healthy subject of 28 years with a weight of 57 kg and a height of 1.6 m with the security conditions and an informed consent letter.
The trajectories for the hip and the knee joint were generated using the approximation of 8th-order trigonometric polynomials as described in Section 3. The walking trajectories are from a healthy patient and are generated in real time by the card as the reference trajectories for the tracking controller.
The system is operated using a graphic interface elaborated in LabView. The graphic interface consists of modulators to vary the amplitude and frequency of the hip and knee trajectories during the walking. The maximum amplitude and frequency are limited to avoid patient injury. Additionally, the prototype movement is restricted by mechanical limits.
The test consists of the following steps: First, the subject is attached to the exoskeleton using a crane. Second, the system is turned on and the therapist sets the values of the walking trajectory in the graphical user interface. Then, the therapist begins the system by clicking on the start button in the graphical interface. The system executes the walking movement until the operator clicks on the end button, and the system finishes the cycle and goes to the initial position. A constant disturbance is added to the prototype to prove the robustness of the controller.
Three experimental tests were performed, so in experimental test A, the exoskeleton was used to follow the desired trajectory for gait rehabilitation (Figure 7A). In test B, the exoskeleton with a healthy user was evaluated to validate the robustness of the controller, considering the user as an unknown uncertainty (Figure 7B). In test C, external disturbances were added to the system and were rejected by the controller, so the robot followed the desired trajectory (Figure 7C).

4.3. Experimental Results

The experimental tests consist of specifying the setup of the walking cycle on the user interface following the steps in Section 4.2. The frequency and the amplitude of the trigonometric functions are established with a frequency of 0.15 Hz and an amplitude of 0.25. Moreover, the experimental values of the neural network are assigned as specified in the simulation section.
Figure 8 shows the tracking trajectories of four walking cycles executed by the lower-limb exoskeleton using the experimental setup. The solid black line is the desired trajectory q d 1 , the dotted red line shows the trajectory tracking of the hip robot joint q 1 for test 1 (only exoskeleton), the dotted blue line shows the trajectory tracking of the hip for test 2 (exoskeleton + user), and the dotted magenta line shows the trajectory tracking for test 3 (exoskeleton + user + disturbances).
Figure 9A and Figure 9B present the tracking error for the hip and the knee, and the means of absolute errors are q ˜ 1 = 1.27 ° and q ˜ 2 = 2.02 ° for the hip and the knee, respectively. The estimated unknown function is shown in Figure 9C. The experimental control inputs for the hip τ 1 and knee τ 2 are shown in Figure 10A,B; the units are the percentage of duty cycle of the pulse width modulation.
The links for the videos of the experimental tests are listed as follows:
  • In test A, the exoskeleton is used to follow the desired trajectory for gait rehabilitation (https://youtu.be/g3pz0G4x_aw accessed on 20 September 2024).
  • In test B, the exoskeleton with a healthy user is tested to validate the robustness of the controller considering the user as an unknown uncertainty (https://youtu.be/r4ciHMjtOJo accessed on 20 September 2024.
  • In test C, external disturbances are added to the system and are rejected by the controller, so the robot follows the desired trajectory (https://youtu.be/kjzzASk1yII accessed on 20 September 2024).

5. Conclusions

A control based on Artificial Neural Networks was applied to a system consisting of a robotic exoskeleton for gait rehabilitation. A leg exoskeleton with elastic joints was used to give comfort to the patient and to sense human intention of motion, as well as a linear feedback with an ANN estimator, which compensated the nonlinear dynamics and disturbances. This compensation provides robustness performance. The linear control guarantied an over-damping response for the tracking trajectory, which is important to smoothly follow the trajectories and avoid hurting the patient, during the therapy when external disturbances were present. The numerical simulation and the experimental results showed good performance of the tracking trajectory and validated the proposed control strategy. The robot joint trajectories converged to the desired trajectories, and tracking errors converged to a bounded region e 2 . 82 given by Equation (40), even when perturbations were included in the experimental setup.
The practical significance of the results is the relevance of the designed low-cost exoskeleton, which is capable of performing walking trajectories for people who have lower-limb disabilities. The test shows that the system is able to reject external disturbances and track the walking trajectories at a time response of less than 2 s. These characteristics guarantee security during rehabilitation even if there exist vibrations.
The further planned work is to sense the human intention of motion to activate the walking cycle and to perform active rehabilitation. Active rehabilitation means that the patient executes part of the effort to move the exoskeleton, and then the exoskeleton provides assistance as the patient needs it. Finally, we planned the validation of the system by a rehabilitation institute and the acceptance of a protocol to perform tests on injured patients.

Author Contributions

Conceptualization, K.I.E.-E. and R.L.-G.; methodology, K.I.E.-E., Y.R.-L., R.L.-G., S.S. and R.L.; software, K.I.E.-E. and Y.R.-L.; validation, S.S. and R.L.; formal analysis, K.I.E.-E.; investigation, K.I.E.-E., Y.R.-L., R.L.-G., S.S. and R.L.; resources, R.L.; data curation, K.I.E.-E.; writing—original draft preparation, K.I.E.-E., Y.R.-L. and R.L.-G.; writing—review and editing, R.L.-G., S.S. and R.L.; visualization, R.L.-G.; supervision, R.L.; project administration, S.S. and R.L.; funding acquisition, R.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Lee, T.; Kim, I.; Baek, Y.S. Design of a 2DoF Ankle Exoskeleton with a Polycentric Structure and a Bi-Directional Tendon-Driven Actuator Controlled Using a PID Neural Network. Actuators 2021, 10, 9. [Google Scholar] [CrossRef]
  2. Dao, Q.T.; Dinh, V.V.; Trinh, M.C.; Tran, V.C.; Nguyen, V.L.; Duong, M.D.; Bui, N.T. Nonlinear Extended Observer-Based ADRC for a Lower-Limb PAM-Based Exoskeleton. Actuators 2022, 11, 369. [Google Scholar] [CrossRef]
  3. Ran, M.; Li, J.; Xie, L. A new extended state observer for uncertain nonlinear systems. Automatica 2021, 131, 109772. [Google Scholar] [CrossRef]
  4. Fareh, R.; Khadraoui, S.; Abdallah, M.Y.; Baziyad, M.; Bettayeb, M. Active disturbance rejection control for robotic systems: A review. Mechatronics 2021, 80, 102671. [Google Scholar] [CrossRef]
  5. Aole, S.; Elamvazuthi, I.; Waghmare, L.; Patre, B.; Meriaudeau, F. Improved Active Disturbance Rejection Control for Trajectory Tracking Control of Lower Limb Robotic Rehabilitation Exoskeleton. Sensors 2020, 20, 3681. [Google Scholar] [CrossRef]
  6. Chen, C.F.; Du, Z.J.; He, L.; Wang, J.Q.; Wu, D.M.; Dong, W. Active Disturbance Rejection with Fast Terminal Sliding Mode Control for a Lower Limb Exoskeleton in Swing Phase. IEEE Access 2019, 7, 72343–72357. [Google Scholar] [CrossRef]
  7. Awad, N.; Humaidi, A.; Sabah, A. Fractional proportional derivative-based active disturbance rejection control of knee exoskeleton device for rehabilitation care. Indones. J. Electr. Eng. Comput. Sci. 2022, 28, 1405–1413. [Google Scholar] [CrossRef]
  8. Guerrero-Castellanos, J.; Rifaï, H.; Arnez-Paniagua, V.; Linares-Flores, J.; Saynes-Torres, L.; Mohammed, S. Robust Active Disturbance Rejection Control via Control Lyapunov Functions: Application to Actuated-Ankle–Foot-Orthosis. Control Eng. Pract. 2018, 80, 49–60. [Google Scholar] [CrossRef]
  9. Zhao, L.; Liu, X.; Zhang, J.; Wang, T. Angle tracking control for double-joint dexterous hand systems based on a piecewise extended state observer. Control Eng. Pract. 2021, 110, 104754. [Google Scholar] [CrossRef]
  10. San-Miguel, A.; Puig, V.; Alenyà, G. Disturbance observer-based LPV feedback control of a N-DoF robotic manipulator including compliance through gain shifting. Control Eng. Pract. 2021, 115, 104887. [Google Scholar] [CrossRef]
  11. Aole, S.; Elamvazuthi, I.; Waghmare, L.; Patre, B.; Bhaskarwar, T.; Meriaudeau, F.; Su, S. Active Disturbance Rejection Control Based Sinusoidal Trajectory Tracking for an Upper Limb Robotic Rehabilitation Exoskeleton. Appl. Sci. 2022, 12, 1287. [Google Scholar] [CrossRef]
  12. Li, H.; Cheng, L.; Li, Z.; Xue, W. Active Disturbance Rejection Control for a Fluid-Driven Hand Rehabilitation Device. IEEE/ASME Trans. Mechatron. 2021, 26, 841–853. [Google Scholar] [CrossRef]
  13. Hagan, M.T.; Demuth, H.B.; Beale, M.H.; Jesús, O.D. Neural Network Design, 2nd ed.; Martin Hagan, Oklahoma UE; 2014; Available online: https://books.google.com.mx/books?id=4EW9oQEACAAJ (accessed on 29 July 2024).
  14. Moura Coelho, R.; Gouveia, J.; Botto, M.A.; Krebs, H.I.; Martins, J. Real-time walking gait terrain classification from foot-mounted Inertial Measurement Unit using Convolutional Long Short-Term Memory neural network. Expert Syst. Appl. 2022, 203, 117306. [Google Scholar] [CrossRef]
  15. Fuentes-Alvarez, R.; Hernandez, J.H.; Matehuala-Moran, I.; Alfaro-Ponce, M.; Lopez-Gutierrez, R.; Salazar, S.; Lozano, R. Assistive robotic exoskeleton using recurrent neural networks for decision taking for the robust trajectory tracking. Expert Syst. Appl. 2022, 193, 116482. [Google Scholar] [CrossRef]
  16. Yu, W.; Rosen, J. Neural PID Control of Robot Manipulators with Application to an Upper Limb Exoskeleton. IEEE Trans. Cybern. 2013, 43, 673–684. [Google Scholar] [CrossRef]
  17. Liu, G.; Sun, N.; Liang, D.; Chen, Y.; Yang, T.; Fang, Y. Neural network-based adaptive command filtering control for pneumatic artificial muscle robots with input uncertainties. Control Eng. Pract. 2022, 118, 104960. [Google Scholar] [CrossRef]
  18. Aly, A.; The Vu, M.; El-Sousy, F.; Hsia, K.H.; Alotaibi, A.; Mousa, G.; Le, D.N.; Mobayen, S. Adaptive Neural Network-Based Fixed-Time Tracking Controller for Disabilities Exoskeleton Wheelchair Robotic System. Mathematics 2022, 10, 3853. [Google Scholar] [CrossRef]
  19. Al-Waeli, K.H.; Ramli, R.; Haris, S.M.; Zulkoffli, Z.B.; Amiri, M.S. Offline ANN-PID Controller Tuning on a Multi-Joints Lower Limb Exoskeleton for Gait Rehabilitation. IEEE Access 2021, 9, 107360–107374. [Google Scholar] [CrossRef]
  20. Han, S.; Wang, H.; Tian, Y.; Christov, N. Time-delay estimation based computed torque control with robust adaptive RBF neural network compensator for a rehabilitation exoskeleton. ISA Trans. 2020, 97, 171–181. [Google Scholar] [CrossRef] [PubMed]
  21. Hao, Z.; Liu, K.; Wei, Q. Adaptive Neural Network Control of Lower Limb Exoskeleton Robots Using Disturbance Observer. In Proceedings of the 2020 5th International Conference on Advanced Robotics and Mechatronics (ICARM), Shenzhen, China, 18–21 December 2020; pp. 246–251. [Google Scholar] [CrossRef]
  22. Peng, Z.; Luo, R.; Huang, R.; Hu, J.; Shi, K.; Cheng, H.; Ghosh, B.K. Data-Driven Reinforcement Learning for Walking Assistance Control of a Lower Limb Exoskeleton with Hemiplegic Patients. In Proceedings of the 2020 IEEE International Conference on Robotics and Automation (ICRA), Paris, France, 31 May–31 August 2020; pp. 9065–9071. [Google Scholar] [CrossRef]
  23. Rosales-Luengas, Y.; Espinosa-Espejel, K.I.; Lopéz-Gutiérrez, R.; Salazar, S.; Lozano, R. Lower Limb Exoskeleton for Rehabilitation with Flexible Joints and Movement Routines Commanded by Electromyography and Baropodometry Sensors. Sensors 2023, 23, 5252. [Google Scholar] [CrossRef] [PubMed]
  24. Liu, J. Radial Basis Function (RBF) Neural Network Control for Mechanical Systems; Springer: Berlin/Heidelberg, Germany, 2013; pp. 134–139. [Google Scholar]
Figure 1. Final prototype and general distribution of actuators.
Figure 1. Final prototype and general distribution of actuators.
Sensors 24 06546 g001
Figure 2. Schematic diagram and frame orientation of the lower-limb exoskeleton.
Figure 2. Schematic diagram and frame orientation of the lower-limb exoskeleton.
Sensors 24 06546 g002
Figure 3. Scheme of the neural network-based control.
Figure 3. Scheme of the neural network-based control.
Sensors 24 06546 g003
Figure 4. Tracking trajectories (two walking cycles) for (A) hip joint and (B) knee joint.
Figure 4. Tracking trajectories (two walking cycles) for (A) hip joint and (B) knee joint.
Sensors 24 06546 g004
Figure 5. (A) Tracking error and (B) unknown parameter estimation.
Figure 5. (A) Tracking error and (B) unknown parameter estimation.
Sensors 24 06546 g005
Figure 6. (A) Control input for hip joint τ 1 ; (B) control input for knee joint τ 2 .
Figure 6. (A) Control input for hip joint τ 1 ; (B) control input for knee joint τ 2 .
Sensors 24 06546 g006
Figure 7. Experimental tests: (A) only exoskeleton, (B) exoskeleton + user, and (C) exoskeleton + user + external disturbance (blue arrow).
Figure 7. Experimental tests: (A) only exoskeleton, (B) exoskeleton + user, and (C) exoskeleton + user + external disturbance (blue arrow).
Sensors 24 06546 g007
Figure 8. Experimental tracking trajectories (four walking cycles) for (A) hip joint and (B) knee joint.
Figure 8. Experimental tracking trajectories (four walking cycles) for (A) hip joint and (B) knee joint.
Sensors 24 06546 g008
Figure 9. (A,B) Experimental tracking error for hip and knee joint and (C) unknown parameter estimation.
Figure 9. (A,B) Experimental tracking error for hip and knee joint and (C) unknown parameter estimation.
Sensors 24 06546 g009
Figure 10. (A) Experimental control input for hip τ 1 ; (B) experimental control input for knee τ 2 .
Figure 10. (A) Experimental control input for hip τ 1 ; (B) experimental control input for knee τ 2 .
Sensors 24 06546 g010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Espinosa-Espejel, K.I.; Rosales-Luengas, Y.; Salazar, S.; Lopéz-Gutiérrez, R.; Lozano, R. Active Disturbance Rejection Control via Neural Networks for a Lower-Limb Exoskeleton. Sensors 2024, 24, 6546. https://doi.org/10.3390/s24206546

AMA Style

Espinosa-Espejel KI, Rosales-Luengas Y, Salazar S, Lopéz-Gutiérrez R, Lozano R. Active Disturbance Rejection Control via Neural Networks for a Lower-Limb Exoskeleton. Sensors. 2024; 24(20):6546. https://doi.org/10.3390/s24206546

Chicago/Turabian Style

Espinosa-Espejel, Karina I., Yukio Rosales-Luengas, Sergio Salazar, Ricardo Lopéz-Gutiérrez, and Rogelio Lozano. 2024. "Active Disturbance Rejection Control via Neural Networks for a Lower-Limb Exoskeleton" Sensors 24, no. 20: 6546. https://doi.org/10.3390/s24206546

APA Style

Espinosa-Espejel, K. I., Rosales-Luengas, Y., Salazar, S., Lopéz-Gutiérrez, R., & Lozano, R. (2024). Active Disturbance Rejection Control via Neural Networks for a Lower-Limb Exoskeleton. Sensors, 24(20), 6546. https://doi.org/10.3390/s24206546

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop