CN116687356B - 基于时频分析的睡眠可持续性检测调节方法、系统和装置 - Google Patents
基于时频分析的睡眠可持续性检测调节方法、系统和装置 Download PDFInfo
- Publication number
- CN116687356B CN116687356B CN202310974382.7A CN202310974382A CN116687356B CN 116687356 B CN116687356 B CN 116687356B CN 202310974382 A CN202310974382 A CN 202310974382A CN 116687356 B CN116687356 B CN 116687356B
- Authority
- CN
- China
- Prior art keywords
- sleep
- sustainability
- curve
- time
- frequency
- 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.)
- Active
Links
- 230000007958 sleep Effects 0.000 title claims abstract description 854
- 238000004458 analytical method Methods 0.000 title claims abstract description 200
- 238000001514 detection method Methods 0.000 title claims abstract description 165
- 238000000034 method Methods 0.000 title claims abstract description 164
- 230000000694 effects Effects 0.000 claims abstract description 132
- 230000008569 process Effects 0.000 claims abstract description 82
- 230000004622 sleep time Effects 0.000 claims abstract description 51
- 238000012545 processing Methods 0.000 claims abstract description 36
- 238000011156 evaluation Methods 0.000 claims abstract description 27
- 230000000638 stimulation Effects 0.000 claims description 59
- 238000001914 filtration Methods 0.000 claims description 57
- 238000012937 correction Methods 0.000 claims description 45
- 230000035790 physiological processes and functions Effects 0.000 claims description 45
- 238000007726 management method Methods 0.000 claims description 34
- 238000000605 extraction Methods 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 24
- 230000009466 transformation Effects 0.000 claims description 24
- 238000005457 optimization Methods 0.000 claims description 16
- 238000011002 quantification Methods 0.000 claims description 13
- 230000000007 visual effect Effects 0.000 claims description 12
- 238000013139 quantization Methods 0.000 claims description 11
- 238000009432 framing Methods 0.000 claims description 10
- 238000007619 statistical method Methods 0.000 claims description 10
- 230000002093 peripheral effect Effects 0.000 claims description 9
- 230000011218 segmentation Effects 0.000 claims description 9
- 238000012544 monitoring process Methods 0.000 claims description 8
- 230000000241 respiratory effect Effects 0.000 claims description 8
- 230000001755 vocal effect Effects 0.000 claims description 7
- 206010062519 Poor quality sleep Diseases 0.000 claims description 6
- 230000001276 controlling effect Effects 0.000 claims description 6
- YHXISWVBGDMDLQ-UHFFFAOYSA-N moclobemide Chemical compound C1=CC(Cl)=CC=C1C(=O)NCCN1CCOCC1 YHXISWVBGDMDLQ-UHFFFAOYSA-N 0.000 claims description 6
- 238000012216 screening Methods 0.000 claims description 6
- 230000036541 health Effects 0.000 claims description 5
- 230000009467 reduction Effects 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 241001123248 Arma Species 0.000 claims description 4
- 238000010220 Pearson correlation analysis Methods 0.000 claims description 4
- 241000728173 Sarima Species 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 4
- 238000013135 deep learning Methods 0.000 claims description 4
- 230000006996 mental state Effects 0.000 claims description 4
- 238000002604 ultrasonography Methods 0.000 claims description 4
- 238000012952 Resampling Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 238000013500 data storage Methods 0.000 claims description 3
- 238000013079 data visualisation Methods 0.000 claims description 3
- 238000012217 deletion Methods 0.000 claims description 3
- 230000037430 deletion Effects 0.000 claims description 3
- 230000004927 fusion Effects 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 230000003862 health status Effects 0.000 claims description 3
- 238000010801 machine learning Methods 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 3
- 230000036385 rapid eye movement (rem) sleep Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 238000012549 training Methods 0.000 claims description 3
- 238000012512 characterization method Methods 0.000 claims 1
- 239000002131 composite material Substances 0.000 claims 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 3
- 241000282414 Homo sapiens Species 0.000 description 2
- 210000001015 abdomen Anatomy 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000003750 conditioning effect Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000004630 mental health Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 230000033764 rhythmic process Effects 0.000 description 2
- 230000004617 sleep duration Effects 0.000 description 2
- 230000035882 stress Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 206010052428 Wound Diseases 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004424 eye movement Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003340 mental effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000004461 rapid eye movement Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000004761 scalp Anatomy 0.000 description 1
- 230000008665 sleep physiology Effects 0.000 description 1
- 230000003860 sleep quality Effects 0.000 description 1
- 230000036578 sleeping time Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/01—Measuring temperature of body parts ; Diagnostic temperature sensing, e.g. for malignant or inflamed tissue
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/23—Updating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/10—Pre-processing; Data cleansing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/10—Pre-processing; Data cleansing
- G06F18/15—Statistical pre-processing, e.g. techniques for normalisation or restoring missing data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H15/00—ICT specially adapted for medical reports, e.g. generation or transmission thereof
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Engineering & Computer Science (AREA)
- Biophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Cardiology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Mathematical Physics (AREA)
- Physiology (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Computational Mathematics (AREA)
- Pulmonology (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Optics & Photonics (AREA)
- Probability & Statistics with Applications (AREA)
Abstract
本发明提供了基于时频分析的睡眠可持续性检测调节方法,包括:初始化睡眠可持续性动态检测调节基本方案,对用户睡眠生理信号进行采集和处理,得到睡眠状态特征曲线和睡眠时相曲线;根据所述睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值并进行时频分析,得到高频时频分析信号并进行去噪声处理,得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线,计算得到睡眠可持续性指数和睡眠可持续性指数曲线,进行趋势预测分析,生成睡眠可持续性动态调节策略并完成用户睡眠过程的动态调节及效果评估,生成用户睡眠可持续性检测调节报告。本发明能够实现对用户睡眠可持续性的高效干预调节。
Description
技术领域
本发明涉及睡眠可持续性检测量化及辅助调节领域,特别涉及基于时频分析的睡眠可持续性检测调节方法、系统和装置。
背景技术
睡眠是人类最重要的生命生理过程之一,稳定且连续性好的睡眠是人们获得精力体力、身心健康的重要保障。现实生活中,由于衰老、睡眠环境、生理疾病及创伤、精神压力等众多因素,整夜睡眠过程中出现多次觉醒、长时间觉醒和睡眠状态急剧震荡等一系列睡眠连续性或可持续性差的事件,极大地影响了人们的睡眠质量和身心体验。长期的低质量的睡眠连续性或可持续性,对人们的身心健康、学习、工作和生活带来重大不良影响。
申请人提出的在先解决方案中国申请CN2023101940795提供了一种睡眠可持续性检测量化及辅助干预的方法、系统和装置,包括以下步骤:对用户睡眠过程的生理状态信号和环境状态信号进行采集监测、信号处理和特征分析,生成生理状态特征和环境状态特征;对所述生理状态特征进行睡眠状态分析、时序成分分析和可持续性量化分析,提取睡眠可持续性指数,生成睡眠可持续性量化日报;重复上述步骤,对用户睡眠过程进行连续监测和跟踪分析,评估睡眠环境对睡眠可持续性的影响,提取最佳睡眠可持续性环境方案并对睡眠环境进行动态优化调整,生成睡眠可持续性量化报告。该技术方案基于睡眠持续期状态特征曲线的时间序列分解,提取了睡眠可持续性指数,实现提出了睡眠可持续性指数作为睡眠连续性或可持续性量化的创新评价指标。但面对睡眠可持续的高效提取、动态检测和动态调节等场景需求,还存在可以进一步的提升空间,主要包括:首先,时间序列分解方法在多场景可持续特征提取中存在适应性局限而带来准确性的不稳定;其次,睡眠持续期状态特征曲线是基于阶梯式睡眠时相分期值平滑得到的,不能够对用户睡眠状态进行细致刻画和精确量化,也进一步限制了睡眠可持续性评价的准确性和灵敏度;最后,除了睡眠环境调节干预外,如何实现用户更快速高效的、长期稳定的且具备高个性化的睡眠可持续性指数检测量化和动态调节,并能够动态持续地提升检测准确性和调节高效性。
如何提取更细致的用户睡眠状态连续变化特征,更准确地量化评估睡眠可持续性;如何实现睡眠可持续性检测量化和动态调节的高效协同,并能够实现检测调节的过程动态优化,持续提高检测量化效率和干预调节效果,是目前国内外产品技术方案和实际应用场景中需要进一步优化或解决的问题。
发明内容
针对现有方法的以上缺陷及改进需求,本发明的目的在于提供一种基于时频分析的睡眠可持续性检测调节方法,通过对睡眠状态特征曲线进行时频分析,得到高频时频分析信号并进行去噪声处理,提取睡眠可持续性时频特征曲线,计算得到睡眠可持续性指数和曲线;对用户睡眠状态进行趋势预测,生成睡眠可持续性动态调节策略并完成动态调节及效果评估,进而动态迭代优化睡眠可持续性的动态检测调节过程方法策略;对检测调节过程进行统计分析,生成用户睡眠可持续性检测调节报告,建立并持续更新用户个性睡眠可持续性数据库和睡眠可持续性通用数据库,最终实现提高检测量化效率和干预调节效果,辅助用户获得更高的睡眠质量和睡眠连续性。本发明还提供了基于时频分析的睡眠可持续性检测调节系统,用于实现上述方法。本发明还提供了基于时频分析的睡眠可持续性检测调节装置,用于实现上述系统。
根据本发明的目的,本发明提出了一种基于时频分析的睡眠可持续性检测调节方法,包括以下步骤:
初始化睡眠可持续性动态检测调节基本方案,对用户睡眠生理信号进行采集和处理,得到睡眠状态特征曲线和睡眠时相曲线;
根据所述睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值并进行时频分析,得到高频时频分析信号并进行去噪声处理,得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线;
根据所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线,计算得到睡眠可持续性指数和睡眠可持续性指数曲线;
对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析,生成睡眠可持续性动态调节策略并完成用户睡眠过程的动态调节及效果评估,动态迭代优化睡眠可持续性的动态检测调节过程方法策略;
生成用户睡眠可持续性检测调节报告;
初始化建立或持续更新用户个性睡眠可持续性数据库,更新睡眠可持续性通用数据库。
更优地,所述初始化睡眠可持续性动态检测调节基本方案,对用户睡眠生理信号进行采集和处理,得到睡眠状态特征曲线和睡眠时相曲线的具体步骤还包括:
获取用户关键生理信息,从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选并初始化所述睡眠可持续性动态检测调节基本方案;
对用户睡眠生理信号进行动态地监测采集和信号处理,得到睡眠生理状态数据;
对所述睡眠生理状态数据进行动态地特征分析和特征选择,生成所述睡眠状态特征曲线;
对所述睡眠生理状态数据进行动态地睡眠时相分析,生成所述睡眠时相曲线。
更优地,所述用户关键生理信息至少包括性别、年龄、生理健康状态和心理精神状态。
更优地,所述睡眠可持续性通用数据库具体为由规模的不同健康状态人群的睡眠可持续性的检测量化及动态调节数据库,至少包括所述用户关键生理信息、睡眠可持续性检测调节周期、睡眠可持续性检测量化过程方法参数、睡眠可持续性动态调节过程方法参数或策略、睡眠可持续性指数曲线、睡眠可持续性调节效果曲线和睡眠时相曲线。
更优地,所述用户个性睡眠可持续性数据库至少包括所述用户关键生理信息、睡眠状态特征曲线、睡眠时相曲线、睡眠可持续性指数曲线、动态调节效果曲线、时相指数分布特征、时相效果分布特征、动态调节综合效果系数、时频分析方法、趋势预测分析方法、睡眠可持续性动态调节策略、所述睡眠可持续性动态检测调节基本方案。
更优地,所述睡眠可持续性动态检测调节基本方案至少包括睡眠可持续性检测量化过程方案、睡眠可持续性动态调节过程方案;若用户为新用户,所述睡眠可持续性动态检测调节基本方案从所述睡眠可持续性通用数据库中获取,反之则从所述用户个性睡眠可持续性数据库中获取。
更优地,所述睡眠生理信号至少包括脑电信号、心电信号、呼吸信号中的任一项。
更优地,所述信号处理至少包括重采样、重参考、去伪迹、信号矫正、降噪、工频陷波、带通滤波、平滑处理和时帧分割;其中,所述时帧分割是指根据信号的采样率,以预设分帧步长对信号数据进行预设分帧时长窗口的连续滑动分割。
更优地,所述特征分析至少包括数值分析、包络分析、时频分析、熵分析、分形分析和复杂度分析中的任一项。
更优地,所述睡眠状态特征曲线具体由不同睡眠生理状态数据经所述特征分析得到不同睡眠生理状态特征进行线性组合融合而得到,用于描述用户在不同睡眠时期、不同睡眠时相、不同睡眠状态下的生理状态特征的连续变化曲线;其中,所述睡眠时期至少包括入睡前、睡眠中和睡眠后,所述睡眠时相至少包括清醒觉醒期、浅睡眠期、深睡眠期和快速眼动睡眠期。
更优地,所述睡眠时相曲线的生成方法具体为:
1)通过机器学习算法对规模睡眠用户样本的所述睡眠生理状态数据及其对应睡眠时相分期数据进行学习训练和数据建模,得到睡眠时相识别模型;
2)将当前用户的所述睡眠生理状态数据输入所述睡眠时相识别模型,得到所对应的睡眠时相分期并按时序生成所述睡眠时相曲线。
更优地,所述根据所述睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值并进行时频分析,得到高频时频分析信号并进行去噪声处理,得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线的具体步骤还包括:
根据所述睡眠状态特征曲线的曲线生成属性,确定所述特征曲线高频阈值;
对所述睡眠状态特征曲线进行动态地时频分析,根据所述特征曲线高频阈值,提取得到所述高频时频分析信号;
对所述高频时频分析信号进行动态地去噪声处理,分离得到所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线。
更优地,所述特征曲线高频阈值由所述睡眠状态特征曲线的睡眠生理状态特征线性组合融合的数据来源决定。
更优地,所述时频分析至少包括时频变换、时域滤波和频域滤波中的任一项;其中,所述时频变换至少包括标准傅里叶变换、快速傅里叶变换、短时傅里叶变换、S变换、Gabor变换、希尔伯特黄变换、Wigner-Ville分布、平滑Wigner-Ville分布、伪Wigner-Ville分布、平滑伪Wigner-Ville分布、Choi-Williams分布中的任一项,所述时域滤波至少包括平均值滤波、中值滤波、卡尔曼滤波和Savitzky-Golay滤波中的任一项,所述频域滤波至少包括低通滤波、高通滤波、带通滤波、带阻滤波中的任一项。
更优地,所述去噪声处理至少包括高斯滤波、均值滤波、傅里叶变换滤波、小波变换滤波和小波包变换滤波中的任一种方法,用于消除由于信号来自体表生理采集和分帧式的信号特征分析对睡眠状态刻画的不连续性,而在所述睡眠状态特征曲线中引入的高斯噪声或白噪声干扰。
更优地,所述根据所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线,计算得到睡眠可持续性指数和睡眠可持续性指数曲线的具体步骤还包括:
根据所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线,动态计算得到所述睡眠可持续性指数;
按照时序拼接所述睡眠可持续性指数,生成所述睡眠可持续性指数曲线。
更优地,所述睡眠可持续性指数的计算方法,具体为:
1)获取当前时帧对应的睡眠时相分期、所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线;
2)计算所述睡眠状态特征曲线和所述睡眠可持续性时频噪声曲线的差值,得到目标睡眠状态特征曲线。
3)分别对所述目标睡眠状态特征曲线和所述睡眠可持续性时频特征曲线先后做取绝对值和求均值运算,得到目标睡眠状态特征均值和睡眠可持续性TF特征均值;
4)根据所述目标睡眠状态特征均值和所述睡眠可持续性TF特征均值的相对变化量,计算得到睡眠可持续性节点特征系数;
5)根据所述睡眠可持续性节点特征系数、时频分析方法修正系数、用户个性修正系数和睡眠时相分期修正系数的均值,计算得到所述睡眠可持续性指数。
更优地,所述睡眠可持续性节点特征系数的一种计算公式,具体为:
;
其中,为所述睡眠可持续性节点特征系数,/>分别为所述目标睡眠状态特征均值和所述睡眠可持续性TF特征均值。
更优地,所述睡眠可持续性指数的一种计算公式,具体为:
;
其中,为所述睡眠可持续性指数,/>分别为所述睡眠可持续性节点特征系数、所述时频分析方法修正系数、所述用户个性修正系数和所述睡眠时相分期修正系数,且/>,/>2。
更优地,所述对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析,生成睡眠可持续性动态调节策略并完成用户睡眠过程的动态调节及效果评估,动态迭代优化睡眠可持续性的动态检测调节过程方法策略的具体步骤还包括:
对所述睡眠状态特征曲线进行动态地趋势预测分析,得到睡眠状态特征预测值;
对所述睡眠可持续性指数曲线进行动态地趋势预测分析,得到睡眠可持续性指数预测值;
根据所述睡眠状态特征预测值和所述睡眠可持续性指数预测值,按照睡眠可持续性检测调节周期和预设睡眠调节知识库,生成所述睡眠可持续性动态调节策略;
根据所述睡眠可持续性动态调节策略,连接并控制睡眠调节外围设备,对用户睡眠过程进行动态调节;
对调节效果进行动态跟踪评估,计算动态调节效果系数并生成动态调节效果曲线;
根据所述动态调节效果曲线和所述动态调节效果系数,动态迭代优化睡眠可持续性的动态检测调节过程方法策略。
更优地,所述趋势预测分析的方法至少包括AR、MA、ARMA、ARIMA、SARIMA、VAR、深度学习中的任一项。
更优地,所述睡眠可持续性动态调节策略至少包括睡眠场景、睡眠时相、调节方式、执行方式、调节方法、调节强度、调节时点、持续时间、目标调节值和装置控制参数;其中,所述调节方式至少包括声乐刺激、超声刺激、光刺激、电刺激、磁刺激、温度刺激、湿度刺激、触觉刺激和浓度调控中的任一方式,所述执行方式至少包括离体式和接触式中的任一方式。
更优地,所述睡眠调节外围设备至少包括声乐刺激设备、超声刺激设备、光刺激设备、电刺激设备、磁刺激设备、温度刺激设备、湿度刺激设备、触觉刺激设备和浓度调控设备中的任一项,并由具体的所述调节方式决定。
更优地,所述动态调节效果系数的一种计算方式,具体为:
;
其中,为所述动态调节效果系数,/>为趋势预测分析方法修正系数且,/>为动态调节后的所述睡眠时相分期修正系数,分别为动态调节前、动态调节后的所述睡眠可持续性指数。
更优地,所述动态调节效果曲线由按照时序拼接所述动态调节效果系数而生成得到。
更优地,所述动态检测调节过程方法策略至少包括去高斯噪声处理方法参数、时频分析方法参数、所述特征曲线高频阈值、趋势预测分析方法参数、所述睡眠可持续性动态调节策略和动态调节效果评估方法参数。
更优地,所述生成用户睡眠可持续性检测调节报告的具体步骤还包括:
根据所述睡眠时相曲线和所述睡眠可持续性指数曲线,计算不同睡眠时相下所述睡眠可持续性指数的分布特征,得到时相指数分布特征;
根据所述睡眠时相曲线和所述动态调节效果曲线,计算不同睡眠时相下所述动态调节效果系数的分布特征,得到时相效果分布特征;
计算所述睡眠可持续性指数曲线和所述动态调节效果曲线的相关性,得到动态调节综合效果系数;
按照预设报告生成周期,生成所述用户睡眠可持续性检测调节报告。
更优地,所述分布特征至少包括平均值、均方根、最大值、最小值、方差、标准差、变异系数、峰度和偏度中的任意一项。
更优地,所述相关性的计算方法至少包括相干性分析、皮尔逊相关分析、杰卡德相似分析、线性互信息分析、线性相关分析、欧氏距离分析、曼哈顿距离分析、切比雪夫距离分析中的任意一项。
更优地,所述用户睡眠可持续性检测调节报告至少包括所述用户关键生理信息、所述睡眠状态特征曲线、所述睡眠时相曲线、所述睡眠可持续性指数曲线、所述动态调节效果曲线、所述时相指数分布特征、所述时相效果分布特征、所述动态调节综合效果系数、睡眠可持续性检测调节总结。
更优地,所述初始化建立或持续更新用户个性睡眠可持续性数据库,更新睡眠可持续性通用数据库的具体步骤还包括:
动态收集检测调节过程数据和分析结果,初始化建立或持续动态更新所述用户个性睡眠可持续性数据库;
按照预设通用数据库更新周期,根据所述用户个性睡眠可持续性数据库和所述用户睡眠可持续性检测调节报告,将更新所述睡眠可持续性通用数据库的用户关键数据信息。
根据本发明的目的,本发明提出了一种基于时频分析的睡眠可持续性检测调节系统,包括以下模块:
状态检测分析模块,用于初始化睡眠可持续性动态检测调节基本方案,对用户睡眠生理信号进行采集和处理,得到睡眠状态特征曲线和睡眠时相曲线;
特征曲线提取模块,用于根据所述睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值并进行时频分析,得到高频时频分析信号并进行去噪声处理,得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线;
指数曲线提取模块,用于根据所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线,计算得到睡眠可持续性指数和睡眠可持续性指数曲线;
动态调节优化模块,用于对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析,生成睡眠可持续性动态调节策略并完成用户睡眠过程的动态调节及效果评估,动态迭代优化睡眠可持续性的动态检测调节过程方法策略;
统计分析报告模块,用于生成用户睡眠可持续性检测调节报告;
数据库更新模块,用于初始化建立或持续更新用户个性睡眠可持续性数据库,更新睡眠可持续性通用数据库;
数据运行管理模块,用于对所述系统的所有数据进行可视化管理、统一存储和运行管理。
更优地,所述状态检测分析模块还包括以下功能单元:
方案初始化单元,用于获取用户关键生理信息,从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选并初始化所述睡眠可持续性动态检测调节基本方案;
信号采集处理单元,用于对用户睡眠生理信号进行动态地监测采集和信号处理,得到睡眠生理状态数据;
信号特征分析单元,用于对所述睡眠生理状态数据进行动态地特征分析和特征选择,生成所述睡眠状态特征曲线;
睡眠时相识别单元,用于对所述睡眠生理状态数据进行动态地睡眠时相分析,生成所述睡眠时相曲线。
更优地,所述特征曲线提取模块还包括以下功能单元:
信号阈值确定单元,用于根据所述睡眠状态特征曲线的曲线生成属性,确定所述特征曲线高频阈值;
高频分量提取单元,用于对所述睡眠状态特征曲线进行动态地时频分析,根据所述特征曲线高频阈值,提取得到所述高频时频分析信号;
特征曲线生成单元,用于对所述高频时频分析信号进行动态地去噪声处理,分离得到所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线。
更优地,所述指数曲线提取模块还包括以下功能单元:
特征指数计算单元,用于根据所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线,动态计算得到所述睡眠可持续性指数;
指数曲线生成单元,用于按照时序拼接所述睡眠可持续性指数,生成所述睡眠可持续性指数曲线。
更优地,所述动态调节优化模块还包括以下功能单元:
状态特征预测单元,用于对所述睡眠状态特征曲线进行动态地趋势预测分析,得到睡眠状态特征预测值;
特征指数预测单元,用于对所述睡眠可持续性指数曲线进行动态地趋势预测分析,得到睡眠可持续性指数预测值;
动态策略生成单元,用于根据所述睡眠状态特征预测值和所述睡眠可持续性指数预测值,按照睡眠可持续性检测调节周期和预设睡眠调节知识库,生成所述睡眠可持续性动态调节策略;
动态调节执行单元,用于根据所述睡眠可持续性动态调节策略,连接并控制睡眠调节外围设备,对用户睡眠过程进行动态调节;
调节效果评估单元,用于对调节效果进行动态跟踪评估,计算动态调节效果系数并生成动态调节效果曲线;
过程方法优化单元,用于根据所述动态调节效果曲线和所述动态调节效果系数,动态迭代优化睡眠可持续性的动态检测调节过程方法策略。
更优地,所述统计分析报告模块还包括以下功能单元:
指数分布分析单元,用于根据所述睡眠时相曲线和所述睡眠可持续性指数曲线,计算不同睡眠时相下所述睡眠可持续性指数的分布特征,得到时相指数分布特征;
效果分布分析单元,用于根据所述睡眠时相曲线和所述动态调节效果曲线,计算不同睡眠时相下所述动态调节效果系数的分布特征,得到时相效果分布特征;
综合效果评估单元,用于计算所述睡眠可持续性指数曲线和所述动态调节效果曲线的相关性,得到动态调节综合效果系数;
用户报告生成单元,用于按照预设报告生成周期,生成所述用户睡眠可持续性检测调节报告;
用户报告管理单元,用于对所述用户睡眠可持续性检测调节报告的格式输出、展现形式进行统一管理。
更优地,所述数据库更新模块还包括以下功能单元:
个性数据库更新单元,用于动态收集检测调节过程数据和分析结果,初始化建立或持续动态更新所述用户个性睡眠可持续性数据库;
通用数据库更新单元,用于按照预设通用数据库更新周期,根据所述用户个性睡眠可持续性数据库和所述用户睡眠可持续性检测调节报告,将更新所述睡眠可持续性通用数据库的用户关键数据信息。
更优地,所述数据运行管理模块还包括以下功能单元:
用户信息管理单元,用于用户基本信息的登记输入、编辑、查询、输出和删除;
数据可视化管理单元,用于对所述系统中所有数据的可视化展现管理;
数据存储管理单元,用于对所述系统中所有数据的统一存储管理;
数据运营管理单元,用于对所述系统中所有数据的备份、迁移和导出。
根据本发明的目的,本发明提出了一种基于时频分析的睡眠可持续性检测调节装置,包括以下模组:
状态检测分析模组,用于初始化睡眠可持续性动态检测调节基本方案,对用户睡眠生理信号进行采集和处理,得到睡眠状态特征曲线和睡眠时相曲线;
特征曲线提取模组,用于根据所述睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值并进行时频分析,得到高频时频分析信号并进行去噪声处理,得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线;
指数曲线提取模组,用于根据所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线,计算得到睡眠可持续性指数和睡眠可持续性指数曲线;
动态调节优化模组,用于对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析,生成睡眠可持续性动态调节策略并完成用户睡眠过程的动态调节及效果评估,动态迭代优化睡眠可持续性的动态检测调节过程方法策略;
统计分析报告模组,用于生成用户睡眠可持续性检测调节报告;
数据库更新模组,用于初始化建立或持续更新用户个性睡眠可持续性数据库,更新睡眠可持续性通用数据库;
数据可视化模组,用于对所述装置中所有过程数据和结果数据的统一可视化展示管理;
数据管理中心模组,用于对所述装置中所有过程数据和结果数据的统一存储和数据运营管理。
本发明在申请人的在先研究基础上进一步优化睡眠可持续性指数的具体量化设计和计算方式,进一步优化睡眠状态特征曲线的提取方式,进一步将去噪声处理、时频分析方法应用于睡眠可持续性特征的提取,使得对用户睡眠可持续性的检测量化更加全面细致、适应度广、准确度高和灵敏度高;还进一步提供了动态调节效果系数的计算方案,进一步提供了用户个性睡眠可持续性数据库和睡眠可持续性通用数据库的建立、更新和应用机制,进一步提供了检测量化到动态调节的反向反馈应用框架,使得对于检测量化和动态调节过程的协同控制有了有力依据,实现检测调节的过程动态优化,持续提高检测量化效率和干预调节效果。本发明能够提供更加科学高效的睡眠可持续性检测量化及动态调节的实施方法和落地方案,能够赋能相关睡眠量化或调节的产品及服务,满足不同用户场景需求,辅助用户睡眠。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明技术方案的进一步理解,并且构成说明书的一部分,与本申请的实施例一起用于解释本发明的技术方案,并不构成对本发明技术方案的限制。
图1是本发明一实施例所提供的一种基于时频分析的睡眠可持续性检测调节方法的流程步骤示意图;
图2是本发明一实施例所提供的一种基于时频分析的睡眠可持续性检测调节系统的模块组成示意图;
图3是本发明一实施例所提供的一种基于时频分析的睡眠可持续性检测调节装置的模组构成示意图。
具体实施方式
为了更清楚地说明本发明的目的和技术方案,下面将结合本发明申请实施例中的附图,对本发明进行进一步介绍说明。显而易见地,下面描述的实施例仅是本发明一部分实施例,而不是全部的实施例。在没有创造性劳动前提下,本领域普通技术人员基于本发明的实施例所得到的其他实施例,都应属于本发明的保护范围。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互任意组合。
申请人发现,人类睡眠是个非平稳时序过程,而睡眠过程中多次觉醒睡眠状态突变震荡是非平稳的重要因素或诱因,如何跟踪多变状态、意外状态和突变状态是睡眠可持续性或连续性检测调节的一个基本挑战。同时由于头皮脑电EEG、心电ECG等体表生理采集技术的原理限制及对应的传统信号特征分析提取过程,睡眠状态量化上存在非常多的计算误差和引入噪声,尤其是面对睡眠可持续性的评价上。本技术方案通过敏感特征的选择、去噪声处理和时频分析等方法,能够提取细致和精确的睡眠状态特征曲线,完成睡眠可持续性的高鲁棒性的量化设计;进而完成睡眠可持续性的动态调节和效果评估,进一步实现检测调节的反向动态优化,使得整体检测调节方案具备更好的场景适应性和技术应用效果。
结合图1所示,本发明实施例提供的一种基于时频分析的睡眠可持续性检测调节方法,包括以下步骤:
P100:初始化睡眠可持续性动态检测调节基本方案,对用户睡眠生理信号进行采集和处理,得到睡眠状态特征曲线和睡眠时相曲线。
第一步、获取用户关键生理信息,从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选并初始化睡眠可持续性动态检测调节基本方案。
本实施例中,用户关键生理信息至少包括性别、年龄、生理健康状态和心理精神状态。
本实施例中,睡眠可持续性通用数据库具体为由规模的不同健康状态人群的睡眠可持续性的检测量化及动态调节数据库,至少包括用户关键生理信息、睡眠可持续性检测调节周期、睡眠可持续性检测量化过程方法参数、睡眠可持续性动态调节过程方法参数或策略、睡眠可持续性指数曲线、睡眠可持续性调节效果曲线和睡眠时相曲线。
本实施例中,用户个性睡眠可持续性数据库至少包括用户关键生理信息、睡眠状态特征曲线、睡眠时相曲线、睡眠可持续性指数曲线、动态调节效果曲线、时相指数分布特征、时相效果分布特征、动态调节综合效果系数、时频分析方法、趋势预测分析方法、睡眠可持续性动态调节策略、睡眠可持续性动态检测调节基本方案。
本实施例中,睡眠可持续性动态检测调节基本方案至少包括睡眠可持续性检测量化过程方案、睡眠可持续性动态调节过程方案;若用户为新用户,睡眠可持续性动态检测调节基本方案从睡眠可持续性通用数据库中获取,反之则从用户个性睡眠可持续性数据库中获取。
在实际应用场景中,首先确定当前用户是新用户还是老用户,来选择睡眠可持续性动态检测调节基本方案的获取来源;当前用户的当前状态,包括生理健康状态和心理精神状态,是睡眠可持续性动态检测调节基本方案选择的重要考量。
第二步、对用户睡眠生理信号进行动态地监测采集和信号处理,得到睡眠生理状态数据。
本实施例中,睡眠生理信号至少包括脑电信号、心电信号、呼吸信号中的任一项。信号处理至少包括重采样、重参考、去伪迹、信号矫正、降噪、工频陷波、带通滤波、平滑处理和时帧分割;其中,时帧分割是指根据信号的采样率,以预设分帧步长对信号数据进行预设分帧时长窗口的连续滑动分割。
本实施例中,采集监测用户睡眠过程的脑电信号和呼吸信号作为睡眠生理信号,来陈述本技术方案的具体实施过程。首先,通过脑电图机对用户睡眠脑电进行采集记录,采样率为1024Hz,记录电极为F3、F4、C3和C4,参考电极M1、M2;对脑电信号进行统一信号处理,包括以M1和M2进行左右交叉重参考、去伪迹、信号矫正、小波降噪、50Hz倍频工频陷波、1.0-80Hz带通滤波,得到纯净脑电信号。通过胸腹压力带对用户胸部腹部的呼吸信号进行采集记录,采样率64Hz;对呼吸信号进行统一信号处理,包括降噪、去伪迹、信号矫正和0.01-10Hz带通滤波,得到纯净呼吸信号。最后,以预设时间窗口长度30s和预设时间平移步长15s对纯净脑电信号和纯净呼吸号进行连续滑动分割,得到睡眠生理状态数据。
本实施例中,预设时间窗口长度30s也是睡眠可持续性检测调节周期。
第三步、对睡眠生理状态数据进行动态地特征分析和特征选择,生成睡眠状态特征曲线。
本实施例中,特征分析至少包括数值分析、包络分析、时频分析、熵分析、分形分析和复杂度分析中的任一项。睡眠状态特征曲线具体由不同睡眠生理状态数据经特征分析得到不同睡眠生理状态特征进行线性组合融合而得到,用于描述用户在不同睡眠时期、不同睡眠时相、不同睡眠状态下的生理状态特征的连续变化曲线;其中,睡眠时期至少包括入睡前、睡眠中和睡眠后,睡眠时相至少包括清醒觉醒期、浅睡眠期、深睡眠期和快速眼动睡眠期。
本实施例中,以对睡眠可持续性敏感度较高的指标特征来生成睡眠状态特征曲线。首先,逐帧对睡眠生理状态数据进行特征提取:脑电数据进行时频分析(频带功率、频带功率占比)、熵分析(近似熵)和复杂度分析(LZC指数);呼吸数据进行呼吸频率和呼吸深度的数值分析(平均值、变异系数)。通过特征选择,将F3-M2通道的δ节律(1-4Hz)频带功率占比、取负后再归一化的近似熵、取负后再归一化的呼吸深度均值进行直接加和,得到睡眠状态特征曲线。
第四步、对睡眠生理状态数据进行动态地睡眠时相分析,生成睡眠时相曲线。
本实施例中,睡眠时相曲线的生成方法具体为:
1)通过机器学习算法对规模睡眠用户样本的睡眠生理状态数据及其对应睡眠时相分期数据进行学习训练和数据建模,得到睡眠时相识别模型;
2)将当前用户的睡眠生理状态数据输入睡眠时相识别模型,得到所对应的睡眠时相分期并按时序生成睡眠时相曲线。
P200:根据所述睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值并进行时频分析,得到高频时频分析信号并进行去噪声处理,得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线。
第一步、根据睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值。
本实施例中,特征曲线高频阈值由睡眠状态特征曲线的睡眠生理状态特征线性组合融合的数据来源决定。
本实施例中,以δ节律频带功率占比、取负后再归一化的近似熵、取负后再归一化的呼吸深度均值进行直接加和得到睡眠状态特征曲线,且分帧处理的参数为预设时间窗口长度30s和预设时间平移步长15s,所以以0.008Hz高通截止频率作为特征曲线高频阈值。
第二步、对睡眠状态特征曲线进行动态地时频分析,根据特征曲线高频阈值,提取得到高频时频分析信号。
本实施例中,时频分析至少包括时频变换、时域滤波和频域滤波中的任一项;其中,时频变换至少包括标准傅里叶变换、快速傅里叶变换、短时傅里叶变换、S变换、Gabor变换、希尔伯特黄变换、Wigner-Ville分布、平滑Wigner-Ville分布、伪Wigner-Ville分布、平滑伪Wigner-Ville分布、Choi-Williams分布中的任一项,时域滤波至少包括平均值滤波、中值滤波、卡尔曼滤波和Savitzky-Golay滤波中的任一项,频域滤波至少包括低通滤波、高通滤波、带通滤波、带阻滤波中的任一项。
本实施例中,通过0.008Hz高通截止频率、高通滤波器对睡眠状态特征曲线进行时频分析,提取得到高频时频分析信号。
第三步、对高频时频分析信号进行动态地去噪声处理,分离得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线。
本实施例中,去噪声处理至少包括高斯滤波、均值滤波、傅里叶变换滤波、小波变换滤波和小波包变换滤波中的任一种方法,用于消除由于信号来自体表生理采集和分帧式的信号特征分析对睡眠状态刻画的不连续性,而在所述睡眠状态特征曲线中引入的高斯噪声或白噪声干扰。
本实施例中,通过高斯滤波从高频时频分析信号中分离出噪声信号和特征信号,即睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线。
P300:根据所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线,计算得到睡眠可持续性指数和睡眠可持续性指数曲线。
第一步、根据睡眠状态特征曲线、睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线,动态计算得到睡眠可持续性指数。
本实施例中,睡眠可持续性指数的计算方法,具体为:
1)获取当前时帧对应的睡眠时相分期、睡眠状态特征曲线、睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线;
2)计算睡眠状态特征曲线和睡眠可持续性时频噪声曲线的差值,得到目标睡眠状态特征曲线。
3)分别对目标睡眠状态特征曲线和睡眠可持续性时频特征曲线先后做取绝对值和求均值运算,得到目标睡眠状态特征均值和睡眠可持续性TF特征均值;
4)根据目标睡眠状态特征均值和睡眠可持续性TF特征均值的相对变化量,计算得到睡眠可持续性节点特征系数;
5)根据睡眠可持续性节点特征系数、时频分析方法修正系数、用户个性修正系数和睡眠时相分期修正系数的均值,计算得到睡眠可持续性指数。
本实施例中,睡眠可持续性节点特征系数的一种计算公式,具体为:
;
其中,为睡眠可持续性节点特征系数,/>分别为目标睡眠状态特征均值和睡眠可持续性TF特征均值。
本实施例中,睡眠可持续性指数的一种计算公式,具体为:
;
其中,为睡眠可持续性指数,/>分别为睡眠可持续性节点特征系数、时频分析方法修正系数、用户个性修正系数和睡眠时相分期修正系数,且/>,/>2。
在实际应用场景中,通常清醒觉醒期、浅睡眠期、深睡眠期和快速眼动睡眠期的修正系数,依次为 1.0、1.6、2.0、1.4;用户个性修正系数/>主要跟用户年龄和性别相关,通常女性比男性要小,年龄越大系数越小;时频分析方法修正系数/>则可以根据睡眠状态特征曲线的生成方式、去噪声处理的方法参数等实际情况设置不同方法的系数。更重要的是,在后续的动态迭代优化睡眠可持续性的动态检测调节过程方法策略中可以进行动态优化调整。
第二步、按照时序拼接睡眠可持续性指数,生成睡眠可持续性指数曲线。
本实施例中,睡眠可持续性指数曲线由按照时序拼接睡眠可持续性指数而生成得到。
P400:对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析,生成睡眠可持续性动态调节策略并完成用户睡眠过程的动态调节及效果评估,动态迭代优化睡眠可持续性的动态检测调节过程方法策略。
第一步、对睡眠状态特征曲线进行动态地趋势预测分析,得到睡眠状态特征预测值。
本实施例中,趋势预测分析的方法至少包括AR、MA、ARMA、ARIMA、SARIMA、VAR、深度学习中的任一项。
本实施例中,用ARIMA方法来完成对睡眠状态特征曲线的预测。
第二步、对睡眠可持续性指数曲线进行动态地趋势预测分析,得到睡眠可持续性指数预测值。
本实施例中,用ARIMA方法来完成对睡眠可持续性指数曲线的预测。
第三步、根据睡眠状态特征预测值和睡眠可持续性指数预测值,按照睡眠可持续性检测调节周期和预设睡眠调节知识库,生成睡眠可持续性动态调节策略。
本实施例中,睡眠可持续性动态调节策略至少包括睡眠场景、睡眠时相、调节方式、执行方式、调节方法、调节强度、调节时点、持续时间、目标调节值和装置控制参数;其中,调节方式至少包括声乐刺激、超声刺激、光刺激、电刺激、磁刺激、温度刺激、湿度刺激、触觉刺激和浓度调控中的任一方式,执行方式至少包括离体式和接触式中的任一方式。
在实际应用场景中,需要根据用户具体情况来制定不同睡眠可持续性动态调节策略。如选择离体式的声乐刺激、光刺激、温度刺激、湿度刺激、浓度调控等方式,能够减少对用户睡眠过程的干扰;但通常情况下,接触式的超声刺激、电刺激、磁刺激、触觉刺激等方式的调节效率更高、效果也更好。
第四步、根据睡眠可持续性动态调节策略,连接并控制睡眠调节外围设备,对用户睡眠过程进行动态调节。
本实施例中,睡眠调节外围设备至少包括声乐刺激设备、超声刺激设备、光刺激设备、电刺激设备、磁刺激设备、温度刺激设备、湿度刺激设备、触觉刺激设备和浓度调控设备中的任一项,并由具体的调节方式决定。
第五步、对调节效果进行动态跟踪评估,计算动态调节效果系数并生成动态调节效果曲线。
本实施例中,动态调节效果系数的一种计算方式,具体为:
;
其中,为动态调节效果系数,/>为趋势预测分析方法修正系数且,/>为动态调节后的睡眠时相分期修正系数,/>分别为动态调节前、动态调节后的睡眠可持续性指数。
本实施例中,动态调节效果曲线由按照时序拼接动态调节效果系数而生成得到。
在实际应用场景中,通常趋势预测分析方法AR、MA、ARMA、ARIMA、SARIMA、VAR、深度学习的方法修正系数,依次为0.70、0.75、0.80、0.85、0.85、0.9、1.0。更重要的是,在后续的动态迭代优化睡眠可持续性的动态检测调节过程方法策略中可以进行动态优化调整。
第六步、根据动态调节效果曲线和动态调节效果系数,动态迭代优化睡眠可持续性的动态检测调节过程方法策略。
本实施例中,动态检测调节过程方法策略至少包括去高斯噪声处理方法参数、时频分析方法参数、特征曲线高频阈值、趋势预测分析方法参数、睡眠可持续性动态调节策略和动态调节效果评估方法参数。
在实际应用场景中,由动态调节效果系数来反向优化睡眠可持续性检测量化和动态调节的方法参数,实现了检测量化和动态调节的协同控制,以调节效果为目标导向,持续提高检测量化和动态调节的效率效果。
P500:生成用户睡眠可持续性检测调节报告。
第一步、根据睡眠时相曲线和睡眠可持续性指数曲线,计算不同睡眠时相下睡眠可持续性指数的分布特征,得到时相指数分布特征。
本实施例中,分布特征至少包括平均值、均方根、最大值、最小值、方差、标准差、变异系数、峰度和偏度中的任意一项。
本实施例中,时相指数分布特征主要包括清醒觉醒期、浅睡眠期、深睡眠期和快速眼动睡眠期的、睡眠可持续性指数的平均值、最大值、最小值的统计分布情况,来观察分析不同睡眠时相下睡眠可持续性指数的主要特征。
第二步、根据睡眠时相曲线和动态调节效果曲线,计算不同睡眠时相下动态调节效果系数的分布特征,得到时相效果分布特征。
本实施例中,时相效果分布特征主要包括清醒觉醒期、浅睡眠期、深睡眠期和快速眼动睡眠期的、动态调节效果的平均值、最大值、最小值的统计分布情况,来观察分析不同睡眠时相下动态调节效果,判断干预难易程度。
第三步、计算睡眠可持续性指数曲线和动态调节效果曲线的相关性,得到动态调节综合效果系数。
本实施例中,相关性的计算方法至少包括相干性分析、皮尔逊相关分析、杰卡德相似分析、线性互信息分析、线性相关分析、欧氏距离分析、曼哈顿距离分析、切比雪夫距离分析中的任意一项。
在实际应用场景中,皮尔逊相关分析能够快速准确地评价睡眠可持续性指数曲线和动态调节效果曲线的相关性,来量化用户睡眠生理对于动态调节的适应性和应激程度,为动态调节策略的制定提供依据。
第四步、按照预设报告生成周期,生成用户睡眠可持续性检测调节报告。
本实施例中,用户睡眠可持续性检测调节报告至少包括用户关键生理信息、睡眠状态特征曲线、睡眠时相曲线、睡眠可持续性指数曲线、动态调节效果曲线、时相指数分布特征、时相效果分布特征、动态调节综合效果系数、睡眠可持续性检测调节总结。
在实际应用场景中,可以根据用户情况设置不同的报告生成周期,来满足对用户睡眠可持续性的观察、分析和总结,可以对用户睡眠可持续性检测调节报告进行多种格式输出和多种形式展现,适应不同场景需要。
P600:初始化建立或持续更新用户个性睡眠可持续性数据库,更新睡眠可持续性通用数据库。
第一步、动态收集检测调节过程数据和分析结果,初始化建立或持续动态更新用户个性睡眠可持续性数据库。
本实施例中,如用户是新用户的首次检测调节,那么需要在第一步检测调节周期结束时就开始初始化建立用户个性睡眠可持续性数据库,为后续检测量化和动态调节的动态过程优化提供基础。
第二步、按照预设通用数据库更新周期,根据用户个性睡眠可持续性数据库和用户睡眠可持续性检测调节报告,将更新睡眠可持续性通用数据库的用户关键数据信息。
本实施例中,睡眠可持续性通用数据库的持续更新,能够保障不同健康状态人群的睡眠可持续性的检测量化及动态调节的关键数据的积累,有利于睡眠可持续性动态检测调节基本方案与新用户的高效匹配,提高新用户首次检测调节的准确性和有效性。
结合图2所示,本发明实施例提供的一种基于时频分析的睡眠可持续性检测调节系统,所述系统被构造以用于执行上述各个方法步骤。所述系统包括如下模块:
状态检测分析模块S100,用于初始化睡眠可持续性动态检测调节基本方案,对用户睡眠生理信号进行采集和处理,得到睡眠状态特征曲线和睡眠时相曲线;
特征曲线提取模块S200,用于根据睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值并进行时频分析,得到高频时频分析信号并进行去噪声处理,得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线;
指数曲线提取模块S300,用于根据睡眠状态特征曲线、睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线,计算得到睡眠可持续性指数和睡眠可持续性指数曲线;
动态调节优化模块S400,用于对睡眠状态特征曲线和睡眠可持续性指数曲线进行趋势预测分析,生成睡眠可持续性动态调节策略并完成用户睡眠过程的动态调节及效果评估,动态迭代优化睡眠可持续性的动态检测调节过程方法策略;
统计分析报告模块S500,用于生成用户睡眠可持续性检测调节报告;
数据库更新模块S600,用于初始化建立或持续更新用户个性睡眠可持续性数据库,更新睡眠可持续性通用数据库;
数据运行管理模块S700,用于对系统的所有数据进行可视化管理、统一存储和运行管理。
本实施例中,状态检测分析模块S100还包括以下功能单元:
方案初始化单元,用于获取用户关键生理信息,从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选并初始化睡眠可持续性动态检测调节基本方案;
信号采集处理单元,用于对用户睡眠生理信号进行动态地监测采集和信号处理,得到睡眠生理状态数据;
信号特征分析单元,用于对睡眠生理状态数据进行动态地特征分析和特征选择,生成睡眠状态特征曲线;
睡眠时相识别单元,用于对睡眠生理状态数据进行动态地睡眠时相分析,生成睡眠时相曲线。
本实施例中,特征曲线提取模块S200还包括以下功能单元:
信号阈值确定单元,用于根据睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值;
高频分量提取单元,用于对睡眠状态特征曲线进行动态地时频分析,根据特征曲线高频阈值,提取得到高频时频分析信号;
特征曲线生成单元,用于对高频时频分析信号进行动态地去噪声处理,分离得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线。
本实施例中,指数曲线提取模块S300还包括以下功能单元:
特征指数计算单元,用于根据睡眠状态特征曲线、睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线,动态计算得到睡眠可持续性指数;
指数曲线生成单元,用于按照时序拼接睡眠可持续性指数,生成睡眠可持续性指数曲线。
本实施例中,动态调节优化模块S400还包括以下功能单元:
状态特征预测单元,用于对睡眠状态特征曲线进行动态地趋势预测分析,得到睡眠状态特征预测值;
特征指数预测单元,用于对睡眠可持续性指数曲线进行动态地趋势预测分析,得到睡眠可持续性指数预测值;
动态策略生成单元,用于根据睡眠状态特征预测值和睡眠可持续性指数预测值,按照睡眠可持续性检测调节周期和预设睡眠调节知识库,生成睡眠可持续性动态调节策略;
动态调节执行单元,用于根据睡眠可持续性动态调节策略,连接并控制睡眠调节外围设备,对用户睡眠过程进行动态调节;
调节效果评估单元,用于对调节效果进行动态跟踪评估,计算动态调节效果系数并生成动态调节效果曲线;
过程方法优化单元,用于根据动态调节效果曲线和动态调节效果系数,动态迭代优化睡眠可持续性的动态检测调节过程方法策略。
本实施例中,统计分析报告模块S500还包括以下功能单元:
指数分布分析单元,用于根据睡眠时相曲线和睡眠可持续性指数曲线,计算不同睡眠时相下睡眠可持续性指数的分布特征,得到时相指数分布特征;
效果分布分析单元,用于根据睡眠时相曲线和动态调节效果曲线,计算不同睡眠时相下动态调节效果系数的分布特征,得到时相效果分布特征;
综合效果评估单元,用于计算睡眠可持续性指数曲线和动态调节效果曲线的相关性,得到动态调节综合效果系数;
用户报告生成单元,用于按照预设报告生成周期,生成用户睡眠可持续性检测调节报告;
用户报告管理单元,用于对用户睡眠可持续性检测调节报告的格式输出、展现形式进行统一管理。
本实施例中,数据库更新模块S600还包括以下功能单元:
个性数据库更新单元,用于动态收集检测调节过程数据和分析结果,初始化建立或持续动态更新用户个性睡眠可持续性数据库;
通用数据库更新单元,用于按照预设通用数据库更新周期,根据用户个性睡眠可持续性数据库和用户睡眠可持续性检测调节报告,将更新睡眠可持续性通用数据库的用户关键数据信息。
本实施例中,数据运行管理模块S700还包括以下功能单元:
用户信息管理单元,用于用户基本信息的登记输入、编辑、查询、输出和删除;
数据可视化管理单元,用于对系统中所有数据的可视化展现管理;
数据存储管理单元,用于对系统中所有数据的统一存储管理;
数据运营管理单元,用于对系统中所有数据的备份、迁移和导出。
结合图3所示,本发明实施例提供的一种基于时频分析的睡眠可持续性检测调节装置,包括以下模组:
状态检测分析模组M100,用于初始化睡眠可持续性动态检测调节基本方案,对用户睡眠生理信号进行采集和处理,得到睡眠状态特征曲线和睡眠时相曲线;
特征曲线提取模组M200,用于根据睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值并进行时频分析,得到高频时频分析信号并进行去噪声处理,得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线;
指数曲线提取模组M300,用于根据睡眠状态特征曲线、睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线,计算得到睡眠可持续性指数和睡眠可持续性指数曲线;
动态调节优化模组M400,用于对睡眠状态特征曲线和睡眠可持续性指数曲线进行趋势预测分析,生成睡眠可持续性动态调节策略并完成用户睡眠过程的动态调节及效果评估,动态迭代优化睡眠可持续性的动态检测调节过程方法策略;
统计分析报告模组M500,用于生成用户睡眠可持续性检测调节报告;
数据库更新模组M600,用于初始化建立或持续更新用户个性睡眠可持续性数据库,更新睡眠可持续性通用数据库;
数据可视化模组M700,用于对装置中所有过程数据和结果数据的统一可视化展示管理;
数据管理中心模组M800,用于对装置中所有过程数据和结果数据的统一存储和数据运营管理。
所述装置被构造以用于对应执行图1的方法中的各个步骤,在此不再赘述。
本发明还提供了可编程的各类处理器(FPGA、ASIC或其他集成电路),所述处理器用于运行程序,其中,所述程序运行时执行上述实施例中的步骤。
本发明还提供了对应的计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述存储器执行所述程序时实现上述实施例中的步骤。
虽然本发明所揭露的实施方式如上,但所述的 内容仅为便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属领域内的技术人员,在不脱离本发明所揭露的精神和原则的前提下,可以在实施的形式及细节上进行任何的修改与变化、等同替换等,这些都属于本发明的保护范围。因此,本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (38)
1.一种基于时频分析的睡眠可持续性检测调节方法,其特征在于,包括以下步骤:
初始化睡眠可持续性动态检测调节基本方案,对用户睡眠生理信号进行采集和处理,得到睡眠状态特征曲线和睡眠时相曲线;
根据所述睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值并进行时频分析,得到高频时频分析信号并进行去噪声处理,得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线;
根据所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线,计算得到睡眠可持续性指数和睡眠可持续性指数曲线;
对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析,生成睡眠可持续性动态调节策略并完成用户睡眠过程的动态调节及效果评估,动态迭代优化睡眠可持续性的动态检测调节过程方法策略;
生成用户睡眠可持续性检测调节报告;
所述睡眠可持续性指数的计算方法,具体为:
1)获取当前时帧对应的睡眠时相分期、所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线;
2)计算所述睡眠状态特征曲线和所述睡眠可持续性时频噪声曲线的差值,得到目标睡眠状态特征曲线;
3)分别对所述目标睡眠状态特征曲线和所述睡眠可持续性时频特征曲线先后做取绝对值和求均值运算,得到目标睡眠状态特征均值和睡眠可持续性TF特征均值;
4)根据所述目标睡眠状态特征均值和所述睡眠可持续性TF特征均值的相对变化量,计算得到睡眠可持续性节点特征系数;
5)根据时频分析方法修正系数、用户个性修正系数和睡眠时相分期修正系数的均值对所述睡眠可持续性节点特征系数进行修正,计算得到所述睡眠可持续性指数;
所述睡眠可持续性节点特征系数的一种计算公式,具体为:
;
其中,为所述睡眠可持续性节点特征系数,/>分别为所述目标睡眠状态特征均值和所述睡眠可持续性TF特征均值;
所述睡眠可持续性指数的一种计算公式,具体为:
;
其中,为所述睡眠可持续性指数,/>分别为所述睡眠可持续性节点特征系数、所述时频分析方法修正系数、所述用户个性修正系数和所述睡眠时相分期修正系数,且/>,/>2。
2.如权利要求1所述的方法,其特征在于,所述初始化睡眠可持续性动态检测调节基本方案,对用户睡眠生理信号进行采集和处理,得到睡眠状态特征曲线和睡眠时相曲线的具体步骤还包括:
获取用户关键生理信息,从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选并初始化所述睡眠可持续性动态检测调节基本方案;
对用户睡眠生理信号进行动态地监测采集和信号处理,得到睡眠生理状态数据;
对所述睡眠生理状态数据进行动态地特征分析和特征选择,生成所述睡眠状态特征曲线;
对所述睡眠生理状态数据进行动态地睡眠时相分析,生成所述睡眠时相曲线。
3.如权利要求2所述的方法,其特征在于:所述用户关键生理信息至少包括性别、年龄、生理健康状态和心理精神状态。
4.如权利要求2或3所述的方法,其特征在于:所述睡眠可持续性通用数据库具体为由规模的不同健康状态人群的睡眠可持续性的检测量化及动态调节数据库,至少包括所述用户关键生理信息、睡眠可持续性检测调节周期、睡眠可持续性检测量化过程方法参数、睡眠可持续性动态调节过程方法参数或策略、睡眠可持续性指数曲线、睡眠可持续性调节效果曲线和睡眠时相曲线。
5.如权利要求2或3所述的方法,其特征在于:所述用户个性睡眠可持续性数据库至少包括所述用户关键生理信息、睡眠状态特征曲线、睡眠时相曲线、睡眠可持续性指数曲线、动态调节效果曲线、时相指数分布特征、时相效果分布特征、动态调节综合效果系数、时频分析方法、趋势预测分析方法、睡眠可持续性动态调节策略、所述睡眠可持续性动态检测调节基本方案。
6.如权利要求2或3所述的方法,其特征在于:所述睡眠可持续性动态检测调节基本方案至少包括睡眠可持续性检测量化过程方案、睡眠可持续性动态调节过程方案;若用户为新用户,所述睡眠可持续性动态检测调节基本方案从所述睡眠可持续性通用数据库中获取,反之则从所述用户个性睡眠可持续性数据库中获取。
7.如权利要求2所述的方法,其特征在于,所述睡眠生理信号至少包括脑电信号、心电信号、呼吸信号中的任一项。
8.如权利要求7所述的方法,其特征在于,所述信号处理至少包括重采样、重参考、去伪迹、信号矫正、降噪、工频陷波、带通滤波、平滑处理和时帧分割;其中,所述时帧分割是指根据信号的采样率,以预设分帧步长对信号数据进行预设分帧时长窗口的连续滑动分割。
9.如权利要求7所述的方法,其特征在于,所述特征分析至少包括数值分析、包络分析、时频分析、熵分析、分形分析和复杂度分析中的任一项。
10.如权利要求7所述的方法,其特征在于:所述睡眠状态特征曲线具体由不同睡眠生理状态数据经所述特征分析得到不同睡眠生理状态特征进行线性组合融合而得到,用于描述用户在不同睡眠时期、不同睡眠时相、不同睡眠状态下的生理状态特征的连续变化曲线;其中,所述睡眠时期至少包括入睡前、睡眠中和睡眠后,所述睡眠时相至少包括清醒觉醒期、浅睡眠期、深睡眠期和快速眼动睡眠期。
11.如权利要求7-10任一项所述的方法,其特征在于:所述睡眠时相曲线的生成方法具体为:
1)通过机器学习算法对规模睡眠用户样本的所述睡眠生理状态数据及其对应睡眠时相分期数据进行学习训练和数据建模,得到睡眠时相识别模型;
2)将当前用户的所述睡眠生理状态数据输入所述睡眠时相识别模型,得到所对应的睡眠时相分期并按时序生成所述睡眠时相曲线。
12.如权利要求1或2或7任一项所述的方法,其特征在于:所述根据所述睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值并进行时频分析,得到高频时频分析信号并进行去噪声处理,得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线的具体步骤还包括:
根据所述睡眠状态特征曲线的曲线生成属性,确定所述特征曲线高频阈值;
对所述睡眠状态特征曲线进行动态地时频分析,根据所述特征曲线高频阈值,提取得到所述高频时频分析信号;
对所述高频时频分析信号进行动态地去噪声处理,分离得到所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线。
13.如权利要求12所述的方法,其特征在于,所述特征曲线高频阈值由所述睡眠状态特征曲线的睡眠生理状态特征线性组合融合的数据来源决定。
14.如权利要求12所述的方法,其特征在于,所述时频分析至少包括时频变换、时域滤波和频域滤波中的任一项;其中,所述时频变换至少包括标准傅里叶变换、快速傅里叶变换、短时傅里叶变换、S变换、Gabor变换、希尔伯特黄变换、Wigner-Ville分布、平滑Wigner-Ville分布、伪Wigner-Ville分布、平滑伪Wigner-Ville分布、Choi-Williams分布中的任一项,所述时域滤波至少包括平均值滤波、中值滤波、卡尔曼滤波和Savitzky-Golay滤波中的任一项,所述频域滤波至少包括低通滤波、高通滤波、带通滤波、带阻滤波中的任一项。
15.如权利要求12所述的方法,其特征在于,所述去噪声处理至少包括高斯滤波、均值滤波、傅里叶变换滤波、小波变换滤波和小波包变换滤波中的任一种方法,用于消除由于信号来自体表生理采集和分帧式的信号特征分析对睡眠状态刻画的不连续性,而在所述睡眠状态特征曲线中引入的高斯噪声或白噪声干扰。
16.如权利要求1或2所述的方法,其特征在于,所述根据所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线,计算得到睡眠可持续性指数和睡眠可持续性指数曲线的具体步骤还包括:
根据所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线,动态计算得到所述睡眠可持续性指数;
按照时序拼接所述睡眠可持续性指数,生成所述睡眠可持续性指数曲线。
17.如权利要求1或2所述的方法,其特征在于,所述对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析,生成睡眠可持续性动态调节策略并完成用户睡眠过程的动态调节及效果评估,动态迭代优化睡眠可持续性的动态检测调节过程方法策略的具体步骤还包括:
对所述睡眠状态特征曲线进行动态地趋势预测分析,得到睡眠状态特征预测值;
对所述睡眠可持续性指数曲线进行动态地趋势预测分析,得到睡眠可持续性指数预测值;
根据所述睡眠状态特征预测值和所述睡眠可持续性指数预测值,按照睡眠可持续性检测调节周期和预设睡眠调节知识库,生成所述睡眠可持续性动态调节策略;
根据所述睡眠可持续性动态调节策略,连接并控制睡眠调节外围设备,对用户睡眠过程进行动态调节;
对调节效果进行动态跟踪评估,计算动态调节效果系数并生成动态调节效果曲线;
根据所述动态调节效果曲线和所述动态调节效果系数,动态迭代优化睡眠可持续性的动态检测调节过程方法策略。
18.如权利要求17所述的方法,其特征在于,所述趋势预测分析的方法至少包括AR、MA、ARMA、ARIMA、SARIMA、VAR、深度学习中的任一项。
19.如权利要求17所述的方法,其特征在于,所述睡眠可持续性动态调节策略至少包括睡眠场景、睡眠时相、调节方式、执行方式、调节强度、调节时点、持续时间、目标调节值和装置控制参数;其中,所述调节方式至少包括声乐刺激、超声刺激、光刺激、电刺激、磁刺激、温度刺激、湿度刺激、触觉刺激和浓度调控中的任一方式,所述执行方式至少包括离体式和接触式中的任一方式。
20.如权利要求17所述的方法,其特征在于:所述睡眠调节外围设备至少包括声乐刺激设备、超声刺激设备、光刺激设备、电刺激设备、磁刺激设备、温度刺激设备、湿度刺激设备、触觉刺激设备和浓度调控设备中的任一项,并由具体的调节方式决定。
21.如权利要求18-20任一项所述的方法,其特征在于:所述动态调节效果系数的一种计算方式,具体为:
;
其中,为所述动态调节效果系数,/>为趋势预测分析方法修正系数且,/>为动态调节后的所述睡眠时相分期修正系数,/>分别为动态调节前、动态调节后的所述睡眠可持续性指数。
22.如权利要求21所述的方法,其特征在于:所述动态调节效果曲线由按照时序拼接所述动态调节效果系数而生成得到。
23.如权利要求18-20任一项所述的方法,其特征在于:所述动态检测调节过程方法策略至少包括去高斯噪声处理方法参数、时频分析方法参数、所述特征曲线高频阈值、趋势预测分析方法参数、所述睡眠可持续性动态调节策略和动态调节效果评估方法参数。
24.如权利要求1或2所述的方法,其特征在于,所述生成用户睡眠可持续性检测调节报告的具体步骤还包括:
根据所述睡眠时相曲线和所述睡眠可持续性指数曲线,计算不同睡眠时相下所述睡眠可持续性指数的分布特征,得到时相指数分布特征;
根据所述睡眠时相曲线和动态调节效果曲线,计算不同睡眠时相下动态调节效果系数的分布特征,得到时相效果分布特征;
计算所述睡眠可持续性指数曲线和所述动态调节效果曲线的相关性,得到动态调节综合效果系数;
按照预设报告生成周期,生成所述用户睡眠可持续性检测调节报告。
25.如权利要求24所述的方法,其特征在于,所述分布特征至少包括平均值、均方根、最大值、最小值、方差、标准差、变异系数、峰度和偏度中的任意一项。
26.如权利要求24所述的方法,其特征在于,所述相关性的计算方法至少包括相干性分析、皮尔逊相关分析、杰卡德相似分析、线性互信息分析、线性相关分析、欧氏距离分析、曼哈顿距离分析、切比雪夫距离分析中的任意一项。
27.如权利要求24所述的方法,其特征在于,所述用户睡眠可持续性检测调节报告至少包括用户关键生理信息、所述睡眠状态特征曲线、所述睡眠时相曲线、所述睡眠可持续性指数曲线、所述动态调节效果曲线、所述时相指数分布特征、所述时相效果分布特征、所述动态调节综合效果系数、睡眠可持续性检测调节总结。
28.如权利要求1或2所述的方法,其特征在于:还包括初始化建立或持续更新用户个性睡眠可持续性数据库,更新睡眠可持续性通用数据库的具体步骤:
动态收集检测调节过程数据和分析结果,初始化建立或持续动态更新所述用户个性睡眠可持续性数据库;
按照预设通用数据库更新周期,根据所述用户个性睡眠可持续性数据库和所述用户睡眠可持续性检测调节报告,将更新所述睡眠可持续性通用数据库的用户关键数据信息。
29.一种基于时频分析的睡眠可持续性检测调节系统,其特征在于,包括以下模块:
状态检测分析模块,用于初始化睡眠可持续性动态检测调节基本方案,对用户睡眠生理信号进行采集和处理,得到睡眠状态特征曲线和睡眠时相曲线;
特征曲线提取模块,用于根据所述睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值并进行时频分析,得到高频时频分析信号并进行去噪声处理,得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线;
指数曲线提取模块,用于根据所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线,计算得到睡眠可持续性指数和睡眠可持续性指数曲线;
动态调节优化模块,用于对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析,生成睡眠可持续性动态调节策略并完成用户睡眠过程的动态调节及效果评估,动态迭代优化睡眠可持续性的动态检测调节过程方法策略;
统计分析报告模块,用于生成用户睡眠可持续性检测调节报告;
数据运行管理模块,用于对所述系统的所有数据进行可视化管理、统一存储和运行管理;
所述睡眠可持续性指数的计算方法,具体为:
1)获取当前时帧对应的睡眠时相分期、所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线;
2)计算所述睡眠状态特征曲线和所述睡眠可持续性时频噪声曲线的差值,得到目标睡眠状态特征曲线;
3)分别对所述目标睡眠状态特征曲线和所述睡眠可持续性时频特征曲线先后做取绝对值和求均值运算,得到目标睡眠状态特征均值和睡眠可持续性TF特征均值;
4)根据所述目标睡眠状态特征均值和所述睡眠可持续性TF特征均值的相对变化量,计算得到睡眠可持续性节点特征系数;
5)根据时频分析方法修正系数、用户个性修正系数和睡眠时相分期修正系数的均值对所述睡眠可持续性节点特征系数进行修正,计算得到所述睡眠可持续性指数;
所述睡眠可持续性节点特征系数的一种计算公式,具体为:
;
其中,为所述睡眠可持续性节点特征系数,/>分别为所述目标睡眠状态特征均值和所述睡眠可持续性TF特征均值;
所述睡眠可持续性指数的一种计算公式,具体为:
;
其中,为所述睡眠可持续性指数,/>分别为所述睡眠可持续性节点特征系数、所述时频分析方法修正系数、所述用户个性修正系数和所述睡眠时相分期修正系数,且/>,/>2。
30.如权利要求29所述的系统,其特征在于,还包括:数据库更新模块,用于初始化建立或持续更新用户个性睡眠可持续性数据库,更新睡眠可持续性通用数据库。
31.如权利要求30所述的系统,其特征在于,所述状态检测分析模块还包括以下功能单元:
方案初始化单元,用于获取用户关键生理信息,从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选并初始化所述睡眠可持续性动态检测调节基本方案;
信号采集处理单元,用于对用户睡眠生理信号进行动态地监测采集和信号处理,得到睡眠生理状态数据;
信号特征分析单元,用于对所述睡眠生理状态数据进行动态地特征分析和特征选择,生成所述睡眠状态特征曲线;
睡眠时相识别单元,用于对所述睡眠生理状态数据进行动态地睡眠时相分析,生成所述睡眠时相曲线。
32.如权利要求29所述的系统,其特征在于,所述特征曲线提取模块还包括以下功能单元:
信号阈值确定单元,用于根据所述睡眠状态特征曲线的曲线生成属性,确定所述特征曲线高频阈值;
高频分量提取单元,用于对所述睡眠状态特征曲线进行动态地时频分析,根据所述特征曲线高频阈值,提取得到所述高频时频分析信号;
特征曲线生成单元,用于对所述高频时频分析信号进行动态地去噪声处理,分离得到所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线。
33.如权利要求32所述的系统,其特征在于,所述指数曲线提取模块还包括以下功能单元:
特征指数计算单元,用于根据所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线,动态计算得到所述睡眠可持续性指数;
指数曲线生成单元,用于按照时序拼接所述睡眠可持续性指数,生成所述睡眠可持续性指数曲线。
34.如权利要求29-33任一项所述的系统,其特征在于,所述动态调节优化模块还包括以下功能单元:
状态特征预测单元,用于对所述睡眠状态特征曲线进行动态地趋势预测分析,得到睡眠状态特征预测值;
特征指数预测单元,用于对所述睡眠可持续性指数曲线进行动态地趋势预测分析,得到睡眠可持续性指数预测值;
动态策略生成单元,用于根据所述睡眠状态特征预测值和所述睡眠可持续性指数预测值,按照睡眠可持续性检测调节周期和预设睡眠调节知识库,生成所述睡眠可持续性动态调节策略;
动态调节执行单元,用于根据所述睡眠可持续性动态调节策略,连接并控制睡眠调节外围设备,对用户睡眠过程进行动态调节;
调节效果评估单元,用于对调节效果进行动态跟踪评估,计算动态调节效果系数并生成动态调节效果曲线;
过程方法优化单元,用于根据所述动态调节效果曲线和所述动态调节效果系数,动态迭代优化睡眠可持续性的动态检测调节过程方法策略。
35.如权利要求34所述的系统,其特征在于,所述统计分析报告模块还包括以下功能单元:
指数分布分析单元,用于根据所述睡眠时相曲线和所述睡眠可持续性指数曲线,计算不同睡眠时相下所述睡眠可持续性指数的分布特征,得到时相指数分布特征;
效果分布分析单元,用于根据所述睡眠时相曲线和所述动态调节效果曲线,计算不同睡眠时相下所述动态调节效果系数的分布特征,得到时相效果分布特征;
综合效果评估单元,用于计算所述睡眠可持续性指数曲线和所述动态调节效果曲线的相关性,得到动态调节综合效果系数;
用户报告生成单元,用于按照预设报告生成周期,生成所述用户睡眠可持续性检测调节报告;
用户报告管理单元,用于对所述用户睡眠可持续性检测调节报告的格式输出、展现形式进行统一管理。
36.如权利要求31所述的系统,其特征在于,所述数据库更新模块还包括以下功能单元:
个性数据库更新单元,用于动态收集检测调节过程数据和分析结果,初始化建立或持续动态更新所述用户个性睡眠可持续性数据库;
通用数据库更新单元,用于按照预设通用数据库更新周期,根据所述用户个性睡眠可持续性数据库和所述用户睡眠可持续性检测调节报告,将更新所述睡眠可持续性通用数据库的用户关键数据信息。
37.如权利要求29-33任一项所述的系统,其特征在于,还包括:所述数据运行管理模块还包括以下功能单元:
用户信息管理单元,用于用户基本信息的登记输入、编辑、查询、输出和删除;
数据可视化管理单元,用于对所述系统中所有数据的可视化展现管理;
数据存储管理单元,用于对所述系统中所有数据的统一存储管理;
数据运营管理单元,用于对所述系统中所有数据的备份、迁移和导出。
38.一种基于时频分析的睡眠可持续性检测调节装置,其特征在于,包括以下模组:
状态检测分析模组,用于初始化睡眠可持续性动态检测调节基本方案,对用户睡眠生理信号进行采集和处理,得到睡眠状态特征曲线和睡眠时相曲线;
特征曲线提取模组,用于根据所述睡眠状态特征曲线的曲线生成属性,确定特征曲线高频阈值并进行时频分析,得到高频时频分析信号并进行去噪声处理,得到睡眠可持续性时频噪声曲线和睡眠可持续性时频特征曲线;
指数曲线提取模组,用于根据所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线,计算得到睡眠可持续性指数和睡眠可持续性指数曲线;
动态调节优化模组,用于对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析,生成睡眠可持续性动态调节策略并完成用户睡眠过程的动态调节及效果评估,动态迭代优化睡眠可持续性的动态检测调节过程方法策略;
统计分析报告模组,用于生成用户睡眠可持续性检测调节报告;
数据库更新模组,用于初始化建立或持续更新用户个性睡眠可持续性数据库,更新睡眠可持续性通用数据库;
数据可视化模组,用于对所述装置中所有过程数据和结果数据的统一可视化展示管理;
数据管理中心模组,用于对所述装置中所有过程数据和结果数据的统一存储和数据运营管理;
所述睡眠可持续性指数的计算方法,具体为:
1)获取当前时帧对应的睡眠时相分期、所述睡眠状态特征曲线、所述睡眠可持续性时频噪声曲线和所述睡眠可持续性时频特征曲线;
2)计算所述睡眠状态特征曲线和所述睡眠可持续性时频噪声曲线的差值,得到目标睡眠状态特征曲线;
3)分别对所述目标睡眠状态特征曲线和所述睡眠可持续性时频特征曲线先后做取绝对值和求均值运算,得到目标睡眠状态特征均值和睡眠可持续性TF特征均值;
4)根据所述目标睡眠状态特征均值和所述睡眠可持续性TF特征均值的相对变化量,计算得到睡眠可持续性节点特征系数;
5)根据时频分析方法修正系数、用户个性修正系数和睡眠时相分期修正系数的均值对所述睡眠可持续性节点特征系数进行修正,计算得到所述睡眠可持续性指数;
所述睡眠可持续性节点特征系数的一种计算公式,具体为:
;
其中,为所述睡眠可持续性节点特征系数,/>分别为所述目标睡眠状态特征均值和所述睡眠可持续性TF特征均值;
所述睡眠可持续性指数的一种计算公式,具体为:
;
其中,为所述睡眠可持续性指数,/>分别为所述睡眠可持续性节点特征系数、所述时频分析方法修正系数、所述用户个性修正系数和所述睡眠时相分期修正系数,且/>,/>2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310974382.7A CN116687356B (zh) | 2023-08-04 | 2023-08-04 | 基于时频分析的睡眠可持续性检测调节方法、系统和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310974382.7A CN116687356B (zh) | 2023-08-04 | 2023-08-04 | 基于时频分析的睡眠可持续性检测调节方法、系统和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116687356A CN116687356A (zh) | 2023-09-05 |
CN116687356B true CN116687356B (zh) | 2024-05-07 |
Family
ID=87843633
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310974382.7A Active CN116687356B (zh) | 2023-08-04 | 2023-08-04 | 基于时频分析的睡眠可持续性检测调节方法、系统和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116687356B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102274022A (zh) * | 2011-05-10 | 2011-12-14 | 浙江大学 | 一种基于脑电信号的睡眠状态监测方法 |
CN105942974A (zh) * | 2016-04-14 | 2016-09-21 | 禅客科技(上海)有限公司 | 一种基于低频脑电的睡眠分析方法及系统 |
GB201803102D0 (en) * | 2018-02-26 | 2018-04-11 | Sentimoto Ltd | Sleep activity detection method and apparatus |
CN115862877A (zh) * | 2023-03-03 | 2023-03-28 | 安徽星辰智跃科技有限责任公司 | 睡眠可持续性检测量化及辅助干预的方法、系统和装置 |
CN116392088A (zh) * | 2023-06-08 | 2023-07-07 | 安徽星辰智跃科技有限责任公司 | 基于时频分析的睡眠趋稳性量化及调节方法、系统和装置 |
CN116392085A (zh) * | 2023-06-06 | 2023-07-07 | 安徽星辰智跃科技有限责任公司 | 基于趋势分析的睡眠趋稳性量化及调节方法、系统和装置 |
CN116509336A (zh) * | 2023-06-27 | 2023-08-01 | 安徽星辰智跃科技有限责任公司 | 基于波形分析的睡眠周期性检测及调节方法、系统和装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100809041B1 (ko) * | 2006-06-20 | 2008-03-03 | 삼성전자주식회사 | 수면 상태 감지 장치 및 그 방법 |
US9655559B2 (en) * | 2014-01-03 | 2017-05-23 | Vital Connect, Inc. | Automated sleep staging using wearable sensors |
-
2023
- 2023-08-04 CN CN202310974382.7A patent/CN116687356B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102274022A (zh) * | 2011-05-10 | 2011-12-14 | 浙江大学 | 一种基于脑电信号的睡眠状态监测方法 |
CN105942974A (zh) * | 2016-04-14 | 2016-09-21 | 禅客科技(上海)有限公司 | 一种基于低频脑电的睡眠分析方法及系统 |
GB201803102D0 (en) * | 2018-02-26 | 2018-04-11 | Sentimoto Ltd | Sleep activity detection method and apparatus |
CN115862877A (zh) * | 2023-03-03 | 2023-03-28 | 安徽星辰智跃科技有限责任公司 | 睡眠可持续性检测量化及辅助干预的方法、系统和装置 |
CN116392085A (zh) * | 2023-06-06 | 2023-07-07 | 安徽星辰智跃科技有限责任公司 | 基于趋势分析的睡眠趋稳性量化及调节方法、系统和装置 |
CN116392088A (zh) * | 2023-06-08 | 2023-07-07 | 安徽星辰智跃科技有限责任公司 | 基于时频分析的睡眠趋稳性量化及调节方法、系统和装置 |
CN116509336A (zh) * | 2023-06-27 | 2023-08-01 | 安徽星辰智跃科技有限责任公司 | 基于波形分析的睡眠周期性检测及调节方法、系统和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN116687356A (zh) | 2023-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116392085B (zh) | 基于趋势分析的睡眠趋稳性量化及调节方法、系统和装置 | |
CN116369866B (zh) | 基于小波变换的睡眠趋稳性量化及调节方法、系统和装置 | |
CN116392083B (zh) | 基于局部变换的睡眠趋稳性量化及调节方法、系统和装置 | |
CN115910351B (zh) | 一种睡眠周期性检测量化及辅助干预的方法、系统和装置 | |
CN115804573A (zh) | 一种睡眠深度量化及干预的方法、系统和装置 | |
CN115862877B (zh) | 睡眠可持续性检测量化及辅助干预的方法、系统和装置 | |
CN116509338B (zh) | 基于模态分析的睡眠周期性检测及调节方法、系统和装置 | |
CN116058804B (zh) | 一种睡眠情绪活跃水平动态调节的方法、系统和装置 | |
CN116058805B (zh) | 一种睡眠记忆活跃水平动态调节的方法、系统和装置 | |
CN116705247B (zh) | 基于局部分解的睡眠可持续性检测调节方法、系统和装置 | |
CN116392088B (zh) | 基于时频分析的睡眠趋稳性量化及调节方法、系统和装置 | |
CN116525063B (zh) | 基于时频分析的睡眠周期性检测及调节方法、系统和装置 | |
CN116509336B (zh) | 基于波形分析的睡眠周期性检测及调节方法、系统和装置 | |
CN116504357B (zh) | 基于小波分析的睡眠周期性检测及调节方法、系统和装置 | |
CN116392087B (zh) | 基于模态分解的睡眠趋稳性量化及调节方法、系统和装置 | |
CN115969330B (zh) | 一种睡眠情绪活跃水平检测量化的方法、系统和装置 | |
CN116687356B (zh) | 基于时频分析的睡眠可持续性检测调节方法、系统和装置 | |
CN116671875B (zh) | 基于小波变换的睡眠可持续性检测调节方法、系统和装置 | |
CN116682535B (zh) | 基于数值拟合的睡眠可持续性检测调节方法、系统和装置 | |
CN116035536B (zh) | 一种睡眠行为活跃水平检测量化的方法、系统和装置 | |
CN116779110B (zh) | 基于模态分解的睡眠可持续性检测调节方法、系统和装置 | |
CN116509337B (zh) | 基于局部分解的睡眠周期性检测及调节方法、系统和装置 | |
CN116509337A (zh) | 基于局部分解的睡眠周期性检测及调节方法、系统和装置 | |
CN116312928B (zh) | 一种眼疲劳恢复进程的辅助调节方法、系统和装置 | |
CN116312972B (zh) | 基于眼刺激的睡眠记忆情感张力调节方法、系统和装置 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |