Next Article in Journal
High-Resolution Hyperspectral Imaging Using Low-Cost Components: Application within Environmental Monitoring Scenarios
Previous Article in Journal
Secure IIoT Information Reinforcement Model Based on IIoT Information Platform Using Blockchain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Knee Extension Force Using Mechanomyography Signals Based on GRA and ICS-SVR

1
Institute of Intelligent Machines, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei 230031, China
2
Department of Science Island, University of Science and Technology of China, Hefei 230026, China
3
School of Electrical and Photoelectric Engineering, West Anhui University, Lu’an 237012, China
*
Authors to whom correspondence should be addressed.
Sensors 2022, 22(12), 4651; https://doi.org/10.3390/s22124651
Submission received: 6 April 2022 / Revised: 14 June 2022 / Accepted: 15 June 2022 / Published: 20 June 2022
(This article belongs to the Topic Human–Machine Interaction)

Abstract

:
During lower-extremity rehabilitation training, muscle activity status needs to be monitored in real time to adjust the assisted force appropriately, but it is a challenging task to obtain muscle force noninvasively. Mechanomyography (MMG) signals offer unparalleled advantages over sEMG, reflecting the intention of human movement while being noninvasive. Therefore, in this paper, based on MMG, a combined scheme of gray relational analysis (GRA) and support vector regression optimized by an improved cuckoo search algorithm (ICS-SVR) is proposed to estimate the knee joint extension force. Firstly, the features reflecting muscle activity comprehensively, such as time-domain features, frequency-domain features, time–frequency-domain features, and nonlinear dynamics features, were extracted from MMG signals, and the relational degree was calculated using the GRA method to obtain the correlation features with high relatedness to the knee joint extension force sequence. Then, a combination of correlated features with high relational degree was input into the designed ICS-SVR model for muscle force estimation. The experimental results show that the evaluation indices of the knee joint extension force estimation obtained by the combined scheme of GRA and ICS-SVR were superior to other regression models and could estimate the muscle force with higher estimation accuracy. It is further demonstrated that the proposed scheme can meet the need of muscle force estimation required for rehabilitation devices, powered prostheses, etc.

1. Introduction

The central nervous system quantitatively controls force production of the skeletal muscles through the successive recruitment of motor units (MUs) [1]. These forces contribute to the generation of the forces required to perform various movements and to interact with the external environment. The skeletal muscles, as the power source of the motor system, work together with the bones and joints to accomplish basic limb movements such as standing, sitting, walking, and jumping, as well as complex movements, under the innervation of the nervous system. Although these behaviors do not directly affect people’s survival, they are directly related to living independently and autonomously.
In the real world, common muscular and neurological disorders, such as stroke and hemiplegia, as well as cumulative diseases such as arthritis, lead to movement problems in the skeletal muscle system and severely affect the body’s ability to move freely [2]. Rehabilitation with assistive devices has received increasing attention for helping these patients regain their ability to live independently and improve their standard of living.
Some rehabilitation auxiliary devices in structured environments are able to obtain kinematic information about the human body during limb movement through sensors or the mechanical principle of the device itself to achieve simple repetitive tasks [3]. However, these devices in unstructured environments should be augmented with force-sensing functions to improve flexibility and interaction capabilities; otherwise, they will fail.
Currently, the force control and interaction capabilities of these assistive devices lag far behind the expected results [4,5] and do not reflect the user’s movement intention, much less operate flexibly. Most patients do not receive satisfactory active training needs, which plays a key role in the rehabilitation process. Therefore, it is necessary to accurately predict and estimate the human muscle force, which helps these auxiliary devices to not only improve their naturalness and flexibility, but also respond to and assist the user in a timely and appropriate manner [6,7]. Accordingly, their users can receive effective and comfortable training.
The knee joint of the lower limbs is one of the most widely studied objects in rehabilitation, and it determines the movement of the human lower limb. Knee joint force plays a crucial role in completing lower-limb movement and interacting with the external environment. The knee joint force is a biomechanical measure used to infer structural loads on the knee joint. The timely and accurate estimation of the lower-limb force information is expected to be used to determine the output power and force information of the auxiliary equipment, as well as to determine more flexible control.
Neural information is not only complex and time-varying, but also has a nonlinear mapping relationship with human movement intention, which poses a challenge in terms of intention decoding and subsequent human–computer interaction [8]. Compared with the traditional muscle force measurement method, human surface biosignals, such as electromyography (sEMG) and mechanomyography (MMG), are easy to measure and do not damage the human body during the measurement process. Therefore, the estimation of muscle force based on human surface biosignals has received more and more attention from researchers and has been widely used in many fields, such as rehabilitation training, intelligent prosthesis control, sports medicine, and human–computer interaction [9,10].
Since the features of sEMG show a nonlinear and complex relationship with muscle force, in recent years, many researchers have started to estimate or predict human kinematic parameters with sEMG using models such as neural network models and support vector machines (SVM), obtaining fruitful results. Peng et al. [11] used two three-layer backpropagation neural networks (BPNNs) to estimate muscle torques at hip and knee joints, on the basis of sEMG, and used them in rehabilitation robots to achieve real-time coordinated active training, obtaining better results. Wang et al. [12] used a cross-model selection (CMS) technology to estimate the handgrip force on the basis of sEMG and obtained a high accuracy. Song et al. [13] used a recurrent artificial neural network (RANN) to estimate the joint torque on the basis of sEMG and also obtained a better prediction accuracy. Even though these studies demonstrated that sEMG obtains outstanding results in estimating or predicting human kinematic parameters, it is worth noting that sEMG sensors can be affected by environmental noise, skin surface cleanliness, and placement, which often lead to problems in terms of the actual test signal instability during sEMG testing. These not only bring difficulties but also further limit their practical application [14,15].
Related studies have shown that, in addition to sEMG signals, MMG signals can be generated during muscle activity, which are generated by the mechanical vibration of muscle fibers [16,17,18]. MMG signals can provide information about the number and firing rates of recruited motor units, reflecting the characteristics of muscle activity [19]. As a result, MMG signals are similar to sEMG and also contain rich information such as muscle state, movement pattern, and movement intention [20]. MMG signals can be detected by various types of sensors that do not need to be placed in precise locations and can also be placed on clothes. Additionally, since MMG signals are mechanical signals, their detection process is also not affected by the change in skin impedance [21]. Compared with sEMG, MMG signals have unparalleled advantages, have attracted extensive attention from many researchers in recent years, and have also been applied in many fields, such as prosthetic hand control [9], muscle fatigue assessment [22], human movement recognition [13], and human kinetic parameter estimation [23,24,25,26,27]. However, because of the characteristics of weak signal strength, low frequency, and strong randomness, it also poses a challenge to decode MMG signals and use them to estimate and predict muscle force. When MMG signals are used for muscle force estimation, there are two critical issues to be addressed: one is to extract the features that can effectively reflect muscle activity; the other is to construct an appropriate model for muscle force estimation.
In recent years, researchers have made valuable attempts in muscle force estimation using MMG signals. Akataki et al. [28] found the MMG/force relationship through the RMS amplitude/force and MPF/force relationships during voluntary isometric ramp contractions of the biceps brachii muscles. Beck et al. [29] found a linear increase in MMG amplitude with isokinetic torque, but no significant change relationship between MMG MPF and isokinetic torque during isokinetic muscle actions of the biceps brachii. With further research, the researchers found the presence of a nonlinear relationship between MMG and muscle force, indicating that MMG can be used to estimate muscle force. Therefore, subsequent researchers also started to extract time-domain features and frequency-domain features from MMG signals and fed them into a machine learning model to achieve a nonlinear MMG–force relationship [25,26,27]. Although some achievements have been made in muscle force estimation, there are still some drawbacks such as fewer features considered and lower accuracy of the estimated models. Recently, to improve the accuracy of muscle force estimation, our research group tried to build an optimized SVR to estimate the knee joint extension force by selecting root mean square (RMS), mean absolute value (MAV), zero crossing (ZC), mean power frequency (MPF), and sample entropy (SE). Although good effects and accuracy were achieved, there are still some shortcomings, such as inefficient models and empirical selection of features that do not fully reflect the muscle activity.
At present, the main features for characterizing MMG signals include time-domain features, frequency-domain features, and time–frequency-domain features. Due to the nonlinear chaotic nature of MMG signals [30], nonlinear dynamics features should also be considered. However, all the features are fed into a machine learning model, which will increase not only the computational complexity but also the computational burden, eventually leading to overfitting of the model and degradation of the estimation accuracy. Therefore, dimension reduction techniques, such as principal component analysis (PCA), independent component analysis (ICA), and linear discriminant analysis (LDA), are widely used to obtain effective features [31,32]. However, these methods tend to destroy the original data structure. Therefore, we need to investigate other methods of feature dimension reduction to avoid destroying the original data structure.
Support vector machine (SVM) is a machine learning method based on statistical learning theory proposed by the research group led by Vapnik in 1995 [33]. With outstanding nonlinear processing capability, high-dimensional mapping ability, and kernel computing technology, the support vector regression (SVR) model has already yielded significant and abundant application results in many fields [34,35]. However, when applying the SVR model in practice, the selection of critical parameters, such as penalty coefficient C and kernel function width parameter σ, can directly affect the prediction accuracy of the model. Usually, trial-and-error and grid search methods are used to select critical parameters, which are not only time-consuming but also ineffective.
Since Beni and Wang first proposed the concept of swarm intelligence in 1989, swarm intelligence algorithms have been continuously developed and used by many researchers to optimize and find the optimal critical parameters in machine learning models such as SVM, GRNN, and BPNN [36,37,38]. Among intelligent algorithms, the cuckoo search (CS) algorithm has attracted the attention of researchers because of its powerful optimization ability and few adjustment parameters [38]. In particular, the application of CS to SVR in many fields has achieved good effects. Yang et al. [39] used CS-SVR to evaluate and monitor the health status of radar equipment and obtained a high recognition rate. Bo et al. [40] used CS-SVM to predict the postoperative recurrence time and recurrence location of patients with liver cancer for studying the mechanism of liver cancer recurrence; the results indicated that CS-SVM could effectively predict the time and location of cancer recurrence, and the prediction results of CS-SVM outperformed those of BPNN. Although the CS algorithm has been applied in many fields with its advantages, it still has some shortcomings, such as uneven population diversity and poor adaptiveness of the search process [41]. Therefore, developing an improved CS (ICS) method to overcome these shortcomings and using ICS to optimize the critical parameters of SVR are beneficial to the solution of practical problems.
As the application of MMG signals in engineering technology continues to be developed, we must consider the influence of effective features reflecting muscle activity and the effects of estimation models when using MMG signals to estimate muscle force. Therefore, in this paper, the MMG feature expression and the adaptability of the model were fully considered, and a scheme based on GRA and the ICS-SVR model was proposed to estimate knee joint extension force. Figure 1 shows the specific process of knee joint extension force estimation based on the proposed scheme, via which the highly accurate estimation of knee joint extension force can be realized. The main contributions of this paper are as follows:
(1)
In this paper, the features reflecting muscle activity were fully considered to avoid the limitations of a single category of features and conventional features. We extracted time-domain features, frequency-domain features, and time–frequency-domain features from MMG signals, as well as nonlinear dynamics features. In order to obtain effective features that are highly correlated with muscle force, GRA was employed for effective feature selection, aiming to achieve high accuracy muscle force estimation using these effective features reflecting muscle activity.
(2)
In this paper, muscle force estimation was based on the SVR model, whose performance depends entirely on critical parameters (C, σ). The cuckoo search (CS) algorithm can optimize the SVR parameters due to the advantages of fast convergence, few parameters, and easy implementation. To obtain better optimization performance, we improved the CS algorithm using a chaotic Tent initialization population and adaptive control parameters. Compared with other optimization algorithms, the optimal global minimum and optimal convergence performance of ICS were obtained in test benchmark functions.
(3)
In this paper, combining the advantages of GRA and ICS-SVR, we designed an MMG–force scheme to effectively obtain muscle activity features and accurately perform muscle force estimation.
The remainder of the paper is arranged as follows: the model building process and related algorithms are given in Section 2; signal acquisition and processing are described in Section 3; experiments and results are given in Section 4; the discussion is given in Section 5; the conclusion is described in Section 6.

2. ICS-SVR

2.1. SVR

SVM, developed by Vapnik [33], is considered one of the most powerful supervised machine learning algorithms based on the structural risk minimization principle. As an effective tool to solve certain problems in complex data processing, SVM can overcome traditional difficulties such as dimensional disaster and overlearning. Compared with neural networks, it has a reliable statistical theory foundation and higher prediction accuracy for high-dimensional small sample data, which is especially suitable for nonlinear problems and has been successfully applied in the fields of classification and prediction.
SVR is a class of support vector machines for fitting regression problems, which has excellent fitting ability, robustness, and fault tolerance for nonlinear and unstable data. On the basis of a previous study, this paper also chose the RBF function. However, the critical parameters of SVR are difficult to select using experience and trial-and-error methods [42]. Traditional parameter search methods such as grid search, particle swarm optimization algorithm (PSO), and genetic optimization algorithm (GA) have the disadvantages of falling into local optimization or slow calculation speed. Therefore, in this paper, we chose the improved swarm intelligence algorithm ICS to search for these critical parameters, so as to avoid the overfitting phenomenon caused by improper parameter selection.

2.2. Cuckoo Search Algorithm

CS is a novel type of heuristic algorithm used to simulate the egg-laying behavior of a cuckoo, which has been extensively applied to solve optimization problems in different fields of engineering [43]. The algorithm has many advantages such as fast convergence, high stability, few parameters, and simple operation, and it is very effective in solving global optimization because it can maintain a balance between local and global random walks using switching parameters. The algorithm is developed from the parasitic reproduction strategy of the cuckoo population itself and Lévy flight behavior. For specific practical applications, the following three ideal assumptions are made [44]:
(1)
Each cuckoo lays one egg in a randomly selected host nest at a time.
(2)
Following the survival of the fittest principle, a strong surviving egg among all the host nests is inherited by the next generation.
(3)
For a fixed number of host nests, the probability of the intruder egg being found by the host bird is Pa ∈ [0, 1].
On the basis of the above assumptions, the position of the cuckoo’s nest can be updated as follows:
X i t + 1 = X i t + α L é v y ( λ ) ,
where X i t + 1 denotes the updated position of the i-th bird’s nest in the t + 1-th iteration, α is the factor of step size used to control the random walk step, the product is a kind of calculation that denotes entry-wise multiplication, and L é v y ( λ ) is a random number, which is drawn from a L é v y distribution,
L é v y ( λ ) ~ μ = t λ ,   1 < λ 3 .
In the implementation, L é v y ( λ ) can be calculated as follows:
L é v y ( λ ) ~ φ μ | v | 1 β ,
where μ and v follow a normal distribution, and β is a distribution factor between 0.3 and 1.99. Here, β = 1.5 ; φ is expressed as follows:
φ = { г ( 1 + β ) s i n ( π β 2 ) г ( 1 + β 2 ) × β × 2 ( β 1 ) / 2 } ,
where г ( · ) is the gamma function.
α is expressed as
α = α 0 ( x i t x b e s t t ) ,
where x b e s t t denotes the optimal bird’s nest of the t-th iteration; α 0 is the control parameter that controls the random step size and usually takes the value of 0.01.
Equation (1) can be reformulated as
X i t + 1 = X i t + α 0 φ μ | v | 1 β ( x i t x b e s t t ) .

2.3. Improved Cuckoo Search Algorithm

Although the CS algorithm is used in many fields with its outstanding advantages, it still has some drawbacks and shortcomings [45]. Firstly, the random initial population lacks uniformity of distribution. Secondly, the optimal solution is greatly affected by the fixed probability P a of the host bird finding the parasitic egg. Thirdly, the L é v y flight mechanism lacks adaptivity due to alternating large and small steps for the global search. Specifically, the lack of adaptive ability in local search will reduce the search efficiency and increase the search time, making it difficult to obtain a balance between precision and search ability.
Consequently, to address these issues, improvements were made in this paper in terms of both population diversity and parameter adaptation, respectively.

2.3.1. Initial Population Chaoticization

Defined as highly unstable/unpredictable motion in finite phase space, chaos often occurs in deterministic nonlinear dynamical systems. According to the change in chaos behavior, chaos has randomness and strong ergodicity. Thus, we proposed a chaotic cuckoo initial population to improve the population diversity through chaotic mapping, while reducing the possibility of falling into local optimum at the initial stage. Tent chaotic mapping was used in this paper, which can be expressed as follows:
x n + 1 = { 2 x n ,   0 x n 1 2 2 ( 1 x n ) ,   1 2 x n 1 .
The specific implementation process is described below.
First, a random x 0 is generated, and then the n nest locations are obtained by iterating using Equation (7).
x k , j ( i ) = [ x k , 1 ( i ) , x k , 2 ( i ) , , x k , n ( i ) ] T ,   k = C , σ ,
where i represents the iteration number, j represents the number of nests, and x ( 0 ,   1 ) .
Then, the chaotic sequence values are mapped to the value range of SVR critical parameters.
c x k , j ( i ) = L b j + x k , j ( i ) ( U b j L b j ) ,
where Lb = 0.001; Ub = 500.

2.3.2. Adaptive Control Parameters

The original CS algorithm adopts a fixed value for both P a and a 0 . Typically, P a is 0.25, and a 0 is 0.01. However, in the face of practical different optimization problems, it is often necessary to adjust these parameters to suit specific needs on the basis of personal experience; hence, so it is a challenging task to choose the appropriate parameters P a and a 0 [44].
Due to the lack of control parameter adaptivity, the algorithm does not guarantee fast convergence. If a 0 is relatively large, it can improve the global search ability of the algorithm, but reduce the local search ability; if a 0 is too small, the local search is better, but the total number of iterations increases, and the convergence may become worse. Therefore, a 0 should gradually decrease with the iterative process of the algorithm. In addition, during the iterative search process, as the quality of cuckoo individuals gradually increases, the intensity of population evolution should be appropriately increased to avoid the algorithm falling into local optimal solutions; thus, P a should also gradually decrease with the iterative process. Therefore, to improve the performance of the CS algorithm, we propose an adaptive control parameter strategy, which is represented as follows:
a 0 = a m i n + 0.5 ( a m a x a m i n ) ( c o s ( l o g ( 1 + ( e 1 ) t T ) π ) + 1 ) ,
P a = P a m i n + ( P a m a x P a m i n ) ( 200 t T ) 2 ,
where t and T are the number of the current iteration and the number of total iterations, respectively. The parameters in the adaptive control parameter strategy were set as follows: a m i n = 0.001; a m a x = 0.2; P a ,   m i n = 0.1; P a , m a x = 0.70; T = 200.

2.4. Architecture of the ICS-SVR Model

The designed ICS can quickly find the optimal solution without searching all the parameter points; accordingly, this method can greatly improve the efficiency of finding the optimal solution in a wide range of SVR parameters C and σ. The optimal parameters obtained can be used to improve the fitting performance of SVR. In this paper, a muscle force estimation model was constructed on the basis of ICS-SVR, and the procedure of ICS-SVR is illustrated in Figure 2.

2.5. Gray Correlation Analysis

Choosing the optimal features is crucial in pattern recognition because it determines the accuracy and generalization ability of classifiers and prediction models [46]. Therefore, it is a challenging task to extract effective features from the MMG signals that can reflect muscle activity. In practical engineering problems, since the data collected by the instrumentation are of great importance, it is not recommended to destroy the original structure of the data during modeling if not absolutely necessary.
Gray system theory was proposed by Deng Julong in 1982 on the basis of the mathematical theory of systems engineering, which is one of the important achievements in the field of uncertainty system research [47]. GRA is an important part of gray system theory, where a gray relationship is an uncertain relationship between two variables. The dominant factors are determined by calculating the correlation degree between multiple factors and the same reference sequence [48]. Therefore, in this paper, GRA was used to select the effective MMG features that are most associated with knee joint extension force. The specific implementation of GRA is described below.
Let F = { F ( k ) | k = 1 , 2 , , n } be the knee joint extension force sequence and X i = { x i ( k ) | k = 1 , 2 , , n ;   i = 1 , 2 , , m }   be the MMG feature sequences.
Relational coefficients between the reference sequence (the knee joint extension force sequence) and comparison sequences (the MMG feature sequences) were calculated as follows:
ξ i ( k ) = min i min k | F ( k ) x i ( k ) | + ρ max i max k | F ( k ) x i ( k ) | | F ( k ) x i ( k ) | + ρ max i max k | F ( k ) x i ( k ) | .
Let Δ i ( k ) = | F ( k ) x i ( k ) | ; then, the difference sequence is
Δ i ( k ) = ( Δ i ( 1 ) , Δ i ( 2 ) , , Δ i ( n ) ) .
After simplifying Equation (12), we can get
ξ i ( k ) = min i min k Δ i ( k ) + ρ max i max k Δ i ( k ) Δ i ( k ) + ρ max i max k Δ i ( k ) .
A small resolution coefficient (ρ) indicates a great resolution. In general, the range of ρ is (0, 1). The resolution is best when ρ ≤ 0.5463; usually, ρ = 0.5.
The relational coefficient represents the value of the correlation degree between the knee joint extension force sequence and the MMG feature sequence.
Usually, the average of the correlation coefficient at each moment is quantitatively expressed as the relational degree between the knee joint extension force sequence and the MMG feature sequence. The relational degree r i is calculated as follows:
r i = 1 n k = 1 n ξ i ( k ) ,   k = 1 , 2 , , n .

2.6. Performance of the Models

The performance evaluation indicators for the estimation models were the root-mean-square error (RMSE), mean absolute percentage error (MAPE), and correlation coefficient (R), which are defined as follows:
RMSE = 1 N i = 1 N ( y ^ i y i ) 2 ,
MAPE = 1 N i = 1 N | y ^ i y i y i | ,
R = C o v ( y ^ i , y i ) D ( y ^ i ) D ( y i ) ,
where y ^ i is the estimated value, y i is the actual value, D ( · ) is the calculated variance, C o v ( · ) is the covariance, and N is the number of samples in the test set. Generally speaking, the closer RMSE and MAPE are to 0 and R is to 1, the closer the estimated value of the model is to the actual value.

3. Signal Acquisition and Preprocessing

3.1. Experimental Procedure and Signal Processing

A total of five healthy male subjects aged 22–25 years without neuromuscular and musculoskeletal diseases were recruited and briefed on the experiment and the potential associated risks. Informed consent forms were signed by the subjects before experimentation. The subjects were asked to sit comfortably in a test chair with their right leg fixed and bent at a 90° angle. The MMG and the force signals were recorded using a DAQ device (VK702, VKinging, Inc., Shenzhen, China) at a sampling frequency of 1 kHz. MMG signals came from the acceleration sensor (ADXL335, Analog Devices, Inc., Norwood, MA, USA), while the force signals came from the force sensor (DYLF-30, 0–300 N, Bengbu Sensor, Inc., Bengbu, China). Three accelerometers were placed on the clothes and bound with elastic bandages to the muscle belly positions of the rectus femoris (RF), vastus lateralis (VL), and vastus medialis (VM) muscles, respectively. The data collection process was completed across two days. On the first day, subjects were familiarized with the experimental process and precautions, and maximal voluntary contraction force (MVC) was collected. On the second day, forces of 10–80% MVC and corresponding MMG signals were collected. Sufficient rest time (2–5 min) was provided between experiments to avoid muscle fatigue.
To attenuate the effects of movement and noise, the MMG signals were filtered by a fourth-order Butterworth filter with a passband of 5–100 Hz [49], and the force signals were filtered by a third-order Butterworth filter with a low passband of 2 Hz. Due to the fluctuation of force in the process of data collection, 6 s of stable muscle force data, whose fluctuation range was less than 5% of the relative value, were selected as the experimental data.

3.2. Feature Extraction

In this paper, sliding windows were used for feature extraction of MMG signals. The window length was 1000 data points, and the overlap length was 100 data points.
Feature extraction plays a key role in pattern recognition, which directly determines the accuracy of prediction and estimation. Therefore, in this paper, we extracted not only time-domain features, frequency-domain features, and time–frequency-domain features, but also nonlinear dynamics features. Specifically, root-mean-square (RMS), kurtosis, standard deviation (SD), slope sign change (SSC), mean absolute values (MAV), zero crossings (ZC), and waveform length (WL) were selected as time-domain features. Mean power frequency (MPF) and median frequency (MDF) were selected as the frequency-domain features. After a three-layer wavelet packet decomposition, wavelet packet energy (WPE) and the energy of every frequency band (WP1 WP2, WP3, WP4, WP5, WP6, WP7, WP8) were calculated as time–frequency-domain features. Lempel–Ziv complexity (LZC), sampling entropy (SampEn), approximate entropy (ApEn), fuzzy entropy (FuzzyEn), distribution entropy (DistEn), box-counting fractal dimension (FD), and largest Lyapunov exponent (LyapExp) were calculated as nonlinear dynamics features. The details of the specific nonlinear dynamics feature extraction method can be found in [50,51,52,53,54]. As a result, we extracted 25 dimensional features from each segment, with a total of 75 dimensional features for the three channels.
Additionally, for eight different strength forces, a sample dataset was obtained for each subject, resulting in a total of 400 features. To further train and test the model, the 400 samples were randomly shuffled; 90% of them were selected as the training set and 10% were selected as the test set.

3.3. Data Normalization

To eliminate the impact of the value range or uneven magnitude of different features, the MMG features and muscle force were normalized using the Z-score method, which can be expressed as follows:
z = x x ¯ σ ,
where x ¯ = 1 n i = 1 n x , and σ = 1 n 1 i = 1 n ( x x ¯ ) . The overfitting of the ICS-SVR model can be avoided to some extent by normalizing the data.

4. Experiments and Results

In this paper, our experiments consisted of three parts. In the first part, test benchmark functions were used to test the performance of different swarm intelligence algorithms. If ICS was shown to have better search performance than other algorithms, it could be used to optimize SVR critical parameters. In the second part, the relational degree between MMG features and knee joint extension force was calculated to obtain different feature combinations, which were input into the ICS-SVR model to obtain the estimated knee extension force. Then, the optimal feature combination sequence was obtained by comparing and analyzing the estimation accuracy under different feature combinations. In the third part, the ICS-SVR model was compared with other classical models in the estimation of the knee joint extension force. If the estimation results of the model were found to be closer to real data than other models, the model could be considered the best.

4.1. Performance Analysis of ICS Algorithm

To verify the effectiveness of the designed ICS algorithm, a comparative simulation experiment was designed to test the optimization performance.
Performance analysis was performed on the ICS algorithm and other algorithms, such as particle swarm optimization (PSO), gray wolf optimizer (GWO), and the original CS algorithm, using the four global optimal test benchmark functions without constraints. Parameter values in the four optimization algorithms are shown in Table 1. The four test benchmark functions are expressed below.
(1)
Rosenbrock function:
f 1 = i = 1 n 1 [ 100 ( x i + 1 x i 2 ) 2 + ( x i 1 ) 2 ] , 2.048 x i 2.048 .
(2)
Griewank function:
f 2 = i = 1 n x i 2 4000 i = 1 n cos ( x i i ) + 1 , 600 x i 600 .
(3)
Cross-in-tray function:
f 3 = 0.0001 ( | sin ( x i ) sin ( x i + 1 ) exp ( | 100 x i 2 + x i + 1 2 π | ) | + 1 ) 0.1 , 10 x i 10 .
(4)
Schaffer function:
f 4 = 0.5 + ( sin x i 2 + x i + 1 2 ) 2 0.5 [ 1 + 0.01 ( x i 2 + x i + 1 2 ) ] 2 , 10 x i 10 .
The Rosenbrock function, also known as the valley or banana function, has a global minimum f 1 ( X * ) = 0 at X * = ( 1 , , 1 ) . The Griewank function has many widespread and regularly distributed local minima and a global minimum f 2 ( X * ) = 0 at X * = ( 0 , , 0 ) . The cross-in-tray function has a global minimum f 3 ( X * ) = 2.06261 at four points of X * = ( ± 1.3491 ,   ± 1.3491 ) . The Schaffer function has a strongly oscillating behavior, with many local minima and a global minimum f 4 ( X * ) = 0 at X * = ( 0 , , 0 ) [45,55,56]. The search space of the four test benchmark functions is shown in Figure 3.
The purpose of the test is to find the global minimum of all four functions. The PSO algorithm, the GWO algorithm, the original CS algorithm, and the ICS algorithm were terminated when the number of iterations reached 200. The convergence curves are shown in Figure 4. As can be seen from Figure 4, the ICS algorithm could converge to or was much closer to the optimal solution in a shorter time and was more accurate than the other algorithms in the four test benchmark functions. That is to say, the proposed new algorithm could effectively improve the convergence speed and convergence accuracy in search attempts.
To eliminate the difference in each experiment, the four algorithms were executed 50 times each. The numerical results of the selected indices, including the optimal solution, worst solution, average solution, and SD, are presented in Table 2.
Comparing the results presented in Table 2, it can be obviously seen that they almost always succeeded in finding the global minimum; however, the ICS algorithm was superior to PSO, GWO, and CS in all indices. In general, it can be concluded that the ICS algorithm contributed to superior performance, while CS and GWO were better than PSO, and PSO performed the worst. Therefore, the designed ICS algorithm was used in this paper to optimize the critical parameters of SVR.

4.2. Feature Combination Sequence Selection with GRA

A total of 25 features were extracted from the MMG signal segment of each channel, and a total of 75 features were obtained from the three channels. Among these features, there were many unrelated or low-correlation features, and inputting all of these features into the model would lead to overfitting and a heavy computational burden. Therefore, in this paper, GRA was used for MMG feature selection to extract features exhibiting a high correlation with knee joint extension force. Firstly, based on GRA, the relational degree between knee joint extension force and various MMG features was calculated, and the results are shown in Figure 5. Usually, the classification of the relational degree is as follows: 0.8–1 indicates a strong correlation, 0.6–0.8 indicates a general correlation, and 0.6–0 indicates a weak correlation. To obtain the optimal features to describe the knee extension force, we selected the features with a correlation of 0.6–1 as feature combination A, features with a correlation of 0.75–1 as feature combination B, features with a correlation of 0.8–1 as feature combination C, features with a correlation of 0.85–1 as feature combination D, and features with a correlation of 0.9–1 as feature combination E, and we input them into the ICS-SVR model for knee joint extension force estimation; the results are presented in Table 3.
Generally, a smaller value of RMSE and MAPE and a larger value of R indicate a better effect of the feature combination of knee joint extension force estimation. From Table 3, it can be seen that feature combination D, i.e., the feature combination with a relational degree greater than 0.85, had the best results for knee joint extension force estimation with the smallest RMSE, the smallest MAPE, and the largest R. Therefore, the best estimation effect was obtained by selecting feature combination D, i.e., the selected MMG feature sequence was 2, 19, 20, 25, 27, 44, 45, 50, 51, 53, 54, 57, 61, and 71. This also indicates that different feature combinations had a certain influence on the estimation of knee joint extension force, and some features did not correlate well with the estimated results.
As a result, in the process of muscle force estimation, the feature sequence of feature combination D was selected in this paper to avoid the overfitting of the model caused by too many input features. Figure 6 shows the results of muscle force estimation for subject S1 with feature combination D. As can be seen in Figure 6, the estimated values were mostly consistent with the actual observed values, indicating that the proposed method is effective.
To further verify the validity of the proposed method, muscle force estimation was performed in five subjects. The results of the knee joint extension force estimation for the five subjects are shown in Table 4. The results indicate that muscle force estimation using feature combination D in the five subjects experiment outperformed the effect of the other feature combinations.

4.3. Comparative Performance of the Proposed Model with Classical Machine Learning Models

To further demonstrate the validity of the ICS-SVR model, a comparison was made with BPNN, ELM, and CS-SVR. The network structures of BPNN and ELM were 14–10–1 and 14–8–1, respectively. The estimated results of subject S1 upon executing the different models 10 times are shown in Table 5. As shown in Table 5, the proposed model outperformed BPNN, ELM, and CS-SVR in terms of estimation effect. Additionally, compared with the fluctuation of their estimated results, the proposed model had the least fluctuation and the most stable estimation results.
Furthermore, different models were used to estimate muscle force for the five subjects. The results of knee joint extension force estimation are shown in Figure 7 and Table 6. As can be seen in Figure 7, the force estimation effects of each subject with ICS-SVR were superior to those with BPNN, ELM, and CS-SVR. In addition, using the proposed model for muscle force estimation, all five subjects obtained outstanding average results with the highest R values and the lowest RMSE and MAPE, as shown in Table 6. In particular, compared with BPNN, ELM, and CS-SVR models, the estimated results of the proposed model had less fluctuation. Meanwhile, in terms of numerical results, the accuracy of the muscle force estimation of the proposed model in this paper was optimal compared with the results of the literature [25,26,27,28,29]. Furthermore, it is shown that our proposed model had unparalleled applicability and prediction accuracy in knee joint extension force estimation.

5. Discussion

The human body often needs to exert force during various movements, which play an important role in daily activities. To achieve flexible human–computer interaction during effective training with rehabilitation devices, this study aimed to decode human movement intentions through muscle interfaces to achieve the estimation of interaction forces during isometric contractions. Muscle force estimation can be used to improve the naturalness and flexibility of human–machine coordinated movement. Many previous studies used sEMG signals for human movement classification and kinetic parameter estimation. However, sEMG detection has some disadvantages compared with MMG signals, which brings many inconveniences to practical applications. Therefore, decoding human intentions using MMG signals can provide a more natural and portable human–computer interaction, which has attracted the attention of a wide range of researchers. Furthermore, MMG signals can objectively reflect muscle activity and are expected to be widely used in rehabilitation equipment, powered prosthesis, etc.
The muscle contraction part of the kinematic chain through the limb effector can generate forces that interact with the environment [57]. Due to the anthropotomy and motion control principle, the excitation–contraction of muscles will be dynamic and coupled. Therefore, it is challenging to estimate muscle force accurately and efficiently. In this paper, the relationship between knee joint extension force and MMG was established through GRA and ICS-SVR. Our proposed combined scheme of GRA and ICS-SVR could better deal with the difficulty of MMG feature selection and achieve accurate knee joint extension force estimation.
For the estimated models, if excessive MMG features are input, the training speed and efficiency of the estimation models would be reduced, and the convergence may not be obtained. Therefore, it is worth noting that a suitable feature sequence for the force estimation is an issue that must be considered. The main principle in selecting MMG features is the effectiveness of the force estimation. The relational degree between the knee joint extension force sequence and the MMG feature sequence calculated by GRA facilitates the easy selection of related features that are highly correlated with the reference sequence, thus effectively avoiding the difficulty of feature selection and the long time needed to test the selected features. Therefore, GRA can be performed on the MMG features to fulfill the above effectiveness, avoiding the destruction of the data and maintaining the physical meaning of the original data. From the experimental results, it can be concluded that the optimal feature combination obtained using GRA had a better estimation effect than all feature combinations, whereby RMSE was reduced by 25.40%, MAPE was reduced by 13.83%, and MAPE was improved by 0.03%.
According to the findings of a previous study, SVR is an effective method for muscle force estimation, but its effect is extremely susceptible to critical parameters. Therefore, we proposed an ICS method to optimize the critical parameters of SVR. Compared with PSO, GWO, and CS, the designed ICS algorithm was optimal in terms of search performance. The adaptability of the CS algorithm was enhanced by using a chaotic optimization initial population and adaptive control parameters, such that the critical parameters searched became more suitable for MMG characteristics. Subsequently, to demonstrate the validity of the model, we compared the proposed ICS-SVR with classical machine learning algorithms, such as BPNN, ELM, and CS-SVR, as a function of specific knee joint extension force estimation cases. The results showed that the proposed ICS-SVR model performed optimally in knee joint extension force estimation. In addition, the overall average results were about the same compared to the pre-IGWO-SVR model results, but obtained a nice performance improvement in RMSE, from 0.2492 to 0.2205. The features in the previous study were selected empirically, while the present study was based on GRA for automatic selection, which is more adaptable to complex time-series MMG signals; this result can perhaps bring immeasurable effect in subsequent muscle force estimation under complex conditions. Further research is needed to apply the proposed scheme to the muscle force estimation under complex conditions based on MMG signals and to the estimation of movement parameters and muscle force based on other time-series signals, such as sEMG and EEG.
MMG can be used as a useful tool to examine various conditions of muscle activity. Consequently, it can be applied in many fields such as muscle function assessment, power prosthesis, and medical rehabilitation [5,10]. Usually, muscle atrophy or muscle damage may occur as a result of injury or disease, and muscle weakness may occur. The main purpose of rehabilitation training is to restore normal function to the atrophied muscles. To facilitate an effective rehabilitation process, the relevant equipment should provide the appropriate training intensity according to the patient’s muscle activity to provide the required assistance. Nevertheless, one of the most important training processes is muscle isometric contraction. To make this process more effective, MMG can be used to monitor the degree of muscle activation and muscle force interaction with the external environment, which can help to provide an understanding of muscle function and recovery effects, and further provide a reasonable training program.
In this paper, a high estimation accuracy was achieved in the knee extension force estimation. In addition to the fact that this result only applies to the muscle isometric contraction case, this study had some limitations, such as limited sample size and single physical health status.

6. Conclusions

In this paper, a combined scheme of GRA and ICS-SVR was developed to accurately estimate muscle force from MMG signals. The method screened out the features reflecting muscle activity by GRA and reduced the input features of the model, thereby effectively improving the efficiency of the muscle force estimation model. Meanwhile, the CS algorithm was improved by a chaotic optimization initial population and adaptive control parameters, resulting in a further improvement of the global search performance. Finally, ICS was applied to the critical parameter search of SVR to obtain excellent muscle force estimation effects. In particular, the mean values of RMSE and MAPE obtained by ICS-SVR were significantly lower than those obtained by BPNN, ELM, and CS-SVR, and the mean value of R was closer to 1. In conclusion, the proposed combined scheme of ICS-SVR and GRA is a feasible, effective, and promising method for the estimation of knee joint extension force, which is also expected to further improve the working performance of related devices. Future research will be focused on experiments on dynamic force generation during isometric muscle contraction and force generation during joint movement to further refine key interaction techniques required for specific rehabilitation and powered prostheses.

Author Contributions

Conceptualization, L.G. and D.W.; methodology, Z.L.; software, Z.L.; validation, L.G., W.L. and D.W.; formal analysis, Z.L. and W.L.; investigation, L.G. and H.C.; writing—original draft preparation, Z.L.; writing—review and editing, Z.L., W.L. and D.W.; visualization, H.C., W.L. and G.Z.; supervision, L.G.; project administration, D.W.; funding acquisition, L.G., H.C. and Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDA22040303), the Key Scientific Research Projects of Anhui Province Higher Education (Grant No. KJ2020A0630), the University Synergy Innovation Program of Anhui Province (Grant No. GXXT-2021-026), and the Key Research Projects supported by the National Natural Science Foundation of China (Grant No. 92067205), as well as partially funded by the HFIPS Director’s Fund (Grant No.YZJJ2021QN25).

Institutional Review Board Statement

Ethical review and approval were approved by the Institutional Review Committee of Hefei Institute of Physical Science, Chinese Academy of Sciences. The approval number was No. SWYX-Y-2021-75. All subjects were students and staff of the research group, and no patients were involved. Moreover, there was no harm to the body in our experiment.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data used to support the findings of this study are included within the article.

Acknowledgments

This work was partly carried out at the Robotic Sensors and Human–Computer Interaction Laboratory, Hefei Institude of Intelligent Machinery, Chinese Academy of Sciences.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rohlén, R.; Stålberg, E.; Stöverud, K.; Yu, J.; Grönlund, C. A Method for Identification of Mechanical Response of Motor Units in Skeletal Muscle Voluntary Contractions Using Ultrafast Ultrasound Imaging—Simulations and Experimental Tests. IEEE Access 2020, 8, 50299–50311. [Google Scholar] [CrossRef]
  2. Hou, J.; Sun, Y.; Sun, L.; Pan, B.; Huang, Z.; Wu, J.; Zhang, Z. A Pilot Study of Individual Muscle Force Prediction during Elbow Flexion and Extension in the Neurorehabilitation Field. Sensors 2016, 16, 2018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Wang, B.; Yang, C.; Xie, Q. Human-machine interfaces based on EMG and Kinect applied to teleoperation of a mobile humanoid robot. In Proceedings of the 10th World Congress on Intelligent Control and Automation (WCICA), Beijing, China, 6–8 July 2012; pp. 3903–3908. [Google Scholar]
  4. Ehrsson, H.H.; Fagergren, A.; Jonsson, T.; Westling, G.; Johansson, R.S.; Forssberg, H. Cortical activity in precision- versus power-grip tasks: An fMRI study. J. Neurophysiol. 2000, 83, 528–536. [Google Scholar] [CrossRef] [PubMed]
  5. Uwamahoro, R.; Sundaraj, K.; Subramaniam, I.D. Assessment of muscle activity using electrical stimulation and mechanomyography: A systematic review. Biomed. Eng. Online 2021, 20, 1. [Google Scholar] [CrossRef]
  6. Lenzi, T.; De Rossi, S.M.M.; Vitiello, N.; Carrozza, M.C. Intention-based EMG control for powered exoskeletons. IEEE Trans. Biomed. Eng. 2012, 59, 2180–2190. [Google Scholar] [CrossRef]
  7. Villa-Parra, A.C.; Delisle-Rodriguez, D.; López-Delis, A.; Bastos-Filho, T.; Sagaró, R.; Frizera-Neto, A. Towards a robotic knee exoskeleton control based on human motion intention through EEG and sEMG signals. Procedia Manuf. 2015, 3, 1379–1386. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, Q.; Fang, L.; Zhang, Q.; Xiong, C. Simultaneous estimation of joint angle and interaction force towards sEMG-driven human-robot interaction during constrained tasks. Neurocomputing 2022, 484, 38–45. [Google Scholar] [CrossRef]
  9. Staudenmann, D.; Roeleveld, K.; Stegeman, D.F.; van Dieën, J.H. Methodological aspects of SEMG recordings for force estimation—A tutorial and review. J. Electromyogr. Kinesiol. 2010, 20, 375–387. [Google Scholar] [CrossRef]
  10. Islam, M.A.; Sundaraj, K.; Ahmad, R.B.; Ahamed, N.U. Mechanomyogram for muscle function assessment: A review. PLoS ONE 2013, 8, e58902. [Google Scholar] [CrossRef]
  11. Peng, L.; Hou, Z.G.; Kasabov, N.; Hu, J.; Peng, L.; Wang, W.Q. sEMG-based torque estimation for robot-assisted lower limb rehabilitation. In Proceedings of the International Joint Conference on Neural Networks (IJCNN), Killarney, Ireland, 12–17 July 2015; pp. 1–5. [Google Scholar]
  12. Wang, K.; Zhang, X.; Ota, J.; Huang, Y. Development of a SEMG-Handgrip Force Model Based on Cross Model Selection. IEEE Sens. J. 2019, 19, 1829–1838. [Google Scholar] [CrossRef]
  13. Song, R.; Tong, K.Y. Using recurrent artificial neural network model to estimate voluntary elbow torque in dynamic situations. Med. Biol. Eng. Comput. 2005, 43, 473–480. [Google Scholar] [CrossRef] [PubMed]
  14. Asheghabadi, A.S.; Moqadam, S.B.; Xu, J. Multichannel Finger Pattern Recognition Using Single-Site Mechanomyography. IEEE Sens. J. 2021, 21, 8184–8193. [Google Scholar] [CrossRef]
  15. Vidovic, M.M.C.; Hwang, H.J.; Amsüss, S.; Hahne, J.M.; Farina, D.; Müller, K.R. Improving the Robustness of Myoelectric Pattern Recognition for Upper Limb Prostheses by Covariate Shift Adaptation. IEEE Trans. Neural Syst. Rehabil. 2016, 24, 961–970. [Google Scholar] [CrossRef] [PubMed]
  16. Oster, G.; Jaffe, J.S. Low frequency sounds from sustained contraction of human skeletal muscle. Biophys. J. 1980, 30, 119–127. [Google Scholar] [CrossRef] [Green Version]
  17. Talib, I.; Sundaraj, K.; Lam, C.K. Choice of mechanomyography sensors for diverse types of muscle activities. IJEECE 2018, 10, 79–82. [Google Scholar]
  18. Beck, T.W.; Housh, T.J.; Johnson, G.O.; Weir, J.P.; Cramer, J.T.; Coburn, J.W.; Malek, M.H. Mechanomyographic amplitude and mean power frequency versus torque relationships during isokinetic and isometric muscle actions of the biceps brachii. J. Electromyogr. Kinesiol. 2004, 14, 555–564. [Google Scholar] [CrossRef]
  19. Beck, T.W.; Housh, T.J.; Johnson, G.O.; Cramer, J.T.; Weir, J.P.; Coburn, J.W.; Malek, M.H. Does the frequency content of the surface mechanomyographic signal reflect motor unit firing rates? A brief review. J. Electromyogr. Kinesiol. 2007, 17, 1–13. [Google Scholar] [CrossRef]
  20. Islam, M.A.; Sundaraj, K.; Ahmad, R.B.; Ahamed, N.U.; Ali, M.A. Mechanomyography Sensor Development, Related Signal Processing, and Applications: A Systematic Review. IEEE Sens. J. 2013, 13, 2499–2516. [Google Scholar] [CrossRef]
  21. Wu, H.; Wang, D.; Huang, Q.; Gao, L. Real-time continuous recognition of knee motion using multi-channel mechanomyography signals detected on clothes. J. Electromyogr. Kinesiol. 2018, 38, 94–102. [Google Scholar] [CrossRef]
  22. Naeem, J.; Hamzaid, N.A.; Islam, M.A.; Azman, A.W.; Bijak, M. Mechanomyography-based muscle fatigue detection during electrically elicited cycling in patients with spinal cord injury. Med. Biol. Eng. Comput. 2019, 57, 1199–1211. [Google Scholar] [CrossRef]
  23. Hondo, N.; Tsuji, T. Torque Estimation of Knee Flexion and Extension Movements From a Mechanomyogram of the Femoral Muscle. IEEE Trans. Neural. Syst. Rehabil. Eng. 2022, 30, 1120–1126. [Google Scholar] [CrossRef] [PubMed]
  24. Ibitoye, M.O.; Hamzaid, N.A.; Abdul Wahab, A.K.; Hasnan, N.; Olatunji, S.O.; Davis, G.M. Estimation of Electrically-Evoked Knee Torque from Mechanomyography Using Support Vector Regression. Sensors 2016, 16, 1115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Youn, W.; Kim, J. Estimation of elbow flexion force during isometric muscle contraction from mechanomyography and electromyography. Med. Biol. Eng. Comput. 2010, 48, 1149–1157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Youn, W.; Kim, J. Feasibility of using an artificial neural network model to estimate the elbow flexion force from mechanomyography. J. Neurosci. Methods 2011, 194, 386–393. [Google Scholar] [CrossRef]
  27. Lei, K.F.; Cheng, S.C.; Lee, M.Y.; Lin, W.Y. Measurement and estimation of muscle contraction strength using mechanomyography based on artificial neural network algorithm. Biomed. Eng. Appl. Basis Commun. 2013, 25, 1350020. [Google Scholar] [CrossRef]
  28. Akataki, K.; Mita, K.; Watakabe, M.; Itoh, K. Mechanomyogram and force relationship during voluntary isometric ramp contractions of the biceps brachii muscle. Eur. J. Appl. Physiol. 2001, 84, 19–25. [Google Scholar] [CrossRef]
  29. Beck, T.W.; Housh, T.J.; Johnson, G.O.; Weir, J.P.; Cramer, J.T.; Coburn, J.W.; Malek, M.H. Mechanomyographic and electromyographic time and frequency domain responses during submaximal to maximal isokinetic muscle actions of the biceps brachii. Eur. J. Appl. Physiol. 2004, 92, 352–359. [Google Scholar] [CrossRef]
  30. Xie, H.B.; Guo, J.Y.; Zheng, Y.P. Uncovering chaotic structure in mechanomyography signals of fatigue biceps brachii muscle. J. Biomech. 2010, 43, 1224–1226. [Google Scholar] [CrossRef]
  31. Li, C.; Ren, J.; Huang, H.; Wang, B.; Zhu, Y.; Hu, H. PCA and deep learning based myoelectric grasping control of a prosthetic hand. Biomed. Eng. Online 2018, 17, 107. [Google Scholar] [CrossRef] [Green Version]
  32. Subasi, A.; Gursoy, M.I. EEG signal classification using PCA, ICA, LDA and support vector machines. Expert Syst. Appl. 2010, 37, 8659–8666. [Google Scholar] [CrossRef]
  33. Cortes, C.; Vapnik, V. Support-Vector Networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  34. Zhang, X.; Ding, S.; Xue, Y. An improved multiple birth support vector machine for pattern classification. Neurocomputing 2017, 225, 119–128. [Google Scholar] [CrossRef]
  35. Hua, X.; Ding, S. Weighted least squares projection twin support vector machines with local information. Neurocomputing 2015, 160, 228–237. [Google Scholar] [CrossRef]
  36. Yi, J.H.; Xu, W.H.; Chen, Y.T. Novel back propagation optimization by Cuckoo Search algorithm. Sci. World J. 2014, 2014, 878262. [Google Scholar] [CrossRef] [Green Version]
  37. Wang, F.; Yin, T.; Lei, C.; Zhang, Y.; Wang, Y.; Liu, J. Prediction of lower limb joint angle using sEMG based on GA-GRNN. In Proceedings of the IEEE International Conference on Cyber Technology in Automation, Control, and Intelligent Systems (CYBER), Shenyang, China, 8–12 June 2015; pp. 1894–1899. [Google Scholar]
  38. Wu, Q.; Zhao, X. Traffic time series prediction based on CS and SVR. In Proceedings of the Chinese Control and Decision Conference (CCDC), Shenyang, China, 9 July 2018; pp. 3427–3432. [Google Scholar]
  39. Yang, Y.; Zhang, M.; Dai, Y. A Fuzzy Comprehensive CS-SVR Model-based health status evaluation of radar. PLoS ONE 2019, 14, e0213833. [Google Scholar] [CrossRef] [Green Version]
  40. Bo, J.; Li, S.; Ma, P.; Zhu, Y.; Zhang, Y. Research on Early Warning Mechanism and Model of Liver Cancer Rehabilitation Based on CS-SVM. J. Healthc. Eng. 2021, 12, 6658776. [Google Scholar]
  41. Yang, X.; Deb, S. Cuckoo search: Recent advances and applications. Neural Comput. Appl. 2014, 24, 169–174. [Google Scholar] [CrossRef] [Green Version]
  42. Zhou, Z.; Zhang, R.; Wang, Y.; Zhu, Z.; Zhang, J. Color difference classification based on optimization support vector machine of improved grey wolf algorithm. Optik 2018, 170, 17–29. [Google Scholar] [CrossRef]
  43. Yang, X.S.; Deb, S. Cuckoo search via Levy flights. In Proceedings of the World Congress on Nature & Biologically Inspired Computing (NaBIC), Coimbatore, India, 9–11 December 2009; pp. 210–214. [Google Scholar]
  44. Sheng, Z.; Wang, J.; Zhou, S.; Zhou, B. Parameter estimation for chaotic systems using a hybrid adaptive cuckoo search with simulated annealing algorithm. Chaos 2014, 24, 013133. [Google Scholar] [CrossRef]
  45. Du, X.; Wang, J.; Jegatheesan, V.; Shi, G. Parameter estimation of activated sludge process based on an improved cuckoo search algorithm. Bioresour. Technol. 2018, 249, 447–456. [Google Scholar] [CrossRef]
  46. Alves, N.; Chau, T. Classification of the mechanomyogram: Its potential as a multifunction access pathway. In Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Minneapolis, MN, USA, 13 November 2009; pp. 2951–2954. [Google Scholar]
  47. Deng, J.L. Control problems of grey systems. Syst. Control Lett. 1982, 1, 288–294. [Google Scholar]
  48. Wang, H.; Wang, L.; Yang, L.; Shi, X.; Wen, Z.; Dong, X. Exploring the relationship between the dielectric properties and viability of human normal hepatic tissues from 10 Hz to 100 MHz based on grey relational analysis and BP neural network. Comput. Biol. Med. 2021, 134, 104494. [Google Scholar] [CrossRef] [PubMed]
  49. Zuniga, J.M.; Housh, T.J.; Camic, C.L.; Hendrix, C.R.; Mielke, M.; Schmidt, R.J.; Johnson, G.O. The effects of accelerometer placement on mechanomyographic amplitude and mean power frequency during cycle ergometry. J. Electromyogr. Kinesiol. 2010, 20, 719–725. [Google Scholar] [CrossRef] [PubMed]
  50. Torres, A.; Fiz, J.A.; Jane, R.; Laciar, E.; Galdiz, J.B.; Gea, J.; Morera, J. Rényi entropy and Lempel-Ziv complexity of mechanomyographic recordings of diaphragm muscle as indexes of respiratory effort. In Proceedings of the 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Vancouver, BC, Canada, 20–25 August 2008; pp. 2112–2115. [Google Scholar]
  51. Lima, C.A.M.; Coelho, A.L.V.; Madeo, R.C.B.; Peres, S.M. Classification of electromyography signals using relevance vector machines and fractal dimension. Neural Comput. Appl. 2016, 27, 791–804. [Google Scholar] [CrossRef]
  52. Li, P.; Karmakar, C.; Yan, C.; Palaniswami, M.; Liu, C. Classification of 5-S Epileptic EEG Recordings Using Distribution Entropy and Sample Entropy. Front. Physiol. 2016, 7, 136. [Google Scholar] [CrossRef] [Green Version]
  53. Xue, Y.; Ji, X.; Zhou, D.; Li, J.; Ju, Z. SEMG-Based Human In-Hand Motion Recognition Using Nonlinear Time Series Analysis and Random Forest. IEEE Access 2019, 7, 176448–176457. [Google Scholar] [CrossRef]
  54. Chen, X.; Chen, J.; Liang, J.; Li, Y.; Courtney, C.A.; Yang, Y. Entropy-Based Surface Electromyogram Feature Extraction for Knee Osteoarthritis Classification. IEEE Access 2019, 7, 164144–164151. [Google Scholar] [CrossRef]
  55. Cuong-Le, T.; Hoang-Le, M.; Samir, K.; Magd, A.W.; Minh, T.T.; Seyedali, M. A novel version of Cuckoo search algorithm for solving optimization problems. Expert Syst. Appl. 2021, 186, 115669. [Google Scholar] [CrossRef]
  56. Jamil, M.; Yang, X.S. A Literature Survey of Benchmark Functions for Global Optimization Problems. Mathematics 2013, 4, 150–194. [Google Scholar]
  57. Artemiadis, P.K.; Kyriakopoulos, K.J. EMG-Based Control of a Robot Arm Using Low-Dimensional Embeddings. IEEE Trans. Rob. 2010, 26, 393–398. [Google Scholar] [CrossRef]
Figure 1. The block diagram of knee joint extension force estimation.
Figure 1. The block diagram of knee joint extension force estimation.
Sensors 22 04651 g001
Figure 2. The procedure of ICS-SVR.
Figure 2. The procedure of ICS-SVR.
Sensors 22 04651 g002
Figure 3. The search space of the four test benchmark functions.
Figure 3. The search space of the four test benchmark functions.
Sensors 22 04651 g003
Figure 4. Comparison of convergence curves of the four algorithms for the four test benchmark functions.
Figure 4. Comparison of convergence curves of the four algorithms for the four test benchmark functions.
Sensors 22 04651 g004
Figure 5. GRA analysis of MMG features of subject S1.
Figure 5. GRA analysis of MMG features of subject S1.
Sensors 22 04651 g005
Figure 6. Results of knee joint extension force estimation for subject S1 with feature combination D.
Figure 6. Results of knee joint extension force estimation for subject S1 with feature combination D.
Sensors 22 04651 g006
Figure 7. Statistical analysis of knee joint extension force estimation with different regression models. (a) The R-value of the force estimation against the actual observed values. (b) The RMSE of the force estimation against the actual observed values. (c) The MAPE of the force estimation against the actual observed values.
Figure 7. Statistical analysis of knee joint extension force estimation with different regression models. (a) The R-value of the force estimation against the actual observed values. (b) The RMSE of the force estimation against the actual observed values. (c) The MAPE of the force estimation against the actual observed values.
Sensors 22 04651 g007
Table 1. Parameter values in the four optimization algorithms with a population size of 20.
Table 1. Parameter values in the four optimization algorithms with a population size of 20.
AlgorithmParameter
PSOc1 = c2 = 1, k = 0.5, wV = 0.9, wP = 0.9
GWOr1, r2 ∈ (0, 1), a ∈ (0, 2)
CS P a = 0.25 ,   β = 1.5 ,   a 0 = 0.01
ICS β = 1.5
PSO = particle swarm optimization algorithm; GWO = gray wolf optimization algorithm; CS = cuckoo search algorithm; ICS = improved cuckoo search algorithm.
Table 2. Results of optimization algorithms for solving four benchmark functions with D = 2.
Table 2. Results of optimization algorithms for solving four benchmark functions with D = 2.
Test FunctionAlgorithmOptimal
Solution
Worst
Solution
Average
Solution
SD
Rosenbrock functionPSO0.01423.96891.07741.0380
GWO4.377 × 10−72.503 × 10−42.383 × 10−54.341 × 10−5
CS1.389 × 10−133.605 × 10−47.822 × 10−65.093 × 10−5
ICS1.816 × 10−139.313 × 10−74.922 × 10−81.541 × 10−7
Griewank functionPSO0.383813.03721.71291.8178
GWO00.02710.00540.0056
CS9.678 × 10−70.00890.00440.0031
ICS2.618 × 10−60.00810.00250.0029
Cross-in-tray functionPSO−2.0622−1.8755−2.02600.0438
GWO−2.0626−2.0626−2.06261.257 × 10−7
CS−2.0626−2.0626−2.06266.519 × 10−11
ICS−2.0626−2.0626−2.06263.352 × 10−11
Schaffer functionPSO0.01130.24430.09400.0391
GWO00.08510.02040.0367
CS4.749 × 10−70.08510.00190.0120
ICS2.113 × 10−90.00856.824 × 10−40.0016
PSO = particle swarm optimization algorithm; GWO = gray wolf optimization algorithm; CS = cuckoo search algorithm; ICS = improved cuckoo search algorithm.
Table 3. Estimated results of different feature combinations of subject S1.
Table 3. Estimated results of different feature combinations of subject S1.
CombinationRMSEMAPER
Feature combination A0.17680.04050.9963
Feature combination B0.74930.03770.9954
Feature combination C0.16110.04160.9946
Feature combination D0.13190.03490.9966
Feature combination E0.42820.03270.9957
RMSE = root-mean-square error; MAPE = mean absolute percentage error; R = correlation coefficient.
Table 4. Results of the knee joint extension force estimation for the five subjects with feature sequences of different feature combinations.
Table 4. Results of the knee joint extension force estimation for the five subjects with feature sequences of different feature combinations.
CombinationRMSE ± SDMAPE ± SDR ± SD
Feature combination A0.7511 ± 0.76450.0550 ± 0.01200.9937 ± 0.0042
Feature combination B0.6012 ± 0.38400.0604 ± 0.01900.9912 ± 0.0046
Feature combination C0.6706 ± 0.52020.0649 ± 0.02410.9920 ± 0.0026
Feature combination D0.2761± 0.23960.0522± 0.02080.9949± 0.0016
Feature combination E0.8214 ± 0.67180.0644 ± 0.04640.9916 ± 0.0059
Values are the mean ± SD.
Table 5. Comparison of the estimation results for S1 using different models.
Table 5. Comparison of the estimation results for S1 using different models.
ModelRMSEMAPER
BPNN0.3952 ± 0.32460.0361 ± 0.00640.9954 ± 0.0012
ELM1.0464 ± 0.66730.1071 ± 0.01900.9681 ± 0.0130
CS-SVR0.1424 ± 0.02740.0358 ± 0.00260.9923 ± 0.0008
ICS-SVR0.1295±0.00210.0349±8.06 × 10−60.9966±1.38 × 10−6
BPNN = backpropagation neural network; ELM = extreme learning machine; CS-SVR = support vector regression optimized by cuckoo search algorithm; ICS-SVR = support vector regression optimized by the improved cuckoo search algorithm.
Table 6. Comparison of the estimation results for the five subjects using different models.
Table 6. Comparison of the estimation results for the five subjects using different models.
ModelRMSEMAPER
BPNN0.4706 ± 0.12990.0575 ± 0.02410.9919 ± 0.0022
ELM2.5252 ± 0.95070.2223 ± 0.08220.9137 ± 0.0341
CS-SVR0.2673 ± 0.10610.0603 ± 0.01930.9932 ± 0.0011
ICS-SVR0.2205±0.08400.0565±0.01850.9945±0.0013
Values are the mean ± SD.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Z.; Gao, L.; Lu, W.; Wang, D.; Cao, H.; Zhang, G. Estimation of Knee Extension Force Using Mechanomyography Signals Based on GRA and ICS-SVR. Sensors 2022, 22, 4651. https://doi.org/10.3390/s22124651

AMA Style

Li Z, Gao L, Lu W, Wang D, Cao H, Zhang G. Estimation of Knee Extension Force Using Mechanomyography Signals Based on GRA and ICS-SVR. Sensors. 2022; 22(12):4651. https://doi.org/10.3390/s22124651

Chicago/Turabian Style

Li, Zebin, Lifu Gao, Wei Lu, Daqing Wang, Huibin Cao, and Gang Zhang. 2022. "Estimation of Knee Extension Force Using Mechanomyography Signals Based on GRA and ICS-SVR" Sensors 22, no. 12: 4651. https://doi.org/10.3390/s22124651

APA Style

Li, Z., Gao, L., Lu, W., Wang, D., Cao, H., & Zhang, G. (2022). Estimation of Knee Extension Force Using Mechanomyography Signals Based on GRA and ICS-SVR. Sensors, 22(12), 4651. https://doi.org/10.3390/s22124651

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