The invention discloses a method for detecting the effect connectivity between EEG signals based on three-dimensional non-linear partial direct coherence function, which includes the following steps: (1) constructing a single-input multi-output non-linear autoregressive model; (2) estimating the coefficients of the model constructed by step (1) using FROLS algorithm; (3) transforming the three-dimensional PDS formally to obtain the P of YJ with the signal Yi described by frequency response function. Definition formula of DC; (4) Using multi-dimensional Fourier transform of Volterra series kernel function to analyze SIMO NARX model in frequency domain and calculate the non-linear frequency response function of the model; (5) Substituting the non-linear frequency response function calculated in step (4) into the PDC definition formula in step (3), the three-dimensional NPDC is obtained, and the cause of one signal to another signal under the condition of considering three-dimensional signal at the same time is obtained. Fruit effect. This method can detect causality between three-dimensional EEG signals.
【技术实现步骤摘要】
基于三维非线性偏直接相干函数的脑电信号间效应连通性检测方法
本专利技术本专利技术属于生物医学领域,具体涉及一种检测脑电信号间效应连通性的方法。
技术介绍
大脑不同区域信号之间的效应连通性对于确定病灶区域具有重要作用。目前对脑电信号间效应连通性的检测包括时域分析和频域分析。PDC(PartialDirectedCoherence,偏直接相干函数)是一种常用的从频域上分析信号之间因果关系的方法,该方法可以辨别出两个信号之间直接或者间接的因果关系。但该方法是基于线性模型的,而真实的EEG(ElectroEncephaloGraphy,脑电图)是具有非线性特性的,PDC算法无法完全检测到信号中非线性因果关系。NPDC(NonlinearPDC,非线性PDC)是基于PDC所改进的适用于NARX(NonlinearAutoRegressiveeXogenous,非线性自回归模型)的方法,它可以检测到信号之间非线性的相互影响。目前常用的为二维NPDC,二维NPDC算法是一种基于二维PDC算法改进的适用于二维NARX模型的大脑效应连通性算法,该算法能够检测信号之间的线性与非线性因果关系。但是二维NPDC算法仅能处理二维信号,处理多维信号时无法区分信号之间的直接和间接因果关系。
技术实现思路
专利技术目的:针对现有技术中的不足,本专利技术提供了一种基于三维NPDC的脑电信号效应连通性检测方法,可以检测三维信号之间的因果关系。技术方案:本专利技术采用如下技术方案:基于三维非线性偏直接相干函数的脑电信号间效应连通性检测方法,包括如下步骤:(1)构造单输入多输出的非线性自回归模型,所 ...
【技术保护点】
1.基于三维非线性偏直接相干函数的脑电信号间效应连通性检测方法,其特征在于,包括如下步骤:(1)构造单输入多输出的非线性自回归模型,所述模型为:
【技术特征摘要】
1.基于三维非线性偏直接相干函数的脑电信号间效应连通性检测方法,其特征在于,包括如下步骤:(1)构造单输入多输出的非线性自回归模型,所述模型为:其中为脑电信号在时刻k的采样值,N为模型的非线性阶数,Nout为输出信号的个数,q、p分别是输入信号u(k-ki)和输出信号的非线性程度,且p+q≤N;ki为信号的延迟值,K为模型的延迟阶数,ki≤K,是u(k-ki)和的线性或非线性组合的模型项,为模型项的系数,为脑电信号的采样值和模型预测值之间的误差项;(2)应用FROLS算法对步骤(1)构造的模型进行系数估计;(3)对三维PDS进行形式变换,得到用频率响应函数描述的信号yi对yj的PDC的定义式;(4)应用Volterra级数核函数的多维傅里叶变换对SIMONARX模型进行频域分析,计算出模型的非线性频率响应函数;(5)将步骤(4)计算出的非线性频率响应函数代入步骤(3)中的PDC定义式,得到三维NPDC,得出在同时考虑三维信号的情况下某一信号对另一信号的因果影响。2.根据权利要求1所述的基于三维非线性偏直接相干函数的脑电信号间效应连通性检测方法,其特征在于,所述步骤(2)具体包括:(2-1)将非线性自回归模型改写为线性参数形式:其中pl(n)是脑电信号y(n-ki)和u(n-ki)的线性或非线性组合的模型项,L为候选项个数,θl为线性模型系数,e(n)为线性模型误差项;(2-2)将步骤(2-1)中的线性参数形式模型转换为正交模型:其中wl(n)是相互正交的,gl为正交模型的系数;(2-3)令D={p1,p2,…,pL}为由L个候选基底所组成的初始字典,pl=[pl(1),pl(2),…,pl(N)]T,令ql=pl且σ=yTy,l=1,2,…,L,计算向量y为由脑电信号采样值构成的向量;令即ERR[m1]=max{ERR(1)[l]:1≤l≤L},则第一个重要的模型项被选出,第一个正交向量可以被选为(2-4)假设在算法的第s-1步,选出了一个子集Ds-1,由s-1个重要的模型项组成,所述s-1个模型项经过正交变换组成了新的正交向量q1,q2,…,qs-1;在算法的第s步,令l≠m1,l≠m2,…,l≠ms-1,对于l=1,2,…,L,计算:令那么第s个重要的模型项就可以被选出来,第s个正交向量可以经过正交变换得到(2-5)重复步骤(2-4),当被筛选出的所有模型项的ERR值的和达到预先设定的阈值,就停止筛选过程;(2-6)设经过上述步骤从L个候选项中选出的L0个重要模型项的线性组合,L0<L,即正交模型为:其等价于其中参数向量可由计算得到。3.根据权利要求1所述的基于三维非线性偏直接相干函数的脑电信号间效应连...
【专利技术属性】
技术研发人员:杨淳沨,杨文琪,刘彦超,伍家松,孔佑勇,姜龙玉,杨冠羽,舒华忠,
申请(专利权)人:东南大学,
类型:发明
国别省市:江苏,32
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。