主权项 |
一种用于解剖学定位的Theta节律信号延时算法,其特征在于,包括下述步骤:(1)通过对大的插值倍数进行分解,形成若干个小的插值倍数的级联形式:首先确定插值倍数Md,并确定数字延时值d=delay×Md/T,T为采样信号周期,然后将插值倍数Md进行分解,形成至少三个小插值倍数的级联形式:Md=M1×M2×…×Mp M1≥M2≥…≥Mp其中,p为≥3的正整数,M1、M2、…、Mp是分解后的小插值倍数;(2)通过数学推导,得到级联延时算法的数学表达式:设与第一级插值倍数M1对应的FIR低通滤波器为hL1(n),系数为[hL1(0),hL1(1),…,hL1(N1‑1)],长度为N1;与第二级插值倍数M2对应的FIR低通滤波器为hL2(n),系数为[hL2(0),hL2(1),…,hL2(N2‑1)],长度为N2;……,与第p级插值倍数Mp对应的FIR低通滤波器为hLp(n),系数为[hLp(0),hLp(1),…,hLp(Np‑1)],长度为Np;由hL1(n),hL2(n),…,hLp(n)可得,延时值为d的p级级联延时算法的数学表达式为: <mrow> <msub> <mi>h</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <msub> <mi>l</mi> <mn>1</mn> </msub> <mo>=</mo> <mi>Hi</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>Lo</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </munderover> <msub> <mi>h</mi> <mrow> <mi>L</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>nM</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>l</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>×</mo> <munderover> <mo>∏</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>2</mn> </mrow> <mrow> <mi>p</mi> <mo>-</mo> <mn>2</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <msub> <mi>l</mi> <mi>i</mi> </msub> <mo>=</mo> <mi>Hi</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>Lo</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </munderover> <msub> <mi>h</mi> <mi>Li</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>l</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <msub> <mi>M</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>l</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <munderover> <mi>Σ</mi> <mrow> <msub> <mi>l</mi> <mrow> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mi>Hi</mi> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>Lo</mi> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> </munderover> <msub> <mi>h</mi> <mi>L</mi> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>l</mi> <mrow> <mi>p</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> <msub> <mi>M</mi> <mrow> <mo>(</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>l</mi> <mrow> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <msub> <mi>h</mi> <mi>Lp</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>l</mi> <mrow> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <msub> <mi>M</mi> <mi>p</mi> </msub> <mo>-</mo> <mi>d</mi> <mo>)</mo> </mrow> </mrow>其中, <mrow> <mi>Hi</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>High</mi> <mo>[</mo> <mfrac> <mi>d</mi> <mrow> <munderover> <mo>∏</mo> <mrow> <mi>j</mi> <mo>=</mo> <mi>i</mi> </mrow> <mi>p</mi> </munderover> <msub> <mi>M</mi> <mi>j</mi> </msub> </mrow> </mfrac> <mo>]</mo> <mo>,</mo> </mrow> <mrow> <mi>Lo</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>Low</mi> <mo>[</mo> <munderover> <mi>Σ</mi> <mrow> <mi>k</mi> <mo>=</mo> <mi>i</mi> </mrow> <mrow> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mfrac> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mi>k</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mrow> <munderover> <mo>∏</mo> <mrow> <mi>j</mi> <mo>=</mo> <mi>i</mi> </mrow> <mi>k</mi> </munderover> <msub> <mi>M</mi> <mi>j</mi> </msub> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mi>p</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>+</mo> <mi>d</mi> <mo>)</mo> </mrow> <mrow> <munderover> <mo>∏</mo> <mrow> <mi>j</mi> <mo>=</mo> <mi>i</mi> </mrow> <mi>p</mi> </munderover> <msub> <mi>M</mi> <mi>j</mi> </msub> </mrow> </mfrac> <mo>]</mo> <mo>,</mo> </mrow>i取正整数;High[A]表示大于A的最小整数;Low[A]表示小于A的最大整数;hd(n)中n的有效取值范围为:Hi(1)≤n≤Lo(1)上式中, <mrow> <mi>Hi</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>High</mi> <mo>[</mo> <mfrac> <mi>d</mi> <mrow> <munderover> <mo>∏</mo> <mrow> <mi>j</mi> <mo>=</mo> <mi>i</mi> </mrow> <mi>p</mi> </munderover> <msub> <mi>M</mi> <mi>j</mi> </msub> </mrow> </mfrac> <mo>]</mo> <mo>,</mo> </mrow> <mrow> <mi>Lo</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>Low</mi> <mo>[</mo> <munderover> <mi>Σ</mi> <mrow> <mi>k</mi> <mo>=</mo> <mi>i</mi> </mrow> <mrow> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mfrac> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mi>k</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mrow> <munderover> <mo>∏</mo> <mrow> <mi>j</mi> <mo>=</mo> <mi>i</mi> </mrow> <mi>k</mi> </munderover> <msub> <mi>M</mi> <mi>j</mi> </msub> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mi>p</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>+</mo> <mi>d</mi> <mo>)</mo> </mrow> <mrow> <munderover> <mo>∏</mo> <mrow> <mi>j</mi> <mo>=</mo> <mi>i</mi> </mrow> <mi>p</mi> </munderover> <msub> <mi>M</mi> <mi>j</mi> </msub> </mrow> </mfrac> <mo>]</mo> </mrow>(3)通过频谱分析方法,消除级联延时算法中的附加时延:根据hd(n)的表达式,引入一个0延时滤波器h0(n),其离散时间傅里叶变换为: <mrow> <msub> <mi>H</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <msup> <mi>e</mi> <mi>jω</mi> </msup> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>n</mi> <mo>=</mo> <mo>-</mo> <mo>∞</mo> </mrow> <mo>∞</mo> </munderover> <msub> <mi>h</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>jωn</mi> </mrow> </msup> </mrow>并且: <mrow> <msub> <mi>H</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <msup> <mi>e</mi> <mi>jω</mi> </msup> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>n</mi> <mo>=</mo> <mo>-</mo> <mo>∞</mo> </mrow> <mo>∞</mo> </munderover> <msub> <mi>h</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>jωn</mi> </mrow> </msup> </mrow>令H‑0(ejω)=H0(e‑jω),则消除了附加时延的级联插值‑滤波‑延时‑抽取算法的滤波器为:Hdd(ejω)=Hd(ejω)×H‑0(ejω)=|H0(ejω)|e‑jωd即:hdd(n)=hd(n)*h‑0(n)。 |