1、精選優(yōu)質(zhì)文檔-傾情為你奉上1、限幅濾波法（又稱程序判斷濾波法）A、方法：根據(jù)經(jīng)驗判斷，確定兩次采樣允許的最大偏差值（設(shè)為A）每次檢測到新值時判斷：如果本次值與上次值之差<=A,則本次值有效如果本次值與上次值之差>A,則本次值無效,放棄本次值,用上次值代替本次值B、優(yōu)點：能有效克服因偶然因素引起的脈沖干擾C、缺點無法抑制那種周期性的干擾平滑度差#define A 10 char value;char filter() char new_value; new_value = get_ad(); if ( ( new_value - value > A ) | ( value -

2、new_value > A ) return value; return new_value; 2、中位值濾波法A、方法：連續(xù)采樣N次（N取奇數(shù)）把N次采樣值按大小排列取中間值為本次有效值B、優(yōu)點：能有效克服因偶然因素引起的波動干擾對溫度、液位的變化緩慢的被測參數(shù)有良好的濾波效果C、缺點：對流量、速度等快速變化的參數(shù)不宜#define N 11 char filter() char value_bufN; char count,i,j,temp; for ( count=0;count<N;count+) value_bufcount = get_ad(); delay(); fo

3、r (j=0;j<N-1;j+) for (i=0;i<N-j;i+) if ( value_buf>value_bufi+1 ) temp = value_buf; value_buf = value_bufi+1; value_bufi+1 = temp; return value_buf(N-1)/2; 3、算術(shù)平均濾波法A、方法：連續(xù)取N個采樣值進行算術(shù)平均運算N值較大時：信號平滑度較高，但靈敏度較低N值較小時：信號平滑度較低，但靈敏度較高N值的選取：一般流量，N=12；壓力：N=4B、優(yōu)點：適用于對一般具有隨機干擾的信號進行濾波這樣信號的特點是有一個平均值，信號在某

4、一數(shù)值范圍附近上下波動C、缺點：對于測量速度較慢或要求數(shù)據(jù)計算速度較快的實時控制不適用比較浪費RAM#define N 12 char filter() int sum = 0; for ( count=0;count<N;count+) sum + = get_ad(); delay(); return (char)(sum/N); 4、遞推平均濾波法（又稱滑動平均濾波法）A、方法：把連續(xù)取N個采樣值看成一個隊列隊列的長度固定為N每次采樣到一個新數(shù)據(jù)放入隊尾,并扔掉原來隊首的一次數(shù)據(jù).(先進先出原則)把隊列中的N個數(shù)據(jù)進行算術(shù)平均運算,就可獲得新的濾波結(jié)果N值的選取：流量，N=12；壓

5、力：N=4；液面，N=412；溫度，N=14B、優(yōu)點：對周期性干擾有良好的抑制作用，平滑度高適用于高頻振蕩的系統(tǒng)C、缺點：靈敏度低對偶然出現(xiàn)的脈沖性干擾的抑制作用較差不易消除由于脈沖干擾所引起的采樣值偏差不適用于脈沖干擾比較嚴重的場合比較浪費RAM#define N 12 char value_bufN;char i=0; char filter() char count; int sum=0; value_bufi+ = get_ad(); if ( i = N ) i = 0; for ( count=0;count<N,count+) sum = value_bufcount; r

6、eturn (char)(sum/N); 5、中位值平均濾波法（又稱防脈沖干擾平均濾波法）A、方法：相當(dāng)于“中位值濾波法”+“算術(shù)平均濾波法”連續(xù)采樣N個數(shù)據(jù)，去掉一個最大值和一個最小值然后計算N-2個數(shù)據(jù)的算術(shù)平均值N值的選取：314B、優(yōu)點：融合了兩種濾波法的優(yōu)點對于偶然出現(xiàn)的脈沖性干擾，可消除由于脈沖干擾所引起的采樣值偏差C、缺點：測量速度較慢，和算術(shù)平均濾波法一樣比較浪費RAM#define N 12 char filter() char count,i,j; char value_bufN; int sum=0; for (count=0;count<N;count+) val

7、ue_bufcount = get_ad(); delay(); for (j=0;j<N-1;j+) for (i=0;i<N-j;i+) if ( value_buf>value_bufi+1 ) temp = value_buf; value_buf = value_bufi+1; value_bufi+1 = temp; for(count=1;count<N-1;count+) sum += valuecount; return (char)(sum/(N-2); 6、限幅平均濾波法A、方法：相當(dāng)于“限幅濾波法”+“遞推平均濾波法”每次采樣到的新數(shù)據(jù)先進行限幅

8、處理，再送入隊列進行遞推平均濾波處理B、優(yōu)點：融合了兩種濾波法的優(yōu)點對于偶然出現(xiàn)的脈沖性干擾，可消除由于脈沖干擾所引起的采樣值偏差C、缺點：比較浪費RAM略 參考子程序1、37、一階滯后濾波法(低通濾波)A、方法：取a=01本次濾波結(jié)果=（1-a）*本次采樣值+a*上次濾波結(jié)果B、優(yōu)點：對周期性干擾具有良好的抑制作用適用于波動頻率較高的場合C、缺點：相位滯后，靈敏度低滯后程度取決于a值大小不能消除濾波頻率高于采樣頻率的1/2的干擾信號#define a 50 char value; char filter() char new_value; new_value = get_ad(); retu

9、rn (100-a)*value + a*new_value; 8、加權(quán)遞推平均濾波法A、方法：是對遞推平均濾波法的改進，即不同時刻的數(shù)據(jù)加以不同的權(quán)通常是，越接近現(xiàn)時刻的數(shù)據(jù)，權(quán)取得越大。給予新采樣值的權(quán)系數(shù)越大，則靈敏度越高，但信號平滑度越低B、優(yōu)點：適用于有較大純滯后時間常數(shù)的對象和采樣周期較短的系統(tǒng)C、缺點：對于純滯后時間常數(shù)較小，采樣周期較長，變化緩慢的信號不能迅速反應(yīng)系統(tǒng)當(dāng)前所受干擾的嚴重程度，濾波效果差#define N 12 char code coeN = 1,2,3,4,5,6,7,8,9,10,11,12;char code sum_coe = 1+2+3+4+5+6+7

10、+8+9+10+11+12; char filter() char count; char value_bufN; int sum=0; for (count=0,count<N;count+) value_bufcount = get_ad(); delay(); for (count=0,count<N;count+) sum += value_bufcount*coecount; return (char)(sum/sum_coe); 9、消抖濾波法A、方法：設(shè)置一個濾波計數(shù)器將每次采樣值與當(dāng)前有效值比較：如果采樣值當(dāng)前有效值，則計數(shù)器清零如果采樣值<>當(dāng)前有效值

11、，則計數(shù)器+1，并判斷計數(shù)器是否>=上限N(溢出)如果計數(shù)器溢出,則將本次值替換當(dāng)前有效值,并清計數(shù)器B、優(yōu)點：對于變化緩慢的被測參數(shù)有較好的濾波效果,可避免在臨界值附近控制器的反復(fù)開/關(guān)跳動或顯示器上數(shù)值抖動C、缺點：對于快速變化的參數(shù)不宜如果在計數(shù)器溢出的那一次采樣到的值恰好是干擾值,則會將干擾值當(dāng)作有效值導(dǎo)入系統(tǒng)#define N 12 char filter() char count=0; char new_value; new_value = get_ad(); while (value !=new_value); count+; if (count>=N) return

12、 new_value; delay(); new_value = get_ad(); return value; 10、限幅消抖濾波法A、方法：相當(dāng)于“限幅濾波法”+“消抖濾波法”先限幅,后消抖B、優(yōu)點：繼承了“限幅”和“消抖”的優(yōu)點改進了“消抖濾波法”中的某些缺陷,避免將干擾值導(dǎo)入系統(tǒng)C、缺點：對于快速變化的參數(shù)不宜第11種方法：IIR 數(shù)字濾波器A. 方法：確定信號帶寬， 濾之。Y(n) = a1*Y(n-1) + a2*Y(n-2) + . + ak*Y(n-k) + b0*X(n) + b1*X(n-1) + b2*X(n-2) + . + bk*X(n-k)B. 優(yōu)點：高通，低通，帶

13、通，帶阻任意。設(shè)計簡單(用matlab）C. 缺點：運算量大。int BandpassFilter4(int InputAD4) int ReturnValue; int ii; RESLO=0; RESHI=0; MACS=*PdelIn; OP2=1068; /FilterCoeff44; MACS=*(PdelIn+1); OP2=8; /FilterCoeff43; MACS=*(PdelIn+2); OP2=-2001;/FilterCoeff42; MACS=*(PdelIn+3); OP2=8; /FilterCoeff41; MACS=InputAD4; OP2=1068; /

14、FilterCoeff40; MACS=*PdelOu; OP2=-7190;/FilterCoeff48; MACS=*(PdelOu+1); OP2=-1973; /FilterCoeff47; MACS=*(PdelOu+2); OP2=-19578;/FilterCoeff46; MACS=*(PdelOu+3); OP2=-3047; /FilterCoeff45; *p=RESLO; *(p+1)=RESHI; mytestmul<<=2; ReturnValue=*(p+1); for (ii=0;ii<3;ii+) DelayInputii=DelayInpu

15、tii+1; DelayOutputii=DelayOutputii+1; DelayInput3=InputAD4; DelayOutput3=ReturnValue; / if (ReturnValue<0) / / ReturnValue=-ReturnValue; / return ReturnValue; 第12種方法：RC濾波的一種實現(xiàn).RcDigital(double & X, double & Y) static int MidFlag; static double Yn_1,Xn_1; double MyGetX=0,MyGetY=0; double A

16、lfa; Alfa=0.7; if(X=0|Y=0) MidFlag=0; Xn_1=0; Yn_1=0; MyGetX=0; MyGetY=0; if(X>0&&Y>0) if(MidFlag=1) MyGetY = (1 - Alfa) * Y + Alfa * Yn_1; MyGetX = (1 - Alfa) * X + Alfa * Xn_1; Xn_1 = MyGetX; Yn_1 = MyGetY; else MidFlag=1; MyGetX = X; MyGetY = Y; Xn_1 = X; Yn_1 = Y; X = MyGetX; Y = M

17、yGetY; 實際上的應(yīng)用并沒有單一的應(yīng)用，例如：實際的AD上存在一定的沖擊干擾，這時就中位值平均濾波，又由于效率的問題對其添加滑動平均濾波法，提高其動態(tài)性。sint32 slAD =0; void filter() static sint32 slaTemp10 = 0; static uint8 k =0 ; static sint32 slFirst = 0; static uint8 ucFirstFlag =0; uint8 queue10=0,1,2,3,4,5,6,7,8,9; uint8 ucTemp =0; uint8 i,j =0; sint32 slTemp =0; slaTempk = get_ad()/2; k+; for(j=0;j<10;j+) for(i=j+1;i&