单片机首要效果是操控外围的器材，并完成必定的通讯和数据处理。但在某些特定场合，不可避免地要用到数学运算，尽管单片机并不拿手完成算法和进行杂乱的运算。下面首要是介绍如何用单片机完成数字滤波。
在单片机进行数据收集时，会遇到数据的随机差错，随机差错是由随机搅扰引起的，其特点是在相同条件下丈量同一量时，其巨细和符号会现无规则的改动而无法猜测，但屡次丈量的成果契合核算规则。为战胜随机搅扰引起的差错，硬件上可选用滤波技能，软件上可选用软件算法完成数字滤波。滤波算法往往是体系测控算法的一个重要组成部分，实时性很强。
选用数字滤波算法战胜随机搅扰的差错具有以下长处：
1、数字滤波无需其他的硬件本钱，只用一个核算进程，可靠性高，不存在阻抗匹配问题。尤其是数字滤波能够对频率很低的信号进行滤波，这是模仿滤波器做不到的。
2、数字滤波运用软件算法完成，多输入通道可共用一个滤波程序，下降体系开支。
3、只需恰当改动滤波器的滤波程序或运算，就能便利地改动其滤波特性，这关于滤除低频搅扰和随机信号会有较大的效果。
4、在单片机体系中常用的滤波算法有限幅滤波法、中值滤波法、算术均匀滤波法、加权均匀滤波法、滑动均匀滤波等。
（1）限幅滤波算法
该运算的进程中将两次相邻的采样相减，求出其增量，然后将增量的绝对值，与两次采样答应的最大差值A进行比较。A的巨细由被测方针的详细情况而定，假如小于或等于答应的最大差值，则本次采样有用;不然取前次采样值作为本次数据的样本。
算法的程序代码如下：
#defineA //答应的最大差值chardata; //上一次的数据char filter（）{chardatanew; //新数据变量datanew=get_data（）; //取得新数据变量if（（datanew-data）》A||（data-datanew》A））return data;elsereturndatanew;}
阐明：限幅滤波法首要用于处理改动较为缓慢的数据，如温度、物体的方位等。运用时，要害要选取适宜的门约束A。一般这可由经历数据取得，必要时可通过试验得到。
（2）中值滤波算法
该运算的进程是对某一参数接连采样N次（N一般为奇数），然后把N次采样的值按从小到大摆放，再取中心值作为本次采样值，整个进程实际上是一个序列排序的进程。
算法的程序代码如下：
#define N11 //界说取得的数据个数char filter（）{charvalue_buff［N］; //界说存储数据的数组char count，i，j，temp;for（count=0;count{value_buf［count］=get_data（）;delay（）; //假如收集数据比较慢，那么就需要延时或中止}for（j=0;j{for（value_buff［i］》value_buff［i+1］{temp=value_buff［i］;value_buff［i］=value_buff［i+1］;value_buff［i+1］=temp;}}returnvalue_buff［（N-1）/2］;}
阐明：中值滤波比较适用于去掉由偶尔要素引起的动摇和采样器不稳定而引起的脉动搅扰。若被丈量值改动比较慢，选用中值滤波法效果会比较好，但假如数据改动比较快，则不宜选用此办法。
（3）算术均匀滤波算法
该算法的基本原理很简单，便是接连取N次采样值后进行算术均匀。
算法的程序代码如下：
char filter（）{int sum=0;for（count=0;count{sum+=get_data（）;delay（）：}return （char）（sum/N）;}
阐明：算术均匀滤波算法适用于对具有随机搅扰的信号进行滤波。这种信号的特点是有一个均匀值，信号在某一数值邻近上下动摇。信号的均匀滑润程度完全到决于N值。当N较大时，滑润度高，灵敏度低;当N较小时，滑润度低，但灵敏度高。为了便利求均匀值，N一般取4、8、16、32之类的2的整数幂，以便在程序顶用移位操作来替代除法。
（4）加权均匀滤波算法
由于前面所说的“算术均匀滤波算法”存在滑润度和灵敏度之间的对立。为了和谐滑润度和灵敏度之间的联系，可选用加权均匀滤波。它的原理是对接连N次采样值别离乘上不同的加权系数之后再求累加，加权系数一般先小后大，以杰出后边若干采样的效果，加强体系对参数改动趋势的知道。各个加权系数均小于1的小数，且满意总和等于1的完毕条件。这样加权运算之后的累加和即为有用采样值。设D为N个采样值的加权均匀值：XN-i为第N-i次采样值;N为采样次数;Ci为加权系数。加权系数Ci表现了各种采样值在均匀值中所占的份额。一般来说采样次数越靠后，取的份额越大，这样可增加新采样在均匀值中所占的比重。加权均匀值滤波法可杰出一部分信号抵抗另一部分信号，以进步采样值改动的灵敏度。
样例程序代码如下：
char codejq［N］={1，2，3，4，5，6，7，8，9，10，11，12}; //code数组为加权系数表，存在程序存储区char codesum_jq=1+2+3+4+5+6+7+8+9+10+11+12;char filter（）{char count;char value_buff［N］;int sum=0;for（count=0;count{value_buff［count］=get_data（）;delay（）;}for（count=0;countsum+=value_buff［count］*jq［count］;return（char）（sum/sum_jq）;}
（5）滑动均匀滤波算法
以上介绍和各种均匀滤波算法有一个共同点，即每获取一个有用采样值有必要接连进行若干次采样，当采速度慢时，体系的实时得不到确保。这儿介绍的滑动均匀滤波算法只采样一次，将一次采样值和曩昔的若干次采样值一同求均匀，得到的有用采样值即可投入运用。假如取N个采样值求均匀，存储区中有必要拓荒N个数据的暂存区。每新收集一个数据便存入暂存区中，一起去掉一个最老数据，保存这N个数据始终是最新更新的数据。选用环型行列结构能够便利地完成这种数据寄存方法。
程序代码如下：
char value_buff［N］;char i=0;char filter（）{char count;int sum=0;value_buff［i++］=get_data（）;if（i==N）i=0;for（count=0;countsum=value_buff［count］;return （char）（sum/N）;}
（6）低通滤波
将一般硬件RC低通滤波器的微分方程用差分方程来表求，变能够选用软件算法来模仿硬件滤波的功用，经推导，低通滤波算法如下：
Yn=a* Xn+（1-a） *Yn-1
式中 Xn——本次采样值
Yn-1——前次的滤波输出值;
，a——滤波系数，其值一般远小于1;
Yn——本次滤波的输出值。
由上式能够看出，本次滤波的输出值首要取决于前次滤波的输出值（留意不是前次的采样值，这和加权均匀滤波是有本质区别的），本次采样值对滤波输出的奉献是比较小的，但多少有些批改效果，这种算法便模仿了详细有教大惯性的低通滤波器功用。滤波算法的截止频率可用以下式核算：
fL=a/2Pit pi为圆周率3.14…
式中 a——滤波系数;
， t——采样间隔时间;
例如：当t=0.5s（即每秒2次），a=1/32时;
fL=（1/32）/（2*3.14*0.5）=0.01Hz
当方针参数为改动很慢的物理量时，这是很有用的。别的一方面，它不能滤除高于1/2采样频率的干搅信号，本例中采样频率为2Hz，故对1Hz以上的干搅信号应选用其他方法滤除，低通滤波算法程序于加权均匀滤波类似，但加权系数只要两个：a和1-a。为核算便利，a取一整数，1-a用256-a，来替代，核算成果舍去最低字节即可，由于只要两项，a和1-a，均以当即数的方式编入程序中，不别的设表格。尽管采样值为单元字节（8位A/D）。为确保运算精度，滤波输出值用双字节表明，其间一个字节整数，一字节小数，不然有或许由于每次舍去尾数而使输出不会改动。
设Yn-1寄存在30H（整数）和31H（小数）两单元中，Yn寄存在32H（整数）和33H（小数）中。
小结一下吧：数字滤波器，说白了，便是屡次采样求均匀值的一个进程，准确一点的，便是再顺序摆放，去掉首位再求均匀值，便是求均匀数！！
责任编辑:pj
