专利名称：线性相位iir滤波器的一种实现方法
技术领域：
本发明属于信号处理领域，特别是线性相位nR滤波器的一种实现方法。
背景技术：
为了避免各种干扰对信号的采集记录过程产生影响，一般需要设置滤波器对这些 干扰进行滤除。在实际滤除过程中，往往希望不改变波形或减少波形失真，因而如果设计的 滤波器有线性相位响应，则这会是它们很理想的特性，会使波形失真很少。在陷波器的设计 过程亦是如此。通常HR滤波器很容易满足线性相位特性，但它的计算量相当的大；用简单 整系数方法设计的窄带工频陷波器也能很好满足线性相位特性，但它有较大延时，并且会 要求采样率是陷波点的整数倍，这对不同频带陷波器的设计也是一个制约；而常规方法设 计的IIR陷波器不能有较精确的线性相位特性，因为它的零点和极点不会镜像对称，如果 要求极点关于零点镜像对称，则极点会在单位圆的外部，造成滤波器不稳定。后来有方法上的改进，比如，通过全通滤波器来实现IIR滤波的相位均等，但它的 计算量有时比传统HR滤波器还有所增加；通过等波纹及最大平坦群延时方法设计IIR滤 波器，但它没有理想的幅频特性；还有报道较多的传统的零相位IIR滤波器，是由时序反转 及多路滤波组合而成，通过递归计算，保持了幅频特性以及线性相位，但是它主要是针对有 限长序列反转，并且反向UR滤波不能实时实现，限制了其在连续信号等实时性要求较高 的场合的应用。
发明内容
为解决上述问题，本发明提供线性相位UR滤波器的一种实现方法，本方法能够 满足实时场合无限连续信号应用，并且在实现线性相位的同时大大的减少计算量，减小波 形失真误差，保证很好的幅频特性。本发明为解决其问题所采用的技术方案是
线性相位HR滤波器的一种实现方法，包括以下步骤
(1)对无限长连续信号采样后得到x(n)= {X(0)，X(l)，d2)，…}，按照时序划分成 长度为L的等长片段，对等长片段依次进行时序反转，反转后得到a (n) = Ix (L-I)，… ,x(0),x (2L-1)，…，χ (L)，χ (3L-1)，…，χ (2L)，χ (4L-1)，…}；
(2)在反转后得到的片段后面进行L个零输入的延长，对延长后所得的2L长度的新片 段利用反向HR滤波器进行反向滤波，反向滤波后输出片段结果为当前片段的前L个滤波 输出响应加上后一个片段的后L个滤波输出响应，即yk=Xk、Xk+1T，其中yk为反向滤波后输出 的第k个片段，xkL为第k个片段的前L个滤波输出响应，xk+1T为第k+Ι个片段后一个片段 的后L个滤波输出响应，反向滤波的片段输出序列为 ·(η) = {、Υι，Υ2，···}；
(3)步骤(2)所得到的片段输出序列再按照步骤(1)的方法进行时序反转；
(4)将步骤(3)所得到的片段序列利用前向HR滤波器进行前向滤波，得到输出结果。其中，所述反向UR滤波器和前向UR滤波器的传递函数相同。进一步优选的是，所述反向HR滤波器和前向UR滤波器采用相同的2阶IIR Butterworth滤波器。步骤(1)中时序反转的方法包括以下步骤
a.预设一长度为L的寄存器，寄存器的两端分别为第一输入输出端和第二输入输出
端；
b.χ (η)第一个片段的L个信号x(0)，···，x(L)从第一输入输出端依次存入寄存器；
c.第一个片段的L个信号从第一输入输出端依次取出得到x(L-l)，…，χ(0)，第二个片 段的L个信号χ (L)，…，χ (2L-1)从第二输入输出端依次存入寄存器；
d.第二个片段的L个信号从第二输入输出端依次取出得到χ(2L-1)，…，χ (L)，第三个 片段的L个信号χ (2L)，…，χ (3L-1)从第一输入输出端依次存入寄存器；
e.依照上述步骤持续循环。本发明将信号划分为等长片段后再进行时序反转，这种反转的方法可以满足无线 信号场合的应用，保证实时性要求；针对分片处理后所导致的截断误差，本发明在反转后得 到的片段后面进行零输入的延长，使得相邻片段反向滤波后的输出的叠加作为当前片段的 反向滤波输出，实现片段滤波截断误差的补偿修正，从而减小波形失真误差；另外，通过对 传递函数进行分析可知，基于本发明的方法制作的滤波器具有良好的幅频特性和线性相位 特性，本发明的方法与HR滤波器的实现方法以及通过全通滤波器来实现IIR滤波的相位 均等的方法相比计算量大大减少。
下面结合附图和实施例对本发明作进一步说明 图1为基于本发明的线性IIR滤波器的系统框图2为划分等长片段进行时序反转的一种实施例的示意图； 图3为反向滤波输出响应的一种实施例的示意图4为一种常规IIR陷波器和本发明的线性UR滤波器幅频特性的对比图； 图5为本发明的线性IIR滤波器的相位特性图6为常规IIR陷波器和本发明的线性IIR滤波器分别对50Hz信号滤波后的信号特 性对比图7为本发明的线性HR滤波器对有工频干扰的心电信号进行陷波的效果。
具体实施例方式本发明的线性相位IIR滤波器的实现方法包括划分等长片段进行时序反转、零输 入延长后进行反向滤波、再时序反转、前向滤波等步骤。下面分别详细描述。步骤(1):划分等长片段进行时序反转。对无限长连续信号采样后得到 X (η) = {χ (0)，χ⑴，χ (2)，···}，按照时序划分成长度为L的等长片段，对等长片段依次进行 时序反转。该时序反转的方法一般采用以下步骤
a.预设一长度为L的寄存器，寄存器的两端分别为第一输入输出端和第二输入输 出端；b. χ (η)第一个片段的L个信号x(0)，···，x(L)从第一输入输出端依次存入寄存 器；c.第一个片段的L个信号从第一输入输出端依次取出得到x(L-l)，…，X(O)，第二个 片段的L个信号x(L)，…，W2L-1)从第二输入输出端依次存入寄存器；d.第二个片段的L个信号从第二输入输出端依次取出得到x(2L-l)，…，X(L)，第三个片段的L个信号 x(2L)，…，x(3L-l)从第一输入输出端依次存入寄存器；e.依照上述步骤持续循环。反转 后得到 a (η) = {χ (L-I)，…，χ (0)，χ (2L-1)，…，χ (L)，χ (3L-1)，…，χ (2L)，χ (4L-1)，...}。利 用上述反转方法进行时序反转的具体实施例参照图2，在此实施例中，寄存器的长度L=5， 经反转后得到的序列为 a(n) = {x(4),x(3),x(2),x(l),x(0),x(9),x(8),x(7),x(6),x(5), x(14)，…}。从上面可以看出，本发明的时序反转方法可以无限的连续不断的对信号片段进 行时序反转，能够满足实时场合无限连续信号应用。步骤(2)零输入延长后进行反向滤波。在反转后得到的片段后面进行L个零输 入的延长，对延长后所得的2L长度的新片段利用反向UR滤波器进行反向滤波。假设反向 IIR滤波器的传递函数为H(ζ)，反向滤波后的输出yk(n) =h (-η) *χ (η)，其中kL彡η彡(k+1) L-1，上式进行下面的形式变换
Yk (η)
=h (-η) (η) =Σ n+L m=nx(m)h(n_m)
=Σ (k+1) L-I m=nx (m) h (n-m) + Σ n+L m= (k+1) Lx (m)h(n-m) =h (-n) *xk (n) +h (_n) *xk+1 (n)
从上式可以看出反向滤波后输出片段结果为当前片段的前L个滤波输出响应加上后 一个片段的后L个滤波输出响应，即yk=Xk、Xk+1T，其中yk为反向滤波后输出的第k个片段， xkL为第k个片段的前L个滤波输出响应，xk+1T为第k+Ι个片段后一个片段的后L个滤波输 出响应。这种采用相邻片段滤波输出叠加的结果作用当前片段的输出的方法可以看作是对 步骤(1)中对信号分段截断而产生的误差的补偿修正，能够很好地减少和消除周期误差。经 过步骤(2)的零输入延长及反向滤波处理后的片段输出序列为f (η) = {、Y1, y2,…}。图3 为一种反向滤波输出响应的实施例的示意图，本实施例中的每个片段的长度L=5，从图中可 以看出，反向滤波输出的每个片段都是由两部分相加得到的。步骤(3)将步骤(2)所得到的片段输出序列再按照步骤(1)的方法进行时序反 转。片段输出序列f(n)H^ypy2,…}，按照时序划分成长度为L的等长片段，对等长片段 依次进行时序反转，反转后得到 b (n) = {y4, y3, y2, Y1, y0, y9, y8, Y7, Y6, Y5, Yi4,…}。步骤(4):将步骤(3)所得到的片段序列利用前向UR滤波器进行前向滤波，得 到输出结果。作为优选的实施方式，该前向UR滤波器和反向UR滤波器的传递函数相 同，在一般的陷波器系统中，该前向UR滤波器和反向UR滤波器优选采用相同的2阶IIR Butterworth 滤波器。参照图1，为基于本发明的线性IIR滤波器的系统框图。信号输入序列为x(n)， 输出序列为y (η)，本发明的步骤(1)于图中的A部分进行，步骤(2)于图中的B部分进 行，步骤(3)于图中的C部分进行，步骤(4)于图中的D部分进行。假设B部分中的反向 IIR滤波器的传递函数为Ha(z)，D部分中的前向HR滤波器的传递函数为Hb(z)，那么整 个滤波器的系统传递函数H^1 (Z)=Hb(z)Ha(P)^(H)/2,系统的传递函数的幅度|Hjz) | = |Hb(Z) M^Uz—1) I，当反向IIR滤波器和前向IIR滤波器的传递函数为相同的H(Z)时，系统 传递函数H^1⑴劝⑴对广)^^2，幅度iHjz) I= IH(ζ) |2，即系统滤波器的衰减深度为 滤波器H(ζ)衰减深度的2倍。
下面以一去除50Hz工频干扰的系统陷波器实例来说明本发明的效果。其中用到 的前向HR滤波器和反向UR滤波器都为2阶IIR Butterworth滤波器，此滤波器采用常 规IIR陷波器设计方法设计，截止频率分别为48Hz和52Hz。参照图4-图6，其中图4为常 规IIR陷波器和系统幅频特性的对比图，线10为常规IIR陷波器H(Z)的幅频特性曲线，线 20为系统的幅频特性曲线，从图中可以看出，系统对50Hz工频的衰减效果比常规IIR陷波 器更好。图5为系统的相位特性图，从图中可以看出，系统的相位特性为线性。图6为常规 UR陷波器和系统滤波器分别对50Hz信号滤波后的信号特性对比图，图中幅度接近0. 02的 部分为利用常规IIR陷波器进行滤波后的信号，幅度大大小于0. 005接近零点的部分为利 用本系统进行滤波后的信号，从图中可以看出，本系统对50Hz工频干扰信号的滤除效果远 远优于常规IIR陷波器。此外，实验进一步表明，本系统与常规IIR陷波器相比，其对方波 进行滤波时引起的方波过冲震荡和波形失真误差会更小。本发明可以应用于对心电信号的处理方面。利用本发明的滤波器对有工频干扰的 心电信号进行陷波的效果可参照图7，其中图中的灰色线M为原始输出信号，黑色线N为经 过滤波后的信号，从图中可知，经过滤波后原始的信号的工频纹波基本被消除，所输出的心 电信号更干净更直观。
权利要求
1.线性相位IIR滤波器的一种实现方法，其特征在于包括以下步骤(1)对无限长连续信号采样后得到X(n)= {x(0)，x(l)，d2)，…}，按照时序划分成 长度为L的等长片段，对等长片段依次进行时序反转，反转后得到a (n) = Ix (L-I)，… ,x(0),x (2L-1)，…，χ (L)，χ (3L-1)，…，χ (2L)，χ (4L-1)，…}；(2)在反转后得到的片段后面进行L个零输入的延长，对延长后所得的2L长度的新片 段利用反向HR滤波器进行反向滤波，反向滤波后输出片段结果为当前片段的前L个滤波 输出响应加上后一个片段的后L个滤波输出响应，即yk=Xk、Xk+1T，其中yk为反向滤波后输出 的第k个片段，xkL为第k个片段的前L个滤波输出响应，xk+1T为第k+Ι个片段后一个片段 的后L个滤波输出响应，反向滤波的片段输出序列为 ·(η) = {、Υι，Υ2，···}；(3)步骤(2)所得到的片段输出序列再按照步骤(1)的方法进行时序反转；(4)将步骤(3)所得到的片段序列利用前向HR滤波器进行前向滤波，得到输出结果。
2.根据权利要求1所述的线性相位HR滤波器的一种实现方法，其特征在于所述反向 UR滤波器和前向UR滤波器的传递函数相同。
3.根据权利要求2所述的线性相位IIR滤波器的一种实现方法，其特征在于所述反向 IIR滤波器和前向UR滤波器采用相同的2阶IIR Butterworth滤波器。
4.根据权利要求1所述的线性相位IIR滤波器的一种实现方法，其特征在于步骤(1) 中时序反转的方法包括以下步骤a.预设一长度为L的寄存器，寄存器的两端分别为第一输入输出端和第二输入输出端；b.χ (η)第一个片段的L个信号x(0)，···，x(L)从第一输入输出端依次存入寄存器；c.第一个片段的L个信号从第一输入输出端依次取出得到χ(L-I)，…，χ (0)，第二个片 段的L个信号χ (L)，…，χ (2L-1)从第二输入输出端依次存入寄存器；d.第二个片段的L个信号从第二输入输出端依次取出得到χ(2L-1)，…，χ (L)，第三个 片段的L个信号W2L)，…，x(3L-l)从第一输入输出端依次存入寄存器；e.依照上述步骤持续循环。
全文摘要
本发明公开了线性相位IIR滤波器的一种实现方法，包括划分等长片段进行时序反转、零输入延长后进行反向滤波、再时序反转、前向滤波等步骤。本发明能够满足实时场合无限连续信号应用，并且在实现线性相位的同时大大的减少计算量，减小波形失真误差，保证很好的幅频特性。
文档编号H03H17/02GK102082559SQ20101056917
公开日2011年6月1日 申请日期2010年12月2日 优先权日2010年12月2日
发明者刘伟 申请人:广东宝莱特医用科技股份有限公司