[發明專利]一種用于解剖學定位的Theta 節律信號延時算法

無效

• 申請號: 201310117626.6
• 申請日: 2013-04-07
• 公開（公告）號: CN103218526A
• 公開（公告）日: 2013-07-24
• 發明（設計）人: 郭偉;羅二平;謝康寧;申廣浩;吳小明;湯池;劉娟;閆一力;孫濤
• 申請（專利權）人: 中國人民解放軍第四軍醫大學
• 主分類號: G06F19/00
• 分類號: G06F19/00
• 代理公司: 西安通大專利代理有限責任公司 61200
• 代理人: 朱海臨

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。