《数字滤波》PPT课件.ppt

《《数字滤波》PPT课件.ppt》由会员分享，可在线阅读，更多相关《《数字滤波》PPT课件.ppt（37页珍藏版）》请在三一办公上搜索。

1、第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,一、数字滤波器基本概念,滤波器的作用,高速继电保护装置都工作在故障发生后的最初瞬变过程中，这时的电压和电流信号由于混有衰减直流分量和复杂的谐波成分而发生严重的畸变。目前大多数保护装置的原理是建立在反映正弦基波或某些整倍数谐波基础上，所以滤波器一直是继电保护装置的关键器件。目前所使用的微机继电保护几乎都采用了数字滤波器。,滤波器,模拟滤波器,数字式滤波器,第二章 微型计算机继电保护的数字信号基础,数字滤波器,2.2 数字滤波器,数字滤波器是将输入模拟信号x(t)经过采样和模数转换变成数字量后，进行某种数学运算（所谓数字滤波器通常是指一

2、种程序或算法（也可以用硬件实现）去掉信号中的无用成分，然后再经过数模转换得到模拟量输出y(t)。,（1）滤波精度高，具有高度的灵活性 通过改变滤波算法或某些参数，可以非常灵活地调整滤波特性，易于适应不同应用场合的要求。,（2）稳定性和可靠性高 模拟器件受环境影响较大，而数字系统受影响要小得多。,数字滤波器的优点：,（3）便于分时复用 采用模拟滤波器时，每一个输入通道都需要一个滤波器，而数字滤波器通过分时复用，一套数字滤波器就可以完成所有通道的滤波任务，保证了各通道的滤波特性完全一致。,按频率特性分类,带阻滤波器,高通滤波器,低通滤波器,带通滤波器,第二章 微型计算机继电保护的数字信号基础,按脉

3、冲响应时间的不同分类,2.2 数字滤波器,按运算结构不同分类,无限长冲击响应滤波器（IIR）,有限长脉冲响应滤波器（FIR）,IIR滤波器的脉冲响应h(n)包括无限个采样值，h(n)在 区间内有无限个不等于零的采样值。,FIR滤波器的脉冲响应h(n)，在 区间内有有限个采样点，n1和n2是两个有限值。,非递归型,递归型,数字滤波器的运算过程可用下述常系数线性差分方程来描述,通过选择滤波器系数ai和bj，可以滤除输入信号序列x(n)中的某些无用频率成分。,第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,当bi全为零时，上式称为非递归滤波器，非递归滤波器，即当前的输出y(n)只与过去

4、和当前的输入值x(n-i)有关，而与过去的输出值y(n-j)无关。,只含有前馈通道而无反馈通道,如果系数bj不全为0，即当前的输出不仅与过去和当前的输入有关，还与过去的输出有关，这就是递归型滤波器。,含有反馈通道,第二章 微型计算机继电保护的数字信号基础,二、简单滤波器及级联滤波器,2.2 数字滤波器,最简单的数字滤波器是利用加、减法运算构成的线性滤波器。这种滤波器是假定输入信号由稳态基波、稳态整次谐波和稳态直流所组成，不考虑暂态过程和其它高频成分。所以这种滤波比较粗糙，但运算量较小。一般用于速度较低的保护中。如过负荷保护、过电流保护和一些后备保护。,减法滤波器（差分滤波器）,差分方程为,y(

5、n)=x(n)-x(n-k),z变换,Y(z)=X(z)(1 z-k),转移函数,将 代入得到它的幅频特性为,通常令fs=N f1，N=1，2，即每基频周期内有N点采样。,第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,如果想消除m次谐波，确定参数k（滤波器的阶次）值，则应当使时,p=1，2，.,第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,例1 已知采样频率fs=600Hz，基波频率f1=50Hz。则,（1）要消除偶次谐波，取m=2p（p=0，1，2，.），代入上式得k=6，可消去直流成分及2、4、6等偶次谐波。滤波器的差分方程为,y(n)=x(n)x(n-6)

6、,（2）要消除三次谐波及其整数倍的谐波，则可取m=3p，得k=4。滤波器的差分方程为：y(n)=x(n)x(x-4)。,（3）要消除直流、工频及其所有整数倍的谐波分量，可取m=p，k=12。,在电力系统稳态时，该滤波器没有输出；在系统发生故障后的一个基波周期内，它只输出故障分量，所以可以用来使实现启动元件，选相元件及其它利用故障分量原理构成的保护。,取前6Ts时刻的采样值,第二章 微型计算机继电保护的数字信号基础,减法滤波器程序入,取当前采样值,滤波结果 存入RAM区,返回,减法程序流程,2.2 数字滤波器,第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,2加法滤波器,差分方程为

7、,y(n)=x(n)+x(n-k),z变换为,H(z)=1+z-k,幅频特性为,p=0，1，2，,如果已知k，可得滤除的谐波次数,第二章 微型计算机继电保护的数字信号基础,例2 已知采样频率fs=600Hz，基波频率f1=50Hz，则N=fs/f1=12，得,2.2 数字滤波器,要消除1、3、5等奇次谐波，可令m=2p+1，得k=6，这种滤波器也不输出基波分量。,要消除3、9、15等次谐波，可令m=3(2p+1)，得k=2。,y(n)=x(n)+x(n-6),y(n)=x(n)+x(n-2),第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,3积分滤波器,差分方程为,y(n)=x(

8、n)+x(n-1)+x(n-2)+x(n-k),z变换为,H(z)=1+z-1+z-2+z-k=,幅频特性为,p=0，1，2，,如果已知k，可得滤除的谐波次数,若p取0时，m也为0，但幅频特性不为0，所以不论fs、k取何值，都不能滤除直流分量（具有低通特性）。,第二章 微型计算机继电保护的数字信号基础,频率特性：,2.2 数字滤波器,例3 已知fs=600Hz，代入上式得,（1）要滤除偶次谐波，可令m=2p，得k=5。,y(n)=x(n)+x(n-1)+x(n-2)+x(n-5),滤波器的差分方程为,（2）若要滤除3的整数倍谐波，可令m=3p，得k=3，滤波器的差分方程为,y(n)=x(n)+

9、x(n-1)+x(n-2)+x(n-3),第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,4加减法滤波器,差分方程为,z变换为,（1）k为奇数时的幅频特性,p=0，1，2，,如果已知k，可得滤除的谐波次数,可见这种滤波器可以滤除N/(k+1)整数倍次谐波。当p=0，m=0时，幅频特性也为零，即无论k、fs取何值，都能滤除直流分量。它的幅频特性曲线见图,第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,例4：已知fs=600Hz，则,若要消除3的整数倍谐波，可令m=3p，可得 k=3。,若要消除偶次谐波，可令m=2p，得 k=5。,第二章 微型计算机继电保护的数字信号基

10、础,2.2 数字滤波器,（2）k为偶数时的幅频特性,p=0，1，2，,如果已知k，可得滤除的谐波次数,第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,例5：已知fs=600Hz，则,若令m=4p+2，则可滤除2、6、10等次谐波，代入上式得k=2。,5简单滤波单元的组合,一个由m个简单滤波单元组成的级联滤波器的转移函数可表示为,幅频特性,相频特性,时延特性,第二章 微型计算机继电保护的数字信号基础,例6 设采样频率为fs=600Hz，要求完全滤除直流分量和2、3、4、6次谐波分量。,2.2 数字滤波器,解：采用两个减法滤波（带通）单元和两个积分滤波（低通),单元级联：,H1(z)

11、=(1-z-2)，0、6次,H2(z)=(1-z-6)，0、2、4次,H3(z)=(1+z-1+z-2)，4次,H4(z)=(1+z-1+z-2+z-3)3、6次,减法滤波单元可滤除直流分量和2、4、6次谐波分量，两个积分单元可滤除3、4、6次谐波，其幅频特性曲线如图所示。,第二章 微型计算机继电保护的数字信号基础,例7 有一个50Hz基频带通滤波器，分别由下列滤波单元组成,2.2 数字滤波器,一个减法滤波器,H1(z)=1-z-6,两个积分滤波器,采样频率为fs=1200Hz。试分析该滤波器能完全滤除多少次谐波？,解：滤波器的转移函数为,H(z)=H1(z)H2(z)H3(z)=(1-z-6

12、)(1-z-8)(1-z-10)/(1-z-1)2,幅频特性,当fs=1200Hz时，fs=24f1，Ts=1/（24f1），上式可写成,第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,幅频特性曲线,当m=4p时，减法滤波器可滤除直流分量和4、8、12次谐波。,当m=3p时，第一个积分滤波器H2(z)可滤除第3、6、9次谐波。,当m=2.4p时，第二个积分滤波器H3(z)可滤除2.4、4.8、7.2、9.6次谐波。,第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,三、用零极点配置法设计数字滤波器,这是通过直接在Z平面上合理的设置零点和极点，以得到符合要求的频率响应特

13、性和时延特性。,（一）零、极点对系统频率响应的影响,N阶数字滤波器的转移传递函数的一般形式为,式中，ci、di分别为Z平面上的零点、极点。,第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,为零向量，,为极向量。,第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,幅频特性,相频特性,零、极点位置对系统频率响应的影响,极点的影响,第二章 微型计算机继电保护的数字信号基础,零点的影响,2.2 数字滤波器,例题：如图所示为离散控制系统框图。试分析其稳定性。其中采样周期T=1s，K=1，H(s)=1。,解：此系统的脉冲传递函数为,其中,第二章 微型计算机继电保护的数字信号基础,

14、2.2 数字滤波器,第二章 微型计算机继电保护的数字信号基础,2.2 数字滤波器,所以,闭环系统的特征方程为,当K=1时，,解得,2.2 数字滤波器,（二）用零、极点配置法设计数字滤波器,1零、极点配置法,因为数字滤波器的差分方程系数必须是实数，这就要求数字滤波器的脉冲传函H(z)的零、极点必须是以共轭的形式成对出现；从稳定性考虑，要求极点的幅值必须小于1，即要位于单位圆内。从提高滤波效果看，系统的零点设置在单位圆上比较好。,H(z)零极点配置及与滤波效果的关系：,1）如果在单位圆上设置一对共轭复零点，如，,则数字滤波器的频率特性将会在f=f1处出现零传输，即可以完全滤除输入信号中的f1频率分

15、量。,2）如果在单位圆内设置一对共轭复极点，如，,则数字滤波器的频率特性将会在f=f2处出现峰值，而且r（极点到圆心的距离，称为极径）越接近1（即极点越接近圆上），峰值越尖锐，如图所示。说明数字滤波器对频率为f2的信号分量具有提取作用。,2.2 数字滤波器,简单零、极点的频率响应,零点的配置：,如果要求数字滤波器滤除输入信号中的第k次谐波分量，则应该在Z平面上 处设置零点（其中，r0，Ts为采样周期，k=fi/f1），若z是复数，为使脉冲传函的系数为实数，该点应是一对共轭复数，因此脉冲传函中的相应因子应设置为：,2.2 数字滤波器,若要完全滤除k次谐波，可令r=1，即在单位圆上设置共轭零点。当

16、零点设置在正实轴上时，有k=0（欲滤除直流），对应的传递函数为,滤除直流,当零点位于负实轴上时，对应的传函为,极点的配置：,这里讨论的是单一滤除某一次谐波时，数字滤波器的脉冲传函H(z)中的零点配置方法。而它的极点配置方法与零点配置方法类似，如果要提取输入信号中的第k次谐波分量，则极点的配置应为,，,在H(z)中的因子是,为保证数字滤波器的稳定性，r的取值小于1。,2.2 数字滤波器,综上所述，零极点配置法设计数字滤波器，实质上是根据工程的需要，提出要滤除的或要提取的谐波频率kf1来选择数字滤波器的脉冲传函H(z)的零极点，根据来构造数字滤波器的H(z)中的因子表达式Hk(z)，然后再利用H(

17、z)与对应的时域差分方程来实现满足要求的数字滤波器算法。,2.2 数字滤波器,2零点和极点配置方法的应用,（1）全零点数字滤波器,在H(z)中只含零点，不含极点的滤波器。它属于一种非递归型数字滤波器。,设采样频率fs是输入信号基波频率f1的N倍，且N为偶数。现要滤除K=m以外的整数次谐波，为满足这一要求，数字滤波器的零点应选为Z平面单位圆上的一些点r=1，则全零点数字滤波器的脉冲传函应选为各零点因子脉冲传函之积（相当于级联），即,，,式中A为放大倍数。幅频特性为,2.2 数字滤波器,如：设N=12，现要求保留基波，完全滤除直流和其他整数次谐波分量，则数字滤波器的脉冲传函为,将上式完全展开成关于

18、z的多项式，再根据时域差分方程与H(z)之间的对应关系，即可求出滤波器的时域差分方程，即时域算法。下表给出了滤除直流分量直到12次谐波分量的脉冲传函H(z)中所对应因子的表达式。,2.2 数字滤波器,可见，有一些脉冲传函是相同的，它们是关于6次谐波滤除因子对称的。这是由于只要谐波频率k1f1与k2f2所对应的脉冲传函零点 与 是共轭复数关系，滤除k1次和k2次谐波的脉冲传函因子Hk1(z)和Hk2(z)就是相同的。同时也说明了同一个数字滤波器可以滤除两种不同频率的谐波信号。,在实际工程中，当需要同时滤除第k次和k+1次谐波信号时，可以将这两个谐波滤波器的脉冲传函相乘，就得到了所要设计的数字滤波

19、器的脉冲传函H(z)，即,H(z)=Hk(z)Hk+1(z),例 设采样频率fs=600Hz，要求保留基波，完全滤除直流分量和其他谐波分量（最高k=5次）。试用零点配置法设计数字滤波器。,解 根据全零点数字滤波器的设计方法得,H0(z)=1-z-1，H2(z)=1-z-1+z-2，H3(z)=1+z-2 H4(z)=1+z-1+z-2，H5(z)=1+1.732z-1+z-2,2.2 数字滤波器,其脉冲传函为,差分方程为,由于需要滤除的谐波次数较多，所以这个滤波器的数据窗较长（10个采样值）。,特点：在输入信号中，只含有直流分量、基波和有限个基波频率整数倍频的谐波分量时，采用这种方法提取基波分

20、量或某一个基波整数倍频的谐波分量是一个非常有实用价值的数字滤波算法。但当要滤除的谐波分量的个数较多时，滤波器的数据量用的较多，使滤波器的响应时间增加。,2.2 数字滤波器,（2）狭窄带通滤波器,由前面分析可知，全零点数字滤波器能够滤除输入信号中不需要的基波及各整数倍次谐波分量，但是在输入信号中一旦出现了无法预知的非整数倍频的谐波成分时，这种滤波器的输出信号中将含有无用的谐波成分。,为了抑制这种基波的非整数倍次谐波，需要在所提取的基频或基波某一个整数倍频处设置一个极点滤波器，以提高有用信号在输出信号中所占的比例。,狭窄带通滤波器是按照幅频特性来定义的。如果在数字滤波器所要提取k次谐波信号频率处设

21、置一个极点滤波器，此极点选择为,当r接近1时，就可以得到很窄的通带和比较陡峭的过渡带。,2.2 数字滤波器,1）极点必须是以共轭形式成对出现，且幅值r小于1，即极点位于单位圆内。设对基频的周期采样点数为N，即采样频率fs=Nf1，采样周期Ts=1/(Nf1)，则,假设在基波频率f=kf1，k=1处设置一个狭窄带通滤波器。,其对应的脉冲传函因子为,当要求完全滤除直流和N/2次谐波分量，提取基频分量时，滤波器的脉冲传函为,2.2 数字滤波器,其时域差分方程为,其幅频特性为,幅频特性曲线,r1r2r3,fs1fs2fs3,r的影响,fs的影响,2.2 数字滤波器,例如，在线路距离保护中，常常需要准确抽取故障信号中的基频分量，为了很好地抑制故障信号中其他非基频噪声分量的影响，可采用以基频频率为中心频率的带通滤波器。并要求该滤波器的通带带宽小、阻带衰耗大、过渡带陡峭。假设微机保护的采样频率为1000Hz（每周波20点)，采用零、极点配置法设计一递归型数字滤波器，其滤波方程为,滤波器的幅频特性曲线,从图中可以看出，采用递归型数字滤被器可以获得相当理想的滤波特性，并且计算简单，便于实时应用。而若采用非递归型滤波器，要想取得同样的滤波效果，滤波器的结构将十分复杂，数据窗也很长，运算量要大得多。,