CN104997508B - 阵列式sEMG信号自动分解方法 - Google Patents
阵列式sEMG信号自动分解方法 Download PDFInfo
- Publication number
- CN104997508B CN104997508B CN201510512579.4A CN201510512579A CN104997508B CN 104997508 B CN104997508 B CN 104997508B CN 201510512579 A CN201510512579 A CN 201510512579A CN 104997508 B CN104997508 B CN 104997508B
- Authority
- CN
- China
- Prior art keywords
- semg
- array
- waveform
- sequence vector
- moving cell
- 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.)
- Expired - Fee Related
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/389—Electromyography [EMG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7221—Determining signal validity, reliability or quality
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- Psychiatry (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Fuzzy Systems (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明提供一种阵列式sEMG信号自动分解方法,首先对阵列式sEMG信号预处理,并且计算初始发放序列向量;其次利用时域减法确定运动单元个数,即聚类个数;然后根据聚类个数,利用最小距离分类器对阵列式sEMG波形聚类;最后重新计算新的运动单元发放序列向量,循环执行程序直至分解完成,并对所有发放序列归类整理,优化结果。该分解方法准确性高,计算快速、实现简单。
Description
技术领域
本发明涉及一种阵列式sEMG信号自动分解方法。
背景技术
表面肌电信号( surface EMG, sEMG) 是利用表面电极从人体体表检测肌电信号,与针电极肌电信号(Needle EMG, NEMG)相比,它具有无创性、易于患者接受的特点,因此应用前景广阔。实验表明,利用阵列式sEMG能够提高运动单元(MU)的检测率,特别是提高小幅值运动单元活动电位(MUAP)的检出和识别效果。sEMG的起源是MUAP,活动电位由给定肌肉收缩过程中所激活的每一个运动单位所释放。在任何一个给定的募集模式,众多的运动单位以非同步的模式被激活,这些运动单位活动的总和构成了肌电信号的强度。阵列式sEMG信号实质上是多个运动单位活动电位的总和,其波幅典型地在1~5000uv之间,频率范围为10~400Hz。在临床上,通过阵列式sEMG可以较全面地了解神经肌肉的功能状态,鉴别神经源性和肌源性疾病,判断神经损伤的部位、程度及恢复状况,并且阵列式sEMG信号的检测分析对康复医学及运动医学也有具有重要意义。
阵列式sEMG分解实质上是对sEMG包含的运动单元发放序列进行分类,提取出运动单元各电位序列。目前,sEMG分类方法主要有:K均值聚类算法、模板匹配法、人工神经网络(ANN)算法、实时线性混叠盲信号分离算法、独立成分分折(ICA)、卷积核补偿算法等方法。K-均值聚类算法需要指定聚类的类别数,而在肌电信号缺少运动单元发放的先验知识,难以对类别进行精确的指定。模板匹配法由于模板获取困难,应用受限。ANN可以解决含更多迭加波形情形和在低信噪比时更好地消除绝对误差,然而,ANN方法一旦训练后,网络就固定不变,当模板的形状发生变化时,神经网络还需重新训练,所以强健性不好。ICA是一种盲信号分解技术, 它假设构成肌电信号的各MUAPT 相互独立,然后把信号分解成若干相互独立的成分。卷积核补偿算法方法是一种盲信号分解方法,该方法在波形叠加严重时效果仍然不理想。阵列式sEMG的信噪比较低,MUAP波形的变异性强且相互间的叠加程度较大,这是导致其分解困难的主要原因。很多研究者将插入式肌电信号的分解方法做一些修改应用到表面肌电信号的分解研究中,但是都不能达到插入式肌电信号的分解效果。总的来看,表面肌电信号分解研究还处于探索阶段,是肌电研究领域的难点之一。
发明内容
鉴于上述问题,本发明的目的在于提供一种阵列式sEMG信号自动分解方法。
为实现上述目的,针对阵列式表面电极肌电信号,提出基于最小距离器的聚类方法,对发放波形进行分类,并且在这个过程中,不求解矩阵,通过卷积核补偿方法得到肌肉运动单元的数目和发放序列,实现阵列式sEMG的分解。由于该方法同时考虑发放波形和发放时刻,相对于其它方法,该方法具有sEMG分解精度高的优点。
本发明公开了一种阵列式sEMG信号自动分解方法。其包括以下步骤:
步骤一:阵列式sEMG信号预处理:对sEMG信号滤波,剔除干扰;
步骤二:计算初始发放序列向量:利用阵列sEMG信号各通道信号的相关性,计算初始发放序列向量,作为提取初值;
步骤三:确定运动单元个数:利用时域减法,从sEMG中获得运动单元个数,作为聚类个数;
步骤四:对阵列式sEMG波形聚类:根据所述聚类个数,利用最小距离分类器,对阵列式sEMG波形聚类;
步骤五:计算新的运动单元发放序列向量:根据所述波形聚类结果,取含最多时刻的类,计算该类时刻波形均值,计算新的发放序列向量;
步骤六:对所有发放序列归类整理:重复步骤二-----步骤五,剔除重复的以及不合理的发放序列向量,优化结果。
优化的技术措施包括:
上述sEMG信号的互相关矩阵表示为:
其中是采样时刻,是第个采样时刻的阵列信号,是第个采样时刻的阵列信号转置,是数序期望。
运动单元发放序列表示为:
其中阵列信号互相关矩阵的逆矩阵。
上述欧氏距离公式表示为:
其中是某个时刻对应的阵列sEMG的值,是K个时刻对应的阵列sEMG值的平均值。
新的运动单元发放序列的计算公式表示为:
其中是N个波形平均值转置。
与现有技术相比,本发明的一种阵列式sEMG信号自动分解方法,在计算初始发放序列时,由于最大峰值往往是由于干扰造成的,导致其值畸高,所以本发明采用的是次大峰值所对应的时刻,得到更为准确的发放序列初值;本发明利用时域减法获得相减次数M,该次数M作为分类个数,实际上提供了运动单元数目,为后续的最小距离器分类提供先验知识,提高了分解准确性;本发明利用欧氏距离作为分类依据,实现方便、简单;由于本发明分类针对的是波形,而互相关矩阵计算的是发放时刻,所以分解过程同时考虑了发放时刻和发放波形,大大提高了分解的准确性;本发明不需要计算运动单元发放序列和阵列sEMG信号间混合矩阵,大大减少计算时间,提高了效率。本分解方法运行过程不需要人工干预,使用方便。
附图说明
图1是本发明流程图。
具体实施方式
以下结合附图实例对本发明作进一步详细描述,本领域技术人员可由本说明书所揭露的内容轻易地实现。
如图1所示为本发明的流程图。
步骤一:阵列式sEMG信号预处理:对sEMG信号滤波,剔除干扰。由于sEMG信号中包含各种干扰信号,预处理首先需要采用带通滤波器,保留10Hz--500Hz频段信号,然后采用陷波滤波器,滤除50Hz工频干扰。
步骤二:计算初始发放序列向量:利用阵列sEMG信号各通道信号的相关性,计算初始发放序列向量,作为提取初值。具体过程为:首先计算阵列sEMG信号互相关矩阵以及互相关矩阵逆矩阵,互相关矩阵表示为:
其中是采样时刻,是第个采样时刻的阵列信号,是第个采样时刻的阵列信号转置,是数序期望。计算互相关矩阵的逆矩阵,即
然后采样时刻取sEMG信号能量的中值,能量按照下式计算:
取能量中值所对应的时刻。最后利用如下公式计算得到运动单元发放序列向量的初值
步骤三:确定运动单元个数:利用时域减法,从sEMG中获得运动单元个数,作为聚类个数。具体过程为:首先从运动单元发放序列向量的初值中,找到次大峰值,记该次大峰值所对应时刻为,其次计算时刻发放序列,计算公式如下:
然后发放序列中找到K个最大峰值所对应的时刻,一般取。取K个时刻所对应的阵列sEMG值,相加并除以K,得到K个时刻的均值。最后在时域上,用阵列sEMG信号逐次减K个时刻波形平均值,直到结果为负数停止,记录相减次数M,该M值即运动单元个数。
步骤四:对阵列式sEMG波形聚类:根据所述聚类个数M,利用最小距离分类器,对阵列式sEMG波形聚类。
最小距离分类器中采用欧氏距离,欧氏距离公式表示为:
其中是某个时刻对应的阵列sEMG的值,是K个时刻对应的阵列sEMG值的平均值。
步骤五:计算新的运动单元发放序列向量:根据上述波形聚类结果,取含最多时刻的类,计算该类时刻波形均值,计算新的发放序列向量。
新的运动单元发放序列的计算公式表示为:
其中是N个波形平均值转置。
步骤六:对所有发放序列归类整理:重复步骤二--步骤五,直到初始发放时刻不能提取,sEMG提取完成。剔除重复的以及不合理的发放序列向量,优化结果。不合理的发放序列指发放时刻小于15毫秒的序列,需剔除。
Claims (3)
1.一种阵列式sEMG信号自动分解方法,其特征是包括以下步骤:
步骤一:阵列式sEMG信号预处理:对sEMG信号滤波,剔除干扰;
步骤二:计算初始发放序列向量:利用阵列sEMG信号各通道信号的相关性,计算次大峰值的初始发放序列向量,作为提取初值;
步骤三:确定运动单元个数:利用时域减法,从sEMG中减去K个最大峰值所对应波形的平均值,获得运动单元个数,作为聚类个数;
步骤四:对阵列式sEMG波形聚类:根据所述聚类个数,利用最小距离分类器,对阵列式sEMG波形聚类;
步骤五:计算新的运动单元发放序列向量:根据所述波形聚类结果,取含最多时刻的类,计算该类时刻波形均值,计算新的发放序列向量;
步骤六:对所有发放序列归类整理:重复步骤二-----步骤五,剔除重复以及不合理的发放序列向量,优化结果。
2.根据权利1要求所述的一种阵列式sEMG信号自动分解方法,其特征是:利用K个时刻波形平均值,通过时域减法确定聚类个数。
3.根据权利2要求所述的一种阵列式sEMG信号自动分解方法,其特征是:在最小距离分类器中,采用欧氏距离,以sEMG波形作为分类目标,提高分解可靠性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510512579.4A CN104997508B (zh) | 2015-08-19 | 2015-08-19 | 阵列式sEMG信号自动分解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510512579.4A CN104997508B (zh) | 2015-08-19 | 2015-08-19 | 阵列式sEMG信号自动分解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104997508A CN104997508A (zh) | 2015-10-28 |
CN104997508B true CN104997508B (zh) | 2018-05-08 |
Family
ID=54370578
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510512579.4A Expired - Fee Related CN104997508B (zh) | 2015-08-19 | 2015-08-19 | 阵列式sEMG信号自动分解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104997508B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105326501B (zh) * | 2015-12-10 | 2021-02-12 | 宁波工程学院 | 基于多通道sEMG的肌肉状态评估方法 |
CN105342610B (zh) * | 2015-12-10 | 2017-12-22 | 宁波工程学院 | 肌肉运动单元深度估计快速方法 |
CN105662336B (zh) * | 2015-12-23 | 2019-03-19 | 黑龙江科技大学 | 一种信号去噪处理方法及装置 |
CN105975917B (zh) * | 2016-04-28 | 2019-06-21 | 宁波工程学院 | 面向强干扰的阵列式表面肌电信号分解方法 |
CN107526952B (zh) * | 2016-06-22 | 2020-09-01 | 宁波工程学院 | 基于多通道表面肌电信号的身份识别方法 |
CN106264517B (zh) * | 2016-09-30 | 2019-05-14 | 浙江大学 | 一种选择心电测量位置的方法及系统 |
CN108403114B (zh) * | 2018-02-11 | 2021-02-02 | 宁波工程学院 | 一种面向恒力的阵列式表面肌电信号分解方法 |
CN108378848A (zh) * | 2018-02-11 | 2018-08-10 | 宁波工程学院 | 肌肉运动单元空间位置估计方法 |
CN110720912B (zh) * | 2019-10-12 | 2022-04-26 | 宁波工程学院 | 一种基于波形相关性的肌肉运动单元提取方法 |
CN110720910B (zh) * | 2019-10-12 | 2022-03-29 | 宁波工程学院 | 一种基于相关性的肌肉运动单元搜索方法 |
CN110720911B (zh) * | 2019-10-12 | 2022-08-30 | 宁波工程学院 | 一种肌肉运动单元提取方法 |
CN111419230A (zh) * | 2020-04-17 | 2020-07-17 | 上海交通大学 | 一种用于运动单元解码的表面肌电信号采集系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101488189B (zh) * | 2009-02-04 | 2012-01-18 | 天津大学 | 基于独立分量自动聚类处理的脑电信号处理方法 |
CN102961203B (zh) * | 2012-12-10 | 2015-04-22 | 杭州电子科技大学 | 基于emd样本熵的表面肌电信号识别方法 |
-
2015
- 2015-08-19 CN CN201510512579.4A patent/CN104997508B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN104997508A (zh) | 2015-10-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104997508B (zh) | 阵列式sEMG信号自动分解方法 | |
Negro et al. | Multi-channel intramuscular and surface EMG decomposition by convolutive blind source separation | |
Nazarpour et al. | Application of higher order statistics to surface electromyogram signal classification | |
CN105997064B (zh) | 一种用于人体下肢表面肌电信号的辨识方法 | |
CN110353702A (zh) | 一种基于浅层卷积神经网络的情感识别方法及系统 | |
CN107260166A (zh) | 一种实用化在线脑电伪迹剔除方法 | |
Knight | Signal fraction analysis and artifact removal in EEG | |
CN110059564B (zh) | 基于功率谱密度和互相关熵谱密度融合的特征提取方法 | |
CN102961203A (zh) | 基于emd样本熵的表面肌电信号识别方法 | |
CN112370017B (zh) | 脑电信号分类模型的训练方法和装置、以及电子设备 | |
CN109009098B (zh) | 一种运动想象状态下的脑电信号特征识别方法 | |
Phinyomark et al. | EMG amplitude estimators based on probability distribution for muscle–computer interface | |
AYDEMİR | Common spatial pattern-based feature extraction from the best time segment of BCI data | |
CN105956547B (zh) | 基于阵列式表面肌电信号平滑的分解方法 | |
CN109620260A (zh) | 心理状态识别方法、设备及存储介质 | |
CN113378737A (zh) | 植入式脑机接口神经元锋电位分类方法 | |
Rong et al. | Classification of surface EMGs using wavelet packet energy analysis and a genetic algorithm-based support vector machine | |
CN105975917B (zh) | 面向强干扰的阵列式表面肌电信号分解方法 | |
Naik et al. | Evaluation of higher order statistics parameters for multi channel sEMG using different force levels | |
Vigneshwari et al. | Analysis of finger movements using EEG signal | |
Assi et al. | Kmeans-ICA based automatic method for ocular artifacts removal in a motorimagery classification | |
KR20140110304A (ko) | 상관 특징 분석 기반 실시간 emg 패턴 인식 방법 | |
Putra et al. | Classification of EMG during walking using principal component analysis and learning vector quantization for biometrics study | |
Ahmed et al. | Control of prosthetic hand by using mechanomyography signals based on support-vector machine classifier | |
Fu et al. | A detection scheme for frontalis and temporalis muscle EMG contamination of EEG data |
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: 20180508 Termination date: 20200819 |
|
CF01 | Termination of patent right due to non-payment of annual fee |