CN106923824B - 基于多空间信号特征的脑电放松度识别方法及装置 - Google Patents
基于多空间信号特征的脑电放松度识别方法及装置 Download PDFInfo
- Publication number
- CN106923824B CN106923824B CN201710187187.4A CN201710187187A CN106923824B CN 106923824 B CN106923824 B CN 106923824B CN 201710187187 A CN201710187187 A CN 201710187187A CN 106923824 B CN106923824 B CN 106923824B
- Authority
- CN
- China
- Prior art keywords
- electroencephalogram
- signal
- processed
- characteristic quantity
- amplitude
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 210000004556 brain Anatomy 0.000 claims abstract description 65
- 238000005070 sampling Methods 0.000 claims abstract description 34
- 238000009826 distribution Methods 0.000 claims abstract description 18
- 230000003321 amplification Effects 0.000 claims abstract description 10
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 10
- 238000012706 support-vector machine Methods 0.000 claims description 43
- 238000000605 extraction Methods 0.000 claims description 33
- 238000003062 neural network model Methods 0.000 claims description 28
- 238000004364 calculation method Methods 0.000 claims description 22
- 230000009467 reduction Effects 0.000 claims description 16
- 239000000284 extract Substances 0.000 abstract description 7
- 238000004422 calculation algorithm Methods 0.000 description 22
- 230000006870 function Effects 0.000 description 15
- 238000012549 training Methods 0.000 description 15
- 238000013528 artificial neural network Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 9
- 238000001914 filtration Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000003044 adaptive effect Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 239000011159 matrix material Substances 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000005315 distribution function Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000015654 memory Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 208000013738 Sleep Initiation and Maintenance disease Diseases 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 206010022437 insomnia Diseases 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000012847 principal component analysis method Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 208000019901 Anxiety disease Diseases 0.000 description 1
- 206010019233 Headaches Diseases 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 208000003443 Unconsciousness Diseases 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000000540 analysis of variance Methods 0.000 description 1
- 230000036506 anxiety Effects 0.000 description 1
- 238000013542 behavioral therapy Methods 0.000 description 1
- 230000001149 cognitive effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 231100000869 headache Toxicity 0.000 description 1
- 230000000147 hypnotic effect Effects 0.000 description 1
- 210000000987 immune system Anatomy 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000007787 long-term memory Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000003183 myoelectrical effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 210000004243 sweat Anatomy 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 208000011580 syndromic disease Diseases 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/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/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/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/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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Psychiatry (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychology (AREA)
- Developmental Disabilities (AREA)
- Educational Technology (AREA)
- Hospice & Palliative Care (AREA)
- Child & Adolescent Psychology (AREA)
- Social Psychology (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种基于多空间信号特征的脑电放松度识别方法及装置,所述方法包括:从待处理脑电序列信号中提取信号波;获取信号波的采样点数及各秒中幅度小于幅度阈值且幅度绝对值最大的点,计算其幅度平均值,根据采样点数及幅度平均值计算幅度概率密度,获得在时域空间的特征量;计算信号波的相空间分布密度,获得在相空间的特征量;计算信号波的能量,并根据能量及各个脑电波的频率范围,计算各个信号波的中心频率,获得在频域空间的特征量;根据在时域空间的特征量、相空间的特征量及在频域空间的特征量进行分类识别,得到脑电放松度。本发明可全面提取各个脑电波的不同空间的特征,从而实现准确的脑电放松度识别。
Description
技术领域
本发明涉及放松治疗领域,尤其涉及一种基于多空间信号特征的脑电放松度识别方法及装置。
背景技术
放松训练是行为疗法中使用最广的技术之一,是在心理学实验的基础上建立和发展起来的咨询和治疗方法,其在治疗焦虑抑郁症、神经性头痛、失眠、高血压病,减轻更年期综合征和转变不良行为模式等方面取得了较好的疗效。
现有的放松训练主要有录音指导、口头指导和生物反馈指导。其中,录音指导方法僵化、没有变化,无法根据受训者的状态变化内容;口头指导则要求对口头指导的对象要求很高,且受到时间、场地限制;生物反馈指导以脑电反馈为主,能够结合前两种方式的优点,因而受到广泛关注。
进行生物反馈指导需要识别用户的放松度,而计算放松度首先需要从用户的脑电信号中提取各个频段的脑电波(包括Delta、Theta、Alpha、Beta、Gamma波),再提取各个脑电波的特征,将这些特征输入到分类器进行分类识别。
现有的特征提取方法一般只能从单一角度提取脑电波的特征,评价方式单一,不能保证分类结果的准确性。而且现有的特征提取算法的计算和处理过程复杂,一方面,增加了对硬件的要求,另一方面,由于计算复杂,也无法及时得到分类结果,进而影响了放松治疗的效果。
发明内容
针对上述问题,本发明的目的在于提供一种基于多空间信号特征的脑电放松度识别方法及装置,可全面的提取出各个脑电波的特征。
本发明提供了一种基于多空间信号特征的脑电放松度识别方法,包括如下步骤:
从接收到的待处理脑电序列信号中提取出对应于各个脑电波的信号波;
获取每个信号波的采样点数及所述信号波的各秒中幅度小于幅度阈值且幅度绝对值最大的点,计算这些点的幅度平均值,并根据所述采样点数及所述幅度平均值计算各个信号波的幅度概率密度,获得所述脑电切片信号在时域空间的特征量;
计算对应于各个脑电波的信号波的相空间分布密度,获得所述待处理脑电序列信号在相空间的特征量;
计算与各个脑电波对应的信号波的能量,并根据各个信号波的能量及各个脑电波的频率范围,计算各个信号波的中心频率,获得所述待处理脑电序列信号在频域空间的特征量;
根据所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行分类识别,得到脑电放松度。
优选地,所述获取每个信号波的采样点数及所述信号波内各秒中幅度小于幅度阈值且幅度绝对值最大的点,计算这些点的幅度平均值,并根据所述采样点数及所述幅度平均值计算各个信号波的幅度概率密度,获得所述脑电切片信号在时域空间的特征量具体包括:
统计每个信号波的采样点数n;
查找所述信号波内各秒中幅度小于幅度阈值并且幅度绝对值最大的点,计算这些点的幅度平均值AmpMaxAverrage;
统计幅度值在(-AmpMaxAverrage×k,+AmpMaxAverrage×k)区间的点的数目Ampnum,其中k为经验值;
根据采样点数n及Ampnum计算各个形变的幅度概率密度,获得所述脑电切片信号在时域空间的特征量。
优选地,所述根据所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行分类识别,得到脑电放松度具体包括:
利用至少一个训练好的支持向量机和至少一个神经网络模型对所述特征量进行分类,得到所述待处理脑电序列信号在各个支持向量机及各个神经网络模型下的分类;
将出现次数最多的分类设置为所述待处理脑电序列信号的分类;
根据所述分类与脑电放松度的对应关系,识别得到与所述待处理脑电序列信号对应的脑电放松度。
优选地,所述计算对应于各个脑电波的信号波的相空间分布密度,获得所述待处理脑电序列信号在相空间的特征量具体包括:
根据所述各个脑电波的信号波及待处理脑电序列信号形成对应的包含有信号的二维图表;
用与二维图表等大小的m*m的网格对应覆盖在各个二维图表上,并统计覆盖有信号的网格数;其中,m为大于1的整数;
根据覆盖有信号的网格数与所述网格的全部网格数,计算各个信号波及待处理脑电序列信号的相空间分布密度,获得待处理脑电序列信号在相空间的特征量。
优选地,在所述根据所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行分类识别,得到脑电放松度之前,还包括:
对所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行特征选取及特征降维。
本发明还提供了一种基于多空间信号特征的脑电放松度识别装置,包括:
信号提取单元,用于从接收到的待处理脑电序列信号中提取出对应于各个脑电波的信号波;
时域特征提取单元,用于获取每个信号波的采样点数及所述信号波内各秒中幅度小于幅度阈值且幅度绝对值最大的点,计算这些点的幅度平均值,并根据所述采样点数及所述幅度平均值计算各个信号波的幅度概率密度,获得所述脑电切片信号在时域空间的特征量;
相空间特征提取单元,用于计算对应于各个脑电波的信号波的相空间分布密度,获得所述待处理脑电序列信号在相空间的特征量;
频域特征提取单元,用于计算与各个脑电波对应的信号波的能量,并根据各个信号波的能量及各个脑电波的频率范围,计算各个信号波的中心频率,获得所述待处理脑电序列信号在频域空间的特征量;
脑电放松度识别单元,用于根据所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行分类识别,得到脑电放松度。
优选地,所述时域特征提取单元具体包括:
采样点数统计模块,用于统计每个信号波的采样点数n。
幅度平均值计算模块,用于查找所述信号波内各秒中幅度小于幅度阈值并且幅度绝对值最大的点,计算这些点的幅度平均值AmpMaxAverrage。
点数统计模块,用于统计幅度值在(-AmpMaxAverrage×k,+AmpMaxAverrage×k)区间的点的数目Ampnum,其中k为经验值。
幅度概率密度计算模块,根据采样点数n及Ampnum计算各个形变的幅度概率密度,获得所述脑电切片信号在时域空间的特征量。
优选地,所述脑电放松度识别单元具体包括:
分类模块,用于利用至少一个训练好的支持向量机和至少一个神经网络模型对所述特征量进行分类,得到所述待处理脑电序列信号在各个支持向量机及各个神经网络模型下的分类。
分类统计模块,用于将出现次数最多的分类设置为所述待处理脑电序列信号的分类;
放松度识别模块,用于根据所述分类与脑电放松度的对应关系,识别得到与所述待处理脑电序列信号对应的脑电放松度。
优选地,所述相空间特征提取单元具体包括:
二维图表生成模块,用于根据所述各个脑电波的信号波及待处理脑电序列信号形成对应的包含有信号的二维图表。
网格数统计模块,用与二维图表等大小的m*m的网格对应覆盖在各个二维图表上,并统计覆盖有信号的网格数;其中,m为大于1的整数。
特征量计算模块,用于根据覆盖有信号的网格数与所述网格的全部网格数,计算各个信号波及待处理脑电序列信号的相空间分布密度,获得待处理脑电序列信号在相空间的特征量。
优选地,还包括:
特征选取降维单元,用于对所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行特征选取及特征降维。
本发明提供的基于多空间信号特征的脑电放松度识别方法及装置,同时从时域空间、相空间、频域空间三个空间提取待处理脑电序列序号的特征量并基于提取得到的特征量进行分类识别得到最终的脑电放松度。相比于单一角度的特征提取,对信号的评价方式更多样化,能更全面的体现信号的特性,避免了单一角度的特征提取容易导致的特征过于片面的问题而影响最终识别精度的问题。本发明可以大大提高识别分类的精度和准确度,为放松治疗提供准确的依据。
附图说明
为了更清楚地说明本发明的技术方案,下面将对实施方式中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的基于多空间信号特征的脑电放松度识别方法的流程示意图。
图2是通过切片得到待处理脑电序列信号的示意图。
图3是本发明实施例提供的对对原始脑电序列信号进行加权移动平均计算的原理图。
图4是自适应滤波器的工作原理图。
图5是SVM的最优超平面分类的示意图。
图6是SVM高维映射的示意图。
图7是本发明实施例提供的基于多空间信号特征的脑电放松度识别装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1,本发明实施例提供了一种基于多空间信号特征的脑电放松度识别方法,其可包括如下步骤:
S101,从接收到的待处理脑电序列信号中提取出对应于各个脑电波的信号波。
在本发明实施例中,所述的各个脑电波可包括Delta波、Theta波、Alpha波、Beta波、Gamma波。其中,一般地,Delta波的频率范围为0.5~3Hz,Theta波的频率范围为3~7Hz、Alpha波的频率范围为8~13Hz、Beta波的频率范围为14~17Hz、Gamma波的频率范围为34~50Hz。
其中,Delta波:深度睡眠脑波状态。
当人的大脑频率处于Delta波时,为深度睡眠、无意识状态。人的睡眠品质好坏与Delta波有非常直接的关系。Delta波睡眠是一种很深沉的睡眠状态,如果在辗转难眠时自己召唤出近似Delta波状态,就能很快地摆脱失眠而进入深沉睡眠。
Theta波:深度放松、无压力的潜意识状态。
当人的大脑频率处于Theta波时,人的意识中断,身体深沉放松,对于外界的信息呈现高度的受暗示状态,即被催眠状态。Theta波对于触发深沉记忆、强化长期记忆等帮助极大,所以Theta波被称为"通往记忆与学习的闸门"。
Alpha波:学习与思考的最佳脑波状态。
当人的大脑频率处于Alpha波时,人的意识清醒,但身体却是放松的,它提供意识与潜意识的“桥梁”。在这种状态下,身心能量耗费最少,相对的脑部获得的能量较高,运作就会更加快速、顺畅、敏锐。Alpha波被认为是人们学习与思考的最佳脑波状态。
Beta波:紧张、压力、脑疲劳时的脑波状态。
人们清醒时,大部分时间大脑频率处于Beta波状态。随着Beta波的增加,身体逐渐呈紧张状态,因而削减了体内免疫系统能力,此时人的能量消耗加剧,容易疲倦,若不充分休息,容易堆积压力。适当的Beta波对注意力提升以及认知行为的发展有积极作用。
在本发明实施例中,在获得所述待处理脑电序列信号后,可根据每个脑电波的频率范围通过滤波(如卡尔曼滤波)、小波变换或者自回归模型等提取算法从所述待处理脑电序列信号中提取出对应于各个脑电波的信号波。其中,提取时可仅用一个算法提取得到对应于各个脑电波的信号波,也可同时通过多个算法进行提取,再对不同算法提取得到的结果进行加权求和,得到最终的信号波。利用多个提取算法进行信号波的提取,可避免单一算法提取出现的误差偏大或稳定性不高的问题。
S102,获取每个信号波的采样点数及所述信号波的各秒中幅度小于幅度阈值且幅度绝对值最大的点,计算这些点的幅度平均值,并根据所述采样点数及所述幅度平均值计算各个信号波的幅度概率密度,获得所述脑电切片信号在时域空间的特征量。
在本发明实施例中,可通过幅度概率密度算法提取各个信号波的幅度概率密度来获取各个信号波的时域特征。
首先,统计每个信号波的采样点数n。
然后,查找所述信号波内各秒中幅度小于幅度阈值并且幅度绝对值最大的点,计算这些点的幅度平均值AmpMaxAverrage;其中,幅度阈值通过全部切片自学习得到,其为经验参数。
接着,统计幅度值在(-AmpMaxAverrage×k,+AmpMaxAverrage×k)区间的点的数目Ampnum,其中k为经验值,例如可选取0.4。
最后,根据采样点数n及Ampnum计算各个形变的幅度概率密度,获得所述脑电切片信号在时域空间的特征量AmpPD,即
在本发明实施例中,分别计算得到每个信号波(Delta、Theta、Alpha、Beta、Gamma波)的幅度概率密度AmpPDP1~AmpPDP5;按照同样的方法可以计算脑电切片信号的幅度概率密度AmpPDEEG,即得到所述脑电切片信号在时域空间的特征量。
S103,计算对应于各个脑电波的信号波的相空间分布密度,获得所述待处理脑电序列信号在相空间的特征量。
具体地:
S1031,根据所述各个脑电波的信号波及待处理脑电序列信号形成对应的包含有信号的二维图表。
其中,二维图表的横坐标为时间,纵坐标为信号的幅值。
S1032,用与二维图表等大小的m*m的网格对应覆盖在各个二维图表上,并统计覆盖有信号的网格数。
其中,m为大于1的整数,且m的取值由信号的长度决定。
S1033,根据覆盖有信号的网格数与所述网格的全部网格数,计算各个信号波及待处理脑电序列信号的相空间分布密度,获得待处理脑电序列信号在相空间的特征量。
即,相空间分布密度=md/m2。其中,md为覆盖有信号的网格数。
在本发明实施例中,在得到各个信号波及所述待处理脑电序列信号的相空间分布密度后,就得到了所述待处理脑电序列信号在相空间的特征量。
S104,计算与各个脑电波对应的信号波的能量,并根据各个信号波的能量及各个脑电波的频率范围,计算各个信号波的中心频率,获得所述待处理脑电序列信号在频域空间的特征量。
在本发明实施例中,可通过中心频率算法提取各个信号波的中心频率,得到所述待处理脑电序列信号在频域空间的特征量。
具体地:
S1041,计算与各个脑电波对应的信号波的能量。
S1042,根据各个信号波的能量及各个脑电波的频率范围,计算各个信号波的中心频率,获得所述待处理脑电序列信号在频域空间的特征量。
如公式1所示,中心频率FC:
其中,fH为对应的脑电波的上限频率,fL为对应的脑电波的下限频率,例如Delta波的上限频率为3Hz,下限频率为0.5Hz。
在本发明实施例中,依次计算各个信号波的中心频率,就得到所述待处理脑电序列信号在频域空间的特征量。
S105,根据所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行分类识别,得到脑电放松度。
在本发明实施例中,在得到所述待处理脑电序列信号在时域空间、相空间及频域空间的特征量之后,将其输入到预先设置好的分类器中,就可以得到当前的脑电放松度。
本发明实施例提供的基于多空间信号特征的脑电放松度识别方法,同时从时域空间、相空间及频域空间三个空间提取各个脑电波的信号波的特征,得到待处理脑电序列序号的特征量,并基于提取得到的特征量进行分类识别得到最终的脑电放松度。相比于单一角度的特征提取,对信号的评价方式更多样化,能更全面的体现信号的特性,避免了单一角度的特征提取容易导致的特征过于片面的问题而影响最终识别精度的问题。本发明可以大大提高识别分类的精度和准确度,为放松治疗提供准确的依据。此外,本发明实施例提供的特征提取方法计算简单、快速,对硬件要求低,因而能及时输出分类结果,便于实时的放松治疗。
优选地,在步骤S104之后,步骤S105之前,还可以包括:
S106,基于主成分分析法对所述待处理脑电序列信号在时域空间的特征量、在相空间的特征量及在频域空间的特征量进行降维处理,获得降维后的特征量。
具体地:
S1061,将所述待处理脑电序列信号在时域空间的特征量、在相空间的特征量及在频域空间的特征量设置为输入样本空间中的特征量,并对所述输入样本空间进行数据标准化处理。
具体地,将待处理脑电序列信号在时域空间的特征量、在相空间的特征量及在频域空间的特征量设置为输入样本空间X中的元素。对该样本空间X进行数据标准化处理具体为:
其中:
其中,Xi'j是标准化后的新数据;Mj、Sj分别表示原始数据某一列的算术平均值和标准(偏)差。
S1062,根据数据标准化处理后的所述输入样本空间得到协方差矩阵。
其中,协方差矩阵D=XTX,即:
其中:
S1063,计算所述协方差矩阵的特征根和与每个特征根对应的特征向量;其中,所述特征根的数量为p个,且所述的p个特征根按照大小排序。
其中,DP=Pλ (7)
当只考虑第j个特征值时,有DPj=Pjλj,即求解|D-λjI|=0。依次解出各个λ,并使其按大小顺序排列,即λ1≥λ2≥…,≥λp≥0;然后可求出各特征值对应的特征向量P,进而特征方程求解完成。
S1064,获取所述的p个特征根中,贡献率之和大于预定阈值的前m个特征根。
其中,每个特征根的贡献率等于所述特征根的值除以全部的p个特征根的值之和。
首先,计算单个主成分的贡献率并进行累计,根据累计贡献率确定主成分的个数m,从而确定所需要选取的主成分。贡献率的计算公式如公式8所述。累计贡献率即前m个贡献率的累积和,如公式9所示。所述阈值Dmax一般取在85%~95%之间。根据上一步骤中的特征根排序可知,λ1≥λ2≥…,≥λp≥0,从前往后(也是从大到小)依次对特征根进行累加,当累计贡献率大于Dmax时,停止计算,此时累计计算的特征根λ的数目为m,则只需要选取前m个主成分即可。
S1065,根据与所述的前m个特征根对应的特征向量及所述输入样本空间,得到主成分得分矩阵。
其中,所述主成分得分矩阵中的特征量为所述降维后的特征量。
其中,所述主成分得分矩阵
其中,该主成分得分矩阵T中的各个元素即为经过降维后的特征量。
需要说明的是,在本发明实施例中,还可计算主成分的载荷,其中,所述主成分载荷主要反映主成分得分与原变量xj的关联程度,计算公式为:得到各主成分的载荷以后,就可以知道选取的每一个主成分分别对应的原始特征,如有需要,可以根据原始特征的量纲转换回去。
在本发明实施例中,在采用主成分分析法筛选出所获得的待处理脑电序列信号在时域空间、相空间及频域空间的特征量中较为重要的特征量后,即可得到降维后的待处理脑电序列信号的特征量。通过对待处理脑电序列信号在时域空间、相空间和频域空间的特征量进行降维处理,使得后续进行脑电放松度识别时的数据计算量减小,从而提高脑电放松度识别的速度,实现脑电放松度的实时识别。
优选地,在步骤S104之后,步骤S105之前,还可以包括:
S107,结合方差分析和F检验计算每个特征量的显著性指标,并选取显著性指标大于预设阈值的特征量。
具体地:
首先,根据所述特征量、由标准设备同步采集得到的脑电放松度及预先拟合的与每个特征量对应的线性拟合曲线进行方差分析,计算得到每个特征量的拟合值与期望的平方和及原值与拟合值的平方和。
在本发明实施例中,在准备阶段,在采集脑电切片信号的同时,还会利用标准设备(如神念设备)同步采集脑电放松度作为标准的脑电放松度。这样当采集一系列脑电切片信号(A1,A2,A3…An)后,就可以同时得到一系列标准脑电放松度(Y1,Y2….Yn)。其中,在采集脑电切片信号时,可从所述脑电切片信号提取多个特征量,例如,对于脑电切片信号A1,其可提取特征量(X11,X11,X11….),对于脑电切片信号A2,其可提取特征量(X21,X21,X21….)。这样,对于同一种特征量(例如,这里设X11,X21,…Xn1是同一种特征量),可拟合出其与脑电放松度的函数关系,即对坐标点(X11,Y1)、(X21,Y2)、…(Xn1,Yn)进行拟合,这里假设放松度与特征量的关系均是线性关系的,因此每种特征量对应的线性拟合曲线可表示为Yi=Xβi+εi。
在选取特征量阶段,在计算得到特征量后,将当前的所述特征量及由标准设备同步采集得到的脑电放松度代入到与每个特征量对应的线性拟合曲线中进行方差分析,即可分别计算得到对应的拟合值与期望的平方和(SSR),原值与拟合值的平方和(SSE)。
然后,根据拟合值与期望的平方和、原值与拟合值的平方和及自由度,计算得到每个特征量的显著性指标。
其中,显著性指标F=MSR/MSE,且MSR=SSR/1;MSE=SSE/n-2;其中,1和n-2表示自由度。
最后,根据每个特征量的显著性指标,选取出显著性指标大于预设阈值的特征量。
在本发明实施中,设定显著性指标的阈值Fmin。选取所述特征量中,F>Fmin的特征量作为需要输入的特征量。
本优选实施例中,结合方差分析与F检验对所述特征量进行选取,选出较为显著的特征,从而可以提高分类的精确度和速度。
优选地,在步骤S101之前,还包括:
S01,对原始脑电信号进行切片,获得时间长度为30秒的待处理脑电序列信号,并对所述待处理脑电序列信号的各个时刻的脑电信号进行加权移动平均计算,得到去除低频直流信息后的待处理脑电序列信号。
如图2所示,在本发明实施例中,所述待处理脑电序列信号可通过对原始脑电信号进行切片获得。其中,所述原始脑电信号可通过脑电电极采集获得。一般地,脑电电极采集的原始脑电信号的持续时间比较长(如数小时甚至更长),因此进行需要对原始脑电信号进行切片,例如,每个切片的片段为30秒,即每段所述待处理脑电序列信号的长度为30秒。
在本优选实施例中,为了保证提取和滤波的效率和准确性,还可对脑电序列信号进行相应的预处理,例如,通过预处理去除待处理脑电序列信号中的低频直流信息,以避免这些低频直流信息对脑电波提取的干扰。
在本优选实施例中,为了去除原始脑电序列信号中的低频直流信息,可基于加权移动平均算法对降采样后的原始脑电序列信号的各个时刻的脑电信号进行计算,得到所述待处理脑电序列信号。具体地:
首先,基于当前的第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个时刻的脑电信号的能量。
接着,根据预设的权值分布函数为获取的N个脑电信号的能量分配权值;其中,N个脑电信号的能量的权值之和为1。
在本优选实施例中,所述权值分布函数为正态分布函数,如可为:其中,w(i)为第i个时刻的脑电信号的权值值,t(i)为第i个时刻的脑电信号的时间,τ表示需要放大的局部信息量。如图3所示,采用这种权值分布,避免了将第j点附近所有点都看成是一样的比重,而是按照距离(时间差)赋予其一个比重,实现局部信息量的放大,减弱了距离太远的信息对当前点的影响。
需要说明的是,在计算得到各个脑电信号的能量的权值后,还需要进行归一化,保证N个脑电信号的能量的权值之和为1。
然后,对所述N个脑电信号的能量根据分配的权值进行加权求和,得到新的第j个时刻的脑电信号的能量。
即:
最后,依次对所述原始脑电序列信号的各个时刻的脑电信号的能量进行加权求和后,根据所有时刻的新的脑电信号的能量,生成待处理脑电序列信号。
在本发明实施例中,一般地,需要将所述待处理脑电序列信号进行再次切片,如切成时间长度为6秒的片段。
本优选实施例中,对脑电信号进行去低频直流信息,避免了这些低频直流信息与脑电波的频率发生重叠而影响提取的效果,而且对待处理脑电序列信号的各个时刻的脑电信号进行加权移动平均计算,可以避免去除低频直流信息时,造成信号失真,保证信号的真实度。
优选地,在步骤S101之前,还可以包括:
S02,以待处理脑电序列信号为原始信号,以与所述待处理脑电序列信号同步采集得到的伪迹序列信号为参考信号,采用经函数链神经网络优化的自适应滤波器对所述原始脑电序列信号进行滤波,得到去除伪迹序列信号后的待处理脑电序列信号。
需要说明的是,在S101之前,可以单独采用步骤S01对待处理脑电序列信号进行处理,也可以单独采用步骤S02对待处理脑电序列信号进行处理,还可以同时采用步骤S01和步骤S02对待处理脑电序列信号进行处理。当同时采用步骤S01和步骤S02对待处理脑电序列信号进行处理时,步骤S01与步骤S02无特定的先后顺序。
在本优选实施例中,考虑到待处理脑电序列信号中还包含有各种伪迹序列信号,如舌电伪迹,出汗伪迹,眼电伪迹,脉搏伪迹以及肌电伪迹等干扰。其中,以眼电伪迹和肌电伪迹难以去除的问题,这主要是由于其伪迹信号的幅值较高,是脑电信号的几倍甚至几十倍,而且与脑电信号在频域有混叠。
本优选实施例提出一种经函数链神经网络优化的自适应滤波器,滤除待处理脑电信号中的各种伪迹信号。
具体地,首先,构造自适应滤波器,其中自适应滤波器的原理框图如图4所示,其由原始信号(即所述待处理脑电序列信号)和参考信号(与所述待处理脑电序列信号同步采集得到的伪迹序列信号,如舌电伪迹,出汗伪迹,眼电伪迹,脉搏伪迹以及肌电伪迹中的任意一种)两个输入组成。滤波时,参考信号经自适应滤波后,与原始信号进行比较,得到所需脑电序列信号的预估信号(较为纯净的脑电序列信号),其中,滤波器不断地自我重新调整其权值,从而使目标误差达到最小。
其次,将函数链神经网络(FLNN)应用到自适应滤波器,利用一组正交基函数将原输入矢量进行维数扩展,将线性参数扩展为非线性,来增强自适应滤波器的非线性处理能力。FLNN由函数扩展和单层感知器两部分组成,函数链神经网络的正交基采用切比雪夫正交多项式,如公式11所示。FLNN的基函数T如公式12所示,网络输出如公式13所示,通过FLNN实现了对输入的非线性扩展,更有助于描述脑电信号的非线性特性。
优选地,所述步骤S105具体包括:
S1051,利用至少一个训练好的支持向量机和至少一个神经网络模型对所述特征量进行分类,得到所述待处理脑电序列信号在各个支持向量机及各个神经网络模型下的分类;
S1052,将出现次数最多的分类设置为所述待处理脑电序列信号的分类;
S1053,根据所述分类与脑电放松度的对应关系,识别得到当前的脑电放松度。
在本发明实施例中,在得到待处理脑电序列信号的特征量后,将其输入到训练好的至少一个支持向量机(Support Vector Machine,SVM)和至少一个神经网络模型内,对所述特征量进行分类,识别得到与所述待处理脑电序列信号对应的脑电放松度。
具体地,支持向量机的基本思想是在样本空间或特征空间中构造出最优超平面,使得超平面与不同类样本集之间的距离最大,从而达到最大的泛化能力,如图5所示。
下面将介绍SVM的原理。
首先,针对给定二分类样本对{(xi,yi),xi∈RN,yi=±1}(以此类推五分类样本对为{(xi,yi),xi∈RN,yi=1,2,3,4,5}),xi为训练样本,x为待判决样本。训练样本集为线性不可分时,需引入非负松驰变量αi,i=1,2,......,l;分类超平面的优化问题转化为公式14所示。其中,2/||w||表示分类间隔,使分类间隔最大等价于使||w||2最小。使||w||2最小的分类就成为最优分类面。C为误差惩罚参数,是SVM中最重要的可调参数之一。
其次,选取径向基(Radial Basis Function RBF)核函数,如公式15所示。其中γ为RBF核函数的宽度,是SVM中另一个重要的可调参数。
Kx,xi=exp(-γ*||x-xi||2) (15)
最后,应用核函数技术,将输入空间中的非线性问题,通过函数映射到高维特征空间中,在高维空间中构造线性判别函数,求解最优超平面,使得超平面与不同类样本集之间的距离最大,从而达到最大的泛化能力,如图6所示。
在本发明实施例中,构造好SVM后,就可以进行训练了,具体地,将提取得到的特征量作为训练SVM的输入样本X,将神念设备同步采集得到的“放松度”作为金标准,也就是SVM的输出Y。(X,Y)共同组成SVM的训练样本对,进行SVM训练。
在训练好SVM后,就可以利用该SVM进行分类,从而实现放松度的分类识别。
需要说明的是,SVM的分类性能受到诸多因素影响,其中误差惩罚参数C和RBF核函数的宽度γ两个因素最为关键。C为误差惩罚参数,是SVM中最重要的可调参数之一,表示对错分样本比例和算法复杂度折衷,即在确定的特征子空间中调节学习机器置信范围和经验风险比例,使学习机器的推广能力最好。核函数及参数的选择也直接影响到SVM分类好坏。
在本发明实施例中,不同的支持向量机的误差惩罚参数和核函数的宽度参数可由不同的参数寻优算法进行优化得到,如此,每个支持向量机都将识别得到一个分类的结果。
其中,优选地,所述参数寻优算法包括:结合交叉验证法与网格搜索算法、结合留一法与遗传算法、结合交叉验证法与遗传算法、结合交叉验证法与粒子群算法。
在本优选实施例中,所述的至少一个神经网络模型可采用Levenberg-Marquart算法、标准BP算法、增加动量项的BP算法、改进的GA-BP算法进行训练得到,其中,每个算法训练得到一个神经网络模型,即对应一个分类器。
训练时,将提取得到的各个特征量作为训练神经网络的输入样本X,将神念设备(或其他标准设备)同步采集得到的“放松度”作为金标准,也就是神经网络的输出Y。(X,Y)共同组成神经网络的训练样本对,再根据训练算法进行神经网络训练。
在本优选实施例中,利用训练得到的神经网络模型,就可以作为分类器进行脑电放松度识别,识别得到脑电放松度的等级,如可以分别为弱、较弱、中、较强、强五级(即1~5级)。当前也可以分为更其他等级数,如4级,7级等,本发明不做具体限定。
在本优选实施例中个,在得到SVM和神经网络模型之后,将各个空间的特征量输入到每个SVM和神经网络模型中,获得这些特征量在每个SVM和神经网络模型中的分类,并统计每个分类出现的次数,再将出现次数最多的分类作为最终的分类,并根据与分类对应的脑电放松度的等级,识别得到脑电放松度。
本优选实施例中,通过支持向量机与神经网络模型的组合进行特征量的识别分类,并将出现最多的分类作为最终的分类来识别脑电放松度。由于同时兼顾了支持向量机与神经网络模型两种分类器,因而这种分类方式同时具有支持向量机与神经网络模型的优点,对于某些情况或特征,单独使用支持向量机或者神经网络模型可能分类效果不佳或不准确,而同时采样这两种分类器来分类,可最大程度的避免一种分类器识别效果不佳对最终分类结果的影响,因而分类效果更稳定,不易出现偏差。
请参阅图7,本发明还提供一种基于多空间信号特征的脑电放松度识别装置100,包括:
信号提取单元10,用于从接收到的待处理脑电序列信号中提取出对应于各个脑电波的信号波;
时域特征提取单元20,用于获取每个信号波的采样点数及所述信号波内各秒中幅度小于幅度阈值且幅度绝对值最大的点,计算这些点的幅度平均值,并根据所述采样点数及所述幅度平均值计算各个信号波的幅度概率密度,获得所述脑电切片信号在时域空间的特征量;
相空间特征提取单元30,用于计算对应于各个脑电波的信号波的相空间分布密度,获得所述待处理脑电序列信号在相空间的特征量;
频域特征提取单元40,用于计算与各个脑电波对应的信号波的能量,并根据各个信号波的能量及各个脑电波的频率范围,计算各个信号波的中心频率,获得所述待处理脑电序列信号在频域空间的特征量;
脑电放松度识别单元50,用于根据所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行分类识别,得到脑电放松度。
优选地,还包括:
加权移动平均单元,用于对原始脑电信号进行切片,获得时间长度为30秒的待处理脑电序列信号,并对所述待处理脑电序列信号的各个时刻的脑电信号进行加权移动平均计算,得到去除低频直流信息后的待处理脑电序列信号;或/及
自适应滤波单元,用于以待处理脑电序列信号为原始信号,以与所述待处理脑电序列信号同步采集得到的伪迹序列信号为参考信号,采用经函数链神经网络优化的自适应滤波器对所述原始脑电序列信号进行滤波,得到去除伪迹序列信号后的待处理脑电序列信号。
优选地,所述时域特征提取单元20具体包括:
采样点数统计模块,用于统计每个信号波的采样点数n。
幅度平均值计算模块,用于查找所述信号波内各秒中幅度小于幅度阈值并且幅度绝对值最大的点,计算这些点的幅度平均值AmpMaxAverrage。
点数统计模块,用于统计幅度值在(-AmpMaxAverrage×k,+AmpMaxAverrage×k)区间的点的数目Ampnum,其中k为经验值。
幅度概率密度计算模块,根据采样点数n及Ampnum计算各个形变的幅度概率密度,获得所述脑电切片信号在时域空间的特征量。
优选地,所述相空间特征提取单元30具体包括:
二维图表生成模块,用于根据所述各个脑电波的信号波及待处理脑电序列信号形成对应的包含有信号的二维图表。
网格数统计模块,用与二维图表等大小的m*m的网格对应覆盖在各个二维图表上,并统计覆盖有信号的网格数;其中,m为大于1的整数。
特征量计算模块,用于根据覆盖有信号的网格数与所述网格的全部网格数,计算各个信号波及待处理脑电序列信号的相空间分布密度,获得待处理脑电序列信号在相空间的特征量。
优选地,所述基于多空间信号特征的脑电放松度识别装置还包括特征选取降维单元,用于对所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行特征选取及特征降维。
其中,可基于主成分分析法进行特征降维,基于方差分析和F检验进行特征选取。
优选地,所述脑电放松度识别单元50具体包括:
分类模块,用于利用至少一个训练好的支持向量机和至少一个神经网络模型对所述特征量进行分类,得到所述待处理脑电序列信号在各个支持向量机及各个神经网络模型下的分类。
分类统计模块,用于将出现次数最多的分类设置为所述待处理脑电序列信号的分类;
放松度识别模块,用于根据所述分类与脑电放松度的对应关系,识别得到与所述待处理脑电序列信号对应的脑电放松度。
本优选实施例中,通过支持向量机与神经网络模型的组合进行特征量的识别分类,并将出现最多的分类作为最终的分类来识别脑电放松度。由于同时兼顾了支持向量机与神经网络模型两种分类器,因而这种分类方式同时具有支持向量机与神经网络模型的优点,对于某些情况或特征,单独使用支持向量机或者神经网络模型可能分类效果不佳或不准确,而同时采样这两种分类器来分类,可最大程度的避免一种分类器识别效果不佳对最终分类结果的影响,因而分类效果更稳定,不易出现偏差。
以上所揭露的仅为本发明一种较佳实施例而已,当然不能以此来限定本发明之权利范围,本领域普通技术人员可以理解实现上述实施例的全部或部分流程,并依本发明权利要求所作的等同变化,仍属于发明所涵盖的范围。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(RandomAccessMemory,RAM)等。
Claims (8)
1.一种基于多空间信号特征的脑电放松度识别方法,其特征在于,包括如下步骤:
从接收到的待处理脑电序列信号中提取出对应于各个脑电波的信号波;
获取每个信号波的采样点数及所述信号波的各秒中幅度小于幅度阈值且幅度绝对值最大的点,计算这些点的幅度平均值,并根据所述采样点数及所述幅度平均值计算各个信号波的幅度概率密度,获得脑电切片信号在时域空间的特征量;
计算对应于各个脑电波的信号波的相空间分布密度,获得所述待处理脑电序列信号在相空间的特征量;
计算与各个脑电波对应的信号波的能量,并根据各个信号波的能量及各个脑电波的频率范围,计算各个信号波的中心频率,获得所述待处理脑电序列信号在频域空间的特征量;
根据所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行分类识别,得到脑电放松度;
其中,所述计算对应于各个脑电波的信号波的相空间分布密度,获得所述待处理脑电序列信号在相空间的特征量具体包括:
根据所述各个脑电波的信号波及待处理脑电序列信号形成对应的包含有信号的二维图表;
用与二维图表等大小的m*m的网格对应覆盖在各个二维图表上,并统计覆盖有信号的网格数;其中,m为大于1的整数;
根据覆盖有信号的网格数与所述网格的全部网格数,计算各个信号波及待处理脑电序列信号的相空间分布密度,获得待处理脑电序列信号在相空间的特征量。
2.根据权利要求1所述的基于多空间信号特征的脑电放松度识别方法,其特征在于,所述获取每个信号波的采样点数及所述信号波内各秒中幅度小于幅度阈值且幅度绝对值最大的点,计算这些点的幅度平均值,并根据所述采样点数及所述幅度平均值计算各个信号波的幅度概率密度,获得所述脑电切片信号在时域空间的特征量具体包括:
统计每个信号波的采样点数n;
查找所述信号波内各秒中幅度小于幅度阈值并且幅度绝对值最大的点,计算这些点的幅度平均值AmpMaxAverrage;
统计幅度值在(-AmpMaxAverrage×k,+AmpMaxAverrage×k)区间的点的数目Ampnum,其中k为经验值;
根据采样点数n及Ampnum计算各个形变的幅度概率密度,获得所述脑电切片信号在时域空间的特征量。
3.根据权利要求1所述的基于多空间信号特征的脑电放松度识别方法,其特征在于,所述根据所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行分类识别,得到脑电放松度具体包括:
利用至少一个训练好的支持向量机和至少一个神经网络模型对所述特征量进行分类,得到所述待处理脑电序列信号在各个支持向量机及各个神经网络模型下的分类;
将出现次数最多的分类设置为所述待处理脑电序列信号的分类;
根据所述分类与脑电放松度的对应关系,识别得到与所述待处理脑电序列信号对应的脑电放松度。
4.根据权利要求1所述的基于多空间信号特征的脑电放松度识别方法,其特征在于,在所述根据所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行分类识别,得到脑电放松度之前,还包括:
对所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行特征选取及特征降维。
5.一种基于多空间信号特征的脑电放松度识别装置,其特征在于,包括:
信号提取单元,用于从接收到的待处理脑电序列信号中提取出对应于各个脑电波的信号波;
时域特征提取单元,用于获取每个信号波的采样点数及所述信号波内各秒中幅度小于幅度阈值且幅度绝对值最大的点,计算这些点的幅度平均值,并根据所述采样点数及所述幅度平均值计算各个信号波的幅度概率密度,获得脑电切片信号在时域空间的特征量;
相空间特征提取单元,用于计算对应于各个脑电波的信号波的相空间分布密度,获得所述待处理脑电序列信号在相空间的特征量;
频域特征提取单元,用于计算与各个脑电波对应的信号波的能量,并根据各个信号波的能量及各个脑电波的频率范围,计算各个信号波的中心频率,获得所述待处理脑电序列信号在频域空间的特征量;
脑电放松度识别单元,用于根据所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行分类识别,得到脑电放松度;
其中,所述相空间特征提取单元具体包括:
二维图表生成模块,用于根据所述各个脑电波的信号波及待处理脑电序列信号形成对应的包含有信号的二维图表;
网格数统计模块,用与二维图表等大小的m*m的网格对应覆盖在各个二维图表上,并统计覆盖有信号的网格数;其中,m为大于1的整数;
特征量计算模块,用于根据覆盖有信号的网格数与所述网格的全部网格数,计算各个信号波及待处理脑电序列信号的相空间分布密度,获得待处理脑电序列信号在相空间的特征量。
6.根据权利要求5所述的基于多空间信号特征的脑电放松度识别装置,其特征在于,所述时域特征提取单元具体包括:
采样点数统计模块,用于统计每个信号波的采样点数n;
幅度平均值计算模块,用于查找所述信号波内各秒中幅度小于幅度阈值并且幅度绝对值最大的点,计算这些点的幅度平均值AmpMaxAverrage;
点数统计模块,用于统计幅度值在(-AmpMaxAverrage×k,+AmpMaxAverrage×k)区间的点的数目Ampnum,其中k为经验值;
幅度概率密度计算模块,根据采样点数n及Ampnum计算各个形变的幅度概率密度,获得所述脑电切片信号在时域空间的特征量。
7.根据权利要求5所述的基于多空间信号特征的脑电放松度识别装置,其特征在于,所述脑电放松度识别单元具体包括:
分类模块,用于利用至少一个训练好的支持向量机和至少一个神经网络模型对所述特征量进行分类,得到所述待处理脑电序列信号在各个支持向量机及各个神经网络模型下的分类;
分类统计模块,用于将出现次数最多的分类设置为所述待处理脑电序列信号的分类;
放松度识别模块,用于根据所述分类与脑电放松度的对应关系,识别得到与所述待处理脑电序列信号对应的脑电放松度。
8.根据权利要求5所述的基于多空间信号特征的脑电放松度识别装置,其特征在于,还包括:
特征选取降维单元,用于对所述待处理脑电序列信号在时域空间的特征量、相空间的特征量及在频域空间的特征量进行特征选取及特征降维。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710187187.4A CN106923824B (zh) | 2017-03-27 | 2017-03-27 | 基于多空间信号特征的脑电放松度识别方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710187187.4A CN106923824B (zh) | 2017-03-27 | 2017-03-27 | 基于多空间信号特征的脑电放松度识别方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106923824A CN106923824A (zh) | 2017-07-07 |
CN106923824B true CN106923824B (zh) | 2019-12-20 |
Family
ID=59424799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710187187.4A Active CN106923824B (zh) | 2017-03-27 | 2017-03-27 | 基于多空间信号特征的脑电放松度识别方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106923824B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109953757A (zh) * | 2017-12-14 | 2019-07-02 | 中国航天员科研训练中心 | 面向持续追踪控制与射击类任务脑力负荷实时监测方法 |
CN108491792B (zh) * | 2018-03-21 | 2022-07-12 | 安徽大学 | 基于眼电信号的办公场景人机交互行为识别方法 |
CN112617772B (zh) * | 2021-01-05 | 2022-12-27 | 上海工程技术大学 | 一种基于脉搏波信号的驾驶疲劳识别方法及系统 |
CN113143291B (zh) * | 2021-05-11 | 2023-02-03 | 燕山大学 | 一种快速序列视觉呈现下的脑电特征提取方法 |
CN114209323B (zh) * | 2022-01-21 | 2024-05-10 | 中国科学院计算技术研究所 | 一种识别情绪的方法以及基于脑电数据的情绪识别模型 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0468999A4 (en) * | 1989-04-12 | 1993-01-27 | Aspect Medical Systems, Inc. | Cerebral biopotential analysis system and method |
CN102113879A (zh) * | 2009-12-30 | 2011-07-06 | 上海东方脑科学研究所 | 一种脑波实时评价系统及其评价方法 |
CN102609090A (zh) * | 2012-01-16 | 2012-07-25 | 中国人民解放军国防科学技术大学 | 采用脑电时频成分双重定位范式的快速字符输入方法 |
CN103989485A (zh) * | 2014-05-07 | 2014-08-20 | 朱晓斐 | 基于脑电波的人体疲劳度评价方法 |
CN104173045A (zh) * | 2014-08-15 | 2014-12-03 | 浙江大学医学院附属第二医院 | 一种癫痫发作预警系统 |
CN104887224A (zh) * | 2015-05-29 | 2015-09-09 | 北京航空航天大学 | 面向癫痫脑电信号的特征提取与自动识别方法 |
CN106175799A (zh) * | 2015-04-30 | 2016-12-07 | 深圳市前海览岳科技有限公司 | 基于脑电波评估人体情绪和疲劳状态的方法及系统 |
CN106491125A (zh) * | 2016-11-04 | 2017-03-15 | 广州视源电子科技股份有限公司 | 一种脑电状态识别方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160220836A1 (en) * | 2015-01-30 | 2016-08-04 | Board Of Trustees Of The University Of Arkansas | Device and method of phase-locking brain stimulation to electroencephalographic rhythms |
-
2017
- 2017-03-27 CN CN201710187187.4A patent/CN106923824B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0468999A4 (en) * | 1989-04-12 | 1993-01-27 | Aspect Medical Systems, Inc. | Cerebral biopotential analysis system and method |
CN102113879A (zh) * | 2009-12-30 | 2011-07-06 | 上海东方脑科学研究所 | 一种脑波实时评价系统及其评价方法 |
CN102609090A (zh) * | 2012-01-16 | 2012-07-25 | 中国人民解放军国防科学技术大学 | 采用脑电时频成分双重定位范式的快速字符输入方法 |
CN103989485A (zh) * | 2014-05-07 | 2014-08-20 | 朱晓斐 | 基于脑电波的人体疲劳度评价方法 |
CN104173045A (zh) * | 2014-08-15 | 2014-12-03 | 浙江大学医学院附属第二医院 | 一种癫痫发作预警系统 |
CN106175799A (zh) * | 2015-04-30 | 2016-12-07 | 深圳市前海览岳科技有限公司 | 基于脑电波评估人体情绪和疲劳状态的方法及系统 |
CN104887224A (zh) * | 2015-05-29 | 2015-09-09 | 北京航空航天大学 | 面向癫痫脑电信号的特征提取与自动识别方法 |
CN106491125A (zh) * | 2016-11-04 | 2017-03-15 | 广州视源电子科技股份有限公司 | 一种脑电状态识别方法及装置 |
Non-Patent Citations (4)
Title |
---|
Enhanced performance by time-frequency-phase feature for EEG-based BCI systems;Xu B.etal;《Scientific World Journal》;20140617;全文 * |
基于多变量相空间重构的癫痫脑电分析;孙卫国等;《医疗卫生装备》;20060531(第05期);全文 * |
脑机接口中脑电信号特征提取与分类方法的研究;孙晓晔;《中国优秀硕士学位论文全文数据库》;20131128;全文 * |
认知脑电信号检测分析软件的研究;许佳;《中国优秀硕士学位论文全文数据库 信息科技辑》;20091231(第6期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN106923824A (zh) | 2017-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106963369B (zh) | 一种基于神经网络模型的脑电放松度识别方法及装置 | |
CN106923824B (zh) | 基于多空间信号特征的脑电放松度识别方法及装置 | |
Wu et al. | Estimation of valence of emotion using two frontal EEG channels | |
Zhang et al. | PNN for EEG-based Emotion Recognition | |
CN109994203B (zh) | 一种基于eeg信号深度多视角特征学习的癫痫检测方法 | |
CN114052735B (zh) | 一种基于深度领域自适应的脑电情绪识别方法及系统 | |
Kumari et al. | Seizure detection in EEG using time frequency analysis and SVM | |
CN106108894A (zh) | 一种提高情绪识别模型时间鲁棒性的情绪脑电识别方法 | |
Subasi et al. | EEG signal classification using PCA, ICA, LDA and support vector machines | |
CN112656427A (zh) | 一种基于维度模型的脑电信号情绪识别方法 | |
CN106923825B (zh) | 基于频域及相空间的脑电放松度识别方法及装置 | |
Biju et al. | Ictal EEG classification based on amplitude and frequency contours of IMFs | |
CN106805969B (zh) | 基于卡尔曼滤波和小波变换的脑电放松度识别方法及装置 | |
Cheng et al. | Emotion recognition algorithm based on convolution neural network | |
CN106974648B (zh) | 基于时域及频域空间的脑电放松度识别方法及装置 | |
Saddam et al. | Classification of alcoholic EEG using wavelet packet decomposition, principal component analysis, and combination of genetic algorithm and neural network | |
CN114218986A (zh) | 基于eeg脑电信号数据的状态分类方法 | |
CN109009098A (zh) | 一种运动想象状态下的脑电信号特征识别方法 | |
CN106682605B (zh) | 一种识别脑电放松度的方法及系统 | |
CN110946575A (zh) | 一种基于差分熵和的l1/2脑电波情感分类识别方法 | |
Bawane et al. | Classification of mental task based on EEG processing using self organising feature map | |
Theresia et al. | The back-propagation neural network classification of EEG signal using time frequency domain feature extraction | |
MILICEVIC et al. | Optimal set of EEG features in infant sleep stage classification | |
Peng et al. | DW-FBCSP: EEG emotion recognition algorithm based on scale distance weighted optimization | |
Liu et al. | Research on feature fusion for emotion recognition based on discriminative canonical correlation analysis |
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 |