CN105326500A - 基于表面肌电信号的动作识别方法和设备 - Google Patents

基于表面肌电信号的动作识别方法和设备 Download PDF

Info

Publication number
CN105326500A
CN105326500A CN201410396422.5A CN201410396422A CN105326500A CN 105326500 A CN105326500 A CN 105326500A CN 201410396422 A CN201410396422 A CN 201410396422A CN 105326500 A CN105326500 A CN 105326500A
Authority
CN
China
Prior art keywords
signal
surface electromyogram
electromyogram signal
active surface
electromyographic 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.)
Granted
Application number
CN201410396422.5A
Other languages
English (en)
Other versions
CN105326500B (zh
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.)
Tsinghua University
Huawei Technologies Co Ltd
Original Assignee
Tsinghua University
Huawei Technologies Co Ltd
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 Tsinghua University, Huawei Technologies Co Ltd filed Critical Tsinghua University
Priority to CN201410396422.5A priority Critical patent/CN105326500B/zh
Priority to PCT/CN2015/073369 priority patent/WO2016023349A1/zh
Priority to EP15832528.2A priority patent/EP3167799A4/en
Publication of CN105326500A publication Critical patent/CN105326500A/zh
Priority to US15/428,230 priority patent/US20170143226A1/en
Application granted granted Critical
Publication of CN105326500B publication Critical patent/CN105326500B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/389Electromyography [EMG]
    • A61B5/397Analysis of electromyograms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/389Electromyography [EMG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7275Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7278Artificial waveform generation or derivation, e.g. synthesising signals from measured signals

Abstract

本发明实施例提供一种基于表面肌电信号的动作识别方法和设备,所述方法通过获取多个通道的表面肌电信号,根据多个通道的表面肌电信号确定有效表面肌电信号;然后,确定有效表面肌电信号的频率,根据有效表面肌电信号的频率确定多个通道的表面肌电信号对应的肢体动作。由于表面肌电信号的频率与信号强度等特征无关,因此,本实施例的方法能显著提高基于表面肌电信号的动作识别的准确率。而且,以频率作为识别特征,用户不需要进行大幅度的动作,给用户带来更好的体验。另外,本发明实施例的方法,还可以将表面肌电信号的频率和幅值特征结合起来识别对应的肢体动作,不仅可以提高表面肌电信号识别准确率,而且能够增加肢体动作的识别种类。

Description

基于表面肌电信号的动作识别方法和设备
技术领域
本发明实施例涉及动作识别技术,尤其涉及一种基于表面肌电信号的动作识别方法和设备。
背景技术
表面肌电信号(Electromyography,简称EMG)是一种与神经肌肉活动相关的生物电信号,表面肌电信号能够反应肌肉的收缩模式以及收缩强度等信息,不同的肢体动作对应不同的表面肌电信号,通过分析表面肌电信号可以判别出该表面肌电信号对应的具体动作,因此,表面肌电信号被广泛应用于临床医学、运动医学、生物医学与康复工程等诸多领域。
现有技术中,以表面肌电信号的幅值特征作为识别参数进行肢体动作的识别,该方法通过将采集到的表面肌电信号取滑动加窗,得到每个滑动窗口的窗口序列,然后,计算每个窗口的窗口序列的幅度得到幅值特征,将该幅值特征与预先训练得到的各肢体动作对应的表面肌电信号幅值特征进行比较以确定该幅值特征对应的肢体动作。
现有技术的方法存在以下问题:由于用户的其他动作也会产生表面肌电信号,当用户进行其他运动(跑步、打字等)时,产生的干扰表面肌电信号与目标表面肌电信号会发生叠加,导致识别准确率低,另外,表面肌电信号的幅值的大小正比于用户肢体动作的幅度,为了达到较高的信噪比,用户需要进行较大的力度的肢体动作,在长时间操作的情况下,会使用户感到吃力。
发明内容
本发明实施例提供一种基于表面肌电信号的动作识别方法和设备,能够提高基于表面肌电信号的动作识别的准确率。
本发明第一方面提供一种基于表面肌电信号的动作识别方法,包括:
获取多个通道的表面肌电信号;
根据所述多个通道的表面肌电信号确定有效表面肌电信号;
确定所述有效表面肌电信号的频率;
根据所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作。
结合本发明第一方面,在本发明第一方面的第一种可能的实现方式中,所述根据所述多个通道的表面肌电信号确定有效表面肌电信号,包括:
将所述多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号;
从所述单通道的表面肌电信号的起始时间开始,在每个滑动时刻对所述单通道的表面肌电信号进行滑动得到所述每个滑动时刻对应的窗口,确定所述每个滑动时刻对应的窗口序列,计算所述每个滑动时刻对应的窗口序列的表面肌电信号平均幅值,其中,所述每个滑动时刻之间相差一个滑动间隔,所述窗口序列的表面肌电信号平均幅值为所述窗口序列内表面肌电信号的幅值的绝对值的平均值,所述滑动时刻对应的窗口序列包括所述滑动时刻对应的窗口以及所述滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数;
若所述每个滑动时刻中的滑动时刻T对应的窗口序列的表面肌电信号平均幅值不小于预设幅值,将所述滑动时刻T对应的窗口序列的起始时间作为所述有效表面肌电信号的起始时间,将所述有效表面肌电信号的起始时间加上预设时间得到所述有效表面肌电信号的截止时间,截取所述有效表面肌电信号的开始时间和截止时间之间的所述多个通道的表面肌电信号作为所述有效表面肌电信号。
结合本发明第一方面的第一种可能的实现方式,在本发明第一方面的第二种可能的实现方式中,所述预设幅值为所述多个通道的表面肌电信号叠加后的表面肌电信号的幅值的绝对值的平均值。
结合本发明第一方面以及本发明第一方面的第一种和第二种可能的实现方式,在本发明第一方面的第三种可能的实现方式中,所述确定所述有效表面肌电信号的频率,包括:
分别计算所述有效表面肌电信号与多个正余弦矩阵的相关系数,其中,所述正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同;
判断所述有效表面肌电信号与所述多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数;
若所述最大相关系数大于所述预设的相关系数,则将所述最大相关系数对应的正余弦矩阵的基频作为所述有效表面肌电信号的频率。
本发明第二方面提供一种基于表面肌电信号的动作识别方法,包括:
获取多个通道的表面肌电信号;
根据所述多个通道的表面肌电信号确定有效表面肌电信号;
确定所述有效表面肌电信号的频率;
提取所述有效表面肌电信号的幅值特征;
根据所述有效表面肌电信号的幅值特征和所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作。
结合本发明第二方面,在本发明第二方面的第一种可能的实现方式中,所述根据所述多个通道的表面肌电信号确定有效表面肌电信号,包括:
将所述多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号;
从所述单通道的表面肌电信号的起始时间开始,在每个滑动时刻对所述单通道的表面肌电信号进行滑动得到所述每个滑动时刻对应的窗口,确定所述每个滑动时刻对应的窗口序列,计算所述每个滑动时刻对应的窗口序列的表面肌电信号平均幅值,其中,所述每个滑动时刻之间相差一个滑动间隔,所述窗口序列的表面肌电信号平均幅值为所述窗口序列内表面肌电信号的幅值的绝对值的平均值,所述滑动时刻对应的窗口序列包括所述滑动时刻对应的窗口以及所述滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数;
当所述每个滑动时刻中的滑动时刻T对应的窗口序列的表面肌电信号平均幅值不小于预设幅值时,将所述滑动时刻T对应的窗口序列的起始时间作为所述有效表面肌电信号的起始时间,将所述有效表面肌电信号的起始时间加上预设时间得到所述有效表面肌电信号的截止时间,截取所述有效表面肌电信号的开始时间和截止时间之间的所述多个通道的表面肌电信号作为所述有效表面肌电信号。
结合本发明第二方面的第一种可能的实现方式,在本发明第二方面的第二种可能的实现方式中,所述预设幅值为所述多个通道的表面肌电信号叠加后的表面肌电信号的幅值的绝对值的平均值。
结合本发明第二方面以及本发明第二方面的第一种和第二种可能的实现方式,在本发明第二方面的第三种可能的实现方式中,所述确定所述有效表面肌电信号的频率,包括:
分别计算所述有效表面肌电信号与多个正余弦矩阵的相关系数,其中,所述正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同;
判断所述有效表面肌电信号与所述多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数;
若所述最大相关系数大于所述预设的相关系数,则将所述最大相关系数对应的正余弦矩阵的基频作为所述有效表面肌电信号的频率。
结合本发明第二方面的第三种可能的实现方式,在本发明第二方面的第四种可能的实现方式中,所述提取所述有效表面肌电信号的幅值特征,包括:
对所述有效表面肌电信号的每个通道的表面肌电信号分别进行滑动窗口处理;
计算所述有效表面肌电信号的每个通道的表面肌电信号的每个滑动窗口的平均幅值,其中,所述每个滑动窗口的平均幅值为所述每个滑动窗口内表面肌电信号的幅值的绝对值的平均值,将所述有效表面肌电信号的每个滑动窗口的平均幅值作为所述有效表面肌电信号的幅值特征。
结合本发明第二方面以及本发明第二方面的第一种和第二种可能的实现方式,在本发明第二方面的第五种可能的实现方式中,所述根据所述有效表面肌电信号的幅值特征和所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作,包括:
根据所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的多个备选肢体动作;
将所述有效表面肌电信号的幅值特征与预先训练得到所述多个备选肢体动作的幅值特征进行匹配,得到与所述有效表面肌电信号的幅值特征匹配的肢体动作,将与所述有效表面肌电信号的幅值特征匹配的肢体动作作为所述多个通道的表面肌电信号对应的肢体动作。
本发明第三方面提供一种基于表面肌电信号的动作识别设备,包括:
获取模块,用于获取多个通道的表面肌电信号;
第一确定模块,用于根据所述多个通道的表面肌电信号确定有效表面肌电信号;
第二确定模块,用于确定所述有效表面肌电信号的频率;
识别模块,用于根据所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作。
结合本发明第三方面,在本发明第三方面的第一种可能的实现方式中,所述第一确定模块具体用于:
将所述多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号;
从所述单通道的表面肌电信号的起始时间开始,在每个滑动时刻对所述单通道的表面肌电信号进行滑动得到所述每个滑动时刻对应的窗口,确定所述每个滑动时刻对应的窗口序列,计算所述每个滑动时刻对应的窗口序列的表面肌电信号平均幅值,其中,所述每个滑动时刻之间相差一个滑动间隔,所述窗口序列的表面肌电信号平均幅值为所述窗口序列内表面肌电信号的幅值的绝对值的平均值,所述滑动时刻对应的窗口序列包括所述滑动时刻对应的窗口以及所述滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数;
若所述每个滑动时刻中的滑动时刻T对应的窗口序列的表面肌电信号平均幅值不小于预设幅值,将所述滑动时刻T对应的窗口序列的起始时间作为所述有效表面肌电信号的起始时间,将所述有效表面肌电信号的起始时间加上预设时间得到所述有效表面肌电信号的截止时间,截取所述有效表面肌电信号的开始时间和截止时间之间的所述多个通道的表面肌电信号作为所述有效表面肌电信号。
结合本发明第三方面的第一种可能的实现方式,在本发明第三方面的第二种可能的实现方式中,所述预设幅值为所述多个通道的表面肌电信号叠加后的表面肌电信号的幅值的绝对值的平均值。
结合本发明第三方面以及第三方面的第一种和第二种可能的实现方式,在本发明第三方面的第三种可能的实现方式中,所述第二确定模块具体用于:
分别计算所述有效表面肌电信号与多个正余弦矩阵的相关系数,其中,所述正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同;
判断所述有效表面肌电信号与所述多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数;
若所述最大相关系数大于所述预设的相关系数,则将所述最大相关系数对应的正余弦矩阵的基频作为所述有效表面肌电信号的频率。
本发明第四方面提供一种基于表面肌电信号的动作识别设备,包括:
获取模块,用于获取多个通道的表面肌电信号;
第一确定模块,用于根据所述多个通道的表面肌电信号确定有效表面肌电信号;
第二确定模块,用于确定所述有效表面肌电信号的频率;
提取模块,用于提取所述有效表面肌电信号的幅值特征;
识别模块,用于根据所述有效表面肌电信号的幅值特征和所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作。
结合本发明第四方面,在本发明第四方面的第一种可能的实现方式中,所述第一确定模块具体用于:
将所述多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号;
从所述单通道的表面肌电信号的起始时间开始,在每个滑动时刻对所述单通道的表面肌电信号进行滑动得到所述每个滑动时刻对应的窗口,确定所述每个滑动时刻对应的窗口序列,计算所述每个滑动时刻对应的窗口序列的表面肌电信号平均幅值,其中,所述每个滑动时刻之间相差一个滑动间隔,所述窗口序列的表面肌电信号平均幅值为所述窗口序列内表面肌电信号的幅值的绝对值的平均值,所述滑动时刻对应的窗口序列包括所述滑动时刻对应的窗口以及所述滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数;
当所述每个滑动时刻中的滑动时刻T对应的窗口序列的表面肌电信号平均幅值不小于预设幅值时,将所述滑动时刻T对应的窗口序列的起始时间作为所述有效表面肌电信号的起始时间,将所述有效表面肌电信号的起始时间加上预设时间得到所述有效表面肌电信号的截止时间,截取所述有效表面肌电信号的开始时间和截止时间之间的所述多个通道的表面肌电信号作为所述有效表面肌电信号。
结合本发明第四方面的第一种可能的实现方式,在本发明第四方面的第二种可能的实现方式中,所述预设幅值为所述多个通道的表面肌电信号叠加后的表面肌电信号的幅值的绝对值的平均值。
结合本发明第四方面以及第四方面的第一种和第二种可能的实现方式,在本发明第四方面的第三种可能的实现方式中,所述第二确定模块具体用于:
分别计算所述有效表面肌电信号与多个正余弦矩阵的相关系数,其中,所述正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同;
判断所述有效表面肌电信号与所述多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数;
若所述最大相关系数大于所述预设的相关系数,则将所述最大相关系数对应的正余弦矩阵的基频作为所述有效表面肌电信号的频率。
结合本发明第四方面的第三种可能的实现方式,在本发明第四方面的第四种可能的实现方式中,所所述提取模块具体用于:
对所述有效表面肌电信号的每个通道的表面肌电信号分别进行滑动窗口处理;
计算所述有效表面肌电信号的每个通道的表面肌电信号的每个滑动窗口的平均幅值,其中,所述每个滑动窗口的平均幅值为所述每个滑动窗口内表面肌电信号的幅值的绝对值的平均值,将所述有效表面肌电信号的每个滑动窗口的平均幅值作为所述有效表面肌电信号的幅值特征。
结合本发明第四方面以及第四方面的第一种和第二种可能的实现方式,在本发明第四方面的第五种可能的实现方式中,所述识别模块具体用于:
根据所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的多个备选肢体动作;
将所述有效表面肌电信号的幅值特征与预先训练得到所述多个备选肢体动作的幅值特征进行匹配,得到与所述有效表面肌电信号的幅值特征匹配的肢体动作,将与所述有效表面肌电信号的幅值特征匹配的肢体动作作为所述多个通道的表面肌电信号对应的肢体动作。
本发明第五方面提供一种基于表面肌电信号的动作识别设备,包括:
处理器、存储器和系统总线,所述处理器和所述存储器之间通过所述系统总线连接并完成相互间的通信;
所述存储器,用于存储计算机执行指令;
所述处理器,用于运行所述计算机执行指令,使所述动作识别设备执行本发明第一方面以及本发明第一方面的第一种至第三种可能的实现方式提供的方法。
本发明第六方面提供一种基于表面肌电信号的动作识别设备,包括:
处理器、存储器和系统总线,所述处理器和所述存储器之间通过所述系统总线连接并完成相互间的通信;
所述存储器,用于存储计算机执行指令;
所述处理器,用于运行所述计算机执行指令,使所述动作识别设备执行本发明第二方面以及本发明第二方面的第一种至第五种可能的实现方式提供的方法。
本发明实施例提供的基于表面肌电信号的动作识别方法和设备,通过获取多个通道的表面肌电信号,根据多个通道的表面肌电信号确定有效表面肌电信号;然后,确定有效表面肌电信号的频率,根据有效表面肌电信号的频率确定多个通道的表面肌电信号对应的肢体动作。由于表面肌电信号的频率与信号强度等特征无关,因此,本实施例的方法能显著提高基于表面肌电信号的动作识别的准确率。而且以频率作为识别特征,用户不需要进行大幅度的动作,给用户带来更好的体验。另外,本发明实施例的方法,还可以将表面肌电信号的频率和幅值特征结合起来识别表面肌电信号对应的肢体动作,不仅可以提高表面肌电信号识别准确率,而且能够增加肢体动作的识别种类。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一提供的基于表面肌电信号的动作识别方法的流程图;
图2为本发明实施例二提供的一种有效表面肌电信号的频率确定方法;
图3为本发明实施例三提供的基于表面肌电信号的动作识别方法的流程图;
图4为本发明实施例四提供的基于表面肌电信号的动作识别方法的整体框图;
图5为本发明实施例五提供的基于表面肌电信号的动作识别设备的结构示意图;
图6为本发明实施例六提供的基于表面肌电信号的动作识别设备的结构示意图;
图7为本发明实施例七提供的基于表面肌电信号的动作识别设备的结构示意图;
图8为本发明实施例八提供的基于表面肌电信号的动作识别设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明实施例一提供的基于表面肌电信号的动作识别方法的流程图,本实施例的方法可以由基于表面肌电信号的动作识别设备执行,该基于表面肌电信号的动作识别设备具体可以为智能手机、平板电脑等终端设备。如图1所示,本实施例的方法可以包括以下步骤:
步骤101、获取多个通道的表面肌电信号。
基于表面肌电信号的动作识别设备获取多个通道的表面肌电信号具体为:接收采集装置发送的多个通道的表面肌电信号,通过将采集装置放置在要采集的肌肉群表面进行表面肌电信号的采集,由于本实施例的方法中,是根据表面肌电信号的频率确定对应的肢体动作,因此,在采集表面肌电信号时要求用户以特定的频率重复进行同一种肢体动作。该肢体动作不仅可以是人体上肢前臂的节律性手势,还可以是下肢腿部的节律性动作,甚至可以是脖子、腹部等躯体的节律性动作,这里节律性动作即按照一定频率重复的动作。
采集装置包括多个传感节点,每个传感节点采集的数据作为一个通道的表面肌电信号,采集装置可以嵌入智能手表、手环等穿戴式设备中。在用户进行某一肢体动作时,可能会带动多个肌肉群运动,因此,需要对多个肌肉群都进行表面肌电信号采集,采集到的多个通道的信号综合反映用户的肢体动作。
步骤102、根据多个通道的表面肌电信号确定有效表面肌电信号。
由于表面肌电信号的幅度小、信噪比低,因此,在确定有效表面肌电信号之前需要对获取的多个通道的表面肌电信号进行预处理,对表面肌电信号进行预处理具体为:对表面肌电信号进行信号放大、工频滤波、高通滤波等处理。
采集装置采集到的表面肌电信号的最开始部分的信号可能有一些干扰信号,例如,采集装置工作后用户没有及时的做相应的肢体动作,此时用户虽然没有做肢体动作,采集装置仍然能够采集到一些微弱的表面肌电信号,这些微弱的表面肌电信号属于干扰信号,因此,需要祛除可能的干扰信号得到有效的表面肌电信号。
具体地,可以根据如下方法确定有效表面肌电信号:
第一步、将预处理后多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号。
第二步、从单通道的表面肌电信号的起始时间开始,在每个滑动时刻对单通道表面肌电信号进行滑动得到每个滑动时刻对应的窗口,确定每个滑动时刻对应的窗口序列,计算每个滑动时刻对应的窗口序列的表面肌电信号平均幅值,其中,每个滑动时刻之间相差一个滑动间隔,窗口序列的表面肌电信号平均幅值为窗口序列内表面肌电信号的幅值的绝对值的平均值,其中,滑动时刻对应的窗口序列包括滑动时刻对应的窗口以及滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数。
其中,单通道表面肌电信号的起始时间是指开始有通道信号的时间,在每个滑动时刻对单通道的表面肌电信号进行滑动,即在每个滑动时刻将单通道表面信号在时间上向后移动一个窗口,相邻两个滑动时刻之间相差的时间为一个滑动间隔。例如,窗口宽度为时间t,滑动间隔为t/2,t=1ms(毫秒),每个窗口序列包括四个窗口,共有50个窗口,窗口的编号为1、2、3……50,那么从单通道表面肌电信号的起始时间开始,第一个滑动时刻对应的窗口为4号窗口,第一个滑动时刻对应的窗口序列为0-4ms,共4个窗口,然后,在第二个滑动时刻将单通道表面肌电信号向后滑动0.5ms,得到第二个窗口序列,第二个窗口序列为0.5-4.5ms,依次类推,得到每个滑动时刻对应的窗口序列。在计算第一个滑动时刻对应的窗口序列的表面肌电信号平均幅值时,首先,对该窗口序列内的单通道的表面肌电信号的值取绝对值,然后,将该窗口序列内的单通道的表面肌电信号的绝对值加起来取平均值,得到该窗口序列的表面肌电信号平均幅值。
第三步、若每个滑动时刻中的滑动时刻T对应的窗口序列的表面肌电信号平均幅值不小于预设幅值,将滑动时刻T对应的窗口序列的起始时间作为有效表面肌电信号的起始时间,将有效表面肌电信号的起始时间加上预设时间得到有效表面肌电信号的截止时间,截取有效表面肌电信号的开始时间和截止时间之间的多个通道预处理后的表面肌电信号作为有效表面肌电信号。
在计算得到每个滑动时刻对应的窗口序列的表面肌电信号平均幅值时,判断该滑动时刻对应的窗口序列的表面肌电信号平均幅值是否小于预设幅值,若滑动时刻T对应的窗口序列的表面肌电信号平均幅值不小于预设幅值,即滑动时刻T对应的窗口序列的表面肌电信号平均幅值大于等于预设幅值,则将滑动时刻T对应的窗口序列的起始时间作为有效表面肌电信号的起始时间,有效表面肌电信号的起始时间加上预设时间得到有效表面肌电信号的截止时间。在得到有效表面肌电信号的开始时间和截止时间后,截取开始时间和截止时间之间多个通道的表面肌电信号作为有效表面肌电信号。其中,预设幅值可以为多个通道的表面肌电信号叠加后的表面肌电信号的幅值的绝对值的平均值。
步骤103、确定有效表面肌电信号的频率。
确定有效表面肌电信号的频率的方法有多种,以下列举几种常用的方法:第一种方式,计算有效表面肌电信号与多个正余弦矩阵的相关系数,假设肢体动作可能出现的频率为f1,…,fn,那么选取以f1,…,fn为基频的正余弦矩阵,分别将有效表面肌电信号与各正余弦矩阵进行典型相关分析(CanonicalCorrelationAnalysis,简称CCA)运算,得到有效表面肌电信号与各正余弦矩阵的最大相关系数,如果最大相关系数大于预设的相关系数,则将最大相关系数对应的正余弦矩阵的基频作为有效表面肌电信号的频率。第二种方法,对有效表面信号进行快速傅里叶变换(FastFourierTtransform,简称FFT),根据变换结果确定有效表面肌电信号的频率分布,从而得到有效表面肌电信号的频率。第三种方法,根据有效表面肌电信号两个过零点的时间间隔计算表面肌电信号的频率。
步骤104、根据有效表面肌电信号的频率确定多个通道的表面肌电信号对应的肢体动作。
预先需要建立表面肌电信号的频率与肢体动作的对应关系,例如,1HZ频率对应的肢体动作为“握拳”,2HZ频率对应的肢体动作为“胳膊肘弯曲”,3HZ频率对应的肢体动作为“ok手势”,那么当用户以1HZ的频率做握拳动作时,根据获取到的有效表面肌电信号的频率确定用户的肢体动作为握拳。
本发明实施例,通过获取多个通道的表面肌电信号,根据多个通道的表面肌电信号确定有效表面肌电信号;然后,确定有效表面肌电信号的频率,最后,根据表面肌电信号的频率确定多个通道的表面肌电信号对应的肢体动作。由于表面肌电信号的频率与信号强度等特征无关,因此,本实施例的方法能显著提高基于表面肌电信号的动作识别的准确率。而且以频率作为识别特征,用户不需要进行大幅度的动作,给用户带来更好的体验。
另外,本实施例的方法,具有很好的抗干扰性能,能有效抵抗无关动作所产生的肌电噪音干扰,用户可在跑步、开车、做家务等非静止条件下使用。并且不受使用过程中皮肤湿度变化、电极接触情况变化、肌肉疲劳程度等因素的影响,稳定性高。用户使用方便,无需事先采集用户训练数据,每次使用前也无需重新训练。
本发明实施例二中将对实施例一中步骤103进行详细说明,图2为本发明实施例二提供的一种有效表面肌电信号的频率确定方法,如图2所示,本实施例的方法包括以下步骤:
步骤201、分别计算有效表面肌电信号与多个正余弦矩阵的相关系数,其中,正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同。
假设肢体动作可能出现的频率为f1,…,fn,那么选取以f1,…,fn为基频的正余弦矩阵为对照。分别将有效表面肌电信号与这些正余弦矩阵进行CCA运算,获取有效表面肌电信号与各正余弦矩阵的相关系数。
记n通道的有效表面肌电信号为x=(x1,x2,x3,...,xn)T,x1,x2,x3,...,xn分别代表每个通道的信号,正余弦矩阵为y=[cos(2πf1t),sin(2πf1t),cos(4πf1t),sin(4πf1t)]T,f1是基频。则CCA可以定义为如下问题:分别寻找向量wx和wy,使得x和y在向量wx和wx上的投影X=xTwx和Y=yTwy之间的相关值最大。也就是使得下式中的ρ最大:
ρ = E [ XY ] E [ X 2 ] E [ Y 2 ] = E [ w x T xy T w y ] E [ w x T xx T w x ] E [ w y T yy T w y ] = w x T E [ xy T ] w y w x T E [ xx T ] w x w y T E [ yy T ] w y = w x T C xy w y w x T C xx w x w y T C yy w y
其中Cxy表示x,y的互相关矩阵,Cxx表示x的自相关矩阵,Cyy表示y的自相关矩阵。寻找向量wx和wy的方法如下:
第一步、通过拉格朗日算法,构建如下拉格朗日算子:
L = w x T C xy w y - λx 2 ( w x T C xx w x - 1 ) - λy 2 ( w y T C yy w y - 1 ) .
第二步、将L分别对wx和wy求偏导,得到:
CxywyxCxxwx=0
CyxwyyCyywy=0,
从而得到公式1和公式2:
w y = C yy - 1 C yx w x λ - - - ( 1 )
λ=λx=λy,Cxx -1CxyCyy -1Cyxwx=λ2wx(2)
根据公式2可以求出wx,根据公式1可以求出wy,通过公式2可以看出根据公式2求wx的过程被转换为特征值分解问题,将wx和wy代入到ρ的定义中得到相关系数。
步骤202、判断有效表面肌电信号与多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数。
在计算出有效表面肌电信号与每个正余弦矩阵的相关系数后,找出最大相关系数,判断该最大相关系数是否大于预设的相关系数。
步骤203、若最大相关系数大于预设的相关系数,则将最大相关系数对应的正余弦矩阵的基频作为有效表面肌电信号的频率。
若最大相关系数大于预设的相关系数,说明有效表面肌电信号具备最大相关系数对应的正余弦矩阵的频率,将最大相关系数对应的正余弦矩阵的基频作为有效表面肌电信号的频率。
图3为本发明实施例三提供的基于表面肌电信号的动作识别方法的流程图,本实施例的方法与实施例一的区别在于:本实施例中,将结合表面肌电信号的频率和幅值特征确定表面肌电信号对应的肢体动作。
步骤301、获取多个通道的表面肌电信号。
步骤302、根据多个通道的表面肌电信号确定有效表面肌电信号。
根据多个通道的表面肌电信号确定有效表面肌电信号,具体为:将多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号;从单通道的表面肌电信号的起始时间开始,在每个滑动时刻对单通道的表面肌电信号进行滑动得到每个滑动时刻对应的窗口,确定每个滑动时刻对应的窗口序列,计算每个滑动时刻对应的窗口序列的肌电信号平均幅值,其中,每个滑动时刻之间相差一个滑动间隔,窗口序列的表面肌电信号平均幅值为窗口序列内表面肌电信号的幅值的绝对值的平均值,滑动时刻对应的窗口序列包括滑动时刻对应的窗口以及滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数。当每个滑动时刻中的滑动时刻T对应的窗口序列的表面肌电信号平均幅值不小于预设幅值时,将滑动时刻T对应的窗口序列的起始时间作为有效表面肌电信号的起始时间,将有效表面肌电信号的起始时间加上预设时间得到有效表面肌电信号的截止时间,截取开始时间和所述截止时间之间的多个通道的表面肌电信号作为有效表面肌电信号。其中,预设幅值可以为多个通道的表面肌电信号叠加后的信号的幅值的绝对值的平均值。
步骤303、确定有效表面肌电信号的频率。
确定有效表面肌电信号的频率,具体为:分别计算有效表面肌电信号与多个正余弦矩阵的相关系数,其中,正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同;判断有效表面肌电信号与多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数;若最大相关系数大于预设的相关系数,则将最大相关系数对应的正余弦矩阵的基频作为有效表面肌电信号的频率。
步骤301-303的具体实现方式可以参照实施例一和实施例二的相关描述,这里不再赘述。
步骤304、提取有效表面肌电信号的幅值特征。
提取有效表面肌电信号的幅值特征,具体可以为:首先,对有效表面肌电信号的每个通道的表面肌电信号分别进行滑动窗口处理,得到多个滑动窗口,例如滑动窗口的宽度为100ms,滑动间隔为100ms。然后,计算有效表面肌电信号的每个通道的表面肌电信号的每个滑动窗口的平均幅值,其中,每个滑动窗口的平均幅值为每个滑动窗口内表面肌电信号的幅值的绝对值的平均值,将有效表面肌电信号的每个滑动窗口的平均幅值作为有效表面肌电信号的幅值特征。
需说明的时,本实施例中,步骤303和步骤304在执行时并没有先后顺序,也可以先执行步骤304,再执行步骤303。
步骤305、根据有效表面肌电信号的幅值特征和有效表面肌电信号的频率确定多个通道的表面肌电信号对应的肢体动作。
具体地,首先,根据有效表面肌电信号的频率确定多个通道的表面肌电信号对应的多个备选肢体动作,本实施例中,每个频率可以对应多个肢体动作,例如1HZ的频率可以分别对应以下三种肢体动作:“握拳”、“OK手势”和“胳膊肘弯曲”。那么,当确定表面肌电信号的频率为1HZ时,表面肌电信号对应上述三个备选肢体动作。然后,将有效表面肌电信号的幅值特征与预先训练得到多个备选肢体动作的幅值特征进行匹配,得到与有效表面肌电信号的幅值特征匹配的肢体动作,将与有效表面肌电信号的幅值特征匹配的肢体动作作为多个通道的表面肌电信号对应的肢体动作。可以将备选肢体动作的幅值特征作为模板,采用线性判别式分析(LinearDiscriminantAnalysis,简称LDA)分类器识别表面肌电信号对应的肢体动作。
本实施例的方法,通过获取表面肌电信号的频率和幅值特征,根据表面肌电信号的频率和幅值特征识别表面肌电信号对应的肢体动作,不仅可以提高基于表面肌电信号的动作识别的准确率,而且能够增加表面肌电信号识别种类,因为同一种肢体动作以不同频率做节律性运动可以视为不同的肢体动作,因此本实施例在识别种类数目上有明显的提升。
将表面肌电信号的频率和幅值特征结合确定表面肌电信号对应的肢体动作时,需要预先对表面肌电信号进行训练获取各种肢体动作的幅值特征。如图4所示,图4为本发明实施例四提供的表面肌电信号的识别方法的整体框图,本实施例的识别方法将表面肌电信号的识别过程分为两个部分:基于幅值特征识别方法和基于频率的识别方法。
在基于幅值特征的识别方法中,先要对表面肌电信号进行训练,得到训练模板。表面肌电信号的训练过程具体包括以下步骤:
第一步,采集各种肢体动作的多个通道的表面肌电信号。
让用户将不同肢体动作(例如“握拳”和“OK”)重复做多遍,通过采集装置采集肢体动作关联的主要肌肉群的表面肌电信号,采集的表面肌电信号由多个通道组成,并记录每遍采集时间的时间。
第二步,对每种肢体动作的多个通道的表面肌电信号进行预处理。
例如,对采集到的多通道的表面肌电信号进行50HZ的工频干扰陷波,并使用FIR滤波器进行高通滤波,得到预处理后的表面肌电信号。
第三步,提取每种肢体动作的多个通道的表面肌电信号的幅值特征。
具体的提取方法可以参照实施例三种步骤304的描述,这里不再赘述。
举例来说,用户将握拳动作重复做30遍,每遍采集开始的时间为[t1,t2,…,t30],采集到的表面肌电信号由8个通道的信号组成。在提取表面肌电信号的幅值特征时,依次在(t1,t1+300ms),(t2,t2+300ms),…(t30,t30+300ms)共30个时间段内,将预处理后的表面肌电信号的各通道的信号按照100ms的滑动窗间隔进行滑动,滑动窗口的宽度也为100ms,计算这30个时间段内各窗口的信号的绝对值的平均值M(n),将M(n)作为握拳肢体动作的幅值特征,按照同样的方法可以得到所有肢体动作的幅值特征。
第四步、训练模板的制作,即建立每种肢体动作与幅值特征的对应关系,然后,制作好的训练模板发送给分类器,由分类器根据训练模板进行肢体动作的识别。
上述四个步骤为训练阶段,在识别阶段,同样要提取表面肌电信号的幅值特征,然后,将幅值特征输入分类器进行识别。
在基于频率的识别方法中,不需要对表面肌电信号进行训练,在识别阶段,第一步,采集多个通道的表面肌电信号;第二步,对多个通道的表面肌电信号进行预处理;第三步,确定有效表面肌电信号的频率;第四步,进行CCA计算,即计算有效表面肌电信号与多个正余弦矩阵的相关系数,得到有效表面肌边信号的最大相关系数;第五步,确定有效表面肌电信号的频率。上述第一步至第五步的具体实现方式可操作实施例一和实施例二描述,这里不再赘述。在确定有效表面肌电信号的频率后,根据有效表面肌电信号的频率确定备选肢体动作。
本实施例中,在确定有效表面肌电信号后,提取有效表面肌电信号的幅值特征,将有效表面肌电信号的幅值特征输入分类器,分类器将备选肢体动作的幅值特征作为模板,根据输入的有效表面肌电信号的幅值特征和备选肢体动作的幅值特征进行特征匹配,得到有效表面肌电信号对应的肢体动作。
本发明各实施例的方法具有广泛的应用场景:(1)作为穿戴设备的启动命令:例如通过多次、按特定频率转动手腕来启动肌电手环,避免误操作;(2)在开车、跑步等非静止状态下使用:例如调节音乐音量、切换歌曲、接听电话等;(3)由于识别的准确度高,因此可以应用于触发控制指令。比如:残疾人通过不同节律的手势来控制轮椅的方向和速度;(4)空中鼠标,例如将手指不同频率的转动赋予相应的鼠标的操作;(5)简单的游戏控制:手指的快,慢2个运动可以分别对应赛车游戏的加速和刹车;(6)用于康复治疗过程中对用户的身体协调和控制能力进行测试。
将所述有效表面肌电信号的幅值特征与预先训练得到所述多个备选肢体动作的幅值特征进行匹配,得到与所述有效表面肌电信号的幅值特征匹配的肢体动作,将与所述有效表面肌电信号的幅值特征匹配的肢体动作作为所述多个通道的表面肌电信号对应的肢体动作。
图5为本发明实施例五提供的基于表面肌电信号的动作识别设备的结构示意图,如图5所示,本实施例的基于表面肌电信号的动作识别设备包括:获取模块11、第一确定模块12、第二确定模块13和识别模块14。
其中,获取模块11,用于获取多个通道的表面肌电信号;
第一确定模块12,用于根据所述多个通道的表面肌电信号确定有效表面肌电信号;
第二确定模块13,用于确定所述有效表面肌电信号的频率;
识别模块14,用于根据所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作。
可选地,第一确定模块12具体用于:
将所述多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号;
从所述单通道的表面肌电信号的起始时间开始,在每个滑动时刻对所述单通道的表面肌电信号进行滑动得到所述每个滑动时刻对应的窗口,确定所述每个滑动时刻对应的窗口序列,计算所述每个滑动时刻对应的窗口序列的表面肌电信号平均幅值,其中,所述每个滑动时刻之间相差一个滑动间隔,所述窗口序列的表面肌电信号平均幅值为所述窗口序列内表面肌电信号的幅值的绝对值的平均值,所述滑动时刻对应的窗口序列包括所述滑动时刻对应的窗口以及所述滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数;
若所述每个滑动时刻中的滑动时刻T对应的窗口序列的表面肌电信号平均幅值不小于预设幅值,将所述滑动时刻T对应的窗口序列的起始时间作为所述有效表面肌电信号的起始时间,将所述有效表面肌电信号的起始时间加上预设时间得到所述有效表面肌电信号的截止时间,截取所述有效表面肌电信号的开始时间和截止时间之间的所述多个通道的表面肌电信号作为所述有效表面肌电信号。其中,预设幅值可以为多个通道的表面肌电信号叠加后的表面肌电信号的幅值的绝对值的平均值。
可选地,第二确定模块13具体用于:首先,分别计算所述有效表面肌电信号与多个正余弦矩阵的相关系数,其中,所述正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同;然后,判断所述有效表面肌电信号与所述多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数;若所述最大相关系数大于所述预设的相关系数,则将所述最大相关系数对应的正余弦矩阵的基频作为所述有效表面肌电信号的频率。
本发明实施例的设备可用于执行方法实施例一和实施例二的方案,具体实现方式和技术效果类似,这里不再赘述。
图6为本发明实施例六提供的基于表面肌电信号的动作识别设备的结构示意图,如图6所示,本实施例提供的基于表面肌电信号的动作识别设备包括:获取模块21、第一确定模块22、第二确定模块23、提取模块24和识别模块25。
其中,获取模块21,用于获取多个通道的表面肌电信号;
第一确定模块22,用于根据所述多个通道的表面肌电信号确定有效表面肌电信号;
第二确定模块23,用于确定所述有效表面肌电信号的频率;
提取模块24,用于提取所述有效表面肌电信号的幅值特征;
识别模块25,用于根据所述有效表面肌电信号的幅值特征和所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作。
可选地,第一确定模块22具体用于:
将所述多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号;
从所述单通道的表面肌电信号的起始时间开始,在每个滑动时刻对所述单通道的表面肌电信号进行滑动得到所述每个滑动时刻对应的窗口,确定所述每个滑动时刻对应的窗口序列,计算所述每个滑动时刻对应的窗口序列的表面肌电信号平均幅值,其中,所述每个滑动时刻之间相差一个滑动间隔,所述窗口序列的表面肌电信号平均幅值为所述窗口序列内表面肌电信号的幅值的绝对值的平均值,所述滑动时刻对应的窗口序列包括所述滑动时刻对应的窗口以及所述滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数;
当所述每个滑动时刻中的滑动时刻T对应的窗口序列的肌电信号平均幅值不小于预设幅值时,将所述滑动时刻T对应的窗口序列的起始时间作为所述有效表面肌电信号的起始时间,将所述有效表面肌电信号的起始时间加上预设时间得到所述有效表面肌电信号的截止时间,截取所述有效表面肌电信号的开始时间和截止时间之间的所述多个通道的表面肌电信号作为所述有效表面肌电信号。其中,预设幅值可以为多个通道的表面肌电信号叠加后的表面肌电信号的幅值的绝对值的平均值。
可选地,第二确定模块23具体用于:首先,分别计算所述有效表面肌电信号与多个正余弦矩阵的相关系数,其中,所述正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同;然后,判断所述有效表面肌电信号与所述多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数;若所述最大相关系数大于所述预设的相关系数,则将所述最大相关系数对应的正余弦矩阵的基频作为所述有效表面肌电信号的频率。
可选地,提取模块24具体用于:对所述有效表面肌电信号的每个通道的表面肌电信号分别进行滑动窗口处理;计算所述有效表面肌电信号的每个通道的表面肌电信号的每个滑动窗口的平均幅值,其中,所述每个滑动窗口的平均幅值为所述每个滑动窗口内表面肌电信号的幅值的绝对值的平均值,将所述有效表面肌电信号的每个滑动窗口的平均幅值作为所述有效表面肌电信号的幅值特征。
可选地,识别模块25具体用于:根据所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的多个备选肢体动作;将所述有效表面肌电信号的幅值特征与预先训练得到所述多个备选肢体动作的幅值特征进行匹配,得到与所述有效表面肌电信号的幅值特征匹配的肢体动作,将与所述有效表面肌电信号的幅值特征匹配的肢体动作作为所述多个通道的表面肌电信号对应的肢体动作。
本实施例的设备可用于执行实施例三提供的技术方案,具体实现方式和技术效果类似,这里不再赘述。
图7为本发明实施例七提供的基于表面肌电信号的动作识别设备的结构示意图,如图7所示,本实施例的基于表面肌电信号的动作识别设备300,包括:处理器31、存储器32和系统总线33,处理器31和存储器32之间通过所述系统总线33连接并完成相互间的通信;存储器32,用于存储计算机执行指令321;处理器31,用于运行计算机执行指令321,执行如下所述的方法:
获取多个通道的表面肌电信号;
根据所述多个通道的表面肌电信号确定有效表面肌电信号;
确定所述有效表面肌电信号的频率;
根据所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作。
可选地,处理器31在根据所述多个通道的表面肌电信号确定有效表面肌电信号时,具体用于:
将所述多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号;
从所述单通道的表面肌电信号的起始时间开始,在每个滑动时刻对所述单通道的表面肌电信号进行滑动得到所述每个滑动时刻对应的窗口,确定所述每个滑动时刻对应的窗口序列,计算所述每个滑动时刻对应的窗口序列的表面肌电信号平均幅值,其中,所述每个滑动时刻之间相差一个滑动间隔,所述窗口序列的表面肌电信号平均幅值为所述窗口序列内表面肌电信号的幅值的绝对值的平均值,所述滑动时刻对应的窗口序列包括所述滑动时刻对应的窗口以及所述滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数;
若所述每个滑动时刻中的滑动时刻T对应的窗口序列的表面肌电信号平均幅值不小于预设幅值,将所述滑动时刻T对应的窗口序列的起始时间作为所述有效表面肌电信号的起始时间,将所述有效表面肌电信号的起始时间加上预设时间得到所述有效表面肌电信号的截止时间,截取所述有效表面肌电信号的开始时间和截止时间之间的所述多个通道的表面肌电信号作为所述有效表面肌电信号。其中,预设幅值可以为所述多个通道的表面肌电信号叠加后的表面肌电信号的幅值的绝对值的平均值。
可选地,处理器31在确定所述有效表面肌电信号的频率时,具体用于:分别计算所述有效表面肌电信号与多个正余弦矩阵的相关系数,其中,所述正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同;判断所述有效表面肌电信号与所述多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数;若所述最大相关系数大于所述预设的相关系数,则将所述最大相关系数对应的正余弦矩阵的基频作为所述有效表面肌电信号的频率。
本实施例的设备可用于执行实施例一和实施例二的技术方案,具体实现方式和技术效果类似,这里不再赘述。
图8为本发明实施例八提供的基于表面肌电信号的动作识别设备的结构示意图,如图8所示,本实施例的基于表面肌电信号的动作识别设备400包括:处理器41、存储器42和系统总线43,所述处理器41和所述存储器42之间通过所述系统总线43连接并完成相互间的通信;所述存储器42,用于存储计算机执行指令421;所述处理器41,用于运行所述计算机执行指令421执行如下所述的方法:
获取多个通道的表面肌电信号;
根据所述多个通道的表面肌电信号确定有效表面肌电信号;
确定所述有效表面肌电信号的频率;
提取所述有效表面肌电信号的幅值特征;
根据所述有效表面肌电信号的幅值特征和所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作。
可选地,处理器41在根据所述多个通道的表面肌电信号确定有效表面肌电信号时,具体用于:
将所述多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号;
从所述单通道的表面肌电信号的起始时间开始,在每个滑动时刻对所述单通道的表面肌电信号进行滑动得到所述每个滑动时刻对应的窗口,确定所述每个滑动时刻对应的窗口序列,计算所述每个滑动时刻对应的窗口序列的表面肌电信号平均幅值,其中,所述每个滑动时刻之间相差一个滑动间隔,所述窗口序列的表面肌电信号平均幅值为所述窗口序列内表面肌电信号的幅值的绝对值的平均值,所述滑动时刻对应的窗口序列包括所述滑动时刻对应的窗口以及所述滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数;
当所述每个滑动时刻中的滑动时刻T对应的窗口序列的肌电信号平均幅值不小于预设幅值时,将所述滑动时刻T对应的窗口序列的起始时间作为所述有效表面肌电信号的起始时间,将所述有效表面肌电信号的起始时间加上预设时间得到所述有效表面肌电信号的截止时间,截取所述有效表面肌电信号的开始时间和截止时间之间的所述多个通道的表面肌电信号作为所述有效表面肌电信号。其中,所述预设幅值可以为所述多个通道的表面肌电信号叠加后的表面肌电信号的幅值的绝对值的平均值。
可选地,处理器41在确定所述有效表面肌电信号的频率时,具体用于:分别计算所述有效表面肌电信号与多个正余弦矩阵的相关系数,其中,所述正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同;判断所述有效表面肌电信号与所述多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数;若所述最大相关系数大于所述预设的相关系数,则将所述最大相关系数对应的正余弦矩阵的基频作为所述有效表面肌电信号的频率。
可选地,处理器41在提取所述有效表面肌电信号的幅值特征时,具体用于:对所述有效表面肌电信号的每个通道的表面肌电信号分别进行滑动窗口处理;计算所述有效表面肌电信号的每个通道的表面肌电信号的每个滑动窗口的平均幅值,其中,所述每个滑动窗口的平均幅值为所述每个滑动窗口内表面肌电信号的幅值的绝对值的平均值,将所述有效表面肌电信号的每个滑动窗口的平均幅值作为所述有效表面肌电信号的幅值特征。
可选地,处理器41在根据所述有效表面肌电信号的幅值特征和所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作时,具体用于:根据所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的多个备选肢体动作;将所述有效表面肌电信号的幅值特征与预先训练得到所述多个备选肢体动作的幅值特征进行匹配,得到与所述有效表面肌电信号的幅值特征匹配的肢体动作,将与所述有效表面肌电信号的幅值特征匹配的肢体动作作为所述多个通道的表面肌电信号对应的肢体动作。
本实施例的设备可用于执行实施例三的技术方案,具体实现方式和技术效果类似,这里不再赘述。
本领域普通技术人员可以理解:实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时,执行包括上述各方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (22)

1.一种基于表面肌电信号的动作识别方法,其特征在于,包括:
获取多个通道的表面肌电信号;
根据所述多个通道的表面肌电信号确定有效表面肌电信号;
确定所述有效表面肌电信号的频率;
根据所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作。
2.根据权利要求1所述的方法,其特征在于,所述根据所述多个通道的表面肌电信号确定有效表面肌电信号,包括:
将所述多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号;
从所述单通道的表面肌电信号的起始时间开始,在每个滑动时刻对所述单通道的表面肌电信号进行滑动得到所述每个滑动时刻对应的窗口,确定所述每个滑动时刻对应的窗口序列,计算所述每个滑动时刻对应的窗口序列的表面肌电信号平均幅值,其中,所述每个滑动时刻之间相差一个滑动间隔,所述窗口序列的表面肌电信号平均幅值为所述窗口序列内表面肌电信号的幅值的绝对值的平均值,所述滑动时刻对应的窗口序列包括所述滑动时刻对应的窗口以及所述滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数;
若所述每个滑动时刻中的滑动时刻T对应的窗口序列的表面肌电信号平均幅值不小于预设幅值,将所述滑动时刻T对应的窗口序列的起始时间作为所述有效表面肌电信号的起始时间,将所述有效表面肌电信号的起始时间加上预设时间得到所述有效表面肌电信号的截止时间,截取所述有效表面肌电信号的开始时间和截止时间之间的所述多个通道的表面肌电信号作为所述有效表面肌电信号。
3.根据权利要求2所述的方法,其特征在于,所述预设幅值为所述多个通道的表面肌电信号叠加后的表面肌电信号的幅值的绝对值的平均值。
4.根据权利要求1-3中任一项所述的方法,其特征在于,所述确定所述有效表面肌电信号的频率,包括:
分别计算所述有效表面肌电信号与多个正余弦矩阵的相关系数,其中,所述正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同;
判断所述有效表面肌电信号与所述多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数;
若所述最大相关系数大于所述预设的相关系数,则将所述最大相关系数对应的正余弦矩阵的基频作为所述有效表面肌电信号的频率。
5.一种基于表面肌电信号的动作识别方法,其特征在于,包括:
获取多个通道的表面肌电信号;
根据所述多个通道的表面肌电信号确定有效表面肌电信号;
确定所述有效表面肌电信号的频率;
提取所述有效表面肌电信号的幅值特征;
根据所述有效表面肌电信号的幅值特征和所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作。
6.根据权利要求5所述的方法,其特征在于,所述根据所述多个通道的表面肌电信号确定有效表面肌电信号,包括:
将所述多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号;
从所述单通道的表面肌电信号的起始时间开始,在每个滑动时刻对所述单通道的表面肌电信号进行滑动得到所述每个滑动时刻对应的窗口,确定所述每个滑动时刻对应的窗口序列,计算所述每个滑动时刻对应的窗口序列的表面肌电信号平均幅值,其中,所述每个滑动时刻之间相差一个滑动间隔,所述窗口序列的表面肌电信号平均幅值为所述窗口序列内表面肌电信号的幅值的绝对值的平均值,所述滑动时刻对应的窗口序列包括所述滑动时刻对应的窗口以及所述滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数;
当所述每个滑动时刻中的滑动时刻T对应的窗口序列的表面肌电信号平均幅值不小于预设幅值时,将所述滑动时刻T对应的窗口序列的起始时间作为所述有效表面肌电信号的起始时间,将所述有效表面肌电信号的起始时间加上预设时间得到所述有效表面肌电信号的截止时间,截取所述有效表面肌电信号的开始时间和截止时间之间的所述多个通道的表面肌电信号作为所述有效表面肌电信号。
7.根据权利要求6所述的方法,其特征在于,所述预设幅值为所述多个通道的表面肌电信号叠加后的表面肌电信号的幅值的绝对值的平均值。
8.根据权利要求5-7任一项所述的方法,其特征在于,所述确定所述有效表面肌电信号的频率,包括:
分别计算所述有效表面肌电信号与多个正余弦矩阵的相关系数,其中,所述正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同;
判断所述有效表面肌电信号与所述多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数;
若所述最大相关系数大于所述预设的相关系数,则将所述最大相关系数对应的正余弦矩阵的基频作为所述有效表面肌电信号的频率。
9.根据权利要求8所述的方法,其特征在于,所述提取所述有效表面肌电信号的幅值特征,包括:
对所述有效表面肌电信号的每个通道的表面肌电信号分别进行滑动窗口处理;
计算所述有效表面肌电信号的每个通道的表面肌电信号的每个滑动窗口的平均幅值,其中,所述每个滑动窗口的平均幅值为所述每个滑动窗口内表面肌电信号的幅值的绝对值的平均值,将所述有效表面肌电信号的每个滑动窗口的平均幅值作为所述有效表面肌电信号的幅值特征。
10.根据权利要求5-7中任一项所述的方法,其特征在于,所述根据所述有效表面肌电信号的幅值特征和所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作,包括:
根据所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的多个备选肢体动作;
将所述有效表面肌电信号的幅值特征与预先训练得到所述多个备选肢体动作的幅值特征进行匹配,得到与所述有效表面肌电信号的幅值特征匹配的肢体动作,将与所述有效表面肌电信号的幅值特征匹配的肢体动作作为所述多个通道的表面肌电信号对应的肢体动作。
11.一种基于表面肌电信号的动作识别设备,其特征在于,包括:
获取模块,用于获取多个通道的表面肌电信号;
第一确定模块,用于根据所述多个通道的表面肌电信号确定有效表面肌电信号;
第二确定模块,用于确定所述有效表面肌电信号的频率;
识别模块,用于根据所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作。
12.根据权利要求11所述的设备,其特征在于,所述第一确定模块具体用于:
将所述多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号;
从所述单通道的表面肌电信号的起始时间开始,在每个滑动时刻对所述单通道的表面肌电信号进行滑动得到所述每个滑动时刻对应的窗口,确定所述每个滑动时刻对应的窗口序列,计算所述每个滑动时刻对应的窗口序列的表面肌电信号平均幅值,其中,所述每个滑动时刻之间相差一个滑动间隔,所述窗口序列的表面肌电信号平均幅值为所述窗口序列内表面肌电信号的幅值的绝对值的平均值,所述滑动时刻对应的窗口序列包括所述滑动时刻对应的窗口以及所述滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数;
若所述每个滑动时刻中的滑动时刻T对应的窗口序列的表面肌电信号平均幅值不小于预设幅值,将所述滑动时刻T对应的窗口序列的起始时间作为所述有效表面肌电信号的起始时间,将所述有效表面肌电信号的起始时间加上预设时间得到所述有效表面肌电信号的截止时间,截取所述有效表面肌电信号的开始时间和截止时间之间的所述多个通道的表面肌电信号作为所述有效表面肌电信号。
13.根据权利要求12所述的设备,其特征在于,所述预设幅值为所述多个通道的表面肌电信号叠加后的表面肌电信号的幅值的绝对值的平均值。
14.根据权利要求11-13中任一项所述的设备,其特征在于,所述第二确定模块具体用于:
分别计算所述有效表面肌电信号与多个正余弦矩阵的相关系数,其中,所述正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同;
判断所述有效表面肌电信号与所述多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数;
若所述最大相关系数大于所述预设的相关系数,则将所述最大相关系数对应的正余弦矩阵的基频作为所述有效表面肌电信号的频率。
15.一种基于表面肌电信号的动作识别设备,其特征在于,包括:
获取模块,用于获取多个通道的表面肌电信号;
第一确定模块,用于根据所述多个通道的表面肌电信号确定有效表面肌电信号;
第二确定模块,用于确定所述有效表面肌电信号的频率;
提取模块,用于提取所述有效表面肌电信号的幅值特征;
识别模块,用于根据所述有效表面肌电信号的幅值特征和所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的肢体动作。
16.根据权利要求15所述的设备,其特征在于,所述第一确定模块具体用于:
将所述多个通道的表面肌电信号叠加在一起,对叠加后的表面肌电信号除以通道数量得到单通道的表面肌电信号;
从所述单通道的表面肌电信号的起始时间开始,在每个滑动时刻对所述单通道的表面肌电信号进行滑动得到所述每个滑动时刻对应的窗口,确定所述每个滑动时刻对应的窗口序列,计算所述每个滑动时刻对应的窗口序列的表面肌电信号平均幅值,其中,所述每个滑动时刻之间相差一个滑动间隔,所述窗口序列的表面肌电信号平均幅值为所述窗口序列内表面肌电信号的幅值的绝对值的平均值,所述滑动时刻对应的窗口序列包括所述滑动时刻对应的窗口以及所述滑动时刻之前的N-1个滑动时刻对应的N-1个窗口共N个连续的窗口,N为大于等于2的正整数;
当所述每个滑动时刻中的滑动时刻T对应的窗口序列的表面肌电信号平均幅值不小于预设幅值时,将所述滑动时刻T对应的窗口序列的起始时间作为所述有效表面肌电信号的起始时间,将所述有效表面肌电信号的起始时间加上预设时间得到所述有效表面肌电信号的截止时间,截取所述有效表面肌电信号的开始时间和截止时间之间的所述多个通道的表面肌电信号作为所述有效表面肌电信号。
17.根据权利要求16所述的设备,其特征在于,所述预设幅值为所述多个通道的表面肌电信号叠加后的表面肌电信号的幅值的绝对值的平均值。
18.根据权利要求15-17中任一项所述的设备,其特征在于,所述第二确定模块具体用于:
分别计算所述有效表面肌电信号与多个正余弦矩阵的相关系数,其中,所述正余弦矩阵由基频和倍频的正弦函数与余弦函数组成,每个正余弦矩阵的基频不同;
判断所述有效表面肌电信号与所述多个正余弦矩阵的相关系数中的最大相关系数是否大于预设的相关系数;
若所述最大相关系数大于所述预设的相关系数,则将所述最大相关系数对应的正余弦矩阵的基频作为所述有效表面肌电信号的频率。
19.根据权利要求18所述的设备,其特征在于,所述提取模块具体用于:
对所述有效表面肌电信号的每个通道的表面肌电信号分别进行滑动窗口处理;
计算所述有效表面肌电信号的每个通道的表面肌电信号的每个滑动窗口的平均幅值,其中,所述每个滑动窗口的平均幅值为所述每个滑动窗口内表面肌电信号的幅值的绝对值的平均值,将所述有效表面肌电信号的每个滑动窗口的平均幅值作为所述有效表面肌电信号的幅值特征。
20.根据权利要求15-17中任一项所述的设备,其特征在于,所述识别模块具体用于:
根据所述有效表面肌电信号的频率确定所述多个通道的表面肌电信号对应的多个备选肢体动作;
将所述有效表面肌电信号的幅值特征与预先训练得到所述多个备选肢体动作的幅值特征进行匹配,得到与所述有效表面肌电信号的幅值特征匹配的肢体动作,将与所述有效表面肌电信号的幅值特征匹配的肢体动作作为所述多个通道的表面肌电信号对应的肢体动作。
21.一种基于表面肌电信号的动作识别设备,其特征在于,包括:
处理器、存储器和系统总线,所述处理器和所述存储器之间通过所述系统总线连接并完成相互间的通信;
所述存储器,用于存储计算机执行指令;
所述处理器,用于运行所述计算机执行指令,使所述动作识别设备执行如权利要求1至4任一所述的方法。
22.一种基于表面肌电信号的动作识别设备,其特征在于,包括:
处理器、存储器和系统总线,所述处理器和所述存储器之间通过所述系统总线连接并完成相互间的通信;
所述存储器,用于存储计算机执行指令;
所述处理器,用于运行所述计算机执行指令,使所述动作识别设备执行如权利要求5至10任一所述的方法。
CN201410396422.5A 2014-08-13 2014-08-13 基于表面肌电信号的动作识别方法和设备 Expired - Fee Related CN105326500B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201410396422.5A CN105326500B (zh) 2014-08-13 2014-08-13 基于表面肌电信号的动作识别方法和设备
PCT/CN2015/073369 WO2016023349A1 (zh) 2014-08-13 2015-02-27 基于表面肌电信号的动作识别方法和设备
EP15832528.2A EP3167799A4 (en) 2014-08-13 2015-02-27 Movement recognition method and apparatus based on surface electromyographic signal
US15/428,230 US20170143226A1 (en) 2014-08-13 2017-02-09 Action Recognition Method and Device Based on Surface Electromyography Signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410396422.5A CN105326500B (zh) 2014-08-13 2014-08-13 基于表面肌电信号的动作识别方法和设备

Publications (2)

Publication Number Publication Date
CN105326500A true CN105326500A (zh) 2016-02-17
CN105326500B CN105326500B (zh) 2018-02-09

Family

ID=55277427

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410396422.5A Expired - Fee Related CN105326500B (zh) 2014-08-13 2014-08-13 基于表面肌电信号的动作识别方法和设备

Country Status (4)

Country Link
US (1) US20170143226A1 (zh)
EP (1) EP3167799A4 (zh)
CN (1) CN105326500B (zh)
WO (1) WO2016023349A1 (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106108897A (zh) * 2016-07-20 2016-11-16 西安中科比奇创新科技有限责任公司 一种基于经验模式分解的肌电信号滤波方法
CN107126303A (zh) * 2017-02-15 2017-09-05 上海术理智能科技有限公司 一种基于手机app的上下肢运动辅助方法
CN108108016A (zh) * 2017-12-07 2018-06-01 浙江大学 手势感知器
CN108268844A (zh) * 2018-01-17 2018-07-10 上海术理智能科技有限公司 基于表面肌电信号的动作模式识别方法及装置
CN108403108A (zh) * 2018-02-11 2018-08-17 宁波工程学院 基于波形优化的阵列式表面肌电信号分解方法
CN109062401A (zh) * 2018-07-11 2018-12-21 北京理工大学 一种基于肌电信号的实时手势识别系统
CN109222968A (zh) * 2017-07-10 2019-01-18 丰田自动车株式会社 康复评估设备、康复评估方法以及康复评估程序
CN109512424A (zh) * 2018-11-16 2019-03-26 福州大学 一种高密度或多通道肌电信号的肌肉激活起点检测方法
CN109948465A (zh) * 2019-02-28 2019-06-28 浙江工业大学 一种基于多目标的表面肌电信号分类方法
CN110399846A (zh) * 2019-07-03 2019-11-01 北京航空航天大学 一种基于多通道肌电信号相关性的手势识别方法
CN113986017A (zh) * 2021-12-27 2022-01-28 深圳市心流科技有限公司 一种肌电手势模板的生成方法、装置及存储介质
CN115299950A (zh) * 2022-08-12 2022-11-08 歌尔股份有限公司 肌电信号采集电路、可穿戴设备及控制方法

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109200469B (zh) * 2018-10-24 2020-05-05 山东大学 一种增强人工耳蜗音调感知的优化编码方法及系统
CN109805923A (zh) * 2019-01-29 2019-05-28 北京京东方光电科技有限公司 可穿戴设备、信号处理方法及装置
CN113625882B (zh) * 2021-10-12 2022-06-14 四川大学 一种基于稀疏多通道相关性特征的肌电手势识别方法
CN115034273B (zh) * 2021-12-27 2023-09-01 驻马店市中心医院 一种基于模式识别的肌电生物反馈设备及系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090326406A1 (en) * 2008-06-26 2009-12-31 Microsoft Corporation Wearable electromyography-based controllers for human-computer interface
CN101874758A (zh) * 2010-05-07 2010-11-03 中国科学院深圳先进技术研究院 辅助运动系统及其控制方法
CN101930285A (zh) * 2009-11-23 2010-12-29 上海交通大学 基于表面肌电信号的手写识别方法
WO2011020446A1 (en) * 2009-08-21 2011-02-24 The Chinese University Of Hong Kong Systems and methods for reproducing body motions via networks
CN102139139A (zh) * 2011-01-13 2011-08-03 中国医学科学院生物医学工程研究所 一种肌电反馈控制的电刺激装置及其控制方法
CN102426651A (zh) * 2011-08-25 2012-04-25 武汉理工大学 人体前臂表面肌电信号采集及模式识别系统
CN102961203A (zh) * 2012-12-10 2013-03-13 杭州电子科技大学 基于emd样本熵的表面肌电信号识别方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE0303061D0 (sv) * 2003-11-19 2003-11-19 Maquet Critical Care Ab Metod och anordning för att bestämma en EMG-signal
CN101587546B (zh) * 2009-06-12 2011-05-11 浙江大学 基于支持向量多元分类的下肢动作模式的识别方法
CN101995947B (zh) * 2009-08-21 2012-07-04 香港中文大学 通过网络感知用户动作的系统和方法
CN103495260B (zh) * 2013-06-04 2016-05-25 东南大学 一种基于肌电信号和传感器信号实现精细实时运动的控制方法
CN103654774B (zh) * 2014-01-02 2016-08-17 北京思睿博创科技有限公司 穿戴式移动手环

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090326406A1 (en) * 2008-06-26 2009-12-31 Microsoft Corporation Wearable electromyography-based controllers for human-computer interface
WO2011020446A1 (en) * 2009-08-21 2011-02-24 The Chinese University Of Hong Kong Systems and methods for reproducing body motions via networks
CN101930285A (zh) * 2009-11-23 2010-12-29 上海交通大学 基于表面肌电信号的手写识别方法
CN101874758A (zh) * 2010-05-07 2010-11-03 中国科学院深圳先进技术研究院 辅助运动系统及其控制方法
CN102139139A (zh) * 2011-01-13 2011-08-03 中国医学科学院生物医学工程研究所 一种肌电反馈控制的电刺激装置及其控制方法
CN102426651A (zh) * 2011-08-25 2012-04-25 武汉理工大学 人体前臂表面肌电信号采集及模式识别系统
CN102961203A (zh) * 2012-12-10 2013-03-13 杭州电子科技大学 基于emd样本熵的表面肌电信号识别方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHRISTOS SAPSANIS 等: "EMG based classification of basic hand movements based on timefrequency features", 《2013 21ST MEDITERRANEAN CONFERENCE ON CONTROL & AUTOMATION》 *
DORI PELEG 等: "Classification of Finger Activation for Use in a Robotic Prosthesis Arm", 《IEEE TRANSACTIONS ON NEURAL SYSTEMS AND REHABILITATION ENGINEERING》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106108897A (zh) * 2016-07-20 2016-11-16 西安中科比奇创新科技有限责任公司 一种基于经验模式分解的肌电信号滤波方法
CN107126303A (zh) * 2017-02-15 2017-09-05 上海术理智能科技有限公司 一种基于手机app的上下肢运动辅助方法
CN109222968A (zh) * 2017-07-10 2019-01-18 丰田自动车株式会社 康复评估设备、康复评估方法以及康复评估程序
CN108108016A (zh) * 2017-12-07 2018-06-01 浙江大学 手势感知器
CN108268844A (zh) * 2018-01-17 2018-07-10 上海术理智能科技有限公司 基于表面肌电信号的动作模式识别方法及装置
CN108403108A (zh) * 2018-02-11 2018-08-17 宁波工程学院 基于波形优化的阵列式表面肌电信号分解方法
CN109062401B (zh) * 2018-07-11 2020-12-15 北京理工大学 一种基于肌电信号的实时手势识别系统
CN109062401A (zh) * 2018-07-11 2018-12-21 北京理工大学 一种基于肌电信号的实时手势识别系统
CN109512424B (zh) * 2018-11-16 2021-07-13 福州大学 一种高密度或多通道肌电信号的肌肉激活起点检测方法
CN109512424A (zh) * 2018-11-16 2019-03-26 福州大学 一种高密度或多通道肌电信号的肌肉激活起点检测方法
CN109948465B (zh) * 2019-02-28 2020-11-13 浙江工业大学 一种基于多目标的表面肌电信号分类方法
CN109948465A (zh) * 2019-02-28 2019-06-28 浙江工业大学 一种基于多目标的表面肌电信号分类方法
CN110399846A (zh) * 2019-07-03 2019-11-01 北京航空航天大学 一种基于多通道肌电信号相关性的手势识别方法
CN113986017A (zh) * 2021-12-27 2022-01-28 深圳市心流科技有限公司 一种肌电手势模板的生成方法、装置及存储介质
CN113986017B (zh) * 2021-12-27 2022-05-17 深圳市心流科技有限公司 一种肌电手势模板的生成方法、装置及存储介质
CN115299950A (zh) * 2022-08-12 2022-11-08 歌尔股份有限公司 肌电信号采集电路、可穿戴设备及控制方法

Also Published As

Publication number Publication date
EP3167799A1 (en) 2017-05-17
CN105326500B (zh) 2018-02-09
EP3167799A4 (en) 2017-08-09
US20170143226A1 (en) 2017-05-25
WO2016023349A1 (zh) 2016-02-18

Similar Documents

Publication Publication Date Title
CN105326500A (zh) 基于表面肌电信号的动作识别方法和设备
Anderson et al. Multivariate autoregressive models for classification of spontaneous electroencephalographic signals during mental tasks
Phinyomark et al. Feature reduction and selection for EMG signal classification
CN102573619B (zh) 用于生成注意力水平表示的设备和方法
CN107273798A (zh) 一种基于表面肌电信号的手势识别方法
CN108268844A (zh) 基于表面肌电信号的动作模式识别方法及装置
Fatima et al. Towards a low cost Brain-computer Interface for real time control of a 2 DOF robotic arm
CN104571504A (zh) 一种基于想象动作的在线脑-机接口方法
Goen et al. Review of surface electromyogram signals: its analysis and applications
Velásquez-Martínez et al. Motor imagery classification for BCI using common spatial patterns and feature relevance analysis
Zhang et al. Activity recognition of the torso based on surface electromyography for exoskeleton control
Altindis et al. Use of topological data analysis in motor intention based brain-computer interfaces
CN106484112A (zh) 基于运动想象脑‑机接口的字符拼写器及测试方法
CN105686827A (zh) 一种基于微控制器的肌电信号处理及特征提取方法
Vargic et al. Human computer interaction using BCI based on sensorimotor rhythm
Sourina et al. Fractal-based brain state recognition from EEG in human computer interaction
CN115299960A (zh) 基于短时变分模态分解的电信号分解方法及脑电信号分解装置
Yu et al. Classification of nine directions using the maximum likelihood estimation based on electromyogram of both forearms
Chen et al. Analysis of motor unit activities during multiple motor tasks by real-time emg decomposition: perspective for myoelectric control
Dyson et al. Sequential classification of mental tasks vs. idle state for EEG based BCIs
Nakaya et al. A practical approach for recognition of hand gesture and distinction of its singularity
Lee et al. Feature mapping using PCA, locally linear embedding and isometric feature mapping for EEG-based brain computer interface
Kamath et al. Modeling QRS complex in dECG
Sudalairaj Spatio-Temporal Analysis of EEG using Deep Learning
Stival Subject-independent modeling of semg signals for the motion of a single robot joint through gmm modelization

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180209

Termination date: 20190813

CF01 Termination of patent right due to non-payment of annual fee