CN114916937A - 一种基于bdpca聚类算法的驾驶员脑电疲劳等级划分方法 - Google Patents

一种基于bdpca聚类算法的驾驶员脑电疲劳等级划分方法 Download PDF

Info

Publication number
CN114916937A
CN114916937A CN202210650658.1A CN202210650658A CN114916937A CN 114916937 A CN114916937 A CN 114916937A CN 202210650658 A CN202210650658 A CN 202210650658A CN 114916937 A CN114916937 A CN 114916937A
Authority
CN
China
Prior art keywords
electroencephalogram
fatigue
bdpca
clustering algorithm
clustering
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
Application number
CN202210650658.1A
Other languages
English (en)
Inventor
吕超
刘云清
颜飞
张琼
孔祥尉
闫超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changchun University of Science and Technology
Original Assignee
Changchun University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changchun University of Science and Technology filed Critical Changchun University of Science and Technology
Priority to CN202210650658.1A priority Critical patent/CN114916937A/zh
Publication of CN114916937A publication Critical patent/CN114916937A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
    • A61B5/18—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state for vehicle drivers or machine operators
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
    • A61B5/168—Evaluating attention deficit, hyperactivity
    • 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/316—Modalities, i.e. specific diagnostic methods
    • A61B5/369—Electroencephalography [EEG]
    • A61B5/372—Analysis of electroencephalograms
    • A61B5/374—Detecting the frequency distribution of signals, e.g. detecting delta, theta, alpha, beta or gamma waves
    • 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/7235—Details of waveform analysis
    • 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/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • 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/7253—Details of waveform analysis characterised by using transforms
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2503/00—Evaluating a particular growth phase or type of persons or animals
    • A61B2503/20—Workers
    • A61B2503/22—Motor vehicles operators, e.g. drivers, pilots, captains

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Psychiatry (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Medical Informatics (AREA)
  • Signal Processing (AREA)
  • Psychology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Developmental Disabilities (AREA)
  • Educational Technology (AREA)
  • Child & Adolescent Psychology (AREA)
  • Social Psychology (AREA)
  • Hospice & Palliative Care (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明公开了一种基于BDPCA聚类算法的驾驶员脑电疲劳等级划分方法，通过小波包分解提取脑电信号的节律能量特征和脑电信号的非线性特征作为特征向量。使用共同邻域参数(CNN)改进的DPCA聚类算法对提取到的特征向量进行分析训练。同时，使用贝叶斯准则(BIC)对类簇中心个数进行辅助判定。实验结果表明，改进后的BDPCA聚类算法可以自动对脑电信号中表征的不同疲劳状态等级进行自动划分，客观的定义了疲劳状态等级。BDPCA聚类算法对不同的疲劳状态等级进行了划分，而又无需主观判断提前对信号所属状态进行判定，是一种能够客观划分疲劳状态等级的可行的算法。

Description

一种基于BDPCA聚类算法的驾驶员脑电疲劳等级划分方法
技术领域
本发明涉及疲劳研究技术领域，更具体地说是一种基于BDPCA聚类算法的驾驶员脑电疲劳等级划分方法。
背景技术
世卫组织的报告中说明了道路交通事故已经成为了全球一项重要的死亡原因。而在频发的交通事故中，由于驾驶员的疲劳驾驶导致的交通事故占据着相当大的比例。因此，如果能够对驾驶员的疲劳状态进行检测，并进行相应的等级划分，就能够充分地掌握驾驶员的精神状态。根据驾驶员的疲劳等级状态进行相应的工作划分，这样便能大大的减少疲劳驾驶所导致的交通事故的发生。近些年，对驾驶员疲劳状态的客观检测研究主要围绕肌电、眼电、心电以及脑电等生物电信号进行。除此之外，还有研究人员使用眨眼频率、方向盘轴承偏移量等行为信号进行驾驶员的疲劳状态分析。但是由于生物电信号具有良好的实时差异性和不可伪装性等优点，其对疲劳状态检测的精确度要高于行为信号，而脑电信号相比于其它生物电信号地采集更加简单方便，对机体状态的表征更为全面。因此使用脑电信号对驾驶员的疲劳状态进行检测被认为是最合适的方法之一。随着人工智能技术的快速发展，基于人工智能的分析方法在数据的处理和分类问题上越来越受到关注。Jianfeng Hu将K-NN和SVM应用到脑电信号的分析处理中。闵建亮使用极限学习机对从脑电信号中提取到的特征值进行分析进行对疲劳状态的检测。Saeedi MaryamSVM和多层感知器应用在脑电信号的分析处理中。上述的机器学习方法虽然在对清醒状态与疲劳状态的区分检测中取到了不错的识别精度。但是，仅仅针对疲劳与清醒进行相应的划分，而忽略了其中两种状态之间的过渡信息，并没有建立一个客观明确的疲劳等级划分机制。由于前期进行训练的数据状态大多是由个体的主观感受作为依据，因此难免会出现误差。
发明内容
为了克服现有技术的上述缺陷，本发明采用无监督学习方式的聚类学习算法直接对数据进行处理。相对于有监督的学习方式，无监督学习的聚类算法自动的将疲劳状态相近的脑电信号划分为同一等级，实现对疲劳等级的客观划分。DPCA算法基于数据的局部密度与数据点之间的距离进行分析，但是DPCA算法无法自动确定类簇中心K的数值。因此本方案引入贝叶斯信息准则(BIC)对DPCA聚类算法进行改进的BDPCA算法对数据进行分析。
为实现上述目的，本发明提供如下技术方案：
一种基于BDPCA聚类算法的驾驶员脑电疲劳等级划分方法，包括以下步骤：
S1、脑电信号的采集
S2、对脑电信号进行预处理消除噪声
首先对原始脑电信号进行128HZ重采样，128HZ重采样不仅降低了信号所占内存的大小并且有利于我们后续使用小波包分解进行能量特征的提取；使用陷波滤波器去除50HZ的工频陷波；使用0.45-35HZ的带通滤波器去除低频噪声和高频干扰；使用MATLAB自带的EEGLAB工具箱进行ICA成分分析，去除较为明显的眼电噪声成分；经过上述预处理步骤后的脑电信号用于后续的特征提取；
S3、脑电信号的特征提取
根据脑电信号随疲劳状态程度的变化特性，采用小波包变换对脑电信号进行分解，提取不同频带范围的节律被能量值和信号作为疲劳状态划分的特征值；将预处理后的脑电信号进行4层的小波包分解，共得到16个节点，
其中，与各节律波有关的节点编号及所属的频带范围下表1所示：
表1小波包节点对应脑电频段
Figure BDA0003687581740000031
根据不同节律波所处的分解子频带的系数通过公式3-1计算求取该节点的能量；
Figure BDA0003687581740000032
其中
Figure BDA0003687581740000033
表示小波包分解中第j层第i个节点的分解系数；通过上述公式可以计算出δ波、α波、θ波、β波的能量值；根据研究，
Figure BDA0003687581740000034
与疲劳状态的相关性最大；因此，最终确定了指标
Figure BDA0003687581740000035
作为划分疲劳状态的特征向量之一；
除此之外选择能反映***复杂度的非线性特征来表征大脑***的变化；最终，本次选择使用样本熵(F2)与节律波能量特征(F1)一同作为提取的特征向量进行后续疲劳状态等级的分析；
S4、采用BDPCA聚类算法对脑电信号的疲劳状态等级进行相应的划分
计算输入数据点的局部密度和聚类中心距离；根据局部密度和聚类中心距离做出簇类中心图；根据簇类中心图确定K取值范围及对应的中心点；计算各K值对应的BIC参数值；BIC参数最小值对应的聚类结果为最终结果。
进一步的，所述步骤S4的具体步骤如下所示：
(6)选取截断距离dc；
(7)按照4-1和4-2求取局部密度ρi；
Figure BDA0003687581740000036
Figure BDA0003687581740000037
(8)按式4-3求取聚类中心距离δi
Figure BDA0003687581740000041
(9)使用求取到的局部密度与聚类中心距离画图；根据所做图进行数据点类型的判定；普通点靠近横轴，异常点靠近纵轴；而聚类中心点同时远离横轴和纵轴；
(10)聚类划分：将非中心点数据点划分至局部密度大于该点中距离最近的点所属的簇中；对每个点进行判断分析后结束算法；
进一步的，使用贝叶斯信息准则(BIC)进行簇类中心点的辅助分析判定；BIC准则的核心思想为最大化当前数据集下的所选模型的概率；在针对脑电信号的研究中，
Figure BDA0003687581740000042
为脑电信号数据集，R为数据集的样本数；Mi为所选的聚类模型，Pi为模型参数的个数；假定所聚类的数据集服从k个高斯分布，因此Pi的值为MK个质心和K个方差估计之和；根据贝叶斯定理可知，在当前数据集的条件下该模型被选择的后验概率为4-4所示；
Figure BDA0003687581740000043
在初始条件下，每个模型被选择的概率相同的条件下，寻找后验概率最大的模型就是寻找边缘概率P(x1,x2,...,xnMi)最大的模型；
Figure BDA0003687581740000044
其中θi为模型的参数向量，L为最大似然函数，gi(θi)为参数向量的概率密度函数。BIC参数和最大似然函数的计算进一步简化为式4-6和4-7；
Figure BDA0003687581740000045
Figure BDA0003687581740000046
其中Rn为第n类簇的样本量，see为残差平方和，M为样本维度，R为样本容量，K为模型参数的个数。上述BIC值表征的是模型与真实数据之间的损失，BIC值越小，所选择的模型越好。
进一步的，使用相似度的概念对ρi和δi的值进行重新定义；采用共同邻居参数(CNN)对相似度的计算进行优化；优化后的计算公式如下：
Figure BDA0003687581740000047
其中CNN(sisj)表示数据点si和sj的dc邻域范围内共同存在的数据点的个数；引入该参数后，后续数据点的类型判定由距离和两点之间的局部密度共同决定；数据点将被划分至局部密度大于其的数据点中距离最近的一个，避免了由于两类不同簇的边缘点距离过近造成的误判现象。
进一步的，所述步骤S1中，噪声包括但不限于：肌电、心电等生物电噪声和工频干扰噪声。
本发明的技术效果和优点：
1、可以自动对脑电信号中表征的不同疲劳状态等级进行自动划分，客观的定义了疲劳状态等级。
2、BDPCA聚类算法对不同的疲劳状态等级进行了划分，而又无需主观判断提前对信号所属状态进行判定，是一种能够客观划分疲劳状态等级的可行的算法。。
附图说明
图1为疲劳状态等级划分机制的流程图；
图2为4层分解小波包树；
图3为波恩数据集簇类中心图；
图4为BDPCA算法工作流程图；
图5为簇类中心图；
图6为相近数据点类别误判的情况示意图；
图7为二维脑电特征分布图；
图8为波恩数据集簇类中心图；
图9为驾驶员A簇类中心图；
图10为驾驶员A特征值变化图；
图11为驾驶员A簇类中心图；
图12为驾驶员B特征值变化图；
图13为驾驶员疲劳状态等级图。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
1.疲劳状态划分机制
脑电信号中蕴含着丰富的节律波信息，当疲劳程度加深时，脑电信号中的慢波逐渐增多，快波逐渐减少，大脑***的复杂性下降。这种大脑***和脑电信号与疲劳状态的变化关系证明了通过脑电信号对疲劳状态进行划分的可能性。通过脑电信号的变化特性进行疲劳状态分析的步骤为：脑电信号的采集、聚类分析算法、脑电信号的预处理、脑电信号的特征提取、疲劳状态划分。研究者对脑电信号的特征提取方式研究通常从频域和两个方面进行着手。时域分析方法简单直观，一般提取信号的的幅值、均值、方差、正则性和同步性等参数。简单的时域特征提取方式并不适用于复杂的脑电信号分析。频域分析方法将信号投射到频率域进行分析，但是脑电信号有着较强的随机性，很多有用的信息都是瞬态的，根据“不确定性原理”，信号处理时无法同时拥有较好的时间与频率分辨率，因此在时域和频域结合起来进行处理是非常必要的。而小波包变化既兼顾了频率分辨率，又保留了时间分辨率，同时相对于小波变换，小波包变换对低频信号的划分更为细致。因此本方案使用小波包变换对采集到的的脑电信号进行特征提取。然后使用改进的BDPCA聚类算法对提取的特征向量进行训练，完成对不同疲劳状态等级的聚类区分。一个完整的疲劳状态等级划分机制的流程图如图1所示。
2.特征提取
2.1.脑电信号预处理
通过实验采集到的原始脑电信号中含有大量的诸如肌电、心电等生物电噪声和工频干扰噪声。因此在使用脑电信号进行特征提取之前需要对脑电信号进行预处理消除上述噪声。首先对原始脑电信号进行128HZ重采样，128HZ重采样不仅降低了信号所占内存的大小并且有利于我们后续使用小波包分解进行能量特征的提取。使用陷波滤波器去除50HZ的工频陷波。使用0.45-35HZ的带通滤波器去除低频噪声和高频干扰。使用MATLAB自带的EEGLAB工具箱进行ICA成分分析，去除较为明显的眼电噪声成分。经过上述预处理步骤后的脑电信号用于后续的特征提取。
2.2.脑电信号的特征提取
特征提取的主要目的是从预处理后的脑电信号中提取到与机体疲劳状态密切相关的信息用于后续的分析。脑电信号中主要包含δ波(0.5-4Hz)、θ波(4-7.5Hz)、α波(8-13Hz)、β波(14-30Hz)、γ波(30-45Hz)等五种不同频带范围的电波分量。当疲劳状态发生时，大脑的反应能力下降，同时大脑***的复杂度降低。EEG信号中γ波较多出现于病变机体的脑电信号中，因此并不适用于对疲劳状态的分析。随着疲劳程度的变化，δ波、α波、θ波、β波均会发生相应的变化，为我们对疲劳状态等级的划分提供了理论支持。根据脑电信号随疲劳状态程度的变化特性，本次研究采用小波包变换对脑电信号进行分解，提取不同频带范围的节律被能量值和信号作为疲劳状态划分的特征值。将预处理后的脑电信号进行4层的小波包分解，共得到16个节点，节点位置及编号如图2所示。
其中，与各节律波有关的节点编号及所属的频带范围如表2-1所示：
表2-1小波包节点对应脑电频段
Figure BDA0003687581740000071
根据不同节律波所处的分解子频带的系数通过公式2-1计算求取该节点的能量。
Figure BDA0003687581740000072
其中
Figure BDA0003687581740000073
表示小波包分解中第j层第i个节点的分解系数。通过上述公式可以计算出δ波、α波、θ波、β波的能量值。根据研究，
Figure BDA0003687581740000074
与疲劳状态的相关性最大。因此，最终确定了指标
Figure BDA0003687581740000081
作为划分疲劳状态的特征向量之一。
除此之外选择能反映***复杂度的非线性特征来表征大脑***的变化。最终，本次研究选择使用样本熵(F2)与节律波能量特征(F1)一同作为提取的特征向量进行后续疲劳状态等级的分析。
3.判别模型的选择
3.1.BDPCA聚类算法流程
针对有监督学习与无监督学习的算法模型的特点。本方案研究采用BDPCA聚类算法对脑电信号的疲劳状态等级进行相应的划分。BDPCA聚类算法由密度峰值聚类算法(DPCA)和贝叶斯信息准则(BIC)融合而成。DPCA聚类算法的核心思想为类簇的中心点有着更大的局部密度值，而不同类簇的中心点之间与非类簇中心点之间有着更远的距离。基于上述核心思想DPCA聚类算法的聚类步骤如下所示：
(11)选取截断距离dc。
(12)按照4-1和4-2求取局部密度ρi。
Figure BDA0003687581740000082
Figure BDA0003687581740000083
(13)按式4-3求取聚类中心距离δi
Figure BDA0003687581740000084
(14)使用求取到的局部密度与聚类中心距离画图。根据所做图进行数据点类型的判定。普通点靠近横轴，异常点靠近纵轴。而聚类中心点同时远离横轴和纵轴。
(15)聚类划分。将非中心点数据点划分至局部密度大于该点中距离最近的点所属的簇中。对每个点进行判断分析后结束算法。
DPCA算法虽然根据簇类中心图提出了簇类中心点的选择原则，但是在实际数据分析应用中，研究者大多不知道数据集的分布形式，因此很容易出现疑似簇类中心的点被选作簇类的中心点。图3为使用德国波恩大学癫痫脑电数据集提取到的特征值所画的簇类中心图。
根据图3所示，对簇类中心图的观察，可能为簇类中心点的数据点共有6个，但其实际数据集仅分为3类。因此，仅根据DPCA算法输出的簇类中心图对疑似簇类中心点的判定会产生一定的困难。因此本方案研究使用贝叶斯信息准则(BIC)进行簇类中心点的辅助分析判定。BIC准则的核心思想为最大化当前数据集下的所选模型的概率。在本方案针对脑电信号的研究中，
Figure BDA0003687581740000091
为脑电信号数据集，R为数据集的样本数。Mi为所选的聚类模型，Pi为模型参数的个数，在本次研究中，我们假定所聚类的数据集服从k个高斯分布，因此Pi的值为MK个质心和K个方差估计之和。根据贝叶斯定理可知，在当前数据集的条件下该模型被选择的后验概率为4-4所示。
Figure BDA0003687581740000092
在初始条件下，每个模型被选择的概率相同的条件下，寻找后验概率最大的模型就是寻找边缘概率P(x1,x2,...,xnMi)最大的模型。
Figure BDA0003687581740000093
其中θi为模型的参数向量，L为最大似然函数，gi(θi)为参数向量的概率密度函数。BIC参数和最大似然函数的计算进一步简化为式4-6和4-7。
Figure BDA0003687581740000094
Figure BDA0003687581740000095
其中Rn为第n类簇的样本量，see为残差平方和，M为样本维度，R为样本容量，K为模型参数的个数。上述BIC值表征的是模型与真实数据之间的损失，BIC值越小，所选择的模型越好。BDPCA聚类算法的工作流程如图4所示。
3.2.BDPCA聚类算法的改进
3.1所述的BDPCA聚类算法主要有如下不足：由于对不同的数据点计算的ρi和δi均为离散值。有可能造成不同数据点求取的局部密度ρi与聚类中心距离δi相等的情况。图5为截取的疲劳测试采集到的脑电信号所产生的簇类中心图。
观察根据提取到的特征向量生成的簇类中心图，可以发现局部密度相同的点有很多，因此局部密度最大的点所产生的疑似簇类中心点也很多，虽然可以通过计算BIC参数确定最佳簇类的个数，但会大大增加运算时间，并且会对后续的数据点所属类别的判定造成困难。
因此，本方案使用相似度的概念对ρi和δi的值进行重新定义。根据4-8和4-9计算局部密度与距离。
Figure BDA0003687581740000096
Figure BDA0003687581740000101
通过计算信号点之间的相似度代替原始的ρi和δi，极大地降低了不同数据点的ρi和δi出现相同的情况。但是，由于在计算相似度的时候单纯的使用两点的距离进行计算，因此当某一簇的边缘数据点与另一簇的边缘数据点相距较近时，容易出现的误判的情况。图6为经典双月数据集的划分结果。
从聚类结果图中可以看出聚类效果并不理想。图7为单通道脑电数据集提取到的二维特征分布图。从图中可以看出提取到的特征值数值大小十分接近，而且其分布状态在二维平面并不明显可分，因此直接使用相似度进行脑电数据点的聚类分析并不可行。
因此，为避免上述情况，本方案研究引入共同邻居参数(CNN)对相似度的计算进行优化。优化后的计算公式如下：
Figure BDA0003687581740000102
其中CNN(sisj)表示数据点si和sj的dc邻域范围内共同存在的数据点的个数。与传统相似度计算方法不同，引入该参数后，后续数据点的类型判定由距离和两点之间的局部密度共同决定。数据点将被划分至局部密度大于其的数据点中距离最近的一个，避免了由于两类不同簇的边缘点距离过近造成的误判现象。
4.BDPCA聚类算法在疲劳等级划分中的应用
BDPCA算法已经应用在一些工程的数据处理研究中，如航空发动机故障数据标记，并取得了较好的仿真结果。由于脑电信号和大多数工程信号一样具有随机性、非平稳性等特点，本方案将BDPCA算法应用在脑电数据的处理分析工作中。
4.1.BDPCA算法应用于脑电信号的可行性分析
本方案使用改进后的BDPCA算法在部分波恩癫痫脑电数据中进行聚类分析。通过聚类结果的准确率进行可行性分析。通过小波包变换提取脑电信号的δ波和θ波作为训练的特征向量。使用改进后的BDPCA算法对提取到的特征向量进行训练分析，生成的聚类中心图如图8所示。
根据簇类中心图确定疑似中心点的个数为2-6个，分别计算簇类中心点的个数为2-6个是的BIC参数。计算结果如表4-1所示。
表4-1不同K值的BIC参数
Figure BDA0003687581740000111
根据BIC参数确定最佳聚类簇数为3类。而选取的数据为健康受试者、癫痫患者未发病状态与癫痫患者的发病状态采集的脑电信号。其实际簇类个数与BIC参数的计算结果相符。BDPCA聚类算法的分析结果如表4-2所示。
表4-2波恩癫痫脑电聚类结果
Figure BDA0003687581740000112
根据聚类结果可以看出BDPCA聚类算法在对发病的癫痫患者的划分的准确率达80％。证明了BDPCA聚类算法应用在脑电信号分析中的可行性。
4.2.BDPCA算法在疲劳状态等级划分中的应用
为保证实验数据与真实的驾驶脑电数据相一致与安全性，本方案使用的数据为志愿者驾驶虚拟驾驶设备时采集到的脑电信号进行分析。为更好的查看个体的差异性与消除实验的偶然性，实验选取身体健康、视力正常、近期无服药病史且有驾照的一男一女共两名志愿者同学参加。
实验选择在易于产生疲劳状态的午后进行。正式进行实验开始前，驾驶员需要进行10min的试驾来熟悉操作过程。正式实验开始后要求驾驶员以固定的速度进行驾驶，驾驶过程中前方屏幕会随机出现障碍物，要求驾驶员操控方向盘进行闪避，整个实验过程共持续150min，整个实验期间驾驶员如感觉身体不适可随时停止，停止后实验随之结。如没有出现特殊要求，驾驶员将没有休息的持续驾驶150min，驾驶过程中，便携式EMOTIV采集设备将一直以256HZ的频率采集驾驶员的脑电信号。实验过程中还配有两名工作人员，其中一名工作人员在侧面每5min记录一次驾驶员的眨眼次数，另一名工作人员记录下障碍物出现的时间以及避障结果。上述两类参数将作为后期划分疲劳状态等级的判定标准之一。
采集到的脑电信号经过2-1的预处理步骤后，将信号划分为20s的数据段，提取数据段中的10s作为特征提取的信号。通过小波包分解提取信号的能量特征与样本熵作为算法的输入向量。分别对提取的两名驾驶员的脑电信号特征向量进行BDPCA聚类分析。
驾驶员A生成的簇类中心图如图9所示。
根据簇类中心图确定疑似中心点的个数为2-6个，分别计算簇类中心点的个数为2-6时的BIC参数。计算结果如表4-2所示。
表4-2波恩癫痫脑电聚类结果
Figure BDA0003687581740000121
根据BIC参数BDPCA算法将脑电信号分为5种不同的状态。不同状态的特征值变化如图10所示。
从图中我们可以看出能量特征值是逐渐增加的，样本熵是逐渐减小的这与疲劳程度加深时特征值理论变化趋势相一致。从图中可以观察到的能量变化更为明显，样本熵虽然有减少的趋势，但并不如能量特征值明显。
驾驶员B生成的簇类中心图如图11所示。
根据簇类中心图确定驾驶员B疑似中心点的个数也为2-6个，计算簇类中心点的个数为2-6时的BIC参数。计算结果如表4-3所示。
表4-3波恩癫痫脑电聚类结果
Figure BDA0003687581740000122
根据BIC参数BDPCA算法将脑电信号分为5种不同的状态。不同状态的特征值变化如图12所示。
根据驾驶员B的特征值变化图我们可以发现其变化趋势与驾驶员A的变化趋势相同，符合理论变化趋势。
在BDPCA聚类算法对数据进行分析训练后，使用实验过程中记录的避障结果、眨眼次数与特征值数值一同进行疲劳状态等级的划分。在通过对每一状态对应的时间段记录的眨眼次数与特征值数值的分析过后，发现驾驶员A与B的状态5时间分布较为普遍，并不集中在某一时间段，进一步观察行为信号发现，在提取该状态的10s时间段信号内，均发生过3次及以上的避障操作，因此将该状态作为多次连续动作造成的特殊状态，并不纳入疲劳状态等级中。对其余四种状态进行分析发现驾驶员B的状态1避障率要高于驾驶员A状态1的避障率，眨眼次数也远小于驾驶员A的状态1次数。继续对比发现驾驶员A的状态1的各种行为信息与驾驶员B的状态2相近，而能量特征的数值也有这种对应关系，证明驾驶员A在前10min预驾驶时已出现了疲劳状态。最终通过避障率、眨眼次数以及特征向量数值分析得到的疲劳状态等级划分如图13所示。
5.试验结果
本发明针对脑疲劳等级划分并不客观的问题，提出了一种基于BDPCA的改进聚类算法用于疲劳状态等级划分的方法。并采集不同个体的脑电信号进行验证研究。并通过避障率与眨眼次数等行为信号进行所划分的等级的疲劳状态划定。结果表明两名驾驶员相同疲劳状态等级的特征值的平均值相近，当驾驶员处在清醒与轻度疲劳状态下(F1<1.4)可以质量较高的操作模拟驾驶设备，处于中度疲劳时(F1<1.75)，驾驶质量下降较为明显，已不适合长时间驾驶，当出现重度疲劳时(F1<2.3)，驾驶质量飞速下降，眨眼次数明显变多，已不适合继续驾驶车辆，当出现极度疲劳时(F1<2.8)，驾驶员闭眼时长增加，处于近乎昏睡的状态。从结果上看，BDPCA聚类算法对不同的疲劳状态等级进行了划分，而又无需主观判断提前对信号所属状态进行判定，是一种能够客观划分疲劳状态等级的可行的算法。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (5)

1.一种基于BDPCA聚类算法的驾驶员脑电疲劳等级划分方法，其特征在于，包括以下步骤：
S1、脑电信号的采集
S2、对脑电信号进行预处理消除噪声
首先对原始脑电信号进行128HZ重采样，128HZ重采样不仅降低了信号所占内存的大小并且有利于我们后续使用小波包分解进行能量特征的提取；使用陷波滤波器去除50HZ的工频陷波；使用0.45-35HZ的带通滤波器去除低频噪声和高频干扰；使用MATLAB自带的EEGLAB工具箱进行ICA成分分析，去除较为明显的眼电噪声成分；经过上述预处理步骤后的脑电信号用于后续的特征提取；
S3、脑电信号的特征提取
根据脑电信号随疲劳状态程度的变化特性，采用小波包变换对脑电信号进行分解，提取不同频带范围的节律被能量值和信号作为疲劳状态划分的特征值；将预处理后的脑电信号进行4层的小波包分解，共得到16个节点，
其中，与各节律波有关的节点编号及所属的频带范围下表1所示：
表1小波包节点对应脑电频段
Figure FDA0003687581730000011
根据不同节律波所处的分解子频带的系数通过公式3-1计算求取该节点的能量；
Figure FDA0003687581730000021
其中
Figure FDA0003687581730000022
表示小波包分解中第j层第i个节点的分解系数；通过上述公式可以计算出δ波、α波、θ波、β波的能量值；根据研究，
Figure FDA0003687581730000023
与疲劳状态的相关性最大；因此，最终确定了指标
Figure FDA0003687581730000024
作为划分疲劳状态的特征向量之一；
除此之外选择能反映***复杂度的非线性特征来表征大脑***的变化；最终，本次选择使用样本熵(F2)与节律波能量特征(F1)一同作为提取的特征向量进行后续疲劳状态等级的分析；
S4、采用BDPCA聚类算法对脑电信号的疲劳状态等级进行相应的划分
计算输入数据点的局部密度和聚类中心距离；根据局部密度和聚类中心距离做出簇类中心图；根据簇类中心图确定K取值范围及对应的中心点；计算各K值对应的BIC参数值；BIC参数最小值对应的聚类结果为最终结果。
2.根据权利要求1所述的一种基于BDPCA聚类算法的驾驶员脑电疲劳等级划分方法，其特征在于：所述步骤S4的具体步骤如下所示：
(1)选取截断距离dc；
(2)按照4-1和4-2求取局部密度ρi；
Figure FDA0003687581730000025
Figure FDA0003687581730000026
(3)按式4-3求取聚类中心距离δi
Figure FDA0003687581730000027
(4)使用求取到的局部密度与聚类中心距离画图；根据所做图进行数据点类型的判定；普通点靠近横轴，异常点靠近纵轴；而聚类中心点同时远离横轴和纵轴；
(5)聚类划分：将非中心点数据点划分至局部密度大于该点中距离最近的点所属的簇中；对每个点进行判断分析后结束算法。
3.根据权利要求1所述的一种基于BDPCA聚类算法的驾驶员脑电疲劳等级划分方法，其特征在于：使用贝叶斯信息准则(BIC)进行簇类中心点的辅助分析判定；BIC准则的核心思想为最大化当前数据集下的所选模型的概率；在针对脑电信号的研究中，
Figure FDA0003687581730000028
为脑电信号数据集，R为数据集的样本数；Mi为所选的聚类模型，Pi为模型参数的个数；假定所聚类的数据集服从k个高斯分布，因此Pi的值为MK个质心和K个方差估计之和；根据贝叶斯定理可知，在当前数据集的条件下该模型被选择的后验概率为4-4所示；
Figure FDA0003687581730000031
在初始条件下，每个模型被选择的概率相同的条件下，寻找后验概率最大的模型就是寻找边缘概率P(x1,x2,...,xnMi)最大的模型；
Figure FDA0003687581730000032
其中θi为模型的参数向量，L为最大似然函数，gi(θi)为参数向量的概率密度函数。BIC参数和最大似然函数的计算进一步简化为式4-6和4-7；
Figure FDA0003687581730000033
Figure FDA0003687581730000034
其中Rn为第n类簇的样本量，see为残差平方和，M为样本维度，R为样本容量，K为模型参数的个数。上述BIC值表征的是模型与真实数据之间的损失，BIC值越小，所选择的模型越好。
4.根据权利要求1所述的一种基于BDPCA聚类算法的驾驶员脑电疲劳等级划分方法，其特征在于：使用相似度的概念对ρi和δi的值进行重新定义；采用共同邻居参数(CNN)对相似度的计算进行优化；优化后的计算公式如下：
Figure FDA0003687581730000035
其中CNN(sisj)表示数据点si和sj的dc邻域范围内共同存在的数据点的个数；引入该参数后，后续数据点的类型判定由距离和两点之间的局部密度共同决定；数据点将被划分至局部密度大于其的数据点中距离最近的一个，避免了由于两类不同簇的边缘点距离过近造成的误判现象。
5.根据权利要求1所述的一种基于BDPCA聚类算法的驾驶员脑电疲劳等级划分方法，其特征在于：所述步骤S1中，噪声包括但不限于：肌电、心电等生物电噪声和工频干扰噪声。
CN202210650658.1A 2022-06-10 2022-06-10 一种基于bdpca聚类算法的驾驶员脑电疲劳等级划分方法 Pending CN114916937A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210650658.1A CN114916937A (zh) 2022-06-10 2022-06-10 一种基于bdpca聚类算法的驾驶员脑电疲劳等级划分方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210650658.1A CN114916937A (zh) 2022-06-10 2022-06-10 一种基于bdpca聚类算法的驾驶员脑电疲劳等级划分方法

Publications (1)

Publication Number Publication Date
CN114916937A true CN114916937A (zh) 2022-08-19

Family

ID=82812815

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210650658.1A Pending CN114916937A (zh) 2022-06-10 2022-06-10 一种基于bdpca聚类算法的驾驶员脑电疲劳等级划分方法

Country Status (1)

Country Link
CN (1) CN114916937A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116807478A (zh) * 2023-06-27 2023-09-29 常州大学 一种驾驶员瞌睡开始状态检测方法、装置及设备

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116807478A (zh) * 2023-06-27 2023-09-29 常州大学 一种驾驶员瞌睡开始状态检测方法、装置及设备

Similar Documents

Publication Publication Date Title
CN104586387B (zh) 一种时、频、空域多参数脑电特征提取与融合方法
CN110772268A (zh) 一种多模脑电信号及1dcnn迁移的驾驶疲劳状态识别方法
Picot et al. On-line automatic detection of driver drowsiness using a single electroencephalographic channel
Cheng et al. A motor imagery EEG feature extraction method based on energy principal component analysis and deep belief networks
CN110811649A (zh) 一种基于生物电及行为特征融合的疲劳驾驶检测方法
CN111523601A (zh) 一种基于知识引导和生成对抗学习的潜在情绪识别方法
CN111951637B (zh) 一种任务情景相关联的无人机飞行员视觉注意力分配模式提取方法
CN113158964B (zh) 一种基于残差学习和多粒度特征融合的睡眠分期方法
CN105320969A (zh) 基于多尺度Renyi熵的心率变异性特征分类方法
CN110390272B (zh) 一种基于加权主成分分析的eeg信号特征降维方法
Peng et al. 3D-STCNN: Spatiotemporal Convolutional Neural Network based on EEG 3D features for detecting driving fatigue
CN114916937A (zh) 一种基于bdpca聚类算法的驾驶员脑电疲劳等级划分方法
Zheng et al. EEG emotion recognition based on hierarchy graph convolution network
Singh et al. Physical and physiological drowsiness detection methods
Guan Application of logistic regression algorithm in the diagnosis of expression disorder in Parkinson's disease
Lv et al. Compact vehicle driver fatigue recognition technology based on EEG signal
CN113397482B (zh) 一种人类行为分析方法及***
Chai et al. Classification of driver fatigue in an electroencephalography-based countermeasure system with source separation module
Dey et al. Automatic detection of drowsiness in EEG records based on time analysis
Murugappan et al. Subtractive fuzzy classifier based driver drowsiness levels classification using EEG
Ma et al. Research on drowsy-driving monitoring and warning system based on multi-feature comprehensive evaluation
CN112617761A (zh) 自适应聚点生成的睡眠阶段分期方法
Jantan et al. A multi-model analysis for driving fatigue detection using EEG signals
Adlakha Single trial EEG classification
Zhang et al. EEG and Forehead EOG-Based Driver Fatigue Classification Using Sparse-Deep Belief Networks

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication