一种新型的数据平滑滤波算法

该算法已经在 2个项目上运行大半年了，根据目前客户反馈 数据平滑效果还是很不错的

思路：将需要平滑的数据第一次的值放大N倍(最好是2的多少次方 需要保证放大后的数据不被溢出)，

加上最近的一个数据N 减去第N个数据中最早的那个数据，缩小到1/N，得到的值即为平滑之后的数据

直接以DS18B20温度滤波数据为例：

#include "ds18b20.h"
#include "delay.h"


short tmp_arg; //ζÈƽ»¬Â˲¨
//¸´Î»DS18B20
void DS18B20_Rst(void)   
{                 
DS18B20_IO_OUT(); //SET PG11 OUTPUT
    DS18B20_DQ_OUT=0; //À­µÍDQ
    delay_us(750);     //À­µÍ750us
    DS18B20_DQ_OUT=1; //DQ=1 
delay_us(15);     //15US
}
//µÈ´ýDS18B20µÄ»ØÓ¦
//·µ»Ø1:δ¼ì²âµ½DS18B20µÄ´æÔÚ
//·µ»Ø0:´æÔÚ
u8 DS18B20_Check(void)   
{   
u8 retry=0;
DS18B20_IO_IN();//SET PG11 INPUT  
    while (DS18B20_DQ_IN&&retry<200)
{
retry++;
delay_us(1);
};  
if(retry>=200)return 1;
else retry=0;
    while (!DS18B20_DQ_IN&&retry<240)
{
retry++;
delay_us(1);
};
if(retry>=240)return 1;   
return 0;
}
//´ÓDS18B20¶ÁÈ¡Ò»¸öλ
//·µ»ØÖµ£º1/0
u8 DS18B20_Read_Bit(void)  
{
    u8 data;
DS18B20_IO_OUT();//SET PG11 OUTPUT
    DS18B20_DQ_OUT=0; 
delay_us(2);
    DS18B20_DQ_OUT=1; 
DS18B20_IO_IN();//SET PG11 INPUT
delay_us(12);
if(DS18B20_DQ_IN)data=1;
    else data=0;  
    delay_us(50);           
    return data;
}
//´ÓDS18B20¶ÁÈ¡Ò»¸ö×Ö½Ú
//·µ»ØÖµ£º¶Áµ½µÄÊý¾Ý
u8 DS18B20_Read_Byte(void)     
{        
    u8 i,j,dat;
    dat=0;
for (i=1;i<=8;i++) 
{
        j=DS18B20_Read_Bit();
        dat=(j<<7)|(dat>>1);
    }    
    return dat;
}
//дһ¸ö×Ö½Úµ½DS18B20
//dat£ºÒªÐ´ÈëµÄ×Ö½Ú
void DS18B20_Write_Byte(u8 dat)     
 {             
    u8 j;
    u8 testb;
DS18B20_IO_OUT();//SET PG11 OUTPUT;
    for (j=1;j<=8;j++) 
{
        testb=dat&0x01;
        dat=dat>>1;
        if (testb) 
        {
            DS18B20_DQ_OUT=0; // Write 1
            delay_us(2);                            
            DS18B20_DQ_OUT=1;
            delay_us(60);             
        }
        else 
        {
            DS18B20_DQ_OUT=0; // Write 0
            delay_us(60);             
            DS18B20_DQ_OUT=1;
            delay_us(2);                          
        }
    }
}
//¿ªÊ¼Î¶Èת»»
void DS18B20_Start(void) 
{                 
    DS18B20_Rst();   
DS18B20_Check(); 
    DS18B20_Write_Byte(0xcc); // skip rom
    DS18B20_Write_Byte(0x44); // convert
} 


//³õʼ»¯DS18B20µÄIO¿Ú DQ ͬʱ¼ì²âDSµÄ´æÔÚ
//·µ»Ø1:²»´æÔÚ
//·µ»Ø0:´æÔÚ      
u8 DS18B20_Init(void)
{
  GPIO_InitTypeDef  GPIO_InitStructure;
 
  RCC_APB2PeriphClockCmd(RCC_APB2Periph_GPIOG, ENABLE);//ʹÄÜPORTG¿ÚʱÖÓ 

  GPIO_InitStructure.GPIO_Pin = GPIO_Pin_11;//PORTG.11 ÍÆÍìÊä³ö
  GPIO_InitStructure.GPIO_Mode = GPIO_Mode_Out_PP; 
  GPIO_InitStructure.GPIO_Speed = GPIO_Speed_50MHz;
  GPIO_Init(GPIOG, &GPIO_InitStructure);


  GPIO_SetBits(GPIOG,GPIO_Pin_11);    //Êä³ö1


DS18B20_Rst();


return DS18B20_Check();
}  
//´Óds18b20µÃµ½Î¶ÈÖµ
//¾«¶È£º0.1C
//·µ»ØÖµ£ºÎ¶ÈÖµ £¨-550~1250£© 
short DS18B20_Get_Temp(void)
{
    u8 temp;
    u8 TL,TH;
short tem;
    DS18B20_Start ();   // ds1820 start convert
    DS18B20_Rst();
    DS18B20_Check();  
    DS18B20_Write_Byte(0xcc); // skip rom
    DS18B20_Write_Byte(0xbe); // convert   
    TL=DS18B20_Read_Byte(); // LSB   
    TH=DS18B20_Read_Byte(); // MSB  
     
    if(TH>7)
    {
        TH=~TH;
        TL=~TL; 
        temp=0; //ζÈΪ¸º  
    }else temp=1; //ζÈΪÕý   
    tem=TH; //»ñµÃ¸ß°Ëλ
    tem<<=8;    
    tem+=TL; //»ñµÃµ×°Ëλ
    tem=(float)tem*0.625; //ת»»     
if(temp)return tem; //·µ»ØζÈÖµ
else return -tem;    
}

short Get_tempAverage(void)

{
u32 temp_sum;
short temp;

temp = DS18B20_Get_Temp();
if(DS18B20_Init() == 0)//温度采集起始条件
{
tmp_arg = temp;
temp_sum = tmp_arg;
temp_sum = temp_sum<<3;//放大8倍
}
temp_sum += temp;//加最近的一次温度值
temp_sum -= tmp_arg;//减去8次温度值中最早的那个值
tmp_arg = temp_sum >> 3;
return tmp_arg;
}

DS18B20.H部分

#ifndef __DS18B20_H
#define __DS18B20_H 
#include "sys.h"   


//IO·½ÏòÉèÖÃ
#define DS18B20_IO_IN()  {GPIOB->CRL&=0xFFFFF0FF;GPIOB->CRL|=8 << 1;}
#define DS18B20_IO_OUT() {GPIOB->CRL&=0xFFFFF0FF;GPIOB->CRL|=7 << 1;}
IO²Ù×÷º¯Êý   
#define DS18B20_DQ_OUT PBout(1) //Êý¾Ý¶Ë¿ÚPB1
#define DS18B20_DQ_IN  PBin(1)  //Êý¾Ý¶Ë¿ÚPB1 


extern short tmp_arg; //ζÈƽ»¬Â˲¨


u8 DS18B20_Init(void);//³õʼ»¯DS18B20
short DS18B20_Get_Temp(void);//»ñȡζÈ
void DS18B20_Start(void);//¿ªÊ¼Î¶Èת»»
void DS18B20_Write_Byte(u8 dat);//дÈëÒ»¸ö×Ö½Ú
u8 DS18B20_Read_Byte(void);//¶Á³öÒ»¸ö×Ö½Ú
u8 DS18B20_Read_Bit(void);//¶Á³öÒ»¸öλ
u8 DS18B20_Check(void);//¼ì²âÊÇ·ñ´æÔÚDS18B20
void DS18B20_Rst(void);//¸´Î»DS18B20  
short Get_tempAverage(void); //ζÈƽ»¬Â˲¨
#endif