Dynamic Parameter Identification of a Lower Extremity Exoskeleton Using RLS-PSO

The lower extremity exoskeleton is a device for auxiliary assistance of human movement. The interaction performance between the exoskeleton and the human is determined by the lower extremity exoskeleton’s controller. The performance of the controller is affected by the accuracy of the dynamic equation. Therefore, it is necessary to study the dynamic parameter identification of lower extremity exoskeleton. The existing dynamic parameter identification algorithms for lower extremity exoskeletons are generally based on Least Square (LS). There are some internal drawbacks, such as complicated experimental processes and low identification accuracy. A dynamic parameter identification algorithm based on Particle Swarm Optimization (PSO) with search space defined by Recursive Least Square (RLS) is developed in this investigation. The developed algorithm is named RLS-PSO. By defining the search space of PSO, RLS-PSO not only avoids the convergence of identified parameters to the local minima, but also improves the identification accuracy of exoskeleton dynamic parameters. Under the same experimental conditions, the identification accuracy of RLS-PSO, PSO and LS was quantitatively compared and analyzed. The results demonstrated that the identification accuracy of RLS-PSO is higher than that of LS and PSO.


Introduction
Over the past decades, advances in electromechanical technology have accelerated the development of lower extremity exoskeletons.Unless otherwise stated, the exoskeleton referred to in this article refers to the lower extremity exoskeleton.The exoskeleton is a device that interacts with the human.The main applications of exoskeletons are physical assistance and rehabilitation [1][2][3][4][5][6][7][8][9][10].During human-machine interaction, under no circumstances should the exoskeleton directly or indirectly cause injury to the human, either in the routine operation or in fault [11].Therefore, exoskeleton research should not be confined to mechanical design.It is also challenging and important to improve the performance of the exoskeleton's controller and to realize compliance control in the process of human-machine interaction [12][13][14].One of the effective methods to enhance the flexibility of the exoskeleton's controller is to build an accurate dynamic model [15][16][17][18][19].
Dynamic parameters are used to build the dynamic model of the exoskeleton.There are three main methods to obtain the dynamic parameters of the exoskeleton.First, obtaining the CAD data directly from the manufacturers.This method is not applicable to non-standard products such as exoskeletons.Second, obtaining the CAD data from 3D design software.This method ignores many factors such as the dynamics of hoses and wires, and the internal dynamics of actuators that affect the dynamics of the exoskeleton.Third, obtaining the dynamic parameters through dynamic parameter identification algorithms.This method has high accuracy and requires a series of experiments for data acquisition [20][21][22].The exoskeleton is made up of many non-standard parts, such as pipes and actuators [6][7][8][9][10], which cause the dynamic parameters of the exoskeleton to fluctuate during working.Obviously, it is difficult to obtain the dynamic parameters from the manufacturers or 3D design software.Therefore, parameter identification algorithms are expected to be an effective way to obtain more accurate dynamic parameters [23,24].
In order to obtain more accurate dynamic parameters, researchers are increasingly paying attention to the application of parameter identification algorithms for the exoskeleton.Researchers have carried out some preliminary research works.Targeting the swing phase [25,26], Justin et al. identified the dynamic parameters of Berkeley Exoskeleton (BLEEX) [21] by means of Least Square (LS) [27,28].The joint friction coefficients of BLEEX were identified by static experiments, and the inertial parameters [29] were identified by dynamic experiments.The method did not make full use of the non-correlation between the parameters and the linear relationship between the parameters and the parameterized dynamic equations [30,31].Therefore, the method only identified one parameter in each experiment.Designing a corresponding identification experiment for each parameter was necessary.The process of parameter identification experiment was very complicated, and the friction parameters on the hip joint lost identifiability because of the limited range of motion.The parameter identification accuracy of this method is low.Justin et al. did not compare the identification accuracy of LS with other algorithms.Lin Cong et al. also used LS for dynamic parameters identification of the exoskeleton [16].However, the identified parameters were only compared with the CAD data calculated by Solidworks.Due to the existence of a large number of non-standard parts [16], it is difficult to prove the accuracy of this method.Jonas [20] proposed a weighted LS method and verified it on an exoskeleton with flexible actuators.The accuracy of the algorithm was also compared with the CAD data of 3D design software.The CAD data could not offer the real value of all the parameters to be identified.Therefore, the accuracy of the method was also not verified.In conclusion, the current parameter identification algorithms for exoskeletons are either complicated, have low accuracy, or lack of quantitative accuracy verification.The identification algorithms are generally based on the LS method.The performance of other parameter identification algorithms in exoskeleton dynamic parameter identification has not been studied.However, Zafer [32] has demonstrated that Particle Swarm Optimization (PSO) has higher identification accuracy than LS on an RX-60 manipulator.Therefore, the introduction of the PSO is expected to improve the accuracy of identification of an exoskeleton's dynamic parameters.
PSO identifies the parameters through the intelligence of particle swarm and iterative optimization processes [33].Each dimension represents a parameter to be identified, and the ranges of parameters to be identified are often unknown.In order to avoid losing the real value of parameters, PSO optimization is often carried out in the infinite domain [34].If the number of iterations is insufficient or the factors are not set properly, the probability that PSO converges to the local minima increases.Taking the local minima as the identification values, the error of the exoskeleton dynamic model will increase.Therefore, it is necessary to reduce the search space of PSO from the infinite domain to the finite domain in advance, to reduce the number of potential local minima within the search space as much as possible.
The fluctuation of exoskeleton parameters provides a possibility to reduce the search space.The approximate range of each parameter can be identified by the recursive identification algorithm before accurate identification of PSO.Recursive least squares (RLS) can identify the fluctuation of parameters.Therefore, the search space of PSO can be defined by RLS.The accuracy of PSO is expected to be further improved with the finite search space.The existing identification algorithms of exoskeleton dynamic parameters are generally based on LS with low accuracy.In this paper, a PSO with search space defined by RLS is designed, named RLS-PSO.The RLS-PSO was verified on a hydraulically driven exoskeleton of the lower extremity.The accuracy of RLS-PSO was compared quantitatively with LS and PSO.

Methods
The dynamic equations of the exoskeleton were derived and converted into the form required by the parameter identification algorithm.The physical parameters to be identified were determined by the dynamic equations.According to the linearization requirements of RLS, the physical parameters were converted into inertial parameters [15,16,23,35].RLS-PSO was designed based on the linearized dynamic equations.

Parametric Dynamics Equation of Exoskeleton
The RLS-PSO was verified on the active power-assist exoskeleton (APAL) [36] as shown in Figure 1a, which consists of trunk, thigh, shank, and foot modules.As shown in Figure 1b, the trunk and left lower extremity were fixed on the pedestal during the experiment.The right lower extremity was suspended.The hip and knee joints can only swing in the sagittal plane.The forefoot of the exoskeleton was attached to the shank by the fastening band in Figure 1b.The ankle has zero degrees of freedom.The hip and knee joints are driven by hydraulic cylinders, and the exoskeleton is driven only by torque and gravity.The hydraulic actuator of hip joint A is fixed in the trunk and will not affect the dynamic parameters which are related to the thigh and hip joint.Hydraulic cylinder CD drives knee joint B and is fixed on AB.Through the four-linkage mechanism, ABCD, the linear motion of the CD is converted to the rotation of B. The stretching and swinging of the CD cause changes in the position of lower extremity's center of mass (COM) [36].Therefore, the dynamic parameters of the thigh and shank fluctuate.The fluctuation was used for defining the search space of PSO.The dynamic equations of the exoskeleton were derived and converted into the form required by the parameter identification algorithm.The physical parameters to be identified were determined by the dynamic equations.According to the linearization requirements of RLS, the physical parameters were converted into inertial parameters [15,16,23,35].RLS-PSO was designed based on the linearized dynamic equations.

Parametric Dynamics Equation of Exoskeleton
The RLS-PSO was verified on the active power-assist exoskeleton (APAL) [36] as shown in Figure 1a, which consists of trunk, thigh, shank, and foot modules.As shown in Figure 1b, the trunk and left lower extremity were fixed on the pedestal during the experiment.The right lower extremity was suspended.The hip and knee joints can only swing in the sagittal plane.The forefoot of the exoskeleton was attached to the shank by the fastening band in Figure 1b.The ankle has zero degrees of freedom.The hip and knee joints are driven by hydraulic cylinders, and the exoskeleton is driven only by torque and gravity.The hydraulic actuator of hip joint A is fixed in the trunk and will not affect the dynamic parameters which are related to the thigh and hip joint.In the sagittal plane, the right leg of the APAL exoskeleton was simplified into a two-linkage model [37], as shown in Figure 2.
Equation ( 1) is the dynamic equation of APAL's lower extremity.Equation (1) was derived by the Lagrange method.
where 1 q  and 2 q  represent the angular velocity of hip and knee joint, respectively, and 1 q  and 2 q  represent the angular acceleration of hip and knee joint, respectively.h T and k T represent the In the sagittal plane, the right leg of the APAL exoskeleton was simplified into a two-linkage model [37], as shown in Figure 2.
Equation ( 1) is the dynamic equation of APAL's lower extremity.Equation (1) was derived by the Lagrange method.
Appl.Sci.2019, 9, 324 4 of 17 where .q 1 and .q 2 represent the angular velocity of hip and knee joint, respectively, and .. q 1 and .. q 2 represent the angular acceleration of hip and knee joint, respectively.T h and T k represent the hip and knee drive torque required by APAL to complete the trajectory of q 1 and q 2 , respectively.When the lower extremity of APAL is vertically downward, q 1 and q 2 are zero.The counterclockwise rotation is positive.
Appl.Sci.2018, 8, x FOR PEER REVIEW 4 of 18 hip and knee drive torque required by APAL to complete the trajectory of 1 q and 2 q , respectively.
When the lower extremity of APAL is vertically downward, 1 q and 2 q are zero.The counterclockwise rotation is positive.
The angle of hip joint q2: The angle of knee joint mt: The mass of thigh ms: The mass of shank It: The moment of inertia of thigh Is: The moment of inertia of shank Lt: The length of OK Ls: The length of KL hGt: The distance from the COM of thigh to HK hGs: The distance from the COM of shank to KL LGt: The length of O 1 O projected on HK LGs: The length of O 2 K projected on KL Considering the friction on the hip and knee joint, friction torques T were introduced.Equation ( 2) is the real joint torque 1 T and 2 T .
where Since there are couplings between the physical parameters t 2) is nonlinear to the above physical parameters.In order to meet the linearization requirements of RLS [38,39], Equation ( 2) was parameterized to Equation (3).The physical parameters Considering the friction on the hip and knee joint, friction torques T f 1 and T f 2 were introduced.Equation ( 2) is the real joint torque T 1 and T 2 . where f c1 and f v1 represent the friction parameters of the hip joint, f c2 and f v2 represent the friction parameter of the knee joint.
Since there are couplings between the physical parameters I t , I s , m t , m s , L t , L s , L Gt , L Gs , h Gt , h Gs in Equation (2), Equation ( 2) is nonlinear to the above physical parameters.In order to meet the linearization requirements of RLS [38,39], Equation ( 2) was parameterized to Equation (3).The physical parameters I t , I s , m t , m s , L t , L s , L Gt , L Gs , h Gt , h Gs in Equation ( 2) were parameterized into a set of inertial parameters M xt , M xs , M yt , M ys , J t , J s that are capable of fully expressing the dynamic characteristics and independent of each other [16,22,23,40].Equation ( 3) is linear for the parameters to-be-identified.
where 3) is the linearization dynamic equation of APAL, where M xt , M xs , M yt , M ys , J t , J s , f c1 , f c2 , f v1 , f v2 are the parameters to-be-identified.

Design of RLS-PSO Parameter Identification Algorithm
The design of the RLS-PSO was based on Equation ( 3).The RLS-PSO consists of two parts.First, defining the search space for each parameter.Second, identifying parameters within the search space.The search spaces of the parameters were defined by RLS.The accurate identification of the parameters was implemented by the PSO within the search space.

Establishment of the RLS
As mentioned earlier, the parameters to-be-identified fluctuated due to the movement of the hydraulic actuators.The RLS can identify the time-variant parameters.Therefore, the RLS was adopted to define the fluctuation ranges of M xt , M xs , M yt , M ys , J t , J s , f c1 , f c2 , f v1 , f v2 .The fluctuation ranges of M xt , M xs , M yt , M ys , J t , J s , f c1 , f c2 , f v1 , f v2 form a 10-dimension search space.In order to facilitate computer programming and calculation, Equation (3) was converted into matrix form as Equation (4). where In Equation ( 4), H h and H k represent motion information vectors.X represents the parameter vector.While the hip and knee joints are swinging, the hip joint generates a torque to counteract the interference caused by the swing of knee joint to complete the target trajectory.Similarly, the knee joint also generates a torque that counteracts the interference caused by the swing of hip joint.Equation ( 5) was derived from Equation (4).Equation ( 5) is the equation of the LS.
X was identified in each iteration of RLS to define the fluctuation range of X.Therefore, Equation ( 5) was converted into a recursive form as Equation ( 6).Equation ( 6) is the RLS.
where kg represents the number of iterations of RLS.The detailed derivation processes of Equation ( 6) and the principle of RLS can be referred to [41].The fluctuation ranges of X given by the RLS is the search space of PSO.

PSO with a Finite Search Space
PSO is an algorithm based on particle swarm intelligence.PSO is mainly used for parameter identification and optimization.In each iteration, PSO adjusts the direction and distance of the flight based on the best solution of the particle's own motion trajectories and the best solution of the whole particles' motion trajectories.Through continuous iterations, PSO finally converges into the neighborhood of the global minima.Under the same parameters, defining the search space is expected to improve the identification accuracy of PSO.
Take the identification of X(1) and X(2) as an example.As shown in Figure 3, when the ranges of X(1) and X(2) are unknown, the search space of PSO is 2-dimensional (n is the number of parameters to-be-identified) infinite domain (the blue rectangle).If the number of iterations is insufficient or the settings are inappropriate, the PSO will converge to the local minima.In order to solve this problem, the X's search space (the yellow dotted line rectangle) defined by the RLS is adopted to replace the infinite domain.The search space of PSO is greatly reduced.Under the same number of iterations, PSO is more likely to converge into a small neighborhood of real values (the red dotted circle).Since a much smaller search space is defined, a large number of local minima are excluded from the new search space.
Appl.Sci.2018, 8, x FOR PEER REVIEW 6 of 18 X was identified in each iteration of RLS to define the fluctuation range of X.Therefore, Equation ( 5) was converted into a recursive form as Equation (6).Equation ( 6) is the RLS.
where kg represents the number of iterations of RLS.The detailed derivation processes of Equation ( 6) and the principle of RLS can be referred to [41].The fluctuation ranges of X given by the RLS is the search space of PSO.

PSO with a Finite Search Space
PSO is an algorithm based on particle swarm intelligence.PSO is mainly used for parameter identification and optimization.In each iteration, PSO adjusts the direction and distance of the flight based on the best solution of the particle's own motion trajectories and the best solution of the whole particles' motion trajectories.Through continuous iterations, PSO finally converges into the neighborhood of the global minima.Under the same parameters, defining the search space is expected to improve the identification accuracy of PSO.
Take the identification of ( ) 1 X and ( ) 2 X as an example.As shown in Figure 3, when the ranges of ( ) are unknown, the search space of PSO is 2-dimensional (n is the number of parameters to-be-identified) infinite domain (the blue rectangle).If the number of iterations is insufficient or the settings are inappropriate, the PSO will converge to the local minima.In order to solve this problem, the X's search space (the yellow dotted line rectangle) defined by the RLS is adopted to replace the infinite domain.The search space of PSO is greatly reduced.Under the same number of iterations, PSO is more likely to converge into a small neighborhood of real values (the red dotted circle).Since a much smaller search space is defined, a large number of local minima are excluded from the new search space.

Identified value
Searching space Local best Linearly decreasing inertia weight was adopted in PSO.At the beginning of the iteration, the global search performance of PSO was good.PSO converged to the vicinity of the global minima Linearly decreasing inertia weight was adopted in PSO.At the beginning of the iteration, the global search performance of PSO was good.PSO converged to the vicinity of the global minima quickly.In the latter part of the iteration, the local search performance of PSO was good, so as to accurately approximate the global minima.
The PSO initialized the position S 1 i and velocity V 1 i of the particles by Equation (7).
where max(S 1 i ) and min(S 1 i ) represent the maximum and minimum value of X(i)'s search space, respectively.max(V 1 i ) and min(V 1 i ) represent the maximum and minimum value of X(i)'s velocity, respectively.To ensure that the particles were always iteratively searched within the search space, once a particle's position or velocity exceeded the given range, the value was reset to its nearest maximum or minimum.The PSO performed G iterations.In the kg(kg > 1, kg ∈ N+) iteration, the PSO updated the positions and velocities of the particles by Equation (8).
where S kg i and V kg i represent the position and velocity of particle X(i) in the kg iteration, respectively.r 1 and r 2 represent the normal distribution random numbers between [0, 1].c 1 and c 2 represent the learning factor of the particle.w(kg) represents the inertia time-variant weight in the kg iteration.p kg−1 gb and p kg−1 i_lb represent the global minima of the kg − 1 iteration and the local minima of the particle X(i), respectively.Equations ( 7) and ( 8) are the PSO with a finite search space.Detailed derivation and principle of PSO can be found in [42].

Flowchart of RLS-PSO
Based on the established RLS and PSO, RLS-PSO was designed.The hip and knee joints of APAL swing according to the hip and knee trajectories q 1 and q 2 , respectively.The sensors measure hip and knee torques represented by T 1_m and T 2_m , respectively.Hip and knee angles respectively represented by q 1_m and q 2_m are also measured by sensors.The hip and knee joints angular velocities respectively represented by .q 1_m and .q 2_m are calculated by forward difference.The hip and knee joints angular accelerations respectively represented by .. q 1_m and .. q 2_m are calculated by quadratic forward difference.Figure 4 is the flowchart of RLS-PSO.
The steps of RLS-PSO are as follows: 1. RLS identifies the fluctuation range of each parameter in X by Equation ( 6), thereby defining the search space of PSO; 2.
Within the search space defined by the RLS, the PSO optimizes X by Equations ( 7) and (8).
The estimated values of the hip and knee torques respectively represented by T 1_est and T 2_est are calculated by substituting X identified in each iteration into Equation (4).T 1_est and T 2_est are subtracted from T 1_m and T 2_m , respectively.The absolute values of the differences are ∆T 1 and ∆T 2 , respectively.In each iteration, the optimization goal of the PSO is: When the iteration reaches G times or ∆T 1 + ∆T 2 < ε (ε = 10 −2 , G = 10000), the PSO stops searching and the global minima X Gbest is the identified parameter vector.The steps of RLS-PSO are as follows: 1. RLS identifies the fluctuation range of each parameter in X by Equation ( 6), thereby defining the search space of PSO; 2. Within the search space defined by the RLS, the PSO optimizes X by Equations ( 7) and (8) T Δ , respectively.In each iteration, the optimization goal of the PSO is: 3. When the iteration reaches G times or Δ Δ < ), the PSO stops searching and the global minima Gbest X is the identified parameter vector.

Data Acquisition and Discussion
In order to verify the validity of the RLS-PSO, the data acquisition of h H , k H , 1 T , and 2 T under given trajectories was completed on the APAL.Based on the identified torque, the identification accuracy of LS, PSO, and RLS-PSO has been compared quantitatively.All the parameter identification algorithms were compiled in Python.

Data Acquisition
The acceleration signal is affected by the vibration of the pedestal [42], resulting in an increase in the identification error.The vibration caused by the simultaneous swing of the hip and knee joint is larger than the vibration of the single knee joint.Therefore, in order to reduce the influence of the vibration of the pedestal, a segmental identification strategy was adopted.First, the hip joint was fixed.The knee joint swung.At this time, the vibration of the pedestal was relatively small.The identification of (6:10)  X was completed first.By then, the hip and knee joints swung together to complete the identification of (1:5) X .The flowchart of segment identification is shown in Figure 5.
There are three steps of the segmental identification strategy:

Data Acquisition and Discussion
In order to verify the validity of the RLS-PSO, the data acquisition of H h , H k , T 1 , and T 2 under given trajectories was completed on the APAL.Based on the identified torque, the identification accuracy of LS, PSO, and RLS-PSO has been compared quantitatively.All the parameter identification algorithms were compiled in Python.

Data Acquisition
The acceleration signal is affected by the vibration of the pedestal [42], resulting in an increase in the identification error.The vibration caused by the simultaneous swing of the hip and knee joint is larger than the vibration of the single knee joint.Therefore, in order to reduce the influence of the vibration of the pedestal, a segmental identification strategy was adopted.First, the hip joint was fixed.The knee joint swung.At this time, the vibration of the pedestal was relatively small.The identification of X(6 : 10) was completed first.By then, the hip and knee joints swung together to complete the identification of X(1 : 5).The flowchart of segment identification is shown in Figure 5.There are three steps of the segmental identification strategy: 1.
The hip joint is fixed.The knee joint follows the target trajectory.T 2 and H k are acquired.Substitute T 2 and H k (6 : 10) into Equation ( 6) to define the range of X(6 : 10) as the search space of the PSO.The identification of X(6 : 10) is completed by Equation ( 7) and ( 8); 2.
Both the hip and knee joints follow the target trajectory.After that, H h is acquired.Substitute H h (6 : 10) and X(6 : 10) into equation ( 9) to calculate T 2 1 .T 2 1 is the torque generated by the hip joint to resist the interference of the swing on the knee joint.T 1 is calculated by equation ( 4).T 1 is the difference between T 1 and T 2 1 ; T 2 1 = H h (6 : 10)X(6 : 10) (9)

3.
Substituting H h and T 1 into Equation ( 6) defines the range of X(1 : 5) as the search space of the PSO.X(1 : 5) is identified by Equations ( 7) and ( 8).

RLS-PSO Multiply Knee joint swing
Hip & knee joint swing  As shown in Figure 6, the APAL's trunk and left lower extremity were fixed to the pedestal and the right lower extremity was suspended.The PID control law was adopted to make the hip and knee joint track the given trajectories in Table 1.The sampling frequency was 1000 Hz.When people walk and run, the frequency of motion is 0.85~5 Hz [43,44].According to the sampling theorem, the measured and calculated data were zero-phase filtered with a 10 Hz low-pass filter.Compared to differential filters, zero-phase filters do not cause signal phase shifts and waveform distortion.Therefore, a zero-phase filter was adopted to reduce the effect of waveform distortion on the accuracy of the identification.
The tracking accuracy and response time of the hip and knee joints in the experiment are shown in Table 2.As shown in Figure 6, the APAL's trunk and left lower extremity were fixed to the pedestal and the right lower extremity was suspended.The PID control law was adopted to make the hip and knee joint track the given trajectories in Table 1.The sampling frequency was 1000 Hz.When people walk and run, the frequency of motion is 0.85~5 Hz [43,44].According to the sampling theorem, the measured and calculated data were zero-phase filtered with a 10 Hz low-pass filter.Compared to differential filters, zero-phase filters do not cause signal phase shifts and waveform distortion.Therefore, a zero-phase filter was adopted to reduce the effect of waveform distortion on the accuracy of the identification.Table 2 shows that APAL has good accuracy and response performance.The hip and knee joints (a) (b)

The Identification of Parameters
In order to verify the accuracy of the RLS-PSO, a quantitative accuracy evaluation method was adopted.In order to verify the effectiveness of RLS-PSO in improving the accuracy of exoskeleton parameter identification, the identification accuracy of LS, PSO, and RLS-PSO was quantitatively The tracking accuracy and response time of the hip and knee joints in the experiment are shown in Table 2. Table 2 shows that APAL has good accuracy and response performance.The hip and knee joints accurately tracked the given trajectories in Table 1.The measured data can be adopted for parameter identification.

The Identification of Parameters
In order to verify the accuracy of the RLS-PSO, a quantitative accuracy evaluation method was adopted.In order to verify the effectiveness of RLS-PSO in improving the accuracy of exoskeleton parameter identification, the identification accuracy of LS, PSO, and RLS-PSO was quantitatively compared.To this end, the method in Figure 7 was designed which can quantitatively evaluate the accuracy of the identification algorithm.

Quantitative Evaluation Method for Identification Accuracy
Since the real values of X are unknown, the accuracy of the identification algorithm cannot be evaluated by comparing the identification values with the real values.The original intention and application of exoskeleton parameter identification is to calculate joint torques.Therefore, using the error between the identified torque and the measured torque to evaluate the identification accuracy is consistent with the original intention and application of exoskeleton parameter identification.The method of evaluating identification accuracy by the error between identified torque and measured torque has been verified and applied [21,32].Ghan et al. evaluated the identification accuracy by the agreement between the identified torque and measured torque [21].Bing et al. used the comparison of identified torque and measured torque for evaluating the identification accuracy and comparing between PSO and LS.In summary, when the parameters are unknown, it is feasible and effective to evaluate identification accuracy by the errors of identified torque and measured torque in the exoskeleton parameter identification whose original intention and application is the calculation of joint torques.Therefore, this method was also adopted in this paper to evaluate identification accuracy as shown in Figure 7.The method can quantitatively evaluate the accuracy of the identification algorithm.Take the experiment in Figure 6a as an example.The identification of X was implemented based on the data measured from trajectory 1. H 2 was acquired in trajectory 2. T 2 was calculated from Equation ( 4).The difference between T 2 and T 2 was the torque identification error T error .The identification of (6:10)  X was implemented by the data acquired from trajectory 1.The search space of (6:10) X in Table 3 was defined by RLS.(6:10) X was identified by LS, PSO, and

Quantitative Analysis of RLS-PSO Identification Accuracy
The identification of X(6 : 10) was implemented by the data acquired from trajectory 1.The search space of X(6 : 10) in Table 3 was defined by RLS.X(6 : 10) was identified by LS, PSO, and RLS-PSO, respectively.The identification values are shown in Table 4.In each iteration of PSO and RLS-PSO, the best particle which obtains the minimum error is selected.The best particle convergence of PSO and RLS-PSO is shown in Figure 8.The LS does not have iteration processes, so LS is not shown in Figure 8.  Trajectory 2 was the reference group, and the identification values of (6:10) X was quantified by the quantitative evaluation method in Figure 7. Trajectory 2 was the reference group, and the identification values of X(6 : 10) was quantified by the quantitative evaluation method in Figure 7. Figure 9 is a comparison of the identified torque, the measured torque, and the torque error T error .The absolute values of T error were averaged.The mean ± std of the absolute torque error of LS as shown in Figure 10 is 1.0427 ± 0.9047 Nm.The mean ± std of the absolute torque error of the PSO is 0.5126 ± 0.4386 Nm.The mean ± std of the absolute torque error of the RLS-PSO is 0.1887 ± 0.2475 Nm.As shown in Figure 10, it is clear that the most accurate X(1 : 5) was identified by the RLS-PSO.Under the same conditions, the accuracy of RLS-PSO is improved by 81.9% compared with LS, and the accuracy of RLS-PSO is improved by 63.19% compared with PSO.Similarly, Table 5 is the search space of X(1 : 5) defined by RLS.The segmental identification strategy was adopted to identify X(1 : 5).Table 6 is the identified values of X(1 : 5).The best particle convergence of PSO and RLS-PSO is shown in Figure 11.     Figure 12 shows the identified torque, measured torque, and torque error T error .The mean ± std of the absolute torque error of LS as shown in Figure 13 11, the number iterations are 120 and 280 for RLS-PSO and PSO, respectively.Obviously, the convergence rate of RLS-PSO is larger than that of PSO.
It can be seen from Figure 10 and Figure 13 that the accuracy of RLS-PSO is significantly higher than LS and PSO.This result shows that under the same conditions, pre-defining the search space of PSO through RLS can effectively improve the identification accuracy of PSO.
Compared with the LS adopted by BLEEX, RLS-PSO does not require static experiments, which not only reduces the complexity of the experiments but also ensures the accuracy of dynamic parameter identification.The hydraulic drive units caused fluctuations in the dynamic parameters during the movement of the exoskeleton.In order to solve this problem, the fluctuation range of each parameter to-be-identified was defined by RLS and used as the search space of PSO.Accurate parameter identification was achieved by PSO in the search space.RLS-PSO converted the problem from the time-variant of the parameters to the optimization of the parameters.It not only solved the time-variant problem of parameters, but also reduced the potential local minima of PSO in the infinite domain.Under the same experimental condition, compared with LS and PSO, RLS-PSO greatly improved the identification accuracy.Comparing Figure 10 and Figure 13, it can be seen that the identification accuracy of (6:10) X is higher than that of (1:5) X .The main reason for this phenomenon is that the large inertia of the thigh and shank has a larger effect on the pedestal than the shank only.As the pedestal's vibration increased, the interference of the angular acceleration signal increased.In order to minimize the errors caused by the vibration, the segmental identification method was adopted.However, the identification error of (6:10) X still accumulated in the identification of (1:5) X .Despite this, the identification accuracy of RLS-PSO is still higher than the accuracy of LS and PSO.Furthermore, Figure 8 and Figure 11 show that the defined RLS search space led to a better convergence performance of RLS-PSO than that of PSO.Thus, the effectiveness of reducing the search space of the PSO from infinite domain to finite domain is verified.In Figure 11, the number iterations are 120 and 280 for RLS-PSO and PSO, respectively.Obviously, the convergence rate of RLS-PSO is larger than that of PSO.

Conclusions
It can be seen from Figures 10 and 13 that the accuracy of RLS-PSO is significantly higher than LS and PSO.This result shows that under the same conditions, pre-defining the search space of PSO through RLS can effectively improve the identification accuracy of PSO.
Compared with the LS adopted by BLEEX, RLS-PSO does not require static experiments, which not only reduces the complexity of the experiments but also ensures the accuracy of dynamic parameter identification.The hydraulic drive units caused fluctuations in the dynamic parameters during the movement of the exoskeleton.In order to solve this problem, the fluctuation range of each parameter to-be-identified was defined by RLS and used as the search space of PSO.Accurate parameter identification was achieved by PSO in the search space.RLS-PSO converted the problem from the time-variant of the parameters to the optimization of the parameters.It not only solved the time-variant problem of parameters, but also reduced the potential local minima of PSO in the infinite domain.Under the same experimental condition, compared with LS and PSO, RLS-PSO greatly improved the identification accuracy.
Comparing Figures 10 and 13, it can be seen that the identification accuracy of X(6 : 10) is higher than that of X(1 : 5).The main reason for this phenomenon is that the large inertia of the thigh and shank has a larger effect on the pedestal than the shank only.As the pedestal's vibration increased, the interference of the angular acceleration signal increased.In order to minimize the errors caused by the vibration, the segmental identification method was adopted.However, the identification error of X(6 : 10) still accumulated in the identification of X(1 : 5).Despite this, the identification accuracy of RLS-PSO is still higher than the accuracy of LS and PSO.Furthermore, Figures 8 and 11 show that the defined RLS search space led to a better convergence performance of RLS-PSO than that of PSO.Thus, the effectiveness of reducing the search space of the PSO from infinite domain to finite domain is verified.

f and 1 vf 2 c f and 2 vf
represent the friction parameters of the hip joint, represent the friction parameter of the knee joint.

Figure 5 .
Figure 5. Flowchart of segmental identification strategy.The trajectories of the hip and knee joints of APAL are shown in Table 1.No.1 and No.2 are the trajectories of knee joint swing experiments.No.3 and No.4 are the trajectories of the hip and knee joints swing experiments.

Figure 6 .
Figure 6.Data acquisition of APAL: (a) Hip joint was fixed and knee joint swung; (b) Hip and knee joints swung.

Figure 6 .
Figure 6.Data acquisition of APAL: (a) Hip joint was fixed and knee joint swung; (b) Hip and knee joints swung.

Figure 7 .
Figure 7. Quantitative evaluation method for identification accuracy.

Figure 8 .
Figure 8.The convergences of X (6:10) in PSO and RLS-PSO: (a) The convergence of best particle in RLS-PSO's iterations; (b) The convergence of best particle in PSO's iterations.

Figure 9 Figure 8 .
Figure 8.The convergences of X (6:10) in PSO and RLS-PSO: (a) The convergence of best particle in RLS-PSO's iterations; (b) The convergence of best particle in PSO's iterations.

Figure 10 .
Figure 10.The mean ± std. of the absolute errors.Figure 10.The mean ± std. of the absolute errors.

Figure 10 .
Figure 10.The mean ± std. of the absolute errors.Figure 10.The mean ± std. of the absolute errors.

Figure 10 .
Figure 10.The mean ± std. of the absolute errors.

Figure 11 .
Figure 11.The convergences of X(1:5) in PSO and RLS-PSO: (a) The convergence of best particle in RLS-PSO's iterations; (b) The convergence of best particle in PSO's iterations.

Figure 12
Figure 12 shows the identified torque, measured torque, and torque error error T .The mean ± std of the absolute torque error of LS as shown in Figure 13 is 4.9922 ± 3.7366 Nm.The mean ± std of absolute torque error of PSO is 4.3516 ± 3.3602 Nm.The mean ± std of absolute torque error of RLS-PSO is 1.2128 ± 3.100 Nm.The RLS-PSO obsessed the best accuracy.Compared with the LS, the accuracy of the RLS-PSO is improved by 75.70%.Compared with the PSO, the accuracy of the RLS-PSO is improved by 72.13%.

Figure 11 .
Figure 11.The convergences of X(1:5) in PSO and RLS-PSO: (a) The convergence of best particle in RLS-PSO's iterations; (b) The convergence of best particle in PSO's iterations.

Figure 11 .
Figure12shows the identified torque, measured torque, and torque error T error .The mean ± std of the absolute torque error of LS as shown in Figure13is 4.9922 ± 3.7366 Nm.The mean ± std

Figure 12 Figure 12 .
Figure 12 shows the identified torque, measured torque, and torque error error T .The mean ± std of the absolute torque error of LS as shown in Figure 13 is 4.9922 ± 3.7366 Nm.The mean ± std of absolute torque error of PSO is 4.3516 ± 3.3602 Nm.The mean ± std of absolute torque error of RLS-PSO is 1.2128 ± 3.100 Nm.The RLS-PSO obsessed the best accuracy.Compared with the LS, the accuracy of the RLS-PSO is improved by 75.70%.Compared with the PSO, the accuracy of the RLS-PSO is improved by 72.13%.

Figure 13 .
Figure 13.The mean ± std of the absolute errors.

Figure 13 .
Figure 13.The mean ± std of the absolute errors.

Table 1 .
Trajectories of hip and knee joints.

Table 2 .
Mean of tracking accuracy and response time of APAL.

Table 1 .
Trajectories of hip and knee joints.

Table 2 .
Mean of tracking accuracy and response time of APAL.

Table 3 .
The search space of X(6 : 10) defined by RLS.

Table 5 .
The search space of X(1 : 5) defined by RLS.