基于脑电波的便携式睡眠质量监测系统

基于脑电波的 便携式睡眠质量监测系统 金旭扬 导师： 华东理工大学信息学院 万永菁 上海中学信息学科组 吴奕明 1 摘要 睡眠是人体重要的生理活动，睡眠质量近年来受到高度关注；本文从脑电波角度探寻睡 眠监测的有效易行方法，从软硬件角度设计了便携式睡眠质量监测系统。研究分析便携式脑 电采集设备采集的数据和 CAP 睡眠脑电数据库，用功率谱分析和 BP 神经网络探究了睡眠分 期的有效算法。实验进行了初步的睡眠分期与质量评估，证明了便携式睡眠质量监测系统的 准确性及利用脑电数据进行睡眠分期的有效性。本课题研究，提出了利用单导连脑电信号进 行睡眠分期的可行性，为之后研究便携式、市场化的睡眠监测设备以及其他应用提供了重要 的实验参考依据。 关键词：脑电；脑机接口；睡眠监测；睡眠分期；BP 神经网络 一、引言 1.1 睡眠质量研究背景及意义 睡眠是一种重要的生理现象。从生到死，人类始终是在觉醒和睡眠中度过。人类通过高 质量的睡眠，可以消除疲劳，更好地恢复精神和体力，使人在睡眠之后保持良好的觉醒状态， 提高工作、学习效率。 人类用于睡眠的时间占人一生中的三分之一。然而迄今我们对这一重要的生理现象的认 识还微乎其微，对睡眠进行科学的研究只有短短的几十年历史。1937 年，Lomis、Harvey 和 Hobart 注意到，睡眠不是处于一种稳定状态，而是要发生一系列非常有规律的周期性变 化。 1 1986 年，Rechtschaffen 等人重新肯定了 Dement 和 Kleitman 的分期标准，并根据十年 来的经验作了一些必要的修改和补充，使之更趋完善。 2 2007 年，美国睡眠医学会基于上述标准进行改进，发布了新的睡眠分期专业标准，其 中规定了各个指标具体的采集标准及判定方法。 3 1.2 脑电信号分析方法综述 随着电子技术的发展，数字处理技术逐步应用到 EEG 的分析中来。经典的 EEG 分析方法 有：以分析 EEG 波形的几何性质，如幅度、均值、峭度等为主的时域分析方法和以分析 EEG 各频率功率、相干等为主的领域方法。早在 70 年代初，W.C.Yeo 和 J.P.Smith4就应用 Walsh 谱分析离线地研究了一个处于睡眠状态的男性的三段脑电图。R.D.Larsen 等 5应用 Walsh 顺序的 Walsh 函数对 EEG 进行展开，并定义了双值自相关函数，尔后讨论了可以按双 值自相关函数来显示各种睡眠 EEG 的特征。 1982 年，美国物理学家 Hopfield 提出了 HNN 模型，从而有力地推动了应用神经网络方 法解释许多复杂生命过程的进展。自八十年代末以来，人工神经网络的应用已涉及到了脑电 分析的各个方面，其中包括自发脑电的睡眠分级及睡眠 EEG 分析。S.Roberts 和 L.Tarassenko6,7把人工神经网络应用于睡眠 EEG 的自动分析。他们采用无监督学习网络对 大量没有经过人工判别的数据进行自组织分类，少量的经过人工判别的标准样本则用来自组 织分类结果做解释和量化，从而在网络中形成了 8 个聚类区。根据 EEG 在 8 个聚类区之间随 时间运动的轨迹可以对一夜的睡眠状况有定性的了解。 8 1.3 脑电监测设备介绍 目前，脑电监测设备大致有二： 2 一为大型的、医院专用的多导睡眠监测系统。这种系统需要测量多导连的脑电图、眼电 图、肌电图、口鼻气流、呼吸运动、血氧饱和度等众多指标，且有严格的判定规则、需要专 业知识。 3 二为便携式脑机接口设备。此类设备通常体积小、使用方便、成本也较低，测量的脑电 图多为单导连，但由于获取的数据用途较为单一，可以很好地完成睡眠监测的任务。 9 1.4 课题研究目标 本课题利用便携式脑电波采集设备实时获取脑电数据，并且与终端设备通讯实时存储、 分析数据。利用 Windows、Android 等移动平台下编写的软件实现此功能，实现人体的睡眠 监控。 二、方法和假设 2.1 系统软硬件平台的基本架构 2.1.1 睡眠质量监测系统的硬件组成 用于采集数据的设备是宏智力公司出品的 Brainlink 意念力头箍，它采用基于 Neurosky 芯片平台的 Thinkgear 芯片，主要用于检测脑电信号。实验采用手机（Android）系统和电 脑（Windows）系统作为采集终端。 图 2-1 睡眠质量监测系统框图 图 2-2 宏智力公司出品的 Brainlink 意念力头箍 2.1.2 睡眠质量监测系统的数据采集方式 NeuroSky 的脑电波采集设备较为轻便，只有前额、左耳垂两个电极（一导连）。设备采 用 AAA 电池供电，根据介绍续航能力有 8 小时（若再并联一颗电池可以更长），没有传统脑 电采集中与脑电频段接近的 50Hz 工频交流干扰信号。设备采用无线蓝牙连接，更有利于睡 眠时数据的传输；耳垂采用导电夹，容易固定；利用心电图电极片改装前额电极，也可以弥 补原本接触不良的缺点。 为了完成单向传输数据的目的，使用的蓝牙模拟串口（发送）芯片能耗低、续航能力强、 编程较为简易。接收端可以是任何蓝牙 4.0 设备，只需一次配对后就可自动连接，对于手机、 电脑硬件的要求不高。 3 初步测试时，采集使用的是 Microsoft Windows 平台，使用 Neurosky 提供的 API 接口， 在 Visual C+上编写简单的程序即可完成数据的存盘。采样频率约为 513.5Hz，远高于脑电 信号的最高有效频率 30Hz 的两倍，符合采样定理。 图 2-3 Windows 7 下的采集、分析软件 利用 Neurosky 提供的 Android API 接口，在 Android 平台下的脑电波预览、采集工作 也得以完成，程序可以在后台运行，并且将采样数据即使存盘，在实际使用过程中更为方便， 也省去了用电脑建立连接、定义接口的繁杂步骤，适合移动平台。 图 2-4 Android 下的采集、预览软件 2.2 基于脑电信号的睡眠质量监测方法 2.2.1 脑电信号预处理方法 脑电波在时域上属于非平稳随机信号，实验中采集的脑电波只有一导连，因此信号不稳 定、噪波严重。需要经过初步的低通数字滤波预处理。为方便起见，频率衰减带上限取到高 于脑电波分析中有效频率 30Hz 的 50Hz。 数字滤波器包括有限冲激响应(FIR)滤波器和无限冲激响应(IIR)滤波器两大类。FIR 滤 波器可以得到严格的线性相位，相比 IIR 需要采用较高的阶数（约是 IIR 的五至十倍），但 软件实现方便。 10 假设 FIR 滤波器的单位冲击响应 h(n)为一个长度为 N 的序列，那么滤波器的系统函数为： (2-1) 4 上式的差分形式为： (2-2) 由于理想滤波器在边界频率处不连续，故其时域信号 hd(n)一定是无限时宽的，无法实 现。因此，需要把具有理想线性相位特性的滤波器曲线用窗函数截取： (2-3) 这种设计思想称为窗函数设计法。 其中，常用的汉明窗（Hamming Window）函数如下： (2-4) 幅值函数为： (2-5)11 使用 Matlab 的 fir1 工具设计 300 点的 FIR 低通滤波器，采用汉明窗，以 512Hz 作为采 样频率，50Hz 作为率减带，得到的滤波器幅频响应曲线如下： 图 2-5 300 点低通滤波器幅频响应曲线 （采用归一化角频率，2 即为实际采样频率的 512Hz） 2.2.2 脑电信号的频域分析 脑电波按频率从高到低划分依次为： 波(1430Hz)， 波(814Hz)， 波(48Hz)， 波(0.54Hz)。 脑电波 频率范围 精神状态 波 1430Hz 运动感觉节律，放松可集中注意力，有协调性，思 考，对于自我和周围环境意识清楚机警，激动 波 814Hz 放松但不困倦，安静，有意识 波 48Hz 直觉的，回忆的，幻想，想象，浅睡 5 波 0.54Hz 深度睡眠，非快动眼睡眠，无意识 表 2-1 脑电波的频段划分以及不同类型脑电波所反映出的脑部精神状态 12 离散时间序列 x(n)的傅立叶（Fourier Transform）变换是： (2-6) 如已知随机信号 x(n)的自相关函数 r(k)，那么功率谱密度函数就定义为： (2-7) 功率谱函数的另一定义是： (2-8) 理论上，离散信号处理方法对有限带宽的信号能做准确分析，但有限带宽信号在时域上 是无限长的，只取其中有限长的一段进行傅立叶变换，相当于在原信号上加了矩形窗运算。 加窗在频域上，对原功率谱起到了平滑的作用。 (2-9)13 其中，w(n)表示窗口函数。常用的窗有三角窗、汉宁窗、汉明窗、布莱克曼窗等。这些 窗的旁瓣电平比矩形窗低，但分辨率也较矩形窗低。 14这里使用的是汉明窗，在 2.2.1 节 已有详细介绍。 2.2.3 睡眠分期判定的改进算法 根据美国睡眠医学会 2007 年的标准，睡眠分期的脑电标准如下： 睡眠阶段 划分规则（仅含脑电，且忽略例外情况） W（觉醒期） 枕区 波含量大于 50% N1（非快速眼动期 1） 波减弱，低幅度、47Hz 的波含量大于 50% N2（非快速眼动期 2） 开始阶段：出现与觉醒无关的 K 复合波或者纺锤波 持续阶段：低幅度、47Hz 的波（不含 K 复合波或纺 锤波） N3（非快速眼动期 3） 0.52Hz 的慢波含量大于 20% R（快速眼动期） 出现低幅度、47Hz 的波（不含 K 复合波或纺锤波） 表 2-2 睡眠分期的脑电标准 3 其中，非快速眼动期睡眠深度从深到浅，依次是：N3、N2、N1。 人工神经网络（Artificial Neural Network）是由大量简单的处理单元广泛连接组成 的复杂网络，用于 6 模拟人类大脑神经网络的结构和行为。它反映了人脑功能的许多基本特性，但它并不是 人脑全部的真实写照，而只是对其作某种简化、抽象和模拟 15。 在各种学习算法中，多层网络的反向传播算法（简称 BP 算法）应用最为广泛。BP 算法 最早是由 Werbos 在 1974 年提出来的， Rumelhart 等人于 1985 年发展了该理论，提出了清 晰而又严格的算法。BP 算法适用于前向网络，它采用有导师学习的训练形式，提供输入矢 量集的同时提供输出矢量集，通过反向传播学习算法，调整网络的连接权值，以使网络输出 在最小均方差意义下，尽量向期望输出接近，反向学习的进程由正向传播和反向传播组成。 在正向传播过程中，输入信息经隐含神经元逐层处理并传向输出层，如果输出层不能得到期 望的输出，则转入反向传播过程，将实际输出与期望输出之间的误差沿原来的连接通路返回， 通过修改各层神经元的连接权值，使误差减小，然后转入正向传播过程，反复循环，直至误 差小于给定的值为止。 设有 N 个训练对组成的训练集，每一个训练对用输入矢量 Xi= (xi1, xi2, , xim)和输 出矢量 Di=(di1, di2, , din)， 1 i N。在前向传播中，把 Xi作为网络的输入，根据 现有的 W 计算网络的输出 Yi=(yi1, yi2, , yin)。比较实际输出 Yi与期望输出 Di之间的差 异，计算每一个输出单元的平方误差(y ij- dij)2， 1 j n。把这些误差进行加总得到误 差函数： (2-10) 所要做的就是通过改变 W 来减小 E，以使得所有的输入矢量都尽量与相应的输出矢量相 匹配。因此学习的过程就转化为定义在权值空间上的目标函数 E 的极小化问题。 在训练过程中总是以尽可能快的减小 E 的方式进行。一般它依赖于在权值空间中是否沿 梯度方向搜索，所以采用梯度下降法来训练权值。每一个权值 wij的变化量w ij按如下方式 计算： (2-11) 其中 Z 为学习率，是控制算法收敛速度的参数。 在第一阶段得到的总误差平方和又在第二阶段被一层一层地反向传播回去，从输出单元 到输入单元。权值的调整决定于传播过程中的每一步。由于 Ii、f i和 E 都是连续可微的，因 此，可以应用以下公式计算E / w ij的值： (2-12) W 的修改可以有两种方式，一是对于每一训练对(X i, Di)都修改一次 W，另一种方式是输 入全部的训练对后再加总w ij并进行修改。训练矢量集中训练对的数目称为一个 epoch。 当 epoch 不是非常大的时候，后一种方式能够加快收敛的速度。因为第一种方式只能针对某 一特定的训练对减小误差函数，而可能增大其它训练对的误差函数；第二种方式总是以减小 总体误差函数为目标的。所以采用第二种方式 16。 图 2-6 BP 神经网络结构图 三、数据记录与分析 7 3.1 系统各部分效果验证 3.1.1 滤波器效果验证 为了验证实时滤波的效果以及硬件性能，在清醒时采集的脑电波中截取了一段约 4 秒的 信号进行验证。经过检验，滤波器效果良好，经过实时滤波的数据已经符合脑电波分析的要 求。 图 3-1 滤波前后的脑电信号波形对比（采样率 512Hz） 图 3-2 滤波前后的脑电信号频谱对比（采样率 512Hz） 3.1.2 加窗频域分析效果验证 由于睡眠深时低频率的脑电波所占功率比例会增强，睡眠浅时会减弱，因此使用单一频 段的波所占功率比例可以简单判读睡眠的深浅程度。 实验对象佩戴脑电采集设备一晚上约九小时的睡眠初步验证，每 30 秒数据、加汉明窗 频域分析了 波（0.54Hz）所占的功率比例，得到了如下的图像。经过比对与参考，该图 像已能大致反映睡眠的深浅程度，符合人体睡眠周期的客观规律。加窗频域分析的效果，得 以验证。 图 3-3 加窗频域分析后所得 波所占的功率比例 3.2 睡眠分期判定方法 8 3.2.1 睡眠深浅的目测方法 实验对象佩戴脑电采集设备，未服用任何辅助药物或干预治疗，进行了连续五晚的睡眠 脑电波采集。受试者身体健康，入睡时间正常且有规律，睡眠周期较为完整。由于睡眠监测 实验的“第一晚效应”，即受试者在第一次佩戴设备时难以入睡、或是睡眠期间易惊醒、没 有完整的睡眠周期的情况，因此只采用了第二晚至第五晚的数据进行分析。采样频率固定为 512Hz，每晚的睡眠连续时长均超过 6 小时。 多次取连续十分钟（600 秒）片段，每 30 秒加汉明窗频域分析得到了各频段的脑电波所 占的功率比例。目测发现，有如下规律： 波和 波频率的峰或谷出现的位置大都重合， 即一种波形的频率处于峰值时，另一种波形的频率处于谷值。在经过数十次的非连续采集片 段分析后，确认了本次实验中上述规律的普遍性。 图 3-4 波和 波频率比例的峰或谷出现的位置重合 3.2.2 BP 神经网络分析 由于条件的限制，难以用本次实验所用的便携式设备获取大量、准确的脑电波数据，也 难以得到由医生给出的专业睡眠分期判断。这部分研究采用从 PhysioNet17获得的 CAP 睡眠 脑电数据库 18进行分析和算法验证。该睡眠监测实验在意大利帕尔马的 Ospedale Maggiore 睡眠障碍研究中心进行。此数据库有 108 例多导睡眠记录，每例至少记录了三导连的脑电信 号（根据 10-20 国际通用系统，电极为： F3 或 F4、C3 或 C4、O1 或 O2，以 A1 或 A2 作为参 考电极）。其中的 16 例由健康的成年受试者完成，这 16 名受试者无神经系统疾病，未使用 会影响中枢神经的药物。受试者为 9 名女性、7 名男性，年龄从 23 岁至 42 岁不等。此外， 在睡眠中心接受过训练的神经病学家，还根据 Rechtschaffen 293297. 5Larsen R.D., et al. Math Biosci 1976; 31: 237253. 6Roberts S, et al. Med Biol Eng Comput 1992; 30: 509517. 7Roberts S, et al. IEEE procedings-F 1992; 139(6): 420425. 8张杰, 王明时. 睡眠脑电的研究. 国外医学生物医学工程分册, 1997 年第 20 卷第 2 期. 9SHAMBROOM, J. R., FBREGAS, S. E. and JOHNSTONE, J. (2012), Validation of an automated wireless system to monitor sleep in healthy adults. Journal of Sleep Research, 21: 221230. doi: 10.1111/j.1365-2869.2011.00944.x. 10铙志强, 叶念渝. FIR 和 IIR 数字滤波器的探讨与实现. 计算机与数字工程, 2005, 33(7). 11史洁玉. MATLAB 信号处理超级学习手册. 人民邮电出版社, 2014. 9. 12陈群. 脑电生物反馈治疗仪与睡眠分期的研究. 广东工业大学, 2013. 13李颖洁, 邱意弘, 朱贻盛. 脑电信号分析方法及其应用. 北京: 科学出版社, 2009. 14Petre Stoica, Randolph L. Moses. Spectral Analysis of Signals. Pearson Prentice Hall, 2005. 15施鸿宝. 神经网络及其应用. 西安: