睡眠状态分析中清醒状态检测方法和系统与流程

技术特征：

1.一种睡眠状态分析中清醒状态检测方法，其特征在于，包括：

在用户开始睡眠过程后，采集用户的实时脑电信号；

对所述实时脑电信号进行小波分解，并根据设定低频段的小波系数进行信号重建得到脑电信号；

计算所述脑电信号的样本熵，将该样本熵与预先计算的样本熵阈值进行比较；

若所述样本熵大于所述样本熵阈值，则判定用户当前处于清醒状态；

所述计算所述脑电信号的样本熵包括，从开始采集脑电信号后，获取设定时间内的脑电信号；将设定时间内的脑电信号分成多个样本，并分别计算各个脑电信号样本的样本熵，得到样本熵集合；根据所述样本熵集合计算样本熵阈值；

所述根据所述样本熵集合计算样本熵阈值的计算公式如下：

sampen_vali＝sampen(y[p_start:p_end])

p_start＝(i-1)*time_length*fs+1

p_end＝p_start+time_length*fs-1

p_end＜T·fs

式中，其中sampen_thre为样本熵阈值，sampen_vali为样本熵集合中第i个样本的样本熵，n为样本熵集合中脑电信号样本的数量，sampen为求样本熵的运算，其输入y[p_start:p_end]为脑电信号y在第p_start点开始到第p_end点为止的部分，time_length为计算样本熵的每个样本的时间长度，fs为脑电信号的采样率，T为开始采集脑电信号后的设定时间，v为设定参数。

2.根据权利要求1所述的睡眠状态分析中清醒状态检测方法，其特征在于，所述参数v的取值如下：

其中，

x＝v

式中，X表示样本熵集合，Xi表示样本熵集合X中第i个元素，Φ(x)表示标准正态分布函数的积分，P(X≤x)表示样本熵的集合X中的取值小于x的概率。

3.根据权利要求2所述的睡眠状态分析中清醒状态检测方法，其特征在于，所述T＝300s，time_length＝30s，v＝2.58。

4.根据权利要求1至3任一项所述的睡眠状态分析中清醒状态检测方法，其特征在于，还包括：

在用户开始睡眠过程后，采集用户的实时眼电信号；

对所述实时眼电信号进行小波分解，并根据设定低频段的小波系数进行信号重建得到眼电信号；

根据同一时刻脑电信号和眼电信号的相关性以及眨眼眼电波形的特征，在眼电信号上检测眨眼活动；

当检测到眨眼活动时，判定用户当前处于清醒状态。

5.根据权利要求4所述的睡眠状态分析中清醒状态检测方法，其特征在于，所述根据同一时刻脑电信号和眼电信号的相关性以及眨眼眼电波形的特征，在眼电信号上检测眨眼活动的步骤包括：

利用具有设定信号幅度范围和时间长度的滑动窗口截取低频眼电信号；

分别计算滑动窗口内眼电信号与同一时刻脑电信号的波形的相关系数，滑动窗口内眼电信号波形尖峰的尖锐程度参数以及尖峰的持续时间；

若所述相关系数、尖锐程度参数和持续时间均分别满足预设的相关系数阈值、尖锐程度参数阈值和持续时间阈值，判断该滑动窗口内眼电信号存在眨眼活动。

6.根据权利要求5所述的睡眠状态分析中清醒状态检测方法，其特征在于，所述计算滑动窗口内眼电信号波形尖峰的尖锐程度参数的步骤包括：

分别计算眼电信号波形在滑动窗口内的上部区域面积和下部区域面积，计算公式如下：

式中，pi为滑动窗口内的眼电信号，pmax为滑动窗口内眼电信号的最大值，pmin为滑动窗口内眼电信号的最小值，areaup表示上部区域面积，areadown表示下部区域面积；

根据所述上部区域面积和下部区域面积计算所述眼电信号波形尖峰的面积，计算公式如下：

式中，blinkarea表示尖峰的面积，if表示满足条件；

根据尖峰面积计算尖锐程度参数，计算公式如下：

blinkratio＝blinkarea/in-blinkarea

式中，blinkratio表示尖锐程度参数，in-blinkarea表示非尖峰部分的面积。

7.根据权利要求6所述的睡眠状态分析中清醒状态检测方法，其特征在于，所述计算滑动窗口内眼电信号波形尖峰的持续时间的步骤包括：

根据所述上部区域面积和下部区域面积计算所述眼电信号波形尖峰的方向，计算公式如下：

式中，blinkdirection表示尖峰方向，标记为1表示尖峰向下，标记为-1表示尖峰向上；

在尖峰向下时，根据尖峰两侧的局部极大值点计算尖峰的持续时间；或者在尖峰向上时，根据尖峰两侧的局部极小值点计算尖峰的持续时间；计算公式如下：

blinkts＝vertexright-vertexleft,

式中，blinkts表示尖峰的持续时间，vertexright表示尖峰的右侧顶点时刻，表示尖峰的左侧顶点时刻，right_min_loca表示尖峰右侧的局部极小值点，right_max_loca表示尖峰右侧的局部极大值点，left_min_loca表示尖峰左侧的局部极小值点，left_max_loca表示尖峰左侧的局部极大值点，s.t.表示约束条件。

8.一种睡眠状态分析中清醒状态检测系统，其特征在于，包括：

脑电信号采集模块，用于在用户开始睡眠过程后，采集用户的实时脑电信号；

低频信号重建模块，用于对所述实时脑电信号进行小波分解，并根据设定低频段的小波系数进行信号重建得到脑电信号；

样本熵比较模块，用于计算所述脑电信号的样本熵，将该样本熵与预先计算的样本熵阈值进行比较；

清醒状态判定模块，用于若所述样本熵大于所述样本熵阈值，则判定用户当前处于清醒状态；

所述计算所述脑电信号的样本熵包括，从开始采集脑电信号后，获取设定时间内的脑电信号；将设定时间内的脑电信号分成多个样本，并分别计算各个脑电信号样本的样本熵，得到样本熵集合；根据所述样本熵集合计算样本熵阈值；

所述根据所述样本熵集合计算样本熵阈值的计算公式如下：

sampen_vali＝sampen(y[p_start:p_end])

p_start＝(i-1)*time_length*fs+1

p_end＝p_start+time_length*fs-1

p_end＜T·fs

式中，其中sampen_thre为样本熵阈值，sampen_vali为样本熵集合中第i个样本的样本熵，n为样本熵集合中脑电信号样本的数量，sampen为求样本熵的运算，其输入y[p_start:p_end]为脑电信号y在第p_start点开始到第p_end点为止的部分，time_length为计算样本熵的每个样本的时间长度，fs为脑电信号的采样率，T为开始采集脑电信号后的设定时间，v为设定参数。