CN107715276A - 闭环路径中睡眠状态反馈的声音睡眠控制系统及其方法 - Google Patents
闭环路径中睡眠状态反馈的声音睡眠控制系统及其方法 Download PDFInfo
- Publication number
- CN107715276A CN107715276A CN201711191664.0A CN201711191664A CN107715276A CN 107715276 A CN107715276 A CN 107715276A CN 201711191664 A CN201711191664 A CN 201711191664A CN 107715276 A CN107715276 A CN 107715276A
- Authority
- CN
- China
- Prior art keywords
- mrow
- sleep
- sound
- msup
- controller
- 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.)
- Pending
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
- A61M21/00—Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis
- A61M21/02—Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis for inducing sleep or relaxation, e.g. by direct nerve stimulation, hypnosis, analgesia
-
- 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]
-
- 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
- A61B5/374—Detecting the frequency distribution of signals, e.g. detecting delta, theta, alpha, beta or gamma waves
-
- 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/4836—Diagnosis combined with treatment in closed-loop systems or methods
-
- 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/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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
- A61M21/00—Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis
- A61M2021/0005—Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis by the use of a particular sense, or stimulus
- A61M2021/0027—Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis by the use of a particular sense, or stimulus by the hearing sense
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
- A61M2230/00—Measuring parameters of the user
- A61M2230/08—Other bio-electrical signals
- A61M2230/10—Electroencephalographic signals
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Psychiatry (AREA)
- Artificial Intelligence (AREA)
- Psychology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Anesthesiology (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Pain & Pain Management (AREA)
- Acoustics & Sound (AREA)
- Hematology (AREA)
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
Abstract
本发明涉及一种闭环路径中睡眠状态反馈的声音睡眠控制系统及其方法,包括传感器、控制器和执行器;传感器为脑波检测装置,实时采集人体睡眠的EEG信号,发送至控制器;控制器为计算机,根据EEG信号辨识出人体睡眠所处时期,生成与下阶段睡眠频谱相似的声音模型,发送至执行器;执行器为催眠声音播放器,播放助眠声音。本发明基于反馈控制思想建立以人体睡眠状态为控制对象,并以脑电信号为反馈量的睡眠控制系统,在研究睡眠过程脑电信号频率特征的基础上根据生物共振原理构建控制声音库和睡眠控制规则。本发明基于睡眠控制规则的睡眠控制系统能够渐进式地控制人体进入深度睡眠,且睡眠质量得到显著改善。
Description
技术领域
本发明涉及人体睡眠控制领域,具体涉及一种闭环路径中睡眠状态反馈的声音睡眠控制系统及其方法。
背景技术
随着社会生活节奏的不断加快,竞争日趋激烈,生存压力的不断增大,使得人们的睡眠质量开始出现问题,根据公布的第3版国际睡眠障碍分类可知,睡眠疾病已由之前的92种上升为102种,其中失眠症状仍占较大比重,据调研我国有38.6%的人群患有不同程度的睡眠疾病,老年群体睡眠疾病的发病率高达56%。失眠患者的身体可能会出现众多不适的症状,长时间的失眠会导致各种机体功能的下降甚至消失,包括注意力不集中、语言功能低下、记忆力变差以及反应迟钝等,严重者可能会导致某些人体功能性疾病。
因此,有效治疗失眠方法的研究迫在眉睫。目前,国内外主要的治疗方法有心理疗法、药物疗法、食物疗法及物理疗法等。心理疗法疗效不稳定,影响因素多,只能起到辅助治疗的作用;药物治疗,虽然疗效较好,但长期服用药物容易产生耐受性、依赖性及戒断性反应以及副作用;与传统治疗方法相比,新兴的物理治疗,如声音疗法、微电流刺激疗法、电磁刺激诱导睡眠疗法等,具有安全、副作用小等优点,因此,睡眠物理疗法的研究成为当前治疗睡障碍眠的主要方向。其中,声音疗法因其易实现、低成本是物理疗法中最有活力的一支。
然而,包括声音疗法在内的整个物理治疗睡眠方法,研究尚处于起步阶段,对于睡眠障碍的治疗效果远未达到预期目标,而近期研究又出现了驻足不前的现象,其原因在于:
第一,研究者多是从应用角度出发,重点探索何种物理方法改善睡眠的效果更突出,很少从理论的角度出发,研究声音等物理介质影响睡眠的机理。因此物理方法的研究仅靠实验摸索前行而缺少方向性和发展后劲;
第二,睡眠作为大脑活动的一部分,是一个复杂的生理过程,包含了若干个独立的睡眠状态,在整个睡眠过程中,大脑交替或循环地处在这几个睡眠状态之间,因此,在未知大脑当前状态的情况下盲目地施加物理催眠作用,效果难以保证;
第三,由于人脑睡眠过程是一个多分支渐进的动态过程,目前的睡眠方法不能按照一定的睡眠“路径”对睡眠过程实施控制,而简单地把该过程臆断为“清醒”和“睡眠”2个状态的潜移。
因此,对于以上的国内外人体睡眠治疗研究的现状,在人体睡眠方面的控制研究发展还有很大的改进与再设计空间。
发明内容
为解决现有技术中的上述不足,弥补人体睡眠控制领域的空白,本发明提供一种闭环路径中睡眠状态反馈的声音睡眠控制系统及其方法,利用声音控制,渐进式的诱导人体进入深度睡眠。
本发明所采用的技术方案为:
闭环路径中睡眠状态反馈的声音睡眠控制系统,其特征在于:
包括传感器、控制器和执行器;
传感器为脑波检测装置,实时采集人体睡眠的EEG信号,发送至控制器;
控制器为计算机,根据EEG信号辨识出人体睡眠所处时期,生成与下阶段睡眠频谱相似的声音模型,发送至执行器;
执行器为催眠声音播放器,播放助眠声音。
执行器以DSP为核心,通过功率放大电路驱动EDIFIER R12U多媒体有源音箱发出助眠声音。
闭环路径中睡眠状态反馈的声音睡眠控制方法,其特征在于:
构建由传感器、控制器和执行器组成的声音睡眠控制系统;
传感器实时采集人体睡眠的EEG信号,控制器通过辨识EEG信号判断人体的睡眠状态,给执行器提供构造催眠声音的数学模型,执行器根据该模型生成相应的控制声音,逐步诱导人体进入深度睡眠期。
通过传感器采集的EEG信号,控制器分析得到人体睡眠状态,并进行睡眠控制路径的规划,得到睡眠路径表,通过在不同时期中播放相应的下一时期的催眠声音,将人体睡眠分步式的按照一定的睡眠过程进行诱导。
控制器对人体睡眠过程中所处的不同睡眠状态特征量进行分析后,得到相应的睡眠特征量频谱,根据睡眠路径将睡眠状态排序后,得到的在不同睡眠时期所应使用的睡眠特征诱导量,形成睡眠控制规则表,从而形成睡眠控制规则。
控制器根据人体睡眠中各个时期的特征量的不同,由不同的脑电波的数学模型,分别生成相应的控制声音,形成睡眠声音库。
控制器对睡眠状态的辨识包括:
(1)EEG信号预处理:
首先对EEG信号进行滤波以消除其中的噪声成分,采用具有最大平坦幅度特性的5阶Butterworth滤波器对EEG信号进行0.5~32Hz的带通滤波,其传递函数为:
式中,ωc为滤波器截止频率;
(2)特征量提取:
以EEG信号的β、a、θ、δ4种短时功率谱为特征量,计算步骤为:
①将采集到的EEG信号分成10秒的数据片段x(n)(n=0,1,2,…,L-1),其中,L为x(n)的长度;
②利用Hamming窗作为观测窗,其函数为:
式中,取a=0.46,代入上式可得:
式中,窗函数的长度M=1280;
③对每个窗内的数据做DFT变换,求取X(m,k)
式中,N为窗函数在x(n)时间方向上每隔N点的一次移动,一共做了L/N个M点的DFT变换;
④计算短时功率谱:
式中,E(m)表示信号在m点的短时能量;
⑤根据不同特征波段的频率分布,获取各个特征波段的能量,从而得到睡眠辨识的主要特征量,β、a、θ、δ4种短时功率谱;
(3)睡眠状态辨识:
将Eδ特征量输入到由径向基神经网络构成的辨识器中,通过算法能够准确辨识出人体睡眠时期;
经过睡眠辨识出相应的睡眠状态后,由相应的算法得到睡眠的特征量,从而在声音生成模块库中,按照睡眠路径播放下一阶段的催眠声音,根据构建的睡眠控制规则表,渐进式的将人体睡眠诱导入深睡状态。
本发明具有以下优点:
本发明提供了一种缓解人体失眠症状的闭环反馈声音睡眠控制系统,弥补了人体睡眠控制领域的空白,利用人体睡眠EEG信号中的特征量,其在各个睡眠时期所表现出来能量比例的差异,按照各个特征量的频谱从而生成相应的声音模型,由浅入深的分步式按照睡眠“路径”诱导人体睡眠,从而达到深度睡眠的效果。该系统分步式的控制增强了对睡眠诱导的目的性和有效性,提高了入睡效率,具有安全、副作用小等优点,并且,声音疗法在睡眠物理疗法中也比较容易实现,成本较低。
附图说明
图1为睡眠控制系统框图。
图2为睡眠路径控制思路图。
图3为脑电采集电极分布图。
图4为不同睡眠特征量生成声音模型频谱图(α波)。
图5为不同睡眠特征量生成声音模型频谱图(θ波)。
图6为不同睡眠特征量生成声音模型频谱图(高δ波)。
图7为不同睡眠特征量生成声音模型频谱图(低δ波)。
图8为粉红噪声的功率谱图。
具体实施方式
下面结合具体实施方式对本发明进行详细的说明。
本发明涉及的一种闭环路径中睡眠状态反馈的声音睡眠控制系统,包括传感器、控制器和执行器。
传感器为脑波检测装置,实时采集人体睡眠的EEG信号,发送至控制器。
控制器为计算机,根据EEG信号辨识出人体睡眠所处时期,生成与下阶段睡眠频谱相似的声音模型,发送至执行器。
执行器为催眠声音播放器,播放助眠声音,执行器以DSP为核心,通过功率放大电路驱动EDIFIER R12U多媒体有源音箱发出助眠声音。
控制器通过辨识EEG信号判断人体的睡眠状态,并根据睡眠控制规则,给执行器提供构造催眠声音的数学模型,执行器根据该模型生成相应的控制声音,逐步诱导人体进入深度睡眠期;执行器是由一个以DSP为核心的智能装置,在CCS3.3开发平台及OpenViBe平台上完成了主要算法的代码编写及调试,使其能够根据各种声音的数学模型生成声音信号,并通过功率放大电路驱动EDIFIER R12U多媒体有源音箱发出控制声音。
根据脑电采集装置所采集并且分析得到的人体睡眠状态进行睡眠控制路径的规划,得到睡眠路径表,通过在不同时期中播放相应的下一时期的催眠声音,将人体睡眠分步式的按照一定的睡眠过程进行诱导。
对人体睡眠过程中所处的不同睡眠状态特征量进行分析后,得到相应的睡眠特征量频谱,而根据一定的睡眠路径将睡眠状态排序后,得到的在不同睡眠时期所应使用的睡眠特征诱导量,形成睡眠控制规则表,从而形成睡眠控制规则。
控制器内包含有声音生成模块,即根据人体睡眠中各个时期的特征量的不同,由不同的脑电波的数学模型,分别生成相应的控制声音,形成睡眠声音库。
如图1所示,本发明的整体组成包括了作为控制对象的人体睡眠状态,作为传感器的脑电检测装置,作为控制器的计算机控制算法和作为执行器的催眠声音播放装置。其中,控制器包括了睡眠状态辨识模块,睡眠状态控制算法模块和催眠声音模型生成模块。
睡眠状态辨识:
(1)EEG信号预处理。首先对EEG信号进行滤波以消除其中的噪声成分,采用具有最大平坦幅度特性的5阶Butterworth滤波器对EEG信号进行0.5~32Hz的带通滤波,其传递函数为:
式中,ωc为滤波器截止频率。
(2)特征量提取。以EEG信号的β、a、θ、δ等4种短时功率谱为特征量,计算步骤为:
①将采集到的EEG信号分成10秒的数据片段x(n)(n=0,1,2,…,L-1),其中,L为x(n)的长度;
②利用Hamming窗作为观测窗,其函数为:
式中,取a=0.46,代入上式可得:
式中,窗函数的长度M=1280;
③对每个窗内的数据做DFT变换,求取X(m,k)
式中,N为窗函数在x(n)时间方向上每隔N点的一次移动,一共做了L/N个M点的DFT变换。
④计算短时功率谱
式中,E(m)表示信号在m点的短时能量。
⑤根据不同特征波段的频率分布,获取各个特征波段的能量,从而得到睡眠辨识的主要特征量,β、a、θ、δ等4种短时功率谱;
(3)睡眠状态辨识。将Eδ等特征量输入到由径向基神经网络构成的辨识器中,通过算法能够准确辨识出人体睡眠时期。
经过睡眠辨识出相应的睡眠状态后,反馈给计算机,由相应的算法得到睡眠的特征量,从而在声音生成模块库中,按照睡眠“路径”播放下一阶段的催眠声音,根据构建的睡眠控制规则表,渐进式的将人体睡眠诱导入深睡状态。
如图2所示,本发明控制人体按照一定的睡眠“路径”渐进式地进入深睡循环。具体来说,系统通过EEG辨识出人体在当前睡眠进程中所处的睡眠状态,并根据睡眠状态发出相应的控制声音。例如,当人体处于REM状态,则施加浅睡催眠声音诱导人体进入浅睡状态,再如,当人体处于深睡状态,则施加深度睡眠声音,阻滞睡眠离开深睡期进入REM期的进程,延长深睡时间。具体睡眠过程控制如图2。
如图3所示,处理EEG首先要实现EEG信号的记录,按照电极放置位置的不同,EEG信号的记录方法可分为深部EEG信号、皮层EEG信号以及头皮EEG信号记录方法,其中,由于头皮EEG信号记录方法中电极放置的位置在头皮表面,因此易于操作而被广泛采用。EEG导联标准虽然到目前为止还没有统一,各个厂家都按照自己的方案设置一些固定的EEG导联,但EEG导联均采用国际联合会规定的10-20系统电极法,这种系统电极方法采用解剖学标志位作为EEG电极的标准安放部位。该方法将电极放置在头部19个标准位置和两个参考位置,图3为电极位置示意图,且它的电极放置应遵循一定的基本原则。本发明主要采用位于头部左、右额,左、右枕的对睡眠作用最相关的四路通道进行分析,充分保证了使用者的舒适性,不会加深睡眠障碍者的失眠影响。
如图4-7所示,作为一种非平稳随机信号,EEG的频率并非均匀分布,而是主要集中在0.5-30Hz几个不连续的区域,通常根据EEG信号主要频率集中的区域,将EEG进行以下分类,本发明根据α波、θ波、δ波等脑电波的数学模型,分别生成α波、θ波、高δ波和低δ波4种声音,通过对这4种声音的功率谱分析,可见,这四种波大体均匀地分布在整个脑电波的频率范围内,其中α波声音、θ波声音和δ波声音分别分布在高、中和低频段,而高、低δ波声音使得该声音在低频段分解地更细。根据生物共振原理,这4种声音分别可以引起人体大脑在各频段的活动加剧或衰减。
如图8所示,粉红噪声是自然界中最常见的一种噪声,上世纪80年代日本学家利用睡眠中的共振现象,进行了粉红噪声对睡眠的影响实验,实验结果说明了粉红噪声具有助眠和延长深度睡眠的作用。图8是粉红噪声的功率谱图,对比图8和深度睡眠过程功率谱图可见粉红噪声与人体深睡状态下的频谱具有相同的分布。根据生物共振原理,粉红噪声声音可能具有诱导人体进入或滞留于深睡状态的功能。
因此,本发明将采用α波,θ波,高、低δ波声音及粉红噪声5种声音作为睡眠控制声音,初步构成睡眠控制声音库。之后,按照睡眠“路径”和睡眠原理分别播放不同的声音,达到控制睡眠的效果。
本发明的内容不限于实施例所列举,本领域普通技术人员通过阅读本发明说明书而对本发明技术方案采取的任何等效的变换,均为本发明的权利要求所涵盖。
Claims (7)
1.闭环路径中睡眠状态反馈的声音睡眠控制系统,其特征在于:
包括传感器、控制器和执行器;
传感器为脑波检测装置,实时采集人体睡眠的EEG信号,发送至控制器;
控制器为计算机,根据EEG信号辨识出人体睡眠所处时期,生成与下阶段睡眠频谱相似的声音模型,发送至执行器;
执行器为催眠声音播放器,播放助眠声音。
2.根据权利要求1所述的闭环路径中睡眠状态反馈的声音睡眠控制系统,其特征在于:
执行器以DSP为核心,通过功率放大电路驱动EDIFIER R12U多媒体有源音箱发出助眠声音。
3.闭环路径中睡眠状态反馈的声音睡眠控制方法,其特征在于:
构建由传感器、控制器和执行器组成的声音睡眠控制系统;
传感器实时采集人体睡眠的EEG信号,控制器通过辨识EEG信号判断人体的睡眠状态,给执行器提供构造催眠声音的数学模型,执行器根据该模型生成相应的控制声音,逐步诱导人体进入深度睡眠期。
4.根据权利要求3所述的闭环路径中睡眠状态反馈的声音睡眠控制方法,其特征在于:
通过传感器采集的EEG信号,控制器分析得到人体睡眠状态,并进行睡眠控制路径的规划,得到睡眠路径表,通过在不同时期中播放相应的下一时期的催眠声音,将人体睡眠分步式的按照一定的睡眠过程进行诱导。
5.根据权利要求3所述的闭环路径中睡眠状态反馈的声音睡眠控制方法,其特征在于:
控制器对人体睡眠过程中所处的不同睡眠状态特征量进行分析后,得到相应的睡眠特征量频谱,根据睡眠路径将睡眠状态排序后,得到的在不同睡眠时期所应使用的睡眠特征诱导量,形成睡眠控制规则表,从而形成睡眠控制规则。
6.根据权利要求3所述的闭环路径中睡眠状态反馈的声音睡眠控制方法,其特征在于:
控制器根据人体睡眠中各个时期的特征量的不同,由不同的脑电波的数学模型,分别生成相应的控制声音,形成睡眠声音库。
7.根据权利要求3所述的闭环路径中睡眠状态反馈的声音睡眠控制方法,其特征在于:
控制器对睡眠状态的辨识包括:
(1)EEG信号预处理:
首先对EEG信号进行滤波以消除其中的噪声成分,采用具有最大平坦幅度特性的5阶Butterworth滤波器对EEG信号进行0.5~32Hz的带通滤波,其传递函数为:
<mrow>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>&omega;</mi>
<mi>c</mi>
<mn>2</mn>
</msubsup>
<mrow>
<msup>
<mi>s</mi>
<mn>5</mn>
</msup>
<mo>+</mo>
<mn>3.261</mn>
<msup>
<mi>s</mi>
<mn>4</mn>
</msup>
<msub>
<mi>&omega;</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<mn>5.2361</mn>
<msup>
<mi>s</mi>
<mn>3</mn>
</msup>
<msubsup>
<mi>&omega;</mi>
<mi>c</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mn>5.2361</mn>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>&omega;</mi>
<mi>c</mi>
<mn>3</mn>
</msubsup>
<mo>+</mo>
<mn>3.2361</mn>
<msubsup>
<mi>s&omega;</mi>
<mi>c</mi>
<mn>4</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&omega;</mi>
<mi>c</mi>
<mn>5</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
式中,ωc为滤波器截止频率;
(2)特征量提取:
以EEG信号的β、a、θ、δ 4种短时功率谱为特征量,计算步骤为:
①将采集到的EEG信号分成10秒的数据片段x(n)(n=0,1,2,…,L-1),其中,L为x(n)的长度;
②利用Hamming窗作为观测窗,其函数为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&omega;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>&alpha;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>k</mi>
</mrow>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>&le;</mo>
<mi>k</mi>
<mo>&le;</mo>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,取a=0.46,代入上式可得:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&omega;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0.54</mn>
<mo>-</mo>
<mn>0.46</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>k</mi>
</mrow>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>&le;</mo>
<mi>k</mi>
<mo>&le;</mo>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,窗函数的长度M=1280;
③对每个窗内的数据做DFT变换,求取X(m,k)
<mrow>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mi>&omega;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>m</mi>
<mi>N</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>k</mi>
<mi>&pi;</mi>
</mrow>
<mi>M</mi>
</mfrac>
<mi>n</mi>
</mrow>
</msup>
</mrow>
式中,N为窗函数在x(n)时间方向上每隔N点的一次移动,一共做了L/N个M点的DFT变换;
④计算短时功率谱:
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
<mi>L</mi>
</mfrac>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>m</mi>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
式中,E(m)表示信号在m点的短时能量;
⑤根据不同特征波段的频率分布,获取各个特征波段的能量,从而得到睡眠辨识的主要特征量,β、a、θ、δ 4种短时功率谱;
<mrow>
<msub>
<mi>E</mi>
<mi>&delta;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>0.5</mn>
</mrow>
<mn>3</mn>
</munderover>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>E</mi>
<mi>&theta;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>4</mn>
</mrow>
<mn>7</mn>
</munderover>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>E</mi>
<mi>&alpha;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>8</mn>
</mrow>
<mn>13</mn>
</munderover>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mi>&beta;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>14</mn>
</mrow>
<mn>30</mn>
</munderover>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
(3)睡眠状态辨识:
将Eδ特征量输入到由径向基神经网络构成的辨识器中,通过算法能够准确辨识出人体睡眠时期;
经过睡眠辨识出相应的睡眠状态后,由相应的算法得到睡眠的特征量,从而在声音生成模块库中,按照睡眠路径播放下一阶段的催眠声音,根据构建的睡眠控制规则表,渐进式的将人体睡眠诱导入深睡状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711191664.0A CN107715276A (zh) | 2017-11-24 | 2017-11-24 | 闭环路径中睡眠状态反馈的声音睡眠控制系统及其方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711191664.0A CN107715276A (zh) | 2017-11-24 | 2017-11-24 | 闭环路径中睡眠状态反馈的声音睡眠控制系统及其方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107715276A true CN107715276A (zh) | 2018-02-23 |
Family
ID=61219214
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711191664.0A Pending CN107715276A (zh) | 2017-11-24 | 2017-11-24 | 闭环路径中睡眠状态反馈的声音睡眠控制系统及其方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107715276A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108572557A (zh) * | 2018-04-13 | 2018-09-25 | 深圳迈睿智能科技有限公司 | 睡眠环境的控制系统和控制方法 |
CN109316170A (zh) * | 2018-11-16 | 2019-02-12 | 武汉理工大学 | 基于深度学习的脑电波辅助睡眠及唤醒系统 |
CN110841169A (zh) * | 2019-11-28 | 2020-02-28 | 中国科学院深圳先进技术研究院 | 一种用于睡眠调节的深度声音刺激系统和方法 |
CN111231870A (zh) * | 2018-11-28 | 2020-06-05 | 上海博泰悦臻电子设备制造有限公司 | 车辆睡眠模式切换方法及装置 |
CN111760160A (zh) * | 2020-07-01 | 2020-10-13 | 北京脑陆科技有限公司 | 一种基于eeg信号的助眠方法 |
CN113195031A (zh) * | 2018-12-17 | 2021-07-30 | 皇家飞利浦有限公司 | 用于递送音频输出的系统和方法 |
WO2021258245A1 (zh) * | 2020-06-22 | 2021-12-30 | 华为技术有限公司 | 更新助眠音频信号的方法及装置 |
CN114247027A (zh) * | 2020-09-24 | 2022-03-29 | 原相科技股份有限公司 | 以声音诱导睡眠的方法及执行该方法的声音播放装置 |
CN115989998A (zh) * | 2022-11-22 | 2023-04-21 | 常州瑞神安医疗器械有限公司 | 一种检测帕金森病患者睡眠阶段的方法 |
CN116150591A (zh) * | 2023-04-04 | 2023-05-23 | 安徽星辰智跃科技有限责任公司 | 一种睡眠感官应激水平检测量化的方法、系统和装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105749402A (zh) * | 2016-02-24 | 2016-07-13 | 陕西科技大学 | 一种基于反馈机理的深度催眠系统及方法 |
-
2017
- 2017-11-24 CN CN201711191664.0A patent/CN107715276A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105749402A (zh) * | 2016-02-24 | 2016-07-13 | 陕西科技大学 | 一种基于反馈机理的深度催眠系统及方法 |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108572557A (zh) * | 2018-04-13 | 2018-09-25 | 深圳迈睿智能科技有限公司 | 睡眠环境的控制系统和控制方法 |
CN109316170B (zh) * | 2018-11-16 | 2020-12-22 | 武汉理工大学 | 基于深度学习的脑电波辅助睡眠及唤醒系统 |
CN109316170A (zh) * | 2018-11-16 | 2019-02-12 | 武汉理工大学 | 基于深度学习的脑电波辅助睡眠及唤醒系统 |
CN111231870A (zh) * | 2018-11-28 | 2020-06-05 | 上海博泰悦臻电子设备制造有限公司 | 车辆睡眠模式切换方法及装置 |
CN113195031B (zh) * | 2018-12-17 | 2024-04-12 | 皇家飞利浦有限公司 | 用于递送音频输出的系统和方法 |
CN113195031A (zh) * | 2018-12-17 | 2021-07-30 | 皇家飞利浦有限公司 | 用于递送音频输出的系统和方法 |
CN110841169B (zh) * | 2019-11-28 | 2020-09-25 | 中国科学院深圳先进技术研究院 | 一种用于睡眠调节的深度学习声音刺激系统和方法 |
WO2021103084A1 (zh) * | 2019-11-28 | 2021-06-03 | 中国科学院深圳先进技术研究院 | 一种用于睡眠调节的深度声音刺激系统和方法 |
CN110841169A (zh) * | 2019-11-28 | 2020-02-28 | 中国科学院深圳先进技术研究院 | 一种用于睡眠调节的深度声音刺激系统和方法 |
WO2021258245A1 (zh) * | 2020-06-22 | 2021-12-30 | 华为技术有限公司 | 更新助眠音频信号的方法及装置 |
CN111760160A (zh) * | 2020-07-01 | 2020-10-13 | 北京脑陆科技有限公司 | 一种基于eeg信号的助眠方法 |
CN114247027A (zh) * | 2020-09-24 | 2022-03-29 | 原相科技股份有限公司 | 以声音诱导睡眠的方法及执行该方法的声音播放装置 |
CN114247027B (zh) * | 2020-09-24 | 2024-05-03 | 达发科技股份有限公司 | 以声音诱导睡眠的方法及执行该方法的声音播放装置 |
CN115989998A (zh) * | 2022-11-22 | 2023-04-21 | 常州瑞神安医疗器械有限公司 | 一种检测帕金森病患者睡眠阶段的方法 |
CN115989998B (zh) * | 2022-11-22 | 2023-11-14 | 常州瑞神安医疗器械有限公司 | 一种检测帕金森病患者睡眠阶段的方法 |
CN116150591A (zh) * | 2023-04-04 | 2023-05-23 | 安徽星辰智跃科技有限责任公司 | 一种睡眠感官应激水平检测量化的方法、系统和装置 |
CN116150591B (zh) * | 2023-04-04 | 2023-06-20 | 安徽星辰智跃科技有限责任公司 | 一种睡眠感官应激水平检测量化的方法、系统和装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107715276A (zh) | 闭环路径中睡眠状态反馈的声音睡眠控制系统及其方法 | |
CN109224242B (zh) | 基于vr交互的心理放松系统及方法 | |
CN106725462B (zh) | 基于脑电信号的声光睡眠干预系统和方法 | |
CN104825159A (zh) | 脑电同频检测方法、脑电同频刺激方法与装置 | |
CN105771087A (zh) | 基于音乐和肌电反馈刺激的康复训练系统 | |
CN103372258A (zh) | 失眠治疗仪及失眠治疗方法 | |
CN111569226B (zh) | 基于音频处理的助眠方法及装置 | |
CN107802938B (zh) | 一种用于镇痛的音乐电刺激的产生方法 | |
CN107463249A (zh) | 基于vr头显视觉诱发电位的脑机接口系统及控制方法 | |
KR20200075711A (ko) | 인공지능 수면개선 비침습적 뇌회로 조절치료시스템 및 방법 | |
CN101927052A (zh) | 一种脑电波生物反馈戒毒仪 | |
Krugliakova et al. | Boosting recovery during sleep by means of auditory stimulation | |
Makada et al. | Enhancing memory retention by increasing alpha and decreasing beta brainwaves using music | |
CN107305597A (zh) | 一种基于大数据分析的心理自我治愈系统 | |
Rakoviĉ et al. | An overview of microwave resonance therapy and EEG correlates of microwave resonance relaxation and other consciousness altering techniques | |
CN103313748B (zh) | 个性化的治愈声音数据库 | |
CN116312531A (zh) | 一种基于脑电信号的语音谱图重构方法 | |
CN106964069B (zh) | 一种具有生物反馈功能的微循环修复系统 | |
Bouafif | Monitoring of anesthesia by bispectral analysis of EEG signals | |
TWI543749B (zh) | 可調節腦內不同神經迴路之雙耳波差導入音樂方法 | |
CN104368073A (zh) | 一种睡眠调理方法 | |
Drozdenko et al. | The Influence of a Low-Frequency Musical Fragment on the Neural Oscillations | |
Indursky et al. | Low-frequency rhythmic electrocutaneous hand stimulation during slow-wave night sleep: Physiological and therapeutic effects | |
CN108309287B (zh) | 一种早期筛查大脑功能窄化的物理检测方法 | |
Rakovic et al. | Electroencephalographic (EEG) correlates of some activities which may alter consciousness: the transcendental meditation technique, musicogenic states, microwave resonance relaxation, healer/healee interaction, and alertness/drowsiness |
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 |