一种基于脑电信号的疲劳状态识别方法

1.一种基于脑电信号的疲劳状态识别方法，其特征在于，方法包括：
（1）通过脑电信号数据采集仪从人体头皮表面记录不同疲劳状态下的脑电信号，通过信号采集模拟电路滤出信号中的干扰成分，再将微弱的电信号程控放大并消除电平飘移，得到的模拟脑电信号经数字电路进行模数转换成数字脑电信号传输给主机进行处理；
（2）主机首先对信号进行预处理，然后提取每通道信号的特征信息，构造特征向量，最后采用基于模糊模式识别的评估方法，根据所得到的脑电特征评估疲劳程度；对所采集的脑电信号进行处理的方法如下：
a、信号预处理
对采集的脑电数据进行预处理，即进行数字滤波，提取0-30Hz的信号频率成分；
b、信号特征提取
对采集的脑电数据进行特征提取，以四种脑电基本节律的相对能量作为特征，它们是：
δ ：0-3.5Hz、θ ：3.5-8Hz、α ：8-13Hz、β ：13-30Hz，由于δ、θ、α、β 节律相对能量之和为1，只选取δ、θ、α 三个节律相对能量作为特征信息；
c、模糊模式识别
精神疲劳程度是一个模糊概念，记为A，设有需要对模糊概念A 进行识别的n 个样本组成的集合，有m 个指标特征值表示样本的整体特性，则有样本集的指标特征值矩阵
(1)
式中 表示样本j 指标i 的特征值，i=1,2, … ,m; j=1,2 … ,n；
设样本集依据m个指标按c 个级别的已知指标标准特征值进行识别，则有指标标准特征值矩阵
(2)
式中 表示状态或级别h指标i标准特征值，h=1，2，…，c； i=1，2，…，m ；
若规定1级精神状态为不疲劳状态，其评价指标i 的标准特征值 对于疲劳程度的相对隶属度 ；c级精神状态为非常疲劳状态，其评价指标 对于疲劳程度的相对隶属度 ，则人体精神疲劳状态标准指标 的相对隶属度 可根据如下线性公式确定
(3)
由此可将矩阵Y 变换为相应的相对隶属度矩阵S
(4)
可根据下式
(5)
将指标值越大疲劳程度越深的样本指标值 变为相应的相对隶属度 ；相反，对于指标值越小疲劳程度越深的指标用下式进行变换；
(6)
由此可将矩阵X 变换为相应的相对隶属度矩阵
(7)
将样本j的m个指标相对隶属度 ， ，… ， 分别与矩阵S的第1，第2，…，第m 行的行向量逐一进行比较，可得样本j 的级别上限值 和级别下限值 （1≤ ≤c，1≤≤c）；
设样本集对模糊概念A 各个级别的相对隶属度矩阵为
(8)
式中 表示样本j 对A 级别h 的相对隶属度，j=1，2，…，n ； h=1，2，…，c ；
由于样本j 在级别区间 ， 范围内，故矩阵U 应满足归一化约束条件
(9)
样本j 与级别h 之间的差异用广义欧式距离表示为
(10)
为了完善地描述样本j 与级别h 之间的差异，定义下式
(11)
以 为权重，称为加权广义欧式距离；
为了求解样本j 对模糊概念A 的级别h 的最优相对隶属度，建立目标函数，
(12)
根据目标函数和约束条件构造拉格朗日函数，求得最优相对隶属矩阵，矩阵中的各元素由下式确定
(13)
最后应当用级别特征值F1,F2,…,Fc，对疲劳程度进行判断，级别特征值公式表示为：
(14)
H 是 包 含n 个 元 素 的 向 量， 为 样 本j 的 级 别 特 征 值，如 果
，则判断样本j隶属于级别k。

一种基于脑电信号的疲劳状态识别方法\n技术领域\n[0001] 本发明属于生物医学、信号处理和模式识别技术领域，具体涉及一种基于脑电信号的疲劳状态识别方法，通过基于脑电信号的检测与识别实现对人体疲劳状态的评估。\n背景技术\n[0002] 人在经过一段时间工作，特别是长时间的精神或体力劳动后，会产生精神疲劳，引起生理与心理状态的变化，包括继续努力时难度的增加和信息分析能力下降趋势的增强，行为组织能力的下降以及心情的变化等。在现代社会中，特别是诸如实时监控、交通运输、高危作业、航空航天等作业情况下，由于疲劳状态下警觉性、持续性注意力、工作记忆力，判断力、决策制定等能力的下降，工作人员非常容易出现随意操作和违章行为，从而引发安全事故。因此，建立起客观、可靠，并且对正常工作无干扰的精神疲劳检测方法和评估系统，对于防范由于精神疲劳引发的安全事故是十分必要的，具有极大的经济和社会价值。\n[0003] 精神疲劳检测评估方法通常分为主观方法和客观方法两类。目前主要采用主观的方法进行评估，根据受试者的身体、神经感受状态来判断疲劳的程度。然而，这些方法易受主观因素的影响，不能客观的进行定量评估，无法反映疲劳程度的变化情况，并且评分标准也不易统一。现阶段研究的客观方法主要是利用基于生理特征分析的方法加以评估，这些生理特征主要包括通过脑电信号（electroencephalogram, EEG）检测的脑活动，通过心电信号检测的心率变异性，通过肌电信号检测的皮肤电反应，通过动态图像观测到的眼动、头动、面部表情等。在这些生理信息来源中，EEG被认为是最有前途、最可预报的和最可靠的。\nEEG能够客观的记录大脑机能状态的连续变化，是人脑思维活动、认知和意识状态的一种外在表现。利用EEG来研究人的精神疲劳程度，采用各种高效的信号处理技术挖掘EEG所蕴涵的丰富信息，实现对疲劳状态的判定与评估，是非常重要而又有前途的途径。\n[0004] 近年来，EEG检测技术的进步和现代信号处理技术的发展，推动了基于EEG的精神疲劳研究。目前利用脑电信号分析精神疲劳主要有两类研究方法：一是利用事件相关电位（event related potential, ERP）,二是利用自发脑电信号。ERP是一种诱发脑电，需要外加特定的刺激才能产生，因此在汽车驾驶、高危作业、航空航天等作业情况下，是不可能同期采集ERP信号的，这就极大的限制了该类方式的应用范围。而自发脑电是大脑皮层神经电活动本身自发产生的持续性的电位变化，可以在作业的同期采集，对作业无干扰。\n发明内容\n[0005] 本发明的目的是提出一种基于脑电信号的疲劳状态识别方法。精神疲劳程度是一个模糊概念，不同疲劳状态之间的边界是不分明的，应用模糊模式识别来判定精神疲劳程度相比传统的分类方法更合理，能提高对精神疲劳的分类识别正确率。\n[0006] 本发明为一种基于自发脑电信号（以下简称EEG）的疲劳状态识别方法，其特征为，该方法通过EEG数据采集仪从人体头皮表面记录不同疲劳状态下的脑电信号，然后对所采集的信号进行信号放大、滤波、转换为数字信号，并且通过无线通信传输给主机，主机软件处理接收到的无线脑电数据，首先对信号进行预处理，然后提取每通道信号的特征信息，构造特征向量，最后采用基于模糊模式识别的评估方法，根据所得到的脑电特征评估疲劳程度。\n[0007] 本发明所提出的精神疲劳评估方法和系统对EEG进行处理的方法如下： [0008] 1、EEG采集\n[0009] EEG的幅度在微伏量级，并且信号中含有很强的干扰。在本方案中，首先对模拟脑电信号进行预处理，包括对信号进行前置放大、高通滤波、低通滤波、后级放大及电平迁移，然后进行A/D采样，输入到FPGA芯片中进行信号处理，包括50 Hz陷波和去基线漂移，然后经无线发射器发送出去。数字脑电信号经无线接收器接收后由微处理器进行处理，再通过USB接口芯片上传PC机进行下一步的处理。\n[0010] 2、信号特征提取\n[0011] PC机中对采集的脑电数据进行特征提取，以delta（δ）：0-3.5Hz、theta（θ）：\n3.5-8Hz、alpha（α）：8-13Hz、beta（β）：13-30Hz四种脑电基本节律的相对能量作为特征，因为δ、θ、α、β 节律相对能量之和为1，为避免冗余，只选取δ、θ、α 三个节律相对能量作为特征信息。\n[0012] 3、模糊模式识别\n[0013] 精神疲劳程度是一个模糊概念，记为A，设有需要对模糊概念A 进行识别的n 个样本组成的集合，有m 个指标特征值表示样本的整体特性，则有样本集的指标特征值矩阵[0014] (1)\n[0015] 式中 表示样本j 指标i 的特征值，i=1,2, … ,m; j=1,2 … ,n。\n[0016] 如样本集依据m个指标按c 个级别的已知指标标准特征值进行识别，则有指标标准特征值矩阵\n[0017] (2)\n[0018] 式中 表示状态或级别h指标i标准特征值，h=1，2，…，c； i=1，2，…，m。\n[0019] 若规定1级精神状态为不疲劳状态，其评价指标i 的标准特征值 对于疲劳程度的相对隶属度 ；c 级精神状态为非常疲劳状态，其评价指标 对于疲劳程度的相对隶属度 ，则人体精神疲劳状态标准指标 的相对隶属度 可根据如下线性公式确定\n[0020] \n(3)\n[0021] 由此可将矩阵Y 变换为相应的相对隶属度矩阵S\n[0022] (4)\n[0023] 类似地，可根据下式\n[0024] (5)\n[0025] 将指标值越大疲劳程度越深的样本指标值 变为相应的相对隶属度 。相反，对于指标值越小疲劳程度越深的指标用下式进行变换。\n[0026] (6)\n[0027] 由此可将矩阵X 变换为相应的相对隶属度矩阵\n[0028] (7)\n[0029] 将样本j 的m 个指标相对隶属度 ， ，… ， 分别与矩阵S的第1，第2，…，第m 行的行向量逐一进行比较，可得样本j的级别上限值 和级别下限值 （1≤ ≤c，\n1≤ ≤c）。\n[0030] 设样本集对模糊概念A 各个级别的相对隶属度矩阵为\n[0031] (8)\n[0032] 式中 表示样本j对A级别h的相对隶属度，j=1，2，…，n； h=1，2，…，c。\n[0033] 由于样本j 在级别区间 ， 范围内，故矩阵U 应满足归一化约束条件[0034] (9)\n[0035] 样本j 与级别h 之间的差异用广义欧式距离表示为\n[0036] (10)\n[0037] 为了完善地描述样本j 与级别h 之间的差异，定义下式\n[0038] (11)\n[0039] 以 为权重，称为加权广义欧式距离。\n[0040] 为了求解样本j 对模糊概念A 的级别h 的最优相对隶属度，建立目标函数，[0041] (12)\n[0042] 根据目标函数和约束条件构造拉格朗日函数，求得最优相对隶属矩阵，矩阵中的各元素由下式确定\n[0043] (13)\n[0044] 最后应当用级别特征值F1,F2,…,Fc，对疲劳程度进行判断，级别特征值公式表示为：\n[0045] \n(14)\n[0046] H 是 包 含n 个 元 素 的 向 量， 为 样 本j 的 级 别 特 征 值，如 果，则判断样本j隶属于级别m。\n[0047] 根据上述方法可以构建基于脑电信号的精神疲劳评估系统，该系统包括信号采集和主处理器，信号采集部分通过信号采集模拟电路滤出信号中的干扰成分，再将微弱的电信号程控放大并消除电平飘移，得到的模拟脑电信号经数字电路进行AD转换成数字脑电信号，然后经无线发射器发送出去；主处理器部分通过无线接收器接收到信号后再通过USB接口芯片上传PC机，在PC机中将完成对EEG的预处理和特征提取，并判断待识别EEG所对应的疲劳程度。\n[0048] 本发明结合了生物医学信号处理技术及模糊模式识别技术，提出了一种客观的可行的精神疲劳评估方法，使基于脑电信号的检测与识别在精神疲劳评估领域的应用得到技术上的极大提升。\n附图说明\n[0049] 图1：是实现本发明方法的精神疲劳评估系统的方框图；\n[0050] 图2：是本发明方法的精神疲劳评估算法框图；\n[0051] 图3：是实现本发明方法的模糊模式识别算法程序流程图。\n具体实施方式\n[0052] 以下结合附图说明该方法及系统的具体实施方式：\n[0053] 图1表示实现本发明方法的精神疲劳评估系统的方框图。利用这个系统可以实现对被检测对象进行精神疲劳评估并输出评估结果。\n[0054] 如图1所示，该系统包括了EEG采集模拟电路、采集数字电路、无线发射和接收电路、微处理器、USB接口芯片和PC机等六大部分。信号采集部分通过电极检测到受试者头皮表面的原始脑电信号，脑电信号采集模拟电路滤出干扰信号并将人脑上的微弱的电信号程控放大并消除电平飘移，得到的模拟脑电信号经数字电路进行AD转换成数字脑电信号，然后经无线发射器发送出去，数字脑电信号经无线接收器接收后由微处理器进行处理，再通过USB接口芯片上传到PC机上。在PC机上首先提取各通道脑电信号进行预处理，提取频带特征信息，然后对这些特征信息进行分类识别，评估受试者所处的精神疲劳状态。\n[0055] 系统中的PC机对采集到的EEG信号进行的信号处理具体步骤如下： [0056] （1）对采集的脑电数据进行预处理，即进行数字滤波，提取0-30Hz的信号频率成分，然后进行特征提取：计算δ、θ、α和β四个节律的能量Eδ，Eθ，Eα，Eβ，计算总能量和相对能量 。因为δ、θ、α、β 节律相对\n能量之和为1，为避免冗余，只选取δ、θ、α 三个节律相对能量作为各通道EEG信号的特征信息。\n[0057] （2）通过模糊模式识别算法判断从（1）中提取的特征所对应的疲劳状态。下面以不疲劳，中度疲劳，非常疲劳三级疲劳状态分级为例，对本发明作进一步的说明。设精神疲劳程度分为3个级别：级别1——不疲劳，级别2——中度疲劳，级别3——非常疲劳，这三个级别的级别特征值F1，F2，F3分别设为：0，0.5，1。级别1对疲劳的相对隶属度为0，对疲劳的相对隶属度随着级别特征值的增大而升高，级别3对疲劳的相对隶属度为1。\n[0058] ①在训练阶段，得到三级疲劳状态标准特征值矩阵。其具体步骤为：记录受试者在这三种状态下的三导（Fp1、Fp2、Fpz）脑电数据，每种状态数据N 组，三种状态共计3×N 组脑电数据。以每导EEG数据的δ、θ、α 三个节律相对能量pδ，pθ，pα为特征，三导脑电数据一共有9个特征值，分别对每种状态各组数据的9个特征值进行平均，作为该状态的9个标准指标，可得三级疲劳状态标准特征值矩阵为\n[0059] \n[0060] ②在测试阶段，记录受试者在某种状态下的三导（Fp1、Fp2、Fpz）脑电数据，提取测试数据的特征值组成样本指标矩阵 ，根据公式（4）至（7）可将矩阵X与\nY变换为相应的相对隶属度矩阵R 和S。\n[0061] ，\n[0062] 将矩阵R，S 中数据，带入模型（13）得到该组测试数据对于各级疲劳程度的相对隶属度矩阵 。最后应用公式（14）计算得到该组测试样本的级别特征值向\n量：\n[0063] \n[0064] 根据级别特征值向量H的值判断测试样本所隶属的疲劳级别。