CN104490390B - 基于神经电生理信号联合分析的人体运动能力判别方法 - Google Patents
基于神经电生理信号联合分析的人体运动能力判别方法 Download PDFInfo
- Publication number
- CN104490390B CN104490390B CN201410841357.2A CN201410841357A CN104490390B CN 104490390 B CN104490390 B CN 104490390B CN 201410841357 A CN201410841357 A CN 201410841357A CN 104490390 B CN104490390 B CN 104490390B
- Authority
- CN
- China
- Prior art keywords
- analysis
- hrv
- signal
- emg
- frequency
- 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
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000007103 stamina Effects 0.000 title abstract 3
- 230000007831 electrophysiology Effects 0.000 title abstract 2
- 238000002001 electrophysiology Methods 0.000 title abstract 2
- 230000033001 locomotion Effects 0.000 claims abstract description 43
- 239000013598 vector Substances 0.000 claims abstract description 28
- 238000000513 principal component analysis Methods 0.000 claims abstract description 21
- 230000000694 effects Effects 0.000 claims abstract description 14
- 238000003909 pattern recognition Methods 0.000 claims abstract description 10
- 238000002567 electromyography Methods 0.000 claims description 33
- 230000003183 myoelectrical effect Effects 0.000 claims description 12
- 238000001228 spectrum Methods 0.000 claims description 10
- 238000012850 discrimination method Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 210000005036 nerve Anatomy 0.000 claims description 6
- 230000000747 cardiac effect Effects 0.000 abstract description 4
- 238000010219 correlation analysis Methods 0.000 abstract description 3
- 230000010354 integration Effects 0.000 abstract description 2
- 230000004069 differentiation Effects 0.000 abstract 1
- 230000001360 synchronised effect Effects 0.000 abstract 1
- 239000011159 matrix material Substances 0.000 description 14
- 210000003205 muscle Anatomy 0.000 description 12
- 230000008859 change Effects 0.000 description 8
- 238000002565 electrocardiography Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 238000012706 support-vector machine Methods 0.000 description 7
- 230000003595 spectral effect Effects 0.000 description 6
- 210000003403 autonomic nervous system Anatomy 0.000 description 5
- 230000001186 cumulative effect Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 230000002889 sympathetic effect Effects 0.000 description 5
- 238000012549 training Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 230000009084 cardiovascular function Effects 0.000 description 4
- 230000008602 contraction Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 210000000748 cardiovascular system Anatomy 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000005486 microgravity Effects 0.000 description 3
- 206010001497 Agitation Diseases 0.000 description 2
- 206010028289 Muscle atrophy Diseases 0.000 description 2
- 230000002411 adverse Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000016507 interphase Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000007659 motor function Effects 0.000 description 2
- 230000020763 muscle atrophy Effects 0.000 description 2
- 201000000585 muscular atrophy Diseases 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000001515 vagal effect Effects 0.000 description 2
- 210000001186 vagus nerve Anatomy 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- 238000012356 Product development Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 208000007502 anemia Diseases 0.000 description 1
- 210000000467 autonomic pathway Anatomy 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 210000001124 body fluid Anatomy 0.000 description 1
- 239000010839 body fluid Substances 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000004706 cardiovascular dysfunction Effects 0.000 description 1
- 210000003169 central nervous system Anatomy 0.000 description 1
- 108091008690 chemoreceptors Proteins 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 208000035475 disorder Diseases 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 239000003792 electrolyte Substances 0.000 description 1
- 230000037080 exercise endurance Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000005021 gait Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003387 muscular Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000007427 paired t-test Methods 0.000 description 1
- 210000005037 parasympathetic nerve Anatomy 0.000 description 1
- 210000001002 parasympathetic nervous system Anatomy 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000036454 renin-angiotensin system Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 230000001020 rhythmical effect Effects 0.000 description 1
- 210000001013 sinoatrial node Anatomy 0.000 description 1
- 210000002027 skeletal muscle Anatomy 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 210000002820 sympathetic nervous system Anatomy 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 210000001364 upper extremity Anatomy 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02405—Determining heart rate variability
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
-
- 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]
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Cardiology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Physiology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种基于神经电生理信号联合分析的人体运动能力判别方法,包括以下步骤:同步采集EMG、ECG信号;对ECG信号进行R波峰值识别得到HRV信号,再分别对各信号进行时域分析、频域分析以及非线性分析;对肌电参数与心率变异性参数进行联合分析,采用相关分析寻找具有显著性的关联参数对,研究运动能力和心脏负荷情况的关联效应;提取EMG特征、HRV特征,以及计算EMG信号积分肌电值、均方根分别与HRV参数的相关系数;利用PCA主成分分析去除冗余信息,得到特征向量;将特征向量作为模式识别的输入,从而进行运动能力的判别。本方法可有效地提高人体运动能力判别模型的准确性和简便性,可为航天员的选拔提供依据。
Description
技术领域
本发明涉及人体运动能力判别领域,尤其涉及一种基于神经电生理信号联合分析的人体运动能力判别方法。
背景技术
人是在地球1G重力环境中进化和生活的,因此人体各生理系统的形态、结构和功能都产生了适应1G条件的变化。当航天员进入太空后,长时间的生活在失重环境中(准确的说是处于微重力状态),流体静压消失,体液头向分布会引发机体各生理系统出现一系列的适应性变化,引起人体生理系统的明显改变,例如:水和电解质代谢紊乱、心血管功能失调、航天贫血症等,尤其是对运动功能的影响,如肌萎缩、肌肉耐力的下降。从而导致航天员工作绩效的下降,进而影响整个航天任务的顺利进行。除此之外,还可能影响到航天员从空间返回到地面的再适应过程,其主要表现是立位耐力不良和运动耐力下降。因此,现阶段急需建立准确的人体运动能力变化的判别模型,提出一套能够科学的评判人体运动能力的指标,进而应用于航天员的选拔工作中,以及通过对人体运动能力分析评估,找到能够降低失重环境对于航天员健康以及工作绩效影响的对抗方法。
人体运动功能主要表现在行走的步态和上肢的典型动作中。描述这些动作特征的信息主要有运动学参数、力学参数和生物电信号三大类,从而进行功能评定。人体运动的动力来自骨骼肌的协同收缩,肌电信号是肌肉自主收缩或电诱发收缩过程中产生的生理电变化,其信息可以用来描述运动功能特征和分析运动的质量。
心率变异性是指人体心脏搏动周期存在的微小变异,通常情况下是指逐次心动周期间体表心电图R波间期的微小涨落。近年来研究表明,脑的高级神经活动、中枢神经系统的自发性节律活动、呼吸以及压力、化学感受器活动等因素对心脏的自主神经系统进行调制结果导致心率变异性(Heart rate variability,HRV)的产生。换言之,HRV信号蕴涵了有关心血管系统神经及体液调节的大量信息,因而可以通过检测测试者的心率数据,得到许多有关自主神经系统平衡状态的信息。心率变异性是衡量自主神经系统功能状态的一个重要指标,适当的体育运动可以增强自主神经系统的调节能力,HRV分析表现出了潜在的令人注目的应用价值。
现阶段科学研究对于肌电和心率变异性的研究都分别取得了较为瞩目的成就,但研究者通常只是单方面的通过肌肉以及肌肉电信号的特征来对运动能力进行分析。由于航天员长期处于失重条件这种工作环境,导致肌肉的长时间废用,不仅会引起肌肉萎缩,还直接影响到肌肉的工作能力和耐力,而且肌肉泵作用降低也可造成心血管功能下降。由此,从心血管功能变化规律入手,与肌电特征联合分析,可以为建立更为全面的运动能力变化判别模型与评价方法提供新的思路。
发明内容
本发明提供了一种基于神经电生理信号联合分析的人体运动能力判别方法,本发明通过寻找具有显著性的关联参数对,研究运动能力和心脏负荷情况的关联效应;利用主成分分析爱(Principle component analysis,PCA)去除信息冗余后得到特征向量,作为后续模式识别的输入,从而进行运动能力判别,详见下文描述:
一种基于神经电生理信号联合分析的人体运动能力判别方法,所述人体运动能力判别方法包括以下步骤:
同步采集EMG、ECG信号;对ECG信号进行R波峰值识别得到HRV信号,再分别对各信号进行时域分析、频域分析以及非线性分析;
对肌电参数与心率变异性参数进行联合分析,采用相关分析寻找具有显著性的关联参数对,研究运动能力和心脏负荷情况的关联效应;
提取EMG特征、HRV特征,以及计算EMG信号积分肌电值、均方根分别与HRV参数的相关系数;
利用PCA主成分分析去除冗余信息,得到特征向量;
将特征向量作为模式识别的输入,从而进行运动能力的判别。
所述提取EMG特征、HRV特征,以及计算EMG信号积分肌电值、均方根分别与HRV参数的相关系数具体为:
1)EMG特征提取:通过对肌电信号进行时频域分析得到积分肌电值、均方根、平均功率频率、中值频率四个特征参数,分别记为特征向量P1,P2,P3,P4;
2)HRV特征提取:计算总功率、极低频、低频和高频的频谱熵,分别记为特征向量P5,P6,P7,P8;
3)EMG信号积分肌电值、均方根分别与HRV参数的相关系数,分别记为特征向量P9,P10,P11,P12,P13,P14,P15,P16。
所述利用PCA主成分分析去除冗余信息,得到特征向量的步骤具体为:
对原始数据Pn*16进行标准化处理,得到矩阵Pn*16 *;然后求其协方差矩阵C16*16;
对协方差矩阵进行特征根分解,得到特征根∧及特征向量U,其中特征向量作为主成分的坐标轴,构成新的矢量空间,每个特征根的大小代表每个主成分蕴含的信息量;
求原始数据在新坐标系中的投影F;
求累积贡献率;每个主成分的特征根大小代表其蕴含信息量的多少;
选定合适的累积贡献率,使前d个主成分Fn*d作为新的数据,替代原始数据Pn*16进行模式识别。
本发明提供的技术方案的有益效果是:本发明解决了人体运动能力划分测量复杂,准确率低的问题,提出利用EMG信号特征参数与HRV频域特征参数联合分析的方法,从肌肉系统与心血管系统两大方面入手,更为准确全面且针对性更强。与利用单一信号相比,信息来源更丰富。并且利用PCA主成分分析去除冗余信息,保留了对人体运动能力划分更为敏感的参数特征,可以更准确,快速地进行划分,可有效地提高人体运动能力判别模型的准确性和简便性。并且作为全面严谨的运动能力判别模型可为航天员的选拔提供依据。并根据其运动能力情况在飞行前制定针对性的防护训练任务,提高航天员在外太空的适应能力;设计对抗运动,减缓飞行中微重力的不利影响,提升航天任务的有效执行性。
附图说明
图1为一种基于神经电生理信号联合分析的人体运动能力判别方法的流程图;
图2为心电信号R波识别示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面对本发明实施方式作进一步地详细描述。
参见图1,该基于多生理参数联合分析的人体运动能力判别方法的流程为:首先同步采集运动状态下的肌电信号(Electromyography,EMG)、心电信号(Electrocardio,ECG),然后对ECG信号进行R波识别,如图2“*”处为识别到的R波,进而得到心率变异性(Heartrate variability,HRV)信号,再分别对各信号进行时域分析、频域分析以及非线性分析,并分别计算EMG特征参数与HRV特征参数的相关系数,寻找显著关联参数对,然后利用PCA主成分分析去除冗余信息,得到特征向量作为后续模式识别的输入,从而进行运动能力的判别,下面对各个步骤进行详细描述:
101:同步采集EMG、ECG信号;对ECG信号进行R波峰值识别得到HRV信号,再分别对各信号进行时域分析、频域分析以及非线性分析;
肌电信号包含大量的人体运动信息,经过分析处理后,可对人体运动能力进行分析。心率变异即逐次窦性心跳间期之间的微小涨落,直接受自主神经系统中的交感神经和副交感神经的双重支配。肌肉泵作用的变化会影响心血管功能,而心率变异性能很好的反应心血管功能变化。
在本方法中,由于该两个参数测量方法成熟,与其他生理参数相比能够分别从心血管系统与肌肉系统两大方面更准确、敏感、全面的反应人体运动状况,所以采用这两种生理参数有效融合来进行人体运动能力分析在实际应用中具有优越性。
提取分析的目标肌肉肌电信号参数包括时域参数、频域参数与非线性参数。
时域参数包括积分肌电值、均方根值。频域参数包括平均功率频率、中值频率等。
积分肌电值(Integral electromyography,IEMG)
IEMG是指所得肌电信号经过全波整流后单位时间内曲线下面积的综合,它可以反映肌电信号随时间进行的强弱变化。肌电积分用于分析肌肉在单位时间内的收缩特性。
均方根(root mean square,RMS)
RMS表示肌电信号振幅大小的指标,即可评价信号的强度。
频谱参数包括中值频率、功率谱密度、平均功率频率等;
平均功率频率(mean power frequency,MPF)
其中,S(f)是信号功率谱;f为频率。
中值频率(median frequency,MF)
中值频率的定义将频谱图均分为两部分相等面积的频率,且此功率左右的频谱面积相等。
非线性参数包括:复杂度、李雅普诺夫指数、关联维度、模糊熵等。
统计分析过程包括对不同运动能力状态下上述参数的对比分析,采用的分析思路为采用配对T检验分析其在不同运动能力状态下的显著性变化及采用单因素方差分析研究时间因素对于参数的影响。
心率变异性参数包括时域波形参数、频谱参数和非线性参数。其中时域波形参数包括心搏间期均值、心搏间期标准差、心搏间期差值的均方根值、相邻RR间期差值>50ms心搏数占总RR间期(相邻两个R波间期)数百分比、心率变异系数等。
正常窦性心搏间期的均值(meanRR)
反应心率变异性的整体变化;RRi为第i个RR间期。
正常窦性心搏间期的标准差(SDNN)
反应交感神经核迷走神经的总张力,一般用来描述RR间期的分散特性。
相邻正常窦性心搏间期差值的均方根值(RMSSD):
反映心率变速变化的程度;RRi+1为第i+1个RR间期。
相邻RR间期差值>50ms的心搏数占总RR间期数的百分比(PNN50):反映了心律间期突然改变的情况。
心率变异系数(CV):是测定正常心搏间期变化的大小及快慢,反应窦房结自律性受自主神经(迷走神经、交感神经)系统的调节作用。
心率变异性的时域分析虽然方便,但该分析的灵敏性较差,无法准确地反映出交感神经系统和副交感神经系统的活动情况,而频域分析在这方面具有优势。
频域参数包括RR间期总功率谱、极低频功率、低频功率、高频功率、低频和高频功率比值等。
频域分析就是将RR间期的变化转为频谱计算功率谱密度(Power SpectralDensity,PSD),即将心率变化信号分解为不同频率成分,并将其相对强度定量为功率,提供各种频率成分的功率谱测定。
人类的心率变异性功率谱范围一般是0-0.4Hz,频谱成分和频段划分如下:
TF:总功率,频率范围≤0.4Hz;
VLF:极低频功率,频率范围≤0.04Hz。产生机理尚不清楚,可能与血容量状态、心泵功能、肾素-血管紧张素系统的活动水平等因素的变化对心率的影响有关。
LF:低频功率,频率范围0.04Hz-0.15Hz,主要反映交感神经兴奋性。
HF:高频功率,频率范围0.15Hz-0.4Hz,主要反映迷走神经兴奋性。
LH:低频功率和高频功率比值,可以反映交感神经和迷走神经活动的均衡性。
非线性参数包括RR间期的李亚普诺夫指数、排列熵、近似熵等。
102:对肌电参数与心率变异性参数进行联合分析,采用相关分析寻找具有显著性的关联参数对,研究运动能力和心脏负荷情况的关联效应。
相关分析
当一个或几个相互联系的变量取一定数值时,与之相对应的另一变量的值虽然不确定,但它仍按某种规律在一定的范围内变化。变量间的这种相互关系,成为具有不确定性的相关关系。它反映现象之间客观存在的、不严格、不确定的相互依存关系。这种关系不能在个别现象体现其关系的规律性,必须在大量现象中才能体现出来。
相关分析是用一个指标来表明现象间相互依存关系的密切程度。
Pearson简单相关系数用来度量正态分布的定距变量间的线性相关关系,其相关系数定义如下:
其中,x和y分别代表两组数据,长度都为n,和分别是他们的平均值。
Pearson相关系数假设检验t值计算公式:
式中r是相关系数,n是样本观测量数,n-2是自由度。零假设为总体中两个变量间的相关系数为0,当t>t0.05(n-2)时,p<0.05拒绝原假设,说明两个变量间有显著关系。
例如,在模拟失重情况下要求8名被试完成推舱门任务,分别同步记录模拟失重前1天、模拟失重后7天、模拟失重后14天以及模拟失重后28天完成任务时候的右臂肱二头肌的EMG信号以及ECG信号,分别计算不同时间点下EMG信号特征参数(如IEMG)与由ECG信号得到的HRV信号的特征参数(如LF)之间的相关性,得到四个相关系数,若相关系数满足上述假设检验条件,则说明IEMG与LF存在显著关系,为显著性关联参数对,并可继续研究这种显著关系随着模拟失重时间的变化规律。
103:提取EMG特征、HRV特征,以及计算EMG信号积分肌电值、均方根分别与HRV参数的相关系数;
1)EMG特征提取:通过对肌电信号进行时频域分析得到积分肌电值、均方根、平均功率频率、中值频率四个特征参数,分别记为特征向量P1,P2,P3,P4。
2)HRV特征提取:计算总功率(TF)、极低频(VLF)、低频(LF)和高频(HF)的频谱熵,分别记为特征向量P5,P6,P7,P8。
3)EMG信号积分肌电值、均方根分别与HRV参数的相关系数,分别记为特征向量P9,P10,P11,P12,P13,P14,P15,P16。
104:利用PCA主成分分析去除冗余信息,得到特征向量;
各个参数所含的信息之间具有一定的重叠性和相关性,若直接将它们用于模式识别,会造成模型参数的过度拟合而降低分类的准确性和可靠性,且会因为数据量过大而降低分类的速度。因此在模式分类之前,本发明首先利用PCA对得到的特征矩阵进行降维处理。
PCA根据方差最大化原理,用一组线性无关且相互正交的新向量表征原来的数据矩阵的行(或列),达到压缩变量个数,剔除冗余信息,最大化保存有效信息的目的。原始向量组为(P1,P2,…,P16),主成分记为(F1,F2,…,Fm),通常m小于16。则主成分与原始向量组的关系为:
其中,F1蕴含信息量最多,具有最大方差,F2,…,Fm依次递减。因此主成分分析的过程可以看作是确定权重系数的过程。令P1,P2,…,P16的相关矩阵为R,(ai1、ai2…ai16)则是相关矩阵R的第i个特征向量。
在本发明中,首先对该16个变量进行n次观测,得到的观测数据可用下面的矩阵表示:
其中,P1i,P2i,…,Pni(i=1,2…,16)为向量PI(I=1,2…,16)里面所的单个元素。
用PCA进行特征降维的求解过程如下:
1)对原始数据Pn*16进行标准化处理,矩阵中的元素减去所在列的均值,然后除以所在列的标准差,使得每个变量的均值为0,方差为1,得到矩阵Pn*16 *。
Pn*16 *=[yij]n*16,i=1,2,…,n;j=1,2,…,16. (12)
其中,
2)然后求其协方差矩阵C16*16,Pn*16 *中任两列之间可以计算两变量间的协方差:
其中,n表示样本数量,x,y分别表示某一个特征。于是得到协方差矩阵:
3)对协方差矩阵C进行特征根分解,得到特征根∧及特征向量U,其中特征向量作为主成分的坐标轴,构成新的矢量空间,每个特征根的大小代表每个主成分蕴含的信息量。
C16=U16*16Λ16*16U16*16' (16)
4)求原始数据P在新坐标系中的投影F:
Fn*16=Pn*16U16*16 (17)
5)求累积贡献率。每个主成分的特征根大小代表其蕴含信息量的多少。求前k个主成分的累积贡献率。
其中,λi是求出的第i个特征根。
6)选定合适的累积贡献率,使前d个主成分Fn*d作为新的数据,替代原始数据Pn*16进行模式识别。
105:将特征向量作为模式识别的输入,从而进行运动能力的判别。
支持向量机(Support Vector Machine,SVM)是Cortes和Vapnik于1995年首先提出的,它在解决小样本、非线性及高维模式识别中表现出许多特有的优势。支持向量机方法是建立在统计学习理论的VC维理论和结构风险最小原理基础上的。
由于样本数据集偏小,特征提取后,使用支持向量机(Support Vector Machine,SVM)分类器对特征进行识别,进行运动能力划分。
用支持向量机做模式识别时,将由PCA去除冗余信息后的特征参数作为训练支持向量机的输入参数,并以运动能力级别为输出,经过训练得到基于EMG、HRV特征参数以及EMG与HRV两者特征的相关系数的运动能力判别模型,从而达到运动能力划分的目的。
本发明的主旨是提出一种基于EMG,HRV以及两者特征的相关系数的多生理参数联合分析的人体运动能力判别方法,通过计算EMG与HRV参数的相关系数,并与所提取的EMG,HRV信号的特征参数融合作为特征矩阵,继而利用PCA主成分分析去除冗余信息,从而准确、客观、简便的进行人体运动能力判别。该项发明可有效地提高人体运动能力判别系统准确性和简便性,加深对于运动能力和心脏负荷情况的关联效应的研究。进而在航天领域的飞行训练、航天控制、着陆返回等阶段制定相应的防护训练任务,提高航天员在外太空的适应能力,减缓飞行中微重力的不利影响,提升航天任务的有效执行性;获得可观的社会效益和经济效益。最佳实施方案拟采用专利转让、技术合作或产品开发。由于该技术操作简单,敏感性强,基于该技术开发的产品还可应用于人体运动能力监测、运动科研等多种场景等。
本领域技术人员可以理解附图只是一个优选实施例的示意图,上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种基于神经电生理信号联合分析的人体运动能力判别方法,其特征在于,所述人体运动能力判别方法包括以下步骤:
同步采集EMG、ECG信号;对ECG信号进行R波峰值识别得到HRV信号,再分别对EMG信号、HRV信号进行时域分析、频域分析以及非线性分析;
对肌电参数与心率变异性参数进行联合分析,采用相关分析寻找具有显著性的关联参数对,研究运动能力和心脏负荷情况的关联效应;
提取EMG特征、HRV特征,以及计算EMG信号积分肌电值、均方根分别与HRV特征的相关系数;
利用PCA主成分分析去除冗余信息,得到特征向量;
将特征向量作为模式识别的输入,从而进行运动能力的判别。
2.根据权利要求1所述的一种基于神经电生理信号联合分析的人体运动能力判别方法,其特征在于,所述提取EMG特征、HRV特征,以及计算EMG信号积分肌电值、均方根分别与HRV特征的相关系数具体为:
1)EMG特征提取:通过对肌电信号进行时频域分析得到积分肌电值、均方根、平均功率频率、中值频率四个特征参数,分别记为特征向量P1,P2,P3,P4;
2)HRV特征提取:计算总功率、极低频、低频和高频的频谱熵,分别记为特征向量P5,P6,P7,P8;
3)EMG信号积分肌电值、均方根分别与HRV特征的相关系数,分别记为特征向量P9,P10,P11,P12,P13,P14,P15,P16。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410841357.2A CN104490390B (zh) | 2014-12-30 | 2014-12-30 | 基于神经电生理信号联合分析的人体运动能力判别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410841357.2A CN104490390B (zh) | 2014-12-30 | 2014-12-30 | 基于神经电生理信号联合分析的人体运动能力判别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104490390A CN104490390A (zh) | 2015-04-08 |
CN104490390B true CN104490390B (zh) | 2017-03-08 |
Family
ID=52931918
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410841357.2A Active CN104490390B (zh) | 2014-12-30 | 2014-12-30 | 基于神经电生理信号联合分析的人体运动能力判别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104490390B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113499060B (zh) * | 2015-10-07 | 2024-08-27 | 普莱柯迪尔公司 | 用于产生指示心脏状况的信息的方法和设备 |
CN105868532B (zh) * | 2016-03-22 | 2019-03-05 | 曾金生 | 一种智能评估心脏衰老程度的方法及系统 |
CN106126959A (zh) * | 2016-07-19 | 2016-11-16 | 北京心量科技有限公司 | 一种运动风险分析方法和装置 |
CN106228005B (zh) * | 2016-07-19 | 2019-07-12 | 北京心量科技有限公司 | 数据处理方法及系统 |
CN106202953B (zh) * | 2016-07-19 | 2019-07-02 | 北京心量科技有限公司 | 一种数据分析方法及系统 |
CN107808694B (zh) * | 2016-09-08 | 2023-04-14 | 深圳先进技术研究院 | 基于gmfcs评判脑瘫患儿行走能力的系统及方法 |
CN107145236B (zh) * | 2017-05-12 | 2020-02-07 | 中国科学技术大学 | 一种基于腕部肌腱压力相关特性的手势识别方法及系统 |
CN108926326B (zh) * | 2017-05-22 | 2020-06-26 | 北京大学 | 一种用于评价人体随运动平台的动态平衡控制能力的方法 |
WO2019016759A1 (en) * | 2017-07-19 | 2019-01-24 | Bloom Technologies NV | MONITORING UTERINE ACTIVITY AND ASSESSING A RISK OF BIRTH BEFORE TERM |
TWI728401B (zh) * | 2018-09-25 | 2021-05-21 | 高雄榮民總醫院 | 以運動測試評估自主神經功能障礙的電腦程式產品及電腦可讀取媒體 |
CN113194811A (zh) * | 2018-12-29 | 2021-07-30 | 深圳迈瑞生物医疗电子股份有限公司 | 院内病人恢复状态的评估方法、装置、系统及存储介质 |
CN109998522A (zh) * | 2019-02-27 | 2019-07-12 | 东南大学 | 一种基于穿戴式动态心电的运动负荷检测装置及方法 |
CN110200617A (zh) * | 2019-06-26 | 2019-09-06 | 慧言科技(天津)有限公司 | 一种利用心率变异性定量评价疼痛的系统 |
CN110490257A (zh) * | 2019-08-21 | 2019-11-22 | 四川长虹电器股份有限公司 | 一种基于去趋势项的生理电信号熵分析方法 |
CN110974212A (zh) * | 2019-11-07 | 2020-04-10 | 曲阜师范大学 | 一种融合心电和肌电特征的康复训练运动状态监测方法及系统 |
CN116978561B (zh) * | 2023-07-17 | 2024-03-22 | 北京师范大学-香港浸会大学联合国际学院 | 一种基于模糊熵的运动风险评估方法、系统、设备及介质 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3551924B2 (ja) * | 1998-02-26 | 2004-08-11 | オムロンヘルスケア株式会社 | 運動機器 |
US6678549B2 (en) * | 2001-03-26 | 2004-01-13 | Cusimano Maryrose | Combined physiological monitoring system |
AU2004290588A1 (en) * | 2003-11-18 | 2005-06-02 | Vivometrics, Inc. | Method and system for processing data from ambulatory physiological monitoring |
CN101785659B (zh) * | 2009-10-16 | 2012-05-16 | 中国科学院合肥物质科学研究院 | 一种便携式通用生物信号获取和分析系统及方法 |
CN102274108B (zh) * | 2011-05-17 | 2014-01-15 | 天津大学 | 步行器助行中稳定性与脑电-肌电相干系数模糊融合方法 |
CN103405228B (zh) * | 2013-08-21 | 2015-12-23 | 中国科学院自动化研究所 | 一种便携式心电及表面肌电测量装置 |
CN103479362A (zh) * | 2013-09-05 | 2014-01-01 | 南京邮电大学 | 一种基于无线体域网的运动功能康复监测系统 |
-
2014
- 2014-12-30 CN CN201410841357.2A patent/CN104490390B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104490390A (zh) | 2015-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104490390B (zh) | 基于神经电生理信号联合分析的人体运动能力判别方法 | |
Nayak et al. | A review on the nonlinear dynamical system analysis of electrocardiogram signal | |
CN102138789B (zh) | 一种动态心电和运动记录与分析系统 | |
WO2019100566A1 (zh) | 基于人工智能自学习的静态心电图分析方法和装置 | |
Vidaurre et al. | Study of on-line adaptive discriminant analysis for EEG-based brain computer interfaces | |
CN103038772A (zh) | 预测患者的存活性的方法 | |
Yayık et al. | Regularized HessELM and inclined entropy measurement for congestive heart failure prediction | |
CN113951903B (zh) | 基于脑电数据测定的高速铁路调度员超负荷状态识别方法 | |
Tiwari et al. | A comparative study of stress and anxiety estimation in ecological settings using a smart-shirt and a smart-bracelet | |
Guan et al. | Sports fatigue detection based on deep learning | |
Kelwade et al. | Prediction of heart abnormalities using particle swarm optimization in radial basis function neural network | |
Hussain et al. | Classification of normal and pathological heart signal variability using machine learning techniques | |
Fatimah et al. | Automatic detection of mental arithmetic task and its difficulty level using EEG signals | |
Padhy et al. | The power of tensor-based approaches in cardiac applications | |
Ahmad et al. | Multi-level stress assessment from ecg in a virtual reality environment using multimodal fusion | |
Raj et al. | Development of handheld cardiac event monitoring system | |
Xi et al. | Pilot workload prediction from ECG using deep convolutional neural networks | |
Ouelli et al. | AR modeling for automatic cardiac arrhythmia diagnosis using QDF based algorithm | |
Zhou et al. | A novel cardiac arrhythmias detection approach for real-time ambulatory ECG diagnosis | |
Gusarova et al. | Analysis of normal and pathological heart rate variability based on electrocardiogram data | |
Saxena et al. | Extraction of various features of ECG signal | |
Batista et al. | Spectral and time domain parameters for the classification of atrial fibrillation | |
Angrisani et al. | Metrological performance of a single-channel Brain-Computer Interface based on Motor Imagery | |
Qi et al. | Mutual information entropy research on dementia EEG signals | |
Rouhani et al. | Neural Networks based Diagnosis of heart arrhythmias using chaotic and nonlinear features of HRV signals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |