Study on the Control Method of Knee Joint Human–Exoskeleton Interactive System

The advantages of exoskeletons based on the Bowden cable include being lightweight and flexible, thus being convenient in assisting humans. However, the performance of an exoskeleton is limited by the structure and human–exoskeleton interaction, which is analyzed from the established mathematical model of the human–exoskeleton system. In order to improve the auxiliary accuracy, corresponding control methods are proposed. The disturbance observer is designed to compensate for disturbances and parameter perturbations in the inner loop. The human–exoskeleton interaction feedforward model is integrated into the admittance control, which overcomes the limitation of the force loading caused by the friction of the Bowden cable and the change in stiffness of the human–exoskeleton interaction. Furthermore, an angle prediction method using the encoder as the signal source is designed to reduce the disturbance of the force loading caused by human motion. Finally, the effectiveness of the design method proposed in this paper is verified through experiments.


Introduction
Exoskeletons have been extensively studied in recent years because they can enhance the athletic ability of humans. According to different structures, the lower-extremity exoskeletons can be divided into rigid exoskeletons [1][2][3] and flexible exoskeletons [4][5][6]. The mechanism of a rigid exoskeleton is parallel to the human limbs, which is convenient for weight bearing and rehabilitation training. The corresponding control methods [7,8] have also been studied. However, due to structural and inertial limitations, rigid exoskeletons are unsatisfactory in assisting healthy individuals to walk. Therefore, a wearable, flexible exoskeleton is proposed. The flexible exoskeleton mainly comprises textiles and the Bowden cable [4]. The textile is the bearing point of the sheath of the Bowden cable. The Bowden cable [9,10] has advantages such as long transmission distance and flexible layout. Therefore, the actuator can be concentrated near the center of body weight. The integration and flexibility of the exoskeleton can be effectively enhanced.
The human-exoskeleton system (HESY) contains two inputs and one output [11]: a force command, active human motion, and force response. In [12], the HESY's parameter perturbation and active movement are regarded as disturbance. In other words, the HESY is regarded as a single-input and single-output system. A disturbance observer (DOB) is designed to compensate for the disturbance caused by parameter perturbation and active motion based on the dynamic model. To further improve the dynamic performance of the force loop, a sliding mode controller based on the nominal model [12,13] is proposed. control to compensate for Bowden cable friction disturbance and the HEI's stiffness change on the force loading. Finally, an angle prediction method that adopts the encoder as the signal source is proposed, effectively reducing the influence of human motion and signal acquisition delay on force loading.
The structure of this paper is as follows. The mathematical model of the HESY is established in Section 2. The corresponding controller design is presented in Section 3. The experimental verification of the theory in this paper is shown in Section 4. The conclusions are presented in Section 5.

Mathematical Modeling of the HESY
The farther the distance between the load and the center of body weight, the more metabolism the body consumes. Therefore, the power supply, controller, and servo motor are integrated near the center of body weight. The fixation mechanism of the knee joint [31] is shown in Figure 1, which is driven by a servo motor through a Bowden cable. Two lightweight carbon fiber boards (LCFB) and hinge mechanisms can realize rotation in the sagittal and coronal planes. The servo motor is set to speed mode to improve the dynamic response of the inner loop. The human-exoskeleton interaction feedforward model (HEIF) is added to the admittance control to compensate for Bowden cable friction disturbance and the HEI's stiffness change on the force loading. Finally, an angle prediction method that adopts the encoder as the signal source is proposed, effectively reducing the influence of human motion and signal acquisition delay on force loading. The structure of this paper is as follows. The mathematical model of the HESY is established in Section 2. The corresponding controller design is presented in Section 3. The experimental verification of the theory in this paper is shown in Section 4. The conclusions are presented in Section 5.

Mathematical Modeling of the HESY
The farther the distance between the load and the center of body weight, the more metabolism the body consumes. Therefore, the power supply, controller, and servo motor are integrated near the center of body weight. The fixation mechanism of the knee joint [31] is shown in Figure 1, which is driven by a servo motor through a Bowden cable. Two lightweight carbon fiber boards (LCFB) and hinge mechanisms can realize rotation in the sagittal and coronal planes. The servo motor is set to speed mode to improve the dynamic response of the inner loop.  The schematic diagram of the HESY is shown in Figure 2. According to Figure 2, the mathematical model of the HESY can be established.
The human-exoskeleton interaction feedforward model (HEIF) is added to the admittance control to compensate for Bowden cable friction disturbance and the HEI's stiffness change on the force loading. Finally, an angle prediction method that adopts the encoder as the signal source is proposed, effectively reducing the influence of human motion and signal acquisition delay on force loading.
The structure of this paper is as follows. The mathematical model of the HESY is established in Section 2. The corresponding controller design is presented in Section 3 The experimental verification of the theory in this paper is shown in Section 4. The conclusions are presented in Section 5.

Mathematical Modeling of the HESY
The farther the distance between the load and the center of body weight, the more metabolism the body consumes. Therefore, the power supply, controller, and servo motor are integrated near the center of body weight. The fixation mechanism of the knee joint [31] is shown in Figure 1, which is driven by a servo motor through a Bowden cable. Two lightweight carbon fiber boards (LCFB) and hinge mechanisms can realize rotation in the sagittal and coronal planes. The servo motor is set to speed mode to improve the dynamic response of the inner loop.  τ m is the output torque of the current loop of the servo motor. J e and b e are the equivalent moment of inertia and equivalent damping of the servo motor, respectively. θ is the angular displacement of the input shaft of the servo motor reducer. K and l F are the comprehensive stiffness and elastic deformation of the transmission mechanism between the servo motor and the knee joint fixation mechanism, respectively. θ h is the angle of the knee joint encoder. I t = I h + I l , where I h and I l are the equivalent inertia moment of the knee joint and fixation mechanism, respectively. b t = b h + b l , where b h and b l are the equivalent damping of the knee joint and fixation mechanism, respectively. k t = k h + k l , where k h and k l are the equivalent stiffness of the knee joint and fixation mechanism, respectively. Figure 3 shows the displacement relationship of the HESY, as shown below where θ m is the angular displacement of the output shaft of the servo motor reducer, r e is the radius of the sheave at the output shaft of the reducer, ∆l is the inelastic deformation during the movement, θ h1 and θ h2 are the angles of the knee joint and the deformation of the HEI, respectively, r h is the radius of the sheave at the knee joint, and n is the reduction ratio of the servo motor.
tween the servo motor and the knee joint fixation mechanism, respectively. h θ is the angle of the knee joint encoder. = + t h l I I I , where h I and l I are the equivalent inertia moment of the knee joint and fixation mechanism, respectively. = + are the equivalent damping of the knee joint and fixation mechanism, respectively. = + t h l k k k , where h k and l k are the equivalent stiffness of the knee joint and fixation mechanism, respectively. Figure 3 shows the displacement relationship of the HESY, as shown below where s is the Laplace operator, L is the inductance coefficient, i K is the current feedback coefficient, R is the resistance coefficient, d k is the amplification factor of the speed loop, u k is the amplification factor, e C is the back EMF coefficient, m k is the torque coefficient, and v k is the feedback coefficient of the speed loop. The interaction force is as follows: The mathematical equation of the servo motor speed mode is as follows: where s is the Laplace operator, L is the inductance coefficient, K i is the current feedback coefficient, R is the resistance coefficient, k d is the amplification factor of the speed loop, k u is the amplification factor, C e is the back EMF coefficient, k m is the torque coefficient, and k v is the feedback coefficient of the speed loop. The interaction force is as follows: where F is the interaction force. The knee joint flexion angle in the sagittal plane can reach 60 • . Therefore, only the dynamic equation in the sagittal plane is considered. The dynamic equation of the lower limb can be simplified as where τ h is the muscle torque of the human, τ is the interaction torque, and τ g is the gravity torque. Combined with Equations (1)-(5), the block diagram of the inner loop [31] can be obtained as follows. It can be seen from Figure 4 that there are three factors affecting the force loading accuracy. First, the speed loop of the servo motor is disturbed by Bowden cable friction and the interaction force. Due to the poor anti-interference ability of the integrated motor, its dynamic performance is affected. Second, τ h cannot be accurately measured. Moreover, the impedance characteristic of the knee joint changes with the gait, which enhances the stability of human motion. Therefore, the mathematical modeling of Equation (5) is difficult. Last, the HEI can be deformed during the movement. In this paper, the controller is designed for dealing with the problems mentioned above.
where τ h is the muscle torque of the human, τ is the interaction torque, and τ g is the gravity torque. Combined with Equations (1)-(5), the block diagram of the inner loop [31] can be obtained as follows. It can be seen from Figure 4 that there are three factors affecting the force loading accuracy. First, the speed loop of the servo motor is disturbed by Bowden cable friction and the interaction force. Due to the poor anti-interference ability of the integrated motor, its dynamic performance is affected. Second, τ h cannot be accurately measured. Moreover, the impedance characteristic of the knee joint changes with the gait, which enhances the stability of human motion. Therefore, the mathematical modeling of Equation (5) is difficult. Last, the HEI can be deformed during the movement. In this paper, the controller is designed for dealing with the problems mentioned above.

Design of the Inner Loop Controller
Ensuring the anti-interference ability and stable dynamic performance of the inner loop are the prerequisites for realizing accurate assistance. Mostly, the DOB can observe and compensate for external disturbance and parameter perturbation. A DOB is adopted to improve the dynamic performance of the inner loop. The design of the DOB requires high modeling accuracy. Therefore, the motor speed loop is adopted instead of the knee joint angular velocity closed-loop. The transfer function of the speed loop can be rewritten as

Design of the Inner Loop Controller
Ensuring the anti-interference ability and stable dynamic performance of the inner loop are the prerequisites for realizing accurate assistance. Mostly, the DOB can observe and compensate for external disturbance and parameter perturbation. A DOB is adopted to improve the dynamic performance of the inner loop. The design of the DOB requires high modeling accuracy. Therefore, the motor speed loop is adopted instead of the knee joint angular velocity closed-loop. The transfer function of the speed loop can be rewritten aŝ whereĜ 2 (s) is the nominal model of the speed loop. The meanings of the coefficients in Equation (6) are the same as those in Equation (2). The block diagram of the inner loop based on the DOB is shown in Figure 5. Q(s) is the filter and ξ is the high-frequency noise of the speed signal.
Sensors 2022, 22, x FOR PEER REVIEW 6 of 18 The speed loop can be regarded as a three-input one-output system and the transfer function is as follows: Fr   The speed loop can be regarded as a three-input one-output system and the transfer function is as follows: . In the low-frequency range, Q(s) = 1. On the contrary, Q(s) = 0 in the high-frequency range. Therefore, the transfer function of the speed loop can be rewritten as It can be seen from Equation (11) that the dynamic performance of the speed loop in the low-frequency band is equivalent to that in the designed nominal model. Both the force disturbance and internal parameter perturbation can be effectively compensated for. However, the bandwidth of the filter is limited by the stability conditions of the DOB. The stability conditions are described as [31] where ∆ 2 (jω) is the modeling error of the speed loop. It should be noted that Equation (12) is essential for designing a DOB based on a servo motor instead of knee joints. The inertia and damping of the speed loop can be compensated based on Equation (11). The formula is as follows: where C f 1 and C f 2 are proportional coefficients, respectively.

Design of the HEIF
Combining Figure 5 and Equation (11), it can be seen that the human motion and the deformation of the HEI can affect the assisting effect. It is necessary to measure the comprehensive stiffness to compensate for the deformation of the HEI. The posture is relatively stable during the initial landing. The motor rotates at a constant speed. During this process, the relationship between the linear displacement of the servo motor and the interaction force is recorded. The measured result is shown in Figure 6. According to the measured data, curve fitting is performed, and the formula is shown follows: where C f d1 and C f d2 are the fitting coefficients, respectively. It can be seen from Figure 6 that the HEI has not only nonlinear elastic deformation but also inelastic deformation during the assisting process, which is an inevitable structural feature of the wearable, flexible exoskeleton.
where 1 f C and 2 f C are proportional coefficients, respectively.

Design of the HEIF
Combining Figure 5 and Equation (11), it can be seen that the human motion and the deformation of the HEI can affect the assisting effect. It is necessary to measure the comprehensive stiffness to compensate for the deformation of the HEI. The posture is relatively stable during the initial landing. The motor rotates at a constant speed. During this process, the relationship between the linear displacement of the servo motor and the interaction force is recorded. The measured result is shown in Figure 6. According to the measured data, curve fitting is performed, and the formula is shown follows: where 1 fd C and 2 fd C are the fitting coefficients, respectively. It can be seen from Figure  6 that the HEI has not only nonlinear elastic deformation but also inelastic deformation during the assisting process, which is an inevitable structural feature of the wearable, flexible exoskeleton. Therefore, to ensure the force loading performance, nonlinear elastic deformation needs to be compensated for as a feedforward term, which is shown below:  Therefore, to ensure the force loading performance, nonlinear elastic deformation needs to be compensated for as a feedforward term, which is shown below: However, the inelastic deformation cannot be compensated by Equation (15). Moreover, the interaction force needs to be generated by the displacement of the Bowden cable. The inaccurate speed response in the inner loop can cause displacement deviation. Therefore, iterative force control is proposed to compensate for the inelastic deformation of the HEI and the displacement error caused by the speed loop. The iterative speed command of the next gait cycle is corrected by detecting the interaction force error in the current gait cycle. In order to prevent accidental factors from causing excessive force errors and ensure the safety of the system, the maximum single correction displacement is limited to 0.5 mm. The error threshold and speed correction command are as follows .
where δ i is the speed correction coefficient, e F i−1 is the maximum force error value of the previous gait cycle, e F 0 is the force error threshold value, .
θ ci is the current gait speed correction value, and T ei is the average value of the previous two gait cycles.

Design of the Angle Prediction Method
It can be seen from Equation (1) and Figure 5 that the human motion will interfere with the interaction force. The usual method is to collect the motion information of the knee joint and compensate for it to the servo motor. However, due to factors such as signal delay and phase lag, motion compensation is not timely. The joint angle during walking is a periodic signal. Therefore, parallel oscillators are introduced in this section, and each oscillator learns a particular frequency component of the input signal. In other words, Fourier decomposition is performed on the input signal. Only the fundamental frequency [32] is learned to simplify the design, the others being multiples of it. Each phase of the adaptive oscillator (AO) is as follows: where φ i (t) is the phase of the i-th oscillator, ω(t) is the oscillator's fundamental frequency, is the response value of the oscillator, and v is the proportional coefficient. The gait in this paper is not entirely consistent. Therefore, the natural frequency of the oscillator is set as a state variable, as shown below: .
According to the Fourier series, the output is as follows: where α i (t) is the amplitude of the output signal of the i-th oscillator, and η is the integral gain. A kernel-based nonlinear filter is introduced based on Equation (20) to achieve angle prediction, as shown below:θ where Ψ i (φ 1 (t)) is the i-th kernel function, h is the width of the kernel function, c i = 2π N i, N is the number of kernel functions, and w i is the weighting factor, which can be obtained through iteration as follows: where P is the inverse covariance matrix and λ is a constant. The input variable in Equation (21) is converted to φ 1 (t). φ 1 (t) will change uniformly in the range 0 − 2π when Equation (21) converges. Therefore, φ 1 (t) + ∆ φ can be used to predict human movement, which is shown as follows:θ Combined with the angle relationship in Figure 3, the feedforward compensation term for human motion is This paper has compensated for the disturbance of the HEI and active motion to the assistance. The relationship between the inner loop and the outer loop is established by admittance control, as follows: where I d e and b d e are virtual inertia and damping, respectively, and F d is the force command. Through the above analysis, the system control block diagram is shown in Figure 7.

Test Setup
The knee exoskeleton in this paper is shown in Figure 8. The designed knee joint exoskeleton is composed of the power supply, control unit, knee joint fixing structure, and sensors. The control system in this paper consists of a host computer and a target computer. The PC can be used as the host computer, through which the program is

Test Setup
The knee exoskeleton in this paper is shown in Figure 8. The designed knee joint exoskeleton is composed of the power supply, control unit, knee joint fixing structure, and sensors. The control system in this paper consists of a host computer and a target computer. The PC can be used as the host computer, through which the program is written, compiled, and downloaded to the target computer. The upper computer can also be used to control and change the variables in the target computer and receive the sensors' values from the target computer. To improve the integration, Sbrio-9636 is adopted as the lower computer. The servo motor is a Hai Tai brushless DC motor made in Yiwu, China. The power, rated torque, rated speed, and peak torque is 200 W, 7.6 Nm, 250 RPM, and 15 Nm, respectively. The transmission ratio is 6:1. The human may be damaged when the exoskeleton severely hinders human movement during walking. Therefore, the emergency stop button is connected to the DI channel of Sbrio-9636. Subjects can press the emergency stop button in their hands when they feel inconsistent with the movement of the exoskeleton. Sbrio-9636 controls the servo motor to the initial position in time by judging that it receives the high-level signal of the emergency stop button.

Test Setup
The knee exoskeleton in this paper is shown in Figure 8. The designed knee joint exoskeleton is composed of the power supply, control unit, knee joint fixing structure, and sensors. The control system in this paper consists of a host computer and a target computer. The PC can be used as the host computer, through which the program is written, compiled, and downloaded to the target computer. The upper computer can also be used to control and change the variables in the target computer and receive the sensors' values from the target computer. To improve the integration, Sbrio-9636 is adopted as the lower computer. The servo motor is a Hai Tai brushless DC motor made in Yiwu, China. The power, rated torque, rated speed, and peak torque is 200 W, 7.6 Nm, 250 RPM, and 15 Nm, respectively. The transmission ratio is 6:1. The human may be damaged when the exoskeleton severely hinders human movement during walking. Therefore, the emergency stop button is connected to the DI channel of Sbrio-9636. Subjects can press the emergency stop button in their hands when they feel inconsistent with the movement of the exoskeleton. Sbrio-9636 controls the servo motor to the initial position in time by judging that it receives the high-level signal of the emergency stop button. The exoskeleton in this paper has the following characteristics. First, the antiinterference ability of the inner loop is poor under the disturbance of the interaction force and the friction of the Bowden cable. Second, the stiffness of the HEI changes significantly The exoskeleton in this paper has the following characteristics. First, the anti-interference ability of the inner loop is poor under the disturbance of the interaction force and the friction of the Bowden cable. Second, the stiffness of the HEI changes significantly during the assisting process. Third, the dynamic performance of the outer loop is limited by the stiffness change and the frictional hysteresis of the Bowden cable. Finally, the force loading is disturbed by human movement. Combined with the above situation, the experimental ideas in this section are as follows.
According to Figure 7, the HESY contains two inputs and one output: force command, active human motion, and force response. In order to verify the influence of different input signals on the output, the human body is first fixed in a particular posture. In other words, the force loading experiment is performed while keeping the active motion signal of the human body at 0, and the constant force loading is performed during walking to verify the disturbance caused by the active movement of the human. Finally, the control method proposed in this paper is used to assist the human during walking to verify its effectiveness. Informed consent for inclusion was given by three subjects before the research, which was conducted according to the Declaration of Helsinki. The protocol of Registering Clinical Trials (ChiECRCT20200319) was approved by the Chinese Ethics Committee. The force command [33] in this section is the statistical curve of the knee joint in the flexion direction when the human is walking at a speed of 4.5 km/h.

PID Control in Both Inner and Outer Loop
The posture is stable during the initial landing. Therefore, the initial landing phase is selected as the test pose. When the PID is used in both the inner and outer loops, the force load response curve and the speed response curve are shown in Figure 9. The shaded areas are the standard deviations of the force response and speed response. It can be seen from Figure 9a that there is a 60 ms delay in force response. The extension speed of the knee joint at the end of the swing phase should be reduced under the torque in the flexion direction. The protective flexion torque is rapidly increased in the initial landing phase to prevent knee hyperextension.
in the flexion direction when the human is walking at a speed of 4.5 km/h.

PID Control in Both Inner and Outer Loop
The posture is stable during the initial landing. Therefore, the initial landing phase is selected as the test pose. When the PID is used in both the inner and outer loops, the force load response curve and the speed response curve are shown in Figure 9. The shaded areas are the standard deviations of the force response and speed response. It can be seen from Figure 9a that there is a 60 ms delay in force response. The extension speed of the knee joint at the end of the swing phase should be reduced under the torque in the flexion direction. The protective flexion torque is rapidly increased in the initial landing phase to prevent knee hyperextension.
However, in 0-10% of the gait cycle, the assistant hardly responds when the force command changes rapidly. On the one hand, the extension of the knee joint at the end of the swing phase is hindered by inappropriate assistance. On the other hand, in order to counteract the improper assistance, the torque in the extension direction of the human is forced to increase, which is not conducive to reducing the metabolism. It can be seen from Figure 9b that the tracking performance of the inner loop is temporarily satisfactory. Therefore, the dynamic performance of force loading is limited by the outer loop. Both phase lag and amplitude attenuation in the outer loop can seriously affect the assistance to the human. Therefore, it is necessary to study the control method further to ensure the force loading performance.  However, in 0-10% of the gait cycle, the assistant hardly responds when the force command changes rapidly. On the one hand, the extension of the knee joint at the end of the swing phase is hindered by inappropriate assistance. On the other hand, in order to counteract the improper assistance, the torque in the extension direction of the human is forced to increase, which is not conducive to reducing the metabolism. It can be seen from Figure 9b that the tracking performance of the inner loop is temporarily satisfactory. Therefore, the dynamic performance of force loading is limited by the outer loop. Both phase lag and amplitude attenuation in the outer loop can seriously affect the assistance to the human. Therefore, it is necessary to study the control method further to ensure the force loading performance.

The HEIF in the Outer Loop and the DOB in the Inner Loop
The frequency sweep method is used to establish the nominal model of the actuator without load. The amplitude of the frequency sweep signal is 45 r/min, and the frequency range is 1-7 Hz. The response curve is shown in Figure 10. It can be seen that the response curve of the nominal model (NM) is basically the same as the dynamic response curve (WAM) of the actuator, which can verify the validity of the nominal model. Furthermore, the DOB can be designed based on nominal models.
Through experiments in previous studies, we have verified that the parameter perturbation of the inner loop and the poor anti-interference ability can deteriorate the dynamic performance of the outer loop [31]. Furthermore, it is concluded that the inner loop based on the DOB can improve the dynamic performance potential of the outer loop. Therefore, this paper does not repeat the DOB-based inner loop experiment. Instead, it verifies the effectiveness of the control method through combining the HEIF and the admittance control based on the DOB-based inner loop.
Through experiments in previous studies, we have verified that the parameter perturbation of the inner loop and the poor anti-interference ability can deteriorate the dynamic performance of the outer loop [31]. Furthermore, it is concluded that the inner loop based on the DOB can improve the dynamic performance potential of the outer loop. Therefore, this paper does not repeat the DOB-based inner loop experiment. Instead, it verifies the effectiveness of the control method through combining the HEIF and the admittance control based on the DOB-based inner loop. When the human is stationary, the force response and speed response curves are shown in Figure 11. The shaded areas are the standard deviations of the force response and speed response. It can be seen from Figure 11b that the inner loop still maintains adequate dynamic performance and anti-interference ability when the amplitude and the frequency of the speed command increase. The outer loop makes full use of the relationship between force and deformation to overcome the limitation of force loading caused by changes in the stiffness of the HEI. In other words, the contradiction between the stability under the maximum stiffness and the rapidity under the minimum stiffness is effectively avoided. On the other hand, the HEIF also avoids the phase lag of the force response caused by force closed-loop control and Bowden cable friction. It can be seen from Figure 11a that the phase lag and amplitude attenuation of the force response are compensated effectively when the human body is stationary. The effectiveness of the HEIF based on the DOB in the inner loop in this paper is verified. Then, we verify the interference of human motion on auxiliary accuracy and the effectiveness of the angle prediction method based on the controller of the HEIF in the outer loop and DOB in the inner loop. When the human is stationary, the force response and speed response curves are shown in Figure 11. The shaded areas are the standard deviations of the force response and speed response. It can be seen from Figure 11b that the inner loop still maintains adequate dynamic performance and anti-interference ability when the amplitude and the frequency of the speed command increase. The outer loop makes full use of the relationship between force and deformation to overcome the limitation of force loading caused by changes in the stiffness of the HEI. In other words, the contradiction between the stability under the maximum stiffness and the rapidity under the minimum stiffness is effectively avoided. On the other hand, the HEIF also avoids the phase lag of the force response caused by force closed-loop control and Bowden cable friction. It can be seen from Figure 11a that the phase lag and amplitude attenuation of the force response are compensated effectively when the human body is stationary. The effectiveness of the HEIF based on the DOB in the inner loop in this paper is verified. Then, we verify the interference of human motion on auxiliary accuracy and the effectiveness of the angle prediction method based on the controller of the HEIF in the outer loop and DOB in the inner loop. Transparency is one of the important indicators of the wearable exoskeleton. The Bowden cable can only be driven in one direction. Therefore, to verify the transparency of the exoskeleton in this paper, the force command and walking speed are set to 5 N and 4.5 km/h, respectively. The corresponding force response curve and knee displacement curve are shown in Figure 12. Transparency is one of the important indicators of the wearable exoskeleton. The Bowden cable can only be driven in one direction. Therefore, to verify the transparency of the exoskeleton in this paper, the force command and walking speed are set to 5 N and 4.5 km/h, respectively. The corresponding force response curve and knee displacement curve are shown in Figure 12.
dynamic curve of force response; (b) the dynamic curve of speed response.
Transparency is one of the important indicators of the wearable exoskeleton. The Bowden cable can only be driven in one direction. Therefore, to verify the transparency of the exoskeleton in this paper, the force command and walking speed are set to 5 N and 4.5 km/h, respectively. The corresponding force response curve and knee displacement curve are shown in Figure 12. It can be seen that the force response is 0 N when the knee is flexed. The transparency in the flexion direction is satisfactory. However, the assistance is disturbed by the phase lag of the actuator and the friction of the Bowden cable, when the knee joint is extended. However, the assistance is disturbed by the phase lag of the actuator and the friction of the Bowden cable when the knee joint is extended. Moreover, with the increase in the extension speed, the disturbance force also increases. The disturbance force reaches 35 N at the maximum extension speed, which severely reduces the transparency of the exoskeleton.
Further, the force command is set as a statistical force curve, and the walking speed is 4.5 km/h. The motion information of the knee joint is compensated to the servo motor. The corresponding force response curve and the speed response curve are shown in Figure 13. The shaded areas are the standard deviations of the force response and speed response. According to Figure 13, the force response at the initial landing is only 15.9 N. Moreover, the force response is significantly disturbed in 70-90% of the gait cycle. Therefore, it is necessary to study the angle prediction method. It can be seen that the force response is 0 N when the knee is flexed. The transparency in the flexion direction is satisfactory. However, the assistance is disturbed by the phase lag of the actuator and the friction of the Bowden cable, when the knee joint is extended. However, the assistance is disturbed by the phase lag of the actuator and the friction of the Bowden cable when the knee joint is extended. Moreover, with the increase in the extension speed, the disturbance force also increases. The disturbance force reaches 35 N at the maximum extension speed, which severely reduces the transparency of the exoskeleton.
Further, the force command is set as a statistical force curve, and the walking speed is 4.5 km/h. The motion information of the knee joint is compensated to the servo motor. The corresponding force response curve and the speed response curve are shown in Figure 13. The shaded areas are the standard deviations of the force response and speed response. According to Figure 13, the force response at the initial landing is only 15.9 N. Moreover, the force response is significantly disturbed in 70-90% of the gait cycle. Therefore, it is necessary to study the angle prediction method.

The HEIF-AO in the Outer Loop and the DOB in the Inner Loop
A constant force of 30 N is applied to the human throughout the gait in order to intuitively verify the effectiveness of the angle prediction method while not excessively affecting human movement. The joint angle predicted by AO during walking is shown in Figure 14. It can be seen that the prediction of the joint angle can be achieved through AO. According to Figure 15, when the motion information collected by the encoder is directly compensated to the servo motor, the force response will periodically change between 10 and 60 N with the motion. According to Figures 12 and 15, the range of interference force is correspondingly increased from 35 N to 50 N when the constant force is set to 5 N and 30 N, respectively. The importance of angle prediction compensation is confirmed. The force response varies between 20 and 40 N when the angle predicted by AO is input to the

The HEIF-AO in the Outer Loop and the DOB in the Inner Loop
A constant force of 30 N is applied to the human throughout the gait in order to intuitively verify the effectiveness of the angle prediction method while not excessively affecting human movement. The joint angle predicted by AO during walking is shown in Figure 14. It can be seen that the prediction of the joint angle can be achieved through AO. According to Figure 15, when the motion information collected by the encoder is directly compensated to the servo motor, the force response will periodically change between 10 and 60 N with the motion. According to Figures 12 and 15, the range of interference force is correspondingly increased from 35 N to 50 N when the constant force is set to 5 N and 30 N, respectively. The importance of angle prediction compensation is confirmed. The force response varies between 20 and 40 N when the angle predicted by AO is input to the servo motor. The experimental results show that the AO under the constant force loading can effectively reduce the disturbance by more than 60%.
A constant force of 30 N is applied to the human throughout the gait in order to intuitively verify the effectiveness of the angle prediction method while not excessively affecting human movement. The joint angle predicted by AO during walking is shown in Figure 14. It can be seen that the prediction of the joint angle can be achieved through AO. According to Figure 15, when the motion information collected by the encoder is directly compensated to the servo motor, the force response will periodically change between 10 and 60 N with the motion. According to Figures 12 and 15, the range of interference force is correspondingly increased from 35 N to 50 N when the constant force is set to 5 N and 30 N, respectively. The importance of angle prediction compensation is confirmed. The force response varies between 20 and 40 N when the angle predicted by AO is input to the servo motor. The experimental results show that the AO under the constant force loading can effectively reduce the disturbance by more than 60%.  Finally, the validity of the control method proposed in this paper is verified when the force command and the walking speed are set to a statistical force curve and 4.5 km/h, respectively. The fundamental frequency and the phase of the i-th oscillator curves obtained by AO are shown in Figure 16 and Figure 17, respectively. It can be seen from Figure 16 that the fundamental frequency is constantly changing between 4.7 and 5.3 rad/s. According to Figure 17, the RMSE of the phase of the i-th oscillator is 1.252 0.04 ± , which shows that the motion of the knee joint is not strictly a periodic signal. Correspondingly, this is also why the local prediction bias is relatively large. Finally, the validity of the control method proposed in this paper is verified when the force command and the walking speed are set to a statistical force curve and 4.5 km/h, respectively. The fundamental frequency and the phase of the i-th oscillator curves obtained by AO are shown in Figures 16 and 17, respectively. It can be seen from Figure 16 that the fundamental frequency is constantly changing between 4.7 and 5.3 rad/s. According to Figure 17, the RMSE of the phase of the i-th oscillator is 1.252 ± 0.04, which shows that the motion of the knee joint is not strictly a periodic signal. Correspondingly, this is also why the local prediction bias is relatively large. the force command and the walking speed are set to a statistical force curve and 4.5 km/h, respectively. The fundamental frequency and the phase of the i-th oscillator curves obtained by AO are shown in Figure 16 and Figure 17, respectively. It can be seen from Figure 16 that the fundamental frequency is constantly changing between 4.7 and 5.3 rad/s. According to Figure 17, the RMSE of the phase of the i-th oscillator is 1.252 0.04 ± , which shows that the motion of the knee joint is not strictly a periodic signal. Correspondingly, this is also why the local prediction bias is relatively large.   the force command and the walking speed are set to a statistical force curve and 4.5 km/h, respectively. The fundamental frequency and the phase of the i-th oscillator curves obtained by AO are shown in Figure 16 and Figure 17, respectively. It can be seen from Figure 16 that the fundamental frequency is constantly changing between 4.7 and 5.3 rad/s. According to Figure 17, the RMSE of the phase of the i-th oscillator is 1.252 0.04 ± , which shows that the motion of the knee joint is not strictly a periodic signal. Correspondingly, this is also why the local prediction bias is relatively large.   The force response curve and the speed response curve are shown in Figure 18. The shaded areas are the standard deviations of the force response and speed response. It can be seen from Figure 18 that the mean value of assistance in the initial landing phase increases from 15.9 N to 32.2 N. Moreover, the assisted advancement phenomenon almost disappears in 70-90% of the gait cycle. The overall force tracking has almost no phase lag and amplitude attenuation, which can meet the assistance accuracy. The experimental results verify the validity of the theory proposed in this paper. The force response curve and the speed response curve are shown in Figure 18. The shaded areas are the standard deviations of the force response and speed response. It can be seen from Figure 18 that the mean value of assistance in the initial landing phase increases from 15.9 N to 32.2 N. Moreover, the assisted advancement phenomenon almost disappears in 70-90% of the gait cycle. The overall force tracking has almost no phase lag and amplitude attenuation, which can meet the assistance accuracy. The experimental results verify the validity of the theory proposed in this paper.

Conclusions
This paper analyzes the influence of Bowden cable friction, HEI, and human motion on the assistance by mathematical modeling. The corresponding control strategy is proposed based on the speed inner loop and admittance control. The experimental results

Conclusions
This paper analyzes the influence of Bowden cable friction, HEI, and human motion on the assistance by mathematical modeling. The corresponding control strategy is proposed based on the speed inner loop and admittance control. The experimental results show that the HEIF based on the DOB inner loop can effectively compensate for the force response's phase lag and amplitude attenuation when the motion is ignored. However, in the initial landing phase of human walking, the average assistance is only 15.9 N. In 70%-90% of the gait cycle, the force response is significantly disturbed. The force response periodically changes between 10 and 60 N when the command is 30 N, which further verifies the disturbance of the motion to the assistance. When the predicted angle compensates for the active motion, the disturbance to the force response can be effectively reduced by more than 60%. The error of the mean assistance at the initial landing is reduced from 23.4 N to 7.1 N during walking. Moreover, in 70%-90% of the gait cycle, the disturbance to the force response is compensated. It can be proven that the theory in this paper can meet the needs of assisting knee joints.