基于采集脑电信号的人体昏厥预警方法.pdf

上传人:sha****007 文档编号:5550006 上传时间:2019-02-11 格式:PDF 页数:12 大小:1.75MB
返回 下载 相关 举报
摘要
申请专利号:

CN201410058482.6

申请日:

2014.02.20

公开号:

CN103845052A

公开日:

2014.06.11

当前法律状态:

授权

有效性:

有权

法律详情:

授权|||实质审查的生效IPC(主分类):A61B 5/0476申请日:20140220|||公开

IPC分类号:

A61B5/0476

主分类号:

A61B5/0476

申请人:

清华大学; 中国人民解放军空军航空医学研究所

发明人:

张涛; 李毅峰; 邓略; 陈勇胜

地址:

100084 北京市海淀区100084-82信箱

优先权:

专利代理机构:

北京清亦华知识产权代理事务所(普通合伙) 11201

代理人:

张大威

PDF下载: PDF下载
内容摘要

本发明提出一种基于采集脑电信号的人体昏厥预警方法。该方法包括:通过脑电图EEG采集仪采集不同G值下载人离心机的M个通道的脑电信号,其中,M为正整数;对M个通道中的脑电信号进行预处理,以获取M个通道中每个通道的低频脑电数据;根据低频脑电数据获取每个脑电信号的频段在每个数据段的参数相关率当参数相关率符合昏厥前预警条件时,进行人体昏厥预警提醒。本发明实施例的方法,可对人体昏厥前的状态进行提前识别和预警,有利于全面了解受训者的体征,指导人体训练具有实际应用价值。

权利要求书

权利要求书
1.  一种基于采集脑电信号的人体昏厥预警方法,其特征在于,包括:
通过脑电图EEG采集仪采集不同G值下载人离心机的M个通道中的脑电信号,其中,M为正整数;
对所述M个通道中的脑电信号进行预处理,以获取所述M个通道中每个通道的低频脑电数据;
根据所述低频脑电数据获取每个所述脑电信号的频段在每个数据段的参数相关率RzB(%);]]>
当所述参数相关率符合昏厥前预警条件时,进行人体昏厥预警提醒。

2.  如权利要求1所述的方法,其特征在于,所述对M个通道中的脑电信号进行预处理具体包括:
通过带通滤波器对所述脑电信号进行滤波,并对滤波后的脑电信号进行伪迹干扰消除处理。

3.  如权利要求1所述的方法,其特征在于,所述根据低频脑电数据获取每个所述脑电信号的频段在每个数据段的参数相关率具体包括:
对所述低频脑电数据进行分析以获取所述脑电信号的特征参数,其中,所述脑电信号的特征参数包括平均幅值平均周期能量比和平均中心频率(F‾z(Hz));]]>
对所述脑电信号的特征参数进行归一化处理,并根据归一化处理后的特征参数获取所述低频脑电数据的每个频段的估计值Bz;
根据所述每个频段的估计值Bz获取每个频段的参数相关率

4.  如权利要求3所述的方法,其特征在于,所述平均幅值通过以下公式计算:
(A‾z(μv))=Σj=116Az(j)/16]]>
其中,Az(j)为第j个通道z频段的脑电信号的幅度,j=1,2,...,16,所述Az(j)可通过以下公式计算:
Az(x)=6Sz(j)]]>
其中,Sz(j)表示第j个通道z频段的脑电信号的能量。

5.  如权利要求3所述的方法,其特征在于,所述平均周期能量比通过以下公式计算:
(D‾z(%))=Σj=116Dz(j)/16]]>
其中,Dz(j)为第j个通道z频段的脑电信号的周期能量比,j=1,2,...,16,所述Dz(j)可通过以下公式计算:
Dz(j)=Sz(j)/ST(j)×100
其中,Sz(j)表示第j个通道z频段的脑电信号的能量,j=1,2,...,16,ST(j)表示第j个通道T(0.5-25Hz)频段的脑电信号的能量。

6.  如权利要求3所述的方法,其特征在于,所述平均中心频率通过以下公式计算:
(F‾z(Hz))=Σj=116Fz(j)/16]]>
其中,Fz(j)为第j个通道z频段脑电信号的中心频率,j=1,2,...,16,所述Fz(j)通过以下公式计算:
(F‾z(Hz))=fzC(j)|Pmax,Pmax=maxflowerfz(x)fupperP(fz(j))]]>
表示第j个通道z频段中最大能量谱的中心频率,flower表示第j个通道z频段的下界;fup表示第j个通道z频段的上界,P(fz(j))表示第j个通道在z频段的能量。

7.  如权利要求3所述的方法,其特征在于,所述对脑电信号的特征参数进行归一化处理,并根据归一化处理后的特征参数获取所述低频脑电数据的每个频段的估计值Bz具体包括:
通过以下公式对所述脑电信号的特征参数进行归一化处理:
ΦzQ(i)=(Qz(i)-minQz(i))/(maxQz(i)-minQz(i))]]>
其中,i表示数据段的序列号,Qz(i)表示或或minQz(i)为Qz(i)中的最小值,maxQz(i)为Qz(i)中的最大值;
在对所述脑电信号的特征参数进行归一化处理后,通过以下公式计算估计值Bz:
Bz(i)=ΦzA‾(i)ΦzD‾(i)ΦzF‾(i).]]>

8.  如权利要求3所述的方法,其特征在于,所述根据估计值Bz获取每个频段的参数相关率具体包括:
通过以下公式计算参数相关率
RzB(i)=Bz(i)/Bz(std)×100,]]>
其中,i表示数据段的序列号,Bz(std)表示数据段为i时,估计值Bz中的标准值。

9.  如权利要求1所述的方法,其特征在于,所述脑电信号的频段包括α频段和β频段,所述α频段和β频段分别为8-13Hz和13-25Hz,所述昏厥前预警条件为:
a.α频段和β频段脑电数据的参数相关率和均大于90%;
b.在满足条件a后的2~10s内,所述参数相关率小于50%,且所述参数相关率小于50%的状态持续段数大于或者等于2。

说明书

说明书基于采集脑电信号的人体昏厥预警方法
技术领域
本发明涉及航空医学及生物医学工程领域,尤其涉及一种基于采集脑电信号的人体昏厥预警方法。
背景技术
高性能战斗机问世以来,高G防护已成为航空医学亟待解决的重大理论和现实问题。重力加速度引起的意识丧失(G-LOC(G-induced loss of consciousness,G-LOC)严重威胁着飞行安全。为此,在飞行和训练中,为了保证飞行员的安全,需要实时监测人体的生理状态以进行判断和预警。
目前,可通过监测飞行员的心电、耳脉信号和脑电信号来了解飞行员的身体状态。其中,心电中的心律失常经常作为训练过程中停机的指标,并且心率是选拔高G特级飞行和储备飞行员的具有价值的重要指标。当心电中的心律失常时,说明人体的心脏已出现严重异常,此时需要立即停机以免对受试者造成更大的伤害。耳脉搏和头水平血压具有一致性,通常耳脉搏可作为判断+Gz耐力终点的客观指标,其中,+Gz表示z轴正方向的重力加速度。通常耳脉搏拉平1~2秒,即可作为黑视或意识丧失的预警。现在普遍认为脑电信号是一个很好的检测指标,目前可通过缺氧、意识丧失、下体负压等各种方法来模拟高G时产生的加速度效果来间接研究高G环境下的脑电变化特征。
目前存在的问题是,由于心率的变化受个体差异影响较大,且和行为及心理活动密切相关,容易受到外界影响,因此,将其作为耐力终点判别意义不大。耳脉反映的是脑颅外动脉压的变化,并不是脑颅内动脉压的变化。此外,由于受温度的影响,所测量的耳脉信号并不准确,且不稳定。脑电信号的研究目前多基于静态脑电数据来研究人体处于昏厥或已经临近昏厥时的脑电变化特征,并且以目测分析暴发性的高幅慢波为特征,当发现肉眼可辨识的“高幅慢波”(比如δ波或θ波)时,脑功能状态可能已处于严重的抑制状态,再进行G-LOC预警已失去意义。此外,由于个体不同产生的差异,目前评价脑电变化的判断依据并不准确,缺乏客观性。
发明内容
本发明旨在至少解决上述技术问题之一。
为此,本发明的目的在于提出一种基于采集脑电信号的人体昏厥预警方法。该方法可对人体大脑昏厥前的状态进行提前识别和预警,从而避免了对受训者造成更大的伤害,对解决高G防护问题具有重要现实意义。
为了实现上述目的,本发明实施例的基于采集脑电信号的人体昏厥预警方法,包括:通过脑电图EEG采集仪采集不同G值下载人离心机的M个通道中的脑电信号,其中,M为正整数;对所述M个通道中的脑电信号进行预处理,以获取所述M个通道中每个通道的低频脑电数据;根据所述低频脑电数据获取每个所述脑电信号的频段在每个数据段的参数相关率当所述参数相关率符合昏厥前预警条件时,进行人体昏厥预警提醒。
本发明实施例的基于采集脑电信号的人体昏厥预警方法,具有以下有益效果:1、通过直接在载人离心机上开展人体实验,研究离心机+Gz下的脑电变化特征,根据脑电变化特征可以充分了解到晕厥前或G-LOC前仅通过肉眼判图而无法了解的有关的脑功能状态的变化信息,并根据参数相关率的变化情况对+Gz引起的晕厥提出预警及判别方法,对于全面了解受训者的体征,指导人体训练具有实际应用价值,对于提醒飞行员在训练和执行飞行任务中及早采取相应的防护措施,解决高G防护问题具有重要现实意义;2、通过参数相关率可以消除个体不同产生的差异,在评价人体的脑电变化时能够给出客观的和较为准确的判断依据。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中,
图1是本发明一个实施例的基于采集脑电信号的人体昏厥预警方法的流程图;
图2是alpha频带的参数相关率随不同G值变化的示意图;以及
图3是beta频带的参数相关率随不同G值变化的示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。相反,本发明的实施例包括落入所附加权利要求书的精神和内涵范围内的所有变化、修改和等同物。
在本发明的描述中,需要理解的是,术语“第一”、“第二”等仅用于描述目的,而不能理解为指示或暗示相对重要性。在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。此外,在本发明的描述中,除非另有说明,“多个”的含义是两个或两个以上。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分,并且本发明的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能,这应被本发明的实施例所属技术领域的技术人员所理解。
下面参考附图描述本发明实施例的基于采集脑电信号的人体昏厥预警方法。
目前,由于在载人离心机上对人体进行试验要求的特殊性以及+Gz环境下脑电信号易受干扰并且脑电信号的测量难度较大。因此,目前对脑电变化特征的研究多通过缺氧、意识丧失、下体负压等各种方法来间接研究高G环境下的脑电变化特征,其主要以目测分析暴发性的高幅慢波为特征。然而,当发现肉眼可辨识的“高幅慢波”时,人体脑功能状态可能已处于严重的抑制状态,此时再对G-LOC进行预警已失去意义。如果可以在不同G值下,直接采集载人离心机中的脑电信号,并对采集到的脑电信号进行综合分析以获取不同G值下人体昏厥前或G-LOC前脑电信号的变化特征,则可以更方便地了解受训者的体征变化,当飞行员或者受训者在训练和执行任务中如果脑电变化符合预设昏厥前则可进行提前预警,并及时采取相应的防护措施,从而避免了对飞行员或者受试者造成更大的伤害。为此,本发明提出了一种基于采集脑电信号的人体昏厥预警方法。
图1是本发明一个实施例的基于采集脑电信号的人体昏厥预警方法的流程图。
如图1所示,该基于采集脑电信号的人体昏厥预警方法包括。
S101,通过脑电图EEG采集仪采集不同G值(即重力加速度值)下载人离心机的M个通道中的脑电信号,其中,M为正整数。
具体而言,在新型三轴向高性能新型载人离心机中进行人体试验。其中,新型三轴向高性能新型载人离心机主臂长度为8m,具有3轴向加速度。每个轴向的加速度变化将会对人体产生不同的影响,Gz方向的加速度对人体脑电信号的影响最大,也就是说,z轴方向的重力加速度对人体脑电信号的影响最大。本发明主要研究+Gz(即,z轴正方向的重力加速度)作用下人体脑电信号的变化情况,也就是说,主要研究受训者在受到z轴正方向的重力加速度而引起的脑电信号的变化情况。因此在设计新型三轴向高性能新型载人离心机的加速度曲 线时,将Gx和Gy方向的加速度值(即x轴和y轴方向的重力加速度)控制地尽可能的小。例如,可将Gx方向的加速度控制在1以下,Gy方向的加速度控制在0.5以下。
进一步而言,每一次运行的加速度曲线的具体设置为:离心机静止时,受试者感受加速度为1G。新型三轴向高性能新型载人离心机启动时先以1G/s的加速度增长率到达基线持续数秒,之后以3G/s的加速度增长率到达而各次运行的设定的最大G值,持续10s-15s后,再以3G/s的加速度增长率降为1G,回到初始状态,离心机停止。其中,直接将受试者曝露于较大G值是危险的,因此每次运行的最大G值从2.5G开始,按照0.5G的幅度递增,直至受试者到达耐力终点或出现停机指标。同时,经验丰富的医生还通过载人离心机的生理信号记录系统对受试者的耳脉和心电信号实时监测并进行记录,每次运行后询问受试者主观对座舱内周边灯和中央灯的视觉感知情况,并结合受试者的表情,对其耐力做出综合判断。
由于动态情况下的脑电信号极其微弱,且易受干扰,因此,动态情况下测量脑电信号的难度相对于静态情况较大。如果电极位置固定不好,或过程中产生松动,则采集到的脑电数据中可用的较少。为了可以准确地采集到动态情况下的脑电数据,可在通过新型三轴向高性能新型载人离心机进行试验时,选用便携式脑电记录仪安装固定在座舱内采集和记录脑电信号。电极安放时,需根据受试者头型的不同大小,为他们佩戴相应型号的紧固网帽,每个电极上加贴医用胶布以加紧固定,防止运转中松动。其中,脑电电极放置采用国际10/20标准,取16路电极,取电极Fp1,F3,C3,P3,O1,Fp2,F4,C4,P4,O2,F7,T3,T5,F8,T4,T6。所有电极的参考电极选用同侧耳垂处的电极A1或A2,选用单极导联的测量方式。具体而言,在不同G值的作用下,通过新型三轴向高性能新型载人离心机中的便携式脑电记录仪可以采集离心机中的l6路中的脑电信号,也就是说通过16路电极即可获取载人离心机中的16个通道中的脑电信号。
在本发明的实施例中,在获取16个通道中的脑电信号之后,脑电信号首先进入座舱内的干扰抑制盒经过初步的干扰处理,之后进入到座舱内的放大记录盒,以特有格式记录到SD卡上进行后续处理。例如,SD卡的采样频率128Hz。其中,干扰抑制盒和放大记录盒的位置可选取+Gz作用下,脑电信号易被采集和记录的合适位置加以固定。
S102,对M个通道中的脑电信号进行预处理,以获取M个通道中每个通道的低频脑电数据。
在本发明的实施例中,对M个通道中的脑电信号进行预处理,具体包括:可通过带通滤波器对脑电信号进行滤波,并可对滤波后的脑电信号进行伪迹干扰消除处理。
具体地,可将采集到的16个通道中的脑电信号通过带通数字滤波器进行预处理,然后选择合适的小波基,设计合理的小波包分解层数,利用小波包分解重构的方法,消除脑电信号中的肌电、工频电源、基线漂移、电极干扰等信号,主要提取出含有alpha (α):8-13Hz、beta(β):13-25Hz、delta(δ):(0.5-4Hz)、theta(θ):(4-8Hz)频带的低频脑电数据。具体而言,动态情况下采集的脑电信号中含有各种伪迹干扰,所以需选用合适的方法对其进行伪迹干扰消除,以提高数据特征提取的性能和效果。小波滤波是时频滤波中的一种方法,其可以使信号在不同部位得到相应最佳的时域和频域分辨率,从而把信号在一系列的不同层次的空间上完成分解和重构。它非常适合于对非平稳信号的瞬态特性和时变特性的分析。也就是说,在低频部分,有着较高的频率分辨率和较低的时间分辨率,而在高频部分,则有着较高的时间分辨率和较低的频率分辨率。而其中新兴的小波包方法,比小波分解更为精细,它在低频和高频部分可同时进行分解。即它能够将频带进行多层次划分,进而对多分辨率分析没有进行细分的高频部分进行进一步分解,而且能够根据被分析信号的特征自适应性地选择相应的频带,使之与信号的频谱进行匹配,进而提高时频上的分辨率。其中,脑电信号中的脑电伪迹主要有肌电图EMG(electromyography)、眼电图EOG(electro-oculogram)、体动、旋转、基线漂移及电源等。受训者在座舱内经历+Gz曝露时,不可避免地会紧张且做一定的对抗动作,因此受肌电的影响较大,该频段主要集中在35.8-51Hz的频段;EOG较难去除,它可能混杂在脑电数据的多个频带之中,在消除过程中很容易造成有用信息的丢失,因此根据实验经验可主要对0.5-1Hz频段内的信号进行去除;交流电的工频集中在50Hz左右;电极固定和基线漂移容易形成0.8Hz和0.2Hz以下的低频慢波;离心机旋转产生的干扰作用较强,不可忽视,根据实验可达到的最大G值,离心机主臂半径,通向心加速度的换算公式可计算出可知离心机旋转对信号的干扰大致主要在0.5Hz频率以下。此外,工频电源、磁场、体动等也会产生不同程度的影响。基于以上综合考虑,滤波器下限取1Hz左右,上限取35Hz左右。根据采样频率和动态脑电信号的特点,选择daubechies5小波对原始信号进行6层分解,其最小分辨率可用下式来估计,式中fs为采样频率通过小波包方法可较好地去除+Gz作用下的脑电信号中的EMG、工频电源、及其它高频干扰,对基线漂移、电极干扰等低频段的慢波干扰也有较明显的效果,去噪后的信噪比大幅度提高。
S103,根据低频脑电数据获取每个脑电信号的频段在每个数据段的参数相关率
在本发明的实施例中,对低频脑电数据进行分析以获取脑电信号的特征参数,其中,脑电信号的特征参数包括脑电信号特征参数包括平均幅值平均周期能量比和平均中心频率对脑电信号的特征参数进行归一化处理,并根据归一化处理后的特征参数获取低频脑电数据的每个频段的估计值Bz;根据每个频段的估计值Bz获取每个频段的参数相关率其中,z表示脑电信号的各个频段,例如z可为α、β频段或者其他频段。
具体而言,在获取低频脑电数据之后,可对低频脑电数据以预设时间间隔进行周期化处理,并对周期化的脑电数据进行快速傅里叶变换FFT(Fast Fourier Transformation)。例如,可将低频脑电数据以2s为一段,分成连续的多段,以及对各个通道每段的脑电数据进行快速傅里叶变换,并做出每个2s数据段的周期图,其中,周期图由16个通道内的以0.5Hz长度分割的频域中的傅里叶成分的平方和构成。以及在获取周期图之后,可通过各周期图参数把各个通道的信号特征表达出来。
在本发明的实施例中,平均幅值可通过以下公式计算:
(A‾z(μv))=Σj=116Az(j)/16]]>
其中,Az(j)为第j个通道z频段的脑电信号的幅度,j=1,2,...,16,其中,Az(j)可通过以下公式计算:
Az(x)=6Sz(j)]]>
其中,Sz(j)表示第j个通道z频段的脑电信号的能量。
在本发明的实施例中,平均周期能量比可通过以下公式计算:
(D‾z(%))=Σj=116Dz(j)/16]]>
其中,Dz(j)为第j个通道z频段的脑电信号的周期能量比,j=1,2,...,16,其中,Dz(j)可通过以下公式计算:
Dz(j)=Sz(j)/ST(j)×100
其中,Sz(j)表示第j个通道z频段的脑电信号的能量,j=1,2,...,16,ST(j)表示第j个通道T(0.5-25Hz)频段的脑电信号的能量。
在本发明的实施例中,平均中心频率可通过以下公式计算:
(F‾z(Hz))=Σj=116Fz(j)/16]]>
其中,Fz(j)为第j个通道z频段脑电信号的中心频率,j=1,2,...,16,其中,Fz(j)可通过以下公式计算:
(F‾z(Hz))=fzC(j)|Pmax,Pmax=maxflowerfz(x)fupperP(fz(j))]]>
表示第j个通道z频段中最大能量谱的中心频率,flower表示第j个通道z频段的下界;fup表示第j个通道z频段的上界,P(fz(j))表示第j个通道在z频段的能量。
在本发明的实施例中,在获取每个频道的脑电信号的平均幅值平均周期能量比和平均中心频率后,可通过以下公式对脑电信号的特征参数进行归一化处理:
ΦzQ(i)=(Qz(i)-minQz(i))/(maxQz(i)-minQz(i))]]>
其中,i表示数据段的序列号,Qz(i)表示或或minQz(i)为Qz(i)中的最小值,maxQz(i)为Qz(i)中的最大值;
具体而言,通过上述归一化公式对脑电信号的平均幅值平均周期能量比和平均中心频率分别进行归一化处理。
在对脑电信号的特征参数进行归一化处理后,可通过以下公式计算估计值Bz:
Bz(i)=ΦzA‾(i)ΦzD‾(i)ΦzF‾(i).]]>
在本发明的实施例中,在获取每个频段的估计值Bz后,可通过以下公式计算参数相关率RzB(%):]]>
RzB(i)=Bz(i)/Bz(std)×100]]>
其中,i表示数据段的序列号,Bz(std)表示数据段为i时,估计值Bz中的标准值。具体而言,数据段为i时,所对应的每个频段的估计值Bz中的最大值通常被定义为该数据段i所对应的标准值。也就是说,不同数据段对应着不同的标准值。通过计算参数相关率可以很好地消除个体不同的差异,因此,在评价脑电变化时能够给出客观的和比较准确的判断依据。
S104,当参数相关率符合昏厥前预警条件时,进行人体昏厥预警提醒。
在本发明的实施例中,通过反复试验发现,不同G值下alpha(α):8-13Hz和beta(β):13-25Hz频段脑电数据的参数相关率变化比较的明显,通过分析参数相关率和是否满足预设昏厥前条件即可实现人体昏厥前预警。
具体而言,图2和图3是受试者在载人离心机中分别经过了4G和5G的重力加速度所引起的脑电信号和耳脉信号的变化图,并且图2和图3中对G值期间载人离心机的生理信号记录系统记录的耳脉信号所对应的耳脉较低和耳脉拉平的状态分别予以标出。对应耳脉的不同变化状态,参数相关率和表现出一定的变化特性。通过分析可以发现当耳脉拉平或者耳脉较低时,参数相关率和表现出较明显的变化特性。由于耳脉拉平持续2s通常被认为人体已经处于晕厥前的状态,因此重点研究对应耳脉拉平时的参数相关率和的变化特性。即在基于脑电信号分析人体的当前状态时,可以基于参数相关率和对人体的当前状态进行分析,并判断参数相关率和是否满足预设昏厥前预警条件。在本发明的实施例中,脑电信号的频段包括α频段和β频段,其中α频段和β频段分别为8-13Hz和13-25Hz,昏厥前预警条件为:a.α频段和β频段脑电数据的参数相关率和均大于90%;b.在满足条件a后的2~10s内,参数相关率小于50%,并且参数相关率小于50%的状态持续段数大于或者等于2。当参数相关率和满足预设昏厥前预警条件,即确定人体处于昏厥前状态,进行昏厥预警提醒,从而可避免对受训者造成 更大的伤害。
本发明实施例的基于采集脑电信号的人体昏厥预警方法,具有以下有益效果:1、通过直接在载人离心机上开展人体实验,研究离心机+Gz下的脑电变化特征,根据脑电变化特征可以充分了解到晕厥前或G-LOC前仅通过肉眼判图而无法了解的有关的脑功能状态的变化信息,并根据参数相关率的变化情况对+Gz引起的晕厥提出预警及判别方法,对于全面了解受训者的体征,指导人体训练具有实际应用价值,对于提醒飞行员在训练和执行飞行任务中及早采取相应的防护措施,解决高G防护问题具有重要现实意义;2、通过参数相关率可以消除个体不同产生的差异,在评价脑电变化时能够给出客观的和较为准确的判断依据。
应当理解,本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中,多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如,如果用硬件来实现,和在另一实施方式中一样,可用本领域公知的下列技术中的任一项或他们的组合来实现:具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路,具有合适的组合逻辑门电路的专用集成电路,可编程门阵列(PGA),现场可编程门阵列(FPGA)等。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同物限定。

基于采集脑电信号的人体昏厥预警方法.pdf_第1页
第1页 / 共12页
基于采集脑电信号的人体昏厥预警方法.pdf_第2页
第2页 / 共12页
基于采集脑电信号的人体昏厥预警方法.pdf_第3页
第3页 / 共12页
点击查看更多>>
资源描述

《基于采集脑电信号的人体昏厥预警方法.pdf》由会员分享,可在线阅读,更多相关《基于采集脑电信号的人体昏厥预警方法.pdf(12页珍藏版)》请在专利查询网上搜索。

1、(10)申请公布号 CN 103845052 A (43)申请公布日 2014.06.11 CN 103845052 A (21)申请号 201410058482.6 (22)申请日 2014.02.20 A61B 5/0476(2006.01) (71)申请人 清华大学 地址 100084 北京市海淀区 100084-82 信箱 申请人 中国人民解放军空军航空医学研究 所 (72)发明人 张涛 李毅峰 邓略 陈勇胜 (74)专利代理机构 北京清亦华知识产权代理事 务所 ( 普通合伙 ) 11201 代理人 张大威 (54) 发明名称 基于采集脑电信号的人体昏厥预警方法 (57) 摘要 本发明。

2、提出一种基于采集脑电信号的人体昏 厥预警方法。 该方法包括 : 通过脑电图EEG采集仪 采集不同 G 值下载人离心机的 M 个通道的脑电信 号, 其中, M为正整数 ; 对M个通道中的脑电信号进 行预处理, 以获取 M 个通道中每个通道的低频脑 电数据 ; 根据低频脑电数据获取每个脑电信号的 频段在每个数据段的参数相关率当参数 相关率符合昏厥前预警条件时, 进行人体 昏厥预警提醒。 本发明实施例的方法, 可对人体昏 厥前的状态进行提前识别和预警, 有利于全面了 解受训者的体征, 指导人体训练具有实际应用价 值。 (51)Int.Cl. 权利要求书 2 页 说明书 7 页 附图 2 页 (19)。

3、中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页 说明书7页 附图2页 (10)申请公布号 CN 103845052 A CN 103845052 A 1/2 页 2 1. 一种基于采集脑电信号的人体昏厥预警方法, 其特征在于, 包括 : 通过脑电图 EEG 采集仪采集不同 G 值下载人离心机的 M 个通道中的脑电信号, 其中, M 为正整数 ; 对所述 M 个通道中的脑电信号进行预处理, 以获取所述 M 个通道中每个通道的低频脑 电数据 ; 根据所述低频脑电数据获取每个所述脑电信号的频段在每个数据段的参数相关率 当所述参数相关率符合昏厥前预警条件时, 进行人体昏厥预警提醒。。

4、 2.如权利要求1所述的方法, 其特征在于, 所述对M个通道中的脑电信号进行预处理具 体包括 : 通过带通滤波器对所述脑电信号进行滤波, 并对滤波后的脑电信号进行伪迹干扰消除 处理。 3. 如权利要求 1 所述的方法, 其特征在于, 所述根据低频脑电数据获取每个所述脑电 信号的频段在每个数据段的参数相关率具体包括 : 对所述低频脑电数据进行分析以获取所述脑电信号的特征参数, 其中, 所述脑电信号 的特征参数包括平均幅值平均周期能量比和平均中心频率 对所述脑电信号的特征参数进行归一化处理, 并根据归一化处理后的特征参数获取所 述低频脑电数据的每个频段的估计值 Bz; 根据所述每个频段的估计值 B。

5、z获取每个频段的参数相关率 4. 如权利要求 3 所述的方法, 其特征在于, 所述平均幅值通过以下公式计算 : 其中, Az(j) 为第 j 个通道 z 频段的脑电信号的幅度, j=1,2, .,16, 所述 Az(j) 可通过 以下公式计算 : 其中, Sz(j) 表示第 j 个通道 z 频段的脑电信号的能量。 5. 如权利要求 3 所述的方法, 其特征在于, 所述平均周期能量比通过以下公式 计算 : 其中, Dz(j) 为第 j 个通道 z 频段的脑电信号的周期能量比, j=1,2, .,16, 所述 Dz(j) 可通过以下公式计算 : Dz(j)=Sz(j)/ST(j)100 其中, S。

6、z(j) 表示第 j 个通道 z 频段的脑电信号的能量, j=1,2, .,16, ST(j) 表示第 j 个通道 T(0.5-25Hz) 频段的脑电信号的能量。 权 利 要 求 书 CN 103845052 A 2 2/2 页 3 6. 如权利要求 3 所述的方法, 其特征在于, 所述平均中心频率通过以下公式 计算 : 其中, Fz(j) 为第 j 个通道 z 频段脑电信号的中心频率, j=1,2, .,16, 所述 Fz(j) 通过 以下公式计算 : 表示第 j 个通道 z 频段中最大能量谱的中心频率, flower表示第 j 个通道 z 频段 的下界 ; fup表示第 j 个通道 z 频。

7、段的上界, P(fz(j) 表示第 j 个通道在 z 频段的能量。 7. 如权利要求 3 所述的方法, 其特征在于, 所述对脑电信号的特征参数进行归一化处 理, 并根据归一化处理后的特征参数获取所述低频脑电数据的每个频段的估计值 Bz具体包 括 : 通过以下公式对所述脑电信号的特征参数进行归一化处理 : 其中, i 表示数据段的序列号, Qz(i) 表示或或minQz(i) 为 Qz(i) 中 的最小值, maxQz(i) 为 Qz(i) 中的最大值 ; 在对所述脑电信号的特征参数进行归一化处理后, 通过以下公式计算估计值 Bz: 8. 如权利要求 3 所述的方法, 其特征在于, 所述根据估计。

8、值 Bz获取每个频段的参数相 关率具体包括 : 通过以下公式计算参数相关率 其中, i 表示数据段的序列号, Bz(std) 表示数据段为 i 时, 估计值 Bz中的标准值。 9. 如权利要求 1 所述的方法, 其特征在于, 所述脑电信号的频段包括 频段和 频 段, 所述 频段和 频段分别为 8-13Hz 和 13-25Hz, 所述昏厥前预警条件为 : a. 频段和 频段脑电数据的参数相关率和均大于 90% ; b. 在满足条件 a 后的 2 10s 内, 所述参数相关率小于 50%, 且所述参数相关率 小于 50% 的状态持续段数大于或者等于 2。 权 利 要 求 书 CN 10384505。

9、2 A 3 1/7 页 4 基于采集脑电信号的人体昏厥预警方法 技术领域 0001 本发明涉及航空医学及生物医学工程领域, 尤其涉及一种基于采集脑电信号的人 体昏厥预警方法。 背景技术 0002 高性能战斗机问世以来, 高 G 防护已成为航空医学亟待解决的重大理论和现实问 题。重力加速度引起的意识丧失 (G LOC(G-induced loss of consciousness, G-LOC) 严 重威胁着飞行安全。为此, 在飞行和训练中, 为了保证飞行员的安全, 需要实时监测人体的 生理状态以进行判断和预警。 0003 目前, 可通过监测飞行员的心电、 耳脉信号和脑电信号来了解飞行员的身体状。

10、态。 其中, 心电中的心律失常经常作为训练过程中停机的指标, 并且心率是选拔高 G 特级飞行 和储备飞行员的具有价值的重要指标。当心电中的心律失常时, 说明人体的心脏已出现严 重异常, 此时需要立即停机以免对受试者造成更大的伤害。耳脉搏和头水平血压具有一致 性, 通常耳脉搏可作为判断 +Gz 耐力终点的客观指标, 其中, +Gz 表示 z 轴正方向的重力加 速度。通常耳脉搏拉平 1 2 秒, 即可作为黑视或意识丧失的预警。现在普遍认为脑电信 号是一个很好的检测指标, 目前可通过缺氧、 意识丧失、 下体负压等各种方法来模拟高 G 时 产生的加速度效果来间接研究高 G 环境下的脑电变化特征。 00。

11、04 目前存在的问题是, 由于心率的变化受个体差异影响较大, 且和行为及心理活动 密切相关, 容易受到外界影响, 因此, 将其作为耐力终点判别意义不大。耳脉反映的是脑颅 外动脉压的变化, 并不是脑颅内动脉压的变化。此外, 由于受温度的影响, 所测量的耳脉信 号并不准确, 且不稳定。脑电信号的研究目前多基于静态脑电数据来研究人体处于昏厥或 已经临近昏厥时的脑电变化特征, 并且以目测分析暴发性的高幅慢波为特征, 当发现肉眼 可辨识的 “高幅慢波” ( 比如 波或 波 ) 时, 脑功能状态可能已处于严重的抑制状态, 再 进行 G-LOC 预警已失去意义。此外, 由于个体不同产生的差异, 目前评价脑电。

12、变化的判断依 据并不准确, 缺乏客观性。 发明内容 0005 本发明旨在至少解决上述技术问题之一。 0006 为此, 本发明的目的在于提出一种基于采集脑电信号的人体昏厥预警方法。该方 法可对人体大脑昏厥前的状态进行提前识别和预警, 从而避免了对受训者造成更大的伤 害, 对解决高 G 防护问题具有重要现实意义。 0007 为了实现上述目的, 本发明实施例的基于采集脑电信号的人体昏厥预警方法, 包 括 : 通过脑电图 EEG 采集仪采集不同 G 值下载人离心机的 M 个通道中的脑电信号, 其中, M 为正整数 ; 对所述M个通道中的脑电信号进行预处理, 以获取所述M个通道中每个通道的低 频脑电数据。

13、 ; 根据所述低频脑电数据获取每个所述脑电信号的频段在每个数据段的参数相 关率当所述参数相关率符合昏厥前预警条件时, 进行人体昏厥预警提醒。 说 明 书 CN 103845052 A 4 2/7 页 5 0008 本发明实施例的基于采集脑电信号的人体昏厥预警方法, 具有以下有益效果 : 1、 通过直接在载人离心机上开展人体实验, 研究离心机 +Gz 下的脑电变化特征, 根据脑电变 化特征可以充分了解到晕厥前或 G-LOC 前仅通过肉眼判图而无法了解的有关的脑功能状 态的变化信息, 并根据参数相关率的变化情况对 +Gz 引起的晕厥提出预警及判别方法, 对 于全面了解受训者的体征, 指导人体训练具。

14、有实际应用价值, 对于提醒飞行员在训练和执 行飞行任务中及早采取相应的防护措施, 解决高 G 防护问题具有重要现实意义 ; 2、 通过参 数相关率可以消除个体不同产生的差异, 在评价人体的脑电变化时能够给出客观的 和较为准确的判断依据。 0009 本发明附加的方面和优点将在下面的描述中部分给出, 部分将从下面的描述中变 得明显, 或通过本发明的实践了解到。 附图说明 0010 本发明上述的和 / 或附加的方面和优点从下面结合附图对实施例的描述中将变 得明显和容易理解, 其中, 0011 图 1 是本发明一个实施例的基于采集脑电信号的人体昏厥预警方法的流程图 ; 0012 图 2 是 alpha。

15、 频带的参数相关率随不同 G 值变化的示意图 ; 以及 0013 图 3 是 beta 频带的参数相关率随不同 G 值变化的示意图。 具体实施方式 0014 下面详细描述本发明的实施例, 所述实施例的示例在附图中示出, 其中自始至终 相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考 附图描述的实施例是示例性的, 仅用于解释本发明, 而不能理解为对本发明的限制。相反, 本发明的实施例包括落入所附加权利要求书的精神和内涵范围内的所有变化、 修改和等同 物。 0015 在本发明的描述中, 需要理解的是, 术语 “第一” 、“第二” 等仅用于描述目的, 而不 能理解为指示或。

16、暗示相对重要性。 在本发明的描述中, 需要说明的是, 除非另有明确的规定 和限定, 术语 “相连” 、“连接” 应做广义理解, 例如, 可以是固定连接, 也可以是可拆卸连接, 或一体地连接 ; 可以是机械连接, 也可以是电连接 ; 可以是直接相连, 也可以通过中间媒介 间接相连。对于本领域的普通技术人员而言, 可以具体情况理解上述术语在本发明中的具 体含义。此外, 在本发明的描述中, 除非另有说明,“多个” 的含义是两个或两个以上。 0016 流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为, 表示包括 一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、 片段或。

17、部 分, 并且本发明的优选实施方式的范围包括另外的实现, 其中可以不按所示出或讨论的顺 序, 包括根据所涉及的功能按基本同时的方式或按相反的顺序, 来执行功能, 这应被本发明 的实施例所属技术领域的技术人员所理解。 0017 下面参考附图描述本发明实施例的基于采集脑电信号的人体昏厥预警方法。 0018 目前, 由于在载人离心机上对人体进行试验要求的特殊性以及 +Gz 环境下脑电信 号易受干扰并且脑电信号的测量难度较大。 因此, 目前对脑电变化特征的研究多通过缺氧、 说 明 书 CN 103845052 A 5 3/7 页 6 意识丧失、 下体负压等各种方法来间接研究高 G 环境下的脑电变化特征。

18、, 其主要以目测分 析暴发性的高幅慢波为特征。 然而, 当发现肉眼可辨识的 “高幅慢波” 时, 人体脑功能状态可 能已处于严重的抑制状态, 此时再对 G-LOC 进行预警已失去意义。如果可以在不同 G 值下, 直接采集载人离心机中的脑电信号, 并对采集到的脑电信号进行综合分析以获取不同 G 值 下人体昏厥前或 G-LOC 前脑电信号的变化特征, 则可以更方便地了解受训者的体征变化, 当飞行员或者受训者在训练和执行任务中如果脑电变化符合预设昏厥前则可进行提前预 警, 并及时采取相应的防护措施, 从而避免了对飞行员或者受试者造成更大的伤害。为此, 本发明提出了一种基于采集脑电信号的人体昏厥预警方法。

19、。 0019 图 1 是本发明一个实施例的基于采集脑电信号的人体昏厥预警方法的流程图。 0020 如图 1 所示, 该基于采集脑电信号的人体昏厥预警方法包括。 0021 S101, 通过脑电图 EEG 采集仪采集不同 G 值 (即重力加速度值) 下载人离心机的 M 个通道中的脑电信号, 其中, M 为正整数。 0022 具体而言, 在新型三轴向高性能新型载人离心机中进行人体试验。 其中, 新型三轴 向高性能新型载人离心机主臂长度为 8m, 具有 3 轴向加速度。每个轴向的加速度变化将会 对人体产生不同的影响, Gz 方向的加速度对人体脑电信号的影响最大, 也就是说, z 轴方向 的重力加速度对。

20、人体脑电信号的影响最大。本发明主要研究 +Gz(即, z 轴正方向的重力加 速度) 作用下人体脑电信号的变化情况, 也就是说, 主要研究受训者在受到 z 轴正方向的重 力加速度而引起的脑电信号的变化情况。 因此在设计新型三轴向高性能新型载人离心机的 加速度曲线时, 将 Gx 和 Gy 方向的加速度值 (即 x 轴和 y 轴方向的重力加速度) 控制地尽可 能的小。例如, 可将 Gx 方向的加速度控制在 1 以下, Gy 方向的加速度控制在 0.5 以下。 0023 进一步而言, 每一次运行的加速度曲线的具体设置为 : 离心机静止时, 受试者感受 加速度为 1G。新型三轴向高性能新型载人离心机启动。

21、时先以 1G/s 的加速度增长率到达基 线持续数秒, 之后以3G/s的加速度增长率到达而各次运行的设定的最大G值, 持续10s-15s 后, 再以 3G/s 的加速度增长率降为 1G, 回到初始状态, 离心机停止。其中, 直接将受试者曝 露于较大G值是危险的, 因此每次运行的最大G值从2.5G开始, 按照0.5G的幅度递增, 直至 受试者到达耐力终点或出现停机指标。同时, 经验丰富的医生还通过载人离心机的生理信 号记录系统对受试者的耳脉和心电信号实时监测并进行记录, 每次运行后询问受试者主观 对座舱内周边灯和中央灯的视觉感知情况, 并结合受试者的表情, 对其耐力做出综合判断。 0024 由于动。

22、态情况下的脑电信号极其微弱, 且易受干扰, 因此, 动态情况下测量脑电信 号的难度相对于静态情况较大。如果电极位置固定不好, 或过程中产生松动, 则采集到的 脑电数据中可用的较少。为了可以准确地采集到动态情况下的脑电数据, 可在通过新型三 轴向高性能新型载人离心机进行试验时, 选用便携式脑电记录仪安装固定在座舱内采集和 记录脑电信号。 电极安放时, 需根据受试者头型的不同大小, 为他们佩戴相应型号的紧固网 帽, 每个电极上加贴医用胶布以加紧固定, 防止运转中松动。其中, 脑电电极放置采用国际 10/20 标准, 取 16 路电极, 取电极 Fp1, F3, C3, P3, O1, Fp2, F。

23、4, C4, P4, O2, F7, T3, T5, F8, T4, T6。所 有电极的参考电极选用同侧耳垂处的电极 A1或 A2, 选用单极导联的测量方式。具体而言, 在不同 G 值的作用下, 通过新型三轴向高性能新型载人离心机中的便携式脑电记录仪可以 采集离心机中的 l6 路中的脑电信号, 也就是说通过 16 路电极即可获取载人离心机中的 16 个通道中的脑电信号。 说 明 书 CN 103845052 A 6 4/7 页 7 0025 在本发明的实施例中, 在获取 16 个通道中的脑电信号之后, 脑电信号首先进入座 舱内的干扰抑制盒经过初步的干扰处理, 之后进入到座舱内的放大记录盒, 以。

24、特有格式记 录到 SD 卡上进行后续处理。例如, SD 卡的采样频率 128Hz。其中, 干扰抑制盒和放大记录 盒的位置可选取 +Gz 作用下, 脑电信号易被采集和记录的合适位置加以固定。 0026 S102, 对M个通道中的脑电信号进行预处理, 以获取M个通道中每个通道的低频脑 电数据。 0027 在本发明的实施例中, 对 M 个通道中的脑电信号进行预处理, 具体包括 : 可通过带 通滤波器对脑电信号进行滤波, 并可对滤波后的脑电信号进行伪迹干扰消除处理。 0028 具体地, 可将采集到的 16 个通道中的脑电信号通过带通数字滤波器进行预处理, 然后选择合适的小波基, 设计合理的小波包分解层。

25、数, 利用小波包分解重构的方法, 消除 脑电信号中的肌电、 工频电源、 基线漂移、 电极干扰等信号, 主要提取出含有 alpha() : 8-13Hz、 beta() : 13-25Hz、 delta() : (0.5-4Hz)、 theta() : (4-8Hz) 频带的低频 脑电数据。具体而言, 动态情况下采集的脑电信号中含有各种伪迹干扰, 所以需选用合适 的方法对其进行伪迹干扰消除, 以提高数据特征提取的性能和效果。小波滤波是时频滤波 中的一种方法, 其可以使信号在不同部位得到相应最佳的时域和频域分辨率, 从而把信号 在一系列的不同层次的空间上完成分解和重构。它非常适合于对非平稳信号的瞬。

26、态特性 和时变特性的分析。也就是说, 在低频部分, 有着较高的频率分辨率和较低的时间分辨率, 而在高频部分, 则有着较高的时间分辨率和较低的频率分辨率。 而其中新兴的小波包方法, 比小波分解更为精细, 它在低频和高频部分可同时进行分解。即它能够将频带进行多层次 划分, 进而对多分辨率分析没有进行细分的高频部分进行进一步分解, 而且能够根据被分 析信号的特征自适应性地选择相应的频带, 使之与信号的频谱进行匹配, 进而提高时频上 的分辨率。其中, 脑电信号中的脑电伪迹主要有肌电图 EMG(electromyography) 、 眼电图 EOG(electro-oculogram) 、 体动、 旋转。

27、、 基线漂移及电源等。受训者在座舱内经历 +Gz 曝 露时, 不可避免地会紧张且做一定的对抗动作, 因此受肌电的影响较大, 该频段主要集中 在 35.8-51Hz 的频段 ; EOG 较难去除, 它可能混杂在脑电数据的多个频带之中, 在消除过程 中很容易造成有用信息的丢失, 因此根据实验经验可主要对 0.5-1Hz 频段内的信号进行去 除 ; 交流电的工频集中在 50Hz 左右 ; 电极固定和基线漂移容易形成 0.8Hz 和 0.2Hz 以下的 低频慢波 ; 离心机旋转产生的干扰作用较强, 不可忽视, 根据实验可达到的最大 G 值, 离心 机主臂半径, 通向心加速度的换算公式可计算出可知离心机。

28、旋转对信号的干扰大致主要在 0.5Hz频率以下。 此外, 工频电源、 磁场、 体动等也会产生不同程度的影响。 基于以上综合考 虑, 滤波器下限取 1Hz 左右, 上限取 35Hz 左右。根据采样频率和动态脑电信号的特点, 选择 daubechies5小波对原始信号进行6层分解, 其最小分辨率可用下式来估计, 式中fs为采样 频率通过小波包方法可较好地去除+Gz作用下的脑电信号中的EMG、 工频 电源、 及其它高频干扰, 对基线漂移、 电极干扰等低频段的慢波干扰也有较明显的效果, 去 噪后的信噪比大幅度提高。 0029 S103, 根据低频脑电数据获取每个脑电信号的频段在每个数据段的参数相关率 。

29、0030 在本发明的实施例中, 对低频脑电数据进行分析以获取脑电信号的特征参数, 其 说 明 书 CN 103845052 A 7 5/7 页 8 中, 脑电信号的特征参数包括脑电信号特征参数包括平均幅值平均周期能量比 和平均中心频率对脑电信号的特征参数进行归一化处理, 并根据归一化 处理后的特征参数获取低频脑电数据的每个频段的估计值Bz; 根据每个频段的估计值Bz获 取每个频段的参数相关率其中, z 表示脑电信号的各个频段, 例如 z 可为 、 频 段或者其他频段。 0031 具体而言, 在获取低频脑电数据之后, 可对低频脑电数据以预设时间间隔 进行周期化处理, 并对周期化的脑电数据进行快速。

30、傅里叶变换 FFT(Fast Fourier Transformation) 。例如, 可将低频脑电数据以 2s 为一段, 分成连续的多段, 以及对各个通 道每段的脑电数据进行快速傅里叶变换, 并做出每个 2s 数据段的周期图, 其中, 周期图由 16个通道内的以0.5Hz长度分割的频域中的傅里叶成分的平方和构成。 以及在获取周期图 之后, 可通过各周期图参数把各个通道的信号特征表达出来。 0032 在本发明的实施例中, 平均幅值可通过以下公式计算 : 0033 0034 其中, Az(j) 为第 j 个通道 z 频段的脑电信号的幅度, j=1,2, .,16, 其中, Az(j) 可通过以下。

31、公式计算 : 0035 0036 其中, Sz(j) 表示第 j 个通道 z 频段的脑电信号的能量。 0037 在本发明的实施例中, 平均周期能量比可通过以下公式计算 : 0038 0039 其中, Dz(j) 为第 j 个通道 z 频段的脑电信号的周期能量比, j=1,2, .,16, 其中, Dz(j) 可通过以下公式计算 : 0040 Dz(j)=Sz(j)/ST(j)100 0041 其中, Sz(j) 表示第 j 个通道 z 频段的脑电信号的能量, j=1,2, .,16, ST(j) 表示 第 j 个通道 T(0.5-25Hz) 频段的脑电信号的能量。 0042 在本发明的实施例中。

32、, 平均中心频率可通过以下公式计算 : 0043 0044 其中, Fz(j) 为第 j 个通道 z 频段脑电信号的中心频率, j=1,2, .,16, 其中, Fz(j) 可通过以下公式计算 : 0045 0046 表示第 j 个通道 z 频段中最大能量谱的中心频率, flower表示第 j 个通道 z 频段的下界 ; fup表示第 j 个通道 z 频段的上界, P(fz(j) 表示第 j 个通道在 z 频段的能量。 说 明 书 CN 103845052 A 8 6/7 页 9 0047 在本发明的实施例中, 在获取每个频道的脑电信号的平均幅值平均周 期能量比和平均中心频率后, 可通过以下公。

33、式对脑电信号的特征参数进 行归一化处理 : 0048 0049 其中, i表示数据段的序列号, Qz(i)表示或或minQz(i)为Qz(i) 中的最小值, maxQz(i) 为 Qz(i) 中的最大值 ; 0050 具体而言, 通过上述归一化公式对脑电信号的平均幅值平均周期能量 比和平均中心频率分别进行归一化处理。 0051 在对脑电信号的特征参数进行归一化处理后, 可通过以下公式计算估计值 Bz: 0052 0053 在本发明的实施例中, 在获取每个频段的估计值 Bz后, 可通过以下公式计算参数 相关率 0054 0055 其中, i 表示数据段的序列号, Bz(std) 表示数据段为 i。

34、 时, 估计值 Bz中的标准值。 具体而言, 数据段为 i 时, 所对应的每个频段的估计值 Bz中的最大值通常被定义为该数据 段 i 所对应的标准值。也就是说, 不同数据段对应着不同的标准值。通过计算参数相关率 可以很好地消除个体不同的差异, 因此, 在评价脑电变化时能够给出客观的和比较 准确的判断依据。 0056 S104, 当参数相关率符合昏厥前预警条件时, 进行人体昏厥预警提醒。 0057 在本发明的实施例中, 通过反复试验发现, 不同 G 值下 alpha () : 8-13Hz 和 beta () : 13-25Hz 频段脑电数据的参数相关率变化比较的明显, 通过分析参数相关率和 是。

35、否满足预设昏厥前条件即可实现人体昏厥前预警。 0058 具体而言, 图 2 和图 3 是受试者在载人离心机中分别经过了 4G 和 5G 的重力加速 度所引起的脑电信号和耳脉信号的变化图, 并且图 2 和图 3 中对 G 值期间载人离心机的生 理信号记录系统记录的耳脉信号所对应的耳脉较低和耳脉拉平的状态分别予以标出。 对应 耳脉的不同变化状态, 参数相关率和表现出一定的变化特性。通过分析可以 发现当耳脉拉平或者耳脉较低时, 参数相关率和表现出较明显的变化特性。 由于耳脉拉平持续 2s 通常被认为人体已经处于晕厥前的状态, 因此重点研究对应耳脉拉 平时的参数相关率和的变化特性。即在基于脑电信号分析。

36、人体的当前状态 时, 可以基于参数相关率和对人体的当前状态进行分析, 并判断参数相关率 和是否满足预设昏厥前预警条件。在本发明的实施例中, 脑电信号的频段包 说 明 书 CN 103845052 A 9 7/7 页 10 括 频段和 频段, 其中 频段和 频段分别为 8-13Hz 和 13-25Hz, 昏厥前预警条件 为 : a. 频段和 频段脑电数据的参数相关率和均大于 90% ; b. 在满足条件 a 后的 2 10s 内, 参数相关率小于 50%, 并且参数相关率小于 50% 的状态持 续段数大于或者等于2。 当参数相关率和满足预设昏厥前预警条件, 即确定人 体处于昏厥前状态, 进行昏厥。

37、预警提醒, 从而可避免对受训者造成更大的伤害。 0059 本发明实施例的基于采集脑电信号的人体昏厥预警方法, 具有以下有益效果 : 1、 通过直接在载人离心机上开展人体实验, 研究离心机 +Gz 下的脑电变化特征, 根据脑电变 化特征可以充分了解到晕厥前或 G-LOC 前仅通过肉眼判图而无法了解的有关的脑功能状 态的变化信息, 并根据参数相关率的变化情况对 +Gz 引起的晕厥提出预警及判别方法, 对 于全面了解受训者的体征, 指导人体训练具有实际应用价值, 对于提醒飞行员在训练和执 行飞行任务中及早采取相应的防护措施, 解决高 G 防护问题具有重要现实意义 ; 2、 通过参 数相关率可以消除个。

38、体不同产生的差异, 在评价脑电变化时能够给出客观的和较为 准确的判断依据。 0060 应当理解, 本发明的各部分可以用硬件、 软件、 固件或它们的组合来实现。在上述 实施方式中, 多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件 或固件来实现。例如, 如果用硬件来实现, 和在另一实施方式中一样, 可用本领域公知的下 列技术中的任一项或他们的组合来实现 : 具有用于对数据信号实现逻辑功能的逻辑门电路 的离散逻辑电路, 具有合适的组合逻辑门电路的专用集成电路, 可编程门阵列 (PGA) , 现场 可编程门阵列 (FPGA) 等。 0061 在本说明书的描述中, 参考术语 “一个实。

39、施例” 、“一些实施例” 、“示例” 、“具体示 例” 、 或 “一些示例” 等的描述意指结合该实施例或示例描述的具体特征、 结构、 材料或者特 点包含于本发明的至少一个实施例或示例中。在本说明书中, 对上述术语的示意性表述不 一定指的是相同的实施例或示例。而且, 描述的具体特征、 结构、 材料或者特点可以在任何 的一个或多个实施例或示例中以合适的方式结合。 0062 尽管已经示出和描述了本发明的实施例, 本领域的普通技术人员可以理解 : 在不 脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、 修改、 替换和变型, 本 发明的范围由权利要求及其等同物限定。 说 明 书 CN 103845052 A 10 1/2 页 11 图 1 图 2 说 明 书 附 图 CN 103845052 A 11 2/2 页 12 图 3 说 明 书 附 图 CN 103845052 A 12 。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 人类生活必需 > 医学或兽医学;卫生学


copyright@ 2017-2020 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备2021068784号-1