10种简单的数字滤波算法(C语言源程序)

假定从8位AD中读取数据（如果是更高位的AD可定义数据类型为int）,子程序为get_ad();1、限副滤波/*A值可根据实际情况调整value为有效值，new_value为当前采样值滤波程序返回有效的实际值*/#defineA10charvalue;charfilter(){charnew_value;new_value=get_ad();if((new_value-valueA)||(value-new_valueA)returnvalue;returnnew_value;}2、中位值滤波法/*N值可根据实际情况调整排序采用冒泡法*/#defineN11charfilter(){charvalue_buf[N];charcount,i,j,temp;for(count=0;countN;count++){value_buf[count]=get_ad();delay();}for(j=0;jN-1;j++){for(i=0;iN-j;i++){if(value_buf[i]value_buf[i+1]){temp=value_buf[i];value_buf[i]=value_buf[i+1];value_buf[i+1]=temp;}}}returnvalue_buf[(N-1)/2];}3、算术平均滤波法/**/#defineN12charfilter(){intsum=0;for(count=0;countN;count++){sum+=get_ad();delay();}return(char)(sum/N);}4、递推平均滤波法（又称滑动平均滤波法）/**/#defineN12charvalue_buf[N];chari=0;charfilter(){charcount;intsum=0;value_buf[i++]=get_ad();if(i==N)i=0;for(count=0;countN,count++)sum=value_buf[count];return(char)(sum/N);}5、中位值平均滤波法（又称防脉冲干扰平均滤波法）/**/#defineN12charfilter(){charcount,i,j;charvalue_buf[N];intsum=0;for(count=0;countN;count++){value_buf[count]=get_ad();delay();}for(j=0;jN-1;j++){for(i=0;iN-j;i++){if(value_buf[i]value_buf[i+1]){temp=value_buf[i];value_buf[i]=value_buf[i+1];value_buf[i+1]=temp;}}}for(count=1;countN-1;count++)sum+=value[count];return(char)(sum/(N-2));}6、限幅平均滤波法/**/略参考子程序1、37、一阶滞后滤波法/*为加快程序处理速度假定基数为100，a=0~100*/#definea50charvalue;charfilter(){charnew_value;new_value=get_ad();return(100-a)*value+a*new_value;}8、加权递推平均滤波法/*coe数组为加权系数表，存在程序存储区。*/#defineN12charcodecoe[N]={1,2,3,4,5,6,7,8,9,10,11,12};charcodesum_coe=1+2+3+4+5+6+7+8+9+10+11+12;charfilter(){charcount;charvalue_buf[N];intsum=0;for(count=0,countN;count++){value_buf[count]=get_ad();delay();}for(count=0,countN;count++)sum+=value_buf[count]*coe[count];return(char)(sum/sum_coe);}9、消抖滤波法#defineN12charfilter(){charcount=0;charnew_value;new_value=get_ad();while(value!=new_value);{count++;if(count=N)returnnew_value;delay();new_value=get_ad();}returnvalue;}10、限幅消抖滤波法/**/