CN116392085B - 基于趋势分析的睡眠趋稳性量化及调节方法、系统和装置 - Google Patents
基于趋势分析的睡眠趋稳性量化及调节方法、系统和装置 Download PDFInfo
- Publication number
- CN116392085B CN116392085B CN202310658830.2A CN202310658830A CN116392085B CN 116392085 B CN116392085 B CN 116392085B CN 202310658830 A CN202310658830 A CN 202310658830A CN 116392085 B CN116392085 B CN 116392085B
- Authority
- CN
- China
- Prior art keywords
- sleep
- trend
- curve
- user
- characteristic curve
- 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 638
- 238000004458 analytical method Methods 0.000 title claims abstract description 162
- 238000000034 method Methods 0.000 title claims abstract description 157
- 238000011002 quantification Methods 0.000 title claims abstract description 70
- 230000008569 process Effects 0.000 claims abstract description 72
- 230000000694 effects Effects 0.000 claims abstract description 67
- 230000035790 physiological processes and functions Effects 0.000 claims abstract description 61
- 230000008667 sleep stage Effects 0.000 claims abstract description 59
- 238000009826 distribution Methods 0.000 claims abstract description 52
- 238000012544 monitoring process Methods 0.000 claims abstract description 17
- 230000006399 behavior Effects 0.000 claims abstract description 14
- 238000011156 evaluation Methods 0.000 claims abstract description 14
- 230000000638 stimulation Effects 0.000 claims description 54
- 238000001514 detection method Methods 0.000 claims description 42
- 238000004364 calculation method Methods 0.000 claims description 38
- 238000012937 correction Methods 0.000 claims description 35
- 238000007726 management method Methods 0.000 claims description 33
- 238000012545 processing Methods 0.000 claims description 29
- 230000001105 regulatory effect Effects 0.000 claims description 26
- 238000013139 quantization Methods 0.000 claims description 20
- 238000000605 extraction Methods 0.000 claims description 16
- 238000001914 filtration Methods 0.000 claims description 15
- 230000004622 sleep time Effects 0.000 claims description 14
- 230000004927 fusion Effects 0.000 claims description 12
- 230000000007 visual effect Effects 0.000 claims description 12
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 9
- 229910052760 oxygen Inorganic materials 0.000 claims description 9
- 239000001301 oxygen Substances 0.000 claims description 9
- 230000002093 peripheral effect Effects 0.000 claims description 9
- 230000011218 segmentation Effects 0.000 claims description 9
- 230000001755 vocal effect Effects 0.000 claims description 8
- 238000013135 deep learning Methods 0.000 claims description 7
- 210000000467 autonomic pathway Anatomy 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 6
- 230000004617 sleep duration Effects 0.000 claims description 6
- 241001123248 Arma Species 0.000 claims description 5
- 238000010219 correlation analysis Methods 0.000 claims description 5
- 230000009467 reduction Effects 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 241000728173 Sarima Species 0.000 claims description 4
- 230000008030 elimination Effects 0.000 claims description 4
- 238000003379 elimination reaction Methods 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 4
- YHXISWVBGDMDLQ-UHFFFAOYSA-N moclobemide Chemical compound C1=CC(Cl)=CC=C1C(=O)NCCN1CCOCC1 YHXISWVBGDMDLQ-UHFFFAOYSA-N 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 4
- 238000010220 Pearson correlation analysis Methods 0.000 claims description 3
- 238000012952 Resampling Methods 0.000 claims description 3
- 239000008280 blood Substances 0.000 claims description 3
- 210000004369 blood Anatomy 0.000 claims description 3
- 230000036760 body temperature Effects 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000013523 data management 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
- 201000010099 disease Diseases 0.000 claims description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 3
- 230000003862 health status Effects 0.000 claims description 3
- 238000002582 magnetoencephalography Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000012549 training Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 238000002604 ultrasonography Methods 0.000 claims description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 claims 1
- 230000001568 sexual effect Effects 0.000 claims 1
- 230000033764 rhythmic process Effects 0.000 description 9
- 230000000875 corresponding effect Effects 0.000 description 8
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000000241 respiratory effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000002567 autonomic effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000001054 cortical effect Effects 0.000 description 1
- 238000013136 deep learning model Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 230000008452 non REM sleep Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000036385 rapid eye movement (rem) sleep Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 208000019116 sleep disease Diseases 0.000 description 1
- 230000008665 sleep physiology Effects 0.000 description 1
- 230000009466 transformation Effects 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
- A61B5/48—Other medical applications
- A61B5/4806—Sleep evaluation
- A61B5/4812—Detecting sleep stages or cycles
-
- 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
- A61B5/02055—Simultaneously evaluating both cardiovascular condition and temperature
-
- 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
- A61B5/14542—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 for measuring blood gases
-
- 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]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
本发明提供了基于趋势分析的睡眠趋稳性量化及调节方法,包括:对用户睡眠过程生理状态进行采集监测处理、特征分析和时相识别,得到睡眠深度特征曲线和睡眠分期曲线;根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,反向提取趋势分量,提取睡眠趋势特征曲线;根据所述睡眠分期曲线、所述睡眠深度特征曲线和所述睡眠趋势特征曲线,动态计算提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征;建立并动态更新用户个性化睡眠趋势性数据库,对用户睡眠行为进行动态预测、干预调节、效果评估,动态更新数据库并优化调节策略。本发明能够实现对用户睡眠趋稳性的高效干预调节。
Description
技术领域
本发明涉及睡眠趋稳性检测量化及辅助调节领域,特别涉及基于趋势分析的睡眠趋稳性量化及调节方法、系统和装置。
背景技术
人类的睡眠过程除了非快速眼动睡眠和快速眼动睡眠的多次交替循环属性外,还具有趋势性或趋稳性,即整体睡眠深度基线水平由深到浅,直至睡眠结束。
申请人提出的在先解决方案中国申请CN202310195993提供了睡眠趋稳性检测量化及辅助干预的方法,包括:采集用户睡眠过程的生理体征数据和环境因素数据,进行数据预处理、时帧处理和时帧特征分析,生成生理体征特征和环境因素特征;对所述生理体征特征进行睡眠状态分析、时序成分分析和趋稳性量化分析,提取睡眠趋稳性指数,生成睡眠趋稳性量化日报;重复上述步骤,对用户睡眠过程进行连续监测和跟踪分析,评估睡眠环境对睡眠趋稳性的影响,提取最佳睡眠趋稳性环境方案并对睡眠环境进行动态优化调整,生成睡眠趋稳性量化报告。该技术方案提出了睡眠趋稳性的创新评价指标,通过时间序列分解方法从睡眠持续期状态特征曲线中提取趋势成分,计算趋势强度并得到睡眠趋稳性指数,已经能够初步解决趋稳性量化的问题。但还存在可以提升的空间,主要包括:首先,时间序列分解方法在计算效率、多场景评价和人群适应度上表现出局限性;其次,睡眠持续期状态特征曲线是基于阶梯式睡眠时相分期值平滑得到的,仅能对用户睡眠状态进行简要描述,而不能细致地刻画和量化,最终了带来了睡眠趋稳性的评价不够细致和灵敏;最后,如何更优地通过睡眠环境的调节来实现对用户睡眠趋稳性的高效干预调节。
如何识别更灵敏的用户睡眠状态表征特征;如何进一步更准确快速地提取睡眠趋稳性特征,得到更精确的、不同场景不同人群的睡眠趋稳性评价;如何实现更精准高效的、多手段方式的、实时动态的用户睡眠趋稳性干预调节;如何更进一步整体提高用户个性化的检测量化效率和干预调节效果等,是目前国内外产品技术方案和实际应用场景中需要进一步优化或解决的问题。
发明内容
针对现有方法的以上缺陷及改进需求,本发明的目的在于提供一种基于趋势分析的睡眠趋稳性量化及调节方法,通过对用户睡眠过程生理状态进行采集监测处理、特征分析和时相识别,得到睡眠深度特征曲线并进行去趋势分析,提取睡眠趋势特征曲线,进一步而提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征,完成用户睡眠趋势性检测量化;对用户睡眠行为趋势进行预测,提取个性时相场景睡眠调节策略,实现多手段方式的、实时动态的用户睡眠趋稳性干预调节;将检测量化和干预调节的关键过程数据纳入数据库,建立并持续更新用户个性化睡眠趋势性数据库,持续提高用户个性化的检测量化效率和干预调节效果。本发明还提供了基于趋势分析的睡眠趋稳性量化及调节系统,用于实现上述方法。本发明还提供了基于趋势分析的睡眠趋稳性量化及调节装置,用于实现上述系统。
根据本发明的目的,本发明提出了一种基于趋势分析的睡眠趋稳性量化及调节方法,包括以下步骤:
对用户睡眠过程生理状态进行采集监测处理、特征分析和时相识别,得到睡眠深度特征曲线和睡眠分期曲线;
根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,反向提取趋势分量,提取睡眠趋势特征曲线;
根据所述睡眠分期曲线、所述睡眠深度特征曲线和所述睡眠趋势特征曲线,动态计算提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征;
建立并动态更新用户个性化睡眠趋势性数据库;
根据所述用户个性化睡眠趋势性数据库,对用户睡眠行为进行动态预测,提取个性时相场景睡眠调节策略并对用户睡眠过程进行干预调节;
进一步动态更新所述用户个性化睡眠趋势性数据库并优化所述个性时相场景睡眠调节策略,按照预设报告周期生成睡眠趋势性量化及调节报告。
更优地,所述对用户睡眠过程生理状态进行采集监测处理、特征分析和时相识别,得到睡眠深度特征曲线和睡眠分期曲线的具体步骤还包括:
对用户睡眠过程生理状态进行采集监测和信号处理,得到用户睡眠生理状态时帧数据;
对所述用户睡眠生理状态时帧数据进行特征分析和特征融合,得到所述睡眠深度特征曲线;
对所述用户睡眠生理状态时帧数据进行时相识别和睡眠分期,得到所述睡眠分期曲线。
更优地,所述信号处理至少包括AD数模转换、重采样、重参考、降噪、去伪迹、工频陷波、低通滤波、高通滤波、带阻滤波、带通滤波、矫正处理和时帧分割;其中,所述矫正处理具体为对信号中的包含伪迹或失真的信号数据片段进行信号矫正、信号预测和平滑处理,所述时帧分割是指根据信号采样率,以预设时间窗口长度和预设时间平移步长对目标信号数据进行连续滑动分割。
更优地,所述用户睡眠生理状态时帧数据至少包括脑中枢状态数据、自主神经状态数据中的任一项;其中,所述脑中枢状态数据至少包括脑电图信号、脑磁图信号、血氧水平依赖信号中的任一项,所述自主神经状态数据至少包括中血氧水平依赖信号、心电信号、脉搏信号、呼吸信号、血氧信号、体温信号和皮肤电信号的任一项。
更优地,所述特征分析至少包括数值分析、包络分析、时频分析、熵分析、分形分析和复杂度分析。
更优地,所述特征融合是指从所述特征分析得到的目标特征集中,选择预设特征数量的目标特征并进行加权计算,生成所述睡眠深度特征曲线。
更优地,所述睡眠深度特征曲线具体为表征用户在预设入睡前时期、睡眠持续期、预设结束睡眠后时期的睡眠深度和时相状态的特征曲线,其计算生成方法具体为:
1)对所述用户睡眠生理状态时帧数据中时帧数据逐个进行特征分析,按时序拼接得到睡眠生理状态时帧特征曲线集;
2)从所述睡眠生理状态时帧特征曲线集中,筛选目标相关特征曲线,并按照预设特征融合权重进行加权计算,生成所述睡眠深度特征曲线。
更优地,所述睡眠分期曲线的生成方法具体为:
1)通过深度学习算法对规模睡眠用户样本的所述用户睡眠生理状态时帧数据及其对应睡眠分期数据进行学习训练和数据建模,得到睡眠时相自动分期模型;
2)将当前用户的所述用户睡眠生理状态时帧数据输入所述睡眠时相自动分期模型,得到所对应的睡眠时相分期并按时序生成所述睡眠分期曲线。
更优地,所述根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,反向提取趋势分量,提取睡眠趋势特征曲线的具体步骤还包括:
根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,得到FA波动分量;
根据所述睡眠深度特征曲线和所述FA波动分量,反向提取趋势分量,生成所述睡眠趋势特征曲线。
更优地,所述窗口尺度的选取取决于所述睡眠深度特征曲线的特征来源和组合、动态调节效果方式。
更优地,所述拟合方式至少包括最小二乘法拟合、多项式拟合、线性拟合中的任意一项,并由所述睡眠深度特征曲线的特征来源和组合方式共同决定。
更优地,所述去趋势分析的方法至少包括去趋势分析DFA、多重分形去趋势分析MFDFA、非对称消除趋势波动分析ADFA中的任意一项。
更优地,所述睡眠趋势特征曲线的计算公式具体为:
;
其中,分别为所述睡眠趋势特征曲线、所述睡眠深度特征曲线和所述FA波动分量。
更优地,所述根据所述睡眠分期曲线、所述睡眠深度特征曲线和所述睡眠趋势特征曲线,动态计算提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征的具体步骤还包括:
根据所述睡眠深度特征曲线和所述睡眠趋势特征曲线,计算得到所述睡眠趋势性指数;
根据所述睡眠分期曲线和所述睡眠趋势特征曲线进行相关性计算,得到所述时相趋势相关系数;
基于所述睡眠分期曲线,对所述睡眠趋势特征曲线进行时相分布统计,得到所述时相趋势分布特征。
更优地,所述睡眠趋势性指数的计算方法具体为:
1)获取所述睡眠深度特征曲线和所述睡眠趋势特征曲线;
2)分别对所述睡眠深度特征曲线和所述睡眠趋势特征曲线进行取平方,得到睡眠深度特征平方曲线和睡眠趋势特征平方曲线;
3)计算所述睡眠深度特征平方曲线和所述睡眠趋势特征平方曲线的样本点比值,得到睡眠趋势性曲线;
4)求取所述睡眠趋势性曲线的平均值,得到睡眠趋势性特征值;
5)计算所述睡眠趋势性特征值、去趋势分析方法对应的预设方法修正系数、与用户生物状态信息相关的预设用户个性修正系数的乘积,生成所述睡眠趋势性指数。
更优地,所述用户生物状态信息至少包括性别、年龄、职业、健康状态、疾病状态和教育程度。
更优地,所述睡眠趋势性指数的计算公式具体为:
;
其中,为所述睡眠趋势性指数,/>分别为预设用户个性修正系数和预设方法修正系数,/>、/>分别为所述睡眠趋势特征曲线和所述睡眠深度特征曲线第i点的值,/>为所述睡眠趋势特征曲线的数据长度。
更优地,所述相关性计算的方法至少包括相干性分析、皮尔逊相关分析、杰卡德相似分析、线性互信息分析、线性相关分析、欧氏距离分析、曼哈顿距离分析、切比雪夫距离分析中的任意一项。
更优地,所述时相分布统计具体为按照所述睡眠分期曲线中的睡眠时相分期,对所述睡眠趋势特征曲线中的睡眠趋势特征值进行数值分布统计分析,得到所述睡眠趋势特征曲线的数值分布统计特征。
更优地,所述时相趋势分布特征至少包括平均值、均方根、最大值、最小值、方差、标准差、变异系数、峰度和偏度中的任意一项。
更优地,所述建立并动态更新用户个性化睡眠趋势性数据库的具体步骤还包括:
建立所述用户个性化睡眠趋势性数据库;
将用户睡眠趋稳性动态检测量化的关键过程数据,动态更新入所述用户个性化睡眠趋势性数据库。
更优地,所述用户个性化睡眠趋势性数据库至少包括所述用户生物状态信息、睡眠场景、所述睡眠分期曲线、所述睡眠深度特征曲线、所述睡眠趋势特征曲线、所述睡眠趋势性指数、所述时相趋势相关系数、所述时相趋势分布特征、调节方式、执行方式、调节方法、目标调节值曲线、调节强度曲线、动态调节效果曲线、综合调节效果指数。
更优地,所述根据所述用户个性化睡眠趋势性数据库,对用户睡眠行为进行动态预测,提取个性时相场景睡眠调节策略并对用户睡眠过程进行干预调节的具体步骤还包括:
根据所述用户个性化睡眠趋势性数据库,对所述睡眠趋势特征曲线进行动态趋势预测,提取睡眠趋势预测特征值,生成睡眠趋势预测特征曲线;
根据所述用户个性化睡眠趋势性数据库、预设睡眠趋稳性调节知识库、所述睡眠趋势预测特征值和当前具体睡眠场景,按照预设动态调节周期,动态生成所述个性时相场景睡眠调节策略;
根据所述个性时相场景睡眠调节策略,连接并调控睡眠趋势性调节外围设备,对用户睡眠过程进行动态干预调节;
对动态干预调节的效果进行动态跟踪评估,提取动态调节效果系数并生成动态调节效果曲线,计算得到综合调节效果指数。
更优地,所述睡眠趋势预测特征值的预测方法至少包括AR、MR、ARMA、ARIMA、SARIMA、VAR、深度学习中的任一项。
更优地,所述个性时相场景睡眠调节策略至少包括睡眠场景、睡眠时相、调节方式、执行方式、调节方法、调节强度、调节时点、持续时间、目标调节值和装置控制参数;其中,所述调节方式至少包括声乐刺激、超声刺激、光刺激、电刺激、磁刺激、温度刺激、湿度刺激、触觉刺激和调控中的任一方式,所述执行方式至少包括离体式和接触式中的任一方式。
更优地,所述调控睡眠趋势性调节外围设备至少包括声乐刺激设备、超声刺激设备、光刺激设备、电刺激设备、磁刺激设备、温度刺激设备、湿度刺激设备、触觉刺激设备和调控设备中的任一项,并由具体的所述调节方式决定。
更优地,所述动态调节效果系数的具体计算公式具体为:
;
其中,为所述动态调节效果系数;/>分别为所述个性时相场景睡眠调节策略中的目标调节值、所述睡眠趋势预测特征值、动态调节后量化得到的所述睡眠趋势特征曲线中的睡眠趋势特征值;/>为预设用户个性修正系数。
更优地,所述动态调节效果系数反向应用于去趋势分析方法参数中所述窗口尺度的动态调整、所述个性时相场景睡眠调节策略的生成。
更优地,所述综合调节效果指数具体为所述动态调节效果曲线的平均值或均方根。
更优地,所述进一步动态更新所述用户个性化睡眠趋势性数据库并优化所述个性时相场景睡眠调节策略,按照预设报告周期生成睡眠趋势性量化及调节报告的具体步骤还包括:
对用户不同睡眠场景下的睡眠趋势性进行动态检测量化和动态干预调节,持续动态更新所述用户个性化睡眠趋势性数据库;
将动态检测量化和动态干预调节的关键过程数据,动态更新入所述用户个性化睡眠趋势性数据库,动态优化所述个性时相场景睡眠调节策略;
按照预设报告周期,生成所述睡眠趋势性量化及调节报告。
更优地,所述睡眠趋势性量化及调节报告至少包括所述用户生物状态信息、睡眠场景、所述睡眠分期曲线、所述睡眠深度特征曲线、所述睡眠趋势特征曲线、所述睡眠趋势性指数、所述时相趋势相关系数、所述时相趋势分布特征、所述综合调节效果指数、所述动态调节效果曲线、用户睡眠趋势性量化及调节总结。
根据本发明的目的,本发明提出了一种基于趋势分析的睡眠趋稳性量化及调节系统,包括以下模块:
状态检测分析模块,用于对用户睡眠过程生理状态进行采集监测处理、特征分析和时相识别,得到睡眠深度特征曲线和睡眠分期曲线;
趋势特征提取模块,用于根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,反向提取趋势分量,提取睡眠趋势特征曲线;
趋势动态量化模块,用于根据所述睡眠分期曲线、所述睡眠深度特征曲线和所述睡眠趋势特征曲线,动态计算提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征;
趋势数据库模块,用于建立并动态更新用户个性化睡眠趋势性数据库;
趋势动态调节模块,用于根据所述用户个性化睡眠趋势性数据库,对用户睡眠行为进行动态预测,提取个性时相场景睡眠调节策略并对用户睡眠过程进行干预调节;
数据动态应用模块,用于动态更新所述用户个性化睡眠趋势性数据库并优化所述个性时相场景睡眠调节策略,按照预设报告周期生成睡眠趋势性量化及调节报告;
数据运行管理模块,用于对所述系统的全部过程数据进行可视化管理、统一存储和运行管理。
更优地,所述状态检测分析模块还包括以下功能单元:
状态检测处理单元,用于对用户睡眠过程生理状态进行采集监测和信号处理,得到用户睡眠生理状态时帧数据;
深度特征提取单元,用于对所述用户睡眠生理状态时帧数据进行特征分析和特征融合,得到所述睡眠深度特征曲线;
睡眠时相识别单元,用于对所述用户睡眠生理状态时帧数据进行时相识别和睡眠分期,得到所述睡眠分期曲线。
更优地,所述趋势特征提取模块还包括以下功能单元:
去趋势分析单元,用于根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,得到FA波动分量;
趋势分量识别单元,用于根据所述睡眠深度特征曲线和所述FA波动分量,反向提取趋势分量,生成所述睡眠趋势特征曲线。
更优地,所述趋势动态量化模块还包括以下功能单元:
趋势指数量化单元,用于根据所述睡眠深度特征曲线和所述睡眠趋势特征曲线,计算得到所述睡眠趋势性指数;
相关性计算单元,用于根据所述睡眠分期曲线和所述睡眠趋势特征曲线进行相关性计算,得到所述时相趋势相关系数;
时相分布统计单元,用于基于所述睡眠分期曲线,对所述睡眠趋势特征曲线进行时相分布统计,得到所述时相趋势分布特征。
更优地,所述趋势数据库模块还包括以下功能单元:
数据库建立单元,用于建立所述用户个性化睡眠趋势性数据库;
数据库更新单元,用于将用户睡眠趋稳性动态检测量化的关键过程数据,动态更新入所述用户个性化睡眠趋势性数据库。
更优地,所述趋势动态调节模块还包括以下功能单元:
睡眠状态预测单元,用于根据所述用户个性化睡眠趋势性数据库,对所述睡眠趋势特征曲线进行动态趋势预测,提取睡眠趋势预测特征值,生成睡眠趋势预测特征曲线;
调节策略生成单元,用于根据所述用户个性化睡眠趋势性数据库、预设睡眠趋稳性调节知识库、所述睡眠趋势预测特征值和当前具体睡眠场景,按照预设动态调节周期,动态生成所述个性时相场景睡眠调节策略;
动态调控执行单元,用于根据所述个性时相场景睡眠调节策略,连接并调控睡眠趋势性调节外围设备,对用户睡眠过程进行动态干预调节;
动态效果评估单元,用于对动态干预调节的效果进行动态跟踪评估,提取动态调节效果系数并生成动态调节效果曲线,计算得到综合调节效果指数。
更优地,所述数据动态应用模块还包括以下功能单元:
趋势数据管理单元,用于对用户不同睡眠场景下的睡眠趋势性进行动态检测量化和动态干预调节,持续动态更新所述用户个性化睡眠趋势性数据库;
趋势数据应用单元,用于将动态检测量化和动态干预调节的关键过程数据,动态更新入所述用户个性化睡眠趋势性数据库,动态优化所述个性时相场景睡眠调节策略;
用户报告生成单元,用于按照预设报告周期,生成所述睡眠趋势性量化及调节报告;
用户报告输出单元,用于对所述睡眠趋势性量化及调节报告的格式输出、展现形式进行统一管理。
更优地,所述数据运行管理模块还包括以下功能单元:
用户信息管理单元,用于用户基本信息的登记输入、编辑、查询、输出和删除;
数据可视化管理单元,用于对所述系统中所有数据的可视化展现管理;
数据存储管理单元,用于对所述系统中所有数据的统一存储管理;
数据运营管理单元,用于对所述系统中所有数据的备份、迁移和导出。
根据本发明的目的,本发明提出了一种基于趋势分析的睡眠趋稳性量化及调节装置,包括以下模组:
状态检测分析模组,用于对用户睡眠过程生理状态进行采集监测处理、特征分析和时相识别,得到睡眠深度特征曲线和睡眠分期曲线;
趋势特征提取模组,用于根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,反向提取趋势分量,提取睡眠趋势特征曲线;
趋势动态量化模组,用于根据所述睡眠分期曲线、所述睡眠深度特征曲线和所述睡眠趋势特征曲线,动态计算提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征;
趋势数据库模组,用于建立并动态更新用户个性化睡眠趋势性数据库;
趋势动态调节模组,用于根据所述用户个性化睡眠趋势性数据库,对用户睡眠行为进行动态预测,提取个性时相场景睡眠调节策略并对用户睡眠过程进行干预调节;
数据动态应用模组,用于动态更新所述用户个性化睡眠趋势性数据库并优化所述个性时相场景睡眠调节策略,按照预设报告周期生成睡眠趋势性量化及调节报告;
数据可视化模组,用于对所述装置中所有过程数据和/或结果数据的统一可视化展示管理;
数据管理中心模组,用于对所述装置中所有过程数据和/或结果数据的统一存储和数据运营管理。
本发明所提供的基于趋势分析的睡眠趋稳性量化及调节方法、系统和装置,通过对用户睡眠过程生理状态进行采集监测处理、特征分析和时相识别,得到睡眠深度特征曲线并进行去趋势分析,提取睡眠趋势特征曲线,进一步而提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征,完成用户睡眠趋势性检测量化;对用户睡眠行为趋势进行预测,提取个性时相场景睡眠调节策略,实现多手段方式的、实时动态的用户睡眠趋稳性干预调节;将检测量化和干预调节的关键过程数据纳入数据库,建立并持续更新用户个性化睡眠趋势性数据库,持续提高用户个性化的检测量化效率和干预调节效果。
本发明在申请人的在先研究基础上进一步优化了趋稳性量化的具体设计,将去趋势分析应用于趋势性信息的提取,其考量了完整睡眠时期的状态特征,更加全面、适应度广;本申请还进一步改进了趋稳性指数的计算方式,提升了评价的细粒度和灵敏度;还提供了对应的效果系数计算方案,使得对于调整过程的控制有了有力依据。本发明能够提供更加科学高效的睡眠趋稳性检测量化及干预调节的实施方法和落地方案。在实际应用场景中,本发明提供的基于模态分解的睡眠趋稳性量化及调节方法、系统和装置,能够赋能相关睡眠量化或调节的产品及服务,满足不同用户场景需求,辅助用户睡眠。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明技术方案的进一步理解,并且构成说明书的一部分,与本申请的实施例一起用于解释本发明的技术方案,并不构成对本发明技术方案的限制。
图1是本发明一实施例所提供的一种基于趋势分析的睡眠趋稳性量化及调节方法的流程步骤示意图;
图2是本发明一实施例所提供的一种基于趋势分析的睡眠趋稳性量化及调节系统的模块组成示意图;
图3是本发明一实施例所提供的一种基于趋势分析的睡眠趋稳性量化及调节装置的模组构成示意图。
具体实施方式
为了更清楚地说明本发明的目的和技术方案,下面将结合本发明申请实施例中的附图,对本发明进行进一步介绍说明。显而易见地,下面描述的实施例仅是本发明一部分实施例,而不是全部的实施例。在没有创造性劳动前提下,本领域普通技术人员基于本发明的实施例所得到的其他实施例,都应属于本发明的保护范围。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互任意组合。
申请人发现在通常情况下,人类睡眠生理过程是个非平稳时序过程,其表征曲线也是个非平稳时序数据曲线。为了提取出其中的信息,申请人将去趋势分析应用于睡眠趋势性信息的提取,进一步优化了量化和调节过程。其中,去趋势分析是一种标度指数计算方法,用于分析时间序列的长程相关性,有效地滤去序列中的不同阶数的趋势成分, 从而检测含有噪声且叠加有多项式趋势信号的长程相关, 适合非平稳时间序列的长程幂律相关分析。除了标准的去趋势分析DFA外,还有多重分形去趋势分析MFDFA、非对称消除趋势波动分析ADFA等改进优化方法,以满足不同场景需求下的时序信号波动特征分析。多重分形去趋势分析MFDFA基于分形理论,实现对时间序列的多尺度重构和分形分析,去除局部趋势对时间序列标度的影响,识别不同时间标度下时间序列的分形特征,来分析提取时间序列的波动特性。非对称消除趋势波动分析ADFA引入了两个新的标度指数hq+和hq-(分别表征时间序列中正的波动趋势和负的波动趋势),来分析时间序列的波动特性。
结合图1所示,本发明实施例提供的基于趋势分析的睡眠趋稳性量化及调节方法,包括以下步骤:
P100:对用户睡眠过程生理状态进行采集监测处理、特征分析和时相识别,得到睡眠深度特征曲线和睡眠分期曲线。
第一步、对用户睡眠过程生理状态进行采集监测和信号处理,得到用户睡眠生理状态时帧数据。
本实施例中,信号处理至少包括AD数模转换、重采样、重参考、降噪、去伪迹、工频陷波、低通滤波、高通滤波、带阻滤波、带通滤波、矫正处理和时帧分割;其中,矫正处理具体为对信号中的包含伪迹或失真的信号数据片段进行信号矫正、信号预测和平滑处理,时帧分割是指根据信号采样率,以预设时间窗口长度和预设时间平移步长对目标信号数据进行连续滑动分割。
本实施例中,用户睡眠生理状态时帧数据至少包括脑中枢状态数据、自主神经状态数据中的任一项;其中,脑中枢状态数据至少包括脑电图信号、脑磁图信号、血氧水平依赖信号中的任一项,自主神经状态数据至少包括中血氧水平依赖信号、心电信号、脉搏信号、呼吸信号、血氧信号、体温信号和皮肤电信号的任一项。
本实施例中,采集监测用户睡眠过程的脑电信号和心电信号作为睡眠生理状态,来陈述本技术方案的具体实施过程。
首先,通过脑电图机对用户睡眠脑电进行采集记录,采样率为1024Hz,记录电极为F3、F4、C3和C4,参考电极M1、M2;对脑电信号进行统一信号处理,包括以M1/2进行重参考、去伪迹、小波降噪、50Hz倍频工频陷波、0.5-80hz带通滤波和信号矫正处理,得到纯净脑电信号。通过便携式单导心电仪对用户心电信号进行采集提取,采集位置在左胸上方,采样率512Hz;对心电信号进行统一信号处理,包括去伪迹、小波降噪、0.5-40hz带通滤波和信号矫正处理,得到纯净心电信号。
其次,依次对纯净脑电信号进行信号频带提取,包括δ节律(0.5-4Hz)、θ节律(4-8Hz)、α节律(8-12Hz)、β节律(12-30Hz)、γ节律(30-80Hz),得到频带脑电信号;进而,以预设时间窗口长度30s和预设时间平移步长15s对纯净脑电信号、频带脑电信号和纯净心电信号进行连续滑动分割,得到用户睡眠生理状态时帧数据。
第二步、对用户睡眠生理状态时帧数据进行特征分析和特征融合,得到睡眠深度特征曲线。
本实施例中,特征分析至少包括数值分析、包络分析、时频分析、熵分析、分形分析和复杂度分析。特征融合是指从特征分析得到的目标特征集中,选择预设特征数量的目标特征并进行加权计算,生成睡眠深度特征曲线。
本实施例中,睡眠深度特征曲线为表征用户在预设入睡前时期、睡眠持续期、预设结束睡眠后时期的睡眠深度和时相状态的特征曲线,其计算生成方法具体为:
1)对用户睡眠生理状态时帧数据中时帧数据逐个进行特征分析,按时序拼接得到睡眠生理状态时帧特征曲线集;
2)从睡眠生理状态时帧特征曲线集中,筛选目标相关特征曲线,并按照预设特征融合权重进行加权计算,生成睡眠深度特征曲线。
本实施例中, 逐帧对用户睡眠生理状态时帧数据的脑电数据进行时频分析(频带功率、频带功率占比)、熵分析(样本熵)和复杂度分析(LZC指数:Lempel-Ziv复杂性指数);逐帧对用户睡眠生理状态时帧数据的心电数据进行数值分析,提取用户心率变化特征(心率平均值、心率变异系数)。进而,选择F4-M1通道的δ-θ(δ节律+θ节律)联合频带功率占比、归一化心率平均值倒数的均值(特征融合过程),作为睡眠深度特征曲线。通常情况下,用户睡眠越深,其δ-θ联合频带功率占比越大,归一化的心率平均值倒数越大(反面,心率平均值倒数越小),用户睡眠状态、皮层电生理和自主神经生理表现越稳定规律。
第三步、对用户睡眠生理状态时帧数据进行时相识别和睡眠分期,得到睡眠分期曲线。
本实施例中,睡眠分期曲线的生成方法具体为:
1)通过深度学习算法对规模睡眠用户样本的用户睡眠生理状态时帧数据及其对应睡眠分期数据进行学习训练和数据建模,得到睡眠时相自动分期模型;
2)将当前用户的用户睡眠生理状态时帧数据输入睡眠时相自动分期模型,得到所对应的睡眠时相分期并按时序生成睡眠分期曲线。
在实际使用场景中,通过用户样本的数据积累和分期模型的深度学习,睡眠时相自动分期模型的准确性会越来越高。
P200:根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,反向提取趋势分量,提取睡眠趋势特征曲线。
第一步、根据睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,得到FA波动分量。
本实施例中,窗口尺度的选取取决于所述睡眠深度特征曲线的特征来源和组合方式、动态调节效果。拟合方式至少包括最小二乘法拟合、多项式拟合、线性拟合中的任意一项,并由所述睡眠深度特征曲线的特征来源和组合方式共同决定。
本实施例中,去趋势分析的方法至少包括去趋势分析DFA、多重分形去趋势分析MFDFA、非对称消除趋势波动分析ADFA中的任意一项。
本实施例中,选择多重分形去趋势分析MFDFA作为去趋势分析的方法。多重分形去趋势分析MFDFA采用不同时间尺度来量化时间序列的相关性,能够有效去除局部趋势和波动噪声对时间序列标度的影响,有效提取不同时间尺度下时间序列的分形特征,非常适合非平稳时间序列长期趋势变化规律的分析和提取。其主要步骤如下:
1)对于长度为的时间序列/>,构造其取均值的和序列,/>提取公式如下
;/>
其中,为/>的均值。
2)将序列划分为长度为/>的/>个非重叠的区间(即改变时间尺度),得到区间集。其中,
;
3)对区间集中的每个区间/>内的/>个点,用最小二乘法进行/>阶多项式拟合,得到局部拟合趋势函数或趋势信号/>,进而得到消除趋势后的波动函数或波动信号,即均方误差/>。其中,
;
;
4)计算阶波动函数/>,其中
;
;
值得一提的是,是关于数据长度/>和分形阶数/>的函数,随着/>的增大,/>呈现幂律关系增加,最后可以提取得到赫斯特指数/>。当/>时,/>就是标准的去趋势分析DFA。在实际使用过程中,为保证序列/>的信息在划分过程中不丢失,可以按照/>升序(由小到大)和降序(由大到小)分别对序列/>进行划分,共得到/>个非重叠的区间。同时,为了保证/>的稳定性,通常/>取值在4和/>之间。
本实施例中,睡眠深度特征曲线由δ-θ(δ节律+θ节律)联合频带功率占比、归一化心率平均值倒数计算而来,选择窗口尺度为15,最小二乘法拟合为拟合方式,通过多重分形去趋势分析MFDFA提取睡眠深度特征曲线的FA波动分量。窗口尺度和拟合方式的选取,需要考虑睡眠深度特征曲线的计算方式和特征来源,不同特征来源和特征组合方式会有最佳的窗口尺度和最佳的拟合方式,以得到最优的分析变换结果。
第二步、根据睡眠深度特征曲线和FA波动分量,反向提取趋势分量,生成睡眠趋势特征曲线。
本实施例中,睡眠趋势特征曲线的计算公式具体为:
;
其中,分别为睡眠趋势特征曲线、睡眠深度特征曲线和FA波动分量。
P300:根据所述睡眠分期曲线、所述睡眠深度特征曲线和所述睡眠趋势特征曲线,动态计算提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征。
第一步、根据睡眠深度特征曲线和睡眠趋势特征曲线,计算得到睡眠趋势性指数。
本实施例中,睡眠趋势性指数的计算方法具体为:
1)获取睡眠深度特征曲线和睡眠趋势特征曲线;
2)分别对睡眠深度特征曲线和睡眠趋势特征曲线进行取平方,得到睡眠深度特征平方曲线和睡眠趋势特征平方曲线;
3)计算睡眠深度特征平方曲线和睡眠趋势特征平方曲线的样本点比值,得到睡眠趋势性曲线;
4)求取睡眠趋势性曲线的平均值,得到睡眠趋势性特征值;
5)计算睡眠趋势性特征值、去趋势分析方法对应的预设方法修正系数、与用户生物状态信息相关的预设用户个性修正系数的乘积,生成睡眠趋势性指数。
本实施例中,用户生物状态信息至少包括性别、年龄、职业、健康状态、疾病状态和教育程度。
本实施例中,睡眠趋势性指数的计算公式具体为:
;
其中,为睡眠趋势性指数,/>分别为预设用户个性修正系数和预设方法修正系数,/>、/>分别为睡眠趋势特征曲线和所述睡眠深度特征曲线第i点的值,/>为睡眠趋势特征曲线的数据长度。
通常情况下,MFDFA预设方法修正系数为0.95,正常健康用户的预设用户个性修正系数为1.0。
第二步、根据睡眠分期曲线和睡眠趋势特征曲线进行相关性计算,得到时相趋势相关系数。
本实施例中,相关性计算的方法至少包括相干性分析、皮尔逊相关分析、杰卡德相似分析、线性互信息分析、线性相关分析、欧氏距离分析、曼哈顿距离分析、切比雪夫距离分析中的任意一项。
本实施例中,先对睡眠分期曲线进行曲线平滑操作;再选择线性相关分析作为相关性计算并得到线性相关系数,作为睡眠分期曲线和睡眠趋势特征曲线的时相趋势相关系数。
第三步、基于睡眠分期曲线,对睡眠趋势特征曲线进行时相分布统计,得到时相趋势分布特征。
本实施例中,时相分布统计具体为按照睡眠分期曲线中的睡眠时相分期,对睡眠趋势特征曲线中的睡眠趋势特征值进行数值分布统计分析,得到睡眠趋势特征曲线的数值分布统计特征。时相趋势分布特征至少包括平均值、均方根、最大值、最小值、方差、标准差、变异系数、峰度和偏度中的任意一项。
P400:建立并动态更新用户个性化睡眠趋势性数据库。
第一步、建立用户个性化睡眠趋势性数据库。
本实施例中,用户个性化睡眠趋势性数据库至少包括用户生物状态信息、睡眠场景、睡眠分期曲线、睡眠深度特征曲线、睡眠趋势特征曲线、睡眠趋势性指数、时相趋势相关系数、时相趋势分布特征、调节方式、执行方式、调节方法、目标调节值曲线、调节强度曲线、动态调节效果曲线、综合调节效果指数。
本实施例中,用户个性化睡眠趋势性数据库用mysql存储。
第二步、将用户睡眠趋稳性动态检测量化的关键过程数据,动态更新入用户个性化睡眠趋势性数据库。
本实施例中,主要将检测量化过程中产生的用户生物状态信息、睡眠场景、睡眠分期曲线、睡眠深度特征曲线、睡眠趋势特征曲线、睡眠趋势性指数、时相趋势相关系数、时相趋势分布特征等数据,实时更新入用户个性化睡眠趋势性数据库。
P500:根据所述用户个性化睡眠趋势性数据库,对用户睡眠行为进行动态预测,提取个性时相场景睡眠调节策略并对用户睡眠过程进行干预调节。
第一步、根据用户个性化睡眠趋势性数据库,对睡眠趋势特征曲线进行动态预测,提取睡眠趋势预测特征值,生成睡眠趋势预测特征曲线。
本实施例中,睡眠趋势预测特征值的预测方法至少包括AR、MR、ARMA、ARIMA、SARIMA、VAR、深度学习中的任一项。
本实施例中,使用ARMA方法对睡眠深度特征曲线进行趋势预测分析,得到睡眠趋势预测特征值,并生成或更新得到睡眠趋势预测特征曲线。
在实际适应场景中,趋势分析和指数预测可能采用AR、MR、ARMA、ARIMA、SARIMA、VAR等等常用的时间按序列预测方法,也可以通过深度学习模型来完成睡眠趋势预测特征值的预测计算。
第二步、根据用户个性化睡眠趋势性数据库、预设睡眠趋稳性调节知识库、睡眠趋势预测特征值和当前具体睡眠场景,按照预设动态调节周期,动态生成个性时相场景睡眠调节策略。
本实施例中,个性时相场景睡眠调节策略至少包括睡眠场景、睡眠时相、调节方式、执行方式、调节方法、调节强度、调节时点、持续时间、目标调节值和装置控制参数;其中,调节方式至少包括声乐刺激、超声刺激、光刺激、电刺激、磁刺激、温度刺激、湿度刺激、触觉刺激和调控中的任一方式,执行方式至少包括离体式和接触式中的任一方式。
在实际使用场景中,在保证调节效果的情况下,通过选择离体的声乐刺激、光刺激、温度刺激、湿度刺激和和调控,对用户睡眠干扰较小且用户体验较好。
第三步、根据个性时相场景睡眠调节策略,连接并调控睡眠趋势性调节外围设备,对用户睡眠过程进行动态干预调节。
本实施例中,调控睡眠趋势性调节外围设备至少包括声乐刺激设备、超声刺激设备、光刺激设备、电刺激设备、磁刺激设备、温度刺激设备、湿度刺激设备、触觉刺激设备和调控设备中的任一项,并由具体的调节方式决定。
第四步、对动态干预调节的效果进行动态跟踪评估,提取动态调节效果系数并生成动态调节效果曲线,计算得到综合调节效果指数。
本实施例中,动态调节效果系数的具体计算公式具体为:
;
其中,为所述动态调节效果系数;/>分别为所述个性时相场景睡眠调节策略中的目标调节值、所述睡眠趋势预测特征值、动态调节后量化得到的所述睡眠趋势特征曲线中的睡眠趋势特征值;/>为预设用户个性修正系数。
本实施例中,动态调节效果系数反向应用于去趋势分析方法参数中窗口尺度的动态调整、个性时相场景睡眠调节策略的生成,以不断优化检测量化和干预调节的闭环循环效率。在实际使用场景中,可以通过窗口尺度的缩短或延长,来不断测试和验证动态调节的效果。
在实际使用场景中,动态干预调节的效果也可以通过睡眠趋势预测特征曲线、睡眠趋势特征曲线和目标调节值曲线的相关性计算、曲线距离特征计算和综合评价来实现;如通过睡眠趋势特征曲线和目标调节值曲线的欧氏距离的平均值,就能够准确地评估动态调节效果。
本实施例中,综合调节效果指数具体为动态调节效果曲线的平均值。
P600:进一步动态更新所述用户个性化睡眠趋势性数据库并优化所述个性时相场景睡眠调节策略,按照预设报告周期生成睡眠趋势性量化及调节报告的具体步骤还包括:
第一步、对用户不同睡眠场景下的睡眠趋势性进行动态检测量化和动态干预调节,持续动态更新用户个性化睡眠趋势性数据库。
在实际使用场景中,根据用户基本情况,选择不同场景组合,对用户在不同睡眠压力、不同睡眠环境、不同健康状态等多场景下的睡眠趋势性进行动态检测量化和动态干预调节,能够得到更全面的用户个性化睡眠趋势性数据。
第三步、将动态检测量化和动态干预调节的关键过程数据,动态更新入用户个性化睡眠趋势性数据库,动态优化个性时相场景睡眠调节策略。
本实施例中,需要将用户关键生理指标曲线和生理事件、睡眠场景、睡眠分期曲线、睡眠深度特征曲线、睡眠趋势特征曲线、睡眠趋势性指数、时相趋势相关系数、时相趋势分布特征、调节方式、执行方式、调节方法、目标调节值曲线、调节强度曲线、动态调节效果曲线等关键过程数据,都实时地更新入用户个性化睡眠趋势性数据库,也为实时生成的个性时相场景睡眠调节策略提供精确数据基础。
在实际使用场景中,随着用户个性化相关数据的不断积累和场景调节反馈,用户个性化睡眠趋势性数据库的数据丰富性增加,这不仅能够进一步准确检测量化用户睡眠趋稳性,还能够不断提高用户睡眠趋稳性动态调节干预的质量和效果。
第三步、按照预设报告周期,生成睡眠趋势性量化及调节报告。
本实施例中,睡眠趋势性量化及调节报告至少包括用户生物状态信息、睡眠场景、睡眠分期曲线、睡眠深度特征曲线、睡眠趋势特征曲线、睡眠趋势性指数、时相趋势相关系数、时相趋势分布特征、综合调节效果指数、动态调节效果曲线、用户睡眠趋势性量化及调节总结。
在实际使用场景中,可以按不同时间周期来生成并输出睡眠趋势性量化及调节报告,来满足不同人群的不同场景需求,为用户睡眠健康管理提供健康数据统计和策略改进依据。
如图2所示,本发明实施例提供的一种基于趋势分析的睡眠趋稳性量化及调节系统,所述系统被构造以用于执行上述各个方法步骤。所述系统包括如下模块:
状态检测分析模块S100,用于对用户睡眠过程生理状态进行采集监测处理、特征分析和时相识别,得到睡眠深度特征曲线和睡眠分期曲线;
趋势特征提取模块S200,用于根据睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,反向提取趋势分量,提取睡眠趋势特征曲线;
趋势动态量化模块S300,用于根据睡眠分期曲线、睡眠深度特征曲线和睡眠趋势特征曲线,动态计算提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征;
趋势数据库模块S400,用于建立并动态更新用户个性化睡眠趋势性数据库;
趋势动态调节模块S500,用于根据用户个性化睡眠趋势性数据库,对用户睡眠行为进行动态预测,提取个性时相场景睡眠调节策略并对用户睡眠过程进行干预调节;
数据动态应用模块S600,用于动态更新用户个性化睡眠趋势性数据库并优化个性时相场景睡眠调节策略,按照预设报告周期生成睡眠趋势性量化及调节报告;
数据运行管理模块S700,用于对系统的全部过程数据进行可视化管理、统一存储和运行管理。
本实施例中,状态检测分析模块S100还包括以下功能单元:
状态检测处理单元,用于对用户睡眠过程生理状态进行采集监测和信号处理,得到用户睡眠生理状态时帧数据;
深度特征提取单元,用于对用户睡眠生理状态时帧数据进行特征分析和特征融合,得到睡眠深度特征曲线;
睡眠时相识别单元,用于对用户睡眠生理状态时帧数据进行时相识别和睡眠分期,得到睡眠分期曲线。
本实施例中,趋势特征提取模块S200还包括以下功能单元:
去趋势分析单元,用于根据睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,得到FA波动分量;
趋势分量识别单元,用于根据睡眠深度特征曲线和FA波动分量,反向提取趋势分量,生成睡眠趋势特征曲线。
本实施例中,趋势动态量化模块S300还包括以下功能单元:
趋势指数量化单元,用于根据睡眠深度特征曲线和睡眠趋势特征曲线,计算得到睡眠趋势性指数;
相关性计算单元,用于根据睡眠分期曲线和睡眠趋势特征曲线进行相关性计算,得到时相趋势相关系数;
时相分布统计单元,用于基于睡眠分期曲线,对睡眠趋势特征曲线进行时相分布统计,得到时相趋势分布特征。
本实施例中,趋势数据库模块S400还包括以下功能单元:
数据库建立单元,用于建立用户个性化睡眠趋势性数据库;
数据库更新单元,用于将用户睡眠趋稳性动态检测量化的关键过程数据,动态更新入用户个性化睡眠趋势性数据库。
本实施例中,趋势动态调节模块S500还包括以下功能单元:
睡眠状态预测单元,用于根据用户个性化睡眠趋势性数据库,对睡眠趋势特征曲线进行动态趋势预测,提取睡眠趋势预测特征值,生成睡眠趋势预测特征曲线;
调节策略生成单元,用于根据用户个性化睡眠趋势性数据库、预设睡眠趋稳性调节知识库、睡眠趋势预测特征值和当前具体睡眠场景,按照预设动态调节周期,动态生成个性时相场景睡眠调节策略;
动态调控执行单元,用于根据个性时相场景睡眠调节策略,连接并调控睡眠趋势性调节外围设备,对用户睡眠过程进行动态干预调节;
动态效果评估单元,用于对动态干预调节的效果进行动态跟踪评估,提取动态调节效果系数并生成动态调节效果曲线,计算得到综合调节效果指数。
本实施例中,数据动态应用模块S600还包括以下功能单元:
趋势数据管理单元,用于对用户不同睡眠场景下的睡眠趋势性进行动态检测量化和动态干预调节,持续动态更新用户个性化睡眠趋势性数据库;
趋势数据应用单元,用于将动态检测量化和动态干预调节的关键过程数据,动态更新入用户个性化睡眠趋势性数据库,动态优化个性时相场景睡眠调节策略;
用户报告生成单元,用于按照预设报告周期,生成睡眠趋势性量化及调节报告;
用户报告输出单元,用于对睡眠趋势性量化及调节报告的格式输出、展现形式进行统一管理。
本实施例中,数据运行管理模块S700还包括以下功能单元:
用户信息管理单元,用于用户基本信息的登记输入、编辑、查询、输出和删除;
数据可视化管理单元,用于对系统中所有数据的可视化展现管理;
数据存储管理单元,用于对系统中所有数据的统一存储管理;
数据运营管理单元,用于对系统中所有数据的备份、迁移和导出。
如图3所示,本发明实施例提供的一种基于趋势分析的睡眠趋稳性量化及调节装置,包括以下模组:
状态检测分析模组M100,用于对用户睡眠过程生理状态进行采集监测处理、特征分析和时相识别,得到睡眠深度特征曲线和睡眠分期曲线;
趋势特征提取模组M200,用于根据睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,反向提取趋势分量,提取睡眠趋势特征曲线;
趋势动态量化模组M300,用于根据睡眠分期曲线、睡眠深度特征曲线和睡眠趋势特征曲线,动态计算提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征;
趋势数据库模组M400,用于建立并动态更新用户个性化睡眠趋势性数据库;
趋势动态调节模组M500,用于根据用户个性化睡眠趋势性数据库,对用户睡眠行为进行动态预测,提取个性时相场景睡眠调节策略并对用户睡眠过程进行干预调节;
数据动态应用模组M600,用于动态更新用户个性化睡眠趋势性数据库并优化个性时相场景睡眠调节策略,按照预设报告周期生成睡眠趋势性量化及调节报告;
数据可视化模组M700,用于对装置中所有过程数据和/或结果数据的统一可视化展示管理;
数据管理中心模组M800,用于对装置中所有过程数据和/或结果数据的统一存储和数据运营管理。
所述装置被构造以用于对应执行图1的方法中的各个步骤,在此不再赘述。
本发明还提供了可编程的各类处理器(FPGA、ASIC或其他集成电路),所述处理器用于运行程序,其中,所述程序运行时执行上述实施例中的步骤。
本发明还提供了对应的计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述存储器执行所述程序时实现上述实施例中的步骤。
虽然本发明所揭露的实施方式如上,但所述的内容仅为便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属领域内的技术人员,在不脱离本发明所揭露的精神和原则的前提下,可以在实施的形式及细节上进行任何的修改与变化、等同替换等,这些都属于本发明的保护范围。因此,本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (35)
1.一种基于趋势分析的睡眠趋稳性量化及调节方法,其特征在于,包括以下步骤:
对用户睡眠过程生理状态进行采集监测处理、特征分析和时相识别,得到睡眠深度特征曲线和睡眠分期曲线;
根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,反向提取趋势分量,提取睡眠趋势特征曲线;
根据所述睡眠分期曲线、所述睡眠深度特征曲线和所述睡眠趋势特征曲线,动态计算提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征;
建立并动态更新用户个性化睡眠趋势性数据库;
根据所述用户个性化睡眠趋势性数据库,对用户睡眠行为进行动态预测,提取个性时相场景睡眠调节策略并对用户睡眠过程进行干预调节;
进一步动态更新所述用户个性化睡眠趋势性数据库并优化所述个性时相场景睡眠调节策略,按照预设报告周期生成睡眠趋势性量化及调节报告;
所述根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,反向提取趋势分量,提取睡眠趋势特征曲线的具体步骤还包括:
根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,得到FA波动分量;
根据所述睡眠深度特征曲线和所述FA波动分量,反向提取趋势分量,生成所述睡眠趋势特征曲线;
所述根据所述睡眠分期曲线、所述睡眠深度特征曲线和所述睡眠趋势特征曲线,动态计算提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征的具体步骤还包括:
根据所述睡眠深度特征曲线和所述睡眠趋势特征曲线,计算得到所述睡眠趋势性指数;
根据所述睡眠分期曲线和所述睡眠趋势特征曲线进行相关性计算,得到所述时相趋势相关系数;
基于所述睡眠分期曲线,对所述睡眠趋势特征曲线进行时相分布统计,得到所述时相趋势分布特征;
所述睡眠趋势性指数的计算方法具体为:
1)获取所述睡眠深度特征曲线和所述睡眠趋势特征曲线;
2)分别对所述睡眠深度特征曲线和所述睡眠趋势特征曲线进行取平方,得到睡眠深度特征平方曲线和睡眠趋势特征平方曲线;
3)计算所述睡眠深度特征平方曲线和所述睡眠趋势特征平方曲线的样本点比值,得到睡眠趋势性曲线;
4)求取所述睡眠趋势性曲线的平均值,得到睡眠趋势性特征值;
5)计算所述睡眠趋势性特征值、去趋势分析方法对应的预设方法修正系数、与用户生物状态信息相关的预设用户个性修正系数的乘积,生成所述睡眠趋势性指数。
2.如权利要求1所述的方法,其特征在于,所述对用户睡眠过程生理状态进行采集处理和特征分析,得到睡眠深度特征曲线和睡眠分期曲线的具体步骤还包括:
对用户睡眠过程生理状态进行采集监测和信号处理,得到用户睡眠生理状态时帧数据;
对所述用户睡眠生理状态时帧数据进行特征分析和特征融合,得到所述睡眠深度特征曲线;
对所述用户睡眠生理状态时帧数据进行时相识别和睡眠分期,得到所述睡眠分期曲线。
3.如权利要求2所述的方法,其特征在于:所述信号处理至少包括AD模数转换、重采样、重参考、降噪、去伪迹、工频陷波、低通滤波、高通滤波、带阻滤波、带通滤波、矫正处理和时帧分割;其中,所述矫正处理具体为对信号中的包含伪迹或失真的信号数据片段进行信号矫正、信号预测和平滑处理,所述时帧分割是指根据信号采样率,以预设时间窗口长度和预设时间平移步长对目标信号数据进行连续滑动分割。
4.如权利要求2所述的方法,其特征在于:所述用户睡眠生理状态时帧数据至少包括脑中枢状态数据、自主神经状态数据中的任一项;其中,所述脑中枢状态数据至少包括脑电图信号、脑磁图信号、血氧水平依赖信号中的任一项,所述自主神经状态数据至少包括血氧水平依赖信号、心电信号、脉搏信号、呼吸信号、血氧信号、体温信号和皮肤电信号的任一项。
5.如权利要求2或4所述的方法,其特征在于:所述特征分析至少包括数值分析、包络分析、时频分析、熵分析、分形分析和复杂度分析。
6.如权利要求2或4所述的方法,其特征在于:所述特征融合是指从所述特征分析得到的目标特征集中,选择预设特征数量的目标特征并进行加权计算,生成所述睡眠深度特征曲线。
7.如权利要求2或4所述的方法,其特征在于,所述睡眠深度特征曲线具体为表征用户在预设入睡前时期、睡眠持续期、预设结束睡眠后时期的睡眠深度和时相状态的特征曲线,其计算生成方法具体为:
1)对所述用户睡眠生理状态时帧数据中时帧数据逐个进行特征分析,按时序拼接得到睡眠生理状态时帧特征曲线集;
2)从所述睡眠生理状态时帧特征曲线集中,筛选目标相关特征曲线,并按照预设特征融合权重进行加权计算,生成所述睡眠深度特征曲线。
8.如权利要求2或4所述的方法,其特征在于,所述睡眠分期曲线的生成方法具体为:
1)通过深度学习算法对规模睡眠用户样本的所述用户睡眠生理状态时帧数据及其对应睡眠分期数据进行学习训练和数据建模,得到睡眠时相自动分期模型;
2)将当前用户的所述用户睡眠生理状态时帧数据输入所述睡眠时相自动分期模型,得到所对应的睡眠时相分期并按时序生成所述睡眠分期曲线。
9.如权利要求1所述的方法,其特征在于:所述窗口尺度的选取取决于所述睡眠深度特征曲线的特征来源和组合、动态调节效果方式。
10.如权利要求1所述的方法,其特征在于:所述拟合方式至少包括最小二乘法拟合、多项式拟合、线性拟合中的任意一项,并由所述睡眠深度特征曲线的特征来源和组合方式共同决定。
11.如权利要求1所述的方法,其特征在于:所述去趋势分析的方法至少包括去趋势分析DFA、多重分形去趋势分析MFDFA、非对称消除趋势波动分析ADFA中的任意一项。
12.如权利要求1所述的方法,其特征在于,所述睡眠趋势特征曲线的计算公式具体为:
;
其中,分别为所述睡眠趋势特征曲线、所述睡眠深度特征曲线和所述FA波动分量。
13.如权利要求1所述的方法,其特征在于,所述用户生物状态信息至少包括性别、年龄、职业、健康状态、疾病状态和教育程度。
14.如权利要求1或13所述的方法,其特征在于,所述睡眠趋势性指数的计算公式具体为:
;
其中,为所述睡眠趋势性指数,/>分别为预设用户个性修正系数和预设方法修正系数,/>、/>分别为所述睡眠趋势特征曲线和所述睡眠深度特征曲线,/>为所述睡眠趋势特征曲线的数据长度。
15.如权利要求1所述的方法,其特征在于,所述相关性计算的方法至少包括相干性分析、皮尔逊相关分析、杰卡德相似分析、线性互信息分析、线性相关分析、欧氏距离分析、曼哈顿距离分析、切比雪夫距离分析中的任意一项。
16.如权利要求1所述的方法,其特征在于,所述时相分布统计具体为按照所述睡眠分期曲线中的睡眠时相分期,对所述睡眠趋势特征曲线中的睡眠趋势特征值进行数值分布统计分析,得到所述睡眠趋势特征曲线的数值分布统计特征。
17.如权利要求16所述的方法,其特征在于,所述时相趋势分布特征至少包括平均值、均方根、最大值、最小值、方差、标准差、变异系数、峰度和偏度中的任意一项。
18.如权利要求1所述的方法,其特征在于,所述建立并动态更新用户个性化睡眠趋势性数据库的具体步骤还包括:
建立所述用户个性化睡眠趋势性数据库;
将用户睡眠趋稳性动态检测量化的关键过程数据,动态更新入所述用户个性化睡眠趋势性数据库。
19.如权利要求18所述的方法,其特征在于:所述用户个性化睡眠趋势性数据库至少包括用户生物状态信息、睡眠场景、所述睡眠分期曲线、所述睡眠深度特征曲线、所述睡眠趋势特征曲线、所述睡眠趋势性指数、所述时相趋势相关系数、所述时相趋势分布特征、调节方式、执行方式、调节方法、目标调节值曲线、调节强度曲线、动态调节效果曲线、综合调节效果指数。
20.如权利要求1所述的方法,其特征在于:所述根据所述用户个性化睡眠趋势性数据库,对用户睡眠行为进行动态预测,提取个性时相场景睡眠调节策略并对用户睡眠过程进行干预调节的具体步骤还包括:
根据所述用户个性化睡眠趋势性数据库,对所述睡眠趋势特征曲线进行动态趋势预测,提取睡眠趋势预测特征值,生成睡眠趋势预测特征曲线;
根据所述用户个性化睡眠趋势性数据库、预设睡眠趋稳性调节知识库、所述睡眠趋势预测特征值和当前具体睡眠场景,按照预设动态调节周期,动态生成所述个性时相场景睡眠调节策略;
根据所述个性时相场景睡眠调节策略,连接并调控睡眠趋势性调节外围设备,对用户睡眠过程进行动态干预调节;
对动态干预调节的效果进行动态跟踪评估,提取动态调节效果系数并生成动态调节效果曲线,计算得到综合调节效果指数。
21.如权利要求20所述的方法,其特征在于:所述睡眠趋势预测特征值的预测方法至少包括AR、MA、ARMA、ARIMA、SARIMA、VAR、深度学习中的任一项。
22.如权利要求20所述的方法,其特征在于:所述个性时相场景睡眠调节策略至少包括睡眠场景、睡眠时相、调节方式、执行方式、调节方法、调节强度、调节时点、持续时间、目标调节值和装置控制参数;其中,所述调节方式至少包括声乐刺激、超声刺激、光刺激、电刺激、磁刺激、温度刺激、湿度刺激、触觉刺激和调控中的任一方式,所述执行方式至少包括离体式和接触式中的任一方式。
23.如权利要求22所述的方法,其特征在于:所述调控睡眠趋势性调节外围设备至少包括声乐刺激设备、超声刺激设备、光刺激设备、电刺激设备、磁刺激设备、温度刺激设备、湿度刺激设备、触觉刺激设备和调控设备中的任一项,并由具体的所述调节方式决定。
24.如权利要求20-23任一项所述的方法,其特征在于:所述动态调节效果系数的具体计算公式具体为:
;
其中,为所述动态调节效果系数;/>分别为所述个性时相场景睡眠调节策略中的目标调节值、所述睡眠趋势预测特征值、动态调节后量化得到的所述睡眠趋势特征曲线中的睡眠趋势特征值;/>为预设用户个性修正系数。
25.如权利要求20所述的方法,其特征在于:所述动态调节效果系数反向应用于去趋势分析方法参数中所述窗口尺度的动态调整、所述个性时相场景睡眠调节策略的生成。
26.如权利要求20所述的方法,其特征在于:所述综合调节效果指数具体为所述动态调节效果曲线的平均值或均方根。
27.如权利要求1所述的方法,其特征在于:所述进一步动态更新所述用户个性化睡眠趋势性数据库并优化所述个性时相场景睡眠调节策略,按照预设报告周期生成睡眠趋势性量化及调节报告的具体步骤还包括:
对用户不同睡眠场景下的睡眠趋势性进行动态检测量化和动态干预调节,持续动态更新所述用户个性化睡眠趋势性数据库;
将动态检测量化和动态干预调节的关键过程数据,动态更新入所述用户个性化睡眠趋势性数据库,动态优化所述个性时相场景睡眠调节策略;
按照预设报告周期,生成所述睡眠趋势性量化及调节报告。
28.如权利要求27所述的方法,其特征在于:所述睡眠趋势性量化及调节报告至少包括用户生物状态信息、睡眠场景、所述睡眠分期曲线、所述睡眠深度特征曲线、所述睡眠趋势特征曲线、所述睡眠趋势性指数、所述时相趋势相关系数、所述时相趋势分布特征、综合调节效果指数、动态调节效果曲线、用户睡眠趋势性量化及调节总结。
29.一种基于趋势分析的睡眠趋稳性量化及调节系统,其特征在于,包括以下模块:
状态检测分析模块,用于对用户睡眠过程生理状态进行采集监测处理、特征分析和时相识别,得到睡眠深度特征曲线和睡眠分期曲线;
趋势特征提取模块,用于根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,反向提取趋势分量,提取睡眠趋势特征曲线;
趋势动态量化模块,用于根据所述睡眠分期曲线、所述睡眠深度特征曲线和所述睡眠趋势特征曲线,动态计算提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征;
趋势数据库模块,用于建立并动态更新用户个性化睡眠趋势性数据库;
趋势动态调节模块,用于根据所述用户个性化睡眠趋势性数据库,对用户睡眠行为进行动态预测,提取个性时相场景睡眠调节策略并对用户睡眠过程进行干预调节;
数据动态应用模块,用于动态更新所述用户个性化睡眠趋势性数据库并优化所述个性时相场景睡眠调节策略,按照预设报告周期生成睡眠趋势性量化及调节报告;
数据运行管理模块,用于对所述系统的全部过程数据进行可视化管理、统一存储和运行管理;
所述趋势特征提取模块还包括以下功能单元:
去趋势分析单元,用于根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,得到FA波动分量;
趋势分量识别单元,用于根据所述睡眠深度特征曲线和所述FA波动分量,反向提取趋势分量,生成所述睡眠趋势特征曲线;
所述趋势动态量化模块还包括以下功能单元:
趋势指数量化单元,用于根据所述睡眠深度特征曲线和所述睡眠趋势特征曲线,计算得到所述睡眠趋势性指数;
相关性计算单元,用于根据所述睡眠分期曲线和所述睡眠趋势特征曲线进行相关性计算,得到所述时相趋势相关系数;
时相分布统计单元,用于基于所述睡眠分期曲线,对所述睡眠趋势特征曲线进行时相分布统计,得到所述时相趋势分布特征;
所述睡眠趋势性指数的计算方法具体为:
1)获取所述睡眠深度特征曲线和所述睡眠趋势特征曲线;
2)分别对所述睡眠深度特征曲线和所述睡眠趋势特征曲线进行取平方,得到睡眠深度特征平方曲线和睡眠趋势特征平方曲线;
3)计算所述睡眠深度特征平方曲线和所述睡眠趋势特征平方曲线的样本点比值,得到睡眠趋势性曲线;
4)求取所述睡眠趋势性曲线的平均值,得到睡眠趋势性特征值;
5)计算所述睡眠趋势性特征值、去趋势分析方法对应的预设方法修正系数、与用户生物状态信息相关的预设用户个性修正系数的乘积,生成所述睡眠趋势性指数。
30.如权利要求29所述的系统,其特征在于,所述状态检测分析模块还包括以下功能单元:
状态检测处理单元,用于对用户睡眠过程生理状态进行采集监测和信号处理,得到用户睡眠生理状态时帧数据;
深度特征提取单元,用于对所述用户睡眠生理状态时帧数据进行特征分析和特征融合,得到所述睡眠深度特征曲线;
睡眠时相识别单元,用于对所述用户睡眠生理状态时帧数据进行时相识别和睡眠分期,得到所述睡眠分期曲线。
31.如权利要求29所述的系统,其特征在于,所述趋势数据库模块还包括以下功能单元:
数据库建立单元,用于建立所述用户个性化睡眠趋势性数据库;
数据库更新单元,用于将用户睡眠趋稳性动态检测量化的关键过程数据,动态更新入所述用户个性化睡眠趋势性数据库。
32.如权利要求29所述的系统,其特征在于,所述趋势动态调节模块还包括以下功能单元:
睡眠状态预测单元,用于根据所述用户个性化睡眠趋势性数据库,对所述睡眠趋势特征曲线进行动态趋势预测,提取睡眠趋势预测特征值,生成睡眠趋势预测特征曲线;
调节策略生成单元,用于根据所述用户个性化睡眠趋势性数据库、预设睡眠趋稳性调节知识库、所述睡眠趋势预测特征值和当前具体睡眠场景,按照预设动态调节周期,动态生成所述个性时相场景睡眠调节策略;
动态调控执行单元,用于根据所述个性时相场景睡眠调节策略,连接并调控睡眠趋势性调节外围设备,对用户睡眠过程进行动态干预调节;
动态效果评估单元,用于对动态干预调节的效果进行动态跟踪评估,提取动态调节效果系数并生成动态调节效果曲线,计算得到综合调节效果指数。
33.如权利要求32所述的系统,其特征在于,所述数据动态应用模块还包括以下功能单元:
趋势数据管理单元,用于对用户不同睡眠场景下的睡眠趋势性进行动态检测量化和动态干预调节,持续动态更新所述用户个性化睡眠趋势性数据库;
趋势数据应用单元,用于将动态检测量化和动态干预调节的关键过程数据,动态更新入所述用户个性化睡眠趋势性数据库,动态优化所述个性时相场景睡眠调节策略;
用户报告生成单元,用于按照预设报告周期,生成所述睡眠趋势性量化及调节报告;
用户报告输出单元,用于对所述睡眠趋势性量化及调节报告的格式输出、展现形式进行统一管理。
34.如权利要求29所述的系统,其特征在于,所述数据运行管理模块还包括以下功能单元:
用户信息管理单元,用于用户基本信息的登记输入、编辑、查询、输出和删除;
数据可视化管理单元,用于对所述系统中所有数据的可视化展现管理;
数据存储管理单元,用于对所述系统中所有数据的统一存储管理;
数据运营管理单元,用于对所述系统中所有数据的备份、迁移和导出。
35.一种基于趋势分析的睡眠趋稳性量化及调节装置,其特征在于,包括以下模组:
状态检测分析模组,用于对用户睡眠过程生理状态进行采集监测处理、特征分析和时相识别,得到睡眠深度特征曲线和睡眠分期曲线;
趋势特征提取模组,用于根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,反向提取趋势分量,提取睡眠趋势特征曲线;
趋势动态量化模组,用于根据所述睡眠分期曲线、所述睡眠深度特征曲线和所述睡眠趋势特征曲线,动态计算提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征;
趋势数据库模组,用于建立并动态更新用户个性化睡眠趋势性数据库;
趋势动态调节模组,用于根据所述用户个性化睡眠趋势性数据库,对用户睡眠行为进行动态预测,提取个性时相场景睡眠调节策略并对用户睡眠过程进行干预调节;
数据动态应用模组,用于动态更新所述用户个性化睡眠趋势性数据库并优化所述个性时相场景睡眠调节策略,按照预设报告周期生成睡眠趋势性量化及调节报告;
数据可视化模组,用于对所述装置中所有过程数据和/或结果数据的统一可视化展示管理;
数据管理中心模组,用于对所述装置中所有过程数据和/或结果数据的统一存储和数据运营管理;
所述根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,反向提取趋势分量,提取睡眠趋势特征曲线的具体步骤还包括:
根据所述睡眠深度特征曲线的特性选取窗口尺度和拟合方式并对其进行去趋势分析,得到FA波动分量;
根据所述睡眠深度特征曲线和所述FA波动分量,反向提取趋势分量,生成所述睡眠趋势特征曲线;
所述根据所述睡眠分期曲线、所述睡眠深度特征曲线和所述睡眠趋势特征曲线,动态计算提取睡眠趋势性指数、时相趋势相关系数和时相趋势分布特征的具体步骤还包括:
根据所述睡眠深度特征曲线和所述睡眠趋势特征曲线,计算得到所述睡眠趋势性指数;
根据所述睡眠分期曲线和所述睡眠趋势特征曲线进行相关性计算,得到所述时相趋势相关系数;
基于所述睡眠分期曲线,对所述睡眠趋势特征曲线进行时相分布统计,得到所述时相趋势分布特征;
所述睡眠趋势性指数的计算方法具体为:
1)获取所述睡眠深度特征曲线和所述睡眠趋势特征曲线;
2)分别对所述睡眠深度特征曲线和所述睡眠趋势特征曲线进行取平方,得到睡眠深度特征平方曲线和睡眠趋势特征平方曲线;
3)计算所述睡眠深度特征平方曲线和所述睡眠趋势特征平方曲线的样本点比值,得到睡眠趋势性曲线;
4)求取所述睡眠趋势性曲线的平均值,得到睡眠趋势性特征值;
5)计算所述睡眠趋势性特征值、去趋势分析方法对应的预设方法修正系数、与用户生物状态信息相关的预设用户个性修正系数的乘积,生成所述睡眠趋势性指数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310658830.2A CN116392085B (zh) | 2023-06-06 | 2023-06-06 | 基于趋势分析的睡眠趋稳性量化及调节方法、系统和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310658830.2A CN116392085B (zh) | 2023-06-06 | 2023-06-06 | 基于趋势分析的睡眠趋稳性量化及调节方法、系统和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116392085A CN116392085A (zh) | 2023-07-07 |
CN116392085B true CN116392085B (zh) | 2023-09-12 |
Family
ID=87018333
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310658830.2A Active CN116392085B (zh) | 2023-06-06 | 2023-06-06 | 基于趋势分析的睡眠趋稳性量化及调节方法、系统和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116392085B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117235678B (zh) * | 2023-11-14 | 2024-01-26 | 山东交通学院 | 一种木材切割过程中刀片磨损程度的检测方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001061820A (ja) * | 1999-08-26 | 2001-03-13 | Matsushita Electric Works Ltd | 睡眠状態判定方法及び睡眠状態判定装置 |
US9471544B1 (en) * | 2012-05-24 | 2016-10-18 | Google Inc. | Anomaly detection in a signal |
CN111631682A (zh) * | 2020-04-23 | 2020-09-08 | 平安国际智慧城市科技股份有限公司 | 基于去趋势分析的生理特征集成方法、装置和计算机设备 |
CN115881305A (zh) * | 2023-03-03 | 2023-03-31 | 安徽星辰智跃科技有限责任公司 | 一种睡眠趋稳性检测量化及辅助干预的方法、系统和装置 |
CN115944275A (zh) * | 2023-02-15 | 2023-04-11 | 中科闪耀(苏州)健康科技有限公司 | 一种通过呼吸率变异性判断睡眠分期的方法 |
CN116058804A (zh) * | 2023-03-27 | 2023-05-05 | 安徽星辰智跃科技有限责任公司 | 一种睡眠情绪活跃水平动态调节的方法、系统和装置 |
CN116058805A (zh) * | 2023-03-29 | 2023-05-05 | 安徽星辰智跃科技有限责任公司 | 一种睡眠记忆活跃水平动态调节的方法、系统和装置 |
CN116098578A (zh) * | 2021-11-09 | 2023-05-12 | 深圳迈瑞生物医疗电子股份有限公司 | 睡眠监测方法及睡眠监测系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014176286A1 (en) * | 2013-04-22 | 2014-10-30 | The Regents Of The University Of California | Fractal index analysis of human electroencephalogram signals |
US11439344B2 (en) * | 2015-07-17 | 2022-09-13 | Origin Wireless, Inc. | Method, apparatus, and system for wireless sleep monitoring |
KR102004052B1 (ko) * | 2017-11-09 | 2019-07-25 | 재단법인대구경북과학기술원 | 신체적 안정 판단 시스템 및 방법 |
-
2023
- 2023-06-06 CN CN202310658830.2A patent/CN116392085B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001061820A (ja) * | 1999-08-26 | 2001-03-13 | Matsushita Electric Works Ltd | 睡眠状態判定方法及び睡眠状態判定装置 |
US9471544B1 (en) * | 2012-05-24 | 2016-10-18 | Google Inc. | Anomaly detection in a signal |
CN111631682A (zh) * | 2020-04-23 | 2020-09-08 | 平安国际智慧城市科技股份有限公司 | 基于去趋势分析的生理特征集成方法、装置和计算机设备 |
CN116098578A (zh) * | 2021-11-09 | 2023-05-12 | 深圳迈瑞生物医疗电子股份有限公司 | 睡眠监测方法及睡眠监测系统 |
CN115944275A (zh) * | 2023-02-15 | 2023-04-11 | 中科闪耀(苏州)健康科技有限公司 | 一种通过呼吸率变异性判断睡眠分期的方法 |
CN115881305A (zh) * | 2023-03-03 | 2023-03-31 | 安徽星辰智跃科技有限责任公司 | 一种睡眠趋稳性检测量化及辅助干预的方法、系统和装置 |
CN116058804A (zh) * | 2023-03-27 | 2023-05-05 | 安徽星辰智跃科技有限责任公司 | 一种睡眠情绪活跃水平动态调节的方法、系统和装置 |
CN116058805A (zh) * | 2023-03-29 | 2023-05-05 | 安徽星辰智跃科技有限责任公司 | 一种睡眠记忆活跃水平动态调节的方法、系统和装置 |
Non-Patent Citations (1)
Title |
---|
基于去趋势波动分析的脑电信号睡眠分期;胡叶容;《生物医学工程学进展》;第36卷(第4期);第201-204页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116392085A (zh) | 2023-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116369866B (zh) | 基于小波变换的睡眠趋稳性量化及调节方法、系统和装置 | |
CN116392083B (zh) | 基于局部变换的睡眠趋稳性量化及调节方法、系统和装置 | |
CN116058804B (zh) | 一种睡眠情绪活跃水平动态调节的方法、系统和装置 | |
CN116058805B (zh) | 一种睡眠记忆活跃水平动态调节的方法、系统和装置 | |
CN115804573A (zh) | 一种睡眠深度量化及干预的方法、系统和装置 | |
CN116509337A (zh) | 基于局部分解的睡眠周期性检测及调节方法、系统和装置 | |
CN115881305B (zh) | 一种睡眠趋稳性检测量化及辅助干预的方法、系统和装置 | |
CN115910351B (zh) | 一种睡眠周期性检测量化及辅助干预的方法、系统和装置 | |
CN116392085B (zh) | 基于趋势分析的睡眠趋稳性量化及调节方法、系统和装置 | |
CN116392088B (zh) | 基于时频分析的睡眠趋稳性量化及调节方法、系统和装置 | |
CN116392087B (zh) | 基于模态分解的睡眠趋稳性量化及调节方法、系统和装置 | |
CN116525063B (zh) | 基于时频分析的睡眠周期性检测及调节方法、系统和装置 | |
CN116705247B (zh) | 基于局部分解的睡眠可持续性检测调节方法、系统和装置 | |
CN116504357A (zh) | 基于小波分析的睡眠周期性检测及调节方法、系统和装置 | |
CN115969330B (zh) | 一种睡眠情绪活跃水平检测量化的方法、系统和装置 | |
CN116013470B (zh) | 一种睡眠行为活跃水平动态调节的方法、系统和装置 | |
CN116092641B (zh) | 一种睡眠感官应激水平动态调节的方法、系统和装置 | |
CN116313029B (zh) | 一种数字针灸动态控制优化的方法、系统和装置 | |
CN115640827B (zh) | 对电刺激数据处理的智能闭环反馈网络方法及系统 | |
CN115966308A (zh) | 一种睡眠记忆活跃水平检测量化的方法、系统和装置 | |
CN116509336B (zh) | 基于波形分析的睡眠周期性检测及调节方法、系统和装置 | |
CN116509338B (zh) | 基于模态分析的睡眠周期性检测及调节方法、系统和装置 | |
CN116671875B (zh) | 基于小波变换的睡眠可持续性检测调节方法、系统和装置 | |
CN116509336A (zh) | 基于波形分析的睡眠周期性检测及调节方法、系统和装置 | |
CN116687356A (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 |