Next Article in Journal
Nonlinear Refractive Index in Rectangular Graphene Quantum Dots
Next Article in Special Issue
Spreadsheets as Interoperability Solution for Business Process Representation
Previous Article in Journal
Combining Subgoal Graphs with Reinforcement Learning to Build a Rational Pathfinder
Previous Article in Special Issue
Background Knowledge Based Multi-Stream Neural Network for Text Classification
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
State Key Laboratory of Robotics and System, Harbin Institute of Technology (HIT), Harbin 150001, China
2
Shenzhen Academy of Aerospace Technology, Shenzhen 518057, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(2), 324; https://doi.org/10.3390/app9020324
Submission received: 18 December 2018 / Revised: 9 January 2019 / Accepted: 11 January 2019 / Published: 17 January 2019
(This article belongs to the Special Issue Applied Sciences Based on and Related to Computer and Control)

Abstract

:
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.

1. 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.

2. 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.

2.1. 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.
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.
T h = [ I s + m s ( ( L s L G s ) 2 h G s 2 ) + I t + m t ( ( L t L G t ) 2 h G t 2 ) + m s L t 2 + 2 L t ( m s ( L s L G s ) cos q 2 m s h G s sin q 2 ) ] q ¨ 1 + [ I s + m s ( ( L s L G s ) 2 h G s 2 ) + L t ( m s ( L s L G s ) cos q 2 m s h G s sin q 2 ) ] q ¨ 2 + L t ( m s ( L s L G s ) sin q 2 + m s h G s cos q 2 ) ( q ˙ 1 2 ( q ˙ 1 + q ˙ 2 ) 2 ) + g [ ( m t ( L t L G t ) + m s L t ) sin q 1 + m t h G t cos q 1 + m s ( L s L G s ) sin ( q 1 + q 2 ) + m s h G s cos ( q 1 + q 2 ) ] T k = [ I s + m s ( ( L s L G s ) 2 h G s 2 ) + L t ( m s ( L s L G s ) cos q 2 m s h G s sin q 2 ) ] q ¨ 1 + [ I s + m s ( ( L s L G s ) 2 h G s 2 ) ] q ¨ 2 + L t ( m s ( L s L G s ) sin q 2 + m s h G s cos q 2 ) q ˙ 1 2 + g [ m s ( L s L G s ) sin ( q 1 + q 2 ) + m s h G s cos ( q 1 + q 2 ) ]
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.
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 .
T 1 = T h + T f 1 T 2 = T k + T f 2
where
T f 1 = f c 1 sgn ( q ˙ 1 ) + f v 1 q ˙ 1 T f 2 = f c 2 sgn ( q ˙ 2 ) + f v 2 q ˙ 2
f c 1 and f v 1 represent the friction parameters of the hip joint, f c 2 and f v 2 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 G t , L G s , h G t , h G s 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 G t , L G s , h G t , h G s in Equation (2) were parameterized into a set of inertial parameters M x t , M x s , M y t , M y s , 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.
T 1 = [ J t + 2 L t ( M x s cos q 2 M y s sin q 2 ) ] q ¨ 1 + [ J s + L t ( M x s cos q 2 M y s sin q 2 ) ] q ¨ 2 + L t ( M x s sin q 2 + M y s cos q 2 ) [ q ˙ 1 2 ( q ˙ 1 + q ˙ 2 ) 2 ] + g [ M x t sin q 1 + M y t cos q 1 + M x s sin ( q 1 + q 2 ) + M y s cos ( q 1 + q 2 ) ] + f c 1 sgn ( q ˙ 1 ) + f v 1 q ˙ 1 T 2 = [ J s + L t ( M x s cos q 2 M y s sin q 2 ) ] q ¨ 1 + J s q ¨ 2 + L t ( M x s sin q 2 + M y s cos q 2 ) q ˙ 1 2 + g [ M x s sin ( q 1 + q 2 ) + M y s cos ( q 1 + q 2 ) ] + f c 2 sgn ( q ˙ 2 ) + f v 2 q ˙ 2
where
M x t = m t ( L t L G t ) + m s L t M y t = m t h G t M x s = m s ( L s L G s ) M y s = m s h G s J t = I s + m s [ ( L s L G s ) 2 h G s 2 ] + I t + m t [ ( L t L G t ) 2 h G t 2 ] + m s L t 2 J s = I s + m s [ ( L s L G s ) 2 h G s 2 ]
Equation (3) is the linearization dynamic equation of APAL, where M x t , M x s , M y t , M y s , J t , J s , f c 1 , f c 2 , f v 1 , f v 2 are the parameters to-be-identified.

2.2. 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.

2.2.1. 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 x t , M x s , M y t , M y s , J t , J s , f c 1 , f c 2 , f v 1 , f v 2 . The fluctuation ranges of M x t , M x s , M y t , M y s , J t , J s , f c 1 , f c 2 , f v 1 , f v 2 form a 10-dimension search space. In order to facilitate computer programming and calculation, Equation (3) was converted into matrix form as Equation (4).
T 1 = H h ( q 1 , q ˙ 1 , q ¨ 1 , q 2 , q ˙ 2 , q ¨ 2 ) X T 2 = H k ( q 1 , q ˙ 1 , q ¨ 1 , q 2 , q ˙ 2 , q ¨ 2 ) X
where
X = [ M x t , M y t , J t , f c 1 , f v 1 , M x s , M y s , J s , f c 2 , f v 2 ] T H h = [ g sin q 1 , g cos q 1 , q ¨ 1 , s g n ( q ˙ 1 ) , q ˙ 1 , 2 L t q ¨ 1 cos q 2 + L t q ¨ 2 cos q 2 + L t sin q 2 [ q ˙ 1 2 ( q ˙ 1 + q ˙ 2 ) 2 ] + g sin ( q 1 + q 2 ) , 2 L t q ¨ 1 sin q 2 L t q ¨ 2 sin q 2 + L t [ q ˙ 1 2 ( q ˙ 1 + q ˙ 2 ) 2 ] cos q 2 + g cos ( q 1 + q 2 ) , q ¨ 2 , s g n ( q ˙ 1 ) , q ˙ 1 ] H k = [ 0 , 0 , 0 , 0 , 0 , L t q ¨ 1 cos q 2 + L t q ˙ 1 2 sin q 2 + g sin ( q 1 + q 2 ) , L t q ¨ 1 sin q 2 + L t q ˙ 1 2 cos q 2 + g cos ( q 1 + q 2 ) , q ¨ 1 + q ¨ 2 , s g n ( q ˙ 2 ) , q ˙ 2 ]
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 = ( H i T H i ) 1 H i T T j ( i = h , k ; j = 1 , 2 )
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.
{ K = P ( k g 1 ) H i / ( H i T P ( k g 1 ) H i + 1 / w e i g h t ) n e w = T j ( k g ) H i T X ( k g ) X ( k g + 1 ) = X ( k g ) + K n e w P ( k g ) = 1 / w e i g h t ( e y e ( 2 ) K H i T ) P ( k g 1 ) e y e ( 2 ) = ( 1 0 0 1 ) ( k g > 1 , k g N + ; i = h , k ; j = 1 , 2 )
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.

2.2.2. 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.
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 i 1 and velocity V i 1 of the particles by Equation (7).
S i 1 = min ( S i 1 ) + r a n d × [ max ( S i 1 ) min ( S i 1 ) ] V i 1 = min ( V i 1 ) + r a n d × [ max ( V i 1 ) min ( V i 1 ) ] r a n d [ 0 , 1 ] , ( i = 1 , 2 , , 10 )
where max ( S i 1 ) and min ( S i 1 ) represent the maximum and minimum value of X ( i ) ’s search space, respectively. max ( V i 1 ) and min ( V i 1 ) 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 ( k g > 1 , k g N + ) iteration, the PSO updated the positions and velocities of the particles by Equation (8).
V i k g = w ( t ) × V i k g 1 + c 1 r 1 ( p g b k g 1 S i k g 1 ) + c 2 r 2 ( p i _ l b k g 1 S i k g 1 ) S i k g = S i k g 1 + V i k g , ( r 1 = r a n d , r 2 = r a n d ) w ( k g ) = w m a x k g G ( w m a x w m i n ) r a n d [ 0 , 1 ] , ( i = 1 , 2 , , 10 )
where S i k g and V i k g represent the position and velocity of particle X ( i ) in the k g 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 ( k g ) represents the inertia time-variant weight in the k g iteration. p g b k g 1 and p i _ l b k g 1 represent the global minima of the k g 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].

2.2.3. 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:
  • RLS identifies the fluctuation range of each parameter in X by Equation (6), thereby defining the search space of PSO;
  • 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 _ e s t and T 2 _ e s t are calculated by substituting X identified in each iteration into Equation (4). T 1 _ e s t and T 2 _ e s t 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: Δ T 1 + Δ T 2 < ε , ε 0 + ;
  • When the iteration reaches G times or Δ T 1 + Δ T 2 < ε ( ε = 10 2 , G = 10000 ), the PSO stops searching and the global minima X G b e s t is the identified parameter vector.

3. 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.

3.1. 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:
  • 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);
  • 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 1 2 . T 1 2 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 1 2 ;
    T 1 2 = H h ( 6 : 10 ) X ( 6 : 10 )
  • 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).
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.
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.
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.

3.2. 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.

3.2.1. 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 e r r o r .

3.2.2. 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 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 e r r o r .
The absolute values of T e r r o r 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 e r r o r . 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%.
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.
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 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, 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.

4. Conclusions

In summary, this paper developed RLS-PSO to improve the identification accuracy of exoskeleton dynamic parameters. The effectiveness of RLS-PSO in improving the identification accuracy was verified on the APAL. The RLS-PSO has three main parts: First, the dynamic equation of the APAL in the sagittal plane was derived. The dynamic equation was parameterized. A linearized dynamic equation and a vector consisting of the parameters to-be-identified were derived. Second, the search space of each parameter in the vector to-be-identified was defined by RLS to eliminate as many potential local minima as possible in the iterations of the PSO. RLS converted the problem from the time-variant of the parameters to the optimization of the parameters. Third, the iterative optimization was performed in the defined search space by PSO. According to the quantitative identification accuracy evaluation index, the accuracy of RLS-PSO was higher than that of LS and PSO.
The above results show that although the hydraulic drive units cause the dynamic parameters of the exoskeleton to be time-variant, the RLS-PSO can utilize the fluctuation of the parameters to define a smaller PSO search space. Thereby, more accurate parameter identification can be achieved. The RLS-PSO proposed in this paper is expected to improve the accuracy of the dynamic model in model-based exoskeleton control.

Author Contributions

W.S. conceived the method and wrote the paper; F.Z. and W.G. helped to modify it; J.D. and S.Q. performed the experiments and analyzed the data; F.Z. and X.W. contributed theoretical justification and reference materials. All authors have read and approved the final manuscript.

Funding

This research was funded by National Natural Science Foundation of China (No.61473015, JCYJ20160425150757025, No. U1713222 and No.61773139), and Shenzhen Peacock Plan (KQTD2016112515134654).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hiroaki, K.; Yoshiyuki, S. Power assist method based on Phase Sequence and muscle force condition for HAL. Adv. Robot. 2005, 19, 717–734. [Google Scholar]
  2. Aguirre-Ollinger, G.; Nagarajan, U.; Goswami, A. An admittance shaping controller for exoskeleton assistance of the lower extremities. Auton. Robot. 2016, 40, 1–28. [Google Scholar] [CrossRef]
  3. Collins, S.H.; Wiggin, M.B.; Sawicki, G.S. Reducing the energy cost of human walking using an unpowered exoskeleton. Nature 2015, 522, 212–215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Lee, J.W.; Kim, H.; Jang, J.; Park, S. Virtual model control of exoskeleton for load carriage inspired by human behavior. Auton. Robot. 2015, 38, 211–223. [Google Scholar] [CrossRef]
  5. Gregorczyk, K.N.; Hasselquist, L.; Schiffman, J.M.; Bensel, C.K.; Obusek, J.P.; Gutekunst, D.J. Effects of a lower-body exoskeleton device on metabolic cost and gait biomechanics during load carriage. Ergonomics 2010, 53, 1263–1275. [Google Scholar] [CrossRef]
  6. Kazerooni, H. The Berkeley Exoskeleton Project. Experimental Robotics IX; Springe: Berlin/Heidelberg, Germany, 2006; pp. 9–15. [Google Scholar]
  7. Mir-Nasiri, N. Efficient Exoskeleton for Human Motion Assistance. Wearable Robotics: Challenges and Trends; Springer International Publishing: New York, NY, USA, 2017. [Google Scholar]
  8. Cao, H.; Zhu, J.; Xia, C.; Zhou, H.; Chen, X.; Wang, Y. Design and Control of a Hydraulic-Actuated Leg Exoskeleton for Load-Carrying Augmentation. Intelligent Robotics and Applications; Springer: Berlin/Heidelberg, Germany, 2010; pp. 590–599. [Google Scholar]
  9. Zhang, X.; Guo, Q.; Zhao, C.; Zhang, Y.; Luo, X. Development of a lower extremity exoskeleton suit actuated by hydraulic. In Proceedings of the 2012 IEEE International Conference on Mechatronics and Automation, Chengdu, China, 5–8 August 2012; pp. 587–591. [Google Scholar]
  10. Kim, H.; Shin, Y.J.; Kim, J. Design and locomotion control of a hydraulic exoskeleton for mobility augmentation. Mechatronics 2017, 46, 32–45. [Google Scholar] [CrossRef]
  11. Bicchi, A.; Tonietti, G. Fast and soft-arm tactics. IEEE Robot. Autom. Mag. 2004, 11, 22–33. [Google Scholar] [CrossRef]
  12. Tucker, M.R.; Olivier, J.; Pagel, A.; Bleuler, H.; Bouri, M.; Lambercy, O.; del R Millán, J.; Riener, R.; Vallery, H.; Gassert, R. Control strategies for active lower extremity prosthetics and orthotics: A review. J. NeuroEng. Rehabil. 2015, 12, 1. [Google Scholar] [CrossRef] [PubMed]
  13. Durandau, G.; Sartori, M.; Bortole, M.; Moreno, J.C.; Pons, J.L.; Farina, D. Real-Time Modeling for Lower Limb Exoskeletons. In Wearable Robotics: Challenges and Trends; Biosystems & Biorobotics; González-Vargas, J., Ibáñez, J., Contreras-Vidal, J., van der Kooij, H., Pons, J., Eds.; Springer: Cham, Germany, 2017; Volume 16. [Google Scholar]
  14. Moren, J.C.; Brunetti, F.; Navarro, E.; Forner-Cordero, A.; Pons, J.L. Analysis of the human interaction with a wearable lower-limb exoskeleton. Appl. Bionics Biomech. 2009, 2, 245–256. [Google Scholar] [CrossRef]
  15. Manns, P.; Sreenivasa, M.; Millard, M.; Mombaur, K. Motion Optimization and Parameter Identification for a Human and Lower Back Exoskeleton Model. IEEE Robot. Autom. Lett. 2017, 2, 1564–1570. [Google Scholar] [CrossRef] [Green Version]
  16. Cong, L.; Wu, D.; Long, Y.; Du, Z.; Dong, W. Parameter Identification Based Sensitivity Amplification Control for Lower Extremity Exoskeleton. In Proceedings of the 2017 International Conference on Artificial Intelligence, Automation and Control Technologies, Wuhan, China, 7–9 April 2017; pp. 1–6. [Google Scholar] [CrossRef]
  17. Bertolini, A. Wearable Robots: A Legal Analysis. Wearable Robotics: Challenges and Trends; Springer International Publishing: New York, NY, USA, 2017. [Google Scholar]
  18. Dollar, A.M.; Herr, H. Exoskeletons and Active Orthoses: Challenges and State-of-the-Art. IEEE Trans. Robot. 2008, 24, 144–158. [Google Scholar] [CrossRef]
  19. Long, Y.; Du, Z.; Cong, L.; Wang, W.; Zhang, Z.; Dong, W. Active disturbance rejection control based human gait tracking for lower extremity rehabilitation exoskelton. ISA Trans. 2017, 67, 389. [Google Scholar] [CrossRef] [PubMed]
  20. Vantilt, J.; Aertbeliën, E.; De Groote, F.; De Schutter, J. Optimal excitation and identification of the dynamic model of robotic systems with compliant actuators. In Proceedings of the IEEE International Conference on Robotics & Automation, Seattle, WA, USA, 26–30 May 2015; pp. 2117–2124. [Google Scholar]
  21. Ghan, J.; Kazerooni, H. System identification for the Berkeley exoskeleton (BLEEX). In Proceedings of the IEEE International Conference on Robotics and Automation, Orlando, FL, USA, 15–19 May 2006; pp. 3477–3484. [Google Scholar]
  22. Steger, R.; Kim, S.H.; Kazerooni, H. Control scheme and networked control architecture for the Berkeley exoskeleton (BLEEX). In Proceedings of the IEEE International Conference on Robotics and Automation, Orlando, FL, USA, 15–19 May 2006; pp. 3469–3476. [Google Scholar]
  23. Gautier, M.; Briot, S. Dynamic Parameter Identification of a 6 DOF Industrial Robot using Power Model. In Proceedings of the IEEE International Conference on Robotics & Automation, Karlsruhe, Germany, 6–10 May 2013. [Google Scholar]
  24. Ogawa, Y.; Venture, G.; Ott, C. Dynamic parameters identification of a humanoid robot using joint torque sensors and/or contact forces. In Proceedings of the IEEE-RAS International Conference on Humanoid Robots, Madrid, Spain, 18–20 November 2014. [Google Scholar]
  25. Mrachacz-Kersting, N.; Lavoie, B.A.; Andersen, J.B.; Sinkjær, T. Characterisation of the quadriceps stretch reflex during the transition from swing to stance phase of human walking. Exp. Brain Res. 2004, 159, 108–122. [Google Scholar] [CrossRef] [PubMed]
  26. Pantall, A.; Gregor, R.J.; Prilutsky, B.I. Stance and swing phase detection during level and slope walking in the cat: Effects of slope, injury, subject and kinematic detection method. J. Biomech. 2012, 45, 1529–1533. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Doranga, S.; Wu, C.Q. Parameter Identification for Nonlinear Dynamic Systems via Multilinear Least Square Estimation. In Special Topics in Structural Dynamics; Springer: Cham, Germany, 2014; Volume 6. [Google Scholar]
  28. Teunissen, P.; Montenbruck, O. Teunissen. Least-Squares Estimation and Kalman Filtering. Springer Handbook of Global Navigation Satellite Systems; Springer International Publishing: New York, NY, USA, 2017. [Google Scholar]
  29. Ding, L.; Shan, W.; Zhou, C.; Xi, W. Dynamic Identification for Industrial Robot Manipulators Based on Glowworm Optimization Algorithm. In Proceedings of the International Conference on Intelligent Robotics and Applications, Wuhan, China, 16–18 August 2017; pp. 789–799. [Google Scholar]
  30. Xiu, W.; Zhang, L.; Ma, O. Experimental study of a momentum-based method for identifying the inertia barycentric parameters of a human body. Multibody Syst. Dyn. 2016, 36, 237–255. [Google Scholar] [CrossRef]
  31. Ayusawa, K.; Venture, G.; Nakamura, Y. Identifiability and identification of inertial parameters using the underactuated base-link dynamics for legged multibody systems. Int. J. Robot. Res. 2014, 33, 446–468. [Google Scholar] [CrossRef]
  32. Bingül, Z.; Karahan, O. Dynamic identification of Staubli RX-60 robot using PSO and LS methods. Expert Syst. Appl. 2011, 38, 4136–4149. [Google Scholar] [CrossRef]
  33. Nickabadi, A.; Ebadzadeh, M.M.; Safabakhsh, R. A Novel Particle Swarm Optimization Algorithm with Adaptive Inertia Weight; Elsevier: Amsterdam, The Netherlands, 2011. [Google Scholar]
  34. Zheng, Y.L.; Ma, L.H.; Zhang, L.Y.; Qian, J.X. On the convergence analysis and parameter selection in particle swarm optimization. In Proceedings of the International Conference on Machine Learning & Cybernetics, Xi’an, China, 5 November 2004. [Google Scholar]
  35. Khalil, W.; Dombre, E. Modeling, Identification and Control of Robots; Taylor & Francis, Inc.: Abingdon, UK, 2003. [Google Scholar]
  36. Deng, J.; Wang, P.; Li, M.; Guo, W.; Zha, F.; Wang, X. Structure design of active power-assist exoskeleton APAL robot. Adv. Mech. Eng. 2017, 9. [Google Scholar] [CrossRef]
  37. Li, M.; Deng, J.; Zha, F.; Qiu, S.; Wang, X.; Chen, F. Towards Online Estimation of Human Joint Muscular Torque with a Exoskeleton Robot. Appl. Sci. 2018, 8, 1610. [Google Scholar] [CrossRef]
  38. Ding, F.; Wang, Y.; Ding, J. Recursive Least Squares Parameter Identification Algorithms for Systems with Colored Noise Using the Filtering Technique and the Auxilary Model; Academic Press, Inc.: Cambridge, MA, USA, 2015. [Google Scholar]
  39. Dolanc, G.; Strmčnik, S. Identification of nonlinear systems using a piecewise-linear Hammerstein model. Syst. Control Lett. 2013, 54, 145–158. [Google Scholar] [CrossRef]
  40. Young, P.C. Recursive Least Squares Estimation. Recursive Estimation and Time-Series Analysis; Springer: Berlin/Heidelberg, Germany, 2011; pp. 29–46. [Google Scholar]
  41. Couceiro, M.; Ghamisi, P. Particle Swarm Optimization. Fractional Order Darwinian Particle Swarm Optimization; Springer International Publishing: New York, NY, USA, 2016; pp. 149–150. [Google Scholar]
  42. Khandelwal, S.; Wickström, N. Evaluation of the performance of accelerometer-based gait event detection algorithms in different real-world scenarios using the MAREA gait database. Gait Posture 2017, 51, 84–90. [Google Scholar] [CrossRef] [PubMed]
  43. Wagenaar, R.C.; van Emmerik, R.E. Resonant frequencies of arms and legs identify different walking patterns. J. Biomech. 2000, 33, 853–861. [Google Scholar] [CrossRef]
  44. Farley, C.T.; González, O. Leg stiffness and stride frequency in human running. J. Biomech. 1996, 29, 181–186. [Google Scholar] [CrossRef]
Figure 1. Active power-assist exoskeleton (APAL): (a) APAL; (b) Schematic diagram of APAL lower extremity mechanical structure.
Figure 1. Active power-assist exoskeleton (APAL): (a) APAL; (b) Schematic diagram of APAL lower extremity mechanical structure.
Applsci 09 00324 g001
Figure 2. Simplified model of APAL’s lower extremity.
Figure 2. Simplified model of APAL’s lower extremity.
Applsci 09 00324 g002
Figure 3. Two-dimensional search space.
Figure 3. Two-dimensional search space.
Applsci 09 00324 g003
Figure 4. Flowchart of the recursive least square–particle swarm optimization (RLS-PSO) method.
Figure 4. Flowchart of the recursive least square–particle swarm optimization (RLS-PSO) method.
Applsci 09 00324 g004
Figure 5. Flowchart of segmental identification strategy.
Figure 5. Flowchart of segmental identification strategy.
Applsci 09 00324 g005
Figure 6. Data acquisition of APAL: (a) Hip joint was fixed and knee joint swung; (b) Hip and knee joints swung.
Figure 6. Data acquisition of APAL: (a) Hip joint was fixed and knee joint swung; (b) Hip and knee joints swung.
Applsci 09 00324 g006
Figure 7. Quantitative evaluation method for identification accuracy.
Figure 7. Quantitative evaluation method for identification accuracy.
Applsci 09 00324 g007
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 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.
Applsci 09 00324 g008
Figure 9. Measured torque, identified torque, and T e r r o r .
Figure 9. Measured torque, identified torque, and T e r r o r .
Applsci 09 00324 g009
Figure 10. The mean ± std. of the absolute errors.
Figure 10. The mean ± std. of the absolute errors.
Applsci 09 00324 g010
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. 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.
Applsci 09 00324 g011aApplsci 09 00324 g011b
Figure 12. Measured torque, identified torque, and T e r r o r .
Figure 12. Measured torque, identified torque, and T e r r o r .
Applsci 09 00324 g012
Figure 13. The mean ± std of the absolute errors.
Figure 13. The mean ± std of the absolute errors.
Applsci 09 00324 g013
Table 1. Trajectories of hip and knee joints.
Table 1. Trajectories of hip and knee joints.
No.Hip TrajectoryKnee Trajectory
1−1045 − 30sin(2πt)
2−1060 − 30sin(2πt)
3−30 − 30sin(2πt)75 − 30sin(2πt)
4−45 − 30sin(2πt)60 + 30sin(2πt)
Table 2. Mean of tracking accuracy and response time of APAL.
Table 2. Mean of tracking accuracy and response time of APAL.
Performance IndexHip AngleKnee Angle
Tracking accuracy (°)0.02170.0358
Response time (s)0.68060.9725
Table 3. The search space of X ( 6 : 10 ) defined by RLS.
Table 3. The search space of X ( 6 : 10 ) defined by RLS.
X ( 6 ) X ( 7 ) X ( 8 ) X ( 9 ) X ( 10 )
0.79~1.79−0.98~0.180.39~0.91 0.01~1.13 0.02~3.15
Table 4. The identified values of X ( 6 : 10 ) .
Table 4. The identified values of X ( 6 : 10 ) .
Type X ( 6 ) X ( 7 ) X ( 8 ) X ( 9 ) X ( 10 )
LS1.3859−0.82670.70190.09430.2907
PSO1.2199−0.40160.59060.33230.1830
RLS-PSO1.1822−0.28080.56550.44060.1020
Table 5. The search space of X ( 1 : 5 ) defined by RLS.
Table 5. The search space of X ( 1 : 5 ) defined by RLS.
X ( 1 ) X ( 2 ) X ( 3 ) X ( 4 ) X ( 5 )
5.38~11.642.34~5.202.01~4.250.10~7.880.01~1.70
Table 6. The identified values of X ( 1 : 5 ) .
Table 6. The identified values of X ( 1 : 5 ) .
X ( 1 ) X ( 2 ) X ( 3 ) X ( 4 ) X ( 5 )
LS6.51242.44812.49961.72850.0245
PSO7.99703.49692.99221.59140.0752
RLS-PSO7.44333.17882.84131.72850.0372

Share and Cite

MDPI and ACS Style

Zha, F.; Sheng, W.; Guo, W.; Qiu, S.; Deng, J.; Wang, X. Dynamic Parameter Identification of a Lower Extremity Exoskeleton Using RLS-PSO. Appl. Sci. 2019, 9, 324. https://doi.org/10.3390/app9020324

AMA Style

Zha F, Sheng W, Guo W, Qiu S, Deng J, Wang X. Dynamic Parameter Identification of a Lower Extremity Exoskeleton Using RLS-PSO. Applied Sciences. 2019; 9(2):324. https://doi.org/10.3390/app9020324

Chicago/Turabian Style

Zha, Fusheng, Wentao Sheng, Wei Guo, Shiyin Qiu, Jing Deng, and Xin Wang. 2019. "Dynamic Parameter Identification of a Lower Extremity Exoskeleton Using RLS-PSO" Applied Sciences 9, no. 2: 324. https://doi.org/10.3390/app9020324

APA Style

Zha, F., Sheng, W., Guo, W., Qiu, S., Deng, J., & Wang, X. (2019). Dynamic Parameter Identification of a Lower Extremity Exoskeleton Using RLS-PSO. Applied Sciences, 9(2), 324. https://doi.org/10.3390/app9020324

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

Article Metrics

Back to TopTop