CN106580319A - 一种基于小波变换的脑电放松度识别方法及装置 - Google Patents
一种基于小波变换的脑电放松度识别方法及装置 Download PDFInfo
- Publication number
- CN106580319A CN106580319A CN201611184948.2A CN201611184948A CN106580319A CN 106580319 A CN106580319 A CN 106580319A CN 201611184948 A CN201611184948 A CN 201611184948A CN 106580319 A CN106580319 A CN 106580319A
- Authority
- CN
- China
- Prior art keywords
- brain
- energy
- electric array
- array signal
- eeg signals
- 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
- 230000009466 transformation Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 33
- 210000004556 brain Anatomy 0.000 claims abstract description 239
- 238000000605 extraction Methods 0.000 claims abstract description 15
- 238000005070 sampling Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 230000005611 electricity Effects 0.000 claims description 68
- 238000009826 distribution Methods 0.000 claims description 18
- 238000000354 decomposition reaction Methods 0.000 claims description 17
- 238000005315 distribution function Methods 0.000 claims description 12
- 238000005096 rolling process Methods 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 238000002592 echocardiography Methods 0.000 claims description 9
- 230000009467 reduction Effects 0.000 claims description 9
- 238000012549 training Methods 0.000 claims description 9
- 238000013528 artificial neural network Methods 0.000 claims description 7
- 238000005457 optimization Methods 0.000 claims description 6
- 230000001360 synchronised effect Effects 0.000 claims description 5
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 230000003321 amplification Effects 0.000 abstract description 3
- 238000003199 nucleic acid amplification method Methods 0.000 abstract description 3
- 239000010410 layer Substances 0.000 description 18
- 230000006870 function Effects 0.000 description 12
- 230000008569 process Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 230000015654 memory Effects 0.000 description 5
- 239000000470 constituent Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 210000003205 muscle Anatomy 0.000 description 3
- 238000000513 principal component analysis Methods 0.000 description 3
- 230000001755 vocal effect Effects 0.000 description 3
- 208000020401 Depressive disease Diseases 0.000 description 2
- 208000013738 Sleep Initiation and Maintenance disease Diseases 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 206010022437 insomnia Diseases 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- 206010000117 Abnormal behaviour Diseases 0.000 description 1
- 208000019901 Anxiety disease Diseases 0.000 description 1
- 206010003830 Automatism Diseases 0.000 description 1
- 206010019233 Headaches Diseases 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 238000003723 Smelting Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000000540 analysis of variance Methods 0.000 description 1
- 230000036506 anxiety Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000013542 behavioral therapy Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 230000001149 cognitive effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 206010016256 fatigue Diseases 0.000 description 1
- 231100000869 headache Toxicity 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 230000007787 long-term memory Effects 0.000 description 1
- 230000003340 mental effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000003534 oscillatory effect Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000002040 relaxant effect Effects 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 230000003860 sleep quality Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 208000011580 syndromic disease Diseases 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- 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/375—Electroencephalography [EEG] using biofeedback
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
- A61B5/165—Evaluating the state of mind, e.g. depression, anxiety
-
- 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/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/7253—Details of waveform analysis characterised by using transforms
-
- 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
-
- 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
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Psychiatry (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychology (AREA)
- Anesthesiology (AREA)
- Mathematical Physics (AREA)
- Fuzzy Systems (AREA)
- Evolutionary Computation (AREA)
- Pain & Pain Management (AREA)
- Acoustics & Sound (AREA)
- Power Engineering (AREA)
- Hematology (AREA)
- Child & Adolescent Psychology (AREA)
- Developmental Disabilities (AREA)
- Educational Technology (AREA)
- Hospice & Palliative Care (AREA)
- Social Psychology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种基于小波变换的脑电放松度识别方法,包括如下步骤:对接收的第一脑电序列信号进行降采样,获得第二脑电序列信号;对所述第二脑电序列信号的各个时刻的脑电信号进行加权移动平均计算,得到第三脑电序列信号;对所述第三脑电序列信号进行小波变换,并根据每个脑电波的频率范围重构得到各个脑电波;基于重构得到的所述各个脑电波进行特征提取,并根据提取得到的特征量进行分类识别,得到脑电放松度。本发明还提供了一种基于小波变换的脑电放松度识别装置,可实现准确的脑电放松度识别。
Description
技术领域
本发明涉及放松治疗领域,尤其涉及一种基于小波变换的脑电放松度识别方法及装置。
背景技术
放松训练是行为疗法中使用最广的技术之一,是在心理学实验的基础上建立和发展起来的咨询和治疗方法,其在治疗焦虑抑郁症、神经性头痛、失眠、高血压病,减轻更年期综合征和转变不良行为模式等方面取得了较好的疗效。
现有的放松训练主要有录音指导、口头指导和生物反馈指导。其中,录音指导方法僵化、没有变化,无法根据受训者的状态变化内容;口头指导则要求对口头指导的对象要求很高,且受到时间、场地限制;生物反馈指导以脑电反馈为主,能够结合前两种方式的优点,因而受到广泛关注。
进行生物反馈指导需要识别用户的放松度,而计算放松度首先需要从用户的脑电信号中提取各个频段的脑电波(包括Delta、Theta、Alpha、Beta、Gamma波),但实际脑电信号极其微弱,且非常容易受到外界工频,磁场扰动,低频直流、舌电伪迹,出汗伪迹,眼电伪迹,脉搏伪迹以及肌电伪迹等因素的干扰。其中,工频干扰相对固定,较易去除;低频直流虽然可以通过滤波器滤除,但考虑到滤波器带来的频带衰减,直接通过滤波器精确滤除低频直流信息会对脑电信号产生影响;此外,在提取脑电波时,由于各个脑电波的频率比较接近,因此如何准确的分离各个脑电波也关系到最终放松度识别的准确度。
发明内容
针对上述问题,本发明的目的在于提供一种基于小波变换的脑电放松度识别方法及装置,可在不干扰脑电信号的情况下去除低频直流信号,并可准确的分离出脑电信号中的各个脑电波。
本发明提供了一种基于小波变换的脑电放松度识别方法,包括如下步骤:
对接收的第一脑电序列信号进行降采样,获得第二脑电序列信号;
对所述第二脑电序列信号的各个时刻的脑电信号进行加权移动平均计算,得到第三脑电序列信号;
对所述第三脑电序列信号进行小波变换,并根据每个脑电波的频率范围重构得到各个脑电波;
基于重构得到的所述各个脑电波进行特征提取,并根据提取得到的特征量进行分类识别,得到脑电放松度。
优选地,所述对所述第二脑电序列信号的各个时刻的脑电信号进行加权移动平均计算,得到第三脑电序列信号,具体包括:
基于当前的第j个时刻的脑电信号,获取所述第二脑电序列信号中的位于第(j-(N-1)/2)个时刻至第(j+(N-1)/2)个时刻之间的N个脑电信号的能量;其中,N为预设的影响数,且N为奇数,j为大于(N+1)/2的整数;
根据预设的权值分布函数为获取的N个脑电信号的能量分配权值;其中,N个脑电信号的能量的权值之和为1;
对所述N个脑电信号的能量根据分配的权值进行加权求和,得到新的第j个时刻的脑电信号的能量;
依次对所述第二脑电序列信号的各个时刻的脑电信号的能量进行加权求和后,根据所有时刻的新的脑电信号的能量,生成第三脑电序列信号。
优选地,所述对所述第二脑电序列信号的各个时刻的脑电信号进行加权移动平均计算,得到第三脑电序列信号,具体包括:
基于当前的第j个时刻的脑电信号,获取所述第二脑电序列信号中的位于第(j-N)个时刻至第(j-1)个时刻之间的N个脑电信号的能量;其中,N为预设的影响数,j为大于N的整数;
根据预设的权值分布函数为获取的N个脑电信号的能量分配权值;其中,N个脑电信号的能量的权值之和为1;
对所述N个脑电信号的能量根据分配的权值进行加权求和,得到新的第j个时刻的脑电信号的能量;
依次对所述第二脑电序列信号的各个时刻的脑电信号的能量进行加权求和后,根据的所有时刻的新的脑电信号的能量,生成第三脑电序列信号。
优选地,所述权值分布函数为正态分布函数。
优选地,所述对所述第三脑电序列信号进行小波变换,并根据每个脑电波的频率范围重构得到各个脑电波具体为:
根据香农-奈奎斯特采样原理及所述第三脑电序列信号的采样频率进行频段分层,计算得到每层的频率范围;
依据所述频段分层中每层的频率范围及每个脑电波的频率范围,确定与每个脑电波对应的小波分解和重构所需的层数;
根据与每个脑电波对应的所述小波分解所需的层数及预先选择的母小波进行信号分解,得到与每个脑电波对应的按频段划分的多层波形;
根据与所述小波重构所需的层数对应的系数及分解得到的与每个脑电波对应的所述多层波形,重构得到各个脑电波。
优选地,所述母小波为coif3小波,且所述coif3小波的中心频率-带宽比经小波熵自适应优化。
优选地,在基于加权移动平均算法对所述第二脑电序列信号的各个时刻的脑电信号进行处理,得到第三脑电序列信号之前,还包括:
以所述第三脑电序列信号为原始信号,以与所述第三脑电序列信号同步采集得到的伪迹序列信号为参考信号,采用经函数链神经网络优化的自适应滤波器对所述第三脑电序列信号进行滤波,得到去除伪迹序列信号后的第三脑电序列信号。
优选地,所述基于所述各个脑电波进行特征提取,并根据提取得到的特征量进行分类识别,得到当前的脑电放松度,具体包括:
计算各个脑电波的能量和能量分布密度,并基于各个脑电波的能量及每个脑电波的等电位线能量比率,计算得到每个脑电波的等电位线能量基线;
根据各个脑电波的能量及等电位线能量基线,统计能量小于所述电位线能量基线的点的个数;
根据各个脑电波的能量、能量分布密度、能量小于电位线能量基线的点的个数及第一脑电序列信号的等电位线能量,得到所述第一脑电序列信号的特征量;
基于预先训练好的分类器对所述特征量进行分类,识别得到与所述第一脑电序列信号对应的脑电放松度。
优选地,在根据各个脑电波的能量、能量分布密度、能量小于电位线能量基线的点的个数及第一脑电序列信号的等电位线能量,得到所述第一脑电序列信号的特征量之后,在基于预先训练好的分类器对所述特征量进行分类,识别得到与所述第一脑电序列信号对应的脑电放松度之前,还包括:
基于主成分分析法对所述特征量进行特征选择和特征降维。
本发明还提供了一种基于小波变换的脑电放松度识别装置,包括:
降采样单元,用于对接收的第一脑电序列信号进行降采样,获得第二脑电序列信号;
加权移动平均单元,用于对所述第二脑电序列信号的各个时刻的脑电信号进行加权移动平均计算,得到第三脑电序列信号;
小波变换单元,用于对所述第三脑电序列信号进行小波变换,并根据每个脑电波的频率范围重构得到各个脑电波;
放松度识别单元,用于基于重构得到的所述各个脑电波进行特征提取,并根据提取得到的特征量进行分类识别,得到脑电放松度。
本发明提供的基于小波变换的脑电放松度识别方法及装置,通过采用加权移动平均去除所述第二脑电序列信号中的低频直流信息,保证了在去低频直流信息的同时不会对所述第二脑袋序列信号产生任何影响,保证所述第二脑袋序列信号的平滑,同时采用小波变换的方法准确提取各个脑电波,保证了基于这些提取的脑电波计算得到的特征量的准确性,进而保证了最终的脑电放松度的准确识别,为准确的生物反馈指导提供了数据基础和依据。
附图说明
为了更清楚地说明本发明的技术方案,下面将对实施方式中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的基于小波变换的脑电放松度识别方法的流程示意图。
图2是通过切片得到第一脑电序列信号的示意图。
图3是本发明实施例提供的Shannon小波熵与中心频率-带宽比的关系图。
图4是本发明实施例提供的对对第二脑电序列信号进行加权移动平均计算的原理图。
图5是自适应滤波器的工作原理图。
图6是本发明实施例提供的基于小波变换的脑电放松度识别装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1,本发明实施例提供了一种基于小波变换的脑电放松度识别方法,其可包括如下步骤:
S101,对接收的第一脑电序列信号进行降采样,获得第二脑电序列信号。
在本发明实施例中,所述第一脑电序列信号可通过脑电电极采集获得,其中,一般地,脑电电极采集的原始脑电信号的持续时间比较长(如数小时甚至更长),因此进行需要对原始脑电信号进行切片,以得到所述第一脑电序列信号,例如,如图2所示,每个切片的片段为30秒,即每段所述第一脑电序列信号的长度为30秒。
在本发明实施例中,一般地,所述第一脑电序列信号的采样频率为500HZ,为了减少后续的小波分解的计算量,需进行降采样(如将采样至100HZ),获得第二脑电序列信号。
在本发明实施例中,还可对所述第二脑电序列信号进行高低通滤波。其中,脑电波的目标频段在0.1-60Hz之间,因此可采用FIR低通滤波器滤除所述第二脑电序列信号中的高于100Hz的脑电信息,抑制高频噪声;同时采用FIR高通滤波器滤除所述第二脑电序列信号中低于0.05Hz脑电息,以抑制0.05Hz以下的基线漂移。
S102,对所述第二脑电序列信号的各个时刻的脑电信号进行加权移动平均计算,得到第三脑电序列信号。
在本发明实施例中,所述第二脑电序列信号一般包括低频直流信息,这些低频直流信息可能与各个脑电波的频段重叠,因此为了保证脑电波提取的准确度,需进行去除这些低频直流信息。
在本发明实施例中,可采用加权移动平均算法来去除所述低频直流信息。其中,使用移动平均法进行预测能平滑掉需求的突然波动对预测结果的影响,因此能起到去除所述低频直流信息的作用。
具体地,在一个实施例中,步骤S102包括:
S1021,基于当前的第j个时刻的脑电信号,获取所述第二脑电序列信号中的位于第(j-(N-1)/2)个时刻至第(j+(N-1)/2)个时刻之间的N个脑电信号的能量;其中,N为预设的影响数,且N为奇数,j为大于(N+1)/2的整数。
例如,假设当前要预测的脑电信号x(j)的时刻为第10个时刻(即j=10),影响数N为5,则对当前要预测的脑电信号有影响的脑电信号为第8个至第12个时刻的脑电信号,即x(8)~x(12)。此时,先获取这5个时刻的脑电信号的能量。
S1022,根据预设的权值分布函数为获取的N个脑电信号的能量分配权值;其中,N个脑电信号的能量的权值之和为1。
在本发明实施例中,优选地,所述权值分布函数为正态分布函数,如可为:其中,w(i)为第i个时刻的脑电信号的权值值,t(i)为第i个时刻的脑电信号的时间,τ表示需要放大的局部信息量。如图3所示,采用这种权值分布,避免了将第j点附近所有点都看成是一样的比重,而是按照距离(时间差)赋予其一个比重,实现局部信息量的放大,减弱了距离太远的信息对当前点的影响。
需要说明的是,在计算得到各个脑电信号的能量的权值后,还需要进行归一化,保证N个脑电信号的能量的权值之和为1。
S1023,对所述N个脑电信号的能量根据分配的权值进行加权求和,得到新的第j个时刻的脑电信号的能量。
即:
S1024,依次对所述第二脑电序列信号的各个时刻的脑电信号的能量进行加权求和后,根据所有时刻的新的脑电信号的能量,生成第三脑电序列信号。
本发明实施例中,对脑电信号进行去低频直流时,同时考虑了位于其时刻前与时刻后的脑电信号对自身的影响,因此,在去除直流信息的同时,还能起到避免时刻后的脑电信号对当前时刻的信号起不到影响的作用,无法体现信号的变化的问题。
在另一个实施例中,步骤S102具体包括:
S1025,基于当前的第j个时刻的脑电信号,获取所述第二脑电序列信号中的位于第(j-N)个时刻至第(j-1)个时刻之间的N个脑电信号的能量;其中,N为预设的影响数,j为大于N的整数。
例如,假设当前要修正的脑电信号x(j)的时刻为第10个时刻(即j=10),影响数N为5,则对当前要修正的脑电信号有影响的脑电信号为第5个至第9个时刻的脑电信号,即x(5)~x(9)。此时,先获取这5个时刻的脑电信号的能量。
S1026,根据预设的权值分布函数为获取的N个脑电信号的能量分配权值;其中,N个脑电信号的能量的权值之和为1。
S1027,对所述N个脑电信号的能量根据分配的权值进行加权求和,得到新的第j个时刻的脑电信号的能量。
S1028,依次对所述第二脑电序列信号的各个时刻的脑电信号的能量进行加权求和后,根据所有时刻的新的脑电信号的能量,生成第三脑电序列信号。
本发明实施例中,对脑电信号进行去低频直流时,仅考虑了其时刻前的脑电信号对自身的影响,因此,保证信号的平滑。
S103,对所述第三脑电序列信号进行小波变换,并根据每个脑电波的频率范围重构得到各个脑电波。
具体地,可包括如下步骤:
S1031,根据香农-奈奎斯特采样原理及所述第三脑电序列信号的采样频率进行频段分层,计算得到每层的频率范围。
根据香农-奈奎斯特采样原理,设所述第三脑电序列信号的采样频率为fs,目标频段为f1-f2(Hz),运用小波变换分解的层数为N,由奈奎斯特定律可知:
f1=(fs/2)/N1 (1)
f2=(fs/2)/N2 (2)
N>N1(N1>N2) (3)
即需要重构的小波的层数为N2~N1层。
S1032,依据所述频段分层中每层的频率范围及每个脑电波的频率范围,确定与每个脑电波对应的小波分解和重构所需的层数。
在本发明实施例中,所述脑电波包括频率范围为0.5~3Hz的Delta波、频率范围为3~7Hz的Theta波、频率范围为8~13Hz的Alpha波、频率范围为14~17Hz的Beta波、频率范围为34~50Hz的Gamma波。
其中,Delta波:深度睡眠脑波状态:
当人们的大脑频率处于Delta波时,为深度睡眠、无意识状态。人的睡眠品质好坏与Delta波有非常直接的关系。Delta波睡眠是一种很深沉的睡眠状态,如果在辗转难眠时自己召唤出近似Delta波状态,就能很快地摆脱失眠而进入深沉睡眠。
Theta波:深度放松、无压力的潜意识状态
当人们的大脑频率处于Theta波时,人的意识中断,身体深沉放松,对于外界的信息呈现高度的受暗示状态,即被催眠状态。Theta波对于触发深沉记忆、强化长期记忆等帮助极大,所以Theta波被称为"通往记忆与学习的闸门"。
Alpha波:学习与思考的最佳脑波状态
当人们的大脑频率处于Alpha波时,人的意识清醒,但身体却是放松的,它提供意识与潜意识的“桥梁”。在这种状态下,身心能量耗费最少,相对的脑部获得的能量较高,运作就会更加快速、顺畅、敏锐。Alpha波被认为是人们学习与思考的最佳脑波状态。
Beta波:紧张、压力、脑疲劳时的脑波状态
人们清醒时,大部分时间大脑频率处于Beta波状态。随着Beta波的增加,身体逐渐呈紧张状态,因而削减了体内免疫系统能力,此时人的能量消耗加剧,容易疲倦,若不充分休息,容易堆积压力。适当的Beta波对注意力提升以及认知行为的发展有积极作用。
在本发明实施例中,假设所述待处理心电信号已被降采样至100Hz,则fs为100Hz,信号最高频率为50Hz,根据公式(1)、(2)、(3)可知,每一层对应的频段如下:
频段 | 频率范围/Hz | 频段 | 频率范围/Hz |
A1 | 0~25 | D1 | 25~50 |
A2 | 1~12.5 | D2 | 12.5~25 |
A3 | 0~6.25 | D3 | 6.25~12.5 |
A4 | 0~3.125 | D4 | 3.125~6.25 |
A5 | 0~1.625 | D5 | 1.625~3.125 |
A6 | 0~0.8125 | D6 | 0.8125~1.625 |
A7 | 0~0.40625 | D7 | 0.40625~0.8125 |
A8 | 0~0.203125 | D8 | 0.203125~0.40625 |
A9 | 0~0.10156 | D9 | 0.10156~0.203125 |
以Delta波为例,其频段范围为0.5~3Hz。因此,选用第5、6、7层的近似系数(D5/D6/D7)来重构信号。
S1033,根据与每个脑电波对应的所述小波分解所需的层数及预先选择的母小波进行信号分解,得到与每个脑电波对应的按频段划分的多层波形。
在本发明实施例中,经验证,coifN小波和dmey小波的提取效果较佳,且优选地,以coif3小波基作为母小波时,具有最佳的提取效果。因而本发明实施例采用coif3小波基作为母小波进行小波分解。当然,可以理解的是,在本发明的其他实施例中,也可选取其他的母小波,如db小波等,本发明不做具体限定。
需要说明的是,当以coif3小波基作为母小波时,中心频率和带宽是影响coif3小波时频分辨率关键因素。改变中心频率-带宽比就可以改变coif3小波变换的时频分辨率。当中心频率-带宽比达到最优时,coif3小波变换的时频分辨率最高。
具体地,优化过程如下:
首先,coif3小波的母小波表达式如公式(4)所示。其中,fc表示母波ψ(t)的特征频率,也是中心频率,σt为高斯窗的标准差,通常取值为1,σf为带宽,通常σf=1/2π·σt。分析coif3小波的母小波可知,小波波形振荡衰减的快慢由带宽σf决定,波形的振荡频率由中心频率fc决定。根据公式1可以计算coif3小波的频率分辨率(公式5)和时间分辨率(公式6),其中,fs为采样频率,fc为中心频率,σf为带宽,fi为信号分析频率。
接着,利用Shannon熵优化coif3小波变换中心频率-带宽比的核心思想,就是用概率分布序列pi来表示小波系数,然后计算pi的值,表达式如公式7所示。其中,pi是一个概率分布序列,通过小波系数转换得到,具有不确定性。其转换公式如公式8所示,X(fi,t)为小波系数。中心频率-带宽比fc/σf和Shannon小波熵之间的曲线关系,如图4所示。本次实验中,当中心频率-带宽比fc/σf=4.43时,基于Shannon小波熵概率最优理论,可知当Shannon小波熵达到最小值时,coif3小波中心频率-带宽比参数达到最优,对应的基小波就是与特征成分最匹配的小波。
S1034,根据与所述小波重构所需的层数对应的系数及分解得到的与每个脑电波对应所述多层波形,重构得到各个脑电波。
由步骤S1032可知,小波重构所需的层数为第5、6、7层,此时,即可根据与所述小波重构所需的层数对应的系数及分解得到的所述多层波形进行信号重构,得到Delta波。
在本发明实施例中,只需要根据各个脑电波的频率范围选择相应的层数即可用重构得到各个脑电波,本发明在此不做赘述。
S104,基于重构得到的所述各个脑电波进行特征提取,并根据提取得到的特征量进行分类识别,得到脑电放松度。
本发明实施例提供的基于小波变换的脑电放松度识别方法,通过采用加权移动平均法去除所述第二脑电序列信号中的低频直流信息,保证了在去低频直流信息的同时不会对所述第二脑袋序列信号产生任何影响,保证所述第二脑袋序列信号的平滑,同时采用小波变换的方法准确提取各个脑电波,保证了基于这些提取的脑电波计算得到的特征量的准确性,进而保证了最终脑电放松度的准确识别,为准确的生物反馈指导提供了数据基础和依据。
优选地,在步骤S102之前,还包括:
以所述第三脑电序列信号为原始信号,以与所述第三脑电序列信号同步采集得到的伪迹序列信号为参考信号,采用经函数链神经网络优化的自适应滤波器对所述第三脑电序列信号进行滤波,得到去除伪迹序列信号后的第三脑电序列信号。
在本优选实施例中,考虑到第三脑电序列信号中还包含有各种伪迹序列信号,如舌电伪迹,出汗伪迹,眼电伪迹,脉搏伪迹以及肌电伪迹等干扰。其中,以眼电伪迹和肌电伪迹难以去除的问题,这主要是由于其伪迹信号的幅值较高,是脑电信号的几倍甚至几十倍,而且与脑电信号在频域有混叠。
本优选实施例提出一种经函数链神经网络优化的自适应滤波器,滤除放松脑电信号中的各种伪迹信号。
具体地,首先,构造自适应滤波器,其中自适应滤波器的原理框图如图4所示,其由原始信号(即所述第三脑电序列信号)和参考信号(与所述第三脑电序列信号同步采集得到的伪迹序列信号,如舌电伪迹,出汗伪迹,眼电伪迹,脉搏伪迹以及肌电伪迹中的任意一种)两个输入组成。滤波时,参考信号经自适应滤波后,与原始信号进行比较,得到所需脑电序列信号的预估信号(较为纯净的脑电序列信号),其中,滤波器不断地自我重新调整其权值,从而使目标误差达到最小。
其次,将函数链神经网络(Function Link Neural Network,FLNN)应用到自适应滤波器,利用一组正交基函数将原输入矢量进行维数扩展,将线性参数扩展为非线性,来增强自适应滤波器的非线性处理能力。FLNN由函数扩展和单层感知器两部分组成,函数链神经网络的正交基采用切比雪夫正交多项式,如公式9所示。FLNN的基函数T如公式10所示,网络输出如公式11所示,通过FLNN实现了对输入的非线性扩展,更有助于描述脑电信号的非线性特性。
优选地,步骤S104具体为:
S1041,计算各个脑电波的能量和能量分布密度,并基于各个脑电波的能量及每个脑电波的等电位线能量比率,计算得到每个脑电波的等电位线能量基线。
S1042,根据各个脑电波的能量及等电位线能量基线,统计能量小于电位线能量基线的点的个数。
S1043,根据各个脑电波的能量、能量分布密度、能量小于电位线能量基线的点的个数及第一脑电序列信号的等电位线能量,得到所述第一脑电序列信号的特征量。
S1044,基于预先训练好的分类器对所述特征量进行分类,识别得到与所述第一脑电序列信号对应的脑电放松度。
其中,根据步骤S103提取得到的Delta、Theta、Alpha、Beta、Gamma波可得到各个脑电波的幅度A,由于其随时间变化,因此幅度可表示为A(t)。Delta、Theta、Alpha、Beta、Gamma波的幅度分别为A1(t)~A5(t)。
(1)计算能量P,如公式12所示,其中积分范围为-T/2到T/2。那么Delta、Theta、Alpha、Beta、Gamma波的能量分别为P1~P5,此时能量P1~P5就代表Delta、Theta、Alpha、Beta、Gamma波信号。
(2)计算能量分布密度S,如公式13所示,其中,ω=2πf=2π/T,f为频率,T=1/f;那么Delta、Theta、Alpha、Beta、Gamma波的能量分布密度分别为S1~S5。
在本发明实施例中,得到各个脑电波的能量之后,就可以计算得到每个脑电波的等电位线能量基线。
首先,自学习得到等电位线能量比率fdepsRate。
其中,fdepsRate通过大量数据参数自学习得到,且针对Delta、Theta、Alpha、Beta、Gamma波的fdepsRate值都不一样。
其次:计算等电位线能量基线fdeps。
fdeps=max(P)*fdepsRate (14)
再次:统计切片内能量p小于电位线能量基线fdeps的点的个数fdcount。
即:
if p<fdeps
flcount=flcount+1; (15)
最后:计算fdcount与数据片段总点数的比值fdcountP,该比值fdcountP就是所需的特征量。
分别计算得到每个波(Delta、Theta、Alpha、Beta、Gamma波)的等电位线能量特征fdcountP1~fdcountP5,并计算脑电原始信号的等电位线能量特征fdcountEEG。
这里,S1~S5,P1~P5,fdcountP1~fdcountP5及fdcountEEG就是所述第一脑电序列信号的特征量,一共有16个。
在本发明实施例中,将所述特征量输入到预先训练好的分类器(如神经网络模型或支持向量机中),就可以识别出其分类,并根据分类得到脑电放松度。在得到脑电放松度后,就可以根据脑电放松度进行放松治疗,如可以根据脑电放松度进行放松引导内容的选择、标记和播放,能够精准地选取最适合用户的放松引导内容,搭配耳机,播放给用户;同时伴随着基于放松度的放松引导内容播放音量调制,帮助使用者放松身心,缓解焦虑抑郁,陶冶情操、改善个性弱点、消除心理行为障碍、保持心理和躯体健康。
优选地,在步骤S1043之后,步骤S1044之前,还包括:
基于主成分分析法对所述特征量进行特征选择和特征降维。
通过上述步骤我们得到的特征量分别为P1~P5、S1~S5、fdcountEEG、fdcountP1~fdcountP5,共16维特征,特征的维度较大且含有线性相关的输入项,因此考虑进行特征选择和降维。
在本发明实施例中,可采用主成分分析算法实现特征选择和降维。
具体地,利用主成分分析算法进行特征降维的原理如下:
首先,对所述特征量进行数据标准化处理。
具体地,
其中:
其中,Xi'j是标准化后的特征量;Mj、Sj分别表示原始数据某一列的算术平均值和标准(偏)差。
然后,根据数据标准化处理后的所述特征量得到协方差矩阵。
其中,协方差矩阵D=XTX,即:
其中:
接着,计算所述协方差矩阵的特征根和与每个特征根对应的特征向量;其中,所述特征根的数量为p个,且所述的p个特征根按照大小排序。
其中,DP=Pλ(21)
当只考虑第j个特征值时,有DPj=Pjλj,即求解|D-λjI|=0。依次解出各个λ,并使其按大小顺序排列,即λ1≥λ2≥…,≥λp≥0;然后可求出各特征值对应的特征向量P,进而特征方程求解完成。
再接着,获取所述的p个特征根中,贡献率之和大于预定阈值的前m个特征根;其中,每个特征根的贡献率等于所述特征根的值除以全部的p个特征根的值之和。
首先,计算单个主成分的贡献率并进行累计,根据累计贡献率确定主成分的个数m,从而确定所需要选取的主成分。贡献率的计算公式如公式(22)所述。累计贡献率即前m个贡献率的累积和,如公式(23)所示。所述阈值Dmax一般取在85%~95%之间。根据上一步骤中的特征根排序可知,λ1≥λ2≥…,≥λp≥0,从前往后(也是从大到小)依次对特征根进行累加,当累计贡献率大于Dmax时,停止计算,此时累计计算的特征根λ的数目为m,则只需要选取前m个主成分即可。
最后,根据与所述的前m个特征根对应的特征向量及所述输入样本空间,得到主成分得分矩阵。
其中,所述主成分得分矩阵
请参阅图6,本发明实施例还提供一种基于下小波变换的脑电放松度识别装置100,包括:
降采样单元10,用于对接收的第一脑电序列信号进行降采样,获得第二脑电序列信号;
加权移动平均单元20,用于对所述第二脑电序列信号的各个时刻的脑电信号进行加权移动平均计算,得到第三脑电序列信号;
小波变换单元30,用于对所述第三脑电序列信号进行小波变换,并根据每个脑电波的频率范围重构得到各个脑电波;
放松度识别单元40,用于基于重构得到的所述各个脑电波进行特征提取,并根据提取得到的特征量进行分类识别,得到脑电放松度。
优选地,在一个实施例中,所述加权移动平均单元20,具体包括:
第一脑电信号获取模块,用于基于当前的第j个时刻的脑电信号,获取所述第二脑电序列信号中的位于第(j-(N-1)/2)个时刻至第(j+(N-1)/2)个时刻之间的N个脑电信号的能量;其中,N为预设的影响数,且N为奇数,j为大于(N+1)/2的整数;
第一权值分配模块,用于根据预设的权值分布函数为获取的N个脑电信号的能量分配权值;其中,N个脑电信号的能量的权值之和为1;
第一加权求和模块,用于对所述N个脑电信号的能量根据分配的权值进行加权求和,得到新的第j个时刻的脑电信号的能量,其中,依次对所述第二脑电序列信号的各个时刻的脑电信号的能量进行加权求和后,根据所有时刻的新的脑电信号的能量,生成第三脑电序列信号。
优选地,在另一个实施例中,所述加权移动平均单元20,具体包括:
第二脑电信号获取模块,基于当前的第j个时刻的脑电信号,获取所述第二脑电序列信号中的位于第(j-N)个时刻至第(j-1)个时刻之间的N个脑电信号的能量;其中,N为预设的影响数,j为大于N的整数;
第二权值分配模块,根据预设的权值分布函数为获取的N个脑电信号的能量分配权值;其中,N个脑电信号的能量的权值之和为1;
第二加权求和模块,对所述N个脑电信号的能量根据分配的权值进行加权求和,得到新的第j个时刻的脑电信号的能量;其中,依次对所述第二脑电序列信号的各个时刻的脑电信号的能量进行加权求和后,根据的所有时刻的新的脑电信号的能量,生成第三脑电序列信号。
优选地,所述权值分布函数为正态分布函数。
优选地,所述小波变换单元30具体为:
频率分层模块,用于根据香农-奈奎斯特采样原理及所述第三脑电序列信号的采样频率进行频段分层,计算得到每层的频率范围;
层数确定模块,用于依据所述频段分层中每层的频率范围及每个脑电波的频率范围,确定与每个脑电波对应的小波分解和重构所需的层数;
分解模块,用于根据与每个脑电波对应的所述小波分解所需的层数及预先选择的母小波进行信号分解,得到与每个脑电波对应的按频段划分的多层波形;
重构模块,用于根据与所述小波重构所需的层数对应的系数及分解得到的与每个脑电波对应的所述多层波形,重构得到各个脑电波。
优选地,还包括:
伪迹去除单元,用于以所述第三脑电序列信号为原始信号,以与所述第一脑电序列信号同步采集得到的伪迹序列信号为参考信号,采用经函数链神经网络优化的自适应滤波器对所述第三脑电序列信号进行滤波,得到去除伪迹序列信号后的第三脑电序列信号。
优选地,所述放松度识别单元40具体包括:
能量基线计算模块,用于计算各个脑电波的能量和能量分布密度,并基于各个脑电波的能量及每个脑电波的等电位线能量比率,计算得到每个脑电波的等电位线能量基线;
统计模块,用于根据各个脑电波的能量及等电位线能量基线,统计能量小于所述电位线能量基线的点的个数;
特征量获取模块,用于根据各个脑电波的能量、能量分布密度、能量小于电位线能量基线的点的个数及第一脑电序列信号的等电位线能量,得到所述第一脑电序列信号的特征量;
识别分类单元,用于基于预先训练好的分类器对所述特征量进行分类,识别得到与所述第一脑电序列信号对应的脑电放松度。
优选地,还包括:
特征降维模块,用于对所述特征量进行特征选择和特征降维;其中,所述特征选择和特征降维通过主成分分析算法或者通过ANOVA分析和F检验结合算法实现。
优选地,所述脑电波包括频率范围为0.5~3Hz的Delta波、频率范围为3~7Hz的Theta波、频率范围为8~13Hz的Alpha波、频率范围为14~17Hz的Beta波、频率范围为34~50Hz的Gamma波。
本发明实施例提供的基于小波变换的脑电放松度识别装置100,通过采用加权移动平均算法去除所述第二脑电序列信号中的低频直流信息,保证了在去低频直流信息的同时不会对所述第二脑袋序列信号产生任何影响,保证所述第二脑袋序列信号的平滑,同时采用小波变换的方法准确提取各个脑电波,保证了基于这些提取的脑电波计算得到的特征量的准确性,进而保证了最终脑电放松度的准确识别,为准确的生物反馈指导提供了数据基础和依据。
以上所揭露的仅为本发明一种较佳实施例而已,当然不能以此来限定本发明之权利范围,本领域普通技术人员可以理解实现上述实施例的全部或部分流程,并依本发明权利要求所作的等同变化,仍属于发明所涵盖的范围。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random AccessMemory,RAM)等。
Claims (10)
1.一种基于小波变换的脑电放松度识别方法,其特征在于,包括如下步骤:
对接收的第一脑电序列信号进行降采样,获得第二脑电序列信号;
对所述第二脑电序列信号的各个时刻的脑电信号进行加权移动平均计算,得到第三脑电序列信号;
对所述第三脑电序列信号进行小波变换,并根据每个脑电波的频率范围重构得到各个脑电波;
基于重构得到的所述各个脑电波进行特征提取,并根据提取得到的特征量进行分类识别,得到脑电放松度。
2.根据权利要求1所述的基于小波变换的脑电放松度识别方法,其特征在于,所述对所述第二脑电序列信号的各个时刻的脑电信号进行加权移动平均计算,得到第三脑电序列信号,具体包括:
基于当前的第j个时刻的脑电信号,获取所述第二脑电序列信号中的位于第(j-(N-1)/2)个时刻至第(j+(N-1)/2)个时刻之间的N个脑电信号的能量;其中,N为预设的影响数,且N为奇数,j为大于(N+1)/2的整数;
根据预设的权值分布函数为获取的N个脑电信号的能量分配权值;其中,N个脑电信号的能量的权值之和为1;
对所述N个脑电信号的能量根据分配的权值进行加权求和,得到新的第j个时刻的脑电信号的能量;
依次对所述第二脑电序列信号的各个时刻的脑电信号的能量进行加权求和后,根据所有时刻的新的脑电信号的能量,生成第三脑电序列信号。
3.根据权利要求1所述的基于小波变换的脑电放松度识别方法,其特征在于,所述对所述第二脑电序列信号的各个时刻的脑电信号进行加权移动平均计算,得到第三脑电序列信号,具体包括:
基于当前的第j个时刻的脑电信号,获取所述第二脑电序列信号中的位于第(j-N)个时刻至第(j-1)个时刻之间的N个脑电信号的能量;其中,N为预设的影响数,j为大于N的整数;
根据预设的权值分布函数为获取的N个脑电信号的能量分配权值;其中,N个脑电信号的能量的权值之和为1;
对所述N个脑电信号的能量根据分配的权值进行加权求和,得到新的第j个时刻的脑电信号的能量;
依次对所述第二脑电序列信号的各个时刻的脑电信号的能量进行加权求和后,根据的所有时刻的新的脑电信号的能量,生成第三脑电序列信号。
4.根据权利要求2或3所述的基于小波变换的脑电放松度识别方法,其特征在于,所述权值分布函数为正态分布函数。
5.根据权利要求1所述的基于小波变换的脑电放松度识别方法,其特征在于,所述对所述第三脑电序列信号进行小波变换,并根据每个脑电波的频率范围重构得到各个脑电波具体为:
根据香农-奈奎斯特采样原理及所述第三脑电序列信号的采样频率进行频段分层,计算得到每层的频率范围;
依据所述频段分层中每层的频率范围及每个脑电波的频率范围,确定与每个脑电波对应的小波分解和重构所需的层数;
根据与每个脑电波对应的所述小波分解所需的层数及预先选择的母小波进行信号分解,得到与每个脑电波对应的按频段划分的多层波形;
根据与所述小波重构所需的层数对应的系数及分解得到的与每个脑电波对应的所述多层波形,重构得到各个脑电波。
6.根据权利要求5所述的基于小波变换的脑电放松度识别方法,其特征在于,所述母小波为coif3小波,且所述coif3小波的中心频率-带宽比经小波熵自适应优化。
7.根据权利要求1所述的基于小波变换的脑电放松度识别方法,其特征在于,在基于加权移动平均算法对所述第二脑电序列信号的各个时刻的脑电信号进行处理,得到第三脑电序列信号之前,还包括:
以所述第三脑电序列信号为原始信号,以与所述第三脑电序列信号同步采集得到的伪迹序列信号为参考信号,采用经函数链神经网络优化的自适应滤波器对所述第三脑电序列信号进行滤波,得到去除伪迹序列信号后的第三脑电序列信号。
8.根据权利要求1所述的基于小波变换的脑电放松度识别方法,其特征在于,所述基于所述各个脑电波进行特征提取,并根据提取得到的特征量进行分类识别,得到当前的脑电放松度,具体包括:
计算各个脑电波的能量和能量分布密度,并基于各个脑电波的能量及每个脑电波的等电位线能量比率,计算得到每个脑电波的等电位线能量基线;
根据各个脑电波的能量及等电位线能量基线,统计能量小于所述电位线能量基线的点的个数;
根据各个脑电波的能量、能量分布密度、能量小于电位线能量基线的点的个数及第一脑电序列信号的等电位线能量,得到所述第一脑电序列信号的特征量;
基于预先训练好的分类器对所述特征量进行分类,识别得到与所述第一脑电序列信号对应的脑电放松度。
9.根据权利要求8所述的基于小波变换的脑电放松度识别方法,其特征在于,在根据各个脑电波的能量、能量分布密度、能量小于电位线能量基线的点的个数及第一脑电序列信号的等电位线能量,得到所述第一脑电序列信号的特征量之后,在基于预先训练好的分类器对所述特征量进行分类,识别得到与所述第一脑电序列信号对应的脑电放松度之前,还包括:
基于主成分分析法对所述特征量进行特征选择和特征降维。
10.一种基于小波变换的脑电放松度识别装置,其特征在于,包括:
降采样单元,用于对接收的第一脑电序列信号进行降采样,获得第二脑电序列信号;
加权移动平均单元,用于基于加权移动平均算法对所述第二脑电序列信号的各个时刻的脑电信号进行计算,得到第三脑电序列信号;
小波变换单元,用于对所述第三脑电序列信号进行小波变换,并根据每个脑电波的频率范围重构得到各个脑电波;
放松度识别单元,用于基于重构得到的所述各个脑电波进行特征提取,并根据提取得到的特征量进行分类识别,得到脑电放松度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611184948.2A CN106580319A (zh) | 2016-12-20 | 2016-12-20 | 一种基于小波变换的脑电放松度识别方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611184948.2A CN106580319A (zh) | 2016-12-20 | 2016-12-20 | 一种基于小波变换的脑电放松度识别方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106580319A true CN106580319A (zh) | 2017-04-26 |
Family
ID=58599931
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611184948.2A Pending CN106580319A (zh) | 2016-12-20 | 2016-12-20 | 一种基于小波变换的脑电放松度识别方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106580319A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108491699A (zh) * | 2018-03-01 | 2018-09-04 | 广东欧珀移动通信有限公司 | 电子装置、脑电波解锁方法及相关产品 |
CN109771940A (zh) * | 2019-01-30 | 2019-05-21 | 浙江强脑科技有限公司 | 游戏控制方法、装置及计算机可读存储介质 |
CN117752345A (zh) * | 2024-02-22 | 2024-03-26 | 南京信息工程大学 | 基于脑电信号的人体兴奋程度评估方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101690659A (zh) * | 2009-09-29 | 2010-04-07 | 华东理工大学 | 脑电波分析方法 |
CN102521505A (zh) * | 2011-12-08 | 2012-06-27 | 杭州电子科技大学 | 用于控制意图识别的脑电和眼电信号决策融合方法 |
CN102920453A (zh) * | 2012-10-29 | 2013-02-13 | 泰好康电子科技(福建)有限公司 | 一种脑电波信号处理方法及其装置 |
CN103735262A (zh) * | 2013-09-22 | 2014-04-23 | 杭州电子科技大学 | 一种双树复小波与共空间模式结合的脑电特征提取方法 |
US20150208940A1 (en) * | 2014-01-28 | 2015-07-30 | Covidien Lp | Non-stationary feature relationship parameters for awareness monitoring |
US20150257700A1 (en) * | 2010-08-02 | 2015-09-17 | Chi Yung Fu | Apparatus for treating a patient |
CN106037722A (zh) * | 2016-05-13 | 2016-10-26 | 广州视源电子科技股份有限公司 | 催眠状态脑电信号提取方法与系统 |
-
2016
- 2016-12-20 CN CN201611184948.2A patent/CN106580319A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101690659A (zh) * | 2009-09-29 | 2010-04-07 | 华东理工大学 | 脑电波分析方法 |
US20150257700A1 (en) * | 2010-08-02 | 2015-09-17 | Chi Yung Fu | Apparatus for treating a patient |
CN102521505A (zh) * | 2011-12-08 | 2012-06-27 | 杭州电子科技大学 | 用于控制意图识别的脑电和眼电信号决策融合方法 |
CN102920453A (zh) * | 2012-10-29 | 2013-02-13 | 泰好康电子科技(福建)有限公司 | 一种脑电波信号处理方法及其装置 |
CN103735262A (zh) * | 2013-09-22 | 2014-04-23 | 杭州电子科技大学 | 一种双树复小波与共空间模式结合的脑电特征提取方法 |
US20150208940A1 (en) * | 2014-01-28 | 2015-07-30 | Covidien Lp | Non-stationary feature relationship parameters for awareness monitoring |
CN106037722A (zh) * | 2016-05-13 | 2016-10-26 | 广州视源电子科技股份有限公司 | 催眠状态脑电信号提取方法与系统 |
Non-Patent Citations (2)
Title |
---|
XUE Q,HU YH,TOMPKINS WJ: "Neural-network-based adaptive matched filtering for QRS detection", 《EEE TRANSACTIONS ON BIOMEDICAL ENGINEERING》 * |
胡静,王春生,吴敏,杜玉晓: "基于FL-BPNN自适应滤波器的脑电信号伪迹去除算法", 《 第三十二届中国控制会议 2013》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108491699A (zh) * | 2018-03-01 | 2018-09-04 | 广东欧珀移动通信有限公司 | 电子装置、脑电波解锁方法及相关产品 |
CN108491699B (zh) * | 2018-03-01 | 2020-09-22 | Oppo广东移动通信有限公司 | 电子装置、脑电波解锁方法及相关产品 |
CN109771940A (zh) * | 2019-01-30 | 2019-05-21 | 浙江强脑科技有限公司 | 游戏控制方法、装置及计算机可读存储介质 |
CN109771940B (zh) * | 2019-01-30 | 2022-06-07 | 浙江强脑科技有限公司 | 游戏控制方法、装置及计算机可读存储介质 |
CN117752345A (zh) * | 2024-02-22 | 2024-03-26 | 南京信息工程大学 | 基于脑电信号的人体兴奋程度评估方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Maity et al. | Multifractal detrended fluctuation analysis of alpha and theta EEG rhythms with musical stimuli | |
Güler et al. | Adaptive neuro-fuzzy inference system for classification of EEG signals using wavelet coefficients | |
CN107961007A (zh) | 一种结合卷积神经网络和长短时记忆网络的脑电识别方法 | |
CN106963369B (zh) | 一种基于神经网络模型的脑电放松度识别方法及装置 | |
CN112656427A (zh) | 一种基于维度模型的脑电信号情绪识别方法 | |
CN101515200B (zh) | 基于瞬态视觉诱发脑电的目标选择方法 | |
CN105809124A (zh) | 基于DWT和Parametric t-SNE的运动想象脑电信号的特征提取方法 | |
CN107515424A (zh) | 一种基于vmd与小波包的微震信号降噪滤波方法 | |
CN106778594A (zh) | 基于lmd熵特征和lvq神经网络的运动想象脑电信号识别方法 | |
CN104809434A (zh) | 一种基于单通道脑电信号眼电伪迹去除的睡眠分期方法 | |
CN103961091B (zh) | 基于双树复小波样本熵的运动想象脑电信号特征提取方法 | |
CN108143409A (zh) | 睡眠阶段分期方法及装置 | |
CN111184509A (zh) | 一种基于传递熵的情绪诱导脑电信号分类方法 | |
CN106580319A (zh) | 一种基于小波变换的脑电放松度识别方法及装置 | |
Kang et al. | A method of denoising multi-channel EEG signals fast based on PCA and DEBSS algorithm | |
CN106805969B (zh) | 基于卡尔曼滤波和小波变换的脑电放松度识别方法及装置 | |
CN106805968A (zh) | 一种脑电放松度识别方法及装置 | |
CN112861625B (zh) | 一种堆叠去噪自编码器模型确定方法 | |
CN113180659B (zh) | 一种基于三维特征和空洞全卷积网络的脑电情感识别方法 | |
CN114757236B (zh) | 基于tqwt与svmd的脑电信号去噪优化方法及系统 | |
Tafreshi et al. | Epileptic seizure detection using empirical mode decomposition | |
Husain et al. | An artificial neural network model for classification of epileptic seizures using Huang-Hilbert transform | |
Fattah et al. | An approach for classifying alcoholic and non-alcoholic persons based on time domain features extracted from EEG signals | |
CN106923825B (zh) | 基于频域及相空间的脑电放松度识别方法及装置 | |
CN106648087A (zh) | 一种基于意识任务的特征eeg处理方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170426 |
|
RJ01 | Rejection of invention patent application after publication |