zoukankan      html  css  js  c++  java
  • 数字滤波器

      滤波是信号处理中的一种基本而重要的技术,利用滤波技术可以在形形色色的信号中提取所需的信号,抑制不需要的干扰信号。按处理信号不同,滤波器可分为模拟滤波器与数字滤波器两大类。模拟滤波器是用来处理模拟信号或连续时间信号,数字滤波器是用来处理离散的数字信号。与模拟滤波器相比,数字滤波器具有诸多优点:可以用软件编程;稳定性高,可预测;不会因温度、湿度的影响产生误差,不需要高精度元件;很高的性能价格比。

      滤波器是频率选择电路,只允许输入信号中的一定频率成分通过,与模拟滤波器相同,按通带来划分可将数字滤波器分为4种类型:

    (1)低通滤波器(Lowpass filter) :通过低频,滤除高频。
    (2)高通滤波器(Highpass filter):通过高频,滤除低频。
    (3)带通滤波器(Bandpass filter):通过固定范围的频率。
    (4)带阻滤波器(Bandstop filter):滤除固定范围的频率。

    以上 4 种滤波器的理想频率响应如下图所示。
      滤波器允许通过的频率范围为通频带理想滤波器通频带的增益为1,所以信号的幅值不截止频带是滤波器不允许通过的频率范围。理论上滤波器在通带应有单位增益(0 dB),在阻带有0增益(-∞ dB),然而在实际中通
    阻带之间有一个过渡范围,如下图所示在这一范围增益在0~1之间。在应用中允许增益在单位增益上下轻微变化,即允许有通带纹波这一点也是实际滤波器与理想滤波器的区别。阻带衰减在实际中不是无限的通带纹波与阻带衰减可表示为

     通带纹波= $20lg(frac{A_o(f)}{A_i(f)})$ 

      其中$A_0(f)$、$A_i(f)$分别是一定频率f下的输出和输入振幅,通带波纹的单位为分贝(dB)。

      例如对于-0.02dB的波纹:$-0.02=20lg(frac{A_o(f)}{A_i(f)}) Rightarrow frac{A_o(f)}{A_i(f)}=10^{-0.001}=0.9977$。这表明,输出和输入的幅值之比近似为1。当阻带有-60dB的衰减,则$frac{A_o(f)}{A_i(f)}=10^{-3}=0.001$,表明输出的幅度是输入的1/1000。

    • IIRFIR滤波器

      滤波器也可以按冲击响应来分类对信号

    $$x(i)=egin{cases}1 & i=0 \0 & i eq 0end{cases}$$

      的响应为冲击响应,如下图所示。由傅里叶变换后的频率响应可知滤波器在不同频率下的增益理想情况下通带增益为1,阻带增益为0;所以频率在通带范围内的信号可完全通过在阻带的信号则不能通过

      如果冲击响应经一段时间后增益为零则滤波器为有限冲击响应滤波器(FIR)或递归数字滤波器,反之则为无限冲击响应滤波器(IIR)或非递归数字滤波器。其最基本的不同是FIR的输出只是由当前以及过去的输入值决定;而IIR的输出不但由当前以及过去的输入值决定还与过去的输出值有关。

      例如一个商店计算收款额,现在所卖商品的价格为$x(k)$,过去所卖商品的价格为$x(k-1)$,N为所有商品的总数总额为$y(k)$,则:

    $$y(k)=x(k)+x(k-1)+x(k-2)+...+x(1)$$

      也可以写为

    $$y(k)=y(k-1)+x(k)$$

      上述两个公式对收款额的描述本质上是一样的,只是前一个公式只涉及到输入,而后一个公式不但有输入,而且还有输出。前一公式是非递归形式,即FIR;而后一个公式则是递归形式,即IIR。如果加上6%的税,则上式可以写为:$y(k)=1.06{x(k)+x(k-1)+x(k-2)+...+x(1)}$或$y(k)=y(k-1)+1.06x(k)$ 。前一个式子中有一个公因子1.06,而后一个式子中y(k-1)的因子是1,x(k)的因子是1.06。这些因子都是滤波器的系数因子。IIR滤波器中与输入相乘的因子是前项因子与输出相乘的是递归因子。 

      对IIR和FIR滤波器来说,前者有较平坦的幅频特性,而且由于它的递归性,可以减少存储需求;而后者可以实现相位不失真。由于前者的设计方法源于传统的模拟滤波器,且使用滤波器时人们主要关心幅频特性,所以IIR滤波器使用较多。数字信号处理中常见的IIR滤波器有:巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器、贝塞尔滤波器等,各个滤波器具有不同的频率特性。如果有各种不同的滤波器及参数设计时就必须考虑选择最适合的滤波器影响滤波器的因素有是否需要线性相位是否允许纹波是否需要一个窄带的过渡带。在具体设计一种滤波器前需对各种滤波器做多次试验。

      巴特沃斯滤波器的幅频特性表达式为:

    $$|H(omega )|=frac{1}{sqrt{1+left(frac{omega }{omega _c} ight){}^{2n}}}$$

      式中$omega_c$为通带截止频率(在截止频率处有3dB衰减),$n$为滤波器阶数(阶数越高,过渡带越窄)。从下图可以看出巴特沃斯滤波器在通带内具有最大平坦的幅频特性,在截止频率后单调下降其频响特性是最平滑的通带中是理想的单位响应阻带中响应为零  

      巴特沃斯滤波器在通带与阻带之间过度缓慢,相比于巴特沃斯滤波器,切比雪夫滤波器在通带与阻带之间能够达到快速的过渡。如下图所示,契比雪夫滤波器的频响特点是在通带响应中有一个等幅的纹波阻带中单位衰减契比雪夫滤波器较巴特沃斯滤波器的优点是,用较少的阶数就能使过渡带很陡从而加快了滤波速度降低了绝对误差 。(分贝:表示功率量之比的一种单位,等于功率强度之比的常用对数的10倍;表示场量之比的一种单位,等于场强幅值之比的常用对数的20倍,即$1dB=10lgP=20lgV$。分贝刻度将宽的幅度范围压缩在一个较小的数集里,在声震动测量及显示频域信息时会经常用到

    如下图所示,相同阶数时:
    巴特沃斯滤波器通带最平坦,阻带下降慢;
    切比雪夫滤波器通带等纹波,阻带下降较快;
    椭圆滤波器在通带等纹波,阻带下降最快;
    贝塞尔滤波器通带等纹波,阻带下降慢(即幅频特性最差,但是贝塞尔滤波器具有最佳的线性相位特性)

      各个滤波器具有不同的频率特性,如果有各种不同的滤波器及参数设计时就必须考虑选择最适合的滤波器影响滤波器的因素有是否需要线性相位是否允许纹波是否需要一个窄带的过渡带。下图表示如何选择滤波器在具体设计一种滤波器前需对各种滤波器做多次试验。 

    • LabVIEW中的滤波器

      在进行信号分析和处理时,分析数据的一般过程是:初始化缓冲区、数据分析、数据输出,再根据缓冲区中的数据块进行分析。这种基于数据块的分析方法难以实现高速实时分析。LabVIEW提供了一类逐点分析节点,可以一个数据点接一个数据点分析,此时数据分析是基于数据点的,可以实现实时处理。使用逐点分析可以与信号同步,数据丢失的可能性更小,对硬件设备采样率的要求更低。逐点分析节点位于函数选板的“信号处理→逐点”

      关于使用逐点滤波器的范例见labviewexamplesSignal ProcessingPoint by Point中的PtByPt and Array Based Filter VI

    注意几点:

      (1)采样频率

      一般软件中数字滤波器模板中的频率都是归一的化频率,归一化的频率通过采样频率这一参数和实际频率对应起来。对于各种类型滤波,采样频率均应设置成滤波器输入信号的采样频率。

      (2)纹波幅度

      切比雪夫数字滤波器通带段幅频特性呈波纹状,需此参数控制纹波幅度,一般取0.1dB。巴特沃斯和贝塞尔滤波器通带段幅频特性曲线比较平坦,不需要此参数。

      (3)截止频率

      高截止频率f2和低截止频率f1fs为采样频率)必须符合条件:0 < f1 < f2 < 0.5fs

      (4) 滤波器阶数

      随着阶数的提高,与理想滤波器的特性越接近。但由于它的延迟也和滤波阶数有关,所以滤波器的性能不能单纯的靠增加滤波器的阶数来提高,因为增加滤波器阶数的同时也会使延迟变大。实际测试中应充分了解测试信号和噪声的特点,然后选取合适的滤波方案

      另外也可以使用Express模块中的滤波器,Express VI将基本函数面向应用做了进一步的打包,为用户提供了更加方便、简洁的编程途径。使用 Express VI,初学者无需面对复杂的连线即可快速入门。对于一般用户,Express VI 也可以起到减少连线,简化框图,突出程序主脉络等作用 。

  • 相关阅读:
    scjp考试准备
    scjp考试准备
    scjp考试准备
    scjp考试准备
    maven学习手记
    maven学习手记
    ExtJS MVC 学习手记3
    调整maven配置文件
    ExtJS MVC学习手记 2
    ExtJS MVC学习手记 1
  • 原文地址:https://www.cnblogs.com/21207-iHome/p/6064267.html
Copyright © 2011-2022 走看看