CN104720748A - 一种睡眠阶段确定方法和系统 - Google Patents

一种睡眠阶段确定方法和系统 Download PDF

Info

Publication number
CN104720748A
CN104720748A CN201310722893.6A CN201310722893A CN104720748A CN 104720748 A CN104720748 A CN 104720748A CN 201310722893 A CN201310722893 A CN 201310722893A CN 104720748 A CN104720748 A CN 104720748A
Authority
CN
China
Prior art keywords
time slice
time
characteristics parameter
ecg characteristics
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201310722893.6A
Other languages
English (en)
Other versions
CN104720748B (zh
Inventor
王俊艳
张志鹏
许利群
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Mobile Communications Group Co Ltd
Original Assignee
China Mobile Communications Group Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Mobile Communications Group Co Ltd filed Critical China Mobile Communications Group Co Ltd
Priority to CN201310722893.6A priority Critical patent/CN104720748B/zh
Publication of CN104720748A publication Critical patent/CN104720748A/zh
Application granted granted Critical
Publication of CN104720748B publication Critical patent/CN104720748B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4806Sleep evaluation
    • A61B5/4812Detecting sleep stages or cycles

Abstract

本发明实施例提供一种睡眠阶段确定方法和系统,可以确定心电信号的第一心电特征参数,通过判断在一个时间片段内,该第一心电特征参数是否满足睡眠阶段为醒或REM时的特征,即满足设定条件,来确定该时间片段是否处于醒或REM睡眠阶段。通过醒或REM睡眠阶段、以及其他类型的睡眠阶段的区分,实现睡眠宏结构的确定。而通过规则判断来实现睡眠阶段的确定,保证了睡眠阶段确定的准确性。另外,还可以结合睡眠阶段分类器实现睡眠阶段的综合判断,利用睡眠阶段分类器对其他类型的睡眠阶段进行浅度睡眠阶段和深度睡眠阶段的分类,确定睡眠的微结构,进一步提高睡眠阶段判断的准确性。

Description

一种睡眠阶段确定方法和系统
技术领域
本发明涉及通信领域,尤其涉及一种睡眠阶段确定方法和系统。
背景技术
人的一生有近1/3的时间在睡眠中度过,睡眠的好坏直接关系到人体的记忆、学习、工作和免疫力等多个方面。随着时代的发展,人们的生活节奏越来越快,压力越来越大,很多人出现了睡眠问题,睡眠监测应运而生。
在睡眠监测中,人的睡眠通常用6个阶段的连续来描述,这6个阶段通常包括:
1、清醒阶段,可以记为醒。在这个阶段,睡眠者的意识处于清醒状态。
2、入睡阶段,可以记为慢速眼球运动睡眠(non-rapid eye movements,NREM)1。在NREM1中,睡眠者的眼睛是闭上的,不过如果在此阶段被唤醒,睡眠者可能感觉他还没有睡着。
3、浅睡阶段,可以记为NREM2。在NREM2中,睡眠者的身体准备进入深睡。
4、深睡阶段,可以记为NREM3,睡眠者进入深度睡眠状态。
5、延续深睡阶段,可以记为NREM4,且睡眠者在NREM4中比在NREM3中睡得更深。
6、快速眼球运动阶段,可以记为快速眼球运动睡眠(rapid eye movements,REM)。在REM中,与NREM1至NREM4的生理状态有所区别,睡眠者的眼球在此阶段时会快速移动。
目前临床上监测睡眠阶段的方法中,典型的方法是采用多导睡眠仪(Polysomnography,PSG)采集睡眠期间的生理信号,包括脑电波(EEG)、眼动(EOG)、肌肉运动(EMG)、心电(ECG)、血氧饱和度(SpO2)和呼吸信号。睡眠阶段判断主要基于对EEG的分析,并可以利用ECG和EEG进行辅助判断。然而PSG操作复杂,仅适用于医院的临床研究,在家庭监测中难以推广。
因此,目前通常将睡眠用三个阶段的连续来描述,仅通过ECG进行睡眠阶段确定,使得即使不使用PSG,也可以进行睡眠阶段判断,达到简化睡眠阶段确定的目的。这三个睡眠阶段包括:
1、醒或REM。由于在利用心电或者心率变异率进行睡眠描述时,醒和REM的表现是几乎一样的,因此,将睡眠用三个阶段的连续来描述时,将清醒阶段和REM划分为一个阶段。
2、浅度睡眠阶段。将睡眠用三个阶段的连续来描述时,将入睡阶段和浅睡阶段划分为一个阶段,并统一称为浅度睡眠阶段。
3、深度睡眠阶段。将睡眠用三个阶段的连续来描述时,将深睡阶段和延续深睡阶段划分为一个阶段,并统一称为深度睡眠阶段。
通过ECG进行睡眠阶段判断时,主要包括心电特征提取和睡眠阶段分类器两个部分。其中,常用的心电特征主要是心率变异率(Heart rate variability,HRV)的特征,包括心率变异率的时域特征和频域特征。常用的睡眠阶段分类器主要是统计分类器,如支持向量机分类器,或隐马尔可夫模型分类器。通过ECG进行睡眠阶段判断的示意图可以如图1所示。首先,在训练模块中,对ECG训练集进行心电特征提取,利用提取出的心电特征进行分类器训练,得到睡眠阶段分类器;然后在使用模块中,提取待处理ECG数据的心电特征,并根据提取出的心电特征,采用通过训练模块得到的睡眠阶段分类器确定睡眠阶段。
目前通过ECG进行睡眠阶段确定的方案中,睡眠阶段的确定结果由睡眠阶段分类器决定。而睡眠阶段分类器完全是利用训练集训练得到的。睡眠阶段分类器的训练是一个有监督学习的过程,需要具有标签的样本进行训练,对于基于统计学习的睡眠阶段分类器,如支持向量机或者隐马尔可夫模型,其分类能力依赖于训练集,这样可能引起下面的问题。每个人睡眠期间的心率状况是不同的,而且,由于睡眠环境、身体状态等不同因素,同一个人的不同日期的睡眠期间的心率状况也是不同的。为了适应不同的用户,需要大量的训练样本才可能覆盖整个样本空间。而睡眠数据的获取是不容易的,而且人工标注的工作量也很大。因此,难以获得足够的训练样本。而在有限的训练样本的情况下,睡眠阶段分类器很难具有好的普适性,且睡眠阶段确定的精度难以得到保证。
发明内容
本发明实施例提供一种睡眠阶段确定方法和系统,用于提高睡眠阶段确定的准确性。
一种睡眠阶段确定方法,所述方法包括:
确定指定时间段内心电信号的第一心电特征参数;
在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,若确定所述第一心电特征参数满足设定条件,确定该时间片段对应的睡眠阶段为醒或快速眼球运动睡眠REM,否则,确定该时间片段对应的睡眠阶段非醒或REM;
其中,所述时间片段是根据不大于设定时长,将所述指定时间段依次划分得到的。
一种睡眠阶段确定系统,所述系统包括:
确定模块,用于确定指定时间段内心电信号的第一心电特征参数;
分析模块,用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,若确定所述第一心电特征参数满足设定条件,确定该时间片段对应的睡眠阶段为醒或快速眼球运动睡眠REM,否则,确定该时间片段对应的睡眠阶段非醒或REM;其中,所述时间片段是根据不大于设定时长,将所述指定时间段依次划分得到的。
根据本发明实施例提供的方案,可以确定心电信号的第一心电特征参数,通过判断在一个时间片段内,该第一心电特征参数是否满足睡眠阶段为醒或REM时的特征,即满足设定条件,来确定该时间片段是否处于醒或REM睡眠阶段。通过区分醒或REM睡眠阶段、其它类型的睡眠阶段,实现了睡眠宏结构的确定。而通过规则判断来识别醒或REM睡眠阶段,保证了睡眠阶段确定的准确性。
附图说明
图1为现有技术提供的通过ECG进行睡眠阶段判断的示意图;
图2为本发明实施例一提供的睡眠阶段确定方法的步骤流程图;
图3为本发明实施例二提供的睡眠阶段确定方法的步骤示意图;
图4为本发明实施例三提供的睡眠阶段确定系统的结构示意图;
图5为本发明实施例四提供的可穿戴式节点的结构示意图。
具体实施方式
本案发明人研究发现,将睡眠用三个阶段的连续来描述时,醒或REM与其他两个睡眠阶段具有明显的区别。其中至少包括:(1)、睡眠阶段为REM和醒的平均心率都比睡眠阶段为NREM1~4时快;(2)、睡眠阶段为REM和醒时心率变化的标准差大于睡眠阶段为NREM1~4时。因此,可以利用这种区别来进行醒或REM这个睡眠阶段的识别。通过醒或REM睡眠阶段、以及其他类型的睡眠阶段的区分,实现睡眠宏结构的确定。进一步的,可以利用睡眠阶段分类器进行浅度睡眠阶段和深度睡眠阶段的分类,结合睡眠阶段分类器实现睡眠阶段的综合判断。通过从其他类型的睡眠阶段中实现深度睡眠阶段和浅度睡眠阶段的区分,实现睡眠微结构的确定。
且本案发明人还发现,睡眠阶段具有周期性的特点,因此完全可以将相邻的、属于同一类型睡眠阶段的时段合并在一起,统一进行该合并后时段的睡眠阶段的判断,从而进一步减少睡眠阶段的判断次数,提高睡眠阶段的判断速度。
以下结合说明书附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。并且在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。
实施例一、
本发明实施例一提供一种睡眠阶段确定方法,该方法的步骤流程可以如图2所示,包括:
步骤101、确定心电信号的第一心电特征参数。
在本步骤中,可以确定指定时间段内心电信号的第一心电特征参数。所述第一心电特征参数可以但不限于理解为心率变异率或心率等信号。
所述心电信号可以是接收到的预先采集的心电信号,也可以是在本步骤中采集到的心电信号。
需要说明的是,为了后续描述所述指定时段对应的睡眠阶段,可以根据不大于设定时长,将所述指定时段依次划分得到多个时间片段。具体的,可以将所述指定时段均匀划分得到多个时间片段。例如,可以将所述指定时段均匀划分为多个时间片段,每个时间片段的长度为30秒。
步骤102、确定睡眠阶段。
在本步骤中,在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,若确定所述第一心电特征参数满足设定条件,确定该时间片段对应的睡眠阶段为醒或REM,否则,确定该时间片段对应的睡眠阶段非醒或REM。
优选的,在所述第一心电特征参数为心率变异率时,可以根据睡眠阶段为REM和醒的平均心率都比睡眠阶段为NREM1~4时快,以及睡眠阶段为REM和醒时心率变化的标准差大于睡眠阶段为NREM1~4时,通过判断在所述指定时段的一个时间片段中,所述心率变异率的均值是否较小,复杂度是否较高来确定该时间片段对应的睡眠阶段是否为醒或REM。
具体的,可以通过在所述指定时段的一个时间片段中,确定所述心率变异率的均值是否小于第一设定值;并在该时间片段中,确定所述心率变异率的方差是否大于第二设定值和/或所述心率变异率的样本熵是否大于第三设定值,来确定该时间片段对应的睡眠阶段是否为醒或REM。
当然,任何可以用于识别醒或REM睡眠阶段的特征均可以单独或组合用于设定条件的设定。例如,在所述第一心电特征参数为心率时,可以根据睡眠阶段为REM和醒的平均心率都比睡眠阶段为NREM1~4时快,通过判断在所述指定时段的一个时间片段中,心率的均值是否较小,来确定该时间片段对应的睡眠阶段是否为醒或REM。
进一步的,根据睡眠阶段具有周期性的特点,因此完全可以将相邻的、属于同一类型睡眠阶段的时间片段进行合并,并在后续针对合并后的时间片段,进行睡眠阶段的判断。因此,在步骤102之前,所述方法还可以进一步包括步骤102’。在图2中,以步骤102’在步骤101之后执行为例进行说明。进行时间片段合并可以是根据心电信号进行的,也可以是根据心电信号的第二心电特征参数进行的,还可以是根据心电信号以及心电信号的第二心电特征参数进行的。
步骤102’、进行时间片段的合并。
以根据心电信号进行时间片段的合并为例,在本步骤中,可以根据一个时间片段内的所述心电信号,确定该时间片段对应的第一数据,根据所述第一数据确定该时间片段的第一特征数据,并将所述第一特征数据差异不大于第一设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第一特征数据差异均大于所述第一设定阈值。
以根据心电信号的第二心电特征参数进行时间片段的合并为例,在本步骤中,可以确定所述心电信号的第二心电特征参数;根据一个时间片段内的第二心电特征参数,确定该时间片段对应的第二数据,根据所述第二数据确定该时间片段的第二特征数据,并将所述第二特征数据差异不大于第二设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第二特征数据差异均大于所述第二设定阈值。
更优的,为了提高时间片段合并的精度,避免将不属于同一睡眠阶段的时间片段划分到一起,在本步骤中,还可以根据心电信号以及心电信号的第二心电特征参数进行时间片段的合并。具体的,可以分别根据心电信号的第二心电特征参数进行时间片段的合并,以及根据心电信号进行时间片段的合并,然后对两组合并结果进行综合,取更细致的合并结果作为最终的合并结果。对合并结果进行综合时,可以确定任意两个相邻时间片段的第一特征数据差异均大于第一设定阈值的每个时间片段的端点时间,以及任意两个相邻时间片段的第二特征数据差异均大于第二设定阈值的每个时间片段的端点时间;将所述端点时间按照从大到小或从小到大排序,利用任意两个相邻端点时间确定一个合并后的时间片段。
需要说明的是,在本实施例中,还可以进一步结合睡眠阶段分类器实现睡眠阶段的确定。在确定该时间片段对应的睡眠阶段非醒或REM时,步骤102之后,还可以进一步包括步骤103:
步骤103、进一步确定睡眠阶段。
在本步骤中,可以利用预先训练得到的睡眠阶段分类器确定所述时间片段对应的睡眠阶段为深度睡眠阶段或浅度睡眠阶段。
需要说明的是,目前的研究发现,随着睡眠时间的增加,心率有变慢的趋势的特点。为了消除心率变异率数据中跟睡眠时间有关的信息,可以对心率变异率数据进行去趋势分析,优化根据心率变异率数据进行的操作,进一步提高睡眠阶段确定的精度。因此,在本实施例中,在第一心电特征参数或第二心电特征参数为心率变异率时,可以通过以下方式实现心率变异率的确定。
提取指定时间段内心电信号的R波位置,根据所述R波位置确定所述心电信号对应的原始心率变异率数据;对所述原始心率变异率数据进行线性拟合,得到所述原始心率变异率数据的变化趋势数据;确定所述原始心率变异率数据与所述变化趋势数据的差值,将所述差值作为确定出的所述心电信号的心率变异率。具体的,在本实施例中,可以但不限于采用最小二乘法或梯度下降法进行线性拟合。
下面以利用心电信号以及心电信号的第二心电特征参数进行时间片段的合并,且第一心电特征参数和第二心电特征参数均为心率变异率为例,对本发明实施例一提供的方案进行说明。
实施例二、
本发明实施例二提供一种睡眠阶段确定方法,该方法的步骤示意图可以如图3所示,包括:
步骤201、采集心电信号,确定心率变异率。
在本步骤中,可以采集指定时间段内的心电信号,并根据采集到的心电信号得到去趋势分析后的心率变异率,消除心率变异率数据中跟睡眠时间有关的信息,以降低后续根据心率变异率进行分段和与设定条件匹配时的复杂度,提高分段和与设定条件匹配的精度,从而提高确定睡眠阶段的精度。
具体的,在本步骤中,可以提取指定时间段内采集到的心电信号的R波位置,根据提取出的R波位置确定所述心电信号对应的原始心率变异率数据,并对原始心率变异率数据进行去趋势分析得到后续使用的心率变异率。
步骤202、进行时间片段的合并。
由于第二心电特征参数和第一心电特征参数均为心率变异率,在本步骤中,可以直接分别根据心电信号和心率变异率进行时间片段的合并。
根据心电信号进行时间片段的合并时:
假设未合并前,所述指定时间段包括n个时间片段,第i个时间片段的第一数据可以表示为vi和ti,其中,ti可以理解为第i个时间片段的时长,vi可以理解为第i个时间片段对应的一个特征参数,则所述指定时间段n个时间片段的第一数据集合NNArr可以表示如下:
NNArr=(v1,t1),(v2,t2),...(vn-i+1,tn-i+1),...,(vn,tn)
第i个时间片段的第一特征数据ui和ki可以根据vi和ti得到。其中,ui可以理解为第i个时间片段中数据的均值,ki可以理解为第i个时间片段中数据的斜率。
第一特征数据差异dij可以设定为dij=|ui-uj|+λ|κij|,其中,λ是斜率影响因子,且λ不小于0。第一设定阈值可以用ThTFA表示,可以利用dij≤ThTFA来进行时间片段的合并,直至任意两个相邻时间片段的dij均大于ThTFA
较优的,第一数据中的vi可以为心电信号稳定性特征参数,确定vi的过程就是对心电信号进行稳定性分析的过程。第i个时间片段的vi可以用相邻的两次心跳的距离的平均值表示:
s ( i ) = Σ j = 1 n - 1 D ( beat i , beat i + 1 ) n - 1
其中D(beati,beati+1)是相邻两次心跳的距离,可以定义为两个点集之间的距离,如hausdorf距离或动态时间匹配(DTW)后的距离。
根据心率变异率进行时间片段的合并的过程与根据心电信号进行时间片段的合并的过程类似,在此不再赘述。
根据上述两种方法得到两组对时间片段的合并结果后,对合并结果进行综合,取更细致的合并结果作为最终的合并结果。
步骤203、确定睡眠阶段。
在本步骤中,可以确定心率变异率的特征参数,与设定条件进行匹配。在本实施例中,心率变异率的特征参数可以为心率变异率的均值和方差,设定条件可以为一个合并后时间片段内,心率变异率的均值M(i)小于与该时间片段直接相邻的任意一个时间片段内,心率变异率的均值(M(i+1)和M(i-1))(可以理解为实施例一中涉及的一个合并后时间片段内,心率变异率的均值小于第一设定值,可以理解为一个合并后时间片段内,心率变异率的均值小于与该时间片段直接相邻的任意一个时间片段内,心率变异率的均值的较小值),且心率变异率的方差Std大于第二设定值Th。
如果确定心率变异率的特征参数,与设定条件匹配,即确定心率变异率满足设定条件,则可以确定该时间片段对应的睡眠阶段为醒或REM,否则,确定该时间片段对应的睡眠阶段非醒或REM。
步骤204、进一步确定睡眠阶段。
假设步骤203中确定出一个时间片段对应的睡眠阶段非醒或REM,则在本步骤中,可以利用睡眠阶段分类器进一步确定该时间片段对应的睡眠阶段为深度睡眠阶段还是浅度睡眠阶段。
需要说明的是,步骤204中涉及的睡眠阶段分类器可以利用步骤203中确定出的心率变异率的特征参数,进行进一步分类。
与本发明实施例一基于同一发明构思,提供以下的系统和装置。
实施例三、
本发明实施例三提供一种睡眠阶段确定系统,该系统的结构示意图可以如图4所示,包括确定模块11和分析模块12,其中:
确定模块11,用于确定指定时间段内心电信号的第一心电特征参数;
分析模块12,用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,若确定所述第一心电特征参数满足设定条件,确定该时间片段对应的睡眠阶段为醒或快速眼球运动睡眠REM,否则,确定该时间片段对应的睡眠阶段非醒或REM;其中,所述时间片段是根据不大于设定时长,将所述指定时间段依次划分得到的。
所述分析模块12,具体用于确定模块11确定出的所述第一心电特征参数为心率变异率时,在所述指定时段的一个时间片段中,确定所述心率变异率的均值是否小于第一设定值;并在该时间片段中,确定所述心率变异率的方差是否大于第二设定值和/或所述心率变异率的样本熵是否大于第三设定值。
所述分析模块12,还用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件之前,根据一个时间片段内的所述心电信号,确定该时间片段对应的第一数据,根据所述第一数据确定该时间片段的第一特征数据,并将所述第一特征数据差异不大于第一设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第一特征数据差异均大于所述第一设定阈值;
所述分析模块12,用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,具体包括:在一个合并后得到的时间片段中,确定所述第一心电特征参数是否满足设定条件。
所述分析模块12,还用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件之前,确定所述心电信号的第二心电特征参数;根据一个时间片段内的第二心电特征参数,确定该时间片段对应的第二数据,根据所述第二数据确定该时间片段的第二特征数据,并将所述第二特征数据差异不大于第二设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第二特征数据差异均大于所述第二设定阈值;
所述分析模块12,用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,具体包括:在一个合并后得到的时间片段中,确定所述第一心电特征参数是否满足设定条件。
所述分析模块12,还用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件之前,确定所述心电信号的第二心电特征参数;
根据一个时间片段内的第二心电特征参数,确定该时间片段对应的第二数据,根据所述第二数据确定该时间片段的第二特征数据,并将所述第二特征数据差异不大于第二设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第二特征数据差异均大于所述第二设定阈值;以及,
根据一个时间片段内的所述心电信号,确定该时间片段对应的第一数据,根据所述第一数据确定该时间片段的第一特征数据,并将所述第一特征数据差异不大于第一设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第一特征数据差异均大于所述第一设定阈值;
确定任意两个相邻时间片段的第一特征数据差异均大于第一设定阈值的每个时间片段的端点时间,以及任意两个相邻时间片段的第二特征数据差异均大于第二设定阈值的每个时间片段的端点时间;
将所述端点时间按照从大到小或从小到大排序,利用任意两个相邻端点时间确定一个合并后的时间片段;
所述分析模块12,用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,具体包括:在一个所述利用任意两个相邻端点时间确定出的一个合并后的时间片段中,确定所述第一心电特征参数是否满足设定条件。
所述分析模块12,还用于在确定一个时间片段对应的睡眠阶段非醒或REM时,利用预先训练得到的睡眠阶段分类器确定所述时间片段对应的睡眠阶段为深度睡眠阶段或浅度睡眠阶段。
所述确定模块11,具体用于提取指定时间段内心电信号的R波位置,根据所述R波位置确定所述心电信号对应的原始心率变异率数据;对所述原始心率变异率数据进行线性拟合,得到所述原始心率变异率数据的变化趋势数据;确定所述原始心率变异率数据与所述变化趋势数据的差值,将所述差值作为确定出的所述心电信号的心率变异率。
进一步的,所述系统还可以进一步包括显示模块13,用于显示分析模块12确定出的所述指定时段中,每个时间片段对应的睡眠阶段,实现睡眠阶段信息的直观显示。
另外,本发明实施例还提供了一套可记录整晚心电信号的可穿戴式节点,具有操作简单、佩戴舒适的特点,可以广泛应用于家庭睡眠监测。下面通过实施例四进行说明。
实施例四、
本发明实施例四提供一种可穿戴式节点,该节点的结构示意图可以如图5所示,包括低通滤波模块21,低噪声放大模块22,处理器模块23,电源模块24,电源控制模块25,存储模块26和读写模块27,其中:
从睡眠者左右胳膊分别采集电信号,可以分别表示为RA和LA,低通滤波模块21对采集到的电信号进行低通滤波,以及合并RLD电极后,低噪声放大模块22对从低通滤波模块21接收到的信号进行低噪声放大,得到心电信号。低功耗的处理器模块23可以进行控制和信号质量分析,将从低噪声放大模块22得到的心电信号经读写模块27写入存储模块26。读写模块27还可以将得到的心电信号经通用串行总线(USB)口输出。电源模块24可以为可充电锂电池,保证该节点的工作时间可以达到一周。电源控制模块25可以包括电源管理电路,实现多级休眠以降低功耗,并可以与读写模块27连接。
进一步的,该节点还可以包括加速度采集模块28。因此,本实施例提供的节点不仅可以采集心电信号,还可以通过加速度采集模块28采集运动加速度信号,采集到的加速度信号可以与心电信号共同用于睡眠阶段的识别。处理器模块23可以将从加速度采集模块28得到的运动加速度信号,经读写模块27写入存储模块26。读写模块27还可以将得到的运动加速度信号经通用串行总线(USB)口输出。
本实施例中,低通滤波模块21可以包括ADS1194芯片,处理器模块23可以包括MSP430单片机,加速度采集模块28可以包括三轴数字加速度传感器ADXL345,存储模块26可以包括安全数码(SD)卡,读写模块27可以包括读卡机控制器GL896。
低噪声放大模块22与处理器模块23之间,加速度采集模块28与处理器模块23之间,处理器模块23与读写模块27之间,均可以通过串行外围设备接口(SPI)进行通信。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本申请的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例做出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本申请范围的所有变更和修改。
显然,本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样,倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内,则本申请也意图包含这些改动和变型在内。

Claims (14)

1.一种睡眠阶段确定方法,其特征在于,所述方法包括:
确定指定时间段内心电信号的第一心电特征参数;
在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,若确定所述第一心电特征参数满足设定条件,确定该时间片段对应的睡眠阶段为醒或快速眼球运动睡眠REM,否则,确定该时间片段对应的睡眠阶段非醒或REM;
其中,所述时间片段是根据不大于设定时长,将所述指定时间段依次划分得到的。
2.如权利要求1所述的方法,其特征在于,所述第一心电特征参数为心率变异率;
在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,具体包括:
在所述指定时段的一个时间片段中,确定所述心率变异率的均值是否小于第一设定值;并
在该时间片段中,确定所述心率变异率的方差是否大于第二设定值和/或所述心率变异率的样本熵是否大于第三设定值。
3.如权利要求1或2所述的方法,其特征在于,在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件之前,所述方法还包括:
根据一个时间片段内的所述心电信号,确定该时间片段对应的第一数据,根据所述第一数据确定该时间片段的第一特征数据,并将所述第一特征数据差异不大于第一设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第一特征数据差异均大于所述第一设定阈值;
在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,具体包括:
在一个合并后得到的时间片段中,确定所述第一心电特征参数是否满足设定条件。
4.如权利要求1或2所述的方法,其特征在于,在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件之前,所述方法还包括:
确定所述心电信号的第二心电特征参数;
根据一个时间片段内的第二心电特征参数,确定该时间片段对应的第二数据,根据所述第二数据确定该时间片段的第二特征数据,并将所述第二特征数据差异不大于第二设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第二特征数据差异均大于所述第二设定阈值;
在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,具体包括:
在一个合并后得到的时间片段中,确定所述第一心电特征参数是否满足设定条件。
5.如权利要求1或2所述的方法,其特征在于,在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件之前,所述方法还包括:
确定所述心电信号的第二心电特征参数;
根据一个时间片段内的第二心电特征参数,确定该时间片段对应的第二数据,根据所述第二数据确定该时间片段的第二特征数据,并将所述第二特征数据差异不大于第二设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第二特征数据差异均大于所述第二设定阈值;以及,
根据一个时间片段内的所述心电信号,确定该时间片段对应的第一数据,根据所述第一数据确定该时间片段的第一特征数据,并将所述第一特征数据差异不大于第一设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第一特征数据差异均大于所述第一设定阈值;
确定任意两个相邻时间片段的第一特征数据差异均大于第一设定阈值的每个时间片段的端点时间,以及任意两个相邻时间片段的第二特征数据差异均大于第二设定阈值的每个时间片段的端点时间;
将所述端点时间按照从大到小或从小到大排序,利用任意两个相邻端点时间确定一个合并后的时间片段;
在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,具体包括:
在一个所述利用任意两个相邻端点时间确定出的一个合并后的时间片段中,确定所述第一心电特征参数是否满足设定条件。
6.如权利要求1或2所述的方法,其特征在于,在确定一个时间片段对应的睡眠阶段非醒或REM时,所述方法还包括:
利用预先训练得到的睡眠阶段分类器确定所述时间片段对应的睡眠阶段为深度睡眠阶段或浅度睡眠阶段。
7.如权利要求1或2所述的方法,其特征在于,所述第一心电特征参数为心率变异率;
则,确定指定时间段内心电信号的第一心电特征参数,具体包括:
提取指定时间段内心电信号的R波位置,根据所述R波位置确定所述心电信号对应的原始心率变异率数据;
对所述原始心率变异率数据进行线性拟合,得到所述原始心率变异率数据的变化趋势数据;
确定所述原始心率变异率数据与所述变化趋势数据的差值,将所述差值作为确定出的所述心电信号的心率变异率。
8.一种睡眠阶段确定系统,其特征在于,所述系统包括:
确定模块,用于确定指定时间段内心电信号的第一心电特征参数;
分析模块,用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,若确定所述第一心电特征参数满足设定条件,确定该时间片段对应的睡眠阶段为醒或快速眼球运动睡眠REM,否则,确定该时间片段对应的睡眠阶段非醒或REM;其中,所述时间片段是根据不大于设定时长,将所述指定时间段依次划分得到的。
9.如权利要求8所述的系统,其特征在于,所述分析模块,具体用于确定模块确定出的所述第一心电特征参数为心率变异率时,在所述指定时段的一个时间片段中,确定所述心率变异率的均值是否小于第一设定值;并在该时间片段中,确定所述心率变异率的方差是否大于第二设定值和/或所述心率变异率的样本熵是否大于第三设定值。
10.如权利要求8或9所述的系统,其特征在于,所述分析模块,还用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件之前,根据一个时间片段内的所述心电信号,确定该时间片段对应的第一数据,根据所述第一数据确定该时间片段的第一特征数据,并将所述第一特征数据差异不大于第一设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第一特征数据差异均大于所述第一设定阈值;
所述分析模块,用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,具体包括:在一个合并后得到的时间片段中,确定所述第一心电特征参数是否满足设定条件。
11.如权利要求8或9所述的系统,其特征在于,所述分析模块,还用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件之前,确定所述心电信号的第二心电特征参数;根据一个时间片段内的第二心电特征参数,确定该时间片段对应的第二数据,根据所述第二数据确定该时间片段的第二特征数据,并将所述第二特征数据差异不大于第二设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第二特征数据差异均大于所述第二设定阈值;
所述分析模块,用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,具体包括:在一个合并后得到的时间片段中,确定所述第一心电特征参数是否满足设定条件。
12.如权利要求8或9所述的系统,其特征在于,所述分析模块,还用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件之前,确定所述心电信号的第二心电特征参数;
根据一个时间片段内的第二心电特征参数,确定该时间片段对应的第二数据,根据所述第二数据确定该时间片段的第二特征数据,并将所述第二特征数据差异不大于第二设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第二特征数据差异均大于所述第二设定阈值;以及,
根据一个时间片段内的所述心电信号,确定该时间片段对应的第一数据,根据所述第一数据确定该时间片段的第一特征数据,并将所述第一特征数据差异不大于第一设定阈值的两个相邻时间片段合并为一个时间片段,重复执行本步骤,直至任意两个相邻时间片段的第一特征数据差异均大于所述第一设定阈值;
确定任意两个相邻时间片段的第一特征数据差异均大于第一设定阈值的每个时间片段的端点时间,以及任意两个相邻时间片段的第二特征数据差异均大于第二设定阈值的每个时间片段的端点时间;
将所述端点时间按照从大到小或从小到大排序,利用任意两个相邻端点时间确定一个合并后的时间片段;
所述分析模块,用于在所述指定时段的一个时间片段中,确定所述第一心电特征参数是否满足设定条件,具体包括:在一个所述利用任意两个相邻端点时间确定出的一个合并后的时间片段中,确定所述第一心电特征参数是否满足设定条件。
13.如权利要求8或9所述的系统,其特征在于,所述分析模块,还用于在确定一个时间片段对应的睡眠阶段非醒或REM时,利用预先训练得到的睡眠阶段分类器确定所述时间片段对应的睡眠阶段为深度睡眠阶段或浅度睡眠阶段。
14.如权利要求8或9所述的系统,其特征在于,所述确定模块,具体用于提取指定时间段内心电信号的R波位置,根据所述R波位置确定所述心电信号对应的原始心率变异率数据;对所述原始心率变异率数据进行线性拟合,得到所述原始心率变异率数据的变化趋势数据;确定所述原始心率变异率数据与所述变化趋势数据的差值,将所述差值作为确定出的所述心电信号的心率变异率。
CN201310722893.6A 2013-12-24 2013-12-24 一种睡眠阶段确定方法和系统 Active CN104720748B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310722893.6A CN104720748B (zh) 2013-12-24 2013-12-24 一种睡眠阶段确定方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310722893.6A CN104720748B (zh) 2013-12-24 2013-12-24 一种睡眠阶段确定方法和系统

Publications (2)

Publication Number Publication Date
CN104720748A true CN104720748A (zh) 2015-06-24
CN104720748B CN104720748B (zh) 2017-06-06

Family

ID=53445590

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310722893.6A Active CN104720748B (zh) 2013-12-24 2013-12-24 一种睡眠阶段确定方法和系统

Country Status (1)

Country Link
CN (1) CN104720748B (zh)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105262910A (zh) * 2015-11-02 2016-01-20 上海斐讯数据通信技术有限公司 通过检测脉搏来控制手机关机的方法和系统
CN105640508A (zh) * 2016-03-30 2016-06-08 安徽华米信息科技有限公司 实时睡眠监测方法及装置、智能可穿戴设备
CN106333652A (zh) * 2016-10-18 2017-01-18 首都医科大学 一种睡眠状态分析方法
CN106388771A (zh) * 2016-08-16 2017-02-15 惠州市德赛工业研究院有限公司 一种自动检测人体生理状态的方法及运动手环
CN107715273A (zh) * 2017-10-12 2018-02-23 西南大学 闹钟实现方法、装置及系统
CN107890339A (zh) * 2017-11-09 2018-04-10 常熟理工学院 一种睡眠分期检测方法及可穿戴式睡眠分期检测装置
CN107961429A (zh) * 2017-11-28 2018-04-27 广州视源电子科技股份有限公司 睡眠辅助方法和系统、睡眠辅助装置
CN107998499A (zh) * 2017-11-28 2018-05-08 广州视源电子科技股份有限公司 睡眠辅助内容的处理方法和系统、睡眠辅助服务器系统
CN107998500A (zh) * 2017-11-28 2018-05-08 广州视源电子科技股份有限公司 睡眠辅助内容的播放方法和系统、睡眠辅助装置
CN108201435A (zh) * 2017-12-06 2018-06-26 深圳和而泰数据资源与云技术有限公司 睡眠分期确定方法、相关设备及计算机可读介质
CN108937867A (zh) * 2018-06-29 2018-12-07 出门问问信息科技有限公司 一种睡眠状态的监测方法及装置
CN108968915A (zh) * 2018-06-12 2018-12-11 山东大学 基于熵特征及支持向量机的睡眠状态分类方法及系统
CN109363669A (zh) * 2018-10-30 2019-02-22 深圳和而泰数据资源与云技术有限公司 眼罩和计算机可读存储介质
CN109745037A (zh) * 2019-02-14 2019-05-14 广州视源电子科技股份有限公司 一种图像处理方法、装置、设备和存储介质
US10376670B2 (en) 2013-07-08 2019-08-13 Resmed Sensor Technologies Limited Methods and systems for sleep management
CN111743513A (zh) * 2020-06-29 2020-10-09 武汉大学 基于混合通道eeg信号的s1与rem睡眠状态检测方法及系统
CN112168139A (zh) * 2019-07-05 2021-01-05 腾讯科技(深圳)有限公司 一种健康监控方法、装置及存储介质
US11648373B2 (en) 2013-07-08 2023-05-16 Resmed Sensor Technologies Limited Methods and systems for sleep management
CN116269258A (zh) * 2023-02-22 2023-06-23 昆山市第一人民医院 一种枕式睡眠检测方法、装置及计算机设备

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI698223B (zh) * 2017-11-20 2020-07-11 金寶電子工業股份有限公司 可辨識打瞌睡階段的穿戴式裝置及其辨識方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5902250A (en) * 1997-03-31 1999-05-11 President And Fellows Of Harvard College Home-based system and method for monitoring sleep state and assessing cardiorespiratory risk
CN102319063A (zh) * 2011-07-26 2012-01-18 深圳市理邦精密仪器股份有限公司 一种提高胎心率数据加速识别准确性的装置和方法
WO2012156427A1 (en) * 2011-05-18 2012-11-22 Alain Gilles Muzet System and method for determining sleep and sleep stages of a person

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5902250A (en) * 1997-03-31 1999-05-11 President And Fellows Of Harvard College Home-based system and method for monitoring sleep state and assessing cardiorespiratory risk
WO2012156427A1 (en) * 2011-05-18 2012-11-22 Alain Gilles Muzet System and method for determining sleep and sleep stages of a person
CN102319063A (zh) * 2011-07-26 2012-01-18 深圳市理邦精密仪器股份有限公司 一种提高胎心率数据加速识别准确性的装置和方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
肖蒙等: "基于小波分解和过零点分析的睡眠期间心率变异性研究", 《航天医学与医学工程》 *

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11364362B2 (en) 2013-07-08 2022-06-21 Resmed Sensor Technologies Limited Methods and systems for sleep management
US10376670B2 (en) 2013-07-08 2019-08-13 Resmed Sensor Technologies Limited Methods and systems for sleep management
US11648373B2 (en) 2013-07-08 2023-05-16 Resmed Sensor Technologies Limited Methods and systems for sleep management
CN105262910A (zh) * 2015-11-02 2016-01-20 上海斐讯数据通信技术有限公司 通过检测脉搏来控制手机关机的方法和系统
CN105640508A (zh) * 2016-03-30 2016-06-08 安徽华米信息科技有限公司 实时睡眠监测方法及装置、智能可穿戴设备
CN106388771A (zh) * 2016-08-16 2017-02-15 惠州市德赛工业研究院有限公司 一种自动检测人体生理状态的方法及运动手环
CN106388771B (zh) * 2016-08-16 2019-06-28 惠州市德赛工业研究院有限公司 一种自动检测人体生理状态的方法及运动手环
CN106333652A (zh) * 2016-10-18 2017-01-18 首都医科大学 一种睡眠状态分析方法
CN106333652B (zh) * 2016-10-18 2019-10-29 首都医科大学 一种睡眠状态分析方法
CN107715273A (zh) * 2017-10-12 2018-02-23 西南大学 闹钟实现方法、装置及系统
CN107890339A (zh) * 2017-11-09 2018-04-10 常熟理工学院 一种睡眠分期检测方法及可穿戴式睡眠分期检测装置
CN107890339B (zh) * 2017-11-09 2020-09-08 常熟理工学院 一种睡眠分期检测方法及可穿戴式睡眠分期检测装置
CN107998500A (zh) * 2017-11-28 2018-05-08 广州视源电子科技股份有限公司 睡眠辅助内容的播放方法和系统、睡眠辅助装置
CN107998499A (zh) * 2017-11-28 2018-05-08 广州视源电子科技股份有限公司 睡眠辅助内容的处理方法和系统、睡眠辅助服务器系统
CN107961429A (zh) * 2017-11-28 2018-04-27 广州视源电子科技股份有限公司 睡眠辅助方法和系统、睡眠辅助装置
CN108201435A (zh) * 2017-12-06 2018-06-26 深圳和而泰数据资源与云技术有限公司 睡眠分期确定方法、相关设备及计算机可读介质
CN108968915A (zh) * 2018-06-12 2018-12-11 山东大学 基于熵特征及支持向量机的睡眠状态分类方法及系统
CN108937867A (zh) * 2018-06-29 2018-12-07 出门问问信息科技有限公司 一种睡眠状态的监测方法及装置
CN109363669A (zh) * 2018-10-30 2019-02-22 深圳和而泰数据资源与云技术有限公司 眼罩和计算机可读存储介质
CN109745037A (zh) * 2019-02-14 2019-05-14 广州视源电子科技股份有限公司 一种图像处理方法、装置、设备和存储介质
CN112168139A (zh) * 2019-07-05 2021-01-05 腾讯科技(深圳)有限公司 一种健康监控方法、装置及存储介质
CN111743513B (zh) * 2020-06-29 2021-05-18 武汉大学 基于混合通道eeg信号的s1与rem睡眠状态检测方法及系统
CN111743513A (zh) * 2020-06-29 2020-10-09 武汉大学 基于混合通道eeg信号的s1与rem睡眠状态检测方法及系统
CN116269258A (zh) * 2023-02-22 2023-06-23 昆山市第一人民医院 一种枕式睡眠检测方法、装置及计算机设备
CN116269258B (zh) * 2023-02-22 2023-11-21 昆山市第一人民医院 一种枕式睡眠检测方法、装置及计算机设备

Also Published As

Publication number Publication date
CN104720748B (zh) 2017-06-06

Similar Documents

Publication Publication Date Title
CN104720748B (zh) 一种睡眠阶段确定方法和系统
CN107595245B (zh) 一种睡眠管理方法、系统及终端设备
CN104720746B (zh) 一种睡眠阶段确定方法和系统
CN109316170B (zh) 基于深度学习的脑电波辅助睡眠及唤醒系统
CN105395192A (zh) 一种基于脑电的可穿戴情感识别方法和系统
CN105827731A (zh) 基于融合模型的智能化健康管理服务器、系统及其控制方法
CN107106085A (zh) 用于睡眠监测的设备和方法
CN110604565A (zh) 一种基于便携式eeg设备的脑部健康筛查方法
WO2020187109A1 (zh) 一种用户睡眠检测方法及系统
CN206045144U (zh) 一种新型智能助眠和自然唤醒的装置
Tian et al. A hierarchical classification method for automatic sleep scoring using multiscale entropy features and proportion information of sleep architecture
CN105496356A (zh) 睡眠状态判定装置、睡眠状态判定方法及睡眠管理系统
CN109833031A (zh) 一种基于lstm利用多生理信号的自动睡眠分期方法
CN110974219A (zh) 一种基于侵入式bci的人脑意念识别系统
CN104615851A (zh) 一种睡眠监控方法及终端
CN110037693A (zh) 一种基于面部表情和eeg的情绪分类方法
CN113303814B (zh) 基于深度迁移学习的单通道耳脑电自动睡眠分期方法
CN107890339A (zh) 一种睡眠分期检测方法及可穿戴式睡眠分期检测装置
CN110706816A (zh) 一种基于人工智能进行睡眠环境调控的方法及设备
CN106730232B (zh) 一种智能唤醒方法及系统
CN110584657B (zh) 一种注意力检测方法及系统
CN109044280A (zh) 一种睡眠分期方法及相关设备
CN109620158A (zh) 睡眠监控的方法、智能终端及存储装置
CN109009098A (zh) 一种运动想象状态下的脑电信号特征识别方法
CN108491792B (zh) 基于眼电信号的办公场景人机交互行为识别方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant