一种融合脑电和心电信号的疲劳驾驶检测方法
【专利摘要】本发明涉及一种融合脑电和心电信号的疲劳驾驶检测方法。本发明通过16通道的gUSBamp放大器获取脑电和心电的原始数据，提取脑电的功率谱特征数据—脑电疲劳指数和、脑电的相位同步特征数据—Pz-Fz和P3-P4在delta频段的MPC、心电的时域特征数据—HR和心电的频域特征数据—LF/HF后，直接将脑电心电特征数据在特征数据层面上直接融合，再利用SVM进行分类。本发明强调决策信息的全面性，融合特征的分类效果总体上比单种特征好，运用两种正交的生理指标来检测驾驶疲劳，有助于提高检测的准确率。
【专利说明】—种融合脑电和心电信号的疲劳驾驶检测方法
【技术领域】
[0001]本发明涉及道路交通中的驾驶员疲劳驾驶检测方法，具体涉及一种在数据的特征层上直接将脑电心电特征融合再利用支持向量机进行分类从而得到比单种特征分类效果好的融合特征方法，并且该方法得到的识别率相较于单信号分类的识别率相对有所提升。
【背景技术】
[0002]在驾驶员疲劳驾驶时，其在生理上会表现为:脑电信号会随着人体的精神负荷的变化而变化，心率会随着疲劳的加深而降低。
[0003]人们经过实验证实目前的疲劳驾驶监测基本沿用的是单一的疲劳检测指标，本身分类错误率就较高，而且传统的分类方法(如参数估计与神经网络等)在样本有限的情况下准确率低、识别效果不稳定，这就造成目前疲劳驾驶监测系统低准确率的事实。而支持向量机(Support Vector Machine,SVM)是在统计学习理论的基础上发展而来,根据有限的样本信息寻找最佳结构方案。因此SVM在近几年受到了广泛的关注。在疲劳驾驶监测过程中，脑电与心电信号均是疲劳检测的重要生理信号，若把两者指标进行有效的数据融合再辅助以SVM分类器来判断疲劳程度，则有更高的准确性和可靠性。
【发明内容】
[0004]本发明针对仅有单一疲劳监测指标的疲劳监测手段存在局限性的问题，提出了一种在数据的特征层上直接将脑电心电特征融合再利用支持向量机进行分类从而得到比单种特征分类效果好的融合特征方法。
[0005]本发明采用的技术方案是:通过16通道的gUSBamp放大器(g.Tec公司)获取脑电和心电的原始数据，提取脑电的功率谱特征数据一脑电疲劳指数F(e + a)/e和Fe/a、脑电的相位同步特征数据一Pz-Fz和P3-P4在delta频段的MPC、心电的时域特征数据一HR和心电的频域特征数据一LF/HF后，直接将脑电心电特征数据在特征数据层面上直接融合，再利用SVM进行分类。SVM的具体方法为:先通过合适的非线性映射函数把脑电和心电数据的特征向量映射到足够的高维，找到其对应的分类超平面；再求出使得分类间隔(两类中在分类面的垂直方向最近的距离，如图2)，最大的最优解；最后把最优分类面问题转化为对偶问题求解。
[0006]在SVM中，一个关键的问题是采用适当的核函数K (Xi, Xj)，目前常用的核函数主要有三类:
[0007]多项式核函数:K(x,y) = [ < x, y > +l]q
[0008]高斯径向基核函数(RBF):
【权利要求】
1.一种融合脑电和心电信号的疲劳驾驶检测方法，其特征在于该方法包括以下步骤: 首先计算脑电的功率谱特征数据和相位同步特征数据，具体是: 通过 gUSBamp 放大器来记录 Fz、Pz、Oz、Fpl、Fp2、F7、F3、F4、F8、C3、C4、P7、P3、P4 和P8这15个通道的原始数据； 脑电数据的处理单位为ls，选取F3、Fz, F4、P3、Pz和P4这6个通道的脑电原始数据X (η) = {χ(0), χ(1),...，X (N-1)},采用Welch法进行功率谱估计,选择汉宁窗函数w (η)进行窗处理，则计算通道w的M频段的功率谱估计为:.MU S
其中归一化因子为:
由此可分别计算得到通道w的α，β和Θ频段功率谱估计为I α (W)，I e (W)和I e (w)，其中 w 为 F3、Fz、F4、P3、Pz、P4 ; 贝U，脑电的功率谱特征数据——脑电疲劳指数为:
W = F3、Fz、F4、P3、Pz、P4
w = F3、Fz、F4、P3、Pz、P4 选取Pz-Fz、P3-P4这2对电极通道的脑电原始数据,则脑电相位同步特征数据一Pz-Fz和P3-P4的delta频段的平均相位相干性为:
其中，At为采样周期，N为样本点数，0H(kAt)为在k个采样周期后的信号的瞬时相位，λ的取值范围是[0，1]，若λ =0，表明信号x(t)和y(t)相位完全不同步，若λ = I则表明信号X(t)和y(t)相位完全同步； 然后计算心电时域特征数据和频域特征数据，具体是: 心电信号的原始数据采集是利用通过16通道的gUSBamp放大器的没有用到的最后一个通道得到；一个典型的心电信号包含着P波、QRS波和T波；其中QRS波是最明显的，包含着许多跟心脏状态有关的信息； 由于R波具有变化快，幅度大，持续时间短特征，采用差分阈值法来进行R波的监测；对于得到的滤波后的心电信号X (i)，i = 1,2,...，计算一阶差分和二阶差分为:r ⑴=[X(i+l)-X(1-l)]/2?" i(i)=[广 a+l)-f (1-l)]/2把一阶差分和二阶差分相结合得: F(i) =f' (I)Xfi max+f" (i)Xf" max 给定的阈值为:Rthr = CXfmax 其中，C为经验参数，fmax为F(i)最大值； 若存在连续的F⑴> Rthr,则在R点附近存在一个QRS波，在R点附近查找振幅绝对值最大值点，即为R波的波峰点； 计算窗口时间10s内所有RR期间的平均值为:
贝丨J，心电的时域特征数据一心率HR为:
1s为窗口前移步长，即重叠窗口为90s，重复计算HR，得到驾驶过程中的心率数据；对于心电数据样本时间序列X(tj), j = I, 2, 3,...,N,采用Lomb-Scargle周期图法,计算频率f的周期信号功率为:
其中，h是样本时间，N是样本总数，τ是时间平移变量； 计算低频LF的周期信号功率为:
其中:低频的频率范围为0.04~0.14Hz ； 计算高频HF的周期信号功率为:
其中:高频的频率范围为0.15~0.4Hz ； 贝U，心电的频域特征数据为:LF/HF ； 在上述计算脑电和心电特征数据时，脑电特征提取的窗口为ls，而心电特征的窗口为10s ;为了便于脑电和心电信号的直接比较和融合，对10s内的脑电特征取其均值作为脑电10s的特征数据，以便与心电特征对应； 最后将得到的脑电和心电特征数据整理为:
(Xi, Yi), i = I,..., η, X e Rd, y e {+I, -1} 其中“+I”表示疲劳时的脑电或心电数据，“-1”表示清醒时的脑电或心电数据，并将得到的脑电和心电特征数据作为SVM进行分类的样本数据； 由于脑电和心电数据的特征向量在输入空间Rd中线性不可分，非线性SVM通过非线性映射函数:
Φ: Rd — H 将特征向量X e Rd映射到高维欧几里得空间H ；在非线性SVM情况下，两分类的最优分类面的方程为: Φ (X).w+b = O 之后将问题转化为以下的最优解问题:
最优分类面问题转化为对偶问题:
判断函数为:
【文档编号】A61B5/18GK104127195SQ201410366036
【公开日】2014年11月5日 申请日期:2014年7月29日 优先权日:2014年7月29日
【发明者】孔万增, 周凌霄, 周慧敏, 徐飞鹏, 周展鹏 申请人:杭州电子科技大学