CN106178222A - 基于催眠术的智能辅助睡眠方法和系统 - Google Patents

Abstract

本发明涉及一种基于催眠术的智能辅助睡眠方法和系统，其中所述方法包括：采集用户在睡眠中的生物电信号；根据所述生物电信号识别用户当前的睡眠状态；在用户当前为清醒状态时，播放预设的催眠引导词对用户进行催眠，并检测用户的催眠深度；当用户当前催眠深度达到预定的催眠等级时，播放睡眠指令引导用户进入睡眠。本发明的技术方案，通过识别用户的睡眠状态，再进行相应干预，避免错误干预影响用户睡眠，通过利用催眠内容辅助用户入睡，在设定催眠深度时引导用户入睡，既避免催眠深度过低导致起不到辅助入睡作用的现象，也避免了催眠深度过高带来的风险的情况，有效地提高了辅助睡眠效果。

Description

基于催眠术的智能辅助睡眠方法和系统

技术领域

本发明涉及辅助睡眠技术领域，特别是涉及一种基于催眠术的智能辅助睡眠方法和系统。

背景技术

在睡眠中，人体进行了自我放松及恢复的过程。因此良好的睡眠是保持身体健康的一项基本条件。但是由于工作压力大、生活作息不规律等原因，导致了部分人群的睡眠质量欠佳，表现为失眠、半夜惊醒等。智能辅助睡眠是一种结合现代科技的睡眠方法，当被试者进入催眠状态后，其受暗示性明显提高，能与催眠师保持密切的感应关系，会不加批判地接受其暗示指示。将催眠术应用于辅助睡眠时，当催眠者被催眠师所催眠后，催眠师发出睡眠指令即可使被催眠者进入睡眠状态。与药物干预(安眠药)相比，基于催眠术的辅助睡眠对身体的副作用较小，比较适合日常应用。

目前，市面上也较为常见的基于声音和光照的睡眠干预设备。光照刺激只适合唤醒，对于入睡不仅没有帮助。常见的声音刺激有白噪声和放松音乐等。催眠的声音干预对于不同人来说，效果并不一样，甚至适得其反，由于无法检测用户睡眠和催眠状态，从而影响了辅助睡眠效果。

发明内容

基于此，有必要针对上述问题，提供一种基于催眠术的智能辅助睡眠方法和系统，有效地提高辅助睡眠效果。

一种基于催眠术的智能辅助睡眠方法，包括：

采集用户在睡眠中的生物电信号；

根据所述生物电信号识别用户当前的睡眠状态；

在用户当前为清醒状态时，播放预设的催眠引导词对用户进行催眠，并检测用户的催眠深度；

当用户当前催眠深度达到预定的催眠等级时，播放睡眠指令引导用户进入睡眠。

一种基于催眠术的智能辅助睡眠系统，包括：

采集模块，用于采集用户在睡眠中的生物电信号；

识别模块，用于根据所述生物电信号识别用户当前的睡眠状态；

催眠模块，用于在用户当前为清醒状态时，播放预设的催眠引导词对用户进行催眠，并检测用户的催眠深度；

引导模块，用于当用户当前催眠深度达到预定的催眠等级时，播放睡眠指令引导用户进入睡眠。

上述基于催眠术的智能辅助睡眠方法和系统，通过识别用户的睡眠状态，再进行相应干预，避免错误干预影响用户睡眠，通过利用催眠内容辅助用户入睡，在设定催眠深度时引导用户入睡，既避免催眠深度过低导致起不到辅助入睡作用的现象，也避免了催眠深度过高带来的风险的情况，有效地提高了辅助睡眠效果。

附图说明

图1为本发明的基于催眠术的智能辅助睡眠方法的流程图；

图2为一个示例的智能辅助睡眠算法流程图；

图3为本发明的基于催眠术的智能辅助睡眠系统的结构示意图。

具体实施方式

下面结合附图阐述本发明的基于催眠术的智能辅助睡眠方法和系统的实施例。

参考图1所示，图1为本发明的基于催眠术的智能辅助睡眠方法的流程图，包括：

S10，采集用户在睡眠中的生物电信号；

在本步骤中，主要是通过相关设备检测用户的生物电信号，例如，可以采集脑电信号、眼电信号等生物电信号。

在采集生物信号时，一般是以30s为一帧进行采集，每帧作为一个样本，后续对每帧脑电信号进行分析处理。

S20，根据所述生物电信号识别用户当前的睡眠状态；

在本步骤中，利用采集的生物电信号判断用户当前的睡眠状态，目的是要通过相关技术手段来识别用户当前是否处于睡着状态，如果睡着状态则不宜进行睡眠干预，否则容易起到反作用，影响用户的睡眠，如果处于清醒状态，则需要进行睡眠干预，引导用户进入睡眠状态。

在一个实施例中，利用脑电信号识别用户当前的睡眠状态，则对于根据所述生物电信号识别用户当前的睡眠状态的步骤，可以包括如下过程：

s201，采集用户的脑电信号；

s202，对所述脑电信号进行滤波；其中，所述滤波包括带通滤波并滤除工频干扰；

s203，从滤波后的脑电信号中提取信号特征；

对于上述信号特征，可以包括如下几种当中的的一种或多种组合：(1)脑电信号中δ波频段，θ波频段，α波频段，β波频段的能量在总能量中的比例；(2)在一帧的时间内，脑电信号中δ波频段，θ波频段，α波频段，β波频段能量最大的时间长度；(3)脑电信号的信号基线变化的幅度；(4)脑电信号进行信号小波分解后，小波系数的均值、方差，峭度和斜度；对于小波分解，其分解的层数可以由信号的采样率决定，为了分解出delta波(0～4Hz)，当信号的采样率为128Hz时，可以选择4层分解，当信号的采样率为256Hz时，则可以进行5层分解；

经过小波分解后，对δ频段、θ频段、α频段和β频段的信号进行重构。为了计算的方便，其频段范围分别设置为[0.5Hz,4Hz](δ频段)，[4Hz,8Hz](θ频段)，[8Hz,16Hz](α频段)和[16Hz,32Hz](β频段)。

考虑到超过β频段的信号对睡眠分析没有帮助，因此可以将信号的总能量(p_total)定义为δ频段、θ频段、α频段和β频段信号的能量之和。这四种频段的信号的能量之比利用如下公式进行计算：

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_β分别代表这四种频段的信号的能量在总能量的比例。

由于临床上对睡眠状态进行分析时，可以统计一帧(一般为30秒)内δ波、θ波、α波和β波的持续时间。因此在此可以利用如下公式来计算在各秒内能量占比最大的频段：

$c_{\mathrm{\δ}}=\underset{i=1}{\overset{30}{\Σ}}{f}_{\mathrm{\δ}}^i,f_{\mathrm{\δ}}^i=\left\{\begin{array}{c}1,if{r}_{\mathrm{\δ}}^i=\max (r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\\ 0,if{r}_{\mathrm{\δ}}^i\≠\max (r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\end{array}\right.$

$c_{\mathrm{\θ}}=\underset{i=1}{\overset{30}{\Σ}}{f}_{\mathrm{\θ}}^i,f_{\mathrm{\θ}}^i=\left\{\begin{array}{c}1,if{r}_{\mathrm{\θ}}^i=\max (r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\\ 0,if{r}_{\mathrm{\θ}}^i\≠\max (r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\end{array}\right.$

$c_{\mathrm{\α}}=\underset{i=1}{\overset{30}{\Σ}}{f}_{\mathrm{\α}}^i,f_{\mathrm{\α}}^i=\left\{\begin{array}{c}1,if{r}_{\mathrm{\α}}^i=\max (r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\\ 0,if{r}_{\mathrm{\α}}^i\≠\max (r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\end{array}\right.$

$c_{\mathrm{\β}}=\underset{i=1}{\overset{30}{\Σ}}{f}_{\mathrm{\β}}^i,f_{\mathrm{\β}}^i=\left\{\begin{array}{c}1,if{r}_{\mathrm{\β}}^i=\max (r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\\ 0,if{r}_{\mathrm{\β}}^i\≠\max (r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\end{array}\right.$

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

s204，根据所述信号特征识别用户当前的睡眠状态；

在训练睡眠状态识别模型来识别清醒状态或睡着状态，如支持向量机模型，在训练模型时，采用Grid-test方法选择最优参数，即利用惩罚因子C，和rbf核的参数σ，参数的取值范围分别设置为C：2^-2～2¹²，σ：2^-2～2¹⁰，同时调节这两个参数，将识别率最高的参数作为最优的参数，进而利用测试参数训练出睡眠状态识别模型。

在一个实施例中，识别清醒状态的方法，可以包括如下步骤：

计算所述脑电信号的样本熵，将该样本熵与预先计算的样本熵阈值进行比较，若所述样本熵大于所述样本熵阈值，则判定用户当前处于清醒状态。

其中，样本熵阈值的计算公式可以如下：

$sampen\_thre=\frac{1}{n}\underset{i=1}{\overset{n}{\Σ}}sampen\_{\mathrm{val}}_i+\frac{v}{n}(\underset{i=1}{\overset{n}{\Σ}}sampen\_{\mathrm{val}}_i^2-\underset{i=1}{\overset{n}{\Σ}}sampen\_{\mathrm{val}}_i)$

sampen_val_i＝sampen(y[p_start:p_end])

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

p_end＝t_start+time_length*fs-1

p_end＜T·fs

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

在上述计算方案中，参数v的取值非常重要，通过参数v可以控制识别准确率；因此，为了提高识别准确率，对于参数v取值可以通过如下公式计算：

假设清醒状态时的样本熵的集合X服从标准正态分布，集合X中第i个元素表示为：

$X_i=\frac{sampen\_{\mathrm{val}}_i-u}{\mathrm{\σ}},i=1,...,n$

其中，

此时，x＝v

根据标准正态分布函数的积分：

其中，P(X≤x)表示样本熵的集合X中的取值小于x的概率，由此可以计算，以T＝300s，time_length＝30s为例，当参数v＝2.58时，样本熵集合X的值小于x的概率为99.5％。

另外，在识别清醒状态的方法时，也可以包括如下技术方案：

在用户开始睡眠过程后，采集用户的实时眼电信号，对所述实时眼电信号进行小波分解，并根据设定低频段的小波系数进行信号重建得到眼电信号，根据同一时刻脑电信号和眼电信号的相关性以及眨眼眼电波形的特征，在眼电信号上检测眨眼活动，当检测到眨眼活动时，判定用户当前处于清醒状态。

S30，在用户当前为清醒状态时，播放预设的催眠引导词对用户进行催眠，并检测用户的催眠深度；

上述处理过程，是基于前述技术手段的识别，在用户当前为清醒状态时，播放预设的催眠引导词对用户进行催眠，引导用户进入睡眠，并检测用户的催眠深度，从而监控到目前用户所处的被催眠状态。

进一步地，当识别用户当前的睡眠状态为睡着状态时，暂停催眠干预；如上所述，如果用户处于睡着状态，则不宜进行睡眠干预，否则容易起到反作用，避免影响睡眠质量。

S40，当用户当前催眠深度达到预定的催眠等级时，播放睡眠指令引导用户进入睡眠；

在此，是通过设定催眠等级，来引导用户达到该催眠等级，由于催眠不同深度对用户睡眠质量具有重要影响，因此，催眠必须要达到并确保在一个合适的深度范围内，在通过技术手段检测到设定的睡眠等级时，播放睡眠指令引导用户进入睡眠状态。

在一个实施例中，若用户当前催眠深度未达到预定的催眠等级，则继续播放催眠引导词，并持续检测用户的催眠深度，直至达到预定的催眠等级时，播放睡眠指令引导用户进入睡眠。

对于上述预定的催眠等级，一般选取第三级催眠深度，第三催眠深度能更好引导用户进入睡眠，第三级催眠深度时可以完全控制用户所有的肌肉系统，产生数字阻滞。

如前面所述，可以选取第三级催眠深度，则播放睡眠指令引导用户进入睡眠的步骤，即需要判断用户当前催眠深度是否达到第三级催眠深度；

所述第三级催眠深度判断方法，可以包括如下：

(a)在播放催眠引导词中多次插入靶刺激信息；

(b)将多次靶刺激出现后检测到的第一预设时间(一般取600ms)内的脑电信号相加求平均；

(c)检测第二预设时间范围(一般取300ms～500ms)内是否出现明显的正向波形；若是，判断用户当前催眠深度达到第三级催眠深度，否则判断用户当前催眠深度未达到第三级催眠深度。

综合上述实施例，本发明的技术方案，可以利用催眠技术来帮助用户入睡，适用范围更广，利用生物电信号识别催眠深度，能准确的识别当前催眠状态；实现了智能化的生物反馈功能；在准确的识别催眠状态的基础上，调用合适的催眠内容，避免打扰用户睡眠；并且在中等催眠深度时引导用户入睡，避免催眠深度过高带来的风险，也避免了催眠深度过低无法帮助入睡的作用。

需要说明的是，本发明的技术，除了用于辅助睡眠外，还可以用于午休和日常放松等场景中。

基于上述提供的基于催眠术的智能辅助睡眠方法，在实际应用中，还可以设计如下智能辅助睡眠算法。

参考图2所示，图2为一个示例的智能辅助睡眠算法流程图。

s1，采集用户的生物电信号；

s2，根据生物电信号对用户进行睡眠检测；

s3，判断用户是否处于睡着状态，若是，执行s4，否则执行s5；

s4，暂停干预；

s5，向用户播放催眠引导词；

s6，检测用户的催眠深度；

s7，判断用户是否处于第三级催眠深度，若是，执行s8，否则返回执行s5；

s8，播放睡眠指令，引导用户入睡，转去执行s2。

上述算法形成循环流程，能够应用于相应智能设备上执行，以辅助用户睡眠。

参考图3所示，图3为本发明的基于催眠术的智能辅助睡眠系统的结构示意图，包括：

采集模块10，用于采集用户在睡眠中的生物电信号；

识别模块20，用于根据所述生物电信号识别用户当前的睡眠状态；

催眠模块30，用于在用户当前为清醒状态时，播放预设的催眠引导词对用户进行催眠，并检测用户的催眠深度；

引导模块40，用于当用户当前催眠深度达到预定的催眠等级时，播放睡眠指令引导用户进入睡眠。

本发明的基于催眠术的智能辅助睡眠系统与本发明的基于催眠术的智能辅助睡眠方法一一对应，在上述基于催眠术的智能辅助睡眠方法的实施例阐述的技术特征及其有益效果均适用于基于催眠术的智能辅助睡眠系统的实施例中，特此声明。

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

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

Claims (10)

1.一种基于催眠术的智能辅助睡眠方法，其特征在于，包括：

采集用户在睡眠中的生物电信号；

根据所述生物电信号识别用户当前的睡眠状态；

在用户当前为清醒状态时，播放预设的催眠引导词对用户进行催眠，并检测用户的催眠深度；

当用户当前催眠深度达到预定的催眠等级时，播放睡眠指令引导用户进入睡眠。

2.根据权利要求1所述的基于催眠术的智能辅助睡眠方法，其特征在于，所述生物电信号包括脑电信号或眼电信号。

3.根据权利要求1所述的基于催眠术的智能辅助睡眠方法，其特征在于，还包括：若用户当前催眠深度未达到预定的催眠等级，则继续播放催眠引导词，并持续检测用户的催眠深度，直至达到预定的催眠等级时，播放睡眠指令引导用户进入睡眠。

4.根据权利要求2所述的基于催眠术的智能辅助睡眠方法，其特征在于，根据所述生物电信号识别用户当前的睡眠状态的步骤包括：

采集用户的脑电信号；

对所述脑电信号进行滤波；其中，所述滤波包括带通滤波并滤除工频干扰；

从滤波后的脑电信号中提取信号特征；

根据所述信号特征识别用户当前的睡眠状态。

5.根据权利要求1所述的基于催眠术的智能辅助睡眠方法，其特征在于，所述信号特征包括如下当中的一种或多种组合：

脑电信号中δ波频段，θ波频段，α波频段，β波频段的能量在总能量中的比例；

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

脑电信号的信号基线变化的幅度；

脑电信号进行信号小波分解后，小波系数的均值、方差，峭度和斜度。

6.根据权利要求5所述的基于催眠术的智能辅助睡眠方法，其特征在于，所述脑电信号中δ波频段，θ波频段，α波频段，β波频段的能量在总能量中的比例表示为如下公式：

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_β分别代表δ频段、θ频段、α频段和β频段的信号的能量在总能量的比例。

7.根据权利要求5所述的基于催眠术的智能辅助睡眠方法，其特征在于，所述在一帧的时间内，脑电信号中δ波频段，θ波频段，α波频段，β波频段能量最大的时间长度表示为如下公式：

$c_{\mathrm{\δ}}=\underset{i=1}{\overset{30}{\Σ}}{f}_{\mathrm{\δ}}^i,f_{\mathrm{\δ}}^i=\left\{\begin{array}{c}1,if{r}_{\mathrm{\δ}}^i=max(r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\\ 0,if{r}_{\mathrm{\δ}}^i\≠max(r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\end{array}\right.$

$c_{\mathrm{\θ}}=\underset{i=1}{\overset{30}{\Σ}}{f}_{\mathrm{\θ}}^i,f_{\mathrm{\θ}}^i=\left\{\begin{array}{c}1,if{r}_{\mathrm{\θ}}^i=max(r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\\ 0,if{r}_{\mathrm{\θ}}^i\≠max(r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\end{array}\right.$

$c_{\mathrm{\α}}=\underset{i=1}{\overset{30}{\Σ}}{f}_{\mathrm{\α}}^i,f_{\mathrm{\α}}^i=\left\{\begin{array}{c}1,if{r}_{\mathrm{\α}}^i=max(r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\\ 0,if{r}_{\mathrm{\α}}^i\≠max(r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\end{array}\right.$

$c_{\mathrm{\β}}=\underset{i=1}{\overset{30}{\Σ}}{f}_{\mathrm{\β}}^i,f_{\mathrm{\β}}^i=\left\{\begin{array}{c}1,if{r}_{\mathrm{\β}}^i=max(r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\\ 0,if{r}_{\mathrm{\β}}^i\≠max(r_{\mathrm{\δ}}^i,r_{\mathrm{\θ}}^i,r_{\mathrm{\α}}^i,r_{\mathrm{\β}}^i)\end{array}\right.$

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

8.根据权利要求4所述的基于催眠术的智能辅助睡眠方法，其特征在于，所述根据所述信号特征识别用户当前的睡眠状态的步骤包括：

计算所述脑电信号的样本熵，将该样本熵与预先计算的样本熵阈值进行比较，若所述样本熵大于所述样本熵阈值，则判定用户当前处于清醒状态。

9.根据权利要求1所述的基于催眠术的智能辅助睡眠方法，其特征在于，所述预定的催眠等级为第三级催眠深度；

所述播放睡眠指令引导用户进入睡眠的步骤前，还包括：判断用户当前催眠深度是否达到第三级催眠深度；

所述判断方法包括：

在播放催眠引导词中多次插入靶刺激信息；

将多次靶刺激出现后检测到的第一预设时间内的脑电信号相加求平均；

检测第二预设时间范围内是否出现明显的正向波形；若是，判断用户当前催眠深度达到第三级催眠深度，否则判断用户当前催眠深度未达到第三级催眠深度。

10.一种基于催眠术的智能辅助睡眠系统，其特征在于，包括：

采集模块，用于采集用户在睡眠中的生物电信号；

识别模块，用于根据所述生物电信号识别用户当前的睡眠状态；

催眠模块，用于在用户当前为清醒状态时，播放预设的催眠引导词对用户进行催眠，并检测用户的催眠深度；

引导模块，用于当用户当前催眠深度达到预定的催眠等级时，播放睡眠指令引导用户进入睡眠。