CN106388818A - 睡眠状态监测模型的特征信息提取方法和系统 - Google Patents

Abstract

本发明涉及一种睡眠状态监测模型的特征信息提取方法和系统，其中所述方法包括：采集用户在睡眠过程中产生的脑电信号；基于频域信号处理技术从所述脑电信号中提取多个频率段的脑电波；分别计算各个频率段的脑电波的能量在脑电信号总能量中的比例系数，以及计算各个频率段的脑电波在当前帧脑电信号内所占的时间比；根据所述比例系数和时间比确定睡眠状态监测模型的识别任务类型对应的特征信息。基于本发明技术方案提取的特征信息，能够有效提高睡眠状态监测模型识别睡眠状态的准确度，而且也一定程度上提高了识别效率。

Description

睡眠状态监测模型的特征信息提取方法和系统

技术领域

本发明涉及辅助睡眠技术领域，特别是涉及一种睡眠状态监测模型的特征信息提取方法和系统。

背景技术

在睡眠中，人体进行了自我放松及恢复的过程，因此良好的睡眠是保持身体健康的一项基本条件；但是由于工作压力大、生活作息不规律等原因，导致了部分人群的睡眠质量欠佳，表现为失眠、半夜惊醒等。

目前市面上已经有一些设备来帮助人们入睡，提高睡眠质量。例如在某一特定睡眠状态下通过声音、光信号等人工干预，避免在熟睡状态下叫醒用户等。对于辅助睡眠的设备而言，为了真正达到提高用户睡眠质量的目的，正确的识别用户的睡眠状态是非常重要的。

而要识别用户的睡眠状态，目前主要是利用多导睡眠图(Polysomnography,PSG)，又称睡眠脑电图，多导睡眠图利用多种生命体征对睡眠进行分析，在这些体征信号中，脑电处于核心地位；利用脑电波4种节律：δ波(1-3Hz)，θ波(4-7Hz)，α波(8-12Hz)，β波(14-30Hz)的频率特性作为特征信息，采用预先训练的睡眠状态监测模型对睡眠状态进行识别，由于脑电信号容易受到干扰，因此，这种方式的准确度难以得到保证。

发明内容

基于此，有必要针对上述问题，提供一种睡眠状态监测模型的特征信息提取方法和系统，有效地提高睡眠状态监测模型的识别准确度。

一种睡眠状态监测模型的特征信息提取方法，包括：

采集用户在睡眠过程中产生的脑电信号；

基于频域信号处理技术从所述脑电信号中提取多个频率段的脑电波；

分别计算各个频率段的脑电波的能量在脑电信号总能量中的比例系数，以及计算各个频率段的脑电波在当前帧脑电信号内所占的时间比；

根据所述比例系数和时间比确定睡眠状态监测模型的识别任务类型对应的特征信息。

一种睡眠状态监测模型的特征信息提取系统，包括：

采集模块，用于采集用户在睡眠过程中产生的脑电信号；

提取模块，用于基于频域信号处理技术从所述脑电信号中提取多个频率段的脑电波；

计算模块，用于分别计算各个频率段的脑电波的能量在脑电信号总能量中的比例系数，以及计算各个频率段的脑电波在当前帧脑电信号内所占的时间比；

设置模块，用于根据所述比例系数和时间比确定睡眠状态监测模型的识别任务类型对应的特征信息。

上述睡眠状态监测模型的特征信息提取方法和系统，利用采集的脑电信号，提取多个频率段的脑电波，分别计算各个频率段的脑电波的能量在脑电信号总能量中的比例系数及其在当前帧脑电信号内所占的时间比，并将该比例系数和时间比确定睡眠状态监测模型的识别任务类型对应的特征信息。基于该方案提取的特征信息，能够有效提高睡眠状态监测模型识别睡眠状态的准确度，而且也一定程度上提高了识别效率。

附图说明

图1为一个实施例的睡眠状态监测模型的特征信息提取方法的流程图；

图2为滤波处理前后的脑电信号示意图；

图3为一个实施例的睡眠状态监测模型的特征信息提取系统结构示意图。

具体实施方式

下面结合附图阐述本发明的睡眠状态监测模型的特征信息提取方法和系统的实施例。

参考图1所示，图1为本发明的睡眠状态监测模型的特征信息提取方法的流程图，包括：

S101，采集用户在睡眠过程中产生的脑电信号；

在本步骤中，如在对用户进行辅助睡眠时，通过用户佩戴相关传感设备，检测用户的脑电信号，可以以30s为一帧进行采集。

在一个实施例中，考虑到脑电图的信号非常微弱(微伏级)，容易被来自其他部位的生物电信号干扰。例如眼电信号叠加在脑电信号上导致的基线漂移的现象，因此，还可以包括去除每帧(如30s为一帧)脑电信号中的眼电伪迹的步骤：

(1)对采集的脑电信号进行经验模态分解，将其分解成若干个本征模函数，并计算各个本征模函数与同一时刻的眼电信号之间的相关系数；经验模态分解可以包括如下公式：

式中，EEG_original表示原始脑电信号，imf_i表示第i个本征模函数，Re表示残差函数；其中，眼电信号通过用户佩戴相关传感设备，检测用户的眼电信号，可以以30s为一帧进行采集。

(2)将相关系数大于预设阈值的本征模函数和相关系数最大的本征模函数删除，并利用未删除的本征模函数重建脑电信号；重建脑电信号可以包括如下公式：

式中，EEG_pure表示重建的脑电信号，corrcoef表示相关系数，imf表示第i个本征模函数，EOG表示眼电信号，corrcoef_max表示最大的相关系数，thre表示预设的相关系数阈值。

上述实施例的技术方案，只去除高幅度眼电带来的类似于基线漂移的伪迹，保留了原始信号的大部分细节信息。因此有利于后续的基于时域的脑电信号滤波处理。

S102，基于频域信号处理技术从所述脑电信号中提取多个频率段的脑电波；

在本步骤中，在频域对脑电信号进行处理，从脑电信号中提取出多个频率段的脑电波。

在一个实施例中，提取信号特征数据前，还可以对所采集的脑电信号进行滤波处理，滤除高频噪声和工频干扰。例如，脑电信号的有用信息多集中在0-100Hz的范围内，在采集过程中会掺入频率在该范围外的噪声，因此，可以通过滤波手段将其滤除。可以同带通滤波器滤除高频噪声，并设计一个陷波器(50/60Hz)来滤除工频干扰。

参考图2所示，图2为滤波处理前后的脑电信号示意图，上图为原始信号，下图为经过滤波处理之后的信号，可以发现大部分的高频噪声已被滤除。

对于上述脑电信号，由于强度很弱，在信号采集时，极易被外界信号所干扰。

根据频率的不同，多导睡眠图的脑电信号可以分为4种节律脑电波：δ波(1-3Hz)，θ波(4-7Hz)，α波(8-12Hz)，β波(14-30Hz)；在此，可以提取出δ波频段，θ波频段，α波频段，β波频段这四种脑电波。

作为一个实施例，上述提取多个频率段的脑电波的方法，可以包括如下：

对脑电信号进行小波分解，在小波重构中提取所述脑电信号的δ波频段，θ波频段，α波频段，β波频段；

为了更好地分解出所述各种频率波形，小波分解的层数与脑电信号的采样频率满足如下关系：f＝2^N+2，其中，f为脑电信号的采样频率，N为小波分解的层数；例如，当信号的降采样率为128Hz时，可以选择4层分解，当信号的采样率为256Hz时，则可以进行5层分解。

S103，分别计算各个频率段的脑电波的能量在脑电信号总能量中的比例系数，以及计算各个频率段的脑电波在当前帧脑电信号内所占的时间比；

在本步骤中，通过多个频率段的脑电波的能量在脑电信号总能量中的比例系数以及其在当前帧的脑电信号内所占的时间比，作为识别任务的特征信息；由此，以脑电波频域信息基础上，进一步以能量比例系数和脑电波占据脑电信号帧时间比作为特征量进行识别的方式，提高了分类器识别的精确度。

作为一个实施例，上述计算比例系数和时间比的方法，可以包括如下：

(1)分别计算δ波频段，θ波频段，α波频段，β波频段的能量在脑电信号总能量中的比例系数，计算方法可以包括如下公式：

r_δ＝∑(y_δ)²/p_total

r_θ＝∑(y_θ)²/p_total

r_α＝∑(y_α)²/p_total

r_β＝∑(y_β)²/p_total

其中p_total＝∑(y_δ)²+∑(y_θ)²+∑(y_α)²+∑(y_β)²，y_δ，y_θ，y_α和y_β分别表示重构后的δ频段、θ频段、α频段和β频段的信号，r_δ，r_θ，r_α和r_β分别代表δ频段、θ频段、α频段和β频段的信号的能量在总能量的比例。

具体的，通过计算δ波频段，θ波频段，α波频段，β波频段的能量比例系数，可以得到4个特征信息。

(2)计算在一帧脑电信号内，δ波频段，θ波频段，α波频段，β波频段能量最大的时间长度；4计算方法包括如下公式：

式中，c_δ，c_θ，c_α和c_β表示δ频段、θ频段、α频段和β频段的信号在当前帧内所占能量比例最大的时间长度，分别表示第i秒内δ频段、θ频段、α频段和β频段的信号的能量在总能量的比例。

例如，在以30s为一帧进行对脑电信号进行采集时，可以计算该帧30s内，各个脑电波所占的时间；通过计δ波频段，θ波频段，α波频段，β波频段的时间比，可以得到4个特征信息。

由此，可以得到8个特征信息，应用于睡眠状态监测模型的识别任务。

S104，根据所述比例系数和时间比确定睡眠状态监测模型的识别任务类型对应的特征信息；

在本步骤中，是将前述识别的比例系数和时间比作为特征信息，应用于睡眠状态监测模型的识别任务中；如用于训练新的睡眠状态监测模型(分类器)或者作为样本数据输入已有的睡眠状态监测模型(分类器)中进行睡眠状态识别。

上述睡眠状态监测模型，可以采用RBF核的SVM(Support Vector Machin，支持向量机)分类器，也可以采用神经网络、决策树的分类器等。

作为一个实施例，基于前述计算的δ波频段，θ波频段，α波频段，β波频段对应的比例系数和时间比，分别其设为睡眠状态监测模型的清醒状态和睡眠状态识别任务的特征信息，即应用于识别用户处于清醒还是睡眠状态。

作为一个实施例，本发明实施例提供的睡眠状态监测模型的特征信息提取方法，应用于睡眠状态监测模型进行睡眠状态识别时，可以用来进行包括清醒，非眼快动睡眠和眼快动睡眠等识别任务。

在一个实施例中，步骤S102中，在对所述脑电信号进行小波分解的步骤前，还可以提取脑电信号的基线，计算所述基线的变化幅度；其中，所述变化幅度为基线最大值减去最小值；

以及在去掉脑电信号的基线后，执行所述对脑电信号进行小波分解的步骤，获得小波系数。

在步骤S103中，根据小波系数计算小波系数的特征参数；其中，所述特征参数包括小波系数的均值、方差、峭度系数和/或斜度系数。

在步骤S104中，则将所述比例系数、时间参数、变化幅度和特征参数设为睡眠状态监测模型的清醒状态、非眼快动睡眠状态和眼快动睡眠状态识别任务的特征信息。

进一步地，还可以用来进行包括非眼快动睡眠的入睡期(S1)，非眼快动睡眠期(S2)，中等睡眠期(S3)和眼快动睡眠期(S4)4个状态的等识别任务。

对于非眼快动睡眠可以分为4个时期：S1期(完全清醒至睡眠之间的过渡阶段，脑电波以θ波为主，不出现纺锤波或K综合波)；S2期(浅睡阶段，脑电波为纺锤波与K综合波，δ波少于20％)；S3期(中等深度睡眠，脑电波δ波占20％～50％)；S4期(深睡，脑电波δ波占50％以上)。

在一个实施例中，在上述实施例的步骤S102基础上，在去掉脑电信号的基线后，进一步对脑电信号进行降采样处理(如原信号采样率为512Hz，降采样后为128Hz)，得到降采样信号。

在步骤S103基础上，计算所述降采样信号的基于非线性动力学的特征量；其中，所述基于非线性动力学的特征量可以包括LZ复杂度、样本熵和/或近似熵。

在步骤S104中，则将所述比例系数、时间参数、变化幅度、特征参数和特征量设为睡眠状态监测模型的清醒状态、S1-S4期非眼快动睡眠状态和眼快动睡眠状态识别任务的特征信息。

作为一个实施例中，在利用上述特征信息基础上，还可以采集用户的血氧浓度参数，并根据所述浓度参数计算血氧饱和度参数，将所述血氧饱和度参数设为睡眠状态监测模型的识别任务类型对应的特征信息。

参考图3所示，图3为一个实施例的睡眠状态监测模型的特征信息提取系统结构示意图，包括：

采集模块101，用于采集用户在睡眠过程中产生的脑电信号；

提取模块102，用于基于频域信号处理技术从所述脑电信号中提取多个频率段的脑电波；

计算模块103，用于分别计算各个频率段的脑电波的能量在脑电信号总能量中的比例系数，以及计算各个频率段的脑电波在当前帧脑电信号内所占的时间比；

设置模块104，用于根据所述比例系数和时间比确定睡眠状态监测模型的识别任务类型对应的特征信息。

本发明的睡眠状态监测模型的特征信息提取系统与本发明的睡眠状态监测模型的特征信息提取方法一一对应，在上述睡眠状态监测模型的特征信息提取方法的实施例阐述的技术特征及其有益效果均适用于睡眠状态监测模型的特征信息提取系统的实施例中，特此声明。

以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种睡眠状态监测模型的特征信息提取方法，其特征在于，包括：

采集用户在睡眠过程中产生的脑电信号；

基于频域信号处理技术从所述脑电信号中提取多个频率段的脑电波；

分别计算各个频率段的脑电波的能量在脑电信号总能量中的比例系数，以及计算各个频率段的脑电波在当前帧脑电信号内所占的时间比；

根据所述比例系数和时间比确定睡眠状态监测模型的识别任务类型对应的特征信息。

2.根据权利要求1所述的睡眠状态监测模型的特征信息提取方法，其特征在于，还包括：

对采集的脑电信号进行经验模态分解，将其分解成若干个本征模函数，并计算各个本征模函数与同一时刻的眼电信号之间的相关系数；

将相关系数大于预设阈值的本征模函数和相关系数最大的本征模函数删除，并利用未删除的本征模函数重建脑电信号。

3.根据权利要求1所述的睡眠状态监测模型的特征信息提取方法，其特征在于，所述经验模态分解包括如下公式：

${\mathrm{EEG}}_{original}=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{imf}}_i+\mathrm{Re}$

式中，EEG_original表示原始脑电信号，imf_i表示第i个本征模函数，Re表示残差函数；

所述重建脑电信号包括如下公式：

${\mathrm{EEG}}_{pure}=\underset{i=1}{\overset{m}{\mathrm{\&Sigma;}}}{\mathrm{imf}}_i,corrcoef({\mathrm{imf}}_i,EOG)<\min (thre,{\mathrm{corrcoef}}_{\max })$

式中，EEG_pure表示重建的脑电信号，corrcoef表示相关系数，imf表示第i个本征模函数，EOG表示眼电信号，corrcoef_max表示最大的相关系数，thre表示预设的相关系数阈值。

4.根据权利要求1所述的睡眠状态监测模型的特征信息提取方法，其特征在于，所述基于频域信号处理技术从所述脑电信号中提取多个频率段的脑电波的步骤包括：

对脑电信号进行小波分解，在小波重构中提取所述脑电信号的δ波频段，θ波频段，α波频段，β波频段；

所述分别计算各个频率段的脑电波的能量在脑电信号总能量中的比例系数，以及计算各个频率段的脑电波在当前帧脑电信号内所占的时间比的步骤包括：

分别计算δ波频段，θ波频段，α波频段，β波频段的能量在脑电信号总能量中的比例系数；

分别计算在一帧脑电信号内，δ波频段，θ波频段，α波频段，β波频段能量最大的时间长度；

所述根据所述比例系数和时间比确定睡眠状态监测模型的识别任务类型对应的特征信息的步骤包括：

分别将所述δ波频段，θ波频段，α波频段，β波频段对应的比例系数和时间比设为睡眠状态监测模型的清醒状态和睡眠状态识别任务的特征信息。

5.根据权利要求4所述的睡眠状态监测模型的特征信息提取方法，其特征在于，所述分别计算δ波频段，θ波频段，α波频段，β波频段的能量在脑电信号总能量中的比例系数的步骤包括如下公式：

r_δ＝Σ(y_δ)²/p_total

r_θ＝Σ(y_θ)²/p_total

r_α＝Σ(y_α)²/p_total

r_β＝∑(y_β)²/p_total

其中p_total＝∑(y_δ)²+∑(y_θ)²+∑(y_α)²+∑(y_β)²，y_δ，y_θ，y_α和y_β分别表示重构后的δ频段、θ频段、α频段和β频段的信号，r_δ，r_θ，r_α和r_β分别代表δ频段、θ频段、α频段和β频段的信号的能量在总能量的比例。

6.根据权利要求2所述的睡眠状态监测模型的特征信息提取方法，其特征在于，所述计算在一帧脑电信号内，δ波频段，θ波频段，α波频段，β波频段能量最大的时间长度在当前帧脑电信号内所占的时间的方法包括如下公式：

$c_{\mathrm{\&delta;}}=\underset{i=1}{\overset{30}{\mathrm{\&Sigma;}}}{f}_{\mathrm{\&delta;}}^i,f_{\mathrm{\&delta;}}^i=\left\{\begin{array}{c}1,if{r}_{\mathrm{\&delta;}}^i=\max (r_{\mathrm{\&delta;}}^i,r_{\mathrm{\&theta;}}^i,r_{\mathrm{\&alpha;}}^i,r_{\mathrm{\&beta;}}^i)\\ 0,if{r}_{\mathrm{\&delta;}}^i\&NotEqual;\max (r_{\mathrm{\&delta;}}^i,r_{\mathrm{\&theta;}}^i,r_{\mathrm{\&alpha;}}^i,r_{\mathrm{\&beta;}}^i)\end{array}\right.$

$c_{\mathrm{\&theta;}}=\underset{i=1}{\overset{30}{\mathrm{\&Sigma;}}}{f}_{\mathrm{\&theta;}}^i,f_{\mathrm{\&theta;}}^i=\left\{\begin{array}{c}1,if{r}_{\mathrm{\&theta;}}^i=\max (r_{\mathrm{\&delta;}}^i,r_{\mathrm{\&theta;}}^i,r_{\mathrm{\&alpha;}}^i,r_{\mathrm{\&beta;}}^i)\\ 0,if{r}_{\mathrm{\&theta;}}^i\&NotEqual;\max (r_{\mathrm{\&delta;}}^i,r_{\mathrm{\&theta;}}^i,r_{\mathrm{\&alpha;}}^i,r_{\mathrm{\&beta;}}^i)\end{array}\right.$

$c_{\mathrm{\&alpha;}}=\underset{i=1}{\overset{30}{\mathrm{\&Sigma;}}}{f}_{\mathrm{\&alpha;}}^i,f_{\mathrm{\&alpha;}}^i=\left\{\begin{array}{c}1,if{r}_{\mathrm{\&alpha;}}^i=\max (r_{\mathrm{\&delta;}}^i,r_{\mathrm{\&theta;}}^i,r_{\mathrm{\&alpha;}}^i,r_{\mathrm{\&beta;}}^i)\\ 0,if{r}_{\mathrm{\&alpha;}}^i\&NotEqual;\max (r_{\mathrm{\&delta;}}^i,r_{\mathrm{\&theta;}}^i,r_{\mathrm{\&alpha;}}^i,r_{\mathrm{\&beta;}}^i)\end{array}\right.$

$c_{\mathrm{\&beta;}}=\underset{i=1}{\overset{30}{\mathrm{\&Sigma;}}}{f}_{\mathrm{\&beta;}}^i,f_{\mathrm{\&beta;}}^i=\left\{\begin{array}{c}1,if{r}_{\mathrm{\&beta;}}^i=\max (r_{\mathrm{\&delta;}}^i,r_{\mathrm{\&theta;}}^i,r_{\mathrm{\&alpha;}}^i,r_{\mathrm{\&beta;}}^i)\\ 0,if{r}_{\mathrm{\&beta;}}^i\&NotEqual;\max (r_{\mathrm{\&delta;}}^i,r_{\mathrm{\&theta;}}^i,r_{\mathrm{\&alpha;}}^i,r_{\mathrm{\&beta;}}^i)\end{array}\right.$

式中，c_δ，c_θ，c_α和c_β表示δ频段、θ频段、α频段和β频段的信号在当前帧内所占能量比例最大的时间长度，分别表示第i秒内δ频段、θ频段、α频段和β频段的信号的能量在总能量的比例。

7.根据权利要求1所述的睡眠状态监测模型的特征信息提取方法，其特征在于，在对所述脑电信号进行小波分解的步骤前，还包括：

提取脑电信号的基线，计算所述基线的变化幅度；

以及在去掉脑电信号的基线后，执行所述对脑电信号进行小波分解的步骤，获得小波系数，并根据小波系数计算小波系数的特征参数；

所述根据所述比例系数和时间比确定睡眠状态监测模型的识别任务类型对应的特征信息的步骤包括：

将所述比例系数、时间参数、变化幅度和特征参数设为睡眠状态监测模型的清醒状态、非眼快动睡眠状态和眼快动睡眠状态识别任务的特征信息。

8.根据权利要求7所述的睡眠状态监测模型的特征信息提取方法，其特征在于，所述变化幅度为基线最大值减去最小值；所述特征参数包括小波系数的均值、方差、峭度系数和/或斜度系数。

9.根据权利要求7所述的睡眠状态监测模型的特征信息提取方法，其特征在于，在去掉脑电信号的基线后，还包括：

对脑电信号进行降采样处理，得到降采样信号；计算所述降采样信号的基于非线性动力学的特征量；

所述根据所述比例系数和时间比确定睡眠状态监测模型的识别任务类型对应的特征信息的步骤包括：

将所述比例系数、时间参数、变化幅度、特征参数和特征量设为睡眠状态监测模型的清醒状态、S1-S4期非眼快动睡眠状态和眼快动睡眠状态识别任务的特征信息。

10.一种睡眠状态监测模型的特征信息提取系统，其特征在于，包括：

采集模块，用于采集用户在睡眠过程中产生的脑电信号；

提取模块，用于基于频域信号处理技术从所述脑电信号中提取多个频率段的脑电波；

计算模块，用于分别计算各个频率段的脑电波的能量在脑电信号总能量中的比例系数，以及计算各个频率段的脑电波在当前帧脑电信号内所占的时间比；

设置模块，用于根据所述比例系数和时间比确定睡眠状态监测模型的识别任务类型对应的特征信息。