基于离散小波变换的脑电信号睡眠分期研究基于离散小波变换的脑电信号睡眠分期研究
提出基于脑电信号(EEG)的睡眠分期研究。利用离散小波变换(DWT)的db8小波分解得到的细节分量作为
信号新的表达,把各个细节分量能量作为特征,建立带高斯径向基核函数(RBF)的非线性支持向量机
(SVM)模型。研究发现,其对睡眠分期研究的方案是可行的,满足模型对泛化能力的要求。
摘摘 要要: 提出基于脑电信号(EEG)的
关键词关键词: 睡眠分期;离散小波变换;支持向量机
0 引言引言
睡眠是人体一种非常重要且不可或缺的生理活动,但工作压力、疾病、不良的生活习惯等因素会影响人的睡眠质量,长期
睡眠质量低下会引发各种心理、生理疾病。依据生理信号数据自动识别睡眠,分期对睡眠状态与睡眠质量的科学评估以及相关
疾病的诊断及治疗有很重要的意义和参考价值。
对睡眠的研究往往是基于相关生理数据。在清醒和不同睡眠阶段,人的脑电信号(Electroencephalogram,EEG)在不
同的脑区会出现不同的节律变化分布,可以表征人体的不同睡眠阶段。参考文献[1-3]提取最小值、最大值、均值、标准差等频
域特征和排列组合熵、样本熵等非线性动力学参数,但其样本较少,模型推广性能不足。参考文献[4]分别使用Choi–Williams
分布、连续小波变换和Hilbert-Huang变换(Hilbert-Huang Transform,HHT)提取特征,利用随机森林分类器说明不同方法
提取特征对分类器的准确率有一定影响,其中使用CWT提取特征的分类器准确率最高。参考文献[5-6]分别建立复发性神经分
类器和支持向量机(Support Vector Machine,SVM)模型,对非快速眼动深睡期没有细分。而参考文献[7]使用EEG各个节
律波频带的相对功率作为特征,所建立的SVM模型准确率高于神经网络(Neural Network,NN)模型,该研究只是将睡眠分
为4个阶段,未具体区分非快速眼动睡眠各个阶段,且清醒期的识别率偏低,并没有突出提高清醒期识别率的基础作用。参考
文献[8]研究了25位受试者睡眠EEG、眼电(Electrooculogram,EOG)和肌电(Electromyogram,EMG)信号,提取特征后
在深度信念网络(Deep Belief Networks,DBN)基础上建立的马尔科夫模型准确率高于未提取特征的模型准确率,说明提取
特征是建立模型的重要步骤。由于样本来自疑似睡眠呼吸障碍者,该模型具有较强的针对性和局限性。
由于人的生理信号存在个体差异,当被试者人数增加时与信号特征离散性增大,而导致模型的推广能力下降。本文使用
Sleep-EDF数据库EEG Fpz-Cz、EEG Pz-Oz这2个通道的数据。采样频率为100 Hz,由专家根据R&K标准[9],基于经验标定
受试者每30 s区间的状态。R&K标准将睡眠分为6个阶段,分别是清醒期(Wake)、非快速眼动期(Non-Rapid Eye
Movement,NREM)和快速眼动期(Rapid Eye Movement,REM),其中NREM又可分为浅睡期(包括NREM1、
NREM2)和深睡期(包括NREM3、NREM4)。
1 离散小波变换离散小波变换
利用EEG将睡眠分期的关键是对能表征各个睡眠阶段的特征进行提取。受试者在长时间的测试过程中,电极松动或移动
会带来基线漂移,而对频谱图的研究发现,低频信号的幅值相对较大,尤其以1 Hz以内的幅值最为显著,据此可以判断,信
号主要存在基线漂移现象。通过离散小波变换处理数据,在消除基线漂移之后提取各个分量的能量作为模型特征。
小波变换是在傅里叶变换的基础上发展起来的一种时频分析方法,在处理非平稳信号方面,具有良好的时频特性。由于实
验采集到的信号大多是经过采集系统采集的离散数据,因此离散小波变换的应用非常广泛[10]。将信号f(t)进行k尺度离散小
波变换,计算公式如下:
原始时域信号经由离散小波变换得到逼近分量 Ai(t)和细节分量Di(t)(i=1,2,...,k)。本文使用db8小波分解,
分别得到细节分量Di(t)(i=1,2,...,8)和逼近分量A8(t)。细节分量Di(t)(i=1,2,...,8)从D1(t)到 D8(t)的
频率逐次降低。其中,D7(t)和D8(t)的频谱主要集中在1 Hz以内,故可以将其作为基线漂移分量,置零即可。图1显示了
一名受试者在REM期的EEG Fpz-Cz通道的30 s信号经离散小波变换得到细节分量Di(t)(i=1,2,…,6)的频域信号。