CN112617761A - 自适应聚点生成的睡眠阶段分期方法 - Google Patents

自适应聚点生成的睡眠阶段分期方法 Download PDF

Info

Publication number
CN112617761A
CN112617761A CN202011642364.1A CN202011642364A CN112617761A CN 112617761 A CN112617761 A CN 112617761A CN 202011642364 A CN202011642364 A CN 202011642364A CN 112617761 A CN112617761 A CN 112617761A
Authority
CN
China
Prior art keywords
sleep
data
stage
staging
state
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
CN202011642364.1A
Other languages
English (en)
Other versions
CN112617761B (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.)
Hunan Zennze Technology Co ltd
Original Assignee
Hunan Dongsheng Nanxiang Intelligent Technology 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 Hunan Dongsheng Nanxiang Intelligent Technology Co ltd filed Critical Hunan Dongsheng Nanxiang Intelligent Technology Co ltd
Priority to CN202011642364.1A priority Critical patent/CN112617761B/zh
Publication of CN112617761A publication Critical patent/CN112617761A/zh
Application granted granted Critical
Publication of CN112617761B publication Critical patent/CN112617761B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • 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
    • 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
    • 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/4815Sleep quality
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Medical Informatics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Computation (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

本发明提供了自适应聚点生成的睡眠阶段分期方法包括如下步骤:包括如下步骤:步骤S1、获取睡眠原始监测数据;步骤S2、数据预处理:去除状态异常数据,并对数据进行统计学分位数处理;步骤S3、特征数据分析,进行特征选取提取;步骤S4、自适应聚类算法寻找睡眠特征的聚点:通过加权聚类算法,提取k个特征属性的聚点Truthk,并获得聚点距离范围内的误差指标值αk;步骤S5、构建睡眠分期识别模型:依据聚点Truthk和误差指标αk规则识别睡眠状态及睡眠分期阶段;步骤S6、睡眠分期修正:通过初步判断的睡眠分期阶段进行归一化加权评分分析,精准更新睡眠分期阶段。本发明的睡眠分类模型计算运行时间少,算法的准确度高。

Description

自适应聚点生成的睡眠阶段分期方法
技术领域
本发明涉及大数据应用及人工智能技术领域,特别涉及自适应聚点生成的睡眠阶段分期方法。
背景技术
睡眠质量与人类身体的健康息息相关,如今大多数人都存在睡眠问题的困扰,有些是短暂性的,有些是长期性的,整夜低质量的睡眠会对人们精力、心情、健康等各方面产生很大的负面影响。人们开始逐渐关注睡眠质量的问题。长期监控人们整夜的睡眠情况,去知晓睡眠阶段分布和睡眠质量,探索人们长期睡眠质量不好的主要原因,并采取及时有效的改善措施,就显得更加必要了。因此,人们发明便捷式、可穿戴式、无接触式等设备来监测整夜的睡眠状况。目前,在国际通用的R&K(Rechtschaffen&Kales)准则中,睡眠过程主要分为清醒期(aWake,W)、快速眼动期(Rapid EyeMovement,REM)和非快速眼动期(Non-RapidEye Movement,NREM)。在非快速眼动期中,随着人体睡眠状态的由浅入深,又将非快速眼动期分为四个阶段,分别为睡眠Ⅰ期(Stage1,S1),睡眠Ⅱ期(Stage2,S2),睡眠III期(Stage3,S3)以及睡眠Ⅳ期(Stage4,S4)。美国睡眠医学学会提出S3和S4期都处在深度睡眠状态,将S3和S4期合并为慢波睡眠期(Slow wave Sleep,SS)。睡眠分期对于睡眠质量评估具有重要价值,只要能基于设备监测的睡眠数据进行分期,就可以实现家庭和日常生活状态下的睡眠质量监测和评估。
1932年以来相继引入傅里叶变换、频域分析等对睡眠脑电图进行了研究,因为脑电本身具有随机性和非平稳性质,近几年非线性系统分析理论也被迅速运用到了睡眠分期中。一些研究人员已经采用复杂度、熵值、奇异谱分析、独立成分分析等方法来区分不同的睡眠时相。睡眠脑电的自动分期经过分析在不同睡眠时期下脑电的变化趋势来划分人的睡眠结构,在这个过程中,提取信号的特征和模式识别是实现准确分期的关键环节。通过提取信号的特征和模式识别不同睡眠分期,大多数研究方法主要集中在监督学习算法,例如常见的神经网络、隐马尔可夫模型等用于对睡眠状态的脑电信号进行建模和分析以及支持向量机也经常被用来实现睡眠阶段的自动分类。其中,分析睡眠阶段的心率的变异性,从心率时间间期序列中提取时域、频域和非线性等86个特征,将多导睡眠图仪的三分类结果(醒、快速眼动期、非快速眼动期)作为“金标准”,并采用支持向量机作为多分类器模型,并实验论证了具有较高的准确性,大样本人群测试结果表明该模型具有较好的普适性。然而,脑电等生物电信号存在复杂差异性,一方面表现在不同个体,另一方面也受到记录环境和记录条件的影响。然监督学习算法的分类性能主要取决于训练过程以及训练样本,对特定范围内的数据集会有较好的分类效果,难以适应实际应用中各种不同睡眠数据。而无监督方法不需要对样本进行训练,是根据样本的原始信息来完成分类的。运用无监督方法来解决不同睡眠阶段的分类判别问题是值得研究的方向,也具备实际应用价值。
现有分类模型召回率、精确率两个指标均较低,分类效果有限。另外,算法的运行时间长,算法效率低。
发明内容
本发明提供了自适应聚点生成的睡眠阶段分期方法,解决现有技术中睡眠分类模型计算运行时间多和算法的准确度低的技术问题。
为了达到上述目的,本发明提供的自适应聚点生成的睡眠阶段分期方法,包括如下步骤:
步骤S1、获取睡眠原始监测数据:通过生命体征参数测量设备,获取生命体征参数;
步骤S2、数据预处理:去除状态异常数据,并对数据进行统计学分位数处理;
步骤S3、特征数据分析,进行特征选取提取:提取睡眠各阶段各指标样本数据,通过分析各阶段的差异,选择能表征睡眠阶段的特征属性;
步骤S4、自适应聚类算法寻找睡眠特征的聚点:通过加权聚类算法,提取特征属性的聚点Truthk,并获得聚点距离范围内的误差指标值αk
步骤S5、构建睡眠分期识别模型:依据聚点Truthk和误差指标αk规则识别睡眠状态及睡眠分期阶段;
步骤S6、睡眠分期修正:通过初步判断的睡眠分期阶段进行归一化加权评分分析,精准更新睡眠分期阶段。
优选地,所述步骤S1中所述生命体征参数测量设备为非接触式生命体征监护雷达设备,所述生命体征参数包括设备与人体之间的检测距离、呼吸率、心率、设备返回信号强度大小、人体异常状态、体动状态和设备所检测数据时间。
优选地,所述步骤S2包括如下步骤:
步骤S21、去除状态异常数据:去除状态识别为“小动作”、“持续动作”、“平静”、“安静”、“呼吸骤停”或“跌倒”的情况,且在状态开始至结束的一分钟内所产生的数据将作为噪声不考虑计算;
步骤S22、对去除状态异常数据之后的数据进行统计学分位数处理:使用统计学中分位数处理方法,将所有数值由小到大排列并分成十等份,取中间第二等份至第九等份数据作为重要参考数据。
优选地,所述步骤S3具体为:根据睡眠各阶段不同特征样本数据,绘制各阶段不同特征曲线走势图,通过分析各阶段的差异选择能表征睡眠阶段的特征属性,所述能表征睡眠阶段的特征属性包括设备与人体之间的检测距离、呼吸率、心率、设备返回信号强度大小。
优选地,所述步骤S4具体包括如下步骤:
步骤S41、确定聚点真值估计值公式:
聚点真值估计值表示为:
Figure BDA0002880727080000031
其中xi表示第i数据源(设备每秒上传的数据)所提供的观测值,ws表示数据源的权重,s表示数据源;
步骤S42、确定数据源权值分配模型:
步骤S43、确定加权聚类算法流程;具体为:根据数据源提供的观测值,初始化聚点真值,并通过损失函数及对数标准化函数引入权重,更新聚点真值;
步骤S44、确定聚点范围内的误差指标值αk
Figure BDA0002880727080000041
式中,
Figure BDA0002880727080000042
表示第k个特征属性的第i个数据源所提供的观测值,i表示聚合x的个数,k表示特征属性,Std(·)表示标准差。
优选地,所述步骤S42具体为:每个睡眠时间点的每个属性值x的权重w:
score=|x-max(x)|/std(x)
Figure BDA0002880727080000043
其中θ=0.00001
其中,std(x)为属性值x的标准差,max(x)为属性值x的最大值,score为获取权重的中间参数,∑score为所有属性值x的中间参数之和,max(score)为所有属性值x的中间参数中的最大值;
所述步骤S43具体为:使用加权平均函数WeightedMedian,计算出单个特征属性下,每个属性值x的权重w的聚类结果聚点值Truth:
Truth=WeightedMedian(x,w)
每个属性值x数据对应一个权重,从属性值x数据开始点累加,直到到达数据总和的1/2,则该属性值x数据点选为这个训练数据集合的聚点值Truth。
优选地,所述步骤S5包括如下步骤:
步骤S51、睡眠状态识别:以1秒时间周期为窗口,依据Truthk和αk规则判断时间戳ti是“睡眠”或“非睡眠”状态:
Figure BDA0002880727080000044
式中,
Figure BDA0002880727080000045
表示时间戳在ti时刻的属性值;
步骤S52、睡眠分期阶段识别:
当状态为“睡眠”,且T时间周期内出现次数大于等于τ1时,判断睡眠分期阶段为“深睡眠”;
当状态为“睡眠”,且T时间周期内出现次数小于τ1,且大于等于τ2时,判断睡眠分期阶段为“浅睡眠”;
当状态为“睡眠”,且T时间周期内出现次数小于τ2,且大于等于τ3时,判断睡眠分期阶段为“快速眼动”;
当状态为“睡眠”,且T时间周期内出现次数小于τ3,且大于等于τ4时,判断睡眠分期阶段为“清醒”;
其中,τ1,τ2,τ3,τ4均为历史训练经验参数。
优选地,所述步骤S6具体为:提取每一阶段睡眠状态的各特征属性数据
Figure BDA0002880727080000051
再使用评分函数,计算该睡眠分期阶段下数据睡眠数据评分,根据数据得分重新更新睡眠状态;具体包括如下步骤:
步骤S61、集合
Figure BDA0002880727080000052
睡眠状态处理得出各状态数据集合
Figure BDA0002880727080000053
Figure BDA0002880727080000054
步骤S62、对集合Yk中四个特征属性值进行归一化处理(k=[“距离”、“呼吸率”、“心率”、“信号强度”])
Figure BDA0002880727080000055
步骤S63、对归一化后的四个属性的标准差求和
Figure BDA0002880727080000056
步骤S61、根据score值调整阈值,对“深睡眠”、“浅睡眠”、“快速眼动”状态进行更新:S1<S2<S3
当score2小于等于S1时,更新睡眠分期阶段为“深睡眠”;
当score2大于等于S1,且小于S2时,更新睡眠分期阶段为“浅睡眠”;
当score2大于等于S2,且小于S3时,更新睡眠分期阶段为“快速眼动”;
其中,S1,S2,S3均为历史训练经验参数。
优选地,还包括如下步骤:
步骤S7、形成完整的睡眠睡觉时间序列线:将得出深睡眠、浅睡眠、快速眼动、清醒的时间周期按时间序列拼接,形成一条完整的睡眠睡觉时间序列线。
采用本发明的自适应聚点生成的睡眠阶段分期方法,通过多数据源可信度的真值发现方法,凭借该方法的快速收敛聚类特性,识别睡眠状态的特征差异,找出各特征的聚点中心。并基于睡眠数据聚点进行睡眠分期阶段识别,将整晚睡眠数据分为清醒、快速眼动、浅睡眠、深睡眠期4个睡眠阶段。本发明的算法具有较高的准确性,大大节省了算法的运行消耗时间,具有高效性。
附图说明
图1为本发明的自适应聚点生成的睡眠阶段分期方法的流程示意图;
图2为本发明的自适应聚点生成的睡眠阶段分期方法的一较佳实施例的睡眠分期阶段识别的框架图;
图3为本发明的自适应聚点生成的睡眠阶段分期方法的效果验证过程中一名受试者的整晚睡眠采用人工专家判读标准的时相序列图;
图4为本发明的自适应聚点生成的睡眠阶段分期方法的效果验证过程中采用本发明的自适应聚点生成的睡眠阶段分期方法算法划分的睡眠分期图。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述。
本发明针对现有的问题,提供了自适应聚点生成的睡眠阶段分期方法,本发明的主要算法流程如下:首先提取历史记录中已知的整夜完整睡眠数据进行训练得到模型所需的聚点,然后通过真实已知睡眠阶段(睡眠开始、结束时间,深睡眠、浅睡眠、快速眼动、清醒状态时间段等睡眠报告指标)的数据对模型反复验证,从而得到最合理和最精准的聚点参数,最后使用所训练好的聚点和模型对设备上传的24小时睡眠监测数据进行处理,生成一条完整的睡眠阶段分期时间序列。
如图1所示,本发明的自适应聚点生成的睡眠阶段分期方法包括如下步骤:
步骤S1、获取睡眠原始监测数据:通过生命体征参数测量设备,获取生命体征参数;
步骤S2、数据预处理:去除状态异常数据,并对数据进行统计学分位数处理;
步骤S3、特征数据分析,进行特征选取提取:提取睡眠各阶段各指标样本数据,通过分析各阶段的差异,选择能表征睡眠阶段的特征属性;
步骤S4、自适应聚类算法寻找睡眠特征的聚点:通过加权聚类算法,提取k个特征属性的聚点Truthk,并获得聚点距离范围内的误差指标值αk
步骤S5、构建睡眠分期识别模型:依据聚点Truthk和误差指标αk规则识别睡眠状态及睡眠分期阶段;
步骤S6、睡眠分期修正:通过初步判断的睡眠分期阶段进行归一化加权评分分析,精准更新睡眠分期阶段;
步骤S7、形成完整的睡眠睡觉时间序列线:将得出深睡眠、浅睡眠、快速眼动、清醒的时间周期按时间序列拼接,形成一条完整的睡眠睡觉时间序列线。
所述步骤S1中所述生命体征参数测量设备为非接触式生命体征监护雷达设备,所述生命体征参数包括设备与人体之间的检测距离、呼吸率、心率、设备返回信号强度大小、人体异常状态、体动状态和设备所检测数据时间。
所述步骤S2包括如下步骤:
步骤S21、去除状态异常数据:去除状态识别为“小动作”、“持续动作”、“平静”、“安静”、“呼吸骤停”或“跌倒”的情况,且在状态开始至结束的一分钟内所产生的数据将作为噪声不考虑计算;
步骤S22、对去除状态异常数据之后的数据进行统计学分位数处理:使用统计学中分位数处理方法,将所有数值由小到大排列并分成十等份,取中间第二等份至第九等份数据作为重要参考数据。
所述步骤S3具体为:根据睡眠各阶段不同特征样本数据,绘制各阶段不同特征曲线走势图,通过分析各阶段的差异选择能表征睡眠阶段的特征属性,所述能表征睡眠阶段的特征属性包括设备与人体之间的检测距离、呼吸率、心率、设备返回信号强度大小。
所述步骤S4具体包括如下步骤:
步骤S41、确定聚点真值估计值公式:
聚点真值估计值表示为:
Figure BDA0002880727080000071
其中xi表示第i数据源(设备每秒上传的数据)所提供的观测值,ws表示数据源的权重,s表示数据源;
步骤S42、确定数据源权值分配模型:具体为:每个睡眠时间点的每个属性值x的权重w:
score=|x-max(x)|/std(x)
Figure BDA0002880727080000081
其中θ=0.00001
其中,std(x)为属性值x的标准差,max(x)为属性值x的最大值,score为获取权重的中间参数,∑score为所有属性值x的中间参数之和,max(score)为所有属性值x的中间参数中的最大值。
步骤S43、确定加权聚类算法流程;具体为:根据数据源提供的观测值,初始化聚点真值,并通过损失函数及对数标准化函数引入权重,更新聚点真值具体为:使用加权平均函数WeightedMedian,计算出单个特征属性下,每个属性值x的权重w的聚类结果聚点值Truth:
Truth=WeightedMedian(x,w)
每个属性值x数据对应一个权重,从属性值x数据开始点累加,直到到达数据总和的1/2,则该属性值x数据点选为这个训练数据集合的聚点值Truth。
第k个特征属性下,属性值x的权重w的聚类结果聚点值即为Truthk
步骤S44、确定聚点范围内的误差指标值αk
Figure BDA0002880727080000082
式中,
Figure BDA0002880727080000083
表示第k个特征属性的第i个数据源所提供的观测值,i表示聚合x的个数,k表示特征属性,Std(·)表示标准差。
如图2所示,所述步骤S5包括如下步骤:
步骤S51、睡眠状态识别:以1秒时间周期为窗口,依据Truthk和αk规则判断时间戳ti是“睡眠”或“非睡眠”状态:
Figure BDA0002880727080000084
式中,
Figure BDA0002880727080000085
表示时间戳在ti时刻的属性值;
步骤S52、睡眠分期阶段识别:
当状态为“睡眠”,且T时间周期内出现次数大于等于τ1时,判断睡眠分期阶段为“深睡眠”;
当状态为“睡眠”,且T时间周期内出现次数小于τ1,且大于等于τ2时,判断睡眠分期阶段为“浅睡眠”;
当状态为“睡眠”,且T时间周期内出现次数小于τ2,且大于等于τ3时,判断睡眠分期阶段为“快速眼动”;
当状态为“睡眠”,且T时间周期内出现次数小于τ3,且大于等于τ4时,判断睡眠分期阶段为“清醒”;
其中,τ1,τ2,τ3,τ4均为历史训练经验参数。
所述步骤S6具体为:提取每一阶段睡眠状态的各特征属性数据
Figure BDA0002880727080000091
Figure BDA0002880727080000092
再使用评分函数,计算该睡眠分期阶段下数据睡眠数据评分,根据数据得分重新更新睡眠状态;具体包括如下步骤:
步骤S61、集合
Figure BDA0002880727080000093
睡眠状态处理得出各状态数据集合
Figure BDA0002880727080000094
Figure BDA0002880727080000095
步骤S62、对集合Yk中四个特征属性值进行归一化处理(k=[“距离”、“呼吸率”、“心率”、“信号强度”])
Figure BDA0002880727080000096
步骤S63、对归一化后的四个属性的标准差求和
Figure BDA0002880727080000097
步骤S61、根据score值调整阈值,对“深睡眠”、“浅睡眠”、“快速眼动”状态进行更新:S1<S2<S3
当score2小于等于S1时,更新睡眠分期阶段为“深睡眠”;
当score2大于等于S1,且小于S2时,更新睡眠分期阶段为“浅睡眠”;
当score2大于等于S2,且小于S3时,更新睡眠分期阶段为“快速眼动”;
其中,S1,S2,S3均为历史训练经验参数。
采用本发明的自适应聚点生成的睡眠阶段分期方法,通过多数据源可信度的真值发现方法,凭借该方法的快速收敛聚类特性,识别睡眠状态的特征差异,找出各特征的聚点中心。并基于睡眠数据聚点进行睡眠分期阶段识别,将整晚睡眠数据分为清醒、快速眼动、浅睡眠、深睡眠期4个睡眠阶段。本发明的算法具有较高的准确性,大大节省了算法的运行消耗时间,具有高效性。
本发明将整晚睡眠数据划分为清醒(W)、快速眼动(REM)、浅睡眠(LS)、深睡眠(DS)期4个睡眠阶段。采用专家人工判读具体睡眠状态数据段作为训练数据集和校验标准,将整夜睡眠过程中的脑电信号分成每60s一段,分别提取其中的特征参数进行睡眠阶段的本发明提出的方法及对比实验方法分析,并将分类结果和专家的人工判读进行了对比。
表1被试及待训练睡眠阶段统计表
Figure BDA0002880727080000101
综合5名受试者的睡眠数据,不同数据集的通过无线传感设备所收集的所处的环境产生不同的差异,这些差异性也会影响并体现在睡眠数据中。本发明选取两类数据集进行实验,一类选取一位健康受试者的睡眠数据作为训练数据集,另一类则选取其他5名受试者的睡眠数据作为测试数据集,其各阶段睡眠数据统计结果如表1。
本发明将睡眠过程分为清醒(W)、快速眼动(REM)、浅睡眠(LS)、深睡眠(DS)期4个睡眠阶段,对4个类别分别计算了召回率、精确率两个指标来评价分类效果。另外,本发明不仅采用两个指标评价算法的准确度,还引入算法的运行时间评估效率。显然,运行时间越少和算法的准确度越高,算法的效率越高。召回率是评价在实际分类中某类样本被正确识别出来的概率;精确率是评价在某类预测的样本中能够被正确识别出来的样本概率。
图3、4是一名受试者的整晚睡眠时相序列图,其中图3是人工专家判读标准,图4是本发明提出聚类算法的结果。图4贴近图3人工判读的结果。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (9)

1.自适应聚点生成的睡眠阶段分期方法,其特征在于,包括如下步骤:
步骤S1、获取睡眠原始监测数据:通过生命体征参数测量设备,获取生命体征参数;
步骤S2、数据预处理:去除状态异常数据,并对数据进行统计学分位数处理;
步骤S3、特征数据分析,进行特征选取提取:提取睡眠各阶段各指标样本数据,通过分析各阶段的差异,选择能表征睡眠阶段的特征属性;
步骤S4、自适应聚类算法寻找睡眠特征的聚点:通过加权聚类算法,提取k个特征属性的聚点Truthk,并获得聚点距离范围内的误差指标值αk
步骤S5、构建睡眠分期识别模型:依据聚点Truthk和误差指标αk规则识别睡眠状态及睡眠分期阶段;
步骤S6、睡眠分期修正:通过初步判断的睡眠分期阶段进行归一化加权评分分析,精准更新睡眠分期阶段。
2.根据权利要求1所述的自适应聚点生成的睡眠阶段分期方法,其特征在于,所述步骤S1中所述生命体征参数测量设备为非接触式生命体征监护雷达设备,所述生命体征参数包括设备与人体之间的检测距离、呼吸率、心率、设备返回信号强度大小、人体异常状态、体动状态和设备所检测数据时间。
3.根据权利要求1所述的自适应聚点生成的睡眠阶段分期方法,其特征在于,所述步骤S2包括如下步骤:
步骤S21、去除状态异常数据:去除状态识别为“小动作”、“持续动作”、“平静”、“安静”、“呼吸骤停”或“跌倒”的情况,且在状态开始至结束的一分钟内所产生的数据将作为噪声不考虑计算;
步骤S22、对去除状态异常数据之后的数据进行统计学分位数处理:使用统计学中分位数处理方法,将所有数值由小到大排列并分成十等份,取中间第二等份至第九等份数据作为重要参考数据。
4.根据权利要求2所述的自适应聚点生成的睡眠阶段分期方法,其特征在于,所述步骤S3具体为:根据睡眠各阶段不同特征样本数据,绘制各阶段不同特征曲线走势图,通过分析各阶段的差异选择能表征睡眠阶段的特征属性,所述能表征睡眠阶段的特征属性包括设备与人体之间的检测距离、呼吸率、心率、设备返回信号强度大小。
5.根据权利要求1所述的自适应聚点生成的睡眠阶段分期方法,其特征在于,所述步骤S4具体包括如下步骤:
步骤S41、确定聚点真值估计值公式:
聚点真值估计值表示为:
Figure FDA0002880727070000021
其中xi表示第i个数据源(设备每秒上传的数据)所提供的观测值,ws表示数据源的权重,s表示数据源;
步骤S42、确定数据源权值分配模型:
步骤S43、确定加权聚类算法流程;具体为:根据数据源提供的观测值,初始化聚点真值,并通过损失函数及对数标准化函数引入权重,更新聚点真值;
步骤S44、确定聚点范围内的误差指标值αk
Figure FDA0002880727070000022
式中,
Figure FDA0002880727070000023
表示第k个特征属性的第i个数据源所提供的观测值,i表示聚合x的个数,k表示特征属性,Std(·)表示标准差。
6.根据权利要求5所述的自适应聚点生成的睡眠阶段分期方法,其特征在于,所述步骤S42具体为:每个睡眠时间点的每个属性值x的权重w:
score=|x-max(x)|/std(x)
Figure FDA0002880727070000024
其中θ=0.00001
其中,std(x)为属性值x的标准差,max(x)为属性值x的最大值,score为获取权重的中间参数,∑score为所有属性值x的中间参数之和,max(score)为所有属性值x的中间参数中的最大值;
所述步骤S43具体为:使用加权平均函数WeightedMedian,计算出单个特征属性下,每个属性值x的权重w的聚类结果聚点值Truth:
Truth=WeightedMedian(x,w)
每个属性值x数据对应一个权重,从属性值x数据开始点累加,直到到达数据总和的1/2,则该属性值x数据点选为这个训练数据集合的聚点值Truth。
7.根据权利要求1所述的自适应聚点生成的睡眠阶段分期方法,其特征在于,所述步骤S5包括如下步骤:
步骤S51、睡眠状态识别:以1秒时间周期为窗口,依据Truthk和αk规则判断时间戳ti是“睡眠”或“非睡眠”状态:
Figure FDA0002880727070000031
式中,
Figure FDA0002880727070000032
表示时间戳在ti时刻的属性值;
步骤S52、睡眠分期阶段识别:
当状态为“睡眠”,且T时间周期内出现次数大于等于τ1时,判断睡眠分期阶段为“深睡眠”;
当状态为“睡眠”,且T时间周期内出现次数小于τ1,且大于等于τ2时,判断睡眠分期阶段为“浅睡眠”;
当状态为“睡眠”,且T时间周期内出现次数小于τ2,且大于等于τ3时,判断睡眠分期阶段为“快速眼动”;
当状态为“睡眠”,且T时间周期内出现次数小于τ3,且大于等于τ4时,判断睡眠分期阶段为“清醒”;
其中,τ1,τ2,τ3,τ4均为历史训练经验参数。
8.根据权利要求1所述的自适应聚点生成的睡眠阶段分期方法,其特征在于,所述步骤S6具体为:提取每一阶段睡眠状态的各特征属性数据
Figure FDA0002880727070000034
Figure FDA0002880727070000035
再使用评分函数,计算该睡眠分期阶段下数据睡眠数据评分,根据数据得分重新更新睡眠状态;具体包括如下步骤:
步骤S61、集合
Figure FDA0002880727070000036
睡眠状态处理得出各状态数据集合
Figure FDA0002880727070000037
Figure FDA0002880727070000038
步骤S62、对集合Yk中四个特征属性值进行归一化处理(k=[“距离”、“呼吸率”、“心率”、“信号强度”])
Figure FDA0002880727070000033
步骤S63、对归一化后的四个属性的标准差求和
Figure FDA0002880727070000041
步骤S61、根据score2值调整阈值,对“深睡眠”、“浅睡眠”、“快速眼动”状态进行更新:S1<S2<S3
当score2小于等于S1时,更新睡眠分期阶段为“深睡眠”;
当score2大于等于S1,且小于S2时,更新睡眠分期阶段为“浅睡眠”;
当score2大于等于S2,且小于S3时,更新睡眠分期阶段为“快速眼动”;
其中,S1,S2,S3均为历史训练经验参数。
9.根据权利要求1所述的自适应聚点生成的睡眠阶段分期方法,其特征在于,还包括如下步骤:
步骤S7、形成完整的睡眠睡觉时间序列线:将得出深睡眠、浅睡眠、快速眼动、清醒的时间周期按时间序列拼接,形成一条完整的睡眠睡觉时间序列线。
CN202011642364.1A 2020-12-31 2020-12-31 自适应聚点生成的睡眠阶段分期方法 Active CN112617761B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011642364.1A CN112617761B (zh) 2020-12-31 2020-12-31 自适应聚点生成的睡眠阶段分期方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011642364.1A CN112617761B (zh) 2020-12-31 2020-12-31 自适应聚点生成的睡眠阶段分期方法

Publications (2)

Publication Number Publication Date
CN112617761A true CN112617761A (zh) 2021-04-09
CN112617761B CN112617761B (zh) 2023-10-13

Family

ID=75290218

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011642364.1A Active CN112617761B (zh) 2020-12-31 2020-12-31 自适应聚点生成的睡眠阶段分期方法

Country Status (1)

Country Link
CN (1) CN112617761B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114159024A (zh) * 2021-11-17 2022-03-11 青岛海信日立空调系统有限公司 一种睡眠分期方法及装置

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101272732A (zh) * 2005-05-10 2008-09-24 索尔克生物学研究所 睡眠和清醒状态的自动检测
US20090131803A1 (en) * 2004-11-02 2009-05-21 University College Dublin-National University Of Ireland, Dublin A sleep monitoring system
CA2779265A1 (en) * 2008-11-14 2010-05-20 Philip Low Methods of identifying sleep and waking patterns and uses
JP2011083393A (ja) * 2009-10-14 2011-04-28 Osaka Bioscience Institute 睡眠ステージ自動判定の装置と方法およびそのためのコンピュータプログラム
US20110190599A1 (en) * 2010-02-02 2011-08-04 Nellcor Puritan Bennett Llc System And Method For Diagnosing Sleep Apnea Based On Results Of Multiple Approaches To Sleep Apnea Identification
CN102908130A (zh) * 2005-11-29 2013-02-06 风险获利有限公司 用于监测人体健康状态的设备
US20150245800A1 (en) * 2012-08-20 2015-09-03 Danmarks Tekniske Universitet Method for Detection Of An Abnormal Sleep Pattern In A Person
WO2016109807A1 (en) * 2015-01-02 2016-07-07 Hello, Inc. Room monitoring device and sleep analysis
WO2016110804A1 (en) * 2015-01-06 2016-07-14 David Burton Mobile wearable monitoring systems
CN106473703A (zh) * 2016-09-21 2017-03-08 广州视源电子科技股份有限公司 睡眠状态分类器的训练方法和系统
KR20180017854A (ko) * 2016-08-11 2018-02-21 고려대학교 산학협력단 착용형 센서를 이용한 일주기 생체리듬의 비정상성 탐지 장치 및 방법
CN109124572A (zh) * 2018-06-15 2019-01-04 四川斐讯信息技术有限公司 一种睡眠状态的判断方法、系统及空气净化器
US20190038216A1 (en) * 2016-02-03 2019-02-07 Nanyang Technological University Methods for detecting a sleep disorder and sleep disorder detection devices
KR102011126B1 (ko) * 2018-04-30 2019-08-14 고려대학교 산학협력단 뇌 연결성을 이용한 의식 수준 기반의 수면 단계 측정 방법 및 장치
CN110192847A (zh) * 2019-05-31 2019-09-03 湖南省顺鸿智能科技有限公司 一种生命体征监测方法及系统
US10425912B1 (en) * 2019-01-17 2019-09-24 Cisco Technology, Inc. Characterizing movement behaviors of wireless nodes in a network
CN111067503A (zh) * 2019-12-31 2020-04-28 深圳安视睿信息技术股份有限公司 一种基于心率变异性的睡眠分期方法
US20200138366A1 (en) * 2008-11-14 2020-05-07 Neurovigil, Inc. Methods of identifying sleep & waking patterns and uses

Patent Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090131803A1 (en) * 2004-11-02 2009-05-21 University College Dublin-National University Of Ireland, Dublin A sleep monitoring system
CN101272732A (zh) * 2005-05-10 2008-09-24 索尔克生物学研究所 睡眠和清醒状态的自动检测
CN102908130A (zh) * 2005-11-29 2013-02-06 风险获利有限公司 用于监测人体健康状态的设备
US20200138366A1 (en) * 2008-11-14 2020-05-07 Neurovigil, Inc. Methods of identifying sleep & waking patterns and uses
CA2779265A1 (en) * 2008-11-14 2010-05-20 Philip Low Methods of identifying sleep and waking patterns and uses
CN102438515A (zh) * 2008-11-14 2012-05-02 索尔克生物学研究所 识别睡眠模式与清醒模式的方法及用途
JP2011083393A (ja) * 2009-10-14 2011-04-28 Osaka Bioscience Institute 睡眠ステージ自動判定の装置と方法およびそのためのコンピュータプログラム
US20110190599A1 (en) * 2010-02-02 2011-08-04 Nellcor Puritan Bennett Llc System And Method For Diagnosing Sleep Apnea Based On Results Of Multiple Approaches To Sleep Apnea Identification
US20150245800A1 (en) * 2012-08-20 2015-09-03 Danmarks Tekniske Universitet Method for Detection Of An Abnormal Sleep Pattern In A Person
WO2016109807A1 (en) * 2015-01-02 2016-07-07 Hello, Inc. Room monitoring device and sleep analysis
WO2016110804A1 (en) * 2015-01-06 2016-07-14 David Burton Mobile wearable monitoring systems
US20190038216A1 (en) * 2016-02-03 2019-02-07 Nanyang Technological University Methods for detecting a sleep disorder and sleep disorder detection devices
KR20180017854A (ko) * 2016-08-11 2018-02-21 고려대학교 산학협력단 착용형 센서를 이용한 일주기 생체리듬의 비정상성 탐지 장치 및 방법
CN106473703A (zh) * 2016-09-21 2017-03-08 广州视源电子科技股份有限公司 睡眠状态分类器的训练方法和系统
KR102011126B1 (ko) * 2018-04-30 2019-08-14 고려대학교 산학협력단 뇌 연결성을 이용한 의식 수준 기반의 수면 단계 측정 방법 및 장치
CN109124572A (zh) * 2018-06-15 2019-01-04 四川斐讯信息技术有限公司 一种睡眠状态的判断方法、系统及空气净化器
US10425912B1 (en) * 2019-01-17 2019-09-24 Cisco Technology, Inc. Characterizing movement behaviors of wireless nodes in a network
CN110192847A (zh) * 2019-05-31 2019-09-03 湖南省顺鸿智能科技有限公司 一种生命体征监测方法及系统
CN111067503A (zh) * 2019-12-31 2020-04-28 深圳安视睿信息技术股份有限公司 一种基于心率变异性的睡眠分期方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
邵学晓: "基于聚类算法的睡眠自动分期方法研究", 基于聚类算法的睡眠自动分期方法研究, vol. 2019, no. 09, pages 3 - 40 *
郭祥迁: "基于生理信号的睡眠特征分析及睡眠分期方法研究", 《中国优秀博硕士学位论文全文数据库(硕士)医药卫生科技辑》, vol. 2018, no. 09, pages 5 - 90 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114159024A (zh) * 2021-11-17 2022-03-11 青岛海信日立空调系统有限公司 一种睡眠分期方法及装置
CN114159024B (zh) * 2021-11-17 2023-10-31 青岛海信日立空调系统有限公司 一种睡眠分期方法及装置

Also Published As

Publication number Publication date
CN112617761B (zh) 2023-10-13

Similar Documents

Publication Publication Date Title
Mendonca et al. A review of obstructive sleep apnea detection approaches
Wang et al. Detection analysis of epileptic EEG using a novel random forest model combined with grid search optimization
Rubin et al. Recognizing abnormal heart sounds using deep learning
Dash et al. Multi-channel EEG based automatic epileptic seizure detection using iterative filtering decomposition and Hidden Markov Model
US20170360362A1 (en) Sleep profiling system with feature generation and auto-mapping
CN111956221B (zh) 一种基于小波散射因子与lstm神经网络模型的颞叶癫痫分类方法
CN116861354A (zh) 电力变压器监测数据的异常识别及模式区分方法及系统
CN111407243A (zh) 一种基于深度学习的脉搏信号压力识别方法
CN107766898A (zh) 基于svm的三分类情绪概率判断方法
Furui et al. Non-Gaussianity detection of EEG signals based on a multivariate scale mixture model for diagnosis of epileptic seizures
CN109222966A (zh) 一种基于变分自编码器的脑电信号情感分类方法
CN112641451A (zh) 基于单通道脑电信号多尺度残差网络睡眠分期方法及系统
Na et al. An extended K nearest neighbors-based classifier for epilepsy diagnosis
CN112617761B (zh) 自适应聚点生成的睡眠阶段分期方法
Zhang et al. Nonlinear features of bark wavelet sub-band filtering for pathological voice recognition
US20220061694A1 (en) Lung health sensing through voice analysis
Raiesdana Automated sleep staging of OSAs based on ICA preprocessing and consolidation of temporal correlations
Zhang et al. Automatic sleep staging using multi-dimensional feature extraction and multi-kernel fuzzy support vector machine
CN116070168A (zh) 一种多模态数据多视角睡眠分期方法
Romaissa et al. Epileptic seizure detection from imbalanced EEG signal
CN115399735A (zh) 基于时频双流增强的多头注意机制睡眠分期方法
Javel et al. Epileptic seizure detection via EEG using tree-based pipeline optimization tool
Kaleem et al. Telephone-quality pathological speech classification using empirical mode decomposition
CN115270847A (zh) 基于小波包分解和卷积神经网络的设计决策脑电识别方法
Gini et al. An improved optimization algorithm for epileptic seizure detection in EEG signals using random forest classifier

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20230913

Address after: 501, Xianglu International Garden Phase II Hotel and Apartment, No. 61 Lufeng Road, High tech Development Zone, Changsha City, Hunan Province

Applicant after: HUNAN ZENNZE TECHNOLOGY CO.,LTD.

Address before: 618, building 1, Xiangyu wisdom, 579 Chezhan North Road, Dongfeng Road Street, Kaifu District, Changsha City, Hunan Province, 410000

Applicant before: Hunan Dongsheng Nanxiang Intelligent Technology Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant