《基于脑电频域特征指标化算法测量警觉性的帽子及方法.pdf》由会员分享，可在线阅读，更多相关《基于脑电频域特征指标化算法测量警觉性的帽子及方法.pdf（9页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 103815900 A (43)申请公布日 2014.05.28 CN 103815900 A (21)申请号 201310594232.X (22)申请日 2013.11.22 A61B 5/0476(2006.01) A61B 5/16(2006.01) (71)申请人 刘志勇 地址 100086 北京市海淀区中关村南大街 2 号数码大厦 A 座 509 室 (72)发明人 刘志勇 (74)专利代理机构 北京中济纬天专利代理有限 公司 11429 代理人 张晓霞 (54) 发明名称 基于脑电频域特征指标化算法测量警觉性的 帽子及方法 (57) 摘要 本发明涉及一种。
2、基于脑电频域特征指标化算 法测量警觉性的帽子及方法。所述帽子还包括脑 电信号采集装置、 IC 电路、 信号传输装置及终端 设备， 所述采集装置对人脑活动产生的模拟脑电 波进行采样、 量化， 变成离散数字信号， 进行后续 的处理 ； 所述信号放大器将量化的数字信号进行 放大， 增强处理和传输过程中的抗干扰能力 ； 所 述信号处理器对采集的信号进行去噪预处理， 增 强脑电波信号的强度， 并从中提取反映人认知状 态变化的计算指标， 评估使用者的状态 ； 所述信 号传输装置连接 IC 电路与终端设， 将 IC 电路提 取的计算指标传输给终端设备 ； 所述终端设备为 一 PC 机， 对上述计算指标进行处。
3、理， 并进行显示 和反馈。本发明测量精度高可以应用于各种工作 环境。 (51)Int.Cl. 权利要求书 2 页 说明书 5 页 附图 1 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页 说明书5页 附图1页 (10)申请公布号 CN 103815900 A CN 103815900 A 1/2 页 2 1. 一种基于脑电频域特征指标化算法测量警觉性的帽子， 所述帽子还包括脑电信号采 集装置、 IC 电路、 信号传输装置及终端设备， 其特征在于 ： 脑电信号采集装置包括脑电信号采集电极、 脑电信号参考电极和信号处理器， 所述采 集装置对人脑活动产生的模拟脑电波进。
4、行采样、 量化， 变成离散数字信号， 进行后续的处 理 ； 所述 IC 电路包括信号放大器及信号处理器， 所述信号放大器为前置信号放大器， 将量 化的数字信号进行放大， 增强处理和传输过程中的抗干扰能力 ； 所述信号处理器对采集的信号进行去噪预处理， 增强脑电波信号的强度， 并从中提取 反映人认知状态变化的计算指标， 评估使用者的状态 ； 所述信号传输装置连接 IC 电路与终端设， 将 IC 电路提取的计算指标传输给终端设 备 ； 所述终端设备为一 PC 机， 对上述计算指标进行处理， 并进行显示和反馈。 2. 根据权利要求 1 所述的设备， 其特征在于 ： 所述 IC 电路内置于帽子内部， 。
5、进一步整 合有电池仓及电源开关。 3. 根据权利要求 1 所述的设备， 其特征在于 ： 所述脑电信号采集电极位于头部， 脑电信 号参考电极夹位于耳部。 4. 根据权利要求 1 所述的设备， 其特征在于 ： 所述信号传输装置为无线信号传输装置。 5. 一种基于脑电频域特征指标化算法的测量警觉性的帽子， 其特征在于， 所述方法， 所 述方法步骤如下 ： （1） 初始化设备硬件， 设定可接受条件 ； （2） 采集脑电数据 ； （3） 对采集的数据滤波、 去噪、 时频分析、 计算脑电指标 ； （4） 将计算的脑电指标传输至终端并进行显示 ； （5） 判断是否符合接受条件， 进而判断使用者是否处于疲劳状。
6、态， 如果判定处于疲劳状 态， 则给出警示 ； 如果判定不处于疲劳状态， 则再次至步骤 （2） 进行相关步骤。 6. 根据权利要求 5 所述的方法， 其特征在于， 所述指标化算法具体如下 ： （1） 预处理 ： 对量化的脑电波进行数字滤波， 去除肌电等干扰噪音 ； 所述滤波器为无限 冲激响应 （IIR） 带通滤波器 ； （2） 特征表达和提取 ： 从预处理之后的脑电信号中提取反映综合认知状态的基本指标， 具体包括 alpha 波 （8-13Hz） 、 beta 波 （13-20Hz） ， delta（1-4Hz） ， theta（4-7Hz） ； 运用时 频分析技术将这些指标从原始时域信号中提。
7、取出来， 在频域上以能量或功率的时间序列形 式加以表达 ； （3） 指标化表示 ： 将上述基本指标进行标准化， 使得不同使用者以及同一个使用者不同 时间的同一指标具有相同的物理含义 ； 所述算法输出警觉性水平和紧张度水平两个指标， 所述警觉性水平指标和和紧张度水平指标具体如下 ： 权 利 要 求 书 CN 103815900 A 2 2/2 页 3 a(t)、 b(t) 和 c(t) 分别表示 alpha、 beta 和 theta 三个频段的时序信号， 他们分别由 原始脑电信号的经由时频分析并选定特定频段的能量累加实现 ； 警觉度指标 ： S1(t)=c(t)/a(t)， 其中 t 表示时间。
8、， a 和 c 分别表示 alpha 和 theta 的能量 ； 紧张度指标 ： S2(t)=b(t)*c(t)， 其中 t 表示时间， b 和 c 分别表示 beta 和 theta 的能量 ； （4） 注意力水平的判定 ： 即以正常使用者在没有困倦、 疲劳状态发生的情况下， 持续保 持注意2分钟的警觉性水平指标和和紧张度水平指标两项指标序列的时序均值的60%作为 判决门限， 低于此门限判决为疲劳发生。 7. 根据权利要求 6 所述的算法， 其特征在于 ： 所述步骤 （1） 中， 带通滤波器的低通起始 频率为 1Hz， 高通截止频率为 35Hz。 8. 根据权利要求 5 所述的算法， 其特征。
9、在于 ： 所述步骤 （2） 特征表达子模块的具体算法如下 ： 采用 Morlat 函数为母小波函数， 对脑电时域信号进行连续小波变换 ； 输入信号为单 导联脑电采集到的离散脑电时间序列， 与以上母小波函数卷积并经过变换后得到一系列不 同频段的复时间序列信号即小波系数， 其中时间与输入信号长度保持一致， 频率范围对为 1-35Hz， 保留其中 1-35Hz 用于提取节律特征 ； 对于特定的时刻和频点， 系数表示信号的时 频分布情况， 对其取模， 即用功率表示 ； 按照频段范围高低， 分别提取 alpha 波 （8-13Hz） 、 beta 波 （13-20Hz） ， delta（1-4Hz） ，。
10、 theta（4-7Hz） 对应频段的复小波系数的模的时间序 列， 即功率表示频段能量的时序变化。 9. 根据权利要求 5 所述的算法， 其特征在于 ： 所述步骤 （3） 指标化表示采用特征归一化算法， 即 ： 将某一频段能量占总能量的比值作为指标 ： 其中， t 表示时间， f 表示频率， P 表示功率， 从而 Pf(t) 表示某一频率 f 范围内的能量 随时间变化的函数， 公式的分母部分则表示将 1 到 35Hz 频率范围内的能量累加总和 ； 根据 以上模型将 Pf(t) 除以总能量归一化以后， 每个频段的相对能量变为 0-1 范围内的数值， 即 用 Pf(t) 表示。 权 利 要 求 书。
11、 CN 103815900 A 3 1/5 页 4 基于脑电频域特征指标化算法测量警觉性的帽子及方法 0001 技术领域 0002 本发明涉及一种基于脑电频域特征指标化算法测量警觉性的帽子及方法。 0003 背景技术 0004 脑电可以反映人的注意力或疲劳与否早已为众人所熟知， 并且实验室环境中相关 的技术也非常的成熟。在许多场合了解被试是否为觉醒态具有重要意义， 如教室和工厂车 间等， 尤其是驾驶环境中， 当使用者处于困倦、 注意力不集中状态时马上可以报警或给出提 示。但是由于脑电是非常微弱的电信号 （微伏级） ， 在实际环境中记录的过程很容易受到干 扰， 能够提取稳定的表征疲劳和困倦的指标。
12、并把它非常方便的应用与实际生活或工作中的 技术， 没有查阅到文献或专利报道。 0005 发明内容 0006 本发明的目的， 提供一种基于脑电频域特征指标化算法测量警觉性的帽子及方 法， 以客观的判断、 评估测试者的状态并及时的给出提示。 0007 本发明的技术方案如下 ： 一种基于脑电频域特征指标化算法测量警觉性的帽子， 所述帽子还包括脑电信号采集 装置、 IC 电路、 信号传输装置及终端设备， 其特征在于 ： 脑电信号采集装置包括脑电信号采集电极、 脑电信号参考电极和信号处理器， 所述采 集装置对人脑活动产生的模拟脑电波进行采样、 量化， 变成离散数字信号， 进行后续的处 理 ； 所述 IC。
13、 电路包括信号放大器及信号处理器， 所述信号放大器为前置信号放大器， 将量 化的数字信号进行放大， 增强处理和传输过程中的抗干扰能力 ； 所述信号处理器对采集的信号进行去噪预处理， 增强脑电波信号的强度， 并从中提取 反映人认知状态变化的计算指标， 评估使用者的状态 ； 所述信号传输装置连接 IC 电路与终端设， 将 IC 电路提取的计算指标传输给终端设 备 ； 所述终端设备为一 PC 机， 对上述计算指标进行处理， 并进行显示和反馈。 0008 进一步的， 所述 IC 电路内置于帽子内部， 进一步整合有电池仓及电源开关。 0009 进一步的， 所述脑电信号采集电极位于头部， 脑电信号参考电极。
14、夹位于耳部。 0010 进一步的， 所述信号传输装置为无线信号传输装置。 0011 一种基于脑电频域特征指标化算法的测量警觉性的帽子， 其特征在于， 所述方法， 所述方法步骤如下 ： 说 明 书 CN 103815900 A 4 2/5 页 5 （1） 初始化设备硬件， 设定可接受条件 ； （2） 采集脑电数据 ； （3） 对采集的数据滤波、 去噪、 时频分析、 计算脑电指标 ； （4） 将计算的脑电指标传输至终端并进行显示 ； （5） 判断是否符合接受条件， 进而判断使用者是否处于疲劳状态， 如果判定处于疲劳状 态， 则给出警示 ； 如果判定不处于疲劳状态， 则再次至步骤 （2） 进行相关步。
15、骤。 0012 6、 根据权利要求 5 所述的方法， 其特征在于， 所述指标化算法具体如下 ： （1） 预处理 ： 对量化的脑电波进行数字滤波， 去除肌电等干扰噪音 ； 所述滤波器为无限 冲激响应 （IIR） 带通滤波器 ； （2） 特征表达和提取 ： 从预处理之后的脑电信号中提取反映综合认知状态的基本指标， 具体包括 alpha 波 （8-13Hz） 、 beta 波 （13-20Hz） ， delta（1-4Hz） ， theta（4-7Hz） ； 运用时 频分析技术将这些指标从原始时域信号中提取出来， 在频域上以能量或功率的时间序列形 式加以表达 ； （3） 指标化表示 ： 将上述基本指。
16、标进行标准化， 使得不同使用者以及同一个使用者不同 时间的同一指标具有相同的物理含义 ； 所述算法输出警觉性水平和紧张度水平两个指标， 所述警觉性水平指标和和紧张度水平指标具体如下 ： a(t)、 b(t) 和 c(t) 分别表示 alpha、 beta 和 theta 三个频段的时序信号， 他们分别由 原始脑电信号的经由时频分析并选定特定频段的能量累加实现 ； 警觉度指标 ： S1(t)=c(t)/a(t)， 其中 t 表示时间， a 和 c 分别表示 alpha 和 theta 的能量 ； 紧张度指标 ： S2(t)=b(t)*c(t)， 其中 t 表示时间， b 和 c 分别表示 bet。
17、a 和 theta 的能量 ； （4） 注意力水平的判定 ： 即以正常使用者在没有困倦、 疲劳状态发生的情况下， 持续保 持注意2分钟的警觉性水平指标和和紧张度水平指标两项指标序列的时序均值的60%作为 判决门限， 低于此门限判决为疲劳发生。 0013 进一步的， 所述步骤 （1） 中， 带通滤波器的低通起始频率为 1Hz， 高通截止频率为 35Hz。 0014 进一步的， 所述步骤 （2） 特征表达子模块的具体算法如下 ： 采用 Morlat 函数为母小波函数， 对脑电时域信号进行连续小波变换 ； 输入信号为单 导联脑电采集到的离散脑电时间序列， 与以上母小波函数卷积并经过变换后得到一系列不。
18、 同频段的复时间序列信号即小波系数， 其中时间与输入信号长度保持一致， 频率范围对为 1-35Hz， 保留其中 1-35Hz 用于提取节律特征 ； 对于特定的时刻和频点， 系数表示信号的时 频分布情况， 对其取模， 即用功率表示 ； 按照频段范围高低， 分别提取 alpha 波 （8-13Hz） 、 beta 波 （13-20Hz） ， delta（1-4Hz） ， theta（4-7Hz） 对应频段的复小波系数的模的时间序 列， 即功率表示频段能量的时序变化。 0015 进一步的， 所述步骤 （3） 指标化表示采用特征归一化算法， 即 ： 将某一频段能量占总能量的比值作为指标 ： 说 明 书。
19、 CN 103815900 A 5 3/5 页 6 其中， t 表示时间， f 表示频率， P 表示功率， 从而 Pf(t) 表示某一频率 f 范围内的能量 随时间变化的函数， 公式的分母部分则表示将 1 到 35Hz 频率范围内的能量累加总和 ； 根据 以上模型将 Pf(t) 除以总能量归一化以后， 每个频段的相对能量变为 0-1 范围内的数值， 即 用 Pf(t) 表示。 0016 本发明的有益效果在于 ： 本发明基于脑电频域特征指标化算法测量警觉性， 测量精度高， 能精确的对脑电信号 进行提取和处理， 并且基于脑电频域特征指标化算法， 可以客观的判断使用者当前的精神 状态， 并给出及时的。
20、提示。可以应用于各种工作环境。 0017 附图说明 0018 图 1 是本发明的设备结构示意图。 0019 图 2 是本发明的方法流程示意图。 0020 其中， 1- 脑电信号采集电极， 2- IC 电路， 3- 脑电信号参考电极， 4- 帽子。 0021 具体实施方式 0022 如图 1 所示， 是本发明的帽子结构示意图， 包括帽子、 脑电信号采集装置、 IC 电路、 信号传输装置及终端设备， 脑电信号采集装置包括脑电信号采集电极和脑电信号参考电 极， 所述采集装置对人脑活动产生的脑电波进行采样、 量化， 变成离散数字信号， 进行后续 的处理 ； 该两个电极与帽子相整合， 为了更为精确的采集。
21、脑电信号， 脑电信号采集电极位于 头部， 脑电信号参考电极夹位于耳部。 0023 IC 电路包括信号放大器及信号处理器， 所述信号放大器为前置信号放大器， 将量 化的数字信号进行放大， 增强处理和传输过程中的抗干扰能力 ； IC 电路也是内置于帽子 中。 信号处理器对采集的信号进行去噪预处理， 增强脑电波信号的强度， 并从中提取反映人 认知状态变化的计算指标， 评估使用者的状态 ； 信号传输装置连接 IC 电路与终端设备， 将 IC 电路提取的计算指标传输给终端设备。帽子还进一步整合有电池仓及电源开关。 0024 终端设备为一 PC 机， 对上述计算指标进行处理， 并进行显示和反馈。一般来说，。
22、 终 端设备与 IC 电路之间采用无线信号传输。 0025 本发明的一种基于脑电频域特征指标化算法的测量警觉性的帽子， 所述方法步骤 如下 ： （1） 初始化设备硬件， 设定报警条件 ； （2） 同时采集脑电数据， 一般采集 1 秒 ； （3） 对采集的数据滤波、 去噪、 时频分析、 计算脑电指标 ； （4） 将计算的脑电指标传输至终端并进行显示 ； （5） 判断是否符合报警条件， 进而判断使用者是否处于疲劳状态， 如果判定处于疲劳状 态， 则给出警示 ； 如果判定不处于疲劳状态， 则再次至步骤 （2） 进行相关步骤。 0026 其中， 指标化算法具体如下 ： （1） 预处理 ： 对量化的脑电。
23、波进行数字滤波， 去除肌电等干扰噪音 ； 所述滤波器为无限 说 明 书 CN 103815900 A 6 4/5 页 7 冲激响应 （IIR） 带通滤波器 ； 带通滤波器的低通起始频率为 1Hz， 高通截止频率为 35Hz。 0027 （2） 特征表达和提取 ： 从预处理之后的脑电信号中提取反映综合认知状态的基本 指标， 具体包括 alpha 波 （8-13Hz） 、 beta 波 （13-20Hz） ， delta（1-4Hz） ， theta（4-7Hz） ； 运 用时频分析技术将这些指标从原始时域信号中提取出来， 在频域上以能量或功率的时间序 列形式加以表达 ； （3） 指标化表示 ： 。
24、将上述基本指标进行标准化， 使得不同使用者以及同一个使用者不同 时间的同一指标具有相同的物理含义 ； 所述算法输出警觉性水平和紧张度水平两个指标， 所述警觉性水平指标和和紧张度水平指标具体如下 ： a(t)、 b(t) 和 c(t) 分别表示 alpha、 beta 和 theta 三个频段的时序信号， 他们分别由 原始脑电信号的经由时频分析并选定特定频段的能量累加实现 ； 警觉度指标 ： S1(t)=c(t)/a(t)， 其中 t 表示时间， a 和 c 分别表示 alpha 和 theta 的能量 ； 紧张度指标 ： S2(t)=b(t)*c(t)， 其中 t 表示时间， b 和 c 分别。
25、表示 beta 和 theta 的能量 ； （4） 注意力水平的判定 ： 即以正常使用者在没有困倦、 疲劳状态发生的情况下， 持续保 持注意2分钟的警觉性水平指标和和紧张度水平指标两项指标序列的时序均值的60%作为 判决门限， 低于此门限判决为疲劳发生。 0028 其中， 特征表达子模块的具体算法如下 ： 采用 Morlat 函数为母小波函数， 对脑电时域信号进行连续小波变换 ； 输入信号为单 导联脑电采集到的离散脑电时间序列， 与以上母小波函数卷积并经过变换后得到一系列不 同频段的复时间序列信号即小波系数， 其中时间与输入信号长度保持一致， 频率范围对为 1-35Hz， 保留其中 1-35H。
26、z 用于提取节律特征 ； 对于特定的时刻和频点， 系数表示信号的时 频分布情况， 对其取模， 即用功率表示 ； 按照频段范围高低， 分别提取 alpha 波 （8-13Hz） 、 beta 波 （13-20Hz） ， delta（1-4Hz） ， theta（4-7Hz） 对应频段的复小波系数的模的时间序 列， 即功率表示频段能量的时序变化。 0029 其中， 所述步骤 （3） 指标化表示采用特征归一化算法， 即 ： 将某一频段能量占总能量的比值作为指标 ： 其中， t 表示时间， f 表示频率， P 表示功率， 从而 Pf(t) 表示某一频率 f 范围内的能量 随时间变化的函数， 公式的分母。
27、部分则表示将 1 到 35Hz 频率范围内的能量累加总和 ； 根据 以上模型将 Pf(t) 除以总能量归一化以后， 每个频段的相对能量变为 0-1 范围内的数值， 即 用 Pf(t) 表示。 0030 大脑协同程度指标表示 ： 在以上两项基本指标的基础上， 我们又根据不同节律表示的内在生理、 心理规律， 提出 利用节律信号之间的同步性来综合反映使用者的大脑协同状态， 从而全面表达大脑活动的 状态是否适合工作。 0031 指标计算流程如下 ： 从原始信号中分别提取出 8-13Hz 和 4-7Hz 频段范围的 alpha 和 theta 波， 用 a(t) 和 说 明 书 CN 103815900。
28、 A 7 5/5 页 8 c(t) 表示， 其中 t 表示时间。 0032 对脑电幅度信号 a(t) 和 b(t) 分别进行 Hilbert 变换， 得到其相位信号 a(t) 和 b(t)， 表示信号相位随时间变化的情况 ； 计算 alpha 和 theta 能量时间序列之间的同步性指标 S， 表示一段时间内 （用表示） ， 两 个频段信号总体的相位差值， 即同步性的好坏， 用于衡量全脑参与维持注意水平的程度， 同 步性越好， 越能调动全脑的认知资源维持较高的注意水平， 能够保证警觉性， 客服疲劳等因 素， 提高对外界刺激的加工能力， 从而保持较好的工作状态。 0033 指标 S 的计算模型如。
29、下 ： 其中， S 表示拟计算的同步性指标， 其中表示选定的一段时间长度， 信号将根据该时间 段的长度从原始信号中逐步计算出， e 表示自然常数， 其值约为 2.71828 ； t 表示选定时间段 内的某一时刻 ； (t) 表示节律的相位信息 ； 对一段时间内的相位的差进行累加可以计算 出总体的相位同步性， 以自然对数形式表示可以保证指标在 0 和 1 之间。 0034 本专利中， 选定时间段长度为1s， 即每1s输出以上指标S一次， 以实时地跟踪注意 力指标的变化， 传输给终端并加以提示。 0035 流程图可以按以下流程进行 ： 预处理 - 特征表达和提取 - 特征归一化 - 特征指标化表示 Alpha Theta Beta 前两个合成 A ： 警觉度 ； 后两个合成 B ： 紧张度 ； 第一个和第三个深度合成 C ： 协同性指标 ； 然后指向门限判别和指标输出。 说 明 书 CN 103815900 A 8 1/1 页 9 图 1 图 2 说 明 书 附 图 CN 103815900 A 9 。