CN111631688B - 睡眠自动分期的算法 - Google Patents
睡眠自动分期的算法 Download PDFInfo
- Publication number
- CN111631688B CN111631688B CN202010591697.XA CN202010591697A CN111631688B CN 111631688 B CN111631688 B CN 111631688B CN 202010591697 A CN202010591697 A CN 202010591697A CN 111631688 B CN111631688 B CN 111631688B
- Authority
- CN
- China
- Prior art keywords
- data
- sleep
- sequence
- segment
- data set
- 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.)
- Expired - Fee Related
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 claims abstract description 25
- 238000012937 correction Methods 0.000 claims abstract description 18
- 239000000284 extract Substances 0.000 claims abstract description 3
- 239000013598 vector Substances 0.000 claims description 32
- 238000012549 training Methods 0.000 claims description 31
- 230000008667 sleep stage Effects 0.000 claims description 15
- 238000013528 artificial neural network Methods 0.000 claims description 12
- 238000012360 testing method Methods 0.000 claims description 11
- 238000012795 verification Methods 0.000 claims description 11
- 238000005312 nonlinear dynamic Methods 0.000 claims description 9
- 238000012546 transfer Methods 0.000 claims description 7
- 238000011156 evaluation Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 230000011218 segmentation Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 4
- 238000012544 monitoring process Methods 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000010924 continuous production Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 230000002457 bidirectional effect Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 8
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000010801 machine learning Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000004622 sleep time Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 210000004373 mandible Anatomy 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
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/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/48—Other medical applications
- A61B5/4806—Sleep evaluation
- A61B5/4815—Sleep quality
-
- 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/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- 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/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- 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/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- 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/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/044—Recurrent networks, e.g. Hopfield networks
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Medical Informatics (AREA)
- Theoretical Computer Science (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Signal Processing (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computational Linguistics (AREA)
- Fuzzy Systems (AREA)
- Anesthesiology (AREA)
- Computing Systems (AREA)
- Software Systems (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
本发明公开了一种睡眠自动分期的算法,包括如下步骤:a.特征层利用多层感知器提取抽象特征合并传统基于专家经验的手工特征作为睡眠的信息表征;b.然后在机器模型层利用双向门控循环单元Bi‑GRU作为网络模型;c.最后在修正层利用条件随机场CRF方法作为时间连续性修正。解决现有睡眠自动分期算法稳定性、通用性、普适性以及实用性差的问题。
Description
技术领域
本发明涉及睡眠算法领域,具体涉及一种睡眠自动分期的算法。
背景技术
在人的睡眠过程中利用生理信息根据国际睡眠医学标准可将不同睡眠时间段划分为WAKE、REM、N1、N2和N3期,共五个正常类别。
近年来利用计算机设备和程序监测人睡眠状态并自动地将其分类的方法开始广泛应用。虽然现代电子信息技术、机器学习理论、生物医学工程等方面飞速发展,但利用机器学习理论的睡眠自动分期方法在科研、医学应用、消费电子等领域中仍没有国际公认的标准。主要原因包括人睡眠的基础机制理论并不透彻完善、临床医生和研究人员对睡眠人工分期国际标准缺乏足够理解和信任导致睡眠状态分类一致性低、睡眠监测***的开发人员缺少专家经验、非正常睡眠模式混杂在正常睡眠模式中等。
凡是睡眠自动分期算法涉及特征工程的,开发人员很难有专家经验深刻理解睡眠的生理过程,专家知识主要来自美国睡眠医学学会的睡眠技术指导手册以及实践中摸索得到的人工分期的准则。
凡是睡眠自动分期算法涉及神经网络的,其算法模型很难拟合好实际数据,造成欠拟合模型情况,机器的判定结果的准确率远低于预期。
总体来说,现有睡眠自动分期的***不足以精确达到区分正常睡眠状态。现有的基于特征工程的机器学习算法和基于神经网络的机器学习算法在少量数据范围内达到60%到90%的准确率,现有的算法都将整夜睡眠的生理信号数据片段作独立性假设没有时间连续性的考量,整个自动分期技术的稳定性、通用性、普适性和实用性等亟待优化。
发明内容
本发明的目的在于提供一种睡眠自动分期的算法,解决现有睡眠自动分期算法稳定性、通用性、普适性以及实用性差的问题。
为解决上述的技术问题,本发明采用以下技术方案:
一种睡眠自动分期的算法,包括如下步骤:a.特征层利用多层感知器提取抽象特征合并传统基于专家经验的手工特征作为睡眠的信息表征;b.然后在机器模型层利用双向门控循环单元Bi-GRU作为网络模型;c.最后在修正层利用条件随机场CRF方法作为时间连续性修正。
作为本发明的进一步优选,所述步骤a中利用多层感知器提取抽象特征合并传统基于专家经验的手工特征作为睡眠的信息表征包括如下步骤:S1、将由人睡眠中采集的生理数据制作数据集,包含脑电信号、眼电信号和下颌肌电信号,C={C1,C2,...,CN},C表示所有人的数据集合,N表示人数(或完成整夜睡眠的次数);S2、将每个人的一整夜连续的数据按时间顺序分成S段,各个段代表着不同睡眠阶段信息,随着时间推移,各段睡眠数据需要用不同算法模型来学习拟合,因此每段单独制作一个数据集,产生S个段数据集,每个段数据集包含N×M个样本数据,M为每段数据包含样本点的数量,一个样本点表示一个睡眠期,其时间跨度通常为30秒,每个样本点包含L个信号采样点,所以C的数据集样本量为S×N×M×L;S3、上述S个段数据集分别进行特征工程,其中方式包括提取抽象特征:自编码器;提取传统特征:时域特征、频域特征和非线性动力学特征,所有特征拼接成向量,每段数据生成M个向量,并顺序排列形成序列,即段序列;所述步骤b中的在机器模型层利用双向门控循环单元Bi-GRU作为网络模型包括如下步骤:S4、在机器模型层,双向门控循环单元Bi-GRU网络模型利用上述段序列作为样本来训练,S个段数据集分别包括N×M个段序列,各个段数据集分别按一定比例划分训练集和测试集,各个段数据集分别训练结束后,将模型保存;S5、用于训练的样本点输入到已经训练完成的双向门控循环单元Bi-GRU网络模型后,在最后得到睡眠阶段的分类标签;S6、将每夜晚睡眠的S段样本的标签序列拼接成一个超长序列,即一整夜数据最后对应一个完整的整夜标签序列,整夜标签序列是一维向量,其维度为T=S×M,含有N条整夜标签序列集依然按照先前划分是三集分开,训练集包含Ct条,验证集包含Cv条,测试集包含C条;所述步骤c中的在修正层利用条件随机场CRF方法作为时间连续性修正包括步骤S7,步骤S7是将在S6步骤中得到的标签训练集序列输入到修正层,利用条件随机场CRF方法在上下文信息转移提取方面的优势,具体指CRF线性链方法为上述序列建模,通过维特比算法解码出最优标签序列路径,修正整夜标签序列,使其不断和专家人工睡眠分期判定结果的标签序列吻合。
作为本发明的进一步优选,所述步骤S1中还包括了如下子步骤:S11、在人睡眠中通过带有生理信号采集功能的多导睡眠仪设备按美国睡眠医学学会标准来监测、记录和存储人体生理信号数据;S12、不同种类的生理信号的原始信号数据经过采样后数字化后,分别进行零相位数字滤波以防止具有非平稳性质的生理信号出现相位失真,去除信号中极低频的基线、工频噪声和高频噪声,信号预处理即完成;S13、将上述生理信号中眼电信号、下颌肌电信号和脑电信号提取,利用三导信号进行原始数据集C={C1,C2,...,CN}制作,来自N个人个体的N条整夜睡眠数据,按人划分各自子集。
作为本本发明的进一步优选,所述步骤S2中还包括了如下子步骤:S21、合理设置切分阶段的模式,单人一整夜睡眠数据是按时间连续的,本发明算法将其整个连续过程划分S个阶段,被切分段的生理信号代表着睡眠过程不同信息,算法模拟专家人工睡眠分期时会对整夜信号数据有若干个趋势阶段的预判经验,更具一个阶段具有多长时间的跨度来设置每段样本点个数M;S22、为算法模型的前期准备制作具体数据集,在上述原始数据集基础上划分出S个段数据集,数据量大小不能出错,各段数据集均包含来自N个人的N×M个样本点,样本点是睡眠分期最小单位,其时间跨度为30秒,一个样本点对应于一个类别标签,类别标签有五个,其包括WAKE、REM、N1、N2和N3;S23、将标签和数据对齐,核对段数据集制作并保存至文件。以上数据集的制作是本发明的关键步骤,深刻影响整个算法的结果和性能。
作为本本发明的进一步优选,所述步骤S3中还包括了如下子步骤:S31、提取传统特征,根据信号计算特征向量xi维度为k。特征向量xi中的m个特征包括:时域特征量、频域特征量和非线性动力学特征量,时域特征量包含统计学特征量和几何学特征量。频域特征量包括功率谱密度特征量和时频特征量,非线性动力学特征量包括分形维数特征量和复杂度特征量,上述每种特征量均由各自参数和计算方式决定若干具体特征量;S32、提取抽象特征。通过人工神经网络领域的自编码器进行抽象特征提取,利用其无监督学习特性,拟合出能够对输入数据高效表示的人工神经网络,不添加额外人工辅助工作,而将输入的信号数据用固定的低维度的向量高效表征,即为自编码,其输出维度一般小于输入的信号维度,即自编码器数据降维特性。本发明采用自编码器,具有多个编码层,编码的复杂程度取决于神经网络堆叠层的层数,堆叠层数适当增加,可以有效地对输入数据进行压缩表示;S33、对上述的各个样本点进行特征工程,拼接成特征向量ξi,(ξi∈Rkξi∈Rk)各段中样本点必须按时间顺序排列成段序列,对应特征向量也排列成段序列seq,实际上段序列构造的整个特征空间为三维张量X∈RN×M×k。
作为本本发明的进一步优选,所述步骤S4中还包括了如下子步骤:S41、将seq生成的特征空间的张量X划分为最终训练数据集、验证数据集和测试数据集;S42、将上述按带有时间顺序的段序列seq输入到双向门控循环单元Bi-GRU网络模型,seq生成的特征空间的张量X∈RN×M×k的应满足Bi-GRU输入层要求,合理设置网络模型的结构、训练模式和初始参数,之后加载数据集开始模型的训练;S43、上述双向门控循环单元Bi-GRU网络模型在达到预设定的终止条件后保存模型。
作为本本发明的进一步优选,所述步骤S5中还包括了如下子步骤:S51、利用上述网络模型将数据全部前向传播,得到数据样本与之对应的睡眠分类阶段的标签,其中数据样本为上述段数据集生成的训练集和验证集;S52、将上述机器模型的分类结果标签与之专家人工分类结果标签比对,并且记录各个数据集的准确率、召回率和F1分数等评价指标,网络模型至此构建完成。
作为本本发明的进一步优选,所述步骤S7中还包括了如下子步骤:S71、构建条件随机场CRF模型,将步骤S6中的标签序列训练集和对应专家人工分期的标签序列输入到条件随机场CRF模型,设置特征函数个数K,迭代训练出最优的参数,进而得到条件随机场的条件概率P(y|x),即模型学习到睡眠阶段时间依赖性转移的上下文信息的,阶段转移和睡眠时间息息相关,也是本发明针对睡眠的特性利用CRF模型的概率转移的关键一步;S72、测试CRF模型修正后的结果,利用上述条件概率P(y|x)和验证集里的标签序列xs来计算出最优的标签序列y*,最后计算准确率等评价指标。统计对比CRF修正模型和上述网络模型的结果。
与现有技术相比,本发明至少能达到以下有益效果中的一项:
1.选择合适的数据段切分,使得算法对于睡眠各阶段的适应性更强;
2.传统特征和抽象特征融合,扩增特征的表达能力,使得算法的精准度更高;
3.通过时间关联的网络和概率修正,有效增强算法在时间连续性上的解释。
附图说明
图1为本发明算法整体框图。
图2为本发明数据集制作流程示意图。
图3为本发明数据集制作示意图。
图4为本发明数据样本结构示意图。
图5为本发明栈式自编码器框图。
图6为本发明栈式自编码器第一层示意图。
图7为本发明栈式自编码器第二层示意图。
图8为本发明栈式自编码器输出层示意图。
图9为本发明栈式自编码器整体示意图。
图10为本发明Bi-GRU网络模型示意图。
图11为本发明GRU节点内部结构。
图12为本发明加入CRF修正的整体机器模型层示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
具体实施例1:
图1、图2、图3、图4、图5、图6、图7、图8、图9、图10、图11、图12示出了一种睡眠自动分期的算法,如图1所示,本发明整体思路按照数据集制作,特征工程和模型分层三大流程进行逐步细化一个睡眠自动分期算法模型的构建。
如图2所示,数据集制作流程分为数据采集、数据预处理和整夜数据分段三个子部分。首先,将由人睡眠中采集的生理数据制作数据集,包含脑电信号(Eletroencephalogram,EEG)、眼电信号(Electrooculogram,EOG)和下颌肌电信号(Electromyography,EMG)等常规生理信号。C={C1,C2,...,CN},C表示所有人的数据集合,N表示人数(或完成整夜睡眠的次数)。
在人睡眠中通过带有生理信号采集功能的多导睡眠仪设备按美国睡眠医学学会标准来监测、记录和存储人体生理信号数据。
不同种类的生理信号的原始信号数据经过采样后数字化后,分别进行零相位数字滤波以防止具有非平稳性质的生理信号出现相位失真。去除信号中极低频的基线、工频噪声和高频噪声,信号预处理即完成。
将上述生理信号中眼电信号、下颌肌电信号和脑电信号提取,利用三导信号进行原始数据集C={C1,C2,...,CN}制作。来自N个人个体的N条整夜睡眠数据,按人划分各自子集。
将数据的结构进行可视化,如图3和图4所示,把每个人的一整夜连续的数据按时间顺序分成S段,各个段代表着不同睡眠阶段信息,随着时间推移,各段睡眠数据需要用不同算法模型来学习拟合。因此每段单独制作一个数据集,产生S个段数据集,每个段数据集包含N×M个样本数据,M为每段数据包含样本点的数量。一个样本点表示一个睡眠期,其时间跨度通常为30秒。每个样本点包含L个信号采样点。所以C的数据集样本量为S×N×M×L。
合理设置切分阶段的模式。单人一整夜睡眠数据是按时间连续的,本发明算法将其整个连续过程划分S个阶段,被切分段的生理信号代表着睡眠过程不同信息。算法模拟专家人工睡眠分期时会对整夜信号数据有若干个趋势阶段的预判经验。更具一个阶段具有多长时间的跨度来设置每段样本点个数M。
为算法模型的前期准备制作具体数据集,在上述原始数据集基础上划分出S个段数据集。数据量大小不能出错,各段数据集均包含来自N个人的N×M个样本点,样本点是睡眠分期最小单位,其时间跨度为30秒,一个样本点对应于一个类别标签。类别标签有五个,其包括WAKE、REM、N1、N2和N3。
将标签和数据对齐,核对段数据集制作并保存至文件。以上数据集的制作是本发明的关键步骤,深刻影响整个算法的结果和性能。
上述S个段数据集分别进行特征工程,其中方式包括提取抽象特征:自编码器;提取传统特征:时域特征、频域特征和非线性动力学特征。所有特征拼接成向量,每段数据生成M个向量,并顺序排列形成序列,即段序列。
其中提取传统特征。根据信号计算特征向量xi维度为k。特征向量xi中的m个特征包括:时域特征量、频域特征量和非线性动力学特征量。时域特征量包含统计学特征量和几何学特征量。频域特征量包括功率谱密度特征量和时频特征量。非线性动力学特征量包括分形维数特征量和复杂度特征量。上述每种特征量均由各自参数和计算方式决定若干具体特征量。
其中提取抽象特征。通过人工神经网络领域的自编码器(Autoencoders)进行抽象特征提取。利用其无监督学习特性,拟合出能够对输入数据高效表示的人工神经网络。不添加额外人工辅助工作,而将输入的信号数据用固定的低维度的向量高效表征,即为自编码。其输出维度一般小于输入的信号维度,即自编码器数据降维特性。本发明采用自编码器(StackedAutoencoders,SA),具有多个编码层,编码的复杂程度取决于神经网络堆叠层的层数,堆叠层数适当增加,可以有效地对输入数据进行压缩表示。
如图5、6、7、8、9展示栈式自编码器的网络结构原理,提取抽象特征。
对上述的各个样本点进行特征工程,拼接成特征向量ξi,(ξi∈Rk)各段中样本点必须按时间顺序排列成段序列,对应特征向量也排列成段序列seq。实际上段序列构造的整个特征空间为三维张量X∈RN×M×k。
经过数据集制作流程和特征工程,数据集结构即达到本发明的设计思想,充分体现特征融合和数据的时间依赖性。
如图10所示,在机器模型层,双向门控循环单元Bi-GRU网络模型利用上述段序列作为样本来训练,S个段数据集分别包括N×M个段序列。
将seq生成的特征空间的张量X划分为最终训练数据集、验证数据集和测试数据集。
交叉验证后,各个段数据集分别训练结束,最后将模型保存。
将上述按带有时间顺序的段序列seq输入到双向门控循环单元Bi-GRU网络模型。seq生成的特征空间的张量X∈RN×M×k的应满足Bi-GRU输入层要求。合理设置网络模型的结构、训练模式和初始参数,之后加载数据集开始模型的训练。
上述双向门控循环单元Bi-GRU网络模型在达到预设定的终止条件后保存模型。
利用上述网络模型将数据全部前向传播,得到数据样本与之对应的睡眠分类阶段的标签。其中数据样本为上述段数据集生成的训练集和验证集。
将上述机器模型的分类结果标签与之专家人工分类结果标签比对,并且记录各个数据集的准确率、召回率和F1分数等评价指标,网络模型至此构建完成。
用于训练的样本点输入到已经训练完成的双向门控循环单元Bi-GRU网络模型后,在最后得到睡眠阶段的分类标签。
将每夜晚睡眠的S段样本的标签序列拼接成一个超长序列,即一整夜数据最后对应一个完整的整夜标签序列。整夜标签序列是一维向量,其维度为T=S×M。含有N条整夜标签序列集依然按照先前划分是三集分开,训练集包含Ct条,验证集包含Cv条,测试集包含C条。
将上述标签训练集序列输入到修正层,利用条件随机场CRF方法在上下文信息转移提取方面的优势,具体指CRF线性链方法为上述序列建模,通过维特比算法解码出最优标签序列路径。修正整夜标签序列,使其不断和专家人工睡眠分期判定结果的标签序列吻合。
构建条件随机场CRF模型。将步骤S6中的标签序列训练集和对应专家人工分期的标签序列输入到条件随机场CRF模型,设置特征函数个数K,迭代训练出最优的参数,进而得到条件随机场的条件概率P(y|x),即模型学习到睡眠阶段时间依赖性转移的上下文信息的,阶段转移和睡眠时间息息相关,也是本发明针对睡眠的特性利用CRF模型的概率转移的关键一步。
测试CRF模型修正后的结果。利用上述条件概率P(y|x)和验证集里的标签序列xs来计算出最优的标签序列y*。最后计算准确率等评价指标。统计对比CRF修正模型和上述网络模型的结果。
尽管这里参照本发明的多个解释性实施例对本发明进行了描述,但是,应该理解,本领域技术人员可以设计出很多其他的修改和实施方式,这些修改和实施方式将落在本申请公开的原则范围和精神之内。更具体地说,在本申请公开、附图和权利要求的范围内,可以对主题组合布局的组成部件和/或布局进行多种变型和改进。除了对组成部件和/或布局进行的变形和改进外,对于本领域技术人员来说,其他的用途也将是明显的。
Claims (7)
1.一种睡眠自动分期的算法,其特征在于:包括如下步骤:a.特征层利用多层感知器提取抽象特征合并传统基于专家经验的手工特征作为睡眠的信息表征;b.然后在机器模型层利用双向门控循环单元Bi-GRU作为网络模型;c.最后在修正层利用条件随机场CRF方法作为时间连续性修正;
所述步骤a中利用多层感知器提取抽象特征合并传统基于专家经验的手工特征作为睡眠的信息表征包括如下步骤:S1、将由人睡眠中采集的生理数据制作数据集,包含脑电信号、眼电信号和下颌肌电信号,C={C 1,C 2,...,C N},C表示所有人的数据集合,N表示人数;S2、将每个人的一整夜连续的数据按时间顺序分成S段,各个段代表着不同睡眠阶段信息,随着时间推移,各段睡眠数据需要用不同算法模型来学习拟合,因此每段单独制作一个数据集,产生S个段数据集,每个段数据集包含N×M个样本数据,M为每段数据包含样本点的数量,一个样本点表示一个睡眠期,其时间跨度通常为30秒,每个样本点包含L个信号采样点,所以C的数据集样本量为S×N×M×L;S3、上述S个段数据集分别进行特征工程,其中方式包括提取抽象特征:自编码器;提取传统特征:时域特征、频域特征和非线性动力学特征,所有特征拼接成向量,每段数据生成M个向量,并顺序排列形成序列,即段序列;所述步骤b中的在机器模型层利用双向门控循环单元Bi-GRU作为网络模型包括如下步骤:S4、在机器模型层,双向门控循环单元Bi-GRU网络模型利用上述段序列作为样本来训练,S个段数据集分别包括N×M个段序列,各个段数据集分别按一定比例划分训练集和测试集,各个段数据集分别训练结束后,将模型保存;S5、用于训练的样本点输入到已经训练完成的双向门控循环单元Bi-GRU网络模型后,在最后得到睡眠阶段的分类标签;S6、将每夜晚睡眠的S段样本的标签序列拼接成一个超长序列,即一整夜数据最后对应一个完整的整夜标签序列,整夜标签序列是一维向量,其维度为T=S×M,含有N条整夜标签序列集依然按照先前划分是三集分开,训练集包含Ct条,验证集包含Cv条,测试集包含Cc条;所述步骤c中的在修正层利用条件随机场CRF方法作为时间连续性修正包括步骤S7,步骤S7是将在S6步骤中得到的标签训练集序列输入到修正层,利用条件随机场CRF方法在上下文信息转移提取方面的优势,具体指CRF线性链方法为上述序列建模,通过维特比算法解码出最优标签序列路径,修正整夜标签序列,使其不断和专家人工睡眠分期判定结果的标签序列吻合。
2.根据权利要求1所述的睡眠自动分期的算法,其特征在于:所述步骤S1中还包括了如下子步骤:S11、在人睡眠中通过带有生理信号采集功能的多导睡眠仪设备按美国睡眠医学学会标准来监测、记录和存储人体生理数据;S12、不同种类的生理信号的原始信号数据经过采样后数字化后,分别进行零相位数字滤波以防止具有非平稳性质的生理信号出现相位失真,去除信号中极低频的基线、工频噪声和高频噪声,信号预处理即完成;S13、将上述生理信号中眼电信号、下颌肌电信号和脑电信号提取,利用三导信号进行原始数据集C={C1,C 2…C N}制作,来自N个人个体的N条整夜睡眠数据,按人划分各自子集。
3.根据权利要求1所述的睡眠自动分期的算法,其特征在于:所述步骤S2中还包括了如下子步骤:S21、合理设置切分阶段的模式,单人一整夜睡眠数据是按时间连续的,算法将其整个连续过程划分S个阶段,被切分段的生理信号代表着睡眠过程不同信息,算法模拟专家人工睡眠分期时会对整夜信号数据有若干个趋势阶段的预判经验,根据一个阶段具有多长时间的跨度来设置每段样本点个数M;S22、为算法模型的前期准备制作具体数据集,在原始数据集基础上划分出S个段数据集,数据量大小不能出错,各段数据集均包含来自N个人的N×M个样本点,样本点是睡眠分期最小单位,其时间跨度为30秒,一个样本点对应于一个类别标签,类别标签有五个,类别标签包括WAKE、REM、N1、N2和N3;S23、将标签和数据对齐,核对段数据集制作并保存至文件。
4.根据权利要求1所述的睡眠自动分期的算法,其特征在于:所述步骤S3中还包括了如下子步骤:S31、提取传统特征,根据信号计算特征向量x i维度为k,特征向量x i中的m个特征包括:时域特征量、频域特征量和非线性动力学特征量,时域特征量包含统计学特征量和几何学特征量,频域特征量包括功率谱密度特征量和时频特征量,非线性动力学特征量包括分形维数特征量和复杂度特征量,上述每种特征量均由各自参数和计算方式决定若干具体特征量;S32、提取抽象特征,通过人工神经网络领域的自编码器进行抽象特征提取,利用其无监督学习特性,拟合出能够对输入数据高效表示的人工神经网络,不添加额外人工辅助工作,而将输入的信号数据用固定的低维度的向量高效表征,即为自编码,其输出维度一般小于输入的信号维度,即自编码器数据降维特性,采用自编码器,具有多个编码层,编码的复杂程度取决于神经网络堆叠层的层数,堆叠层数适当增加,可以有效地对输入数据进行压缩表示;S33、对上述的各个样本点进行特征工程,拼接成特征向量ξi,(ξi∈R k)各段中样本点必须按时间顺序排列成段序列,对应特征向量也排列成段序列seq,实际上段序列构造的整个特征空间为三维张量X ∈ RN× M × k。
5.根据权利要求1所述的睡眠自动分期的算法,其特征在于:所述步骤S4中还包括了如下子步骤:S41、将seq生成的特征空间的三维张量X划分为最终训练数据集、验证数据集和测试数据集;S42、将按带有时间顺序的段序列seq输入到双向门控循环单元Bi-GRU网络模型,seq生成的特征空间的张量X ∈ RN× M × k应满足Bi-GRU输入层要求,合理设置网络模型的结构、训练模式和初始参数,之后加载数据集开始模型的训练;S43、上述双向门控循环单元Bi-GRU网络模型在达到预设定的终止条件后保存模型。
6.根据权利要求1所述的睡眠自动分期的算法,其特征在于:所述步骤S5中还包括了如下子步骤:S51、利用上述网络模型将数据全部前向传播,得到数据样本与之对应的睡眠分类阶段的标签,其中数据样本为上述段数据集生成的训练集和验证集;S52、将上述机器模型的分类结果标签与之专家人工分类结果标签比对,并且记录各个段数据集的评价指标,包括准确率、召回率和F1分数,网络模型至此构建完成。
7.根据权利要求1所述的睡眠自动分期的算法,其特征在于:所述步骤S7中还包括了如下子步骤:S71、构建条件随机场CRF模型,将步骤S6中的标签序列训练集和对应专家人工分期的标签序列输入到条件随机场CRF模型,设置特征函数个数K,迭代训练出最优的参数,进而得到条件随机场的条件概率P(y|x);S72、测试CRF模型修正后的结果,利用上述条件概率P(y|x)和验证集里的标签序列x s来计算出最优的标签序列y*,最后计算评价指标,包括准确率,统计对比CRF修正模型和上述网络模型的结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010591697.XA CN111631688B (zh) | 2020-06-24 | 2020-06-24 | 睡眠自动分期的算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010591697.XA CN111631688B (zh) | 2020-06-24 | 2020-06-24 | 睡眠自动分期的算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111631688A CN111631688A (zh) | 2020-09-08 |
CN111631688B true CN111631688B (zh) | 2021-10-29 |
Family
ID=72323139
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010591697.XA Expired - Fee Related CN111631688B (zh) | 2020-06-24 | 2020-06-24 | 睡眠自动分期的算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111631688B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112263218A (zh) * | 2020-10-12 | 2021-01-26 | 上海大学 | 睡眠分期方法及装置 |
CN112294342A (zh) * | 2020-10-30 | 2021-02-02 | 哈尔滨理工大学 | 一种基于深度残差Mask-CCNN的睡眠分期方法 |
CN113080864B (zh) * | 2021-04-07 | 2022-02-01 | 电子科技大学 | 一种通过自动睡眠分期结果的常见睡眠疾病检测方法 |
CN113576410B (zh) * | 2021-07-20 | 2022-09-02 | 电子科技大学 | 一种入睡过程动态连续分析方法 |
CN113397562A (zh) * | 2021-07-20 | 2021-09-17 | 电子科技大学 | 基于深度学习的睡眠纺锤波检测方法 |
CN113456030A (zh) * | 2021-08-05 | 2021-10-01 | 成都云卫康医疗科技有限公司 | 一种基于心率监测数据的睡眠分期方法 |
CN113662551A (zh) * | 2021-08-24 | 2021-11-19 | 深圳康佳电子科技有限公司 | 基于智能眼罩的睡眠监测处理方法、装置、智能眼罩 |
CN114707561B (zh) * | 2022-05-25 | 2022-09-30 | 清华大学深圳国际研究生院 | Psg数据自动分析方法、装置、计算机设备以及存储介质 |
CN116898455B (zh) * | 2023-07-06 | 2024-04-16 | 湖北大学 | 一种基于深度学习模型的睡眠脑电信号检测方法及*** |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107106028A (zh) * | 2014-12-18 | 2017-08-29 | 皇家飞利浦有限公司 | 用于心肺睡眠阶段分类的***和方法 |
CN107495962A (zh) * | 2017-09-18 | 2017-12-22 | 北京大学 | 一种单导联脑电的睡眠自动分期方法 |
CN108062756A (zh) * | 2018-01-29 | 2018-05-22 | 重庆理工大学 | 基于深度全卷积网络和条件随机场的图像语义分割方法 |
CN109602410A (zh) * | 2018-11-16 | 2019-04-12 | 青岛真时科技有限公司 | 一种可穿戴设备及其脉搏监测方法 |
CN109864750A (zh) * | 2019-01-31 | 2019-06-11 | 华南理工大学 | 基于经颅刺激的精神状态评估与调节***及其工作方法 |
CN110801221A (zh) * | 2019-12-09 | 2020-02-18 | 中山大学 | 基于无监督特征学习的睡眠呼吸暂停片段检测方法及设备 |
CN110890155A (zh) * | 2019-11-25 | 2020-03-17 | 中国科学技术大学 | 一种基于导联注意力机制的多类心律失常检测方法 |
CN111091116A (zh) * | 2019-12-31 | 2020-05-01 | 华南师范大学 | 一种用于判断心律失常的信号处理方法及*** |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7509163B1 (en) * | 2007-09-28 | 2009-03-24 | International Business Machines Corporation | Method and system for subject-adaptive real-time sleep stage classification |
US10706329B2 (en) * | 2018-11-13 | 2020-07-07 | CurieAI, Inc. | Methods for explainability of deep-learning models |
-
2020
- 2020-06-24 CN CN202010591697.XA patent/CN111631688B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107106028A (zh) * | 2014-12-18 | 2017-08-29 | 皇家飞利浦有限公司 | 用于心肺睡眠阶段分类的***和方法 |
CN107495962A (zh) * | 2017-09-18 | 2017-12-22 | 北京大学 | 一种单导联脑电的睡眠自动分期方法 |
CN108062756A (zh) * | 2018-01-29 | 2018-05-22 | 重庆理工大学 | 基于深度全卷积网络和条件随机场的图像语义分割方法 |
CN109602410A (zh) * | 2018-11-16 | 2019-04-12 | 青岛真时科技有限公司 | 一种可穿戴设备及其脉搏监测方法 |
CN109864750A (zh) * | 2019-01-31 | 2019-06-11 | 华南理工大学 | 基于经颅刺激的精神状态评估与调节***及其工作方法 |
CN110890155A (zh) * | 2019-11-25 | 2020-03-17 | 中国科学技术大学 | 一种基于导联注意力机制的多类心律失常检测方法 |
CN110801221A (zh) * | 2019-12-09 | 2020-02-18 | 中山大学 | 基于无监督特征学习的睡眠呼吸暂停片段检测方法及设备 |
CN111091116A (zh) * | 2019-12-31 | 2020-05-01 | 华南师范大学 | 一种用于判断心律失常的信号处理方法及*** |
Non-Patent Citations (2)
Title |
---|
基于BiGRU深度神经网络的心肌梗死检测;张行进 等;《计算机应用与软件》;20200229;第37卷(第2期);第48-52页 * |
睡眠分期算法研究;何垣谛;《中国优秀硕士学位论文全文数据集(电子期刊) 医药卫生科技辑》;20200131;第E060-1018页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111631688A (zh) | 2020-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111631688B (zh) | 睡眠自动分期的算法 | |
Jindal et al. | An adaptive deep learning approach for PPG-based identification | |
CN112656427B (zh) | 一种基于维度模型的脑电信号情绪识别方法 | |
CN108776788B (zh) | 一种基于脑电波的识别方法 | |
CN110801221B (zh) | 基于无监督特征学习的睡眠呼吸暂停片段检测设备 | |
CN109497996B (zh) | 一种微状态eeg时域特征的复杂网络构建及分析方法 | |
CN108968915A (zh) | 基于熵特征及支持向量机的睡眠状态分类方法及*** | |
Li et al. | Robust ECG biometrics using GNMF and sparse representation | |
CN108596069A (zh) | 基于深度3d残差网络的新生儿疼痛表情识别方法及*** | |
CN111407243B (zh) | 一种基于深度学习的脉搏信号压力识别方法 | |
CN110584596B (zh) | 基于双输入卷积神经网络的睡眠阶段分类方法及应用 | |
CN112509696A (zh) | 基于卷积自编码器高斯混合模型的健康数据检测方法 | |
De Falco et al. | Deep neural network hyper-parameter setting for classification of obstructive sleep apnea episodes | |
CN114732409A (zh) | 一种基于脑电信号的情绪识别方法 | |
CN115530847A (zh) | 一种基于多尺度注意力的脑电信号自动睡眠分期方法 | |
Wang et al. | Identification of Depression with a Semi-supervised GCN based on EEG Data | |
Liang et al. | Obstructive sleep apnea detection using combination of CNN and LSTM techniques | |
CN113303770A (zh) | 睡眠分期方法及装置 | |
Liu et al. | Automatic sleep arousals detection from polysomnography using multi-convolution neural network and random forest | |
CN114145745B (zh) | 基于图的多任务自监督情绪识别方法 | |
Guan | Application of logistic regression algorithm in the diagnosis of expression disorder in Parkinson's disease | |
Liao et al. | Recognizing diseases with multivariate physiological signals by a DeepCNN-LSTM network | |
Pimentel et al. | Human mental state monitoring in the wild: Are we better off with deeperneural networks or improved input features? | |
Li et al. | Tfformer: A time frequency information fusion based cnn-transformer model for osa detection with single-lead ecg | |
CN115736840A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20211029 |