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

1、实用标准文档基于脑电波的便携式睡眠质量监测系统金旭扬导师：华东理工大学信息学院万永菁上海中学信息学科组吴奕明摘要睡眠是人体重要的生理活动，睡眠质量近年来受到高度尖注；本文从脑电波角度探寻睡眠监测的有 效易行方法，从软硬件角度设计了便携式睡眠质量监测系统。研究分析便携式脑电采集设备采集的数据 和CAP睡眠脑电数据库，用功率谱分析和BP神经网络探究了睡眠分 期的有效算法。实验进行了初步 的睡眠分期与质量评估，证明了便携式睡眠质量监测系统的准确性及利用脑电数据进行睡眠分期的有效 性。本课题研究，提出了利用单导连脑电信号进行睡眠分期的可行性，为之后研究便携式、市场化的睡 眠监测设备以及其他应用提供了重

2、要的实验参考依据。尖键词：脑电；脑机接口 ；睡眠监测；睡眠分期；BP神经网络文案大全实用标准文档一、引言1.1睡眠质量研究背景及意义睡眠是一种重要的生理现象。从生到死，人类始终是在觉醒和睡眠中度过。人类通过高 质量的睡 眠，可以消除疲劳，更好地恢复精神和体力，使人在睡眠之后保持良好的觉醒状态，提高工作、学习效 率。人类用于睡眠的时间占人一生中的三分之一。然而迄今我们对这一重要的生理现象的认 识还微乎其 微，对睡眠进行科学的研究只有短短的几十年历史。1937年，Lomis、Harvey和Hobart注意到，睡眠不是处于一种稳定状态，而是要发生一系列非常有规律的周期性变化。in1986年，Rech

3、tschaffe n等人重新肯定了 Deme nt和Kleitman的分期标准，并根据十年 来的经 验作了一些必要的修改和补充，使之更趋完善。2007年，美国睡眠医学会基于上述标准进行改进，发布了新的睡眠分期专业标准，其中规定了各 个指标具体的采集标准及判定方法。1.2脑电信号分析方法综述随着电子技术的发展，数字处理技术逐步应用到EEG的分析中来。经典的EEG分析方法有：以分析EEG波形的几何性质，如幅度、均值、哨度等为主的时域分析方法和以分析EEG各频率功率、相干等为主的领域方法。早在70年代初，W.C.Yeo和J.P.Smith就应用Walsh谱分析离线地研究了一个处于睡眠状态的男性的三段

4、脑电图。R.D.Larsen等(5)应用Walsh顺序的Walsh函数对EEG进行展开，并定义了双值自相矢函数，尔后讨论了可以按双 值自相矢 函数来显示各种睡眠EEG的特征。1982年，美国物理学家Hopfield提出了 HNN模型，从而有力地推动了应用神经网络方法解释许 多复杂生命过程的进展。自八十年代末以来，人工神经网络的应用已涉及到了脑电分析的各个方面，其 中包括自发脑电的睡眠分级及睡眠EEG分析o S. Roberts和L.Tarassenko阴】把人工神经网络应用于 睡眠EEG的自动分析。他们采用无监督学习网络对大量没有经过人工判别的数据进行自组织分类，少量的经过人工判别的标准样本则

5、用来自组织分类结果 做解释和量化，从而在网络中形成了8个聚类区。根据EEG在8个聚类区之间随时间运动的轨迹可以对一夜的睡眠状况有定性的了解。(6)1.3脑电监测设备介绍目前,脑电监测设备大致有二：文案大全实用标准文档一为大型的、医院专用的多导睡眠监测系统。这种系统需要测量多导连的脑电图、眼电图、肌电 图、口鼻气流、呼吸运动、血氧饱和度等众多指标，且有严格的判定规则、需要专业知识。（8）二为便携式脑机接口设备。此类设备通常体积小、使用方便、成本也较低，测量的脑电图多为单导 连，但由于获取的数据用途较为单一，可以很好地完成睡眠监测的任务。（8）1.4课题研究目标本课题利用便携式脑电波采集设备实时获

6、取脑电数据，并且与终端设备通讯实时存储、分析数据。 利用Win dows An droid等移动平台下编写的软件实现此功能，实现人体的睡眠监控。二、方法和假设2.1系统软硬件平台的基本架构2.1.1睡眠质量监测系统的硬件组成用于采集数据的设备是宏智力公司出品的Brai nli nk意念力头箍，它采用基于Neurosky芯片平台 的Thi nkgear芯片，主要用于检测脑电信号。实验采用手机（An droid ）系统和电脑（Windows系 统作为采集终端。广卩也斷幣井图2-2宏智力公司出品的Brainlink意念力头箍2.1.2睡眠质量监测系统的数据采集方式NeuroSky的脑电波采集设备较为

7、轻便，只有前额、左耳垂两个电极（一导连）。设备采用AAA电 池供电，根据介绍续航能力有8小时（若再并联一颗电池可以更长），没有传统脑电采集中与脑电频段 接近的50Hz工频交流干扰信号。设备采用无线蓝牙连接，更有利于睡眠时数据的传输；耳垂采用导电 夹，容易固定；利用心电图电极片改装前额电极，也可以弥补原本接触不良的缺点。为了完成单向传输数据的目的，使用的蓝牙模拟串口（发送）芯片能耗低、续航能力强、编程较为 简易。接收端可以是任何蓝牙4.0设备，只需一次配对后就可自动连接，对于手机、电脑硬件的要求不 高。图2-1睡眠质量监测系统框图文案大全信号的最咼有效频率30Hz的两倍，符合采样定理In什 H

8、“Zd T WE? QMOcignin实用标准文档初步测试时，采集使用的是Microsoft Windows平台，使用Neurosky提供的API接口，在Visual C+上编写简单的程序即可完成数据的存盘。采样频率约为513.5Hz,远高于脑电和Q 睜S UM BI CCM51*1 CMSelect Set图2-3 Windows 7下的采集、分析软件利用Neurosky提供的An droid API接口，在An droid平台下的脑电波预览、采集工作也得以完 成，程序可以在后台运行，并且将采样数据即使存盘，在实际使用过程中更为方便，也省去了用电脑建 立连接、定义接口的繁杂步骤，适合移动平台

9、。图2-4 An droid下的采集、预览软件2.2基于脑电信号的睡眠质量监测方法2.2.1脑电信号预处理方法脑电波在时域上属于非平稳随机信号，实验中采集的脑电波只有一导连，因此信号不稳定、噪波严重。 需要经过初步的低通数字滤波预处理。为方便起见，频率衰减带上限取到高于脑电波分析中有效频率30Hz 的 50Hz。数字滤波器包括有限冲激响应(FIR)滤波器和无限冲激响应(IIR)滤波器两大类。FIR滤波器可以得到 严格的线性相位，相比HR需要采用较高的阶数(约是HR的五至十倍)，但软件实现方便。10假设FIR滤波器的单位冲击响应h(n)为一个长度为N的序列，那么滤波器的系统函数为：Logging

10、 Now 1 $ 53A7SET FILENAMEH=OH(z)=龙&何州(2-1)文案大全实用标准文档上式的差分形式为：文案大全实用标准文档y (门丿二工瓜恥)巩痒一 m I(2-2)由于理想滤波器在边界频率处不连续，故其时域信号hd(n)定是无限时宽的，无法实现。因此，需要把具有理想线性相位特性的滤波器曲线用窗函数截取：:杜掳二选字阮帧曲(2-3)这种设计思想称为窗函数设计法。其中，常用的汉明窗(Hamming WindoW函数如下：Wm() = 0.54 一 0.46 咖JRN)幅值函数为：H血3 = (创叭3)+仏23你3右寸)+(畑心锂I 4(2-5)使用Matlab的fir1工具设

11、计300点的FIR (ii 低通滤波器，采用汉明窗，以512Hz作为采样频率，50Hz作为率减带，得到的滤波器幅频响应曲线如下：隔伽I血世AmiueM( (:44*1图2-5 300点低通滤波器幅频响应曲线(采用归一化角频率，2 n即为实际采样频率的512Hz)2.2.2脑电信号的频域分析脑电波按频率从高到低划分依次为：B波(1430Hz)，口波(814Hz)波(48Hz)波(0.54Hz)脑电波频率范围精神状态B波14-30HZ运动感觉节律，放松可集中注意力，有协调性，思考，对于 自我和周围环境意识清楚机警，激动a波814Hz放松但不困倦，安静，有意识B波48Hz直觉的，回忆的，幻想，想象，

12、浅睡文案大全实用标准文档s波0.54Hz深度睡眠，非快动眼睡眠，无意识表2-1脑电波的频段划分以及不同类型脑电波所反映出的脑部精神状态11离散时间序列x(n)的傅立叶(FourierTransform )变换是：二工(对-扣粹n=J(2-6)如已知随机信号x(n)的自相矢函数r(k)，那么功率谱密度函数就定义为：I+00尸)二工 rk)ek(2-7)功率谱函数的另一定义是：亠心1丨(2-8)理论上，离散信号处理方法对有限带宽的信号能做准确分析，但有限带宽信号在时域上是无限长 的，只取其中有限长的一段进行傅立叶变换，相当于在原信号上加了矩形窗运算。加窗在频域上，对原 功率谱起到了平滑的作用。11

13、 NF(w)二可刀皿5)工(切5A n=l(2-9)【口其中,w(n)表示窗口函数。常用的窗有三角窗、汉宁窗、汉明窗、布莱克曼窗等。这些窗的旁瓣电 平比矩形窗低，但分辨率也较矩形窗低。阿这里使用的是汉明窗，在2.2.1节已有详细介绍。其中，非快速眼动期睡眠深度从深到浅，依次是：N3 N2、N1。2.2.3睡眠分期判定的改进算法 根据美国睡眠医学会2007年的标准，睡眠分期的脑电标准如下：睡眠阶段划分规则(仅含脑电，且忽略例外情况)W(觉醒期)枕区a波含量大于50%N1 (非快速眼动期1)a波减弱，低幅度、47Hz的波含量大于50%N2 (非快速眼动期2)开始阶段：出现与觉醒无尖的K复合波或者纺

14、锤波持续阶段： 低幅度、47Hz的波(不含K复合波或纺锤波)N3(非快速眼动期3)0.52Hz的慢波含量大于20%R (快速眼动期)出现低幅度、47Hz的波(不含K复合波或纺锤波)表22睡眠分期的脑电标准文案大全w的变化量厶Wj按如下方式(2-11)(2-12)实用标准文档人工神经网络（Artificial Neural Network ）是由大量简单的处理单元广泛连接组成的复杂网络，用于 模拟人类大脑神经网络的结构和行为。它反映了人脑功能的许多基本特，性，但它并不是人脑全部的 真实写照，而只是对其作某种简化、抽象和模拟问o在各种学习算法中，多层网络的反向传播算法（简称BP算法）应用最为广泛。

15、BP算法最早是由Werbos在1974年提出来的，Rumelhart等人于1985年发展了该理论，提出了清 晰而又严 格的算法。BP算法适用于前向网络，它采用有导师学习的训练形式，提供输入矢量集的同时提供输出矢 量集，通过反向传播学习算法，调整网络的连接权值，以使网络输出在最小均方差意义下，尽量向期望 输出接近，反向学习的进程由正向传播和反向传播组成。在正向传播过程中，输入信息经隐含神经元逐 层处理并传向输出层，如果输出层不能得到期望的输出，则转入反向传播过程，将实际输出与期望输出 之间的误差沿原来的连接通路返回，通过修改各层神经元的连接权值，使误差减小，然后转入正向传播 过程，反复循环，直至

16、误差小于给定的值为止。设有N个训练对组成的训练集，每一个训练对用输入矢量X= （X h, x已，x研和输出矢量D=（dii, d已，d in） 1 i N o在前向传播中，把X作为网络的输入，根据现有 的W计算网络的畅出丫=（y u, y已，y in）比较实际输出丫与期望输出D之间的差异，计算每一个输也单元阳平方误筆（,X 2）钿、卜、口主、形-亦斗曰巩【、口辛了舲（yr d j），I J n。把这些误差进仃加总得到误差函数7工（女D-必i=lj=/2/C d c所要做的就是通过改变w来减小E,以使得所有的输入矢量都尽量与相应的输出矢量相 匹配。因此学习的过程就转化为定义在权值空间上的目标函数

17、E的极小化问题在训练过程中总是以尽可能快的减小E的方式进行。一般它依赖于在权值空间中是否沿梯度方向搜索，所以采用梯度下降法来训练权值。每一个权值 计算：-Z竺dwij其中z为学习率，是控制算法收敛速度的参数在第一阶段得到的总误差平方和又在第二阶段被一层一层地反向传播回去，从输出单元到输入单 元。权值的调整决定于传播过程中的每一步。由于h、和E都是连续可微的，因此，可以应用以下公式计算厶E/A Wj的值：E _dE df. dhAwij dfi dli dwijW的修改可以有两种方式，一是对于每一训练对（Xi, D0都修改一次W另一种方式是输入全部的训练对后再加总 w并进行修改。训练矢量集中训练

18、对的数目称为一个epoch。 当epoch不是非常大的时候，后一种方式能够加快收敛的速度。因为第一种方式只能针对某一特定的训练对减小误差函数，而可能增大其它训练对的误差函数；第二种方式总是以减小总体误差函数 为目标的o所以采用第二种方式【冏。图26 BP神经网络结构图文案大全实用标准文档三、数据记录与分析3.1系统各部分效果验证3.1.1滤波器效果验证为了验证实时滤波的效果以及硬件f生能，在清醒时采集的脑电波中截取了一段约4 信号进行验证秒的求。巾轴恂1 W/i jV S AA!、Z/ilUiliyi1r Vw wi*r I U.1 A-JL- -m其中n为隐含节点数r191。由于本文的样本数

19、为60个，故n取6，即隐含 层有6个节点。隐含层采用对数S形转移函数(Logarithmic sigmoid transfer fun cti on)1+&(3-1)输出层采用线性函数：f(x) = X(3-2)使用Matlab的newff工具，采用梯度下降自适应学习率训练函数创建BP神经网络。学习率定为 0.01，目标误差0.01，最大迭代次数500 o 60个样本中，随机选取50个用于训练，另外未经过训练的 10个用于验证。经过验证，用于验证的样本中有4个判断错误，神经网络的效果并不十分理想。但这四组判断错误的数据中，有一组把W误判为N1,有一组把N1误判为REM误差并未影响对睡 眠深浅度

20、的 判断。此外，对于十组检验样本中的N3(深度睡眠)均为判断错误，可见BP神经网络还是有着一定的准确程度。文案大全实用标准文档四、结论本文主要探究了基于脑电波的便携式睡眠质量监测系统的可行性以及其硬件、软件系统 的技术路线 和实现方法，注重考虑了硬件系统的便携性与成本，在矢注算法的有效性同时，探究其简易和可操作程 度。本文的实验证明，便携式脑电波采集设备具有传统设备无可比拟的移动性，适合个人、家庭用户使用； 利用它可以获得较高质量的脑电信号、用于分析，结合神经网络等算法也可以对于睡眠质量进行可靠的监 测与评估，还能够初步实现对于睡眠的分期。但若要根据医学上严格的睡眠分期标准进行评估，需要采集

21、多导连的脑电以及其他心电、呼吸等信号综合评估，且要求较高的职业技能与素养，只凭便携式睡眠监测系统难以满足要求。本 文的实验结果，将为人们更好地研究便携的睡眠监测系统提供有意义的实验依据与参考。探究单导连脑 电信号与人体睡眠的相矢性，以及脑电波的现代高级分析算法，是对本文实验结果进行探讨的重要理论 基础。此外，探究脑电波的诱发、治疗理论，和便携式脑电波采集设备的准确性、可靠性及其市场化后 的诊断、评估等应用价值，也可以作为本课题 后续的发展方向。参考文献1蔡文英，钟龙云，张作生睡眠脑电波的计算机分析中国科学技术大学学报Vol.20, No.2,Jun.,1990.2 Rechtschaffen,

22、 A. and Kales, A. A Manual of Standardized Terminology, Tech niq ues, and Scori ng System for Sleep Stages of Huma n Subjects. Un iversity of California, Brain In formatio n Service/Brain Research In stitute, Los An geles, CA,1968.3 Iber C, An coli-lsrael S, Chess on A, and Qua n SF for the America

23、n Academy of Sleep Medici ne. The AASM Manual for the Scori ng of Sleep and Associated Eve nts: Rules, Termino logy and Tech ni cal Specificati ons, 1st ed.: Westchester, Illi no is: America n Academy of Sleep Medici ne, 2007.4 Yeo W.C., et al. Naval Res Catholic Univ of Amer 1972; 293297. Larsen R.

24、D., et al. Math Biosci 1976; 31: 237253. Roberts S, et al. Med Biol Eng Comput 1992; 30: 509517.7 Roberts S, et al. IEEE procedi ngs-F 1992; 139(6): 420-425. 张杰，王明时睡眠脑电的研究国外医学生物医学工程分册,1997年第20卷第2期.9 SHAMBROOM, J. R., F a BREGAS, S. E. and JOHNSTONE, J. (2012), Validation of an automated wireless sys

25、tem to monitor sleep in healthy adults. Journal of Sleep Research, 21:221 230. doi: 10.1111/j.1365-2869.2011.00944.x.10志强，叶念渝FIR和IIR数字滤波器的探讨与实现计算机与数字工程，2005, 33(7)11史洁玉.MATLAB信号处理超级学习手册人民邮电出版社，2014. 9.12陈群脑电生物反馈治疗仪与睡眠分期的研究广东工业大学，2013.13李颖洁，邱意弘，朱贻盛脑电信号分析方法及其应用北京：科学出版社，2009.14 Petre Stoica, Ran dolph L. Moses. Spectral An a