CN107280663A - 一种基于不同实验难度的疲劳脑电特征研究的方法 - Google Patents
一种基于不同实验难度的疲劳脑电特征研究的方法 Download PDFInfo
- Publication number
- CN107280663A CN107280663A CN201710551773.2A CN201710551773A CN107280663A CN 107280663 A CN107280663 A CN 107280663A CN 201710551773 A CN201710551773 A CN 201710551773A CN 107280663 A CN107280663 A CN 107280663A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- difficulty
- fatigue
- eeg signals
- 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
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/369—Electroencephalography [EEG]
-
- 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/25—Bioelectric electrodes therefor
- A61B5/279—Bioelectric electrodes therefor specially adapted for particular uses
- A61B5/291—Bioelectric electrodes therefor specially adapted for particular uses for electroencephalography [EEG]
-
- 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
-
- 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/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Signal Processing (AREA)
- Psychiatry (AREA)
- Artificial Intelligence (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Power Engineering (AREA)
- Psychology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种基于不同实验难度的疲劳脑电特征研究的方法。精神疲劳是一种逐渐的、累积的现象,一般表现在人体生理活动的减弱。本发明通过设置两组不同难度的Stroop实验来诱发不同状态下的疲劳脑电,使用小波包分解和样本熵算法对两组实验中清醒和疲劳状态下的脑电信号进行特征分析。实验结果表明:从清醒到疲劳状态下,α波和θ波相对能量增加,β波相对能量显著降低(P<0.05),参数α/β、(α+θ)/β的比值(P<0.05)随着疲劳的加深逐渐增大,各脑区的样本熵值呈降低趋势;且相较于较高难度的实验组,低难度任务中的参数变化更为明显,因此,参数α/β、(α+θ)/β可以作为衡量脑力疲劳的潜在指标,同时可以验证适当增加实验难度可以在某种程度上对抗脑力疲劳的产生。
Description
技术领域
本发明设计一种基于不同实验难度的疲劳脑电的特征研究方法。具体指,设计两组不同难度的心理学实验来诱发人体的疲劳脑电,分析不同实验难度下疲劳脑电的变化趋势,该发明属于认知神经科学和信息技术领域的结合,属于数字信号处理技术领域。
背景技术
生理性精神疲劳现已成为城市中的重要问题,严重威胁着人们的身体健康和生命财产安全。在交通驾驶、航空航天活动、人机***监控等工作中,操控人员精神疲劳引起的瞬间的注意力分散、反应迟缓或协调性不够,都可能导致极为严重的事故。因此,疲劳脑电的分析和预防工作就变得尤为重要。
目前精神精神疲劳的判定方法大致分为两类:主观评定法和客观评定法。主观评定法主要通过调查问卷的形式进行,这种评定方法可提供精神疲劳的多种信息,如精神疲劳是否产生、产生时间、产生的原因以及疲劳程度,常用的有Piper疲乏量表、Ep-worth嗜睡量表和Stanford疲乏量表等。但主观评定法易受研究者主观判断能力的影响,在检测疲劳状态时有一定的局限性,并且评分标准不易统一,导致检测***不能始终如一正确地报告驾驶员疲劳状态。客观评定法一般使用心理学指标评定法,如大脑信号的测量—脑电图(Electroenceophalograpy,EEG)、眼睛运动信号的测量—眼电图(Electroculography,EOG)、以及心率和心率变异性信号(HRV)的测量—心电图(Electrocardiograph,ECG)。在这些生理信号参数中,脑电信号由于直反映大脑的相关活动,而被认为是分析脑电信号的最可靠的指标,广泛应用于精神疲劳的各项研究。EEG信号是通过收集放置在头皮上的电极的电位变化来反映大脑皮层的活动的重要生理信号。由于与认知刺激的显著关系,EEG信号已被评估为检测认知任务引起的精神疲劳的最合适的方法之一。
脑电信号是反映大脑活动的最佳指标之一,而分析脑电信号的最关键在于脑电信号的特征提取。目前脑电信号的特征提取方法有很多,主要分为时域分析法、频域分析法和非线性分析法。快速傅里叶变换FFT(Fast Fourier Transform)、自回归模型AR(Autoregressive)、功率谱密度PSD(Power Spectrum Density)都属于频域分析方法,通过把幅度随时间变化的脑电信号变换成脑电功率随频率变化的谱图,从而提取出脑电信号的频域特征。但是这些方法的缺点在于只适合分析平稳信号,对于分析脑电信号这类非平稳信号具有很大的局限性。小波变换属于时域分析方法,继承和发扬了短时傅里叶变化局部化的思想,具有多分辨率特性,能很好地分析非平稳信号的局部信息。样本熵是一种典型的非线性分析方法,通过度量脑电信号的复杂度来反映它的非线性特征,采用一个非负数表示时间序列的复杂性,反应时间序列中新信息的产生率,时间序列越复杂,样本熵值越大。
发明内容
针对现有多任务视觉信息下对认知特性研究的需要,本发明使用小波包分解分析两组不同难度Stroop实验下清醒状态到疲劳状态的α波、β波、θ波的能量比值参数α/β、(α+θ)/β的变化趋势,使用样本熵分析不同难度实验下脑电信号的复杂度变化趋势,从而探究疲劳诱发任务中利用EEG特性评估脑电信号的可行性,以及不同任务难度对EEG参数和疲劳产生的影响。
为了实现上述目的,本发明采用如下技术方案:
一种基于不同实验难度的疲劳脑电的特征研究方法,包括如下步骤:
(1)在每个实验的实施间隙,记录被试者的主观疲劳程度;
(2)对实验中被试者按键反应的反应时间和反应正确率进行计算分析;
(3)对采集到的脑电信号进行预处理,利用低通滤波器进行线性矫正,去除基线漂移和工频干扰;
(4)对(3)中预处理后的脑电信号采用独立分量分析算法,去除眼电伪迹;
(5)对去除眼电伪迹后的脑电信号进行小波包分解和重构,计算α,β,θ频带的能量,并分析能量参数F值和R值;
(6)对去除眼电伪迹的脑电信号进行样本熵分析,分析脑电复杂度的变化。
本发明进一步的改进在于:步骤(1)中被试者的主观疲劳评定采用主观疲乏量表,用0-10分值来表示主观疲劳的程度;
本发明进一步的改进在于:步骤(1)中所实施的实验方案为两组不同难度的stroop实验,利用E-Prime软件编写两组实验的相关程序;
本发明进一步的改进在于:步骤(2)中实验具体实现细节是:在每组实验中,电脑显示器上跳动显示不同颜色的字体,被试者通过字体颜色来按下相应的按键,计算机同时记录被试者的反应时间和反应正确率;
本发明进一步的改进在于:利用SPSS13.0中分析两组不同难度实验下,从清醒到疲劳状态时被试者的反应时间和反应正确率,对其求平均并进行配对T检验;
本发明进一步的改进在于:大脑头皮电极的安置采用国际脑电图学会标定的10-20电极导联定位标准,用双耳垂接法,注射导电膏以增强电极的导电性,正确佩戴脑电帽。实验中通过Neuroscan64脑电设备对脑电信号进行实时采集,并通过放大器进行放大,模数转换,输出到计算机中;
本发明进一步的改进在于:步骤(3)中的预处理是利用Matlab软件中的EEGLAB软件包对原始脑电进行相关预处理,原始脑电含有大量噪声和伪迹干扰,通过低通滤波器滤除工频干扰,并去除基线漂移;
本发明进一步的改进在于:步骤(4)中对于原始脑电中的大量眼电伪迹,采用独立分量分析算法进行相应的去除;
本发明进一步的改进在于:步骤(5)中的处理方法是:预处理后的脑电信号采用小波包分解算法进行分解和重构,计算α,β,θ频带的能量,并分析高低两组不同难度实验下,从清醒到疲劳状态时能量参数F值(α/β)和R值((α+θ)/β)的变化趋势;
本发明进一步的改进在于:步骤(6)中的处理方法是:利用样本熵分析算法计算高低两种不同难度实验下,从清醒到疲劳状态时被试者的脑电复杂度变化趋势。
有益效果
现有技术中,脑电信号诱发采用的是单一心理学任务,或者不同形式的多任务实验,未充分考虑同一类型任务中,不同任务难度对疲劳诱发的影响,本发明的有益效果是:对于单一类型的Stroop实验,设计了两组不同的难度实验—低难度组和高难度组。在低难度组中,字体和颜色和字义是相一致的;在高难度组中,字体的颜色和字义是相矛盾的。通过执行不同难度的两组实验来分析实验难度的不同对疲劳诱发的影响。
脑电信号分析部分也结合了小波包分解和样本熵分析两种算法,从不同分析角度来客观描述在不同难度的实验里,从清醒到疲劳状态下脑电信号的整体变化趋势。
附图说明
图1为疲劳脑电采集和分析的整体框架图。
图2为脑电信号的通道导联示意图。
图3位两组Stroop实验中字体显示样式图。
图4为被试者主观疲劳值的变化趋势。
图5为被试者按键的反应时间和反应正确率变化示意图,(a)为被试者反应时间随
时间变化趋势图(b)为被试者平均反应正确率随时间变化趋势图。
图6为预处理前后脑电信号波形图;(a)为清醒状态下原始脑电信号EEG(b)为清醒状态下基线校正后的脑电信号EEG。
图7为小波包分解示意图。
图8为高低难度实验里清醒到疲劳状态下比值参数F值的变化示意图;(a)为低难度组实验清醒到疲劳状态下F值的变化趋势图(b)为高难度组实验清醒到疲劳状态下F值的变化趋势图。
图9为高低难度实验里清醒到疲劳状态下比值参数R值的变化示意图;(a)为低难度组实验清醒到疲劳状态下R值的变化趋势图(b)为高难度组实验清醒到疲劳状态下R值的变化趋势图。
具体实施方式
为使本发明的目的,技术方案和优点更加清楚,下面将结合附图对本发明实施方案作进一步详细阐述。
如图1所示,该图为本实验的整体方案,实验设计部分使用的是经典心理学实验—Stroop实验,使用E-Prime软件设计两种不同难度的Stroop实验,运行在计算机上。实验过程中,被试者分别进行高低难度组的stroop实验,软件通过按键反馈记录被试者反应时间和反应正确率,在每个实验间隙,记录被试者主观疲劳值。Neuroscan64设备用来采集实验整个过程中的脑电信号用于后续分析。实验数据分析部分:主观疲劳数据和行为数据(反应时间和反应正确率)使用SPSS13.0软件配对T检验,原始脑电信号输入到Matlab软件中,使用EEGLAB工具包处理:先进行低通滤波去基线漂移,在使用独立分量分析算法去除相关眼电伪迹;预处理后的脑电信号分别进行小波包分解和样本熵分析,研究两组不同难度实验里从清醒到疲劳状态下疲劳脑电的变化趋势。
如图2所示,该图为大脑头皮电位放置方法。本发明采用国际脑电图学会标定的10-20电极导联定位标准,用双耳垂接法,注射导电膏以增强电极的导电性,正确佩戴脑电帽。选择不同脑区的FP1、FP2、F7、F8、T7、T8、C3、C4、O1、O2十个通道放置电极,以采集不同脑区的脑电信号,通道M1,M2当做参考电极。
图3为两种不同难度的Stroop实验的卡片图,上四个为低难度组卡片,其中卡片的字体颜色和字义是一致的:即“红”字的颜色为红,“黄”字的颜色为黄,“蓝”字的颜色为蓝,“绿”字的颜色为绿;下四个为高难度组的卡片,其中字体的颜色和字义相矛盾:即“红”字的颜色为黄、绿、蓝中任一种,即“黄”字的颜色为红、绿、蓝中任一种,即“绿”字的颜色为黄、红、蓝中任一种,即“蓝”字的颜色为黄、绿、红中任一种。实验实施过程中,被试者根据界面显示卡片的颜色来按下相对应的按键,而非通过字义来按键。
如图4所示:该折线图显示了被试者在整个实验过程中,主观疲劳值的平均变化趋势。每个实验时间大概为一个小时,共分为12段,每段时间为5分钟,在每个实验间隙,记录被试者的主观疲劳值,最后统计所有被试者的平均主观疲劳值变化趋势。图5所示的为被试者行为数据的变化折线。实验过程中,E-Prime软件根据按键反馈,记录被试者的按键反应时间和反应正确率,该图显示了两组不同难度实验过程中,反应时间和反应正确率的变化趋势。
如图6所示,脑电信号采集设备采集到的脑电信号需要使用Matlab中的EEGLAB工具包进行相关特征分析。图中为EEGLAB里的脑电信号显示波形。该图中(a)图形为未处理的原始脑电信号,可看到有明显的基线漂移现象,(b)图为使用了低通滤波器(0-0.5HZ)滤波后的脑电信号波形,可看出,基线漂移现象得到去除,波形得到了很好的线性矫正。去除基线漂移的脑电信号还需用独立分量分析算法来去除相应的眼电伪迹以及EMG等噪声信号,从而获得较好的脑电信号。
如图7所示:由于EEG信号是时变非平稳信号,选择一种合适的方法来更好的获取反映大脑活动和生理状态的有效信息是进行EEG分析的一个重要前提。本发明采用小波包分解技术对信号进行特征提取。该方法是小波分解的推广,具有任意多尺度的特点,并且克服了小波分解中频率分辨率会随着信号频率的升高而降低的缺陷,能够对信号进行更精细的分析。本发明对预处理的脑电信号进行3层小波包分解,最后得到8个子带频段。
如图8和图9所示,正常人的脑电主要由四种频率成分组成:α波:8-13HZ,β波:14-30HZ,θ波:4-8HZ,δ波:0.5-4HZ。α波是一种慢波,是确定脑电快慢的基准波,一般出现在平静、闭眼清醒的状态下,当睁眼、困倦或者睡眠时,α波会减少;β波表示人的大脑皮层处于兴奋状态,属于脑电信号中的快波;θ波和δ波都属于慢波,在人体处于睡眠状态时出现。一般认为从正常清醒状态向疲劳状态变化的过程中,快波逐渐减少,慢波逐渐增多。在本发明两组不同难度的实验中,从清醒到疲劳状态下,α,β,θ波段的能量发生变化,其中比值参数F值(α/β)和R值((α+θ)/β)在不同脑区的整体变化柱状图可在图中获得。三个波段的能量值可由下列算法得出:
以f(t)表示源信号,经过小波分解后,在第i分解层得到2i个子频带,故源信号f(t)可表示为:
式中:j=0,1,....2i-1,fi,j(tj)为小波包分解在第i层节点(i,j)上的重构信号。根据Parseval定理及式(1)可计算得到信号f(t)小波包分解的能量谱如下:
式中:Ei,j(tj)为f(t)小波包分解到节点(i,j)上的频带能量;xj,k(j=0,1,2....2i-1,k=1,2....m)为重构信号fi,j(tj)的离散点幅值;m为信号采样点数。
得到各节点的频带能量值,便可计算小波包熵值:
其中Pj为各节点能量值在能量总和中所占的比例值,WE为各频带的小波熵值。
本发明除了使用小波包分解算法之外,也使用了样本熵算法分析两组不同难度实验里,从清醒到疲劳状态下脑电复杂度的变化。样本熵是由Richman提出的一种新的时间序列复杂度测度方法,是近似熵的一种改进算法,比近似熵需要更少的数据量。具体算法实现步骤如下:
①设给定一维离散时间序列,构造成一组m维矢量,从Xm(1)到Xm(N-m+1),其中:
Xm(i)=[ui,ui+1,ui+2....ui+m-1](i=1~N-m+1)
②定义任意两个m维向量之间的距离为:
d[Xm(i),Xm(j)]=max|ui+k-uj+k|,0≤k≤m-1;i,j=1~N-m+1,i≠j
③给定阈值r,对每个i值,统计d[X(i),X(j)]<r的个数然后计算其与距离总数(N-m)的比值(称为模板匹配数),记:
其中,SD是一维离散时间序列的标准差;
④计算:
⑤对于m+1点矢量,按以上步骤,可得:
⑥理论上此序列的样本熵值为:SampEn(m,r)=lim{-ln[Bm+1(r)/Bm(r)]},当序列长度为有限值时,其样本熵估计值为:SampEn(m,r)=-ln[Bm+1(r)/Bm(r)]
其中参数m是预先选定的模式维数,r是预先选定的相似容限。
具体样本值变化参见下表:
低难度组样本熵平均值变化
高难度组样本熵平均值变化
上两张表显示了高低难度两组实验里清醒态到疲劳态的样本熵值变化,可以看出被试者的脑电复杂度变化情形。
本发明从不同角度分析了两组不同难度实验里,从清醒到疲劳下的疲劳脑电变化趋势:可看出,相比较高难度组实验,低难度组实验更易诱发人体疲劳的产生,增加实验难度可在一定程度上抑制疲劳的产生。本发明虽考虑到了单一任务的不同难度设置,但目前仅设置了两个难度等级,未设置多个难度等级,后续实验可充分考虑多个难度等级对疲劳脑电诱发的影响。
虽然本发明已以较佳实施例公开如上,但实施例并不是用来限定本发明的。在不脱离本发明之精神和范围内,所做的任何等效变换或润饰,同样属于本发明之保护范围。因此本发明的保护范围应当以本申请的权利要求所界定的内容为标准。
Claims (7)
1.一种基于不同实验难度的疲劳脑电特征研究的方法,其特征在于,包括以下步骤:
(1)各脑区的头皮电极放置:利用置于大脑头皮的各部位电极来采集相应脑区的脑电信号,并对其进行放大,模数转换,以数字信号存储在计算机中;
(2)受试者在视觉刺激实验中完成相关按键任务,记录在不同实验难度里受试者从清醒到疲劳状态下的脑电信号,以及相应的主观数据和行为数据;
(3)对记录的主观数据和行为数据求平均和配对T检验,对采集到的脑电信号进行预处理,利用低通滤波器去除基线漂移,独立分量分析ICA去除眼电伪迹;
(4)对预处理后的脑电信号进行小波包分解,计算各频段的能量,并分析其比值参数F和R值;
(5)对预处理后的脑电信号进行样本熵分析,计算从清醒到疲劳状态下脑电信号的脑电复杂度变化趋势;
(6)分析处理后的主观数据、行为数据以及F值(α/β)、R值((α+θ)/β)、样本熵值,来分析高低两种不同的实验难度对疲劳诱发的影响,分析其相应规律。
2.根据权利要求1所述的方法,其特征在于,步骤(1)中各脑区头皮电极的具体放置方法是:采用国际脑电图学会标定的10-20电极导联定位标准,用双耳垂接法,记录导联FP1、FP2、F7、F8、T7、T8、C3、C4、O1、O2,其中参考电极选取M1、M2,采样频率为125HZ,各通道导联阻抗均小于5kΩ;使用Neuroscan64设备采集脑电信号并进行放大和模数转换,再输入到计算机中。
3.根据权利要求1所述的方法,其特征在于,步骤(2)中的视觉刺激实验设计策略是:使用E-Prime软件设计不同难度的视觉刺激实验—Stroop实验,被试者进行两组实验的测试,记录其主观数据和相应行为数据;
所述的视觉刺激实验主要是:通过设计两组不同难度的Stroop刺激实验来诱发人体的疲劳脑电。低难度实验为红黄蓝绿四种颜色的卡片在屏幕上滚动,字体的颜色和字义保持一致;高难度组实验为字体颜色和字义相矛盾的红黄蓝绿四种颜色卡片在屏幕上滚动。被试者通过观察屏幕上滚动的卡片,根据字体的颜色来分别按下与四种颜色相对应的键盘按键。
4.根据权利要求1所述的方法,其特征在于,步骤(3)中的数据分析具体方法是:使用SPSS13.0软件对主观数据和行为数据求平均并进行配对t检验,对采集的原始脑电信号进行低通滤波去除基线漂移,线性矫正,独立分量分析去除眼电伪迹。
5.根据权利要求1所述的方法,其特征在于,步骤(4)中预处理后的脑电信号进行小波包分解,具体分解算法如下:
以f(t)表示源信号,经过小波分解后,在第i分解层得到2i个子频带,故源信号f(t)可表示为:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<msup>
<mn>2</mn>
<mi>i</mi>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mo>...</mo>
<mo>...</mo>
<mo>+</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<msup>
<mn>2</mn>
<mi>i</mi>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mrow>
<msup>
<mn>2</mn>
<mi>i</mi>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
式中:j=0,1,....2i-1,fi,j(tj)为小波包分解在第i层节点(i,j)上的重构信号;根据Parseval定理及式(1)可计算得到信号f(t)小波包分解的能量谱如下:
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mo>&Integral;</mo>
<mi>&Gamma;</mi>
</munder>
<mo>|</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mi>d</mi>
<mi>t</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<mo>|</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>j</mi>
<mo>.</mo>
<mi>k</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
式中:Ei,j(tj)为f(t)小波包分解到节点(i,j)上的频带能量;xj,k(j=0,1,2....2i-1,k=1,2....m)为重构信号fi,j(tj)的离散点幅值;m为信号采样点数;
得到各节点的频带能量值,便可计算小波包熵值:
<mrow>
<msub>
<mi>P</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>E</mi>
<mi>j</mi>
</msub>
<mi>E</mi>
</mfrac>
</mrow>
<mrow>
<mi>W</mi>
<mi>E</mi>
<mo>=</mo>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mi>j</mi>
</munder>
<msub>
<mi>P</mi>
<mi>j</mi>
</msub>
<msub>
<mi>log</mi>
<mn>2</mn>
</msub>
<msub>
<mi>P</mi>
<mi>j</mi>
</msub>
</mrow>
其中Pj为各节点能量值在能量总和中所占的比例值,WE为各频带的小波熵值。
6.根据权利要求1所述的方法,其特征在于,步骤(5)中预处理后脑电信号的脑电复杂度由下列样本熵算法得出:
①设给定一维离散时间序列,构造成一组m维矢量,从Xm(1)到Xm(N-m+1),其中:
Xm(i)=[ui,ui+1,ui+2....ui+m-1](i=1~N-m+1)
②定义任意两个m维向量之间的距离为:
d[Xm(i),Xm(j)]=max|ui+k-uj+k|,0≤k≤m-1;i,j=1~N-m+1,i≠j
③给定阈值r,对每个i值,统计d[X(i),X(j)]<r的个数然后计算其与距离总数(N-m)的比值(称为模板匹配数),记:
其中,SD是一维离散时间序列的标准差;
④计算:
⑤对于m+1点矢量,按以上步骤,可得:
<mrow>
<msup>
<mi>B</mi>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mi>m</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mi>m</mi>
</mrow>
</munderover>
<msubsup>
<mi>B</mi>
<mi>i</mi>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
⑥理论上此序列的样本熵值为:SampEn(m,r)=lim{-ln[Bm+1(r)/Bm(r)]},当序列长度为有限值时,其样本熵估计值为:SampEn(m,r)=-ln[Bm+1(r)/Bm(r)]
其中参数m是预先选定的模式维数,r是预先选定的相似容限。
7.根据权利要求1所述的方法,其特征在于,步骤(6)中两组实验下各项数据的对比分析:分析高低两种不同实验难度下的主观数据、行为数据、F值、R值以及样本熵值变化趋势,可看出:低难度实验下人体疲劳的产生快于高难度实验下人体疲劳的产生,实验难度的提高一定程度上可抑制疲劳的产生。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710551773.2A CN107280663A (zh) | 2017-07-07 | 2017-07-07 | 一种基于不同实验难度的疲劳脑电特征研究的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710551773.2A CN107280663A (zh) | 2017-07-07 | 2017-07-07 | 一种基于不同实验难度的疲劳脑电特征研究的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107280663A true CN107280663A (zh) | 2017-10-24 |
Family
ID=60100058
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710551773.2A Pending CN107280663A (zh) | 2017-07-07 | 2017-07-07 | 一种基于不同实验难度的疲劳脑电特征研究的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107280663A (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108236464A (zh) * | 2017-12-29 | 2018-07-03 | 重庆邮电大学 | 基于脑电信号的特征提取方法及其检测提取*** |
CN108681391A (zh) * | 2018-03-19 | 2018-10-19 | 南京邮电大学 | 一种基于多模式的脑电信号虚拟键盘设计方法 |
CN109032347A (zh) * | 2018-07-06 | 2018-12-18 | 昆明理工大学 | 一种基于眼电信号控制鼠标的方法 |
CN109443528A (zh) * | 2018-11-16 | 2019-03-08 | 国网江苏省电力有限公司盐城供电分公司 | 一种基于振动信号分析的变压器故障诊断***及其诊断方法 |
CN109480871A (zh) * | 2018-10-30 | 2019-03-19 | 北京机械设备研究所 | 一种面向rsvp脑-机接口的疲劳检测方法 |
CN109875552A (zh) * | 2019-02-01 | 2019-06-14 | 五邑大学 | 一种疲劳检测方法、装置及其存储介质 |
CN110074783A (zh) * | 2019-05-17 | 2019-08-02 | 杭州师范大学 | 经颅磁刺激诱发信号的大脑皮层兴奋性及成像和量化方法 |
CN110742602A (zh) * | 2019-10-15 | 2020-02-04 | 武汉理工大学 | 一种基于脑电与车辆行驶数据的激进驾驶状态识别方法 |
WO2020186915A1 (zh) * | 2019-03-15 | 2020-09-24 | 华为技术有限公司 | 一种注意力检测方法及*** |
CN112221003A (zh) * | 2020-10-14 | 2021-01-15 | 电子科技大学 | 一种基于视听刺激快速诱发疲劳状态的模拟方法及*** |
CN112270235A (zh) * | 2020-10-20 | 2021-01-26 | 西安工程大学 | 一种改进的svm脑电信号情绪识别方法 |
CN112426162A (zh) * | 2020-11-23 | 2021-03-02 | 重庆邮电大学 | 一种基于脑电信号节律熵的疲劳检测方法 |
CN112450879A (zh) * | 2020-10-31 | 2021-03-09 | 复旦大学 | 一种环境光调制大脑认知功能的近红外光谱探测方法 |
CN114209325A (zh) * | 2021-12-23 | 2022-03-22 | 东风柳州汽车有限公司 | 驾驶员疲劳行为监控方法、装置、设备及存储介质 |
CN114224359A (zh) * | 2021-12-18 | 2022-03-25 | 西南交通大学 | 一种针对高铁调度员的脑电疲劳指标判定方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5549118A (en) * | 1993-09-22 | 1996-08-27 | New York University | Simultaneous testing of two motor capabilities of a human subject |
KR20120109160A (ko) * | 2011-03-28 | 2012-10-08 | 한국과학기술원 | 3차원 영상의 시각 피로감 측정 장치 및 유발 요인 추출 장치, 그리고 그 방법 |
KR20140075049A (ko) * | 2012-12-07 | 2014-06-19 | 상명대학교서울산학협력단 | 집중력 평가 방법 및 이를 적용하는 장치 |
CN103961091A (zh) * | 2014-04-15 | 2014-08-06 | 杭州电子科技大学 | 基于双树复小波样本熵的运动想象脑电信号特征提取方法 |
CN104814735A (zh) * | 2015-05-22 | 2015-08-05 | 京东方科技集团股份有限公司 | 判断大脑是否疲劳的方法和装置 |
CN106371610A (zh) * | 2016-09-23 | 2017-02-01 | 重庆金瓯科技发展有限责任公司 | 一种基于脑电信号的驾驶疲劳的检测方法 |
CN106913333A (zh) * | 2015-12-28 | 2017-07-04 | 西南交通大学 | 一种持续性注意水平的敏感性特征指标的选取方法 |
-
2017
- 2017-07-07 CN CN201710551773.2A patent/CN107280663A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5549118A (en) * | 1993-09-22 | 1996-08-27 | New York University | Simultaneous testing of two motor capabilities of a human subject |
KR20120109160A (ko) * | 2011-03-28 | 2012-10-08 | 한국과학기술원 | 3차원 영상의 시각 피로감 측정 장치 및 유발 요인 추출 장치, 그리고 그 방법 |
KR20140075049A (ko) * | 2012-12-07 | 2014-06-19 | 상명대학교서울산학협력단 | 집중력 평가 방법 및 이를 적용하는 장치 |
CN103961091A (zh) * | 2014-04-15 | 2014-08-06 | 杭州电子科技大学 | 基于双树复小波样本熵的运动想象脑电信号特征提取方法 |
CN104814735A (zh) * | 2015-05-22 | 2015-08-05 | 京东方科技集团股份有限公司 | 判断大脑是否疲劳的方法和装置 |
CN106913333A (zh) * | 2015-12-28 | 2017-07-04 | 西南交通大学 | 一种持续性注意水平的敏感性特征指标的选取方法 |
CN106371610A (zh) * | 2016-09-23 | 2017-02-01 | 重庆金瓯科技发展有限责任公司 | 一种基于脑电信号的驾驶疲劳的检测方法 |
Non-Patent Citations (2)
Title |
---|
乔晓艳 等: "基于少次相干平均和样本熵的视听诱发脑电特征提取", 《测试技术学报》 * |
范晓丽 等: "基于EEG的脑力疲劳特征研究", 《北京航空航天大学学报》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108236464B (zh) * | 2017-12-29 | 2021-02-23 | 重庆邮电大学 | 基于脑电信号的特征提取方法及其检测提取*** |
CN108236464A (zh) * | 2017-12-29 | 2018-07-03 | 重庆邮电大学 | 基于脑电信号的特征提取方法及其检测提取*** |
CN108681391A (zh) * | 2018-03-19 | 2018-10-19 | 南京邮电大学 | 一种基于多模式的脑电信号虚拟键盘设计方法 |
CN109032347A (zh) * | 2018-07-06 | 2018-12-18 | 昆明理工大学 | 一种基于眼电信号控制鼠标的方法 |
CN109480871A (zh) * | 2018-10-30 | 2019-03-19 | 北京机械设备研究所 | 一种面向rsvp脑-机接口的疲劳检测方法 |
CN109443528A (zh) * | 2018-11-16 | 2019-03-08 | 国网江苏省电力有限公司盐城供电分公司 | 一种基于振动信号分析的变压器故障诊断***及其诊断方法 |
CN109875552A (zh) * | 2019-02-01 | 2019-06-14 | 五邑大学 | 一种疲劳检测方法、装置及其存储介质 |
CN109875552B (zh) * | 2019-02-01 | 2021-09-28 | 五邑大学 | 一种疲劳检测方法、装置及其存储介质 |
WO2020186915A1 (zh) * | 2019-03-15 | 2020-09-24 | 华为技术有限公司 | 一种注意力检测方法及*** |
CN110074783B (zh) * | 2019-05-17 | 2021-07-20 | 杭州师范大学 | 经颅磁刺激诱发信号的大脑皮层兴奋性及成像和量化方法 |
CN110074783A (zh) * | 2019-05-17 | 2019-08-02 | 杭州师范大学 | 经颅磁刺激诱发信号的大脑皮层兴奋性及成像和量化方法 |
CN110742602A (zh) * | 2019-10-15 | 2020-02-04 | 武汉理工大学 | 一种基于脑电与车辆行驶数据的激进驾驶状态识别方法 |
CN112221003A (zh) * | 2020-10-14 | 2021-01-15 | 电子科技大学 | 一种基于视听刺激快速诱发疲劳状态的模拟方法及*** |
CN112221003B (zh) * | 2020-10-14 | 2021-07-20 | 电子科技大学 | 一种基于视听刺激快速诱发疲劳状态的模拟方法及*** |
CN112270235A (zh) * | 2020-10-20 | 2021-01-26 | 西安工程大学 | 一种改进的svm脑电信号情绪识别方法 |
CN112450879A (zh) * | 2020-10-31 | 2021-03-09 | 复旦大学 | 一种环境光调制大脑认知功能的近红外光谱探测方法 |
CN112450879B (zh) * | 2020-10-31 | 2022-03-18 | 复旦大学 | 一种环境光调制大脑认知功能的近红外光谱探测方法 |
CN112426162A (zh) * | 2020-11-23 | 2021-03-02 | 重庆邮电大学 | 一种基于脑电信号节律熵的疲劳检测方法 |
CN114224359A (zh) * | 2021-12-18 | 2022-03-25 | 西南交通大学 | 一种针对高铁调度员的脑电疲劳指标判定方法 |
CN114224359B (zh) * | 2021-12-18 | 2023-06-30 | 西南交通大学 | 一种针对高铁调度员的脑电疲劳指标判定方法 |
CN114209325A (zh) * | 2021-12-23 | 2022-03-22 | 东风柳州汽车有限公司 | 驾驶员疲劳行为监控方法、装置、设备及存储介质 |
CN114209325B (zh) * | 2021-12-23 | 2023-06-23 | 东风柳州汽车有限公司 | 驾驶员疲劳行为监控方法、装置、设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107280663A (zh) | 一种基于不同实验难度的疲劳脑电特征研究的方法 | |
CN102715903B (zh) | 基于定量脑电图的脑电特征提取方法 | |
CN109770924B (zh) | 一种基于广义一致性构建脑功能网络与相关向量机的疲劳分类方法 | |
Sharma et al. | EEG and cognitive biomarkers based mild cognitive impairment diagnosis | |
Vaid et al. | EEG signal analysis for BCI interface: A review | |
CN103399639B (zh) | 基于ssvep和p300联合脑机接口方法及装置 | |
CN108577834B (zh) | 一种用于癫痫间期棘波自动检测的方法 | |
CN106691474A (zh) | 融合脑电信号与生理信号的疲劳检测*** | |
CN103610447A (zh) | 一种基于前额脑电信号的脑力负荷在线检测方法 | |
CN112137616B (zh) | 一种多感觉脑-体联合刺激的意识检测装置 | |
CN105942974A (zh) | 一种基于低频脑电的睡眠分析方法及*** | |
CN112426162A (zh) | 一种基于脑电信号节律熵的疲劳检测方法 | |
Zhang et al. | An EEG based pervasive depression detection for females | |
CN102824171A (zh) | 脑卒中后抑郁症psd患者脑电特征提取方法 | |
Djamal et al. | Significant variables extraction of post-stroke EEG signal using wavelet and SOM kohonen | |
CN105286860A (zh) | 一种基于双树复小波能量差的运动想象脑电信号识别方法 | |
CN108836301A (zh) | 一种基于奇异谱分析和稀疏表示的单次诱发电位提取方法 | |
CN110367975A (zh) | 一种基于脑机接口的疲劳驾驶检测预警方法 | |
Peng et al. | An improved EEG de-noising approach in electroencephalogram (EEG) for home care | |
CN114668373A (zh) | 一种基于动态模态分解的睡眠阶段分类方法及*** | |
Xu et al. | Fatigue EEG feature extraction based on tasks with different physiological states for ubiquitous edge computing | |
Qi et al. | Driver fatigue assessment based on the feature fusion and transfer learning of EEG and EMG | |
Kaur et al. | Drowsiness detection based on EEG signal analysis using EMD and trained neural network | |
Zammouri et al. | Brain waves-based index for workload estimation and mental effort engagement recognition | |
Smitha et al. | Classifying subjective emotional stress response evoked by multitasking using EEG |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171024 |