9个软件滤波的C程序

1-限副滤波:

1. /* A值可根据实际情况调整
   
2. value为有效值，new_value为当前采样值
   
3. 滤波程序返回有效的实际值 */
   
4. #define A 10
   
5. char value;
   
6. char filter()
   
7. {
   
8. char new_value;
   
9. new_value = get_ad();
   
10. if ( ( new_value - value > A ) || ( value - new_value > A )
    
11. return value;
    
12. return new_value;   
    
13. }

复制代码

2-中位值滤波法

1. 
   
2. /* N值可根据实际情况调整
   
3. 排序采用冒泡法*/
   
4. #define N 11
   
5. 
   
6. char filter()
   
7. {
   
8. char value_buf[N];
   
9. char count,i,j,temp;
   
10. for ( count=0;count<N;count++)
    
11. {
    
12. value_buf[count] = get_ad();
    
13. delay();
    
14. }
    
15. for (j=0;j<N-1;j++)
    
16. {
    
17. for (i=0;i<N-j;i++)
    
18. {
    
19.    if ( value_buf>value_buf[i+1] )
    
20.    {
    
21.       temp = value_buf;
    
22.       value_buf = value_buf[i+1];
    
23.          value_buf[i+1] = temp;
    
24.    }
    
25. }
    
26. }
    
27. return value_buf[(N-1)/2];
    
28. }

复制代码

3-算术平均滤波法

1. 
   
2. /*
   
3. */
   
4. 
   
5. #define N 12
   
6. 
   
7. char filter()
   
8. {
   
9. int sum = 0;
   
10. for ( count=0;count<N;count++)
    
11. {
    
12. sum + = get_ad();
    
13. delay();
    
14. }
    
15. return (char)(sum/N);
    
16. }

复制代码

4-递推平均滤波法（又称滑动平均滤波法）

1. 
   
2. /*
   
3. */
   
4. #define N 12
   
5. 
   
6. char value_buf[N];
   
7. char i=0;
   
8. 
   
9. char filter()
   
10. {
    
11. char count;
    
12. int sum=0;
    
13. value_buf[i++] = get_ad();
    
14. if ( i == N ) i = 0;
    
15. for ( count=0;count<N,count++)
    
16. sum = value_buf[count];
    
17. return (char)(sum/N);
    
18. }
    

复制代码

5-中位值平均滤波法（又称防脉冲干扰平均滤波法）

1. 
   
2. /*
   
3. */
   
4. #define N 12
   
5. 
   
6. char filter()
   
7. {
   
8. char count,i,j;
   
9. char value_buf[N];
   
10. int sum=0;
    
11. for (count=0;count<N;count++)
    
12. {
    
13. value_buf[count] = get_ad();
    
14. delay();
    
15. }
    
16. for (j=0;j<N-1;j++)
    
17. {
    
18. for (i=0;i<N-j;i++)
    
19. {
    
20.    if ( value_buf>value_buf[i+1] )
    
21.    {
    
22.       temp = value_buf;
    
23.       value_buf = value_buf[i+1];
    
24.          value_buf[i+1] = temp;
    
25.    }
    
26. }
    
27. }
    
28. for(count=1;count<N-1;count++)
    
29. sum += value[count];
    
30. return (char)(sum/(N-2));
    
31. }
    

复制代码

6-一阶滞后滤波法

1. 
   
2. /* 为加快程序处理速度假定基数为100，a=0~100 */
   
3. 
   
4. #define a 50
   
5. 
   
6. char value;
   
7. 
   
8. char filter()
   
9. {
   
10. char new_value;
    
11. new_value = get_ad();
    
12. return (100-a)*value + a*new_value;
    
13. }

复制代码

7-加权递推平均滤波法

1. /* coe数组为加权系数表，存在程序存储区。*/
   
2. 
   
3. #define N 12
   
4. 
   
5. char code coe[N] = {1,2,3,4,5,6,7,8,9,10,11,12};
   
6. char code sum_coe = 1+2+3+4+5+6+7+8+9+10+11+12;
   
7. 
   
8. char filter()
   
9. {
   
10. char count;
    
11. char value_buf[N];
    
12. int sum=0;
    
13. for (count=0,count<N;count++)
    
14. {
    
15. value_buf[count] = get_ad();
    
16. delay();
    
17. }
    
18. for (count=0,count<N;count++)
    
19. sum += value_buf[count]*coe[count];
    
20. return (char)(sum/sum_coe);
    
21. }

复制代码

8-消抖滤波法

1. #define N 12
   
2. 
   
3. char filter()
   
4. {
   
5. char count=0;
   
6. char new_value;
   
7. new_value = get_ad();
   
8. while (value ！=new_value);
   
9. {
   
10. count++;
    
11. if (count>=N) return new_value;
    
12.    delay();
    
13. new_value = get_ad();
    
14. }
    
15. return value;
    
16. }

复制代码

9-IIR滤波例子

1. 
   
2. int BandpassFilter4(int InputAD4)
   
3. {
   
4. int ReturnValue;
   
5. int ii;
   
6. RESLO=0;
   
7. RESHI=0;
   
8. MACS=*PdelIn;
   
9. OP2=1068; //FilterCoeff4[4];
   
10. MACS=*(PdelIn+1);
    
11. OP2=8; //FilterCoeff4[3];
    
12. MACS=*(PdelIn+2);
    
13. OP2=-2001;//FilterCoeff4[2];
    
14. MACS=*(PdelIn+3);
    
15. OP2=8; //FilterCoeff4[1];
    
16. MACS=InputAD4;
    
17. OP2=1068; //FilterCoeff4[0];
    
18. MACS=*PdelOu;
    
19. OP2=-7190;//FilterCoeff4[8];
    
20. MACS=*(PdelOu+1);
    
21. OP2=-1973; //FilterCoeff4[7];
    
22. MACS=*(PdelOu+2);
    
23. OP2=-19578;//FilterCoeff4[6];
    
24. MACS=*(PdelOu+3);
    
25. OP2=-3047; //FilterCoeff4[5];
    
26. *p=RESLO;
    
27. *(p+1)=RESHI;
    
28. mytestmul<<=2;
    
29. ReturnValue=*(p+1);
    
30. for (ii=0;ii<3;ii++)
    
31. {
    
32. DelayInput[ii]=DelayInput[ii+1];
    
33. DelayOutput[ii]=DelayOutput[ii+1];
    
34. }
    
35. DelayInput[3]=InputAD4;
    
36. DelayOutput[3]=ReturnValue;
    
37. 
    
38. // if (ReturnValue<0)
    
39. // {
    
40. // ReturnValue=-ReturnValue;
    
41. // }
    
42. return ReturnValue;
    
43. }

复制代码

然而这只是两个～

谢谢大神。谢谢大神。谢谢大神。谢谢大神。谢谢大神。谢谢大神。谢谢大神。谢谢大神。谢谢大神。谢谢大神。

这个数不够呀                              

感谢分享