CN114063787B - Deep learning processing analysis method based on EMG and IMU data - Google Patents

Deep learning processing analysis method based on EMG and IMU data Download PDF

Info

Publication number
CN114063787B
CN114063787B CN202111396749.9A CN202111396749A CN114063787B CN 114063787 B CN114063787 B CN 114063787B CN 202111396749 A CN202111396749 A CN 202111396749A CN 114063787 B CN114063787 B CN 114063787B
Authority
CN
China
Prior art keywords
data
emg
value
signals
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111396749.9A
Other languages
Chinese (zh)
Other versions
CN114063787A (en
Inventor
王鹤翔
许德新
张建强
张原野
冯帝
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN202111396749.9A priority Critical patent/CN114063787B/en
Publication of CN114063787A publication Critical patent/CN114063787A/en
Application granted granted Critical
Publication of CN114063787B publication Critical patent/CN114063787B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F3/00Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
    • G06F3/01Input arrangements or combined input and output arrangements for interaction between user and computer
    • G06F3/017Gesture based interaction, e.g. based on a set of recognized hand gestures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F3/00Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
    • G06F3/01Input arrangements or combined input and output arrangements for interaction between user and computer
    • G06F3/011Arrangements for interaction with the human body, e.g. for user immersion in virtual reality
    • G06F3/014Hand-worn input/output arrangements, e.g. data gloves
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02Preprocessing
    • G06F2218/04Denoising
    • G06F2218/06Denoising by applying a scale-space analysis, e.g. using wavelet analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12Classification; Matching

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Human Computer Interaction (AREA)
  • Image Analysis (AREA)

Abstract

The invention provides a deep learning processing analysis method based on EMG and IMU data, which comprises the following steps: obtaining EMG and IMU data signals corresponding to each gesture action; performing bias removal processing and wavelet noise reduction processing; filtering operation is carried out on ACC signals, and data segmentation operation is carried out on the processed signals; fusing the processed EMG signal and ACC signal data, sequentially storing gesture actions corresponding to each row of the data obtained after data fusion into a matrix in sequence, wherein the obtained matrix is in the shape of n rows and 1 column, and n is the number of rows of the data fusion array; normalizing the elements in the data fusion array; and selecting CNN as a recognition model of the data, and modifying the original CNN. After the operations are all completed, the data can be put into the model for training, and the data analyzed through the method is processed through experimental verification, so that the higher recognition rate is ensured on the premise of improving the recognition speed of the model.

Description

Deep learning processing analysis method based on EMG and IMU data
Technical Field
The invention belongs to the technical field of data processing, and particularly relates to a deep learning processing analysis method based on EMG and IMU data.
Background
In order to realize functions such as remote control, a machine system which is based on deep learning and can realize man-machine interaction through gesture recognition is designed, and the machine system has the characteristics of being novel, efficient, flexible, easy to use and the like. Can be applied to medical systems to assist the movement of patients, and is also beneficial to exploring areas where people are inconvenient to operate.
The core of gesture recognition is the selection and design of a gesture recognition model, and a plurality of existing model algorithms are used for recognition aiming at gesture recognition. For example, a network model algorithm such as LSTM, RNN, SVM can be selected for gesture recognition through electromyographic signals, and an image processing model such as CNN can be selected for classification through gesture recognition.
For the traditional text classification model, the processing time of the model such as RNN, LSTM and the like is long, which is unfavorable for the model to process real-time data. The recognition rate of the machine learning model SVM is lower than that of the deep neural network under the condition of processing a large amount of data, so that the problem of low recognition rate can be caused by selecting the SVM. Whereas CNN network models are generally applicable to image data, they have a lower recognition rate for recognizing text data than LSTM models.
Disclosure of Invention
Aiming at the technical problems, the invention provides a deep learning processing analysis method based on EMG and IMU data, so as to ensure higher recognition rate on the premise of improving the recognition speed of a model.
The image processing model CNN is selected to process text data, but if CNN data is directly selected, the problem of recognition rate is caused. Therefore, the method is improved on the traditional CNN model, and a new algorithm is introduced for improvement, so that the recognition rate of the model is improved under the condition of ensuring that the recognition speed is unchanged.
The specific technical scheme is as follows:
an EMG and IMU data-based deep learning processing analysis method comprises the following steps:
step 1: the wearable data acquisition device executes corresponding gesture actions to obtain EMG and IMU data signals corresponding to each gesture action;
step 2: performing unbiasing processing on the acquired EMG data signals, namely subtracting the average value of any 200 sampling points of the EMG signals when each section is in a resting state from the whole section of EMG signals;
step 3: the EMG signal after the offset removal treatment is subjected to wavelet noise reduction treatment, and the specific operation method comprises the following steps:
(3.1) performing discrete wavelet decomposition on the data; namely, selecting a proper wavelet base, and decomposing an EMG signal into 4 layers through discrete wavelet transformation; obtaining a low-frequency wavelet coefficient and a high-frequency wavelet coefficient from the data after discrete wavelet transformation processing; the discrete wavelet transform formula is as follows:
(3.2) selecting a threshold lambda; the threshold value can be determined by a fixed threshold formula:
wherein, media (|s (i) |) is the median of the wavelet coefficient absolute values;
(3.3) modifying the wavelet decomposition coefficients; if the absolute value of the wavelet decomposition coefficient is smaller than the threshold lambda, making the absolute value of the wavelet decomposition coefficient be 0; when the value is larger than lambda, the value is kept unchanged;
(3.4) performing discrete wavelet inverse transformation on the unmodified low-frequency wavelet coefficient and the modified high-frequency wavelet xishu1 to obtain a denoised signal;
step 4: filtering the ACC signal, namely replacing the maximum value and the minimum value of a certain piece of data with the average value of the maximum value and the minimum value of the certain piece of data;
step 5: performing data segmentation operation on the processed signals;
the specific operation method of the step is as follows:
(5.1) summing the three channel EMG signals to yield sum_emg:
wherein a (i) is the maximum value of the EMG signal of each channel;
(5.2) moving average processing is performed on the absolute value of the sum_emg signal to obtain S (k):
(5.3) establishing a comparison threshold th:
th=μ+1.6σ
in the method, in the process of the invention, μ is the mean value of SUM_EMG in the resting state; sigma is the standard deviation of SUM_EMG in a resting state;
(5.4) comparing S (k) with th, recording sequence numbers corresponding to sequences from the beginning of S (k) being greater than th to the beginning of S (k) being less than th, extracting the sequence, and successfully segmenting the EMG signals;
(5.5) dividing the ACC data according to the obtained EMG serial number to successfully divide the ACC signal;
step 6: fusing the processed EMG signal and ACC signal data, namely fusing the EMG signal and the ACC signal into a large array, wherein each row of the array represents a gesture action, and the first 60% of elements of each row are EMG signals and the second 30% are IMU signals;
step 7: sequentially storing gesture actions corresponding to each row of data obtained after data fusion into a matrix according to a sequence, wherein the obtained matrix is in the shape of n rows and 1 column, and n is the number of rows of the data fusion array;
step 8: normalizing the elements in the data fusion array to ensure that the absolute values of the elements are all numbers in the range of 0-1:
s in 1 For the element value after normalization treatment, s is the untreated original element value, s mean Is the mean value of s, s var Is the variance of s;
step 9: selecting CNN as a recognition model of data, and modifying the original CNN as follows:
(1) Modifying the size of the convolution kernel, setting the height of the convolution kernel to be 1, and adapting the width of the convolution kernel to the size of specific input data;
(2) The number of layers of the network is properly increased according to the size of the data, and the method is suitable for stacking small convolution kernels to replace large convolution kernels so as to reduce parameters required by model training;
(3) A residual block is introduced, i.e. the input of a certain layer and the outputs of the following two layers are added together as the input of the latest layer.
After the operations are all completed, the data can be put into the model for training, and the data analyzed through the method is processed through experimental verification, so that the higher recognition rate is ensured on the premise of improving the recognition speed of the model.
Drawings
FIG. 1 is a schematic flow chart of the present invention;
FIG. 2 is a schematic diagram of the large array of step 6;
fig. 3 is a flowchart illustrating the modification of the original CNN in step 9.
Detailed Description
The specific technical scheme of the invention is described by combining the embodiments.
As shown in fig. 1, a deep learning processing analysis method based on EMG and IMU data includes the following steps:
step 1: the wearable data acquisition device executes corresponding gesture actions to obtain EMG and IMU data signals corresponding to each gesture action;
step 2: performing unbiasing processing on the acquired EMG data signals, namely subtracting the average value of any 200 sampling points of the EMG signals when each section is in a resting state from the whole section of EMG signals;
step 3: the EMG signal after the offset removal treatment is subjected to wavelet noise reduction treatment, and the specific operation method comprises the following steps:
(3.1) performing discrete wavelet decomposition on the data; namely, selecting a proper wavelet base, and decomposing an EMG signal into 4 layers through discrete wavelet transformation; obtaining a low-frequency wavelet coefficient and a high-frequency wavelet coefficient from the data after discrete wavelet transformation processing; the discrete wavelet transform formula is as follows:
(3.2) selecting a threshold lambda; the threshold value can be determined by a fixed threshold formula:
wherein, media (|s (i) |) is the median of the wavelet coefficient absolute values;
(3.3) modifying the wavelet decomposition coefficients; if the absolute value of the wavelet decomposition coefficient is smaller than the threshold lambda, making the absolute value of the wavelet decomposition coefficient be 0; when the value is larger than lambda, the value is kept unchanged;
(3.4) performing discrete wavelet inverse transformation on the unmodified low-frequency wavelet coefficient and the modified high-frequency wavelet xishu1 to obtain a denoised signal;
step 4: filtering the ACC signal, namely replacing the maximum value and the minimum value of a certain piece of data with the average value of the maximum value and the minimum value of the certain piece of data;
step 5: performing data segmentation operation on the processed signals;
the specific operation method of the step is as follows:
(5.1) summing the three channel EMG signals to yield sum_emg:
wherein a (i) is the maximum value of the EMG signal of each channel;
(5.2) moving average processing is performed on the absolute value of the sum_emg signal to obtain S (k):
(5.3) establishing a comparison threshold th:
th=μ+1.6σ
in the method, in the process of the invention, μ is the mean value of SUM_EMG in the resting state; sigma is the standard deviation of SUM_EMG in a resting state;
(5.4) comparing S (k) with th, recording sequence numbers corresponding to sequences from the beginning of S (k) being greater than th to the beginning of S (k) being less than th, extracting the sequence, and successfully segmenting the EMG signals;
(5.5) dividing the ACC data according to the obtained EMG serial number to successfully divide the ACC signal;
step 6: fusing the processed EMG signal and ACC signal data, namely fusing the EMG signal and the ACC signal into a large array, wherein each row of the array represents a gesture action, and the first 60% of the elements of each row are EMG signals and the second 30% of the elements of each row are IMU signals as shown in FIG. 2;
step 7: sequentially storing gesture actions corresponding to each row of data obtained after data fusion into a matrix according to a sequence, wherein the obtained matrix is in the shape of n rows and 1 column, and n is the number of rows of the data fusion array;
step 8: normalizing the elements in the data fusion array to ensure that the absolute values of the elements are all numbers in the range of 0-1:
s in 1 For the element value after normalization treatment, s is the untreated original element value, s mean Is the mean value of s, s var Is the variance of s;
step 9: CNN was selected as the recognition model for the data, as shown in fig. 3, and the following modifications were made to the original CNN:
(1) Modifying the size of the convolution kernel, setting the height of the convolution kernel to be 1, and adapting the width of the convolution kernel to the size of specific input data;
(2) The number of layers of the network is properly increased according to the size of the data, and the method is suitable for stacking small convolution kernels to replace large convolution kernels so as to reduce parameters required by model training;
(3) A residual block is introduced, i.e. the input of a certain layer and the outputs of the following two layers are added together as the input of the latest layer.
In order to verify the actual effect of the model, 4 gesture actions of fist making, left swing, right swing and hand opening are respectively designed, and 4 healthy adults aged 19-21 are selected as experimental objects, wherein the experimental objects comprise two men and two women. In order to increase the collection quantity of data, each person wears the data collection device to execute specific gesture actions 150 times. Considering that repeated execution of a certain action causes muscle soreness and affects the recognition rate of gestures, the test subjects rest for 5 minutes every 10 actions to relieve the muscle soreness. Thus 600 sets of experimental data are available for each gesture. The obtained model data is put into the designed model after pretreatment such as noise reduction and segmentation. In order to verify that introducing a residual network does improve the overall performance of the model, these data are put into a common CNN network and a modified design network, respectively, with the following results:
it can be seen that the improved network is superior to the CNN network model, both in terms of identification rate and identification speed. And for the MYO arm ring, the identification rate and the stability of the invention are better.
The user firstly collects and executes data corresponding to different gestures through the myoelectric sensor and the IMU, the collected data firstly carries out data preprocessing and data segmentation operation, then carries out data processing operation, and after the processing is completed, EMG signal data and IMU signal data are fused together, so that the processing is completed. The improved model can be used for training the processed data by modifying the size of the CNN convolution kernel and introducing a residual block. Multiple experiments prove that the recognition rate of the model to the data is about 84%.
For data processing methods, there are many other functional formulas that can be implemented for noise reduction, filtering, feature segmentation, the key is the problem of selection and combination of these formulas, i.e. how to select the appropriate formula for each data processing step to optimize the overall process
For the selection aspect of the recognition model, other algorithms can be selected as classification recognition models, such as LSTM (long short term memory neural network), SVM (support vector machine) and KNN (K nearest neighbor classifier), but the recognition rate and the recognition speed are not as good as the models designed by adopting the project.

Claims (2)

1. The deep learning processing analysis method based on the EMG and IMU data is characterized by comprising the following steps of:
step 1: the wearable data acquisition device executes corresponding gesture actions to obtain EMG and IMU data signals corresponding to each gesture action;
step 2: performing unbiasing processing on the acquired EMG data signals, namely subtracting the average value of any 200 sampling points of the EMG signals when each section is in a resting state from the whole section of EMG signals;
step 3: performing wavelet noise reduction on the EMG signal subjected to the offset removal treatment;
step 4: filtering the ACC signal, namely replacing the maximum value and the minimum value of a certain piece of data with the average value of the maximum value and the minimum value of the certain piece of data;
step 5: performing data segmentation operation on the processed signals;
the specific operation method of the step 5 is as follows:
(5.1) summing the three channel EMG signals to yield sum_emg:
wherein a (i) is the maximum value of the EMG signal of each channel;
(5.2) moving average processing is performed on the absolute value of the sum_emg signal to obtain S (k):
(5.3) establishing a comparison threshold th:
th=μ+1.6σ
in the method, in the process of the invention, μ is the mean value of SUM_EMG in the resting state; sigma is the standard deviation of SUM_EMG in a resting state;
(5.4) comparing S (k) with th, recording sequence numbers corresponding to sequences from the beginning of S (k) being greater than th to the beginning of S (k) being less than th, extracting the sequence, and successfully segmenting the EMG signals;
(5.5) dividing the ACC data according to the obtained EMG serial number to successfully divide the ACC signal;
step 6: fusing the processed EMG signal and ACC signal data, namely fusing the EMG signal and the ACC signal into a large array, wherein each row of the array represents a gesture action, and the first 60% of elements of each row are EMG signals and the second 30% are IMU signals;
step 7: sequentially storing gesture actions corresponding to each row of data obtained after data fusion into a matrix according to a sequence, wherein the obtained matrix is in the shape of n rows and 1 column, and n is the number of rows of the data fusion array;
step 8: normalizing the elements in the data fusion array to ensure that the absolute values of the elements are all numbers in the range of 0-1:
s in 1 For the element value after normalization treatment, s is the untreated original element value, s mean Is the mean value of s, s var Is the variance of s;
step 9: selecting CNN as a recognition model of the data, and modifying the original CNN;
in step 9, the original CNN is modified as follows:
(1) Modifying the size of the convolution kernel, setting the height of the convolution kernel to be 1, and adapting the width of the convolution kernel to the size of specific input data;
(2) The number of layers of the network is properly increased according to the size of the data, and the method is suitable for stacking small convolution kernels to replace large convolution kernels so as to reduce parameters required by model training;
(3) A residual block is introduced, i.e. the input of a certain layer and the outputs of the following two layers are added together as the input of the latest layer.
2. The method for deep learning processing analysis based on EMG and IMU data according to claim 1, wherein the specific operation method in step 3 is as follows:
(3.1) performing discrete wavelet decomposition on the data; namely, selecting a proper wavelet base, and decomposing an EMG signal into 4 layers through discrete wavelet transformation; obtaining a low-frequency wavelet coefficient and a high-frequency wavelet coefficient from the data after discrete wavelet transformation processing; the discrete wavelet transform formula is as follows:
(3.2) selecting a threshold lambda; the threshold value can be determined by a fixed threshold formula:
wherein, media (|s (i) |) is the median of the wavelet coefficient absolute values;
(3.3) modifying the wavelet decomposition coefficients; if the absolute value of the wavelet decomposition coefficient is smaller than the threshold lambda, making the absolute value of the wavelet decomposition coefficient be 0; when the value is larger than lambda, the value is kept unchanged;
(3.4) performing inverse discrete wavelet transform on the unmodified low-frequency wavelet coefficient and the modified high-frequency wavelet xishu1 to obtain a denoised signal.
CN202111396749.9A 2021-11-23 2021-11-23 Deep learning processing analysis method based on EMG and IMU data Active CN114063787B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111396749.9A CN114063787B (en) 2021-11-23 2021-11-23 Deep learning processing analysis method based on EMG and IMU data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111396749.9A CN114063787B (en) 2021-11-23 2021-11-23 Deep learning processing analysis method based on EMG and IMU data

Publications (2)

Publication Number Publication Date
CN114063787A CN114063787A (en) 2022-02-18
CN114063787B true CN114063787B (en) 2023-09-19

Family

ID=80275597

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111396749.9A Active CN114063787B (en) 2021-11-23 2021-11-23 Deep learning processing analysis method based on EMG and IMU data

Country Status (1)

Country Link
CN (1) CN114063787B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20170030139A (en) * 2015-09-08 2017-03-17 경희대학교 산학협력단 System and method of controlling mobile robot using inertia measurement unit and electromyogram sensor-based gesture recognition
WO2020151075A1 (en) * 2019-01-23 2020-07-30 五邑大学 Cnn-lstm deep learning model-based driver fatigue identification method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20170030139A (en) * 2015-09-08 2017-03-17 경희대학교 산학협력단 System and method of controlling mobile robot using inertia measurement unit and electromyogram sensor-based gesture recognition
WO2020151075A1 (en) * 2019-01-23 2020-07-30 五邑大学 Cnn-lstm deep learning model-based driver fatigue identification method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于多传感器融合的动态手势识别研究分析;马正华;李雷;乔玉涛;戎海龙;曹海婷;;计算机工程与应用(第17期);全文 *

Also Published As

Publication number Publication date
CN114063787A (en) 2022-02-18

Similar Documents

Publication Publication Date Title
WO2021143353A1 (en) Gesture information processing method and apparatus, electronic device, and storage medium
Arvaneh et al. Optimizing spatial filters by minimizing within-class dissimilarities in electroencephalogram-based brain–computer interface
CN110555468A (en) Electroencephalogram signal identification method and system combining recursion graph and CNN
Zeng et al. Optimizing single-trial EEG classification by stationary matrix logistic regression in brain–computer interface
Kumar et al. OPTICAL+: a frequency-based deep learning scheme for recognizing brain wave signals
CN112656431A (en) Electroencephalogram-based attention recognition method and device, terminal equipment and storage medium
Duan et al. Recognition of motor imagery tasks for BCI using CSP and chaotic PSO twin SVM
CN114781441B (en) EEG motor imagery classification method and multi-space convolution neural network model
CN112450885B (en) Epileptic electroencephalogram-oriented identification method
CN111951958A (en) Pain data evaluation method based on self-coding and related components
CN111797674A (en) MI electroencephalogram signal identification method based on feature fusion and particle swarm optimization algorithm
Hsu et al. EEG-Channel-Temporal-Spectral-Attention correlation for motor imagery EEG classification
CN113143261B (en) Myoelectric signal-based identity recognition system, method and equipment
CN114305453A (en) Multi-source manifold electroencephalogram feature transfer learning method
CN114063787B (en) Deep learning processing analysis method based on EMG and IMU data
Ma et al. Iterative outlier removal clustering based time-frequency-spatial feature selection for binary EEG motor imagery decoding
CN113642528B (en) Hand movement intention classification method based on convolutional neural network
CN116522106A (en) Motor imagery electroencephalogram signal classification method based on transfer learning parallel multi-scale filter bank time domain convolution
Zhang et al. Low‐Rank Linear Dynamical Systems for Motor Imagery EEG
Gabardi et al. A multi-artifact EEG denoising by frequency-based deep learning
CN112818315B (en) Electrocardiosignal identity recognition method and system fusing multi-feature information
Li et al. A semi-supervised SVM learning algorithm for joint feature extraction and classification in brain computer interfaces
CN114190884B (en) Longitudinal analysis method, system and device for brain disease data
CN114587384A (en) Motor imagery electroencephalogram signal feature extraction method combining low-rank representation and manifold learning
CN113159205A (en) Sparse time-frequency block common space mode feature extraction method based on optimal channel

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant