线性相位ＩＩＲ滤波器的一种实现方法技术

本发明专利技术公开了线性相位ＩＩＲ滤波器的一种实现方法，包括划分等长片段进行时序反转、零输入延长后进行反向滤波、再时序反转、前向滤波等步骤。本发明专利技术能够满足实时场合无限连续信号应用，并且在实现线性相位的同时大大的减少计算量，减小波形失真误差，保证很好的幅频特性。

【技术实现步骤摘要】

本专利技术属于信号处理领域，特别是线性相位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)从上式可以看出反向滤波后输出本文档来自技高网...

【技术保护点】
线性相位ＩＩＲ滤波器的一种实现方法，其特征在于包括以下步骤：（１）对无限长连续信号采样后得到ｘ（ｎ）＝｛ｘ（０），ｘ（１），ｘ（２），…｝，按照时序划分成长度为Ｌ的等长片段，对等长片段依次进行时序反转，反转后得到ａ（ｎ）＝｛ｘ（Ｌ－１），…，ｘ（０），ｘ（２Ｌ－１），…，ｘ（Ｌ），ｘ（３Ｌ－１），…，ｘ（２Ｌ），ｘ（４Ｌ－１），…｝；（２）在反转后得到的片段后面进行Ｌ个零输入的延长，对延长后所得的２Ｌ长度的新片段利用反向ＩＩＲ滤波器进行反向滤波，反向滤波后输出片段结果为当前片段的前Ｌ个滤波输出响应加上后一个片段的后Ｌ个滤波输出响应，即ｙ↓［ｋ］＝ｘ↓［ｋ］↑［Ｌ］＋ｘ↓［ｋ＋１］↑［Ｔ］，其中ｙ↓［ｋ］为反向滤波后输出的第ｋ个片段，ｘ↓［ｋ］↑［Ｌ］为第ｋ个片段的前Ｌ个滤波输出响应，ｘ↓［ｋ＋１］↑［Ｔ］为第ｋ＋１个片段后一个片段的后Ｌ个滤波输出响应，反向滤波的片段输出序列为ｆ（ｎ）＝｛ｙ↓［０］，ｙ↓［１］，ｙ↓［２］，…｝；（３）步骤（２）所得到的片段输出序列再按照步骤（１）的方法进行时序反转；（４）将步骤（３）所得到的片段序列利用前向ＩＩＲ滤波器进行前向滤波，得到输出结果。

【技术特征摘要】

【专利技术属性】
技术研发人员：刘伟，
申请(专利权)人：广东宝莱特医用科技股份有限公司，
类型：发明
国别省市：44[]