基于脑电、心率变异及两者相干性的自动睡眠分期方法
【专利摘要】本发明公开了一种基于脑电、心率变异及两者相干性的自动睡眠分期方法，所述方法包括以下步骤：提取被测试者的脑电信号和心率变异信号；从脑电信号、心率变异信号中共提取9个特征向量；对9个特征向量进行主成分分析；特征提取后，使用支持向量机分类器对特征进行识别，进行自动睡眠分期。该方法可有效地提高自动睡眠分期系统准确性和简便性，并获得可观的社会效益和经济效益。最佳实施方案拟采用专利转让、技术合作或产品开发。由于该技术操作简单，敏感性强，基于该技术开发的产品可应用于睡眠监测、睡眠科研等多种场景等。
【专利说明】基于脑电、心率变异及两者相干性的自动睡眠分期方法
【技术领域】
[0001]本发明涉及自动睡眠分期领域，特别涉及一种基于脑电、心率变异及两者相干性的自动睡眠分期方法。
【背景技术】
[0002]睡眠是人体重要的生理活动，人类生命约有1/3的时间在睡眠中度过。睡眠是机体自我修复和自我恢复的过程，而良好的睡眠是维持身心健康的基本需要。进行睡眠分期是探讨睡眠期间生理调节变化规律、进行睡眠质量评价的基础和前提，更重要的是，作为睡眠监护仪最核心的技术模块，睡眠分期为针对睡眠的多模式生理学研究提供了一种可靠的量化指标。
[0003]睡眠分期就是根据睡眠期间人体生理信号的不同将睡眠分为不同的阶段。主要分为非快速眼动期(NREM)和快速眼动期(REM)。非快速眼动期从夜间入睡开始，又细分为4期，NREM-1为入睡期，NREM-1I为浅睡期，NREM-1II为中度睡眠期，NREM-1V为深度睡眠期。在这个阶段人的呼吸变浅、变慢，心率变慢、血压下降，全身肌肉松弛，无明显的眼球运动。快速眼动期的特征是眼球快速转动。人体的感觉功能进一步减退，肌肉也更加松弛，肌腱反射消失。这个阶段，体内各种代谢功能都显著增加，以保证大脑组织蛋白的合成和消耗物质的补充，使神经系统正常发育，并为第二天的活动积蓄能量。图1是正常人的睡眠周期图。
[0004]目前,基于各种生理信号的睡眠监测方法主要有多导睡眠图(PSG)、EEG (脑电信号)监测法、ECG (心电信号)监测法、血氧饱和度监测法和身体活动变化记录法等。(I)PSG是进行睡眠医学研究的经典方法，也称金方法。它同时记录EEG、E0G (眼电信号)和EMG (肌电信号)，在此基础上根据R&K规则[1]，进行睡眠状态和睡眠阶段的判定，但此方法需要佩戴至少十枚电极，被试者经常感觉不适，对正常睡眠干扰较大。而且一般采用对脑电波形进行人工判读的方法，但这需要很大的工作量，且有赖于判读者的个人经验，往往夹带有主观因素，缺乏一致性客观标准。(2) EEG监测法、ECG监测法、血氧饱和度监测法和身体活动变化记录法等多采用单一参数或将各生理参数简单融合作为自动睡眠分期的指标，正确率较低。
【发明内容】
[0005]本发明提供了一种基于脑电、心率变异及两者相干性的自动睡眠分期方法，本发明提高了睡眠的质量，且提高了自动睡眠分期的质量，详见下文描述:
[0006]一种基于脑电、心率变异及两者相干性的自动睡眠分期方法，所述方法包括以下步骤:
[0007](I)提取被测试者的脑电信号和心率变异信号；
[0008](2)从脑电信号、心率变异信号中共提取9个特征向量；
[0009](3)对9个特征向量进行主成分分析；
[0010](4)特征降维后，使用支持向量机分类器对特征进行识别，进行自动睡眠分期。[0011]所述从脑电信号、心率变异信号中共提取9个特征向量的操作具体为:
[0012]1)通过对脑电信号进行小波变换、希尔伯特黄变换和奇异值分解获取脑电信号的小波熵、Hilbert-Huang熵，以及奇异值第一主成分，分别记为特征向量P1, P2，P3 ；
[0013]2)计算极低频VLF、低频LF和高频HF的频谱熵，并通过小波变换计算心率变异信号的分形维数，分别记为特征向量p4，P5, P6, P7 ；
[0014]3)对脑电信号以及心率变异信号进行功率谱密度分析，并进行归一化处理，然后分别计算脑电信号delta频段与心率变异信号LF、HF的相干系数，分别记为特征向量P8，P90
[0015]所述对脑电信号以及心率变异信号进行功率谱密度分析，并进行归一化处理，然后分别计算脑电信号delta频段与心率变异信号LF、HF的相干系数，分别记为特征向量P8, P9的操作具体为:
[0016]I)截取同一时间段的EEG、HRV信号，进行预处理，包括变参考、降采样、去噪声与干扰，预处理后的EEG信号记为X，HRV信号记为y ；
[0017]2)采用Welch算法计算EEG、HRV的功率谱密度以及两者的互功率谱密度，分别记为 Pxx (fi), Pyy (f2), Pxy (^, f2)；
[0018]3)用改进的相干函数计算EEG信号delta频段与HRV信号LF、HF的相干系数。
[0019]所述对9个特征向量进行主成分分析的操作具体为:
[0020](I)对原始数据P_进行标准化处理，矩阵中的元素减去所在列的均值，然后除以所在列的标准差，使得每个变量的均值为0，方差为1，得到矩阵PnJ ；
[0021](2)然后求其协方差矩阵C_，Pn*9*中任两列之间可以计算两变量间的协方差；
[0022](3)对协方差矩阵C_进行特征根分解，得到特征根矩阵及特征向量；
[0023](4)求原始数据P_在新的矢量空间中的投影，即主成分向量组F_:
[0024](5)求累积贡献率；
[0025](6)选定预设的累积贡献率，使前d个主成分Fn?作为新的数据进行模式识别。
[0026]本发明提供的技术方案的有益效果是:本发明通过计算EEG信号delta频段与HRV信号LF、HF的相干系数，与提取的EEG、HRV多种频域特征和非线性特征融合作为睡眠分期的特征矩阵，PCA去冗余信息后作为后续模式识别的输入，从而准确、客观的进行自动睡眠分期。该项发明可有效地提高睡眠分期的准确性和简便性，并获得可观的社会效益和经济效益；为睡眠监护仪等提供了最核心的技术支持，给实际应用带来便捷，并且可应用于多种工作场景。
【专利附图】
【附图说明】
[0027]图1为睡眠周期图；
[0028]图2为自动睡眠分期方法流程图。
【具体实施方式】
[0029]为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
[0030]睡眠作为一种生理状态，必然要在身体的多个部位及多个生理过程中反映出来。EEG是中枢神经信号，是睡眠分期最敏感最准确的生理信号，心电信号是自主神经信号，提取的心率变异信号(HRV)反映交感神经系统和副交感神经系统的活动水平，而自主神经系统的调控能力随着睡眠程度的不同而改变。有研究表明，自主神经活动与睡眠在皮层下有共同的调节中枢，心电信号作为主要的自主神经信号，与脑电信号的有效融合，势必能够提高自动睡眠分期的正确率。
[0031]101:提取被测试者的脑电信号和心率变异信号；
[0032]EEG包含大量的生理信息，经过分析处理后，可对睡眠进行分期。相对于其他生理参数，EEG反映睡眠各期特征最为明显，是目前睡眠分期的“金标准”，所以近年来通过脑电信号自动睡眠分期的准确性和客观性研究一直受到人们的关注，研究也最为成熟。EEG通常包含 5 种基本节律:delta 波(0.5_4Hz), theta 波(4_8Hz), aleph 波(8_13Hz), beta 波(13_30Hz ), gamma 波(30_50Hz )。
[0033]心率变异即逐次窦性心跳间期之间的微小涨落，直接受自主神经系统中的交感神经和副交感神经的双重支配。心率变异性是自主神经系统功能最直接、有效的反映，而自主神经系统的调控能力随着睡眠程度的不同而改变，研究表明HRV是进行自动睡眠分期的很好指标。HRV典型的频谱可有三个峰，大致位于0.04Hz以下，0.05-0.15Hz和>0.15Hz，分别称为极低频(VLF)、低频(LF)和高频(HF)峰。
[0034]在本方法中，由于该两个参数测量方法成熟，与其他生理参数相比能更准确、敏感的反应各睡眠阶段，并且心脑电信号在不同的睡眠阶段相关程度不同，所以采用这两种生理参数有效融合来进行自动睡眠分期在实际应用中具有优越性。
[0035]102:从脑电信号、心率变异信号中共提取9个特征向量；
[0036]1)EEG特征提取:通过对脑电信号进行小波变换、希尔伯特黄(Hilbert-Huang)变换和奇异值分解获取脑电信号的小波熵[2]、Hilbert-Huang熵[3]，以及奇异值第一主成分(奇异谱中的最大值)M，分别记为特征向量P1, P2, P3。
[0037]2) HRV特征提取:计算极低频(VLF)、低频(LF)和高频(HF)的频谱熵，并通过小波变换计算HRV的分形维数M Μ，分别记为特征向量P4，P5, P6, P7。
[0038]3) EEG信号delta频段与HRV参数的相干系数:对EEG信号以及HRV信号进行功率谱密度分析，并进行归一化处理，然后利用改进的相干函数分别计算EEG信号delta频段与HRV信号LF、HF的相干系数，分别记为特征向量P8，P9。
[0039]相干分析是在频域上描述两个信号相关程度的实值函数。频域相干性是对两个信号在不同频率下的相位一致性的量度，当两者的某一频率成分fi相位一致时，
[0040]
【权利要求】
1.一种基于脑电、心率变异及两者相干性的自动睡眠分期方法，其特征在于，所述方法包括以下步骤:(1)提取被测试者的脑电信号和心率变异信号；(2)从脑电信号、心率变异信号中共提取9个特征向量；(3 )对9个特征向量进行主成分分析；(4)特征降维后，使用支持向量机分类器对特征进行识别，进行自动睡眠分期。
2.根据权利要求1所述的一种基于脑电、心率变异及两者相干性的自动睡眠分期方法，其特征在于，所述从脑电信号、心率变异信号中共提取9个特征向量的操作具体为:1)通过对脑电信号进行小波变换、希尔伯特黄变换和奇异值分解获取脑电信号的小波熵、Hilbert-Huang熵，以及奇异值第一主成分，分别记为特征向量P1, P2, P3 ；2)计算极低频VLF、低频LF和高频HF的频谱熵，并通过小波变换计算心率变异信号的分形维数，分别记为特征向量P4，P5, P6, P7 ；3)对脑电信号以及心率变异信号进行功率谱密度分析，并进行归一化处理，然后分别计算心率变异信号delta频段与心率变异信号LF、HF的相干系数，分别记为特征向量P8，P90
3.根据权利要求2所述的一种基于脑电、心率变异及两者相干性的自动睡眠分期方法，其特征在于，所述对脑电信号以及心率变异信号进行功率谱密度分析，并进行归一化处理，然后分别计算脑电信号delta频段与心率变异信号LF、HF的相干系数，分别记为特征向量P8，P9的操作具体为:1)截取同一时间段 的EEG、HRV信号，进行预处理，包括变参考、降采样、去噪声与干扰，预处理后的EEG信号记为X，HRV信号记为y ；2)采用Welch算法计算EEG、HRV的功率谱密度以及两者的互功率谱密度，分别记为Pxx (f) , Pyy (f2), Pxy (f1； f2)；3)用改进的相干函数计算EEG信号delta频段与HRV信号LF、HF的相干系数。
4.根据权利要求1所述的一种基于脑电、心率变异及两者相干性的自动睡眠分期方法，其特征在于，所述对9个特征向量进行主成分分析的操作具体为:(1)对原始数据p_进行标准化处理，矩阵中的元素减去所在列的均值，然后除以所在列的标准差，使得每个变量的均值为0，方差为1，得到矩阵PnJ ；(2)然后求其协方差矩阵C_，Pn*9*中任两列之间可以计算两变量间的协方差；(3)对协方差矩阵Cw9进行特征根分解，得到特征根矩阵及特征向量；(4)求原始数据P_在新的矢量空间中的投影，即主成分向量组F_:(5)求累积贡献率；(6)选定预设的累积贡献率，使前d个主成分Fn?作为新的数据进行模式识别。
【文档编号】A61B5/00GK103584840SQ201310608769
【公开日】2014年2月19日 申请日期:2013年11月25日 优先权日:2013年11月25日
【发明者】刘爽, 明东, 綦宏志, 赵欣, 周鹏, 张力新, 万柏坤 申请人:天津大学