CN109492546B - 一种融合小波包与互信息的生物信号特征提取方法 - Google Patents
一种融合小波包与互信息的生物信号特征提取方法 Download PDFInfo
- Publication number
- CN109492546B CN109492546B CN201811241761.0A CN201811241761A CN109492546B CN 109492546 B CN109492546 B CN 109492546B CN 201811241761 A CN201811241761 A CN 201811241761A CN 109492546 B CN109492546 B CN 109492546B
- Authority
- CN
- China
- Prior art keywords
- signals
- biological
- mutual information
- biological signal
- algorithm
- 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
- 238000000605 extraction Methods 0.000 title claims abstract description 15
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 27
- 238000012549 training Methods 0.000 claims abstract description 19
- 238000012360 testing method Methods 0.000 claims abstract description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 37
- 239000013598 vector Substances 0.000 claims description 17
- 238000005070 sampling Methods 0.000 claims description 16
- 238000000034 method Methods 0.000 claims description 15
- 238000001514 detection method Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 8
- 238000013528 artificial neural network Methods 0.000 claims description 5
- 230000009977 dual effect Effects 0.000 claims description 5
- 230000011218 segmentation Effects 0.000 claims description 4
- 206010041235 Snoring Diseases 0.000 claims description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 3
- 239000008280 blood Substances 0.000 claims description 3
- 210000004369 blood Anatomy 0.000 claims description 3
- 210000004556 brain Anatomy 0.000 claims description 3
- 210000001715 carotid artery Anatomy 0.000 claims description 3
- 230000000295 complement effect Effects 0.000 claims description 3
- 230000005021 gait Effects 0.000 claims description 3
- 210000000056 organ Anatomy 0.000 claims description 3
- 229910052760 oxygen Inorganic materials 0.000 claims description 3
- 239000001301 oxygen Substances 0.000 claims description 3
- 230000010349 pulsation Effects 0.000 claims description 3
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 3
- 210000001525 retina Anatomy 0.000 claims description 3
- 238000007635 classification algorithm Methods 0.000 claims description 2
- 238000003066 decision tree Methods 0.000 claims description 2
- 238000003064 k means clustering Methods 0.000 claims description 2
- 238000007477 logistic regression Methods 0.000 claims description 2
- 238000000638 solvent extraction Methods 0.000 claims description 2
- 210000002784 stomach Anatomy 0.000 claims description 2
- 238000012706 support-vector machine Methods 0.000 claims description 2
- 241000221931 Hypomyces rosellus Species 0.000 claims 1
- 238000001914 filtration Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 210000000887 face Anatomy 0.000 description 2
- 210000000554 iris Anatomy 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006806 disease prevention Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000008451 emotion Effects 0.000 description 1
- 230000001037 epileptic effect Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000002496 gastric effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明提供一种融合小波包与互信息的生物信号特征提取方法，包括以下步骤：S1)、获取生物信号，并去除低频场电位和高频噪声；S2)、进行生物信号的分割；S3)、对生物信号进行小波包分解；S4)、随机选取少数用小波包分解系数表示的生物信号作为训练数据；S5)、运用最大互信息准则对训练生物信号数据进行特征选择；S6)、对测试生物信号进行特征提取，得到测试生物信号的低维时‑频特征；S7)、应用分类器对所述生物信号的低维时‑频特征进行分类，进而得到目标生物信号和多种非目标生物信号的多分类结果，本发明所提取的脑电信号特征经过聚类处理结果显示大大提高了正常与癫痫信号的分类精度。
Description
技术领域
本发明涉及生物信号特征提取技术领域，尤其是一种融合小波包与互信息的生物信号特征提取方法。
背景技术
生物信号在医学、生命科学等领域具有重要的研究价值，目前生物信号处理在医疗保健、疾病预防、治疗以及医疗仪器产业中展现出其重要作用，由于生物信号来自于生命体，受环境、生理、心理的影响，信号具有多样性和复杂性等特点。并随着科学技术的发展，所采集的生物信号在维度、广度和数量上也呈现出爆炸式的增长趋势。
为了很好的利用所采集的生物信号，对其进行特征提取是很有必要的。现有的针对特定的生物信号的进行特征提取的技术有很多，比如：北京理工大学的申请号为201410543404.5、发明名称为“多尺度特征融合的心电图心电信号分类方法”所针对的是心电信号进行特征提取的方法，东南大学的申请号为201310168156.6、发明名称为“基于巴氏距离最优小波包分解的语音情感特征提取方法”所针对的是语音信号进行特征提取的方法等等。这些方法针对的是特定生物信号进行特征提取的方法，具有一定的针对性，同时缺乏普适性。
本发明保留了生物信号经过小波包分解后的原始时-频域特征，相对于利用小波包分解后的特征重构新特征或者计算频带能量的方式具有特征分辨率更强的优势。同时利用互信息的良好表达能力，有效剔除冗余的小波包分解特征，使得该特征提取方法具有良好的普适性，针对生物信号具有良好的特征提取能力。
发明内容
针对现有技术的不足，本发明提供一种融合小波包与互信息的生物信号特征提取方法。
本发明的技术方案为：一种融合小波包与互信息的生物信号特征提取方法，包括以下步骤：
S1)、读取数据库中所有的生物信号，并使用滤波器去除所述生物信号中的低频场电位和高频噪声；
S2)、进行生物信号的分割，首先根据阈值检测得到每个生物信号的中心发放时间，并根据该时间向前采样K个点、向后采样K-1个点作为一个周期为2K个采样点的生物信号，从而得到所有单独的生物信号；
S3)、对生物信号进行小波包分解，得到每个生物信号的小波包分解系数；
S4)、随机选取少数用小波包分解系数表示的生物信号作为训练数据，每个用于训练的生物信号数据包括多维的特征，即小波包分解系数，以及一维的类别标签；
S5)、运用最大互信息准则对训练生物信号数据进行特征选择，得到相应的系数索引；
S6)、对经过步骤S1)至步骤S3)处理的测试生物信号，采用步骤S5)得到的系数索引进行特征提取，得到测试生物信号的低维时-频特征；
S7)、应用分类器对所述生物信号的低维时-频特征进行分类，进而得到目标生物信号和多种非目标生物信号的多分类结果。
进一步的，上述技术方案中，步骤S1)中，所述的生物信号为心电信号、脑波信号、肌电图、胃电信号、眼视网膜信号、眼电信号、非电生理信号、人体光电信号、耳声发射信号、人体生理特征信号中的多种，其中，非电生理信号是指由器官运动产生的心音、脉搏、颈动脉搏动、呼吸、鼾声和柯氏音，压力信号；人体光电信号是指光电脉波、血氧饱和度信号；人体生理特征信号是指指纹、掌形、面部、虹膜和步态等相关信号。
进一步的，上述技术方案中，步骤S1)中，所述的滤波器包括低通滤波器、高通滤波器、带通滤波器和带阻滤波器中的一种或多种。通过滤波器对所有的生物信号进行滤波处理，滤除生物信号中的低频场电位和高频噪声信号，得到处于特定频率范围内的生物信号。
进一步的，上述技术方案中，步骤S2)中，生物信号的分割采用低阈值检测、高阈值检测或双阈值检测中的任意一种。
进一步的，上述技术方案中，所述的双阈值检测，具体如下：
将步骤S1)中的生物信号通过双阈值检测，当信号中的电位值超过一定阈值4μ，并且在该电位附近的K个样本内最低电位值有低于另一阈值-2μ，则记该超过阈值4μ的采样位置为一个生物信号中心发放时间点，其中，μ为上述信号电位的均值；
接着根据该中心发放时间向前采样K个点，向后采样K-1个点作为一个周期为2K个采样点的生物信号。
进一步的，上述技术方案中，步骤S3)中，对生物信号进行小波包分解采用哈尔小波族、多贝西小波族、Symlet小波族、Coiflet小波族、双正交小波族、反向双正交小波族、Fejer-Korovkin filters小波族、离散逼近的Meyer小波族中的任意一种，所采用的小波包分解级数为所有可能值。
进一步的，上述技术方案中，步骤S5)中，运用最大互信息准则对训练生物信号数据进行特征选择，采用基于互信息的Min-Redundancy特征集度量准则和基于近似条件互信息的特征集度量准则进行特征选择，其中，基于近似条件互信息的特征集度量准则对训练生物信号进行特征选择，具体为：
S501)、对训练生物信号的每一维特征的熵和特定二维特征组的联合熵进行估计，令某度量空间中一个有值的随机变量为X，在该空间下X的两个观测值x和x′的有距离函数||x-x′||，熵值H(X)通过使用
得到其对应的无偏估计量，即
其中，M为表示样本数，
是使用k近邻算法对logμ(x)得到的估计值，已知样本xa和它的第k个近邻点xa+k，得到：
其中，ψ(x)为双伽玛函数，即伽玛函数Γ(x)的对数导数ψ(x)＝Γ(x)′/Γ(x)，d表示X的维度，cd表示d维单元球的体积，E(logε)表示对距离ε的对数值求期望；
考虑二维联合随机变量Z＝(X,Y)，并有两个观测值z′和z″，使用最大范式，
||z′-z″||＝max{||x′-x″,||y′-y″||}，
在M个样本中取一个点za，该点与其第k个近邻点间的距离为ε(a)/2，并且该点在X,Y空间的投影表示为εx(a)/2，εy(a)/2，所以ε(a)＝max{εx(a),εy(a)}，记nx(a)为与样本xa的距离严格小于εx(a)/2的样本个数，定义ny(a)为与样本ya的距离严格小于εy(a)/2的样本个数，用dZ＝dX+dY，
za＝(xa,ya)代替d、cd、xa，得到基于k近邻估计的联合熵为：
其中，dZ、dX、dY分别表示样本点za的维度以及该点在X,Y空间投影的维度；
分别表示表示dZ、dX、dY维单元球的体积；
S502)、由估计得到的熵
和联合熵
计算特征与类别标签间的相应互信息和条件互信息，互信息定义了随机变量间相互依赖的衡量方式，给定连续型生物信号特征向量X,Y以及离散的类别标签向量Z，互信息I(Z；X)和条件互信息I(X；Y|Z)的估计方式表示为：
S503)、利用生物信号的类别标签向量、每一维特征向量或特定二维特征向量的互信息或条件互信息迭代选择特征k，并记为S[k]，具体为：
首先计算每一特征{fi|i＝1,2.....,m}，则类别标签向量Z的互信息I(Z；fi)，然后选择使得互信息最大的一个特征fj作为第一个特征S[1]，并将其归入最优特征集S中，同时更新剩余特征集
其中，
表示表示在全部特征集中S的补集，即，
此时，|S|＝1，|NS|＝m-1，然后在此基础上迭代步骤S501)，S502)，并使用条件互信息：
然后进行后续N-1个特征的选择，其中N为最优特征个数。
进一步的，上述技术方案中，步骤S7)中，采用的分类器为基于分类算法相关的分类器和基于聚类算法相关的分类器，其中，所述的基于分类算法是指决策树算法、朴素贝叶斯算法、逻辑回归算法、人工神经网络算法、k-近邻算法、支持向量机算法、关联规则算法、集成学习算法，所述的基于聚类算法是指基于划分的k-均值聚类算法及其变种、基于层次的聚类算法及其变种、基于密度的DBSCAN聚类算法及其变种、基于网络的STING聚类算法及其变种、基于神经网络的自组织神经网络SOM算法及其变种、基于统计学的COBWEB算法及其变种。
本发明的有益效果为：
1、本发明提出了方法适用于普适的各类生物信号的特征提取；
2、本发明保留了生物信号经过小波包分解后的原始时-频域特征，并利用互信息的特性有效剔除冗余特征，使得本发明提取的生物信号特征不受限于所选择的分类系统而具有良好的普适性，针对各类生物信号具有良好的特征提取能力，更为有效的提取生物信号中所蕴含的重要特征；
3、对于脑电信号，本发明所提取的特征经过聚类处理结果显示大大提高了正常与癫痫信号的分类精度。
附图说明
图1为本发明的流程示意图；
图2为本发明基于近似条件互信息的特征集度量准则为例对训练生物信号进行特征选择的流程图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明：
如图1所示，一种融合小波包与互信息的生物信号特征提取方法，包括以下步骤：
S1)、读取数据库中所有的生物信号，并通过带通滤波器对所有的生物信号进行带通滤波处理，滤除生物信号中的低频场电位和高频噪声信号，得到处于特定频率范围内的生物信号，其中，所述的生物信号为心电信号、脑波信号、肌电图、胃电信号、眼视网膜信号、眼电信号、非电生理信号、人体光电信号、耳声发射信号、人体生理特征信号中的多种，其中，非电生理信号是指由器官运动产生的心音、脉搏、颈动脉搏动、呼吸、鼾声和柯氏音，压力信号；人体光电信号是指光电脉波、血氧饱和度信号；人体生理特征信号是指指纹、掌形、面部、虹膜和步态等相关信号；
S2)、进行生物信号的分割，首先根据双阈值检测得到每个生物信号的中心发放时间，并根据该时间向前采样K个点、向后采样K-1个点作为一个周期为2K个采样点的生物信号，从而得到所有单独的生物信号；具体如下：
S201)、将步骤S1)中的生物信号通过双阈值检测，当信号中的电位值超过一定阈值4μ，并且在该电位附近的K个样本内最低电位值有低于另一阈值-2μ，则记该超过阈值4μ的采样位置为一个生物信号中心发放时间点，其中，μ为上述信号电位的均值；
S202)、接着根据该中心发放时间向前采样K个点，向后采样K-1个点作为一个周期为2K个采样点的生物信号；
S3)、采用多贝西函数db4小波母函数对生物信号进行小波包分解，得到每个生物信号的小波包分解系数，选择最大的分解级数对所有的生物信号逐一进行小波包分解，得到每个生物信号的小波包分解系数，并将其作为生物信号的原始特征进行后续的特征选择；
S4)、随机选取少数用小波包分解系数表示的生物信号作为训练数据，每个用于训练的生物信号数据包括多维的特征，即小波包分解系数，以及一维的类别标签；
S5)、运用基于近似条件互信息的特征集度量准则对训练生物信号进行特征选择，得到相应的系数索引，如图2所示，具体为：
S501)、对训练生物信号的每一维特征的熵和特定二维特征组的联合熵进行估计，令某度量空间中一个有值的随机变量为X，在该空间下X的两个观测值x和x′的有距离函数||x-x′||，熵值H(X)通过使用
得到其对应的无偏估计量，即
其中，M为表示样本数，
是使用k近邻算法对logμ(x)得到的估计值，已知样本xa和它的第k个近邻点xa+k，得到：
其中，ψ(x)为双伽玛函数，即伽玛函数Γ(x)的对数导数ψ(x)＝Γ(x)′/Γ(x)，d表示X的维度，cd表示d维单元球的体积，E(logε)表示对距离ε的对数值求期望；
考虑二维联合随机变量Z＝(X,Y)，并有两个观测值z′和z″，使用最大范式，
||z′-z″||＝max{||x′-x″,||y′-y″||}，
在M个样本中取一个点za，该点与其第k个近邻点间的距离为ε(a)/2，并且该点在X,Y空间的投影表示为εx(a)/2，εy(a)/2，所以ε(a)＝max{εx(a),εy(a)}，记nx(a)为与样本xa的距离严格小于εx(a)/2的样本个数，定义ny(a)为与样本ya的距离严格小于εy(a)/2的样本个数，用dZ＝dX+dY，
za＝(xa,ya)代替d、cd、xa，得到基于k近邻估计的联合熵为：
其中，dZ、dX、dY分别表示样本点za的维度以及该点在X,Y空间投影的维度；
分别表示表示dZ、dX、dY维单元球的体积；
S502)、由估计得到的熵
和联合熵
计算计算特征与类别标签间的相应互信息和条件互信息，互信息定义了随机变量间相互依赖的衡量方式，给定连续型生物信号特征向量X,Y以及离散的类别标签向量Z，互信息I(Z；X)和条件互信息I(X；Y|Z)的估计方式表示为：
S503)、利用生物信号的类别标签向量、每一维特征向量或特定二维特征向量的互信息或条件互信息迭代选择特征k，并记为S[k]，具体为：
首先计算每一特征{fi|i＝1,2.....,m}，则类别标签向量Z的互信息I(Z；fi)，然后选择使得互信息最大的一个特征fj作为第一个特征S[1]，并将其归入最优特征集S中，同时更新剩余特征集
其中，
表示表示在全部特征集中S的补集，即，
此时，|S|＝1，|NS|＝m-1，然后在此基础上迭代步骤S501)，S502)，并使用条件互信息：
然后进行后续N-1个特征的选择，其中N为最优特征个数；
S6)、对经过步骤S1)至步骤S3)处理的测试生物信号，采用步骤S5)得到的系数索引进行特征提取，得到测试生物信号的低维时-频特征，
S7)、应用基于模糊C均值聚类的分类器对所述生物信号的低维时-频特征进行分类，将低维时-频特征表示的生物信号C作为模糊C均值聚类的输入进行聚类处理，进而得到目标生物信号和多种非目标生物信号的多分类结果。
上述实施例和说明书中描述的只是说明本发明的原理和最佳实施例，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。
Claims (6)
1.一种融合小波包与互信息的生物信号特征提取方法，其特征在于，包括以下步骤：
S1)、读取数据库中所有的生物信号，并使用滤波器去除所述生物信号中的低频场电位和高频噪声；
S2)、进行生物信号的分割，首先根据阈值检测得到每个生物信号的中心发放时间，并根据该时间向前采样K个点、向后采样K-1个点作为一个周期为2K个采样点的生物信号，从而得到所有单独的生物信号；
S3)、对生物信号进行小波包分解，得到每个生物信号的小波包分解系数；
S4)、随机选取少数用小波包分解系数表示的生物信号作为训练数据，每个用于训练的生物信号数据包括多维的特征，即小波包分解系数，以及一维的类别标签；
S5)、运用最大互信息准则对训练生物信号数据进行特征选择，得到相应的系数索引，采用基于互信息的Min-Redundancy特征集度量准则和基于近似条件互信息的特征集度量准则进行特征选择，其中，基于近似条件互信息的特征集度量准则对训练生物信号进行特征选择，具体为：
S501)、对训练生物信号的每一维特征的熵和特定二维特征组的联合熵进行估计：令某度量空间中一个有值的随机变量X并且存在概率密度函数μ(x)，在该空间下X的任意两个观测值x′和x″有距离函数：||x′-x″||，熵值H(X)通过使用
得到其对应的无偏估计量
即
其中，M为表示样本数，
是使用k近邻算法对logμ(x)得到的估计值，已知样本xa和它的第k个近邻点xa+k，得到：
其中，ψ(x)为双伽玛函数，即伽玛函数Γ(x)的对数导数ψ(x)＝Γ(x)′/Γ(x)，d表示X的维度，cd表示d维单元球的体积，E(logε)表示对距离ε的对数值求期望；
考虑二维联合随机变量Z＝(X,Y)，并有两个观测值z′和z″，使用最大范式，||z′-z″||＝max{||x′-x″,||y′-y″||}，
在M个样本中取一个点za，该点与其第k个近邻点间的距离为ε(a)/2，并且该点在X,Y空间的投影表示为εx(a)/2，εy(a)/2，所以ε(a)＝max{εx(a),εy(a)}，记nx(a)为与样本xa的距离严格小于εx(a)/2的样本个数，定义ny(a)为与样本ya的距离严格小于εy(a)/2的样本个数，用dZ＝dX+dY，
za＝(xa,ya)代替d、cd、xa，得到基于k近邻估计的联合熵为：
其中，dZ、dX、dY分别表示样本点za的维度以及该点在X,Y空间投影的维度；
分别表示dZ、dX、dY维单元球的体积；
S502)、由估计得到的熵
和联合熵
计算特征与类别标签间的相应互信息和条件互信息，互信息定义了随机变量间相互依赖的衡量方式，给定连续型生物信号特征向量X,Y以及离散的类别标签向量Z，互信息I(Z；X)和条件互信息I(X；Y|Z)的估计方式表示为：
S503)、利用生物信号的类别标签向量、每一维特征向量或特定二维特征向量的互信息或条件互信息迭代选择特征k，并记为S[k]，具体为：
首先计算每一特征{fi|i＝1,2.....,m}，则类别标签向量Z的互信息I(Z；fi)，然后选择使得互信息最大的一个特征fj作为第一个特征S[1]，并将其归入最优特征集S中，同时更新剩余特征集
其中，
表示在全部特征集中S的补集，即，
此时，|S|＝1，|NS|＝m-1，然后在此基础上迭代步骤S501)，S502)，并使用条件互信息：
然后进行后续N-1个特征的选择，其中N为最优特征个数；
S6)、对经过步骤S1)至步骤S3)处理的测试生物信号，采用步骤S5)得到的系数索引进行特征提取，得到测试生物信号的低维时-频特征；
S7)、应用分类器对所述生物信号的低维时-频特征进行分类，进而得到目标生物信号和多种非目标生物信号的多分类结果。
2.根据权利要求1所述的一种融合小波包与互信息的生物信号特征提取方法，其特征在于：步骤S1)中，所述的生物信号为心电信号、脑波信号、肌电图、胃电信号、眼视网膜信号、眼电信号、非电生理信号、人体光电信号、耳声发射信号、人体生理特征信号中的多种，其中，非电生理信号是指由器官运动产生的心音、脉搏、颈动脉搏动、呼吸、鼾声和柯氏音，压力信号；人体光电信号是指光电脉波、血氧饱和度信号；人体生理特征信号是指指纹、掌形、面部、虹膜和步态相关的信号。
3.根据权利要求1所述的一种融合小波包与互信息的生物信号特征提取方法，其特征在于：步骤S1)中，所述的滤波器包括低通滤波器、高通滤波器、带通滤波器和带阻滤波器中的一种或多种， 通过滤波器对所有的生物信号进行滤波处理，滤除生物信号中的低频场电位和高频噪声信号，得到处于特定频率范围内的生物信号。
4.根据权利要求1所述的一种融合小波包与互信息的生物信号特征提取方法，其特征在于：步骤S2)中，所述的阈值检测为低阈值检测、高阈值检测或双阈值检测中的任意一种，其中，使用双阈值检测，具体如下：
S201)、将步骤S1)中的生物信号通过双阈值检测，当信号中的电位值超过一定阈值4μ，并且在该电位附近的K个样本内最低电位值有低于另一阈值-2μ，则记该超过阈值4μ的采样位置为一个生物信号中心发放时间点，其中，μ为上述信号电位的均值；
S202)、接着根据该中心发放时间向前采样K个点，向后采样K-1个点作为一个周期为2K个采样点的生物信号。
5.根据权利要求1所述的一种融合小波包与互信息的生物信号特征提取方法，其特征在于：步骤S3)中，对生物信号进行小波包分解采用哈尔小波族、多贝西小波族、Symlet小波族、Coifle小波族、双正交小波族、反向双正交小波族、Fejer-Korovkin filters小波族、离散逼近的小波族中的任意一种，所采用的小波包分解级数为所有可能值。
6.根据权利要求1所述的一种融合小波包与互信息的生物信号特征提取方法，其特征在于：步骤S7)中，采用的分类器为基于分类算法相关的分类器和基于聚类算法相关的分类器，其中，所述的基于分类算法是指决策树算法、朴素贝叶斯算法、逻辑回归算法、人工神经网络算法、k-近邻算法、支持向量机算法、关联规则算法、集成学习算法，所述的基于聚类算法是指基于划分的k-均值聚类算法及其变种、基于层次的聚类算法及其变种、基于密度的DBSCAN聚类算法及其变种、基于网络的STING聚类算法及其变种、基于神经网络的自组织神经网络SOM算法及其变种、基于统计学的COBWEB算法及其变种。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811241761.0A CN109492546B (zh) | 2018-10-24 | 2018-10-24 | 一种融合小波包与互信息的生物信号特征提取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811241761.0A CN109492546B (zh) | 2018-10-24 | 2018-10-24 | 一种融合小波包与互信息的生物信号特征提取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109492546A CN109492546A (zh) | 2019-03-19 |
| CN109492546B true CN109492546B (zh) | 2022-08-12 |
Family
ID=65692647
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811241761.0A Active CN109492546B (zh) | 2018-10-24 | 2018-10-24 | 一种融合小波包与互信息的生物信号特征提取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109492546B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110123368B (zh) * | 2019-05-28 | 2021-11-26 | 南京邮电大学 | 一种基于pca串并融合的多路心音特征表征方法 |
| CN111110221A (zh) * | 2019-12-26 | 2020-05-08 | 天津大学 | 一种时-频-非线性多维体表胃电特征提取方法 |
| CN111000553B (zh) * | 2019-12-30 | 2022-09-27 | 山东省计算中心（国家超级计算济南中心） | 一种基于投票集成学习的心电数据智能分类方法 |
| CN111904415B (zh) * | 2020-09-07 | 2023-06-16 | 天津市人民医院 | 一种基于脑电信号评估大脑左右半球信息交互强度的方法 |
| CN112451857A (zh) * | 2020-12-10 | 2021-03-09 | 深圳先进技术研究院 | 视网膜神经节细胞激活状态的判别方法、存储介质和设备 |
| CN116975588B (zh) * | 2023-09-22 | 2023-12-19 | 太原理工大学 | 用于带式输送机的故障诊断方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101832471A (zh) * | 2010-04-19 | 2010-09-15 | 哈尔滨工程大学 | 一种信号识别分类方法 |
| CN103236258A (zh) * | 2013-05-06 | 2013-08-07 | 东南大学 | 基于巴氏距离最优小波包分解的语音情感特征提取方法 |
| CN104367317A (zh) * | 2014-10-15 | 2015-02-25 | 北京理工大学 | 多尺度特征融合的心电图心电信号分类方法 |
| CN105095675A (zh) * | 2015-09-07 | 2015-11-25 | 浙江群力电气有限公司 | 一种开关柜故障特征选择方法及装置 |
| CN105183813A (zh) * | 2015-08-26 | 2015-12-23 | 山东省计算中心（国家超级计算济南中心） | 基于互信息的用于文档分类的并行特征选择方法 |
| CN105426842A (zh) * | 2015-11-19 | 2016-03-23 | 浙江大学 | 基于支持向量机的表面肌电信号多类手部动作识别方法 |
| CN105446484A (zh) * | 2015-11-19 | 2016-03-30 | 浙江大学 | 一种基于隐马尔科夫模型的肌电信号手势识别方法 |
| CN106127191A (zh) * | 2016-07-06 | 2016-11-16 | 东南大学 | 基于小波包分解及逻辑回归的脑电分类方法 |
| CN106503670A (zh) * | 2016-11-03 | 2017-03-15 | 太原理工大学 | 基于相关性分析特征选择的j波检测分类方法 |
-
2018
- 2018-10-24 CN CN201811241761.0A patent/CN109492546B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101832471A (zh) * | 2010-04-19 | 2010-09-15 | 哈尔滨工程大学 | 一种信号识别分类方法 |
| CN103236258A (zh) * | 2013-05-06 | 2013-08-07 | 东南大学 | 基于巴氏距离最优小波包分解的语音情感特征提取方法 |
| CN104367317A (zh) * | 2014-10-15 | 2015-02-25 | 北京理工大学 | 多尺度特征融合的心电图心电信号分类方法 |
| CN105183813A (zh) * | 2015-08-26 | 2015-12-23 | 山东省计算中心（国家超级计算济南中心） | 基于互信息的用于文档分类的并行特征选择方法 |
| CN105095675A (zh) * | 2015-09-07 | 2015-11-25 | 浙江群力电气有限公司 | 一种开关柜故障特征选择方法及装置 |
| CN105426842A (zh) * | 2015-11-19 | 2016-03-23 | 浙江大学 | 基于支持向量机的表面肌电信号多类手部动作识别方法 |
| CN105446484A (zh) * | 2015-11-19 | 2016-03-30 | 浙江大学 | 一种基于隐马尔科夫模型的肌电信号手势识别方法 |
| CN106127191A (zh) * | 2016-07-06 | 2016-11-16 | 东南大学 | 基于小波包分解及逻辑回归的脑电分类方法 |
| CN106503670A (zh) * | 2016-11-03 | 2017-03-15 | 太原理工大学 | 基于相关性分析特征选择的j波检测分类方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109492546A (zh) | 2019-03-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109492546B (zh) | 一种融合小波包与互信息的生物信号特征提取方法 | |
| Al-Saegh et al. | Deep learning for motor imagery EEG-based classification: A review | |
| Chatterjee et al. | Ensemble learning approach to motor imagery EEG signal classification | |
| Singh et al. | Efficient wavelet families for ECG classification using neural classifiers | |
| Kumar et al. | OPTICAL+: a frequency-based deep learning scheme for recognizing brain wave signals | |
| Bhanumathi et al. | CNN based training and classification of MRI brain images | |
| Khatibi et al. | Proposing feature engineering method based on deep learning and K-NNs for ECG beat classification and arrhythmia detection | |
| CN111387975B (zh) | 基于机器学习的脑电信号识别方法 | |
| Khan et al. | Detection of dilated cardiomyopathy using pulse plethysmographic signal analysis | |
| Mahmoudi et al. | Multi-class EEG classification of motor imagery signal by finding optimal time segments and features using SNR-based mutual information | |
| Ahmad et al. | Deep learning based on CNN for emotion recognition using EEG signal | |
| Pan et al. | Recognition of human inner emotion based on two-stage FCA-ReliefF feature optimization | |
| Al-Qaysi et al. | Wavelet-based Hybrid learning framework for motor imagery classification | |
| Singh et al. | Short and noisy electrocardiogram classification based on deep learning | |
| Kanna et al. | Smart Detection and Removal of Artifacts in Cognitive Signals Using Biomedical Signal Intelligence Applications | |
| Phinyomark et al. | Applications of variance fractal dimension: A survey | |
| CN112450885B (zh) | 一种面向癫痫脑电的识别方法 | |
| Mondal et al. | EEG Signal Classification with Machine Learning model using PCA feature selection with Modified Hilbert transformation for Brain-Computer Interface Application | |
| Tripathi et al. | Permuted layer-based CNN for Emotion Detection with Multi-Modality Physiological Signals | |
| Khan et al. | Esad: Expert system for apnea detection using enhanced dwt feature extraction and machine learning algorithms | |
| Ali et al. | Hybridization of Wavelet Decomposition and Machine Learning for Brain Waves based Emotion Recognition | |
| Hashmi et al. | A Comparative Study of Machine Learning Algorithms for EEG Signal Classification | |
| Donmez et al. | Classification of Emotional States from Multi-channel EEG Signals by Convolutional Neural Networks | |
| Fouad et al. | Attempts towards the first brain-computer interface system in INAYA Medical College | |
| Saranya et al. | Early Detection of Heartbeat from Multimodal Data Using RPA Learning with KDNN-SAE. |
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 |