CN105943042A - 基于肌电信号的操作者手部行为感知系统 - Google Patents

基于肌电信号的操作者手部行为感知系统 Download PDF

Info

Publication number
CN105943042A
CN105943042A CN201610396781.XA CN201610396781A CN105943042A CN 105943042 A CN105943042 A CN 105943042A CN 201610396781 A CN201610396781 A CN 201610396781A CN 105943042 A CN105943042 A CN 105943042A
Authority
CN
China
Prior art keywords
signal
electromyographic
operator
electromyographic signal
unit
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.)
Pending
Application number
CN201610396781.XA
Other languages
English (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201610396781.XA priority Critical patent/CN105943042A/zh
Publication of CN105943042A publication Critical patent/CN105943042A/zh
Pending legal-status Critical Current

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
    • A61B5/389Electromyography [EMG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/1118Determining activity level
    • 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/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6813Specially adapted to be attached to a specific body part
    • A61B5/6824Arm or wrist
    • 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/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • 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/7225Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • 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/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2503/00Evaluating a particular growth phase or type of persons or animals
    • A61B2503/20Workers
    • A61B2503/22Motor vehicles operators, e.g. drivers, pilots, captains

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Power Engineering (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

一种基于肌电信号的操作者手部行为感知系统,包括肌电腕表装置和接收处理终端;所述肌电腕表装置包括弹性腕带以及设置在弹性腕带上的各功能模块,功能模块包括肌电电极、信号调理单元、信号处理单元、无线发送单元和供电单元;所述接收处理终端包括无线接收单元以及信号分析模块;信号调理单元对肌电电极采集的肌电信号进行放大和滤波处理后输出给信号处理单元进行模数转换并缓存,处理后的肌电信号数据经无线发送单元将发送给接收处理终端,接收处理终端能够对肌电信号数据进行处理和分析,得到操作者的操作手势。本发明能够有效识别操作者的手势信息,并且具有信噪比高、能耗低等优点。

Description

基于肌电信号的操作者手部行为感知系统
技术领域
本发明涉及在轨服务技术领域,具体说涉及一种基于肌电信号的操作者手部行为感知系统。
背景技术
随着空间研究、开发与应用能力的不断提高,各国相继研制并发射了大量面向各种任务要求的航天器,航天器的结构、组成日趋复杂,性能、技术水平不断提高。在这种情况下,宇航员的工作时间越来越强,工作强度也越来越大,因此保证宇航员的安全是空间任务的重中之重。
宇航员在太空中,各种生理信号需要在线实时监控,其中肌电信号是一种重要的生理信号,它既可以通过肌电信号分析得到手部动作,又可以反映宇航员肌肉的疲劳程度和身体健康程度,甚至还可以从信号中分析宇航员的操作动作正确与否。
而现有的肌电信号采集仪大多没有功耗和体积指标,而且由于体积太大导致不能腕表在宇航服内,而且采用有线连接方式不便实时数据处理,而航天应用对功耗和体积要求非常严格,而且对系统的抗噪声干扰能力也要求非常高,因此目前还没有适合航天员应用的、体积小、易佩戴、功能高度集成的操作者手部行为感知系统。
发明内容
为了克服现有技术存在的不足,本发明提供一种基于肌电信号的操作者手部行为感知系统,适应航天应用场合的操作者手部行为感知。
本发明的技术方案是:
一种基于肌电信号的操作者手部行为感知系统,包括肌电腕表装置和接收处理终端;所述肌电腕表装置包括弹性腕带以及设置在弹性腕带上的各功能模块,这些功能模块包括肌电电极、信号调理单元、信号处理单元、无线发送单元、供电单元;所述接收处理终端包括无线接收单元以及信号分析模块;
所述肌电电极用于采集表面肌电信号;所述信号调理单元对肌电电极采集的肌电信号进行放大和滤波处理后输出给信号处理单元;所述信号处理单元将信号进行模数转换并缓存,同时将处理后的肌电信号数据输出给无线发送单元;所述无线发送单元将接收到的肌电信号数据以无线的方式发送给无线接收单元;所述供电单元为上述肌电电极、信号调理单元、信号处理单元、无线发送单元提供电源;所述无线接收单元能够接收来自无线发送单元的肌电信号数据并提交给信号分析模块;所述信号分析模块对肌电信号数据进行处理和分析,得到操作者的操作手势。
本发明中:所述肌电电极为多个,各肌电电极分布在腕带的不同位置,用于采集对应位置的肌电信号。进一步地,所述的肌电电极采用贴片式Ag/AgCl电极,降低了电阻对肌电信号采集的影响,提高了信噪比。
本发明中:所述的信号调理单元为一个两级信号放大模块,具有信号放大功能,此外它还具有陷波、低通和高通等滤波功能,放大倍数为20000倍,长宽高为450mm×60mm×4mm,功耗小于25mW。
本发明中:所述的信号处理单元包括微处理器和模数转换器(ADC),其中微处理器采用超低功耗的MSP430F2418,主晶振采用16MHz晶振,降低了能耗。模数转换器(ADC)采用微处理器自带的功能模块,基准采用低功耗的外部电源基准max1633,它初始精度为0.04%,具有低温漂系数,低压降,低静态电流和低噪声的优点。
本发明中:所述的无线发送单元采用A7105芯片。
本发明中:所述的供电单元包括锂电池和电源模块,其中电源模块采用低压降的ADP122,它电源纹波小,可以使得系统只需要一节锂电池供电即可正常输出3.3V电压,并且利用磁珠实现了数字电源和模拟电源的隔离。
本发明中:所述的无线接收单元采用2.4GHz的无线路由器。
本发明中:所述的信号分析模块是一台PC机,加载了手势识别算法。
本发明中:所述腕带采用有弹性的树脂材料制作,可以用于适应不同大小的手腕佩戴;所述腕带的两端头设置有佩戴锁紧件,佩戴锁紧件采用小型公母搭扣方式实现。
本发明提供一种手势识别算法,包括以下步骤:
S1:手部行为发起检测,设某个肌电信号采集窗口的时间长度为T,500Hz采集频率下共采集到N个肌电信号值X
S1.1:零均值处理
其中,X(n)为该肌电信号采集窗口中采集到的肌电信号的第n个采样值,为计算得到的肌电信号均值;
S1.2:计算TKE算子:
其中,为步骤S1.1中计算得到的肌电信号均值,ψ(n)为计算得到的TKE算子(该算子参见参考文献:Teager H M,Teager S M.Evidence for Nonlinear SoundReduction Mechanisms in the Vocal Tract.Kluwer Acad.Publ.,1990:241-261.)。
S1.3:计算均值和标准差:
μ 0 = 1 N Σ n = 1 N ψ ( n )
δ 0 = { 1 N - 1 Σ n = 1 N ( ψ ( n ) - μ 0 ) 2 } 1 / 2
其中,ψ(n)为步骤S1.2计算得到的TKE算子,μ0为计算得到的均值,δ0为计算得到的标准差。
S1.4:计算阈值Th:Th=μ0+j·δ0
其中μ0为步骤S1.3计算得到的均值,δ0为步骤S1.3计算得到的标准差,j为事先设定的阈值乘子,Th为计算得到的阈值;
S1.5:计算动作序列s(n):s(n)=sign(ψ(n)-Th)
其中,sign(·)为符号函数,即当ψ(n)-Th>0时s(n)=1,否则s(n)=0。
S1.6:动作序列s(n)滤波:设定T1为手部行为中出现的非活动时长,即在正常的肌肉活动中表现出的静息状态的时间;T2则为手部静止时,由于尖峰噪声等影响而出现的活动时长。
在得到的一系列s(n)中,若相邻两个s(n)=1时刻之差小于T1,将两个s(n)=1时刻之间的s(n)的值全部更改为“1”,这样做的目的就是,避免在肌肉活动时由于收缩过快等原因偶尔出现的静息状态;
在得到的一系列s(n)中,相邻两个s(n)=0时刻之差小于T2,将两个s(n)=0时刻之间的s(n)的值全部更改为“0”,这样做的目的就是,在正常无动作情况下肌肉静息时偶尔出现的噪声等尖峰信号带来的影响。
S1.7:计算动作发起时间:
S2:建立常用手势库;
所述的常用手势库,包括以下8种手势,但并不仅仅限于此8种手势,相似的改进和修饰应视作与该手势库相同:(1)用五指抓住并托举书本,且大拇指与其余四指方向相对;(2)仅使用大拇指、食指和中指,三指抓住盛满水的500mL易拉罐;(3)用五指抓住盛满水的500ml易拉罐,且大拇指与其余四指方向相对;(4)用五指抓住球状物;(5)仅使用大拇指与食指拿捏CD盒;(6)用大拇指、食指和中指摘掉记号笔的笔盖;(7)用五指打开文具盒;(8)用五指拿起铅笔。
S3:预分类,包括基于积分特征值的特征提取和基于线性判别函数的模式识别两个步骤。积分特征值计算如下:
I E M G = 1 N Σ n = 0 N - 1 | E M G ( n ) |
其中,EMG(n)(n=0,1,2…,N-1)为一长度为N的时间序列,在这里为一列肌电信号。
基于线性判别函数的模式识别中的识别模型可表述为:
d ( x ) = d e f ω ′ x
其中x=(x1,x2,…,xn,1)′,ω=(ω12,…,ωnn+1)′,x为增广特征矢量,这里
x1,x2,…,xn分别指各通道肌电信号的积分值IEMG,ω称为增广权矢量。ω的学习过程可描述为:
其中,ω(k)为第k步迭代得到的增广权矢量ω,ρ为预先设定的学习率,为一个正实数。
当学习集的模式全部用作训练完毕后,对其进行检验,若仍有模式被错判,则需要对其开展再次迭代,进行第二轮训练,调整建立新的权矢量。以此类推,直到对于学习集中的模式分类正确率达到98%以上为止,或者迭代M次为止,此时得到的权矢量ω值即为所求。
当确定好判别函数d(x)后,即ω确定后,将测试集的特征矢量xt,即各通道肌电信号的积分值IEMG,并在该特征矢量的基础上增加一个常数维度,构成增广特征矢量x,代入判别函数中,运用如下判别规则进行判断,达到分类的效果:
d ( x ) = &omega; x &GreaterEqual; 0 &RightArrow; x &Element; A < 0 &RightArrow; x &Element; B
其中,d(x)为学习得到的判别函数,ω为学习得到的增广权矢量,x为增广特征矢量;
S4:细分类,包括基于AR模型的特征提取和基于BP神经网络分类器的模式识别两个步骤。其中AR模型可表述如下:
E M G ( n ) + &Sigma; k = 1 p a k E M G ( n - k ) = w ( n )
其中EMG(n)为一系列采集到的肌电信号的第n个采样值;w(n)为输入信号,假定其为白噪声,p为AR模型的阶数。ak则是AR模型的第k个系数,这便是所要求取的特征参数。
这里的BP神经网络分类器其输入为AR模型的特征参数,输出为手势的类别。通过训练集的学习过程得到的BP神经网络权值,输入AR模型的特征参数从而可得到手势类别。
本发明主要技术指标为:
(1)供电电压:一节3.7V锂电池;
(2)肌电信号放大倍数:20000倍;
(3)系统功耗:小于200mW;
(4)信号调理模块几何参数:长宽厚为450mm×60mm×4mm
(5)常用手势库中的8种手势的识别正确率>80%
本发明的优点在于:提供的可佩带式的低功耗表面肌电信号采集仪,采用低功耗的信号放大以及采集模块,处理器采用超低功耗的MSP430F2418,减小了体积和功耗。另外采用了低压降的电源模块以减小电源功耗和减小电源纹波,并考虑数字电源和模拟电源的隔离,改善了系统的信噪比,增强了系统的稳定性和可靠性。本发明能够有效识别操作者的手势信息,并且具有信噪比高、能耗低等优点,适应航天应用场合的操作者手部行为感知。
附图说明
图1是基于肌电信号的操作者手部行为感知系统结构示意图。
图2是基于肌电信号的操作者手部行为感知系统原理方框图
图3电源转换电路原理图
图4电源基准电路原理图
图5手势识别算法示意图
图6操作手势的预分类识别率
图7操作手势细分类的混淆矩阵
图8操作手势总的识别率
图9手部行为发起检测的流程
具体实施方式
为了使本发明所解决的技术问题、技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1所示,本发明提供一种基于肌电信号的操作者手部行为感知系统,其特征在于,包括:肌电电极2、信号调理单元4、信号处理单元6、无线发送单元5、腕带3、佩戴锁紧件1、供电单元7、无线接收单元8和信号分析模块9。
图1所示的腕带3采用有弹性的树脂材料制作,可以用于适应不同大小的手腕佩戴,佩戴锁紧件采用小型公母搭扣方式实现。
图2给出了基于肌电信号的操作者手部行为感知系统的工作原理,肌电电极2采用贴片式Ag/AgCl电极,用于采集表面肌电信号,信号调理单元对肌电电极采集的肌电信号进行放大和滤波处理后输出给信号处理单元;所述信号处理单元将信号进行模数转换并缓存,同时将处理后的肌电信号数据输出给无线发送单元;所述无线发送单元将接收到的肌电信号数据以无线的方式发送给无线接收单元。无线接收单元能够接收来自无线发送单元的肌电信号数据并提交给信号分析模块;所述信号分析模块内加载有手势识别算法能够对肌电信号数据进行处理和分析,得到操作者的操作手势。
信号处理单元包括微处理器和模数转换器,微处理器采用低功耗MSP430F2418,发送模块采用低功耗的nrf24l01,微处理器处理器自带的ADC。所述的电压基准采用外部低功耗基准电源,电源由一节锂电池经过电源转换芯片后提供。
图1中所述的信号调理单元4为一个两级信号放大模块,它还有陷波、低通和高通的功能,放大倍数为20000倍,长宽高为450mm×60mm×4mm,功耗小于25mW,只需要提供+3.3V电源。
图3和图4分别描述了电源模块的电路图。所述的电源模块采用低压降的ADP122,它电源纹波小,可以使得系统只需要一节锂电池供电即可正常输出3.3V电压,C2和C3是退耦电容,L1和L2是磁珠,用于隔离数字电源和模拟电源,磁珠L3是用于隔离数字地和模拟地的作用,R1和R2用于对电池分压后检测电池电量。
信号处理单元包括微处理器和ADC。所述的微处理器采用超低功耗的MSP430F2418,主晶振采用16MHz晶振,它带有1个SPI接口、两个URAT接口,10通道的12bit的ADC接口,以及带有RTC模块和看门狗模块,具有多种低功耗模式,可以关闭不需要的外设模块以节约电量,并且有多个外部IO可以使用。所述的ADC采用微处理器自带的功能模块,基准采用低功耗的外部电源基准max1633,它初始精度为0.04%,具有低温漂系数,低压降,低静态电流和低噪声的优点。
如图5所示,本发明提供一种手势识别算法,包括以下步骤:
S1:手部行为发起检测;
S2:建立常用手势库;
所述的常用手势库,包括以下8种手势,但并不仅仅限于此8种手势,相似的改进和修饰应视作与该手势库相同:(1)用五指抓住并托举书本,且大拇指与其余四指方向相对;(2)仅使用大拇指、食指和中指,三指抓住盛满水的500mL易拉罐;(3)用五指抓住盛满水的500ml易拉罐,且大拇指与其余四指方向相对;(4)用五指抓住球状物;(5)仅使用大拇指与食指拿捏CD盒;(6)用大拇指、食指和中指摘掉记号笔的笔盖;(7)用五指打开文具盒;(8)用五指拿起铅笔。
S3:预分类,包括基于积分特征值的特征提取和基于线性判别函数的模式识别两个步骤。积分特征值计算如下:
I E M G = 1 N &Sigma; n = 0 N - 1 | E M G ( n ) |
其中,EMG(n)(n=0,1,2…,N-1)为一长度为N的时间序列,在这里为一列肌电信号。
基于线性判别函数的模式识别中的识别模型可表述为:
d ( x ) = d e f &omega; &prime; x
其中x=(x1,x2,…,xn,1)′,ω=(ω12,…,ωnn+1)′,x为增广特征矢量,这里x1,x2,…,xn分别指各肌电信号通路的积分值IEMG,ω称为增广权矢量。ω的学习过程可描述为:
当学习集的模式全部用作训练完毕后,对其进行检验,若仍有模式被错判,则需要对其开展再次迭代,进行第二轮训练,调整建立新的权矢量;以此类推,直到对于学习集中的模式分类正确率达到98%以上为止,或者迭代M次为止,此时得到的权矢量ω值即为所求。
当确定好判别函数d(x)后,即ω确定后,将测试集的特征矢量x即各通道肌电信号的积分值IEMG,代入判别函数中,运用如下判别规则进行判断,达到分类的效果:
d ( x ) = &omega; &prime; x &GreaterEqual; 0 &RightArrow; x &Element; A < 0 &RightArrow; x &Element; B
S4:细分类,包括基于AR模型的特征提取和基于BP神经网络分类器的模式识别两个步骤。其中AR模型可表述如下:
E M G ( n ) + &Sigma; k = 1 p a k E M G ( n - k ) = w ( n )
其中EMG(n)为线性系统的输出信号,其表示为一系列采集到的肌电信号的第n个采样值;w(n)为输入信号,假定其为白噪声,p为AR模型的阶数。ak则是AR模型的第k个系数,这便是所要求取的特征参数。
这里的BP神经网络分类器其输入为AR模型的特征参数,输出为手势的类别。通过训练集的学习过程得到的BP神经网络权值,输入AR模型的特征参数从而可得到手势类别。
图6给出了操作手势的预分类识别率,图7给出了操作手势细分类的混淆矩阵,两者相乘可得到图8所示的操作手势总的识别率。常用手势库中的手势识别率达到了82.8%以上,证明了本发明所提手势识别方法的有效性。
手部行为发起检测的流程如图9所示,具体包括以下步骤:
S1.1:零均值处理
S1.2:计算TKE算子:
S1.3:计算均值和标准差:
&mu; 0 = 1 N &Sigma; n = 1 N &psi; ( n )
&delta; 0 = { 1 N - 1 &Sigma; n = 1 N ( &psi; ( n ) - &mu; 0 ) 2 } 1 / 2
S1.4:计算阈值Th:Th=μ0+j·δ0
S1.5:计算动作序列s(n):s(n)=sign(ψ(n)-Th)
S1.6:动作序列s(n)滤波:设定T1为手部行为中,出现的非活动时长,即在正常的肌肉活动中表现出的静息状态的时间;T2则为手部静止时,由于尖峰噪声等影响而出现的活动时长。
在得到的一系列s(n)中,若相邻两个s(n)=1时刻之差小于T1,将两个s(n)=1时刻之间的s(n)的值全部更改为“1”,这样做的目的就是,避免在肌肉活动时由于收缩过快等原因偶尔出现的静息状态;
在得到的一系列s(n)中,相邻两个s(n)=0时刻之差小于T2,将两个s(n)=0时刻之间的s(n)的值全部更改为“0”,这样做的目的就是,在正常无动作情况下肌肉静息时偶尔出现的噪声等尖峰信号带来的影响。
S1.7:计算动作发起时间:
以上所述仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应该提出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (9)

1.一种基于肌电信号的操作者手部行为感知系统,其特征在于:包括肌电腕表装置和接收处理终端;所述肌电腕表装置包括弹性腕带以及设置在弹性腕带上的各功能模块,功能模块包括肌电电极、信号调理单元、信号处理单元、无线发送单元和供电单元;所述接收处理终端包括无线接收单元以及信号分析模块;
所述肌电电极用于采集表面肌电信号;所述信号调理单元对肌电电极采集的肌电信号进行放大和滤波处理后输出给信号处理单元;所述信号处理单元将信号进行模数转换并缓存,同时将处理后的肌电信号数据输出给无线发送单元;所述无线发送单元将接收到的肌电信号数据以无线的方式发送给无线接收单元;所述供电单元为上述肌电电极、信号调理单元、信号处理单元、无线发送单元提供电源;所述无线接收单元能够接收来自无线发送单元的肌电信号数据并提交给信号分析模块;所述信号分析模块内加载有手势识别算法能够对肌电信号数据进行处理和分析,得到操作者的操作手势。
2.根据权利要求1所述的基于肌电信号的操作者手部行为感知系统,其特征在于:所述肌电电极为多个,各肌电电极分布在腕带的不同位置,用于采集对应位置的肌电信号。
3.根据权利要求2所述的基于肌电信号的操作者手部行为感知系统,其特征在于:所述的肌电电极采用贴片式Ag/AgCl电极。
4.根据权利要求1所述的基于肌电信号的操作者手部行为感知系统,其特征在于:信号处理单元包括微处理器和模数转换器,其中微处理器采用MSP430F2418,主晶振采用16MHz晶振;模数转换器采用微处理器自带的功能模块,基准采用外部电源基准max1633,其初始精度为0.04%。
5.根据权利要求1所述的基于肌电信号的操作者手部行为感知系统,其特征在于:所述无线发送单元采用A7105芯片。
6.根据权利要求1所述的基于肌电信号的操作者手部行为感知系统,其特征在于:所述供电单元包括锂电池和电源模块,其中电源模块采用ADP122。
7.根据权利要求1所述的基于肌电信号的操作者手部行为感知系统,其特征在于:所述无线接收单元采用2.4GHz的无线路由器。
8.根据权利要求1所述的基于肌电信号的操作者手部行为感知系统,其特征在于:所述腕带采用有弹性的树脂材料制作;所述腕带的两端头设置有佩戴锁紧件。
9.一种手势识别算法,其特征在于,包括以下步骤:
S1:手部行为发起检测
设某个肌电信号采集窗口的时间长度为T,500Hz采集频率下共采集到N个肌电信号值x;
S1.1:零均值处理
其中,X(n)为该肌电信号采集窗口中采集到的肌电信号的第n个采样值,为计算得到的肌电信号均值;
S1.2:计算TKE算子:
其中,为步骤S1.1中计算得到的肌电信号均值,ψ(n)为计算得到的TKE算子;
S1.3:计算均值和标准差:
&mu; 0 = 1 N &Sigma; n = 1 N &psi; ( n )
&delta; 0 = { 1 N - 1 &Sigma; n = 1 N ( &psi; ( n ) - &mu; 0 ) 2 } 1 / 2
其中,ψ(n)为步骤S1.2计算得到的TKE算子,μ0为计算得到的均值,δ0为计算得到的标准差;
S1.4:计算阈值Th:Th=μ0+j·δ0
其中μ0为步骤S1.3计算得到的均值,δ0为步骤S1.3计算得到的标准差,j为事先设定的阈值乘子,Th为计算得到的阈值;
S1.5:计算动作序列s(n):s(n)=sign(ψ(n)-Th)
其中,sign(·)为符号函数,即当ψ(n)-Th>0时s(n)=1,否则s(n)=0;
S1.6:动作序列s(n)滤波:设定T1为手部行为中出现的非活动时长,即在正常的肌肉活动中表现出的静息状态的时间;T2则为手部静止时,由于尖峰噪声等影响而出现的活动时长;
在得到的一系列s(n)中,若相邻两个s(n)=1时刻之差小于T1,将两个s(n)=1时刻之间的s(n)的值全部更改为“1”;
在得到的一系列s(n)中,相邻两个s(n)=0时刻之差小于T2,将两个s(n)=0时刻之间的s(n)的值全部更改为“0”;
S1.7:计算动作发起时间:
S2:建立常用手势库;
S3:预分类,包括基于积分特征值的特征提取和基于线性判别函数的模式识别两个步骤;
积分特征值计算如下:
I E M G = 1 N &Sigma; n = 0 N - 1 | E M G ( n ) |
其中,EMG(n)(n=0,1,2…,N-1)为一长度为N的一列肌电信号;
基于线性判别函数的模式识别中的识别模型可表述为:
d ( x ) = d e f &omega; &prime; x
其中x=(x1,x2,…,xn,1)′,ω=(ω12,…,ωnn+1)′,x为增广特征矢量,x1,x2,…,xn分别指各通道肌电信号的积分值IEMG;ω称为增广权矢量;ω的学习过程可描述为:
其中,ω(k)为第k步迭代得到的增广权矢量ω,ρ为预先设定的学习率,为一个正实数;
当学习集的模式全部用作训练完毕后,对其进行检验,若仍有模式被错判,则需要对其开展再次迭代,进行第二轮训练,调整建立新的权矢量;以此类推,直到对于学习集中的模式分类正确率达到98%以上为止,或者迭代M次为止,此时得到的权矢量ω值即为所求;
当确定好判别函数d(x)后,即ω确定后,测试集的特征矢量xt,即各通道肌电信号的积分值IEMG,并在该特征矢量的基础上增加一个常数维度,构成增广特征矢量x,代入判别函数中,运用如下判别规则进行判断,达到分类的效果:
d ( x ) = &omega; x &GreaterEqual; 0 &RightArrow; x &Element; A < 0 &RightArrow; x &Element; B
其中,d(x)为学习得到的判别函数,ω为学习得到的增广权矢量,x为增广特征矢量;
S4:细分类,包括基于AR模型的特征提取和基于BP神经网络分类器的模式识别两个步骤;
AR模型可表述如下:
E M G ( n ) + &Sigma; k = 1 p a k E M G ( n - k ) = w ( n )
其中EMG(n)为一系列采集到的肌电信号的第n个采样值;w(n)为输入信号,假定其为白噪声,p为AR模型的阶数;ak则是AR模型的第k个系数即待求取的特征参数;
BP神经网络分类器其输入为AR模型的特征参数,输出为手势的类别;通过训练集的学习过程得到的BP神经网络权值,输入AR模型的特征参数从而可得到手势类别。
CN201610396781.XA 2016-06-07 2016-06-07 基于肌电信号的操作者手部行为感知系统 Pending CN105943042A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610396781.XA CN105943042A (zh) 2016-06-07 2016-06-07 基于肌电信号的操作者手部行为感知系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610396781.XA CN105943042A (zh) 2016-06-07 2016-06-07 基于肌电信号的操作者手部行为感知系统

Publications (1)

Publication Number Publication Date
CN105943042A true CN105943042A (zh) 2016-09-21

Family

ID=56907803

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610396781.XA Pending CN105943042A (zh) 2016-06-07 2016-06-07 基于肌电信号的操作者手部行为感知系统

Country Status (1)

Country Link
CN (1) CN105943042A (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108717327A (zh) * 2018-05-21 2018-10-30 中北大学 一种用于反恐作战手势识别的智能感知手套及感知方法
CN109062401A (zh) * 2018-07-11 2018-12-21 北京理工大学 一种基于肌电信号的实时手势识别系统
CN109077726A (zh) * 2018-10-29 2018-12-25 博睿康科技(常州)股份有限公司 一种肌电检测装置与肌电识别方法
CN110720910A (zh) * 2019-10-12 2020-01-24 宁波工程学院 一种基于相关性的肌肉运动单元搜索方法
CN111723717A (zh) * 2020-06-11 2020-09-29 天津大学 一种无声语音识别方法及系统
CN111857349A (zh) * 2020-07-28 2020-10-30 中国科学技术大学 一种带自修复和自校准功能的腕带式手势识别设备及方法
CN112401906A (zh) * 2020-11-10 2021-02-26 河北省科学院应用数学研究所 一种基于幅值的自动脑电去伪迹方法
CN116172584A (zh) * 2022-12-14 2023-05-30 博灵脑机(杭州)科技有限公司 一种自适应的多通道肌电信号实时质量评估与检测方法
US11705748B2 (en) 2018-02-28 2023-07-18 Russell Wade Chan Wearable gesture recognition device for medical screening and associated operation method and system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102846317A (zh) * 2012-07-18 2013-01-02 上海交通大学 佩戴式肌电信号采集系统
CN103750836A (zh) * 2014-01-14 2014-04-30 优尔美有限公司 一种可穿戴式肌电设备
CN103777752A (zh) * 2013-11-02 2014-05-07 上海威璞电子科技有限公司 基于手臂肌肉电流检测和运动传感器的手势识别设备
CN103853333A (zh) * 2014-03-21 2014-06-11 上海威璞电子科技有限公司 一种手势控制玩具的方案
KR20140107741A (ko) * 2013-02-28 2014-09-05 주식회사메디아나 근전도 신호를 이용한 정보 입력 장치 및 그 방법

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102846317A (zh) * 2012-07-18 2013-01-02 上海交通大学 佩戴式肌电信号采集系统
KR20140107741A (ko) * 2013-02-28 2014-09-05 주식회사메디아나 근전도 신호를 이용한 정보 입력 장치 및 그 방법
CN103777752A (zh) * 2013-11-02 2014-05-07 上海威璞电子科技有限公司 基于手臂肌肉电流检测和运动传感器的手势识别设备
CN103750836A (zh) * 2014-01-14 2014-04-30 优尔美有限公司 一种可穿戴式肌电设备
CN103853333A (zh) * 2014-03-21 2014-06-11 上海威璞电子科技有限公司 一种手势控制玩具的方案

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DAPENG YANG ET AL: "EMG Onset Detection Based on Teager-Kaiser Energy Operator and Morphological Close Operation", 《INTERNATIONAL CONFERENCE ON INTELLIGENT ROBOTICS AND APPLICATIONS》 *
赵漫丹 等: "基于肌电信号层级分类的手部动作识别方法", 《北京生物医学工程》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11705748B2 (en) 2018-02-28 2023-07-18 Russell Wade Chan Wearable gesture recognition device for medical screening and associated operation method and system
CN108717327B (zh) * 2018-05-21 2021-01-01 中北大学 一种用于反恐作战手势识别的智能感知手套及感知方法
CN108717327A (zh) * 2018-05-21 2018-10-30 中北大学 一种用于反恐作战手势识别的智能感知手套及感知方法
CN109062401A (zh) * 2018-07-11 2018-12-21 北京理工大学 一种基于肌电信号的实时手势识别系统
CN109077726A (zh) * 2018-10-29 2018-12-25 博睿康科技(常州)股份有限公司 一种肌电检测装置与肌电识别方法
CN110720910A (zh) * 2019-10-12 2020-01-24 宁波工程学院 一种基于相关性的肌肉运动单元搜索方法
CN110720910B (zh) * 2019-10-12 2022-03-29 宁波工程学院 一种基于相关性的肌肉运动单元搜索方法
CN111723717A (zh) * 2020-06-11 2020-09-29 天津大学 一种无声语音识别方法及系统
CN111857349A (zh) * 2020-07-28 2020-10-30 中国科学技术大学 一种带自修复和自校准功能的腕带式手势识别设备及方法
CN111857349B (zh) * 2020-07-28 2022-03-01 中国科学技术大学 一种带自修复和自校准功能的腕带式手势识别设备及方法
CN112401906A (zh) * 2020-11-10 2021-02-26 河北省科学院应用数学研究所 一种基于幅值的自动脑电去伪迹方法
CN112401906B (zh) * 2020-11-10 2021-12-14 河北省科学院应用数学研究所 一种基于幅值的自动脑电去伪迹方法
CN116172584A (zh) * 2022-12-14 2023-05-30 博灵脑机(杭州)科技有限公司 一种自适应的多通道肌电信号实时质量评估与检测方法
CN116172584B (zh) * 2022-12-14 2024-01-23 博灵脑机(杭州)科技有限公司 一种自适应的多通道肌电信号实时质量评估与检测方法

Similar Documents

Publication Publication Date Title
CN105943042A (zh) 基于肌电信号的操作者手部行为感知系统
CN104398263B (zh) 一种基于近似熵和互近似熵的帕金森患者震颤症状量化评测方法
Akmandor et al. Smart, secure, yet energy-efficient, Internet-of-Things sensors
Shi et al. Mobile human airbag system for fall protection using MEMS sensors and embedded SVM classifier
Wu et al. Pilots’ fatigue status recognition using deep contractive autoencoder network
US8301575B2 (en) Method and device for the recognition of the position or movement of a device or a person
CN104569694B (zh) 一种面向航天器飞行过程的电信号特征提取和识别系统
CN106803081A (zh) 一种基于多分类器集成的脑电分类方法
CN103126657B (zh) 一种监测指腕疲劳与工作时间的装置与方法
You et al. Finger motion decoding using EMG signals corresponding various arm postures
CN110123314A (zh) 基于脑电信号判断大脑专注放松状态的方法
Patel et al. Low power real-time seizure detection for ambulatory EEG
CN104490391A (zh) 一种基于脑电信号的战斗人员状态监测系统
Zhang et al. Real-time implementation of a self-recovery EMG pattern recognition interface for artificial arms
Imran et al. Hharnet: Taking inspiration from inception and dense networks for human activity recognition using inertial sensors
Xue et al. SEMG-based human in-hand motion recognition using nonlinear time series analysis and random forest
CN103500342A (zh) 一种基于加速度计的人体行为识别方法
CN106991409A (zh) 一种运动想象脑电信号特征提取与分类系统及方法
CN106951852A (zh) 一种有效的人体行为识别系统
Chang et al. A hierarchical hand motions recognition method based on IMU and sEMG sensors
Ding et al. Cascaded convolutional neural network with attention mechanism for mobile EEG-based driver drowsiness detection system
CN102438522A (zh) 用于确定移动元素的活动的系统和方法
Lu et al. Compressed signal processing on Nyquist-sampled signals
CN105138997A (zh) 一种无线体域网中基于压缩分类的动作识别方法
Mathan et al. Neurophysiological estimation of interruptibility: Demonstrating feasibility in a field context

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20160921