智能仪器8.3.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,8.3,数字滤波技术,当存在随机干扰使被测信号中混入了噪声成分时，可以采用滤波器滤掉信号中的噪声成分、提高信号质量。,模拟滤波器在低频和甚低频时实现是比较困难的，而数字滤波器则不存在这些问题。由于它具有高精度、高可靠性和高稳定性的特点，因此在智能仪器中被广泛用于克服随机误差。,2,采用数字滤波算法克服随机误差具有如下优点：,数字滤波是由软件程序实现的，不需要硬件，因此不存在阻抗匹配问题。,对于多路信号输入通道，可以共用一个软件“滤波器”，从而降低仪表的硬件成本。,适当改变滤波器程序或运算参数，就能方便池改变滤波特性。,3,限幅滤波,生产经验告诉人们，很多物理量的变化是需要一定时间的，相邻两采样值之间的变化也有一个限度，这就是程序判断法的客观依据。,限幅滤波又称为程序判断滤波，由于测控系统存在随机脉冲干扰，或由于变送器不可靠而将尖脉冲干扰引入输入端，从而造成测量信号的失真。对于这种随机干扰，限幅滤波是一种十分有效的方法。其基本方法是比较相邻的两个,4,5,6,7,通过,A,送出,;,补码运算,取反加,1,得真值,8,中位值滤波,中位值滤波就是对某一被测参数连续采样,n,次,(,一般取,n,为奇数,),，然后把,n,次采样值按大小排队，取中间值为本次采样值。,中位值滤波能有效地克服因偶然因素引起的波动或仪表不稳定引起的误码所造成的脉冲干扰。,对温度、液位等缓慢变化的被测参数，采用这种方法能收到良好的效果；但对于流量、压力等快速变化的参数，一般不采用中位值滤波算法。,9,采样次数,N,越大,滤波效果越好,但采样次数,N,太大则会影响仪器的测试速度,所以,N,一般取,3,或,5,。对于变化很慢的参数,有时也可增加次数,例如,15,次。,中值滤波程序主要由数据排列和取中间值两部分组成。数据排列可采用几种常规的排序方法,如冒泡法、沉底法等。,中位值滤波子程序,设,SWAP,为存放采样值的内存单元首地址，,DATA,为存放滤波值的内存单元地址，,N,为采样值个数。,10,11,12,;A=SWAP+N/2,13,算术平均滤波,算术平均滤波法就是连续取,n,个采样值进行平均，其数学表达式为,算术平均滤波法用于对一般具有随机干扰的信号进行滤波。这种信号的特点是围绕着一个平均值，在某一范围附近作上下波动，因此，仅取一个采样值作为滤波值是不准确的。,算术平均滤波法对信号的平滑程度完全取决于,N,。理论上讲，在无系统误差的情况下，当,N,，其平均值趋于最大期望值，但实际上,N,是有限的。当,N,较大时，平滑度高，但灵敏度低；当,N,较小时，平滑度低，但灵敏度高。应根据具体情况选取,N,，既保证滤波效果，又尽量减少计算时间。,14,程序,8.13,算术平均滤波子程序,设,N,为采样值个数。,10,位的采样值存放在,SWAP,开始的单元中，采样值数据高,2,位放于高字节，低,8,位放于低字节，且假定,N,个采样值之和不超过,16,位。滤波结果放在,SWAP,和,SWAP+1,单元中。,DIV21,为双字节除以单字节子程序，,(R7,，,R6),为被除数，,(R5),为除数，商在,R7,，,R6,中。,15,16,算术平均滤波法每计算一次数据需测量,N,次，故这种方法不适于测量速度较快或要求数据计算速度较高的实时系统。,17,去极值平均滤波,算术平均滤波不能将明显的脉冲干扰消除，只是将其影响削弱。因明显干扰使采样值远离真实值，可以比较容易地将其剔除，不参加平均值计算，从而使平均滤波的输出更接近真值。,算法原理：连续采样,N,次，将其累加求和，同时找出其中的最大值与最小值，再从累加和中减去最大值和最小值，按,N-2,个采样值求平均，即得有效采样值为,18,为使平均滤波方便，,N-2,应为,2,、,4,、,8,、,16,。,具体方法有两种：对于快变参数，先连续采,N,次，然后再处理，但要在,RAM,中开辟,N,个数据暂存区；对于慢变参数，可一边采样，一边处理，而不必在,RAM,中开辟数据暂存区。,下面以,N,10,为例，介绍边采样边计算的程序设计方法。程序流程如下图所示。,19,20,去极值平均滤波子程序,设,R4,保存采样最大值，,R5,保存采样最小值，,R6,保存采样值，,R2,保存采样和的进位，,R3,保存采样值之和，按上图所示流程进行处理，结果保存在寄存器,A,中。,21,22,23,;,准备计算,(R2:R3)-R4,;,交换后,R2,保存和的低字节,;,R3,保存和的高字节,;,准备计算,(R3:R2)-R5,;,交换后,R3,保存和的低字节,;,R2,保存和的高字节,24,25,递推（滑动）平均滤波,该方法是先在,RAM,中建立一个数据缓冲区,依顺序存放,N,次采样数据（即把,N,个测量数据看成一个队列,队列的长度固定为,N,）。然后每采进一个新的数据,就将新数据存入队尾,同时将缓冲区中最早采集（队首）的一个数据去掉。,再求出当前,RAM,缓冲区中的,N,个数据的算术平均值。这样,每进行一次采样,就可计算出一个新的平均值,即测量数据取一丢一,测量一次便计算一次平均值,大大加快了数据处理能力。其数学表达式为,26,式中,:,y,n,第,n,次采样值经滤波后的输出；,y,n-i,未经滤波的第,n-i,次采样值；,N,递推平均项数。,通常采用环形队列结构实现数据的存放和平均值的计算，具体方法举例如下：,假设环形队列的地址为,40H,4FH,共,16,个单元,用,R,0,作队尾指示,并且,INPUT,为新采样数据处理子程序,子程序已将新数据置入累加器,A,中,其流程图如下图所示。,27,28,29,30,递推平均滤波法对周期性干扰有良好的抑制作用,平滑度高,灵敏度低；但对偶然出现的脉冲性干扰的抑制作用差,不易消除由于脉冲性干扰引起的采样值偏差。因此,它不适用于脉冲干扰比较严重的场合。通过观察不同,N,值下递推平均的输出响应来选取,N,值,既能少占用计算机时间,又能达到最好的滤波效果。,可以看出,递推平均滤波法与算术平均滤波法在数学处理上是完全相似的,只是,N,个数据的实际意义不同而已。,31,加权递推平均滤波,算术平均滤波和递推平均滤波方法的主要缺点是：为了提高对干扰的抑制效果,必须增大平均范围,N,；但增大,N,将会引起有用信号的失真,特别是会引起有用信号中高频分量丰富的峰值部位的失真。下图表示了递推平均滤波中峰值失真、噪声幅度与平均次数,N,的关系。,为协调三者的关系,可以采用加权平均滤波法。,32,所谓加权平均,是指参加平均运算的各采样值按不同的比例进行相加并求取平均值。加权系数一般先小后大,以突出后若干次采样值的效果,加强系统对参数变化趋势的辨识。,N,项加权平均滤波的算法为,式中,C,0,C,1,C,N-1,为常数,它们的选取有多种方法,但应满足,C,0,+,C,1,+,C,N-1,=1,。,33,一阶惯性滤波,假设滤波器的输入电压为,x(t,),，输出为,y(t,),，它们之间存在如下关系：,数字化时应用其采样值，即,y,n,=,y(nt,),，,x,n,=,x(nt,),如果采样间隔,t,足够小，则上式的离散值近似为：,即,34,令,则可化为,y,n,=ax,n,+by,n-1,可见，系数,a+b,=1,对于直流，,y,n,=,x,n,，即该滤波器的直流增益为,1,。若取采样间隔,t,足够小，则,at,/RC,，滤波器的截止频率为：,系数,a,越大，滤波器的截止频率越高。若取,t,=50s,，,a=1/16,，则截止频率为：,35,当采用模拟滤波器来抑制低频干扰时，要求滤波器有大的时间常数和高精度的,RC,网络。增大时间常数要求增大,R,值，其漏电流也随之增大，从而使,RC,网络的误差增大，降低了滤波效果。而采用数字滤波算法来实现动态的,RC,滤波，则能很好地克服上述模拟滤波器的缺点。在滤波常数要求较大的场合，这种方法更为适用。,一阶惯性滤波算法对于周期干扰具有良好的抑制作用，其不足之处是带来了相位滞后，灵敏度低。滞后程度取决于,b,值的大小。,同时它不能滤除频率高于采样频率,l,2(,奈奎斯特频率,),的干扰信号，例如采样频率为,100Hz,，则它不能滤除,50Hz,以上的干扰信号。对于高于奈奎斯持频率的干扰信号，应该采用模拟滤波器。,36,一阶惯性滤波子程序,设,y,n-1,在,DATA1,为首址的单元中，,x,n,在,DATA2,为首地址的单元中，均为双字节。取,a=0.25,，,b=0.75,，滤波结果在,R2,、,R3,中。,FOF,：,MOV R0,，,DATA1,MOV R1,，,DATA2,CLR C,;0.5Yn-1,存入,(R2:R3),INC R0,MOV A,，,R0,RRC A,MOV R2,，,A,DEC R0,MOV A,，,R0,RRC A,MOV R3,，,A,MOV A,，,R0,;Yn-1,+,Xn,存入,(A:R7),ADD A,，,R1,MOV R7,，,A,INC R0,37,INC R1,MOV A,，,R0,ADDC A,，,R1,CLR C,；,(Xn+Yn-1)0.5,存入,(R6:R7),RRC A,MOV R6,，,A,MOV A,，,R7,RRC A,MOV R7,，,A,CLR C,；,(Xn+Yn-1)0.25,存入,(R6:A),MOV A,，,R6,RRC A,MOV R6,，,A,MOV A,，,R7,RRC A,ADD A,，,R3,；,0.25(Xn+Yn-1)+0.5Yn-1,存于,(R2:R3),MOV R3,，,A,MOV A,，,R6,ADDC A,，,R2,MOV R2,，,A,RET,38,