CN116779110A - 基于模态分解的睡眠可持续性检测调节方法、系统和装置 - Google Patents

基于模态分解的睡眠可持续性检测调节方法、系统和装置 Download PDF

Info

Publication number
CN116779110A
CN116779110A CN202310984123.2A CN202310984123A CN116779110A CN 116779110 A CN116779110 A CN 116779110A CN 202310984123 A CN202310984123 A CN 202310984123A CN 116779110 A CN116779110 A CN 116779110A
Authority
CN
China
Prior art keywords
sleep
sustainability
user
curve
adjustment
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202310984123.2A
Other languages
English (en)
Other versions
CN116779110B (zh
Inventor
何将
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Anhui Xingchen Zhiyue Technology Co ltd
Original Assignee
Anhui Xingchen Zhiyue Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Anhui Xingchen Zhiyue Technology Co ltd filed Critical Anhui Xingchen Zhiyue Technology Co ltd
Priority to CN202310984123.2A priority Critical patent/CN116779110B/zh
Publication of CN116779110A publication Critical patent/CN116779110A/zh
Application granted granted Critical
Publication of CN116779110B publication Critical patent/CN116779110B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/70—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mental therapies, e.g. psychological therapy or autogenous training
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
    • 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]
    • A61B5/346—Analysis of electrocardiograms
    • 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/369—Electroencephalography [EEG]
    • A61B5/372—Analysis of electroencephalograms
    • 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
    • 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/4809—Sleep detection, i.e. determining whether a subject is asleep or not
    • 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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • 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/7271—Specific aspects of physiological measurement analysis
    • A61B5/7275—Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/90—Details of database functions independent of the retrieved data types
    • G06F16/95—Retrieval from the web
    • G06F16/953—Querying, e.g. by the use of web search engines
    • G06F16/9535—Search customisation based on user profiles and personalisation
    • 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/20—Analysing
    • 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
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00—Machine learning
    • 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/70—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Data Mining & Analysis (AREA)
  • Veterinary Medicine (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Evolutionary Computation (AREA)
  • Psychiatry (AREA)
  • Evolutionary Biology (AREA)
  • Primary Health Care (AREA)
  • Physiology (AREA)
  • Epidemiology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Software Systems (AREA)
  • Cardiology (AREA)
  • Psychology (AREA)
  • Signal Processing (AREA)
  • Hospice & Palliative Care (AREA)
  • Developmental Disabilities (AREA)
  • Child & Adolescent Psychology (AREA)
  • Anesthesiology (AREA)
  • Computing Systems (AREA)
  • Pulmonology (AREA)

Abstract

本发明提供了基于模态分解的睡眠可持续性检测调节方法，包括：获取用户关键生理信息，根据睡眠可持续性动态检测调节基本方案，对用户睡眠生理信号进行采集和处理，得到睡眠状态特征曲线和睡眠时相曲线；对所述睡眠状态特征曲线进行去噪声处理和模态分解，根据特征曲线高频阈值提取高频模态分解信号，得到睡眠可持续性特征信号，计算睡眠可持续性指数和曲线；对睡眠可持续性进行趋势预测分析，生成睡眠可持续性动态调节策略，进行动态调节及效果评估；建立或更新用户个性睡眠可持续性数据库，动态迭代优化睡眠可持续性的动态检测调节过程方法策略，生成用户睡眠可持续性检测调节报告。本发明能够实现对用户睡眠可持续性的高效干预调节。

Description

基于模态分解的睡眠可持续性检测调节方法、系统和装置
技术领域
本发明涉及睡眠可持续性检测量化及辅助调节领域，特别涉及基于模态分解的睡眠可持续性检测调节方法、系统和装置。
背景技术
睡眠是人类最重要的生命生理过程之一，稳定且连续性好的睡眠是人们获得精力体力、身心健康的重要保障。现实生活中，由于衰老、睡眠环境、生理疾病及创伤、精神压力等众多因素，整夜睡眠过程中出现多次觉醒、长时间觉醒和睡眠状态急剧震荡等一系列睡眠连续性或可持续性差的事件，极大地影响了人们的睡眠质量和身心体验。长期的低质量的睡眠连续性或可持续性，对人们的身心健康、学习、工作和生活带来重大不良影响。
申请人提出的在先解决方案中国申请CN2023101940795提供了一种睡眠可持续性检测量化及辅助干预的方法、系统和装置，包括以下步骤：对用户睡眠过程的生理状态信号和环境状态信号进行采集监测、信号处理和特征分析，生成生理状态特征和环境状态特征；对所述生理状态特征进行睡眠状态分析、时序成分分析和可持续性量化分析，提取睡眠可持续性指数，生成睡眠可持续性量化日报；重复上述步骤，对用户睡眠过程进行连续监测和跟踪分析，评估睡眠环境对睡眠可持续性的影响，提取最佳睡眠可持续性环境方案并对睡眠环境进行动态优化调整，生成睡眠可持续性量化报告。该技术方案基于睡眠持续期状态特征曲线的时间序列分解，提取了睡眠可持续性指数，实现提出了睡眠可持续性指数作为睡眠连续性或可持续性量化的创新评价指标。但面对睡眠可持续的高效提取、动态检测和动态调节等场景需求，还存在可以进一步的提升空间，主要包括：首先，时间序列分解方法在多场景可持续特征提取中存在适应性局限而带来准确性的不稳定；其次，睡眠持续期状态特征曲线是基于阶梯式睡眠时相分期值平滑得到的，不能够对用户睡眠状态进行细致刻画和精确量化，也进一步限制了睡眠可持续性评价的准确性和灵敏度；最后，除了睡眠环境调节干预外，如何实现用户更快速高效的、长期稳定的且具备高个性化的睡眠可持续性指数检测量化和动态调节，并能够动态持续地提升检测准确性和调节高效性。
如何提取更细致的用户睡眠状态连续变化特征，更准确地量化评估睡眠可持续性；如何实现睡眠可持续性检测量化和动态调节的高效协同，并能够实现检测调节的过程动态优化，持续提高检测量化效率和干预调节效果，是目前国内外产品技术方案和实际应用场景中需要进一步优化或解决的问题。
发明内容
针对现有方法的以上缺陷及改进需求，本发明的目的在于提供一种基于模态分解的睡眠可持续性检测调节方法，通过对用户睡眠生理信号的采集处理和特征分析得到睡眠状态特征曲线和睡眠时相曲线，进一步通过模态分解从睡眠状态特征曲线中提取高频模态分解信号，计算得到睡眠可持续性指数；通过对用户睡眠可持续状态的趋势预测分析，生成睡眠可持续性动态调节策略并完成动态调节及效果评估；对检测调节过程进行统计分析，建立或更新用户个性睡眠可持续性数据库，动态迭代优化睡眠可持续性的动态检测调节过程方法策略，生成用户睡眠可持续性检测调节报告，并更新所述睡眠可持续性通用数据库，最终实现提高检测量化效率和干预调节效果，辅助用户获得更高的睡眠质量和睡眠连续性。本发明还提供了基于模态分解的睡眠可持续性检测调节系统，用于实现上述方法。本发明还提供了基于模态分解的睡眠可持续性检测调节装置，用于实现上述系统。
根据本发明的目的，本发明提出了一种基于模态分解的睡眠可持续性检测调节方法，包括以下步骤：
获取用户关键生理信息，从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选并初始化睡眠可持续性动态检测调节基本方案；
根据所述睡眠可持续性动态检测调节基本方案，对用户睡眠生理信号进行采集和处理，得到睡眠状态特征曲线和睡眠时相曲线；
对所述睡眠状态特征曲线进行去噪声处理和模态分解，根据特征曲线高频阈值提取高频模态分解信号，得到睡眠可持续性特征信号，计算睡眠可持续性指数和睡眠可持续性指数曲线；
对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析，生成睡眠可持续性动态调节策略，对用户睡眠过程进行动态调节及效果评估；
建立或更新所述用户个性睡眠可持续性数据库；
根据所述用户个性睡眠可持续性数据库，动态迭代优化睡眠可持续性的动态检测调节过程方法策略，生成用户睡眠可持续性检测调节报告，并更新所述睡眠可持续性通用数据库。
更优地，所述获取用户关键生理信息，从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选并初始化睡眠可持续性动态检测调节基本方案的具体步骤还包括：
获取所述用户关键生理信息，并从所述睡眠可持续性通用数据库或所述用户个性睡眠可持续性数据库中筛选与用户个性化生理匹配的所述睡眠可持续性动态检测调节基本方案；
根据所述睡眠可持续性动态检测调节基本方案，初始化用户睡眠可持续性的检测量化和动态调节的过程方法或过程策略。
更优地，所述用户关键生理信息至少包括性别、年龄、生理健康状态和心理精神状态。
更优地，所述睡眠可持续性通用数据库具体为由规模的不同健康状态人群的睡眠可持续性的检测量化及动态调节数据库，至少包括所述用户关键生理信息、睡眠可持续性检测调节周期、睡眠可持续性检测量化过程方法参数、睡眠可持续性动态调节过程方法参数或策略、睡眠可持续性指数曲线、睡眠可持续性调节效果曲线和睡眠时相曲线。
更优地，所述用户个性睡眠可持续性数据库至少包括所述用户关键生理信息、睡眠状态特征曲线、睡眠时相曲线、睡眠可持续性指数曲线、动态调节效果曲线、时相指数分布特征、时相效果分布特征、动态调节综合效果系数、模态分解方法、趋势预测分析方法、睡眠可持续性动态调节策略、所述睡眠可持续性动态检测调节基本方案。
更优地，所述睡眠可持续性动态检测调节基本方案至少包括睡眠可持续性检测量化过程方案、睡眠可持续性动态调节过程方案；若用户为新用户，所述睡眠可持续性动态检测调节基本方案从所述睡眠可持续性通用数据库中获取，反之则从所述用户个性睡眠可持续性数据库中获取。
更优地，所述根据所述睡眠可持续性动态检测调节基本方案，对用户睡眠生理信号进行采集和处理，得到睡眠状态特征曲线和睡眠时相曲线的具体步骤还包括：
对用户睡眠生理信号进行动态地监测采集和信号处理，得到睡眠生理状态数据；
对所述睡眠生理状态数据进行动态地特征分析和特征选择，生成所述睡眠状态特征曲线；
对所述睡眠生理状态数据进行动态地睡眠时相分析，生成所述睡眠时相曲线。
更优地，所述睡眠生理信号至少包括脑电信号、心电信号、呼吸信号中的任一项。
更优地，所述信号处理至少包括重采样、重参考、去伪迹、信号矫正、降噪、工频陷波、带通滤波、平滑处理和时帧分割；其中，所述时帧分割是指根据信号的采样率，以预设分帧步长对信号数据进行预设分帧时长窗口的连续滑动分割。
更优地，所述特征分析至少包括数值分析、包络分析、时频分析、熵分析、分形分析和复杂度分析中的任一项。
更优地，所述睡眠状态特征曲线具体由不同睡眠生理状态数据经所述特征分析得到不同睡眠生理状态特征进行线性组合融合而得到，用于描述用户在不同睡眠时期、不同睡眠时相、不同睡眠状态下的生理状态特征的连续变化曲线；其中，所述睡眠时期至少包括入睡前、睡眠中和睡眠后，所述睡眠时相至少包括清醒觉醒期、浅睡眠期、深睡眠期和快速眼动睡眠期。
更优地，所述睡眠时相曲线的生成方法具体为：
1)通过机器学习算法对规模睡眠用户样本的所述睡眠生理状态数据及其对应睡眠时相分期数据进行学习训练和数据建模，得到睡眠时相识别模型；
2)将当前用户的所述睡眠生理状态数据输入所述睡眠时相识别模型，得到所对应的睡眠时相分期并按时序生成所述睡眠时相曲线。
更优地，所述对所述睡眠状态特征曲线进行去噪声处理和模态分解，根据特征曲线高频阈值提取高频模态分解信号，得到睡眠可持续性特征信号，计算睡眠可持续性指数和睡眠可持续性指数曲线的具体步骤还包括：
对所述睡眠状态特征曲线进行动态地去噪声处理，得到目标睡眠状态特征曲线；
对所述目标睡眠状态特征曲线进行动态地模态分解，得到睡眠状态特征IMF分解信号集；
对所述睡眠状态特征IMF分解信号集进行动态地频谱分析，根据特征曲线高频阈值提取高频模态分解信号，生成所述睡眠可持续性特征信号；
根据所述睡眠时相曲线、所述目标睡眠状态特征曲线和所述睡眠可持续性特征信号，动态计算得到所述睡眠可持续性指数并生成所述睡眠可持续性指数曲线。
更优地，所述去噪声处理至少包括高斯滤波、均值滤波、傅里叶变换滤波、小波变换滤波和小波包变换滤波中的任一种方法，用于消除由于信号来自体表生理采集和分帧式的信号特征分析对睡眠状态刻画的不连续性，而在所述睡眠状态特征曲线中引入的高斯噪声或白噪声干扰。
更优地，所述模态分解至少包括经验模态分解、变分模态分解中的任一项。
更优地，所述经验模态分解的方法至少包括EMD、EEMD、CEEMD、 CEEMDAN、ICEEMDAN、ESMD中的任意一项。
更优地，所述变分模态分解的方法至少包括VMD、改进型VMD中的任意一项。
更优地，所述特征曲线高频阈值由所述睡眠状态特征曲线的睡眠生理状态特征线性组合融合的数据来源决定。
更优地，所述睡眠可持续性指数的计算方法，具体为：
1)获取当前时帧对应的睡眠时相分期、所述目标睡眠状态特征曲线和所述睡眠可持续性特征信号；
2)分别对所述目标睡眠状态特征曲线和所述睡眠可持续性特征信号做平方运算并求和，得到目标睡眠状态平方特征值和睡眠可持续性IMF平方特征值；
3)根据所述目标睡眠状态平方特征值和所述睡眠可持续性IMF平方特征值的相对变化量，计算得到睡眠可持续性节点特征系数；
4)根据所述睡眠可持续性节点特征系数、模态分解方法修正系数、用户个性修正系数和睡眠时相分期修正系数的均值和乘积，计算得到所述睡眠可持续性指数。
更优地，所述睡眠可持续性节点特征系数的一种计算公式，具体为：
；
其中，为所述睡眠可持续性节点特征系数，/>分别为所述目标睡眠状态平方特征值和所述睡眠可持续性IMF平方特征值。
更优地，所述睡眠可持续性指数的一种计算公式，具体为：
；
其中，为所述睡眠可持续性指数，/>分别为所述睡眠可持续性节点特征系数、所述模态分解方法修正系数、所述用户个性修正系数和所述睡眠时相分期修正系数，且/>，/>2。
更优地，所述睡眠可持续性指数曲线由按照时序拼接所述睡眠可持续性指数而生成得到。
更优地，所述对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析，生成睡眠可持续性动态调节策略，对用户睡眠过程进行动态调节及效果评估的具体步骤还包括：
对所述睡眠状态特征曲线进行动态地趋势预测分析，得到睡眠状态特征预测值；
对所述睡眠可持续性指数曲线进行动态地趋势预测分析，得到睡眠可持续性指数预测值；
根据所述睡眠状态特征预测值和所述睡眠可持续性指数预测值，按照睡眠可持续性检测调节周期和预设睡眠调节知识库，生成所述睡眠可持续性动态调节策略；
根据所述睡眠可持续性动态调节策略，连接并控制睡眠调节外围设备，对用户睡眠过程进行动态调节；
对调节效果进行动态跟踪评估，计算动态调节效果系数并生成动态调节效果曲线。
更优地，所述趋势预测分析的方法至少包括AR、MA、ARMA、ARIMA、SARIMA、VAR、深度学习中的任一项。
更优地，所述睡眠可持续性动态调节策略至少包括睡眠场景、睡眠时相、调节方式、执行方式、调节方法、调节强度、调节时点、持续时间、目标调节值和装置控制参数；其中，所述调节方式至少包括声乐刺激、超声刺激、光刺激、电刺激、磁刺激、温度刺激、湿度刺激、触觉刺激和浓度调控中的任一方式，所述执行方式至少包括离体式和接触式中的任一方式。
更优地，所述睡眠调节外围设备至少包括声乐刺激设备、超声刺激设备、光刺激设备、电刺激设备、磁刺激设备、温度刺激设备、湿度刺激设备、触觉刺激设备和浓度调控设备中的任一项，并由具体的所述调节方式决定。
更优地，所述动态调节效果系数的一种计算方式，具体为：
；
其中，为所述动态调节效果系数，/>为趋势预测分析方法修正系数且，/>为动态调节后的所述睡眠时相分期修正系数，分别为动态调节前、动态调节后的所述睡眠可持续性指数。
更优地，所述动态调节效果曲线由按照时序拼接所述动态调节效果系数而生成得到。
更优地，所述建立或更新用户个性睡眠可持续性数据库的具体步骤还包括：
根据所述睡眠时相曲线和所述睡眠可持续性指数曲线，计算不同睡眠时相下所述睡眠可持续性指数的分布特征，得到时相指数分布特征；
根据所述睡眠时相曲线和所述动态调节效果曲线，计算不同睡眠时相下所述动态调节效果系数的分布特征，得到时相效果分布特征；
计算所述睡眠可持续性指数曲线和所述动态调节效果曲线的相关性，得到动态调节综合效果系数；
动态收集检测调节过程数据和分析结果，初始化建立或持续动态更新所述用户个性睡眠可持续性数据库。
更优地，所述分布特征至少包括平均值、均方根、最大值、最小值、方差、标准差、变异系数、峰度和偏度中的任意一项。
更优地，所述相关性的计算方法至少包括相干性分析、皮尔逊相关分析、杰卡德相似分析、线性互信息分析、线性相关分析、欧氏距离分析、曼哈顿距离分析、切比雪夫距离分析中的任意一项。
更优地，所述根据所述用户个性睡眠可持续性数据库，动态迭代优化睡眠可持续性的动态检测调节过程方法策略，生成用户睡眠可持续性检测调节报告，并更新所述睡眠可持续性通用数据库的具体步骤还包括：
按照预设检测调节优化周期，根据所述用户个性睡眠可持续性数据库，动态迭代优化所述睡眠可持续性的所述动态检测调节过程方法策略，持续提高检测调节的效率效果；
按照预设报告生成周期，生成所述用户睡眠可持续性检测调节报告；
按照预设通用数据库更新周期，根据所述用户个性睡眠可持续性数据库和所述用户睡眠可持续性检测调节报告，将更新所述睡眠可持续性通用数据库的用户关键数据信息。
更优地，所述动态检测调节过程方法策略至少包括去高斯噪声处理方法参数、模态分解方法参数、所述特征曲线高频阈值、趋势预测分析方法参数、所述睡眠可持续性动态调节策略和动态调节效果评估方法参数。
更优地，所述用户睡眠可持续性检测调节报告至少包括所述用户关键生理信息、所述睡眠状态特征曲线、所述睡眠时相曲线、所述睡眠可持续性指数曲线、所述动态调节效果曲线、所述时相指数分布特征、所述时相效果分布特征、所述动态调节综合效果系数、睡眠可持续性检测调节总结。
根据本发明的目的，本发明提出了一种基于模态分解的睡眠可持续性检测调节系统，包括以下模块：
方案初始化模块，用于获取用户关键生理信息，从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选并初始化所述睡眠可持续性动态检测调节基本方案；
状态采集分析模块，用于根据所述睡眠可持续性动态检测调节基本方案，对用户睡眠生理信号进行采集和处理，得到睡眠状态特征曲线和睡眠时相曲线；
指数模态分析模块，用于对所述睡眠状态特征曲线进行去噪声处理和模态分解，根据特征曲线高频阈值提取高频模态分解信号，得到睡眠可持续性特征信号，计算睡眠可持续性指数和睡眠可持续性指数曲线；
睡眠动态调节模块，用于对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析，生成睡眠可持续性动态调节策略，对用户睡眠过程进行动态调节及效果评估；
数据统计更新模块，用于建立或更新用户个性睡眠可持续性数据库；
检测调节优化模块，用于根据所述用户个性睡眠可持续性数据库，动态迭代优化睡眠可持续性的动态检测调节过程方法策略，生成用户睡眠可持续性检测调节报告，并更新所述睡眠可持续性通用数据库；
数据运行管理模块，用于对所述系统的所有数据进行可视化管理、统一存储和运行管理。
更优地，所述方案初始化模块还包括以下功能单元：
基础方案选择单元，用于获取所述用户关键生理信息，并从所述睡眠可持续性通用数据库或所述用户个性睡眠可持续性数据库中筛选与用户个性化生理匹配的所述睡眠可持续性动态检测调节基本方案；
检测调节初始化单元，用于根据所述睡眠可持续性动态检测调节基本方案，初始化用户睡眠可持续性的检测量化和动态调节的过程方法或过程策略。
更优地，所述状态采集分析模块还包括以下功能单元：
状态采集处理单元，用于对用户睡眠生理信号进行动态地监测采集和信号处理，得到睡眠生理状态数据；
状态特征提取单元，用于对所述睡眠生理状态数据进行动态地特征分析和特征选择，生成所述睡眠状态特征曲线；
睡眠时相分析单元，用于对所述睡眠生理状态数据进行动态地睡眠时相分析，生成所述睡眠时相曲线。
更优地，所述指数模态分析模块还包括以下功能单元：
高斯噪声处理单元，用于对所述睡眠状态特征曲线进行动态地去噪声处理，得到目标睡眠状态特征曲线；
信号模态分解单元，用于对所述目标睡眠状态特征曲线进行动态地模态分解，得到睡眠状态特征IMF分解信号集；
分解信号提取单元，用于对所述睡眠状态特征IMF分解信号集进行频谱分析，根据特征曲线高频阈值提取高频模态分解信号，生成睡眠可持续性特征信号；
睡眠指数计算单元，用于根据所述睡眠时相曲线、所述目标睡眠状态特征曲线和所述睡眠可持续性特征信号，动态计算得到所述睡眠可持续性指数并生成所述睡眠可持续性指数曲线。
更优地，所述睡眠动态调节模块还包括以下功能单元：
状态趋势预测单元，用于对所述睡眠状态特征曲线进行动态地趋势预测分析，得到睡眠状态特征预测值；
指数趋势预测单元，用于对所述睡眠可持续性指数曲线进行动态地趋势预测分析，得到睡眠可持续性指数预测值；
动态策略生成单元，用于根据所述睡眠状态特征预测值和所述睡眠可持续性指数预测值，按照睡眠可持续性检测调节周期和预设睡眠调节知识库，生成所述睡眠可持续性动态调节策略；
睡眠动态调节单元，用于根据所述睡眠可持续性动态调节策略，连接并控制睡眠调节外围设备，对用户睡眠过程进行动态调节；
效果动态评估单元，用于对调节效果进行动态跟踪评估，计算动态调节效果系数并生成动态调节效果曲线。
更优地，所述数据统计更新模块还包括以下功能单元：
时相指数分析单元，用于根据所述睡眠时相曲线和所述睡眠可持续性指数曲线，计算不同睡眠时相下所述睡眠可持续性指数的分布特征，得到时相指数分布特征；
时相效果分析单元，用于根据所述睡眠时相曲线和所述动态调节效果曲线，计算不同睡眠时相下所述动态调节效果系数的分布特征，得到时相效果分布特征；
综合效果评价单元，用于计算所述睡眠可持续性指数曲线和所述动态调节效果曲线的相关性，得到动态调节综合效果系数；
个性数据库更新单元，用于动态收集检测调节过程数据和分析结果，初始化建立或持续动态更新所述用户个性睡眠可持续性数据库。
更优地，所述检测调节优化模块还包括以下功能单元：
过程动态优化单元，用于按照预设检测调节优化周期，根据所述用户个性睡眠可持续性数据库，动态迭代优化所述睡眠可持续性的所述动态检测调节过程方法策略，持续提高检测调节的效率效果；
用户报告生成单元，用于按照预设报告生成周期，生成所述用户睡眠可持续性检测调节报告；
用户报告管理单元，用于对所述用户睡眠可持续性检测调节报告的格式输出、展现形式进行统一管理；
通用数据库更新单元，用于按照预设通用数据库更新周期，根据所述用户个性睡眠可持续性数据库和所述用户睡眠可持续性检测调节报告，将更新所述睡眠可持续性通用数据库的用户关键数据信息。
更优地，所述数据运行管理模块还包括以下功能单元：
用户信息管理单元，用于用户基本信息的登记输入、编辑、查询、输出和删除；
数据可视化管理单元，用于对所述系统中所有数据的可视化展现管理；
数据存储管理单元，用于对所述系统中所有数据的统一存储管理；
数据运营管理单元，用于对所述系统中所有数据的备份、迁移和导出。
根据本发明的目的，本发明提出了一种基于模态分解的睡眠可持续性检测调节装置，包括以下模组：
方案初始化模组，用于获取用户关键生理信息，从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选并初始化所述睡眠可持续性动态检测调节基本方案；
状态采集分析模组，用于根据所述睡眠可持续性动态检测调节基本方案，对用户睡眠生理信号进行采集和处理，得到睡眠状态特征曲线和睡眠时相曲线；
指数模态分析模组，用于对所述睡眠状态特征曲线进行去噪声处理和模态分解，根据特征曲线高频阈值提取高频模态分解信号，得到睡眠可持续性特征信号，计算睡眠可持续性指数和睡眠可持续性指数曲线；
睡眠动态调节模组，用于对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析，生成睡眠可持续性动态调节策略，对用户睡眠过程进行动态调节及效果评估；
数据统计更新模组，用于建立或更新用户个性睡眠可持续性数据库；
检测调节优化模组，用于根据所述用户个性睡眠可持续性数据库，动态迭代优化睡眠可持续性的动态检测调节过程方法策略，生成用户睡眠可持续性检测调节报告，并更新所述睡眠可持续性通用数据库；
数据可视化模组，用于对所述装置中所有过程数据和结果数据的统一可视化展示管理；
数据管理中心模组，用于对所述装置中所有过程数据和结果数据的统一存储和数据运营管理。
本发明在申请人的在先研究基础上进一步优化睡眠可持续性指数的具体量化设计和计算方式，进一步优化睡眠状态特征曲线的提取方式，进一步将去噪声处理、经验模态分解和变分模态分解方法应用于睡眠可持续性特征的提取，使得对用户睡眠可持续性的检测量化更加全面细致、适应度广、准确度高和灵敏度高；还进一步提供了动态调节效果系数的计算方案，进一步提供了用户个性睡眠可持续性数据库和睡眠可持续性通用数据库的建立、更新和应用机制，进一步提供了检测量化到动态调节的反向反馈应用框架，使得对于检测量化和动态调节过程的协同控制有了有力依据，实现检测调节的过程动态优化，持续提高检测量化效率和干预调节效果。本发明能够提供更加科学高效的睡眠可持续性检测量化及动态调节的实施方法和落地方案，能够赋能相关睡眠量化或调节的产品及服务，满足不同用户场景需求，辅助用户睡眠。
本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明技术方案的进一步理解，并且构成说明书的一部分，与本申请的实施例一起用于解释本发明的技术方案，并不构成对本发明技术方案的限制。
图1是本发明一实施例所提供的一种基于模态分解的睡眠可持续性检测调节方法的流程步骤示意图；
图2是本发明一实施例所提供的一种基于模态分解的睡眠可持续性检测调节系统的模块组成示意图；
图3是本发明一实施例所提供的一种基于模态分解的睡眠可持续性检测调节装置的模组构成示意图。
具体实施方式
为了更清楚地说明本发明的目的和技术方案，下面将结合本发明申请实施例中的附图，对本发明进行进一步介绍说明。显而易见地，下面描述的实施例仅是本发明一部分实施例，而不是全部的实施例。在没有创造性劳动前提下，本领域普通技术人员基于本发明的实施例所得到的其他实施例，都应属于本发明的保护范围。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。
申请人发现，人类睡眠是个非平稳时序过程，而睡眠过程中多次觉醒睡眠状态突变震荡是非平稳的重要因素或诱因，如何跟踪多变状态、意外状态和突变状态是睡眠可持续性或连续性检测调节的一个基本挑战。同时由于头皮脑电EEG、心电ECG等体表生理采集技术的原理限制及对应的传统信号特征分析提取过程，睡眠状态量化上存在非常多的计算误差和引入噪声，尤其是面对睡眠可持续性的评价上。本技术方案通过敏感特征的选择、去噪声处理和模态分解等方法，能够提取细致和精确的睡眠状态特征曲线，完成睡眠可持续性的高鲁棒性的量化设计；进而完成睡眠可持续性的动态调节和效果评估，进一步实现检测调节的反向动态优化，使得整体检测调节方案具备更好的场景适应性和技术应用效果。
结合图1所示，本发明实施例提供的一种基于模态分解的睡眠可持续性检测调节方法，包括以下步骤：
P100：获取用户关键生理信息，从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选并初始化所述睡眠可持续性动态检测调节基本方案。
第一步、获取用户关键生理信息，并从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选与用户个性化生理匹配的睡眠可持续性动态检测调节基本方案。
本实施例中，用户关键生理信息至少包括性别、年龄、生理健康状态和心理精神状态。
本实施例中，睡眠可持续性通用数据库具体为由规模的不同健康状态人群的睡眠可持续性的检测量化及动态调节数据库，至少包括用户关键生理信息、睡眠可持续性检测调节周期、睡眠可持续性检测量化过程方法参数、睡眠可持续性动态调节过程方法参数或策略、睡眠可持续性指数曲线、睡眠可持续性调节效果曲线和睡眠时相曲线。
本实施例中，用户个性睡眠可持续性数据库至少包括用户关键生理信息、睡眠状态特征曲线、睡眠时相曲线、睡眠可持续性指数曲线、动态调节效果曲线、时相指数分布特征、时相效果分布特征、动态调节综合效果系数、模态分解方法、趋势预测分析方法、睡眠可持续性动态调节策略、睡眠可持续性动态检测调节基本方案。
本实施例中，睡眠可持续性动态检测调节基本方案至少包括睡眠可持续性检测量化过程方案、睡眠可持续性动态调节过程方案；若用户为新用户，睡眠可持续性动态检测调节基本方案从睡眠可持续性通用数据库中获取，反之则从用户个性睡眠可持续性数据库中获取。
在实际应用场景中，首先确定当前用户是新用户还是老用户，来选择睡眠可持续性动态检测调节基本方案的获取来源；当前用户的当前状态，包括生理健康状态和心理精神状态，是睡眠可持续性动态检测调节基本方案选择的重要考量。
第二步、根据睡眠可持续性动态检测调节基本方案，初始化用户睡眠可持续性的检测量化和动态调节的过程方法或过程策略。
根据睡眠可持续性动态检测调节基本方案，来初始化检测量化和动态调节的初始使用方法、方法参数、初始调节策略、检测调节周期等等。
P200：根据所述睡眠可持续性动态检测调节基本方案，对用户睡眠生理信号进行采集和处理，得到睡眠状态特征曲线和睡眠时相曲线。
第一步、对用户睡眠生理信号进行动态地监测采集和信号处理，得到睡眠生理状态数据。
本实施例中，睡眠生理信号至少包括脑电信号、心电信号、呼吸信号中的任一项。信号处理至少包括重采样、重参考、去伪迹、信号矫正、降噪、工频陷波、带通滤波、平滑处理和时帧分割；其中，时帧分割是指根据信号的采样率，以预设分帧步长对信号数据进行预设分帧时长窗口的连续滑动分割。
本实施例中，以脑电信号和心电信号作为用户睡眠生理信号，来陈述本技术方案的具体实施过程。通过脑电图机对用户睡眠脑电进行动态地采集监测，采样率为1024Hz，记录电极为F3、F4、C3、C4、T3和T4，参考电极M1、M2；对脑电信号进行动态地统一信号处理，包括以M1和M2进行左右交叉重参考、去伪迹、信号矫正、小波降噪、50Hz倍频工频陷波、1.0-80Hz带通滤波，得到纯净脑电信号。通过便携式单导心电仪对用户心电信号进行采集监测，采集位置在左胸上方，采样率1024Hz；对心电信号进行动态地统一信号处理，包括去伪迹、信号处理、小波降噪、0.5-40Hz带通滤波，得到纯净心电信号。以预设时间窗口长度20s和预设时间平移步长10s对纯净脑电信号和纯净心电信号进行连续滑动分割，得到睡眠生理状态数据。
本实施例中，预设时间窗口长度20s也是睡眠可持续性检测调节周期。
第二步、对睡眠生理状态数据进行动态地特征分析和特征选择，生成睡眠状态特征曲线。
本实施例中，特征分析至少包括数值分析、包络分析、时频分析、熵分析、分形分析和复杂度分析中的任一项。睡眠状态特征曲线具体由不同睡眠生理状态数据经特征分析得到不同睡眠生理状态特征进行线性组合融合而得到，用于描述用户在不同睡眠时期、不同睡眠时相、不同睡眠状态下的生理状态特征的连续变化曲线；其中，睡眠时期至少包括入睡前、睡眠中和睡眠后，睡眠时相至少包括清醒觉醒期、浅睡眠期、深睡眠期和快速眼动睡眠期。
本实施例中，以对睡眠可持续性敏感度较高的指标特征来生成睡眠状态特征曲线。首先，逐帧对睡眠生理状态数据进行特征提取：脑电数据进行时频分析（频带功率占比）、熵分析（SVD熵）和复杂度分析（LZC指数）；心电数据进行数值分析，提取用户心率变化特征（心率均值、心率变异系数）。其次，通过特征选择，将F4-M1通道的δ节律（1-4Hz）频带功率占比、取负后再归一化的SVD熵、取负后再归一化的LZC指数，以及取负后再归一化的心率均值进行直接加和，得到睡眠状态特征曲线。
第三步、对睡眠生理状态数据进行动态地睡眠时相分析，生成睡眠时相曲线。
本实施例中，睡眠时相曲线的生成方法具体为：
1)通过机器学习算法对规模睡眠用户样本的睡眠生理状态数据及其对应睡眠时相分期数据进行学习训练和数据建模，得到睡眠时相识别模型；
2)将当前用户的睡眠生理状态数据输入睡眠时相识别模型，得到所对应的睡眠时相分期并按时序生成睡眠时相曲线。
P300：对所述睡眠状态特征曲线进行去噪声处理和模态分解，根据特征曲线高频阈值提取高频模态分解信号，得到睡眠可持续性特征信号，计算睡眠可持续性指数和睡眠可持续性指数曲线。
第一步、对睡眠状态特征曲线进行动态地去噪声处理，得到目标睡眠状态特征曲线。
本实施例中，去噪声处理至少包括高斯滤波、均值滤波、傅里叶变换滤波、小波变换滤波和小波包变换滤波中的任一种方法，用于消除由于信号来自体表生理采集和分帧式的信号特征分析对睡眠状态刻画的不连续性，而在所述睡眠状态特征曲线中引入的高斯噪声或白噪声干扰。
本实施例中，选择高斯滤波来完成去噪声处理。
第二步、对目标睡眠状态特征曲线进行动态地模态分解，得到睡眠状态特征IMF分解信号集。
本实施例中，模态分解至少包括经验模态分解、变分模态分解中的任一项。经验模态分解的方法至少包括EMD、EEMD、CEEMD、 CEEMDAN、ICEEMDAN、ESMD中的任意一项；其中，ESMD具体为极点对称模态分解方法，其借鉴了EMD思想，将外部包络线插值改为内部极点对称插值，借用“最小二乘”思想来优化最后剩余模态使其成为整个数据的“自适应全局均线”，以确定最佳筛选次数。变分模态分解的方法至少包括VMD、改进型VMD中的任意一项；其中，改进型VMD具体是为适应不同场景的信号分析和信号分解需求进行信号处理算法与VMD的有机组合。
本实施例中，选择标准变分模态分解（VMD）对目标睡眠状态特征曲线进行变分模态分解，得到睡眠状态特征IMF分解信号集。相比于EMD经验模态分解，VMD变分模态分解具有较好抗噪能力，克服经验模态分解的端点效应、模态分量/频率混叠等问题，能够降低复杂度高和非线性强的时间序列非平稳性，将时间序列数据分解为一系列具有有限带宽的本征模态函数（IMF），可自适应更新各IMF的最优中心频率和带宽，适用于非平稳性的序列。其主要步骤如下：
1)将原始信号进行多次低通滤波，得到多个频带信号。
2)对每个频带信号进行变分推断，得到该频带信号的局部振动模式。
3)将所有频带信号对应的局部振动模式相加，得到原始信号的 VMD 分解。
本实施例中，通过VMD提取10个IMF分解信号，生成睡眠状态特征IMF分解信号集。
第三步、对睡眠状态特征IMF分解信号集进行动态地频谱分析，根据特征曲线高频阈值提取高频模态分解信号，生成睡眠可持续性特征信号。
本实施例中，特征曲线高频阈值由睡眠状态特征曲线的睡眠生理状态特征线性组合融合的数据来源决定。
本实施例中，以δ节律频带功率占比、取负后再归一化的SVD熵、取负后再归一化的LZC指数，以及取负后再归一化的心率均值来组合生成睡眠状态特征曲线，且分帧处理的参数为预设时间窗口长度20s和预设时间平移步长10s，所以0.013Hz高通截止频率作为特征曲线高频阈值。最后，通过Welch功率谱密度估计对睡眠状态特征IMF分解信号集进行频谱分析，筛选中心频率大于0.013Hz的IMF分解信号，加和生成睡眠可持续性特征信号。
第四步、根据睡眠时相曲线、目标睡眠状态特征曲线和睡眠可持续性特征信号，动态计算得到睡眠可持续性指数并生成睡眠可持续性指数曲线。
本实施例中，睡眠可持续性指数的计算方法，具体为：
1)获取当前时帧对应的睡眠时相分期、目标睡眠状态特征曲线和睡眠可持续性特征信号；
2)分别对目标睡眠状态特征曲线和睡眠可持续性特征信号做平方运算并求和，得到目标睡眠状态平方特征值和睡眠可持续性IMF平方特征值；
3)根据目标睡眠状态平方特征值和睡眠可持续性IMF平方特征值的相对变化量，计算得到睡眠可持续性节点特征系数；
4)根据睡眠可持续性节点特征系数、模态分解方法修正系数、用户个性修正系数和睡眠时相分期修正系数的均值和乘积，计算得到睡眠可持续性指数。
本实施例中，睡眠可持续性节点特征系数的一种计算公式，具体为：
；
其中，为睡眠可持续性节点特征系数，/>分别为目标睡眠状态平方特征值和睡眠可持续性IMF平方特征值。
本实施例中，睡眠可持续性指数的一种计算公式，具体为：
；
其中，为睡眠可持续性指数，/>分别为睡眠可持续性节点特征系数、模态分解方法修正系数、用户个性修正系数和睡眠时相分期修正系数，且/>，/>2。
本实施例中，睡眠可持续性指数曲线由按照时序拼接睡眠可持续性指数而生成得到。
在实际应用场景中，通常清醒觉醒期、浅睡眠期、深睡眠期和快速眼动睡眠期的修正系数，依次为 1.0、1.6、2.0、1.4；用户个性修正系数/>主要跟用户年龄和性别相关，通常女性比男性要小，年龄越大系数越小；模态分解方法修正系数/>则可以根据睡眠状态特征曲线的生成方式、去噪声处理的方法参数等实际情况设置不同方法的系数。更重要的是，在后续的动态迭代优化睡眠可持续性的动态检测调节过程方法策略中可以进行动态优化调整。
P400：对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析，生成睡眠可持续性动态调节策略，对用户睡眠过程进行动态调节及效果评估。
第一步、对睡眠状态特征曲线进行动态地趋势预测分析，得到睡眠状态特征预测值。
本实施例中，趋势预测分析的方法至少包括AR、MA、ARMA、ARIMA、SARIMA、VAR、深度学习中的任一项。
本实施例中，用AR方法来完成对睡眠状态特征曲线的预测。
第二步、对睡眠可持续性指数曲线进行动态地趋势预测分析，得到睡眠可持续性指数预测值。
本实施例中，用AR方法来完成对睡眠可持续性指数曲线的预测。
第三步、根据睡眠状态特征预测值和睡眠可持续性指数预测值，按照睡眠可持续性检测调节周期和预设睡眠调节知识库，生成睡眠可持续性动态调节策略。
本实施例中，睡眠可持续性动态调节策略至少包括睡眠场景、睡眠时相、调节方式、执行方式、调节方法、调节强度、调节时点、持续时间、目标调节值和装置控制参数；其中，调节方式至少包括声乐刺激、超声刺激、光刺激、电刺激、磁刺激、温度刺激、湿度刺激、触觉刺激和浓度调控中的任一方式，执行方式至少包括离体式和接触式中的任一方式。
在实际应用场景中，需要根据用户具体情况来制定不同睡眠可持续性动态调节策略。如选择离体式的声乐刺激、光刺激、温度刺激、湿度刺激、浓度调控等方式，能够减少对用户睡眠过程的干扰；但通常情况下，接触式的超声刺激、电刺激、磁刺激、触觉刺激等方式的调节效率更高、效果也更好。
第四步、根据睡眠可持续性动态调节策略，连接并控制睡眠调节外围设备，对用户睡眠过程进行动态调节。
本实施例中，睡眠调节外围设备至少包括声乐刺激设备、超声刺激设备、光刺激设备、电刺激设备、磁刺激设备、温度刺激设备、湿度刺激设备、触觉刺激设备和浓度调控设备中的任一项，并由具体的调节方式决定。
第五步、对调节效果进行动态跟踪评估，计算动态调节效果系数并生成动态调节效果曲线。
本实施例中，动态调节效果系数的一种计算方式，具体为：
；
其中，为动态调节效果系数，/>为趋势预测分析方法修正系数且，/>为动态调节后的睡眠时相分期修正系数，/>分别为动态调节前、动态调节后的睡眠可持续性指数。
本实施例中，动态调节效果曲线由按照时序拼接动态调节效果系数而生成得到。
在实际应用场景中，通常趋势预测分析方法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还包括以下功能单元：
高斯噪声处理单元，用于对睡眠状态特征曲线进行动态地去噪声处理，得到目标睡眠状态特征曲线；
信号模态分解单元，用于对目标睡眠状态特征曲线进行动态地模态分解，得到睡眠状态特征IMF分解信号集；
分解信号提取单元，用于对睡眠状态特征IMF分解信号集进行频谱分析，根据特征曲线高频阈值提取高频模态分解信号，生成睡眠可持续性特征信号；
睡眠指数计算单元，用于根据睡眠时相曲线、目标睡眠状态特征曲线和睡眠可持续性特征信号，动态计算得到睡眠可持续性指数并生成睡眠可持续性指数曲线。
本实施例中，睡眠动态调节模块S400还包括以下功能单元：
状态趋势预测单元，用于对睡眠状态特征曲线进行动态地趋势预测分析，得到睡眠状态特征预测值；
指数趋势预测单元，用于对睡眠可持续性指数曲线进行动态地趋势预测分析，得到睡眠可持续性指数预测值；
动态策略生成单元，用于根据睡眠状态特征预测值和睡眠可持续性指数预测值，按照睡眠可持续性检测调节周期和预设睡眠调节知识库，生成睡眠可持续性动态调节策略；
睡眠动态调节单元，用于根据睡眠可持续性动态调节策略，连接并控制睡眠调节外围设备，对用户睡眠过程进行动态调节；
效果动态评估单元，用于对调节效果进行动态跟踪评估，计算动态调节效果系数并生成动态调节效果曲线。
本实施例中，数据统计更新模块S500还包括以下功能单元：
时相指数分析单元，用于根据睡眠时相曲线和睡眠可持续性指数曲线，计算不同睡眠时相下睡眠可持续性指数的分布特征，得到时相指数分布特征；
时相效果分析单元，用于根据睡眠时相曲线和动态调节效果曲线，计算不同睡眠时相下动态调节效果系数的分布特征，得到时相效果分布特征；
综合效果评价单元，用于计算睡眠可持续性指数曲线和动态调节效果曲线的相关性，得到动态调节综合效果系数；
个性数据库更新单元，用于动态收集检测调节过程数据和分析结果，初始化建立或持续动态更新用户个性睡眠可持续性数据库。
本实施例中，检测调节优化模块S600还包括以下功能单元：
过程动态优化单元，用于按照预设检测调节优化周期，根据用户个性睡眠可持续性数据库，动态迭代优化睡眠可持续性的动态检测调节过程方法策略，持续提高检测调节的效率效果；
用户报告生成单元，用于按照预设报告生成周期，生成用户睡眠可持续性检测调节报告；
用户报告管理单元，用于对用户睡眠可持续性检测调节报告的格式输出、展现形式进行统一管理；
通用数据库更新单元，用于按照预设通用数据库更新周期，根据用户个性睡眠可持续性数据库和用户睡眠可持续性检测调节报告，将更新睡眠可持续性通用数据库的用户关键数据信息。
本实施例中，数据运行管理模块S700还包括以下功能单元：
用户信息管理单元，用于用户基本信息的登记输入、编辑、查询、输出和删除；
数据可视化管理单元，用于对系统中所有数据的可视化展现管理；
数据存储管理单元，用于对系统中所有数据的统一存储管理；
数据运营管理单元，用于对系统中所有数据的备份、迁移和导出。
结合图3所示，本发明实施例提供的一种基于模态分解的睡眠可持续性检测调节装置，包括以下模组：
方案初始化模组M100，用于获取用户关键生理信息，从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选并初始化睡眠可持续性动态检测调节基本方案；
状态采集分析模组M200，用于根据睡眠可持续性动态检测调节基本方案，对用户睡眠生理信号进行采集和处理，得到睡眠状态特征曲线和睡眠时相曲线；
指数模态分析模组M300，用于对睡眠状态特征曲线进行去噪声处理和模态分解，根据特征曲线高频阈值提取高频模态分解信号，得到睡眠可持续性特征信号，计算睡眠可持续性指数和睡眠可持续性指数曲线；
睡眠动态调节模组M400，用于对睡眠状态特征曲线和睡眠可持续性指数曲线进行趋势预测分析，生成睡眠可持续性动态调节策略，对用户睡眠过程进行动态调节及效果评估；
数据统计更新模组M500，用于建立或更新用户个性睡眠可持续性数据库；
检测调节优化模组M600，用于根据用户个性睡眠可持续性数据库，动态迭代优化睡眠可持续性的动态检测调节过程方法策略，生成用户睡眠可持续性检测调节报告，并更新睡眠可持续性通用数据库；
数据可视化模组M700，用于对装置中所有过程数据和结果数据的统一可视化展示管理；
数据管理中心模组M800，用于对装置中所有过程数据和结果数据的统一存储和数据运营管理。
所述装置被构造以用于对应执行图1的方法中的各个步骤，在此不再赘述。
本发明还提供了可编程的各类处理器（FPGA、ASIC或其他集成电路），所述处理器用于运行程序，其中，所述程序运行时执行上述实施例中的步骤。
本发明还提供了对应的计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述存储器执行所述程序时实现上述实施例中的步骤。
虽然本发明所揭露的实施方式如上，但所述的 内容仅为便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属领域内的技术人员，在不脱离本发明所揭露的精神和原则的前提下，可以在实施的形式及细节上进行任何的修改与变化、等同替换等，这些都属于本发明的保护范围。因此，本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。

Claims (45)

1.一种基于模态分解的睡眠可持续性检测调节方法，其特征在于，包括以下步骤：
获取用户关键生理信息，根据用户关键生理信息筛选并初始化睡眠可持续性动态检测调节基本方案；
根据所述睡眠可持续性动态检测调节基本方案，对用户睡眠生理信号进行采集和处理，得到睡眠状态特征曲线和睡眠时相曲线；
对所述睡眠状态特征曲线进行去噪声处理和模态分解，根据特征曲线高频阈值提取高频模态分解信号，得到睡眠可持续性特征信号，计算睡眠可持续性指数和睡眠可持续性指数曲线；
对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析，生成睡眠可持续性动态调节策略，对用户睡眠过程进行动态调节及效果评估；
所述睡眠可持续性指数的计算方法，具体为：
1)获取当前时帧对应的睡眠时相分期、目标睡眠状态特征曲线和所述睡眠可持续性特征信号；
2)分别对所述目标睡眠状态特征曲线和所述睡眠可持续性特征信号做平方运算并求和，得到目标睡眠状态平方特征值和睡眠可持续性IMF平方特征值；
3)根据所述目标睡眠状态平方特征值和所述睡眠可持续性IMF平方特征值的相对变化量，得到睡眠可持续性节点特征系数；
4)根据所述睡眠可持续性节点特征系数、模态分解方法修正系数、用户个性修正系数和睡眠时相分期修正系数，计算得到所述睡眠可持续性指数；其中，所述可持续性节点特征系数、所述睡眠时相分期修正系数均与所述睡眠可持续性指数正相关，所述模态分解方法修正系数和所述用户个性修正系数用于修正所述睡眠可持续性指数。
2.如权利要求1所述的方法，其特征在于,所述获取用户关键生理信息，根据用户关键生理信息筛选并初始化睡眠可持续性动态检测调节基本方案的具体步骤还包括：
获取所述用户关键生理信息，并从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选与用户个性化生理信息匹配的所述睡眠可持续性动态检测调节基本方案；
根据所述睡眠可持续性动态检测调节基本方案，初始化用户睡眠可持续性的检测量化和动态调节的过程方法或过程策略。
3.如权利要求2所述的方法，其特征在于：所述用户关键生理信息至少包括性别、年龄、生理健康状态和心理精神状态。
4.如权利要求2或3所述的方法，其特征在于：所述睡眠可持续性通用数据库具体为由规模的不同健康状态人群的睡眠可持续性的检测量化及动态调节数据库，至少包括所述用户关键生理信息、睡眠可持续性检测调节周期、睡眠可持续性检测量化过程方法参数、睡眠可持续性动态调节过程方法参数或策略、睡眠可持续性指数曲线、睡眠可持续性调节效果曲线和睡眠时相曲线。
5.如权利要求2或3所述的方法，其特征在于：所述用户个性睡眠可持续性数据库至少包括所述用户关键生理信息、睡眠状态特征曲线、睡眠时相曲线、睡眠可持续性指数曲线、动态调节效果曲线、时相指数分布特征、时相效果分布特征、动态调节综合效果系数、模态分解方法、趋势预测分析方法、睡眠可持续性动态调节策略、所述睡眠可持续性动态检测调节基本方案。
6.如权利要求2或3所述的方法，其特征在于：所述睡眠可持续性动态检测调节基本方案至少包括睡眠可持续性检测量化过程方案、睡眠可持续性动态调节过程方案；若用户为新用户，所述睡眠可持续性动态检测调节基本方案从所述睡眠可持续性通用数据库中获取，反之则从所述用户个性睡眠可持续性数据库中获取。
7.如权利要求1或2所述的方法，其特征在于，所述根据所述睡眠可持续性动态检测调节基本方案，对用户睡眠生理信号进行采集和处理，得到睡眠状态特征曲线和睡眠时相曲线的具体步骤还包括：
对用户睡眠生理信号进行动态地监测采集和信号处理，得到睡眠生理状态数据；
对所述睡眠生理状态数据进行动态地特征分析和特征选择，生成所述睡眠状态特征曲线；
对所述睡眠生理状态数据进行动态地睡眠时相分析，生成所述睡眠时相曲线。
8.如权利要求7所述的方法，其特征在于，所述睡眠生理信号至少包括脑电信号、心电信号、呼吸信号中的任一项。
9.如权利要求8所述的方法，其特征在于：所述信号处理至少包括重采样、重参考、去伪迹、信号矫正、降噪、工频陷波、带通滤波、平滑处理和时帧分割；其中，所述时帧分割是指根据信号的采样率，以预设分帧步长对信号数据进行预设分帧时长窗口的连续滑动分割。
10.如权利要求7所述的方法，其特征在于：所述特征分析至少包括数值分析、包络分析、时频分析、熵分析、分形分析和复杂度分析中的任一项。
11.如权利要求10所述的方法，其特征在于：所述睡眠状态特征曲线具体由不同睡眠生理状态数据经所述特征分析得到不同睡眠生理状态特征进行线性组合融合而得到，用于描述用户在不同睡眠时期、不同睡眠时相、不同睡眠状态下的生理状态特征的连续变化曲线；其中，所述睡眠时期至少包括入睡前、睡眠中和睡眠后，所述睡眠时相至少包括清醒觉醒期、浅睡眠期、深睡眠期和快速眼动睡眠期。
12.如权利要求8或11所述的方法，其特征在于，所述睡眠时相曲线的生成方法具体为：
1)通过机器学习算法对规模睡眠用户样本的所述睡眠生理状态数据及其对应睡眠时相分期数据进行学习训练和数据建模，得到睡眠时相识别模型；
2)将当前用户的所述睡眠生理状态数据输入所述睡眠时相识别模型，得到所对应的睡眠时相分期并按时序生成所述睡眠时相曲线。
13.如权利要求1或2所述的方法，其特征在于，所述对所述睡眠状态特征曲线进行去噪声处理和模态分解，根据特征曲线高频阈值提取高频模态分解信号，得到睡眠可持续性特征信号，计算睡眠可持续性指数和睡眠可持续性指数曲线的具体步骤还包括：
对所述睡眠状态特征曲线进行动态地去噪声处理，得到目标睡眠状态特征曲线；
对所述目标睡眠状态特征曲线进行动态地模态分解，得到睡眠状态特征IMF分解信号集；
对所述睡眠状态特征IMF分解信号集进行动态地频谱分析，根据特征曲线高频阈值提取高频模态分解信号，生成所述睡眠可持续性特征信号；
根据所述睡眠时相曲线、所述目标睡眠状态特征曲线和所述睡眠可持续性特征信号，动态计算得到所述睡眠可持续性指数并生成所述睡眠可持续性指数曲线。
14.如权利要求13所述的方法，其特征在于，所述去噪声处理至少包括高斯滤波、均值滤波、傅里叶变换滤波、小波变换滤波和小波包变换滤波中的任一种方法，用于消除由于信号来自体表生理采集和分帧式的信号特征分析对睡眠状态刻画的不连续性，而在所述睡眠状态特征曲线中引入的高斯噪声或白噪声干扰。
15.如权利要求13所述的方法，其特征在于，所述模态分解至少包括经验模态分解、变分模态分解中的任一项。
16.如权利要求15所述的方法，其特征在于，所述经验模态分解的方法至少包括EMD、EEMD、CEEMD、 CEEMDAN、ICEEMDAN、ESMD中的任意一项。
17.如权利要求15所述的方法，其特征在于，所述变分模态分解的方法至少包括VMD、改进型VMD中的任意一项。
18.如权利要求13所述的方法，其特征在于，所述特征曲线高频阈值由所述睡眠状态特征曲线的睡眠生理状态特征线性组合融合的数据来源决定。
19.如权利要求1或2所述的方法，其特征在于，所述睡眠可持续性节点特征系数的一种计算公式，具体为：
；
其中，为所述睡眠可持续性节点特征系数，/>分别为所述目标睡眠状态平方特征值和所述睡眠可持续性IMF平方特征值。
20.如权利要求19所述的方法，其特征在于：所述睡眠可持续性指数的一种计算公式，具体为：
；
其中，为所述睡眠可持续性指数，/>分别为所述睡眠可持续性节点特征系数、所述模态分解方法修正系数、所述用户个性修正系数和所述睡眠时相分期修正系数，且/>，/>2。
21.如权利要求19所述的方法，其特征在于，所述睡眠可持续性指数曲线由按照时序拼接所述睡眠可持续性指数而生成得到。
22.如权利要求1或2所述的方法，其特征在于：所述对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析，生成睡眠可持续性动态调节策略，对用户睡眠过程进行动态调节及效果评估的具体步骤还包括：
对所述睡眠状态特征曲线进行动态地趋势预测分析，得到睡眠状态特征预测值；
对所述睡眠可持续性指数曲线进行动态地趋势预测分析，得到睡眠可持续性指数预测值；
根据所述睡眠状态特征预测值和所述睡眠可持续性指数预测值，按照睡眠可持续性检测调节周期和预设睡眠调节知识库，生成所述睡眠可持续性动态调节策略；
根据所述睡眠可持续性动态调节策略，连接并控制睡眠调节外围设备，对用户睡眠过程进行动态调节；
对调节效果进行动态跟踪评估，计算动态调节效果系数并生成动态调节效果曲线。
23.如权利要求22所述的方法，其特征在于：所述趋势预测分析的方法至少包括AR、MA、ARMA、ARIMA、SARIMA、VAR、深度学习中的任一项。
24.如权利要求22所述的方法，其特征在于：所述睡眠可持续性动态调节策略至少包括睡眠场景、睡眠时相、调节方式、执行方式、调节强度、调节时点、持续时间、目标调节值和装置控制参数；其中，所述调节方式至少包括声乐刺激、超声刺激、光刺激、电刺激、磁刺激、温度刺激、湿度刺激、触觉刺激和浓度调控中的任一方式，所述执行方式至少包括离体式和接触式中的任一方式。
25.如权利要求22所述的方法，其特征在于：所述睡眠调节外围设备至少包括声乐刺激设备、超声刺激设备、光刺激设备、电刺激设备、磁刺激设备、温度刺激设备、湿度刺激设备、触觉刺激设备和浓度调控设备中的任一项，并由具体的调节方式决定。
26.如权利要求23-25任一项所述的方法，其特征在于，所述动态调节效果系数的一种计算方式，具体为：
；
其中，为所述动态调节效果系数，/>为趋势预测分析方法修正系数且，/>为动态调节后的所述睡眠时相分期修正系数，/>分别为动态调节前、动态调节后的所述睡眠可持续性指数。
27.如权利要求26所述的方法，其特征在于，所述动态调节效果曲线由按照时序拼接所述动态调节效果系数而生成得到。
28.如权利要求22所述的方法，其特征在于：还包括建立或更新用户个性睡眠可持续性数据库的具体步骤：
根据所述睡眠时相曲线和所述睡眠可持续性指数曲线，计算不同睡眠时相下所述睡眠可持续性指数的分布特征，得到时相指数分布特征；
根据所述睡眠时相曲线和所述动态调节效果曲线，计算不同睡眠时相下所述动态调节效果系数的分布特征，得到时相效果分布特征；
计算所述睡眠可持续性指数曲线和所述动态调节效果曲线的相关性，得到动态调节综合效果系数；
动态收集检测调节过程数据和分析结果，初始化建立或持续动态更新所述用户个性睡眠可持续性数据库。
29.如权利要求28所述的方法，其特征在于：所述分布特征至少包括平均值、均方根、最大值、最小值、方差、标准差、变异系数、峰度和偏度中的任意一项。
30.如权利要求28或29所述的方法，其特征在于：所述相关性的计算方法至少包括相干性分析、皮尔逊相关分析、杰卡德相似分析、线性互信息分析、线性相关分析、欧氏距离分析、曼哈顿距离分析、切比雪夫距离分析中的任意一项。
31.如权利要求28所述的方法，其特征在于，所述方法还包括：根据所述用户个性睡眠可持续性数据库，动态迭代优化睡眠可持续性的动态检测调节过程方法策略，生成用户睡眠可持续性检测调节报告，并更新所述睡眠可持续性通用数据库。
32.如权利要求31所述的方法，其特征在于：所述根据所述用户个性睡眠可持续性数据库，动态迭代优化睡眠可持续性的动态检测调节过程方法策略，生成用户睡眠可持续性检测调节报告，并更新所述睡眠可持续性通用数据库的具体步骤还包括：
按照预设检测调节优化周期，根据所述用户个性睡眠可持续性数据库，动态迭代优化所述睡眠可持续性的所述动态检测调节过程方法策略，持续提高检测调节的效率效果；
按照预设报告生成周期，生成所述用户睡眠可持续性检测调节报告；
按照预设通用数据库更新周期，根据所述用户个性睡眠可持续性数据库和所述用户睡眠可持续性检测调节报告，将更新所述睡眠可持续性通用数据库的用户关键数据信息。
33.如权利要求32所述的方法，其特征在于：所述动态检测调节过程方法策略至少包括去高斯噪声处理方法参数、模态分解方法参数、所述特征曲线高频阈值、趋势预测分析方法参数、所述睡眠可持续性动态调节策略和动态调节效果评估方法参数。
34.如权利要求32所述的方法，其特征在于：所述用户睡眠可持续性检测调节报告至少包括所述用户关键生理信息、所述睡眠状态特征曲线、所述睡眠时相曲线、所述睡眠可持续性指数曲线、所述动态调节控效果曲线、所述时相指数分布特征、所述时相效果分布特征、所述动态调节综合效果系数、睡眠可持续性检测调节总结。
35.一种基于模态分解的睡眠可持续性检测调节系统，其特征在于，包括以下模块：
方案初始化模块，用于获取用户关键生理信息，根据用户关键生理信息筛选并初始化所述睡眠可持续性动态检测调节基本方案；
状态采集分析模块，用于根据所述睡眠可持续性动态检测调节基本方案，对用户睡眠生理信号进行采集和处理，得到睡眠状态特征曲线和睡眠时相曲线；
指数模态分析模块，用于对所述睡眠状态特征曲线进行去噪声处理和模态分解，根据特征曲线高频阈值提取高频模态分解信号，得到睡眠可持续性特征信号，计算睡眠可持续性指数和睡眠可持续性指数曲线；
睡眠动态调节模块，用于对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析，生成睡眠可持续性动态调节策略，对用户睡眠过程进行动态调节及效果评估；
数据运行管理模块，用于对所述系统的所有数据进行可视化管理、统一存储和运行管理；
所述睡眠可持续性指数的计算方法，具体为：
1)获取当前时帧对应的睡眠时相分期、目标睡眠状态特征曲线和所述睡眠可持续性特征信号；
2)分别对所述目标睡眠状态特征曲线和所述睡眠可持续性特征信号做平方运算并求和，得到目标睡眠状态平方特征值和睡眠可持续性IMF平方特征值；
3)根据所述目标睡眠状态平方特征值和所述睡眠可持续性IMF平方特征值的相对变化量，得到睡眠可持续性节点特征系数；
4)根据所述睡眠可持续性节点特征系数、模态分解方法修正系数、用户个性修正系数和睡眠时相分期修正系数，计算得到所述睡眠可持续性指数；其中，所述可持续性节点特征系数、所述睡眠时相分期修正系数均与所述睡眠可持续性指数正相关，所述模态分解方法修正系数和所述用户个性修正系数用于修正所述睡眠可持续性指数。
36.如权利要求35所述的系统，其特征在于，所述方案初始化模块还包括以下功能单元：
基础方案选择单元，用于获取所述用户关键生理信息，并从睡眠可持续性通用数据库或用户个性睡眠可持续性数据库中筛选与用户个性化生理匹配的所述睡眠可持续性动态检测调节基本方案；
检测调节初始化单元，用于根据所述睡眠可持续性动态检测调节基本方案，初始化用户睡眠可持续性的检测量化和动态调节的过程方法或过程策略。
37.如权利要求35所述的系统，其特征在于，所述状态采集分析模块还包括以下功能单元：
状态采集处理单元，用于对用户睡眠生理信号进行动态地监测采集和信号处理，得到睡眠生理状态数据；
状态特征提取单元，用于对所述睡眠生理状态数据进行动态地特征分析和特征选择，生成所述睡眠状态特征曲线；
睡眠时相分析单元，用于对所述睡眠生理状态数据进行动态地睡眠时相分析，生成所述睡眠时相曲线。
38.如权利要求35所述的系统，其特征在于，所述指数模态分析模块还包括以下功能单元：
高斯噪声处理单元，用于对所述睡眠状态特征曲线进行动态地去噪声处理，得到目标睡眠状态特征曲线；
信号模态分解单元，用于对所述目标睡眠状态特征曲线进行动态地模态分解，得到睡眠状态特征IMF分解信号集；
分解信号提取单元，用于对所述睡眠状态特征IMF分解信号集进行频谱分析，根据特征曲线高频阈值提取高频模态分解信号，生成睡眠可持续性特征信号；
睡眠指数计算单元，用于根据所述睡眠时相曲线、所述目标睡眠状态特征曲线和所述睡眠可持续性特征信号，动态计算得到所述睡眠可持续性指数并生成所述睡眠可持续性指数曲线。
39.如权利要求35-38任一项所述的系统，其特征在于，所述睡眠动态调节模块还包括以下功能单元：
状态趋势预测单元，用于对所述睡眠状态特征曲线进行动态地趋势预测分析，得到睡眠状态特征预测值；
指数趋势预测单元，用于对所述睡眠可持续性指数曲线进行动态地趋势预测分析，得到睡眠可持续性指数预测值；
动态策略生成单元，用于根据所述睡眠状态特征预测值和所述睡眠可持续性指数预测值，按照睡眠可持续性检测调节周期和预设睡眠调节知识库，生成所述睡眠可持续性动态调节策略；
睡眠动态调节单元，用于根据所述睡眠可持续性动态调节策略，连接并控制睡眠调节外围设备，对用户睡眠过程进行动态调节；
效果动态评估单元，用于对调节效果进行动态跟踪评估，计算动态调节效果系数并生成动态调节效果曲线。
40.如权利要求39所述的系统，其特征在于，还包括：数据统计更新模块，用于建立或更新用户个性睡眠可持续性数据库；
检测调节优化模块，用于根据所述用户个性睡眠可持续性数据库，动态迭代优化睡眠可持续性的动态检测调节过程方法策略，生成用户睡眠可持续性检测调节报告，并更新所述睡眠可持续性通用数据库。
41.如权利要求40所述的系统，其特征在于，所述数据统计更新模块还包括以下功能单元：
时相指数分析单元，用于根据所述睡眠时相曲线和所述睡眠可持续性指数曲线，计算不同睡眠时相下所述睡眠可持续性指数的分布特征，得到时相指数分布特征；
时相效果分析单元，用于根据所述睡眠时相曲线和所述动态调节效果曲线，计算不同睡眠时相下所述动态调节效果系数的分布特征，得到时相效果分布特征；
综合效果评价单元，用于计算所述睡眠可持续性指数曲线和所述动态调节效果曲线的相关性，得到动态调节综合效果系数；
个性数据库更新单元，用于动态收集检测调节过程数据和分析结果，初始化建立或持续动态更新所述用户个性睡眠可持续性数据库。
42.如权利要求40所述的系统，其特征在于，所述检测调节优化模块还包括以下功能单元：
过程动态优化单元，用于按照预设检测调节优化周期，根据所述用户个性睡眠可持续性数据库，动态迭代优化所述睡眠可持续性的所述动态检测调节过程方法策略，持续提高检测调节的效率效果；
用户报告生成单元，用于按照预设报告生成周期，生成所述用户睡眠可持续性检测调节报告；
用户报告管理单元，用于对所述用户睡眠可持续性检测调节报告的格式输出、展现形式进行统一管理；
通用数据库更新单元，用于按照预设通用数据库更新周期，根据所述用户个性睡眠可持续性数据库和所述用户睡眠可持续性检测调节报告，将更新所述睡眠可持续性通用数据库的用户关键数据信息。
43.如权利要求35-38任一项所述的系统，其特征在于，所述数据运行管理模块还包括以下功能单元：
用户信息管理单元，用于用户基本信息的登记输入、编辑、查询、输出和删除；
数据可视化管理单元，用于对所述系统中所有数据的可视化展现管理；
数据存储管理单元，用于对所述系统中所有数据的统一存储管理；
数据运营管理单元，用于对所述系统中所有数据的备份、迁移和导出。
44.一种基于模态分解的睡眠可持续性检测调节装置，其特征在于，包括以下模组：
方案初始化模组，用于获取用户关键生理信息，根据用户关键生理信息筛选并初始化所述睡眠可持续性动态检测调节基本方案；
状态采集分析模组，用于根据所述睡眠可持续性动态检测调节基本方案，对用户睡眠生理信号进行采集和处理，得到睡眠状态特征曲线和睡眠时相曲线；
指数模态分析模组，用于对所述睡眠状态特征曲线进行去噪声处理和模态分解，根据特征曲线高频阈值提取高频模态分解信号，得到睡眠可持续性特征信号，计算睡眠可持续性指数和睡眠可持续性指数曲线；
睡眠动态调节模组，用于对所述睡眠状态特征曲线和所述睡眠可持续性指数曲线进行趋势预测分析，生成睡眠可持续性动态调节策略，对用户睡眠过程进行动态调节及效果评估；
数据可视化模组，用于对所述装置中所有过程数据和结果数据的统一可视化展示管理；
数据管理中心模组，用于对所述装置中所有过程数据和结果数据的统一存储和数据运营管理；
所述睡眠可持续性指数的计算方法，具体为：
1)获取当前时帧对应的睡眠时相分期、目标睡眠状态特征曲线和所述睡眠可持续性特征信号；
2)分别对所述目标睡眠状态特征曲线和所述睡眠可持续性特征信号做平方运算并求和，得到目标睡眠状态平方特征值和睡眠可持续性IMF平方特征值；
3)根据所述目标睡眠状态平方特征值和所述睡眠可持续性IMF平方特征值的相对变化量，得到睡眠可持续性节点特征系数；
4)根据所述睡眠可持续性节点特征系数、模态分解方法修正系数、用户个性修正系数和睡眠时相分期修正系数，计算得到所述睡眠可持续性指数；其中，所述可持续性节点特征系数、所述睡眠时相分期修正系数均与所述睡眠可持续性指数正相关，所述模态分解方法修正系数和所述用户个性修正系数用于修正所述睡眠可持续性指数。
45.如权利要求44所述的装置，其特征在于，还包括：数据统计更新模组，用于建立或更新用户个性睡眠可持续性数据库；
检测调节优化模组，用于根据所述用户个性睡眠可持续性数据库，动态迭代优化睡眠可持续性的动态检测调节过程方法策略，生成用户睡眠可持续性检测调节报告，并更新所述睡眠可持续性通用数据库。
CN202310984123.2A 2023-08-07 2023-08-07 基于模态分解的睡眠可持续性检测调节方法、系统和装置 Active CN116779110B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310984123.2A CN116779110B (zh) 2023-08-07 2023-08-07 基于模态分解的睡眠可持续性检测调节方法、系统和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310984123.2A CN116779110B (zh) 2023-08-07 2023-08-07 基于模态分解的睡眠可持续性检测调节方法、系统和装置

Publications (2)

Publication Number Publication Date
CN116779110A true CN116779110A (zh) 2023-09-19
CN116779110B CN116779110B (zh) 2024-05-31

Family

ID=87986081

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310984123.2A Active CN116779110B (zh) 2023-08-07 2023-08-07 基于模态分解的睡眠可持续性检测调节方法、系统和装置

Country Status (1)

Country Link
CN (1) CN116779110B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118592924A (zh) * 2024-08-08 2024-09-06 北京中成康富科技股份有限公司 基于毫米波雷达的呼吸睡眠监测仪

Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1317943A1 (fr) * 2001-12-07 2003-06-11 Ela Medical Dispositif médical actif comprenant des moyens perfectionnés de discrimination des phases d'éveil et de sommeil
CN106473705A (zh) * 2016-09-21 2017-03-08 广州视源电子科技股份有限公司 用于睡眠状态监测的脑电信号处理方法和系统
WO2018180219A1 (ja) * 2017-03-27 2018-10-04 公立大学法人名古屋市立大学 睡眠状態の評価
WO2019053071A1 (en) * 2017-09-12 2019-03-21 Moona SLEEP METHOD AND APPARATUS FOR SLEEP REGULATION
KR20190035361A (ko) * 2017-09-26 2019-04-03 (주)코랄디자인 스마트 폰 앱 기반의 수면 모니터링 시스템 및 그 방법
CN111493822A (zh) * 2020-03-23 2020-08-07 济南国科医工科技发展有限公司 一种基于睡眠脑电的快速眼动期睡眠行为障碍分类方法
AU2020102907A4 (en) * 2020-10-21 2020-12-17 Ahamad, Danish MR Novel automated machine learning algorithms based system for sleep staging features analysis
CN115171850A (zh) * 2022-09-07 2022-10-11 深圳市心流科技有限公司 一种睡眠方案生成方法、装置、终端设备及存储介质
US20230043406A1 (en) * 2019-12-31 2023-02-09 South China Normal University Noninvasive method and system for sleep apnea detection
CN116072265A (zh) * 2023-04-06 2023-05-05 北京理工大学 基于时间自注意力与动态图卷积的睡眠分期分析系统及方法
CN116369866A (zh) * 2023-06-05 2023-07-04 安徽星辰智跃科技有限责任公司 基于小波变换的睡眠趋稳性量化及调节方法、系统和装置
CN116392087A (zh) * 2023-06-06 2023-07-07 安徽星辰智跃科技有限责任公司 基于模态分解的睡眠趋稳性量化及调节方法、系统和装置
CN116392083A (zh) * 2023-06-01 2023-07-07 安徽星辰智跃科技有限责任公司 基于局部变换的睡眠趋稳性量化及调节方法、系统和装置
CN116504357A (zh) * 2023-06-28 2023-07-28 安徽星辰智跃科技有限责任公司 基于小波分析的睡眠周期性检测及调节方法、系统和装置
CN116525063A (zh) * 2023-06-28 2023-08-01 安徽星辰智跃科技有限责任公司 基于时频分析的睡眠周期性检测及调节方法、系统和装置
CN116509337A (zh) * 2023-06-27 2023-08-01 安徽星辰智跃科技有限责任公司 基于局部分解的睡眠周期性检测及调节方法、系统和装置
CN116509338A (zh) * 2023-06-29 2023-08-01 安徽星辰智跃科技有限责任公司 基于模态分析的睡眠周期性检测及调节方法、系统和装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2002951984A0 (en) * 2002-10-10 2002-10-31 Compumedics Limited Sleep quality and auto cpap awakening
JP7301681B2 (ja) * 2019-09-05 2023-07-03 住友理工株式会社 睡眠状態推定方法および睡眠状態推定装置
CN114067435A (zh) * 2021-11-15 2022-02-18 山东大学 一种基于伪3d卷积网络与注意力机制的睡眠行为检测方法和系统

Patent Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1317943A1 (fr) * 2001-12-07 2003-06-11 Ela Medical Dispositif médical actif comprenant des moyens perfectionnés de discrimination des phases d'éveil et de sommeil
CN106473705A (zh) * 2016-09-21 2017-03-08 广州视源电子科技股份有限公司 用于睡眠状态监测的脑电信号处理方法和系统
WO2018180219A1 (ja) * 2017-03-27 2018-10-04 公立大学法人名古屋市立大学 睡眠状態の評価
WO2019053071A1 (en) * 2017-09-12 2019-03-21 Moona SLEEP METHOD AND APPARATUS FOR SLEEP REGULATION
KR20190035361A (ko) * 2017-09-26 2019-04-03 (주)코랄디자인 스마트 폰 앱 기반의 수면 모니터링 시스템 및 그 방법
US20230043406A1 (en) * 2019-12-31 2023-02-09 South China Normal University Noninvasive method and system for sleep apnea detection
CN111493822A (zh) * 2020-03-23 2020-08-07 济南国科医工科技发展有限公司 一种基于睡眠脑电的快速眼动期睡眠行为障碍分类方法
AU2020102907A4 (en) * 2020-10-21 2020-12-17 Ahamad, Danish MR Novel automated machine learning algorithms based system for sleep staging features analysis
CN115171850A (zh) * 2022-09-07 2022-10-11 深圳市心流科技有限公司 一种睡眠方案生成方法、装置、终端设备及存储介质
CN116072265A (zh) * 2023-04-06 2023-05-05 北京理工大学 基于时间自注意力与动态图卷积的睡眠分期分析系统及方法
CN116392083A (zh) * 2023-06-01 2023-07-07 安徽星辰智跃科技有限责任公司 基于局部变换的睡眠趋稳性量化及调节方法、系统和装置
CN116369866A (zh) * 2023-06-05 2023-07-04 安徽星辰智跃科技有限责任公司 基于小波变换的睡眠趋稳性量化及调节方法、系统和装置
CN116392087A (zh) * 2023-06-06 2023-07-07 安徽星辰智跃科技有限责任公司 基于模态分解的睡眠趋稳性量化及调节方法、系统和装置
CN116509337A (zh) * 2023-06-27 2023-08-01 安徽星辰智跃科技有限责任公司 基于局部分解的睡眠周期性检测及调节方法、系统和装置
CN116504357A (zh) * 2023-06-28 2023-07-28 安徽星辰智跃科技有限责任公司 基于小波分析的睡眠周期性检测及调节方法、系统和装置
CN116525063A (zh) * 2023-06-28 2023-08-01 安徽星辰智跃科技有限责任公司 基于时频分析的睡眠周期性检测及调节方法、系统和装置
CN116509338A (zh) * 2023-06-29 2023-08-01 安徽星辰智跃科技有限责任公司 基于模态分析的睡眠周期性检测及调节方法、系统和装置

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
何敏;: "基于Hilbert-Huang变换获得睡眠脑电信号特征", 电脑知识与技术, no. 20, 15 July 2015 (2015-07-15) *
卢福兰;徐维超;: "基于经验模式分解的自动睡眠分期", 工业控制计算机, no. 06 *
张俊民;王蓓;王行愚;: "基于特征融合的ARMA短时睡眠状态分析", 华东理工大学学报(自然科学版), no. 03 *
肖余粮;和卫星;陈晓平;吉奕;吴秋明;: "小波变换和小波熵在睡眠脑电信号变化特性研究中的应用价值", 中国临床康复, no. 25 *
陈晓平;和卫星;姚静;: "基于小波变换与近似熵的睡眠分期研究", 中国组织工程研究与临床康复, no. 09 *
陈科;单姗;郑红梅;: "具有自适应性的实时睡眠信号处理算法研究", 电子测量与仪器学报, no. 10, 15 October 2016 (2016-10-15) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118592924A (zh) * 2024-08-08 2024-09-06 北京中成康富科技股份有限公司 基于毫米波雷达的呼吸睡眠监测仪

Also Published As

Publication number Publication date
CN116779110B (zh) 2024-05-31

Similar Documents

Publication Publication Date Title
CN116369866B (zh) 基于小波变换的睡眠趋稳性量化及调节方法、系统和装置
CN116392083B (zh) 基于局部变换的睡眠趋稳性量化及调节方法、系统和装置
CN116392085B (zh) 基于趋势分析的睡眠趋稳性量化及调节方法、系统和装置
CN116392087B (zh) 基于模态分解的睡眠趋稳性量化及调节方法、系统和装置
CN116509337B (zh) 基于局部分解的睡眠周期性检测及调节方法、系统和装置
CN116779110B (zh) 基于模态分解的睡眠可持续性检测调节方法、系统和装置
CN115910351B (zh) 一种睡眠周期性检测量化及辅助干预的方法、系统和装置
CN115804573A (zh) 一种睡眠深度量化及干预的方法、系统和装置
CN116509338B (zh) 基于模态分析的睡眠周期性检测及调节方法、系统和装置
CN116058805B (zh) 一种睡眠记忆活跃水平动态调节的方法、系统和装置
CN116705247B (zh) 基于局部分解的睡眠可持续性检测调节方法、系统和装置
CN116058804A (zh) 一种睡眠情绪活跃水平动态调节的方法、系统和装置
CN115969330B (zh) 一种睡眠情绪活跃水平检测量化的方法、系统和装置
CN116035536B (zh) 一种睡眠行为活跃水平检测量化的方法、系统和装置
CN116392088B (zh) 基于时频分析的睡眠趋稳性量化及调节方法、系统和装置
CN116509336B (zh) 基于波形分析的睡眠周期性检测及调节方法、系统和装置
CN116525063B (zh) 基于时频分析的睡眠周期性检测及调节方法、系统和装置
CN116504357B (zh) 基于小波分析的睡眠周期性检测及调节方法、系统和装置
CN116671875B (zh) 基于小波变换的睡眠可持续性检测调节方法、系统和装置
CN116687356B (zh) 基于时频分析的睡眠可持续性检测调节方法、系统和装置
CN116682535B (zh) 基于数值拟合的睡眠可持续性检测调节方法、系统和装置
CN116312928B (zh) 一种眼疲劳恢复进程的辅助调节方法、系统和装置
CN117747118B (zh) 睡眠呼吸周期评价及辅助调节方法、系统和装置
CN116312972B (zh) 基于眼刺激的睡眠记忆情感张力调节方法、系统和装置
CN117747117B (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