Command Filter Backstepping Sliding Model Control for Lower-Limb Exoskeleton

A command filter adaptive fuzzy backstepping control strategy is proposed for lower-limb assisting exoskeleton. Firstly, the humanrobotmodel is established by taking the humanbody as a passive part, and a coupling torque is introduced to describe the interaction between the exoskeleton and human leg. Then, Vicon motion capture system is employed to obtain the reference trajectory. For the purpose of obviating the “explosion of complexity” in conventional backstepping, a second-order command filter is introduced into the sliding mode control strategy. The fuzzy logic systems (FLSs) are also applied to handle with the chattering problem by estimating the uncertainties and disturbances. Furthermore, the stability of the closed-loop system is proved based on the Lyapunov theory. Finally, simulation results are presented to illustrate the effectiveness of the control strategy.


Introduction
Recently, the exoskeleton is increasingly used for powerassisting in industrial [1,2], medical [3][4][5], and military [6,7] areas. The human operator provides locomotion intention and a small muscle strength while the robot can help human complete the desired action with a suitable torque through replicating the operator's movements [2]. The robot in such applications has been actively researched since 1990s, and some of them have been applied for various environments [8]. However, there are challenging problems in exoskeleton research such as the establishment of coupling model of human-robot system and designing of control strategy.
From the modeling point of view, it has been proved that some methods used now are effective such as Newton-Euler equations and Lagrange dynamics [9,10]. Based on aforementioned work, series of improvements have been proposed to increase the accuracy of human-robot system model. The inverse dynamic model is established and the singular points are avoided using damped least squares in [2]. In order to approximate the actual situation of human legs, a variety of musculoskeletal models are developed [11,12]. The comparison between Hill-type and proportional model for human muscle is illustrated in [12], and Hill-type models are proved to be more appropriate.
Nowadays, there have been several published papers on control strategies of human-robot cooperative control [13][14][15][16][17][18]. For instance, a sensitivity amplification control is proposed in [19] which could track the desired trajectory by minimizing the interaction torque. In [20], a robust sliding mode controller is proposed to guarantee the stability in disturbance situation, and the boundary layer is introduced to reduce the chattering problem of sliding model control. Additionally, an adaptive sliding model control based on state observer is proposed in [21], which could update the controller parameters online to improve the safety of the system.
However, most of the control methods mentioned have limitations. In modeling, taking the human effects as disturbance [22], widely used in exoskeleton researches, is unreasonable. Human body is a part of the system obviously. Additional, the system has high-order features when taking the human body as a passive part [19]. Hence, the control strategies above are not available because of the "explosion of complexity" [23].
In this paper, a command filtered backstepping sliding model control equipped with FLSs is proposed based on 2 Mathematical Problems in Engineering a human-robot model. Compared with the analogous literature, the main contributions are summarized as follows: (i) The proposed human-robot model is established by taking the human leg as a passive part and the coupling torque between human and robot is introduced. Compared with the model built in [19], this paper converts the transfer to the state space form and introduces uncertainties and disturbances which are more actual and complicated.
(ii) For the purpose of testing the controller performance, an experiment is implanted to obtain the actual trajectory of the hip-joint. Compared with the sine curve, the upper bound increases rapidly after a few derivative operations, which may cause the system uncontrollable.
(iii) The command filter backstepping sliding model control is proposed. By using the command filter, the analytical derivate is unnecessary and the "explosion of complexity" in the controller design process is avoided [24][25][26]. Besides, the problem mentioned in (ii) is solved. Furthermore, the FLSs are used to approximate the uncertainties and disturbances and provide real-time compensations for the system.
The paper is organized as follows: in Section 2, the high-order human-robot model is introduced based on the linearized models of human leg and exoskeleton. Then the reference trajectory of hip-joint is obtained through experiments. In Section 3, considering the demands of controller design, a model in state space form with uncertainties and disturbances is obtained. Based on that work, the humanrobot controller is designed with a command filter adaptive backstepping sliding mode and, by utilizing the Lyapunov theory, closed-loop system stability is analyzed. Simulation results are discussed in Section 4 and the conclusions are provided in Section 5. Figure 1, the human-robot system can be expressed as a person wearing a lower-limb exoskeleton which provides a back support to assist hip-joint motion. By passing the leg's gravity to the waist, the muscle effort needed for human walking could be reduced to a low level.

Dynamics of the Human-Robot System. As shown in
In order to describe the system, an elementary model is used in this paper which consists of linearized onedegree-of-freedom (1-DOF) models for the human leg and the exoskeleton. In the process of modeling the physical interaction between the human and exoskeleton, a coupling torque, expressed as combination of a linear spring and a damper, is introduced.
Then, the ideal dynamics of the human-robot system are given as follows [19]:  where ℎ , ℎ , ℎ are, respectively, the moment of inertia, joint damping coefficient, and joint stiffness coefficient of the human leg; ℎ ( ) is the hip-joint angle; ℎ ( ) is the net muscle torque acting on the joint. , , are, respectively, the moment of inertia, joint damping coefficient, and joint stiffness coefficient of the exoskeleton; ( ) is the exoskeleton joint angle; ( ) is the actuator torque.
And the coupling torque is defined as follows: where , are, respectively, the equivalent damping coefficient and stiffness coefficient of the interaction torque; is the coupling torque.

Trajectory Generation.
The main way that the exoskeleton helps human complete the locomotion is tracking the human gait cycle. So a reasonable desired position trajectory is an essential factor for testing the model and controller. An experimenter (girl aged 25 years with mass of 52 kg and stature of 165 cm) volunteered to participate in the gait experiment with Vicon motion capture system, device provided by National Research Center for Rehabilitation Technical Aids ( Figure 2). Special trackers are fixed at the particular marks on the experimenter which can be captured by cameras distributed  Hip-joint degree in reasonable location in experiment space. The information from different cameras are combined, and the actual hipjoint angle signals are obtained through data processing. From the considerable data obtained, the most reliable, representative, and authentic data is selected and reorganized.
To ensure the smoothness of the trajectory, the Fourier series is introduced to describe the actual gait cycle. Note that the human locomotion satisfies the smooth characteristics. Mathematical expression of the desired trajectory can be expressed as follows [14]: where 0 is the initial value of 2 ; and are the sine and cosine amplitudes of Fourier series; is the fundamental frequency; is the harmonic order.
The parameters of the trajectory are obtained by the curve fit toolbox of MATLAB. In order to approximate the actual curve and simplify the calculation process, fundamental frequency to third harmonics, that is, = 3 in (3), of the Fourier series is used in this paper. The parameters are illustrated in Table 1.
According to (3) and the parameters given in Table 1, a reliable hip-joint trajectory (expressed as 1 ) can be obtained and shown in Figure 3.
The control objectives for the human-robot system are illustrated as follows: (i) An adaptive controller for high-order human-robot system is designed, such that the position of human leg 1 and exoskeleton 3 can track the actual trajectory 1 obtained from experiment. (ii) The prescribed output tracking error̃1 = 1 − 1 is always bounded. Besides when uncertainties and disturbances exist in the system, the tracking error can converge to a neighborhood of the origin in a short time.

System Description and Control Strategy.
For the exoskeleton, the human leg is a passive part and fulfills the locomotion with the interaction torque between the human and robot. For the dynamics shown in (1), when replacing with (2), the system is given by Note that all the states and the torques are time-varying variables and time flags are omitted for convenience.
Taking the lump uncertainties, parametric/unmodeled uncertainties as well as the external disturbances, into account, (4) can be transformed into state space aṡ 4 =̇are the state vector of the system; the lump uncertainties in human-robot system are Note that the specific parameters of human leg are hard to be measured and the actuator of the exoskeleton includes mechanical errors which cannot be described precisely. Hence, the uncertainties of the system existed and are inevitable.

Mathematical Problems in Engineering
Considering that the system has high-order features, a backstepping sliding model method is introduced to solve the complex problem with a recursive form [27]. To avoid the "explosion of complexity" in the controller design process, a second-order nonlinear command filter is employed in this paper. The command filter ensures that the desired command and its derivative satisfy the same magnitude and rate constrains [24].
In order to handle the system chattering problem caused by the uncertainties and disturbances, fuzzy logic systems (FLSs) are equipped to estimate the upper bounds of the lump uncertainties. FLSs provide real-time compensations for the human-robot system to reduce the switching items of the sliding model control.
Being equipped with second-order command filter and FLSs, a backstepping sliding mode control strategy for the human-robot system with uncertainties and disturbances is proposed in this paper.

Basic Assumptions and FLS.
Some reasonable and useful assumptions are given at first which ensure the stability of the system. Assumption 1. There exist constants 1 , 2 > 0, such that the inequality |Δ 1 | < 1 , |Δ 2 | < 2 holds.
The second step is defuzzification. Center average defuzzification operator is applied in this paper which can be expressed aŝ( where = max ∈ ( ), and ( ), are the membership functions.
Define the fuzzy basis vector as The optimal parameter * can be defined by * = arg min Lemma 3 (Wang [28]). For ∀ ( ), which is continuous function and defined over a compact Ω, for any a constant ∈ + , there exist an FLS and a parameter * such that

Controller Design.
Consider the characteristics of the human-robot system, a backstepping sliding model control with second-order command filter is proposed in this paper. The controller design process is shown in this section. The output tracking errors and compensated tracking errors of the subsystems are defined, respectively, as where and represent the system states and filtered commands of the th subsystem, respectively.
The ( = 1, 2, 3, 4) signals can be obtained bẏ with (0) = 0. (14) is used to achieve the filtering value which are designed to compensate the errors caused by the command filters. Note that they can be computed with integrating processes to avoid the differential operations.

Remark 4. Equation
For the purpose of eliminating the "explosion of complexity," a second-order nonlinear command filter is designed to calculate ,̇. So the command filter is shown as follows: where is the natural frequency of the filter and typically satisfies > , = (1, 2, 3, 4) to ensure the tracking accuracy. ∈ (0, 1] is the damping ratio of the filter system. Mathematical Problems in Engineering 5 , ( = 2, 3, 4) are the virtual control signals and 1 is the desired trajectory.
The filter initial conditions are (0) = 0,̇(0) = 0. Every command filter is designed to compute the filtered commands without differential operation. Furthermore, will track by choosing the suitable parameters. In order to find the optimal parameters of the FLSs, the adaptive laws are chosen for̂1 and̂2 as follows: where 1 , 2 are the positive adaptive coefficients.
Considering the compensating errors and the closed-loop dynamics, the virtual law can be defined as follows: wherê1,̂2 are the estimated values of optimal parameters and 1 , 2 are the basis vectors of FLSs. > 0, ( = 1, 2, 3, 4) are control gains specified by the designer. 1 , 2 are constants that ensure that the inequality 1 > | 1 |, 2 > | 2 | holds. The sats( ) is the switching function that satisfies

Stability Analysis
Theorem 5. For the system illustrated in (5), there exist a range of values for the gains 1 , 2 , 3 , 4 , and the adaptive coefficients 1 , 2 , such that the tracking error̃1 and compensation errors ( = 1, 2, 3, 4) can converge to zero with the compensations provided by FLSs.
Proof. The tracking error̃and the compensated tracking error are given first.
Due to the fact that =̃− for = 1, 2, 3, 4, the dynamics of compensated tracking error can be expressed as follows: Mathematical Problems in Engineering Remark 6. Note the fact that the command will directly output to control plant, so the conclusion is easily obtained that = .
The block diagram of controller design is shown in Figure 4.

Simulation Results
In this section, a simulation of an 1-DOF lower-limb exoskeleton is established. All the parameters of the 1-DOF exoskeleton are illustrated in Table 2.
The PARM is short for parameter, and the notation is omitted for convenience through the paper as long as special notation is not required.
Remark 7. All the parameters are cited in [19]. The parameters, just for calculating, are obtained from real experiment and useful for controller simulation.
The lump uncertainties are chosen as follows: To ensure the stability of the system, the specific parameters of the controller are tuned in a trial-and-error procedure and shown in Table 3.
The simulation results are given in Figures 5-10. It can be seen that the controller designed can guarantee the uncertainties and external disturbances.
x 1d x 3d x 4d x 3d x 2d x 1 x 3 x 4 x 2 x 1cẋ 4ċ x 3cẋ 2c x 1c x 3c x 2c x 4c  Position tracking (rad) Time (s) Figure 5: Position tracking of human-robot system. Figures 5 and 6 illustrate the position tracking and the tracking errors of the human-robot system. In Figure 5, the desired trajectory 1 is represented as solid line, the human leg position 1 is shown as dash-dotted line, and the position of exoskeleton 3 is described as dashed line. Just as it is shown, the controlled plant can steadily track the hip curve in a satisfactory way. Figure 6 shows the errors of the trajectory tracking. It is easy to get that the proposed command filtered fuzzy adaptive backstepping controller could make the errors However there are still differences between the two errors because of the filtered errors. Figures 7 and 8 show the virtual control signals as solid line and the filtered commands , which are produced by passing the virtual control signals through the secondorder filter, as the dashed line. As expected, the control signals of all the subsystems are smooth and bounded. Moreover, the virtual control signals and the corresponding filtered commands satisfy the same magnitude, rate, and bandwidth constrains, which is different from the traditional first-order linear filter. The closed-loop tracking errors could converge to a tiny range around the origin in a very short time (about 0.1 s) just as shown in Figure 9. x 4d x 3c x 3d x 2c x 2d x 1c   are represented as dashed lines. It can be seen that the estimated value (dash) reaching the real value (solid) in less than 0.2 s. Therefrom, upper bounds of the switch terms in sliding model control could be much smaller because of the compensations provided by FLSs.
All the signals remain bounded in a reasonable range during the process. Obviously, the proposed control strategy with command filters and FLSs can be a suitable method for lower-limb exoskeleton.

Conclusion
A human-robot cooperative control strategy based on a convincing high-order model is proposed for a lower-limb assisting exoskeleton. A second-order command filter backstepping method is employed to determine the time derivatives of virtual control signals without differential operations. The FLSs are used to approximate the uncertainties and disturbances and compensate the system timely. In addition, the stability of the system is proved based on the Lyapunov theory. Finally, simulation results are presented to verify the effectiveness of the proposed command filter adaptive fuzzy control strategy.
Future work will focus on the performance of the control strategy in actual experiment and the filtering errors should be proved to converge rigorously. Besides, further research on the state constrain control for high-order nonlinear systems [30] is also needed.

Conflicts of Interest
The authors declare that they have no conflicts of interest.