基于最优核时频分布可视图的癫痫脑电信号识别方法技术

一种基于自适应最优核时频分布可视图的癫痫脑电信号识别法，包括：获得原始脑电数据，计算预处理后的脑电信号的自适应最优核时频分布；进行癫痫状态的分类：对于每一个得到的时间能量序列，进行可视化分析，得到能量时间序列可视图复杂网络，提取获得能量时间序列可视图复杂网络的特征指标，提取自适应最优核时频分布指标，结合支持向量机分类癫痫EEG信号。本发明专利技术的基于最优核时频分布可视图的癫痫脑电信号识别方法，结合自适应最优核时频分布理论与可视图思想构建能量时间序列可视图复杂网络，提取复杂网络指标，实现对癫痫脑电信号的识别。

【技术实现步骤摘要】

本专利技术涉及一种癫痫脑电信号识别法。特别是涉及一种基于最优核时频分布可视图的癫痫脑电信号识别方法。
技术介绍
癫痫是指大脑神经元突发性异常放电，导致短暂的大脑功能障碍的一种慢性疾病。在临床上的表现为发作性运动、感觉、自主神经、意识及精神障碍。癫痫作为一种慢性疾病，虽然短期内对患者没有太大的影响，但是长期频繁的发作可导致患者的身心、智力产生严重影响。脑电信号是脑神经细胞电生理活动在大脑皮层或头皮表面的总体反映。脑电信号中包含了大量的生理与疾病信息，分析癫痫患者与健康人的脑电信号，可以更深入的认识癫痫疾病并获得大量的功能与疾病信息，从而可以根据这些信息对癫痫疾病诊断提供有效的方法。但由于脑电信号是不具备各态历经性的非平稳随机信号，而且其背景噪声很强，因此脑电信号的分析和处理一直是非常吸引人但又是具有相当难度的研究课题。可视图建网从其被创立之初至今已在多领域得到广泛应用，为时间序列分析提供了一个新视角。自适应最优核时频分布在时域内实现了对信号特征紧密的“跟踪”，可有效抑制交叉项产生的干扰，具有很强的信号细节刻画能力。本专利技术的一种基于自适应最优核时频分布可视图的癫痫脑电信号识别方法将时频分析与可视图理论结合，实现了对不同脑电信号的准确分类。
技术实现思路
本专利技术所要解决的技术问题是，提供一种基于最优核时频分布可视图的癫痫脑电信号识别法及应用。本专利技术所采用的技术方案是：一种基于自适应最优核时频分布可视图的癫痫脑电信号识别法，包括如下步骤：1)获得原始脑电数据，计算预处理后的脑电信号的自适应最优核时频分布，包括：(1)自适应最优核时频分布PAOK(t,f)表示为：PAOK(t,f)＝∫∫A(t,τ,υ)Φ(τ,υ)e-j2π(tυ+τf)dτdυ其中，υ为频偏，τ为时延，t为时间，f为频率，Φ(τ,υ)为用于产生自适应最优核时频分布的核函数，形式为二维高斯函数形式：其中，r2＝τ2+υ2，ψ＝arctan(τ/θ),σ(ψ)为控制核函数在ψ方向上的扩展，A(t,τ,υ)为自适应最优核时频分布的短时模糊函数，以时间t为中心，表示为：其中，z为所要计算的自适应最优核时频分布的脑电信号，ω为以时间t为中心，长度为T的对称窗函数，自适应最优核时频分布的短时模糊函数只在[t-T,t+T]范围内计算，u表示时间变量，自适应最优核核函数在极坐标下表示为Φ(r,ψ)，所述核函数由以下的最优化条件来决定：约束条件为：式中，α为核参数，A(r,ψ)和Φ(r,ψ)分别是极坐标下的短时模糊函数和核函数；(2)根据第(1)步所得到的自适应最优核时频分布PAOK(t,f)，得到脑电信号的能量时间序列所述的预处理包括：对原始脑电数据进行带通滤波，去除高低频干扰成分；手动剔除剩余伪迹数据，最终得到预处理后的脑电数据用于后续的分析；2)进行癫痫状态的分类，包括：(1)对于每一个得到的时间能量序列进行可视化分析，得到能量时间序列可视图复杂网络，具体是：如果能量时间序列内的任意两个点(ta,Ea)和(tb,Eb)可视，也就是两点之间有连边，则这两点间的任意一点(tc,Ec)都要满足:(2)提取获得能量时间序列可视图复杂网络的特征指标：全局聚集系数C：局部聚集系数均值聚集系数熵EC：节点度的均值其中，τ△表示的是一个可视图网络中闭三元组的数目，一个三角由三个封闭的三元组构成，τ表示一个可视图网络中的开三元组数目，开三元组即为三个节点由两条连边连接，τi表示一个可视图网络中以节点i为中心的开三元组数目，τi,△表示的是一个可视图网络中以节点i为中心的闭三元组数目，Ci表示节点i的局部聚集系数，ki是节点i的度，即节点i的连边数，N表示网络的节点数；(3)提取自适应最优核时频分布指标：时频分布偏差值Pstd：其中，m和n分别是时频分布的行数和列数，Pmean表示时频分布均值，自适应最优核时频分布的频率范围为(f0-F)；(4)结合支持向量机分类癫痫EEG信号，包括：组合自适应最优核时频分布指标和能量时间序列可视图复杂网络的特征指标，构成三维指标向量，分别采用留一法交叉验证和支持向量机分类器，对所有的三维指标向量进行分类，以及采用十折交叉验证和支持向量机分类器，对所有的三维指标向量进行分类。步骤2)中的第(4)步所述的采用留一法交叉验证和支持向量机分类器，对所有的三维指标向量进行分类，是将步骤2)中得到的每一个多尺度水平有限穿越可视图复杂网络的节点平均聚集系数和聚集系数熵EC构成一个三维指标向量，采用留一法交叉验证和支持向量机分类器，对所有的三维指标向量进行分类，即依次将一个三维指标向量作为测试集，将其余的指标向量作为训练集，并提供对应的训练集标签，遍历所有的三维指标向量作为测试集，其余的指标向量作为训练集后，从而得到所有的三维指标向量的分类结果。步骤2)中的第(4)步所述的采用十折交叉验证和支持向量机分类器，对所有的三维指标向量进行分类，包括：(1)将所有的三维指标向量随机的分为等长的十份子集；(2)将十份子集中的一个子集作为测试集，剩余的九份子集作为训练集，提供对应的训练集标签，得到作为测试集的子集中的三维指标向量的分类结果；(3)重复第(2)步，直至分别得到十份子集的三维指标向量的分类结果，构成一次的十折交叉验证分类结果；(4)重复第(1)步～第(3)步十次，得到十次的十折交叉验证分类结果；(5)对十次的十折交叉验证分类结果取均值作为所有三维指标向量的分类结果。本专利技术的基于最优核时频分布可视图的癫痫脑电信号识别方法，结合自适应最优核时频分布理论与可视图思想构建能量时间序列可视图复杂网络，提取复杂网络指标，实现对癫痫脑电信号的识别。附图说明图1是本专利技术基于最优核时频分布可视图的癫痫脑电信号识别法的流程图；图2是本专利技术应用于区分癫痫患者癫痫发作时期和健康人脑电信号的二维指标分布图；图3是本专利技术应用区分癫痫患者癫痫发作时期与正常时期脑电信号的二维指标分布图。具体实施方式下面结合实施例和附图对本专利技术的基于最优核时频分布可视图的癫痫脑电信号识别方法做出详细说明。本专利技术的基于最优核时频分布可视图的癫痫脑电信号识别方法，通过计算不同大脑状态下的脑电信号的自适应最优核时频分布，得到时间能量序列，对能量时间序列构建可视图复杂网络，提取复杂网络特征指标，结合机器学习实现对癫痫发作状态的辨别与分析。本专利技术的基于最优核时频分布可视图的癫痫脑电信号识别法，如图1所示，包括如下步骤：1)获得原始脑电数据，计算预处理后的脑电信号的自适应最优核时频分布，包括：(1)自适应最优核时频分布PAOK(t,f)表示为：PAOK(t,f)＝∫∫A(t,τ,υ)Φ(τ,υ)e-j2π(tυ+τf)dτdυ其中，υ为频偏，τ为时延，t为时间，f为频率，Φ(τ,υ)为用于产生自适应最优核时频分布的核函数，形式为二维高斯函数形式：其中，r2＝τ2+υ2，ψ＝arctan(τ/θ),σ(ψ)为控制核函数在ψ方向上的扩展，A(t,τ,υ)为自适应最优核时频分布的短时模糊函数，以时间t为中心，表示为：其中，z为所要计算的自适应最优核时频分布的脑电信号，ω为以时间t为中心，长度为T的对称窗函数，自适应最优核时频分布的短时模糊函数只在[t-T,t+T]范围内计算，u表示时间变量，自适应最优核核函数本文档来自技高网...

【技术保护点】
一种基于自适应最优核时频分布可视图的癫痫脑电信号识别法，其特征在于，包括如下步骤：1)获得原始脑电数据，计算预处理后的脑电信号的自适应最优核时频分布，包括：(1)自适应最优核时频分布PAOK(t,f)表示为：PAOK(t,f)＝∫∫A(t,τ,υ)Φ(τ,υ)e‑j2π(tυ+τf)dτdυ其中，υ为频偏，τ为时延，t为时间，f为频率，Φ(τ,υ)为用于产生自适应最优核时频分布的核函数，形式为二维高斯函数形式：Φ(τ,υ)=e-r2/2σ2(ψ)]]>其中，r2＝τ2+υ2，ψ＝arctan(τ/θ),σ(ψ)为控制核函数在ψ方向上的扩展，A(t,τ,υ)为自适应最优核时频分布的短时模糊函数，以时间t为中心，表示为：A(t,υ,τ)=∫z*(uτ2)·ω*(u-t-τ2)·z(u+τ2)·ω(u-t+τ2)·ejυudu]]>其中，z为所要计算的自适应最优核时频分布的脑电信号，ω为以时间t为中心，长度为T的对称窗函数，自适应最优核时频分布的短时模糊函数只在[t‑T,t+T]范围内计算，u表示时间变量，自适应最优核核函数在极坐标下表示为Φ(r,ψ)，所述核函数由以下的最优化条件来决定：maxΦ∫∫|A(r,ψ)·Φ(r,ψ)|2rdrdψ]]>约束条件为：Φ(r,ψ)=e-r2/2σ2(ψ)]]>12π∫∫|Φ(r,ψ)|2rdrdψ=12π∫σ2(ψ)dψ≤α,α≥0]]>式中，α为核参数，A(r,ψ)和Φ(r,ψ)分别是极坐标下的短时模糊函数和核函数；(2)根据第(1)步所得到的自适应最优核时频分布PAOK(t,f)，得到脑电信号的能量时间序列{Et}t=0T=Σf=f0FPAOK(t,f)]]>所述的预处理包括：对原始脑电数据进行带通滤波，去除高低频干扰成分；手动剔除剩余伪迹数据，最终得到预处理后的脑电数据用于后续的分析；2)进行癫痫状态的分类，包括：(1)对于每一个得到的时间能量序列进行可视化分析，得到能量时间序列可视图复杂网络，具体是：如果能量时间序列内的任意两个点(ta,Ea)和(tb,Eb)可视，也就是两点之间有连边，则这两点间的任意一点(tc,Ec)都要满足:Ec<Eb+(Ea-Eb)tb-tctb-ta]]>(2)提取获得能量时间序列可视图复杂网络的特征指标：全局聚集系数C：C=τΔτ]]>局部聚集系数均值Ci=τi,Δτi]]>C‾=1NΣi=1NCi]]>聚集系数熵EC：EC=-Σi=1N(Ci/Σi=1NCi)log(Ci/Σi=1NCi)]]>节点度的均值k‾=1NΣi=1Nki]]>其中，τ△表示的是一个可视图网络中闭三元组的数目，一个三角由三个封闭的三元组构成，τ表示一个可视图网络中的开三元组数目，开三元组即为三个节点由两条连边连接，τi表示一个可视图网络中以节点i为中心的开三元组数目，τi,△表示的是一个可视图网络中以节点i为中心的闭三元组数目，Ci表示节点i的局部聚集系数，ki是节点i的度，即节点i的连边数，N表示网络的节点数；(3)提取自适应最优核时频分布指标：时频分布偏差值Pstd：Pstd=Σt=0TΣf=f0F(PAOK(t,f)-Pmean)2mn-1]]>Pmean=Σt=0TΣf=f0FPAOK(t,f)mn]]>其中，m和n分别是时频分布的行数和列数，Pmean表示时频分布均值，自适应最优核时频分布的频率范围为(f0‑F)；(4)结合支持向量机分类癫痫EEG信号，包括：组合自适应最优核时频分布指标和能量时间序列可视图复杂网络的特征指标，构成三维指标向量，分别采用留一法交叉验证和支持向量机分类器，对所有的三维指标向量进行分类，以及采用十折交叉验证和支持向量机分类器，对所有的三维指标向量进行分类。...

【技术特征摘要】
1.一种基于自适应最优核时频分布可视图的癫痫脑电信号识别法，其特征在于，包括如下步骤：1)获得原始脑电数据，计算预处理后的脑电信号的自适应最优核时频分布，包括：(1)自适应最优核时频分布PAOK(t,f)表示为：PAOK(t,f)＝∫∫A(t,τ,υ)Φ(τ,υ)e-j2π(tυ+τf)dτdυ其中，υ为频偏，τ为时延，t为时间，f为频率，Φ(τ,υ)为用于产生自适应最优核时频分布的核函数，形式为二维高斯函数形式：Φ(τ,υ)=e-r2/2σ2(ψ)]]>其中，r2＝τ2+υ2，ψ＝arctan(τ/θ),σ(ψ)为控制核函数在ψ方向上的扩展，A(t,τ,υ)为自适应最优核时频分布的短时模糊函数，以时间t为中心，表示为：A(t,υ,τ)=∫z*(uτ2)·ω*(u-t-τ2)·z(u+τ2)·ω(u-t+τ2)·ej&upsi...

【专利技术属性】
技术研发人员：高忠科，蔡清，张珊珊，杨宇轩，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津;12