[发明专利]一种用于解剖学定位的Theta 节律信号延时算法

权利要求书

1.一种用于解剖学定位的Theta节律信号延时算法，其特征在于，包括下述步骤：

（1）通过对大的插值倍数进行分解，形成若干个小的插值倍数的级联形式：

首先确定插值倍数M_d，并确定数字延时值d=dela_y×M_d/T，T为采样信号周期，然后将插值倍数M_d进行分解，形成至少三个小插值倍数的级联形式：

M_d=M₁×M₂×…×M_p    M₁≥M₂≥…≥M_p

其中，p为≥3的正整数，M1、M2、…、M_p是分解后的小插值倍数；

（2）通过数学推导，得到级联延时算法的数学表达式：

设与第一级插值倍数M₁对应的FIR低通滤波器为h_L1(n)，系数为[h_L1(0),h_L1(1),…,h_L1(N₁-1)]，长度为N₁；与第二级插值倍数M₂对应的FIR低通滤波器为h_L2(n)，系数为[h_L2(0),h_L2(1),…,h_L2(N₂-1)]，长度为N₂；……，与第p级插值倍数M_p对应的FIR低通滤波器为h_Lp(n)，系数为[h_Lp(0),h_Lp(1),…,h_Lp(N_p-1)]，长度为N_p；由h_L1(n)，h_L2(n)，…，h_Lp(n)可得，延时值为d的p级级联延时算法的数学表达式为：

hd(n)=Σl1=Hi(2)Lo(2)hL1(nM1-l1)×∏i=2p-2Σli=Hi(i+1)Lo(i+1)hLi(li-1Mi-li)×Σlp-1=Hi(p)Lo(p)hL(p-1)(lp-2M(p-1)-lp-1)hLp(lp-1Mp-d)]]>

其中，Hi(i)=High[d∏j=ipMj],]]>Lo(i)=Low[Σk=ip-1(Nk-1)∏j=ikMj+(Np-1+d)∏j=ipMj],]]>i取正整数；

High[A]表示大于A的最小整数；Low[A]表示小于A的最大整数；

h_d(n)中n的有效取值范围为：

Hi(1)≤n≤Lo(1)

上式中，Hi(i)=High[d∏j=ipMj],]]>Lo(i)=Low[Σk=ip-1(Nk-1)∏j=ikMj+(Np-1+d)∏j=ipMj]]]>

（3）通过频谱分析方法，消除级联延时算法中的附加时延：

根据h_d(n)的表达式，引入一个0延时滤波器h₀(n)，其离散时间傅里叶变换为：

H0(ejω)=Σn=-∞∞h0(n)e-jωn]]>

并且：

Hd(ejω)=Σn=-∞∞hd(n)e-jωn]]>

令H_-0(e^jω)=H₀(e^-jω)，则消除了附加时延的级联插值-滤波-延时-抽取算法的滤波器为：

H_dd(e^jω)=H_d(e^jω)×H_-0(e^jω)=|H₀(e^jω)|e^-jωd

即：

h_dd(n)=h_d(n)*h_-0(n)。

2.如权利要求1所述的用于解剖学定位的Theta节律信号延时算法，其特征在于，所述p的取值范围为3～9。