数字滤波算法

        在仪表自动化工作中,经常需要对大量的数据进行处理,这些数据往往是一个时间序列或空间序列,这时常会用到数字滤波技术对数据进行预处理。数字滤波是指利用数学的方法对原始数据进行处理,去掉原始数据中掺杂的噪声数据,获得最具有代表性的数据集合。数据采样是一种通过间接方法取得事物状态的技术,如将事物的温度、压力、流量等属性通过一定的转换技术将其转换为电信号,然后再将电信号转换为数字化的数据。在多次转换中由于转换技术客观原因或主观原因造成采样数据中掺杂少量的噪声数据,影响了最终数据的准确性。为了防止噪声对数据结果的影响,除了采用更加科学的采样技术外,我们还要采用一些必要的技术手段对原始数据进行整理、统计,数字滤波技术是最基本的处理方法,它可以剔除数据中的噪声,提高数据的代表性。       

为消除或减弱外界干扰，提高系统的测量精度和可靠性，通过软件滤波的方法即数字滤波算法对测量数据进行处理。 数字滤波方法有很多种，每种方法有其不同的特点和使用范围。从大的范围可分为3类。  

1、克服大脉冲干扰的数字滤波法

（1）限幅滤波法

（2）中值滤波法

2、抑制小幅度高频噪声的平均滤波法

（1）算数平均

（2）滑动平均

（3）加权滑动平均

（4）一阶滞后滤波法

 3、复合滤波法

（1）中位值平均滤波法

（2）限幅平均滤波法

（一）克服大脉冲干扰的数字滤波法

    克服由仪器外部环境偶然因素引起的突变性扰动或仪器内部不稳定引起误码等造成的尖脉冲干扰，是仪器数据处理的第一步。通常采用简单的非线性滤波法。

1、限幅滤波法

        限幅滤波是通过程序判断被测信号的变化幅度，从而消除缓变信号中的尖脉冲干扰。

A 方法： 根据经验判断，确定两次采样允许的最大偏差值（设为δ）每次检测到新值时判断：如果本次值与上次值之差大于δ,则本次值无效,放弃本次值,用上次值代替本次值；

B 优点：能有效克服因偶然因素引起的脉冲干扰；

C 缺点：无法抑制那种周期性的干扰 平滑度差；

D 适用范围：变化比较缓慢的被测量值。

#define A 10
char value;
char filter(int value)
{
   char new_value;
   new_value = get_ad();
   if(( new_value - value > δ)||( value - new_value >δ)
      return value;
   return new_value;
}

2、中值滤波法

    中位值滤波是一种典型的非线性滤波器，它运算简单，在滤除脉冲噪声的同时可以很好地保护信号的细节信息。

A 方法：连续采样N次（N取奇数） 把N次采样值按大小排列（多采用冒泡法）取中间值为本次有效值；

B 优点：能有效克服因偶然因素引起的波动(脉冲)干扰；

C 缺点：对流量、速度等快速变化的参数不宜；

D 适用范围：对温度、液位的变化缓慢的被测参数有良好的滤波效果。

#define N 11

char filter()
{
  char value_buf[N];
  char count,i,j,temp;
  for(count=0;count<N;count++)
  {
    value_buf[count] = get_ad();
    delay();
	}
  for(j=0;j<N-1;j++)
  {
    for(i=0;i<N-j;i++)
    {
      if( value_buf>value_buf[i+1] )
      {
        temp = value_buf;
        value_buf = value_buf[i+1];
        value_buf[i+1] = temp;
      }
    }
 }
 return value_buf[(N-1)/2];
}

（二）抑制小幅度高频噪声的平均滤波法

   小幅度高频电子噪声：电子器件热噪声、A/D量化噪声等。通常采用具有低通特性的线性滤波器：算数平均滤波法、加权平均滤波法、滑动加权平均滤波法一阶滞后滤波法等。

 1、算术平均滤波法

    算术平均滤波法是对N个连续采样值相加，然后取其算术平均值作为本次测量的滤波值。

A 方法：连续取N个采样值进行算术平均运算N值较大时：信号平滑度较高，但灵敏度较低N值较小时：信号平滑度较低，但灵敏度较高 N值的选取：一般流量，N=12；压力：N=4；

B 优点：对滤除混杂在被测信号上的随机干扰信号非常有效。被测信号的特点是有一个平均值，信号在某一数值范围附近上下波动；

C 缺点：不易消除脉冲干扰引起的误差。对于采样速度较慢或要求数据更新率较高的实时系统，算术平均滤法无法使用的。

#define N 12
char filter()
{
	int sum = 0;
	for( count=0;count<N;count++)
  {
    sum + = get_ad();
    delay();
  }
return (char)(sum/N);
}

2、递推平均滤波法

    又称滑动平均滤波法，对于采样速度较慢或要求数据更新率较高的实时系统，算术平均滤波法无法使用的。滑动平均滤波法把N个测量数据看成一个队列，队列的长度固定为N，每进行一次新的采样，把测量结果放入队尾，而去掉原来队首的一个数据，这样在队列中始终有N个“最新”的数据。

A 方法：把连续取N个采样值看成一个队列，队列的长度固定N，每次采样到一个新数据放入队尾,并扔掉原来队首的一次数据.(先进 先出原则) 把队列中的N个数据进行算术平均运算,就可获得新的滤波结果N值的选取：流量，N=12；压力：N=4；液面，N=4~12；温度，N=1~4；

B 优点：对周期性干扰有良好的抑制作用，平滑度高 适用于高频振荡的系统；

C 缺点：灵敏度低，对偶然出现的脉冲性干扰的抑制作用较差不易消除由于脉冲干扰所引起的采样值偏差 不适用于脉冲干扰比较严重的场合。

#define N 12 
char value_buf[N];
char i=0;
char filter()
{
   char count;
   int  sum=0;
   value_buf[i++] = get_ad();
   if( i == N )   i = 0;
   for( count=0;count<N,count++)
      sum += value_buf[count];
   return (char)(sum/N);
}

3、加权递推平均滤波法

A 方法：是对递推平均滤波法的改进，即不同时刻的数据加以不同的权通常是，越接近现时刻的数据，权取得越大。给予新采样值的权系数越大，则灵敏度越高，但信号平滑度低；

B 优点：适用于有较大纯滞后时间常数的对象 和采样周期较短的系统；

C 缺点：对于纯滞后时间常数较小，采样周期较长，变化缓慢的信号 不能迅速反应系统当前所受干扰的严重程度，滤波效果差。

#define N 12
char code coe[N] = {1,2,3,4,5,6,7,8,9,10,11,12};
char code sum_coe = 1+2+3+4+5+6+7+8+9+10+11+12;
char filter()
{
   char count;
   char value_buf[N];
   int  sum=0;
   for (count=0,count<N;count++)
   {
      value_buf[count] = get_ad();
      delay();
   }
   for(count=0,count<N;count++)
      sum += value_buf[count] * coe[count];
   return (char)(sum/sum_coe);
}

4、一阶滞后滤波法

    一阶低通数字滤波器是用软件的方法实现硬件的RC滤波，以抑制干扰信号。在模拟量输入通道中，常用一阶滞后RC模拟滤波器来抑制干扰。用此种方法来实现对低频干扰时，首先遇到的问题是要求滤波器有大的时间常数（时间常数=RC）和高精度的RC网络。时间常数越大，要求RC值越大，其漏电流也必然增大，从而使RC网络精度下降。采用一阶滞后的数字滤波方法，能很好的克服这种模拟量滤波器的缺点，在滤波常数要求较大的场合，此法更适合。

A 方法： a = T f /( T f +T) T f 为滤波时间常数。T为采样周期，本次滤波结果=（1-a）* 本次采样值+a * 上次滤波结果；

B 优点： 对周期性干扰具有良好的抑制作用 适用于波动频率较高的场合；

C 缺点： 相位滞后，灵敏度低 滞后程度取决于a值大小 不能消除滤波频率高于采样频率的1/2的干扰信号。

#define a 50
char value;
char filter()
{
   char  new_value;
   new_value = get_ad();
   return (1-a)*value + a*new_value;  //a为滤波系数， 其值<<1
}

（三）复合滤波法

   在实际应用中，有时既要消除大幅度的脉冲干扰，有要做到数据平滑。因此常把前面介绍的两种以上的方法结合起来使用，形成复合滤波。去极值平均滤波算法：先用中值滤波算法滤除采样值中的脉冲性干扰，然后把剩余的各采样值进行平均滤波。连续采样N次，剔除其最大值和最小值，再求余下N－2个采样的平均值。显然，这种方法既能抑制随机干扰，又能滤除明显的脉冲干扰。

1、中位值平均滤波法

    又称防脉冲干扰平均滤法，中位值平均滤波法相当于“中位值滤波法”+“算术平均滤波法”

A 方法： 连续采样N个数据，去掉一个最大值和一个最小值 然后计算N-2个数据的算术平均值N值的选取：3~14；

B 优点： 融合了两种滤波法的优点这种方法既能抑制随机干扰，又能滤除明显的脉冲干扰；

C 缺点：测量速度较慢，和算术平均滤波法一样，数据量较大，比较浪费RAM。

2、限幅平均滤波法

    在脉冲干扰较严重的场合，如采用一般的平均值法，则干扰会平均到结果中去。限幅平均滤波法相当于“限幅滤波法”+“递推平均滤波法”

A 方法： 每次采样到的新数据先进行限幅处理，再送入队列进行递推平均滤波处理；

B 优点： 融合了两种滤波法的优点， 对于偶然出现的脉冲性干扰，可消除由于脉冲干扰所引起的采样值偏差；

C 缺点：数据量较大，比较浪费RAM；

D 适用范围： 缓变信号。