C51 COMPILER V7.06   AUTO                                                                  07/20/2008 16:34:55 PAGE 1   


C51 COMPILER V7.06, COMPILATION OF MODULE AUTO
OBJECT MODULE PLACED IN E:\溫度表\AUTO.OBJ
COMPILER INVOKED BY: D:\Keil\C51\BIN\C51.EXE E:\溫度表\AUTO.C DB SB OE

stmt level    source

   1          #include <MATH.H>
   2          #include <intrins.h>
   3          #include <float.h>
   4          extern unsigned int AD_Conver(unsigned char channel);
   5          
   6          /*void RMS_AND_AUTO()
   7          {
   8                  long float temp;
   9                  temp=filter();
  10          
  11          
  12          
  13          }
  14          */
  15          /***************************************************************************************************
  16          函數(shù)名稱：float CalculateRMS(const int x[], unsigned char N)
  17          函數(shù)入口：
  18          函數(shù)出口：
  19          功能描述：使用均方根的方法計算有效值。
  20          ****************************************************************************************************/
  21          /*float CalculateRMS(const int x[], unsigned char N)
  22          {
  23              unsigned long int   RMS = 0;          //有效值變量
  24              unsigned char   k   = 0;
  25          
  26              while (k < N)
  27              {
  28                  RMS  += (unsigned long int)x[k] * (unsigned long int)x[k];
  29                  ++k;
  30              }
  31          
  32              RMS  /= N;
  33          
  34              return sqrt(RMS);
  35          }*/
  36          /*十種軟件濾波方法(轉載)2007年06月12日 星期二 22:20十種軟件濾波方法*/
  37          /*1、限幅濾波法（又稱程序判斷濾波法）
  38               A、方法：
  39                   根據(jù)經驗判斷，確定兩次采樣允許的最大偏差值（設為A）
  40                   每次檢測到新值時判斷：
  41                   如果本次值與上次值之差<=A,則本次值有效
  42                   如果本次值與上次值之差>A,則本次值無效,放棄本次值,用上次值代替本次值
  43               B、優(yōu)點：能有效克服因偶然因素引起的脈沖干擾
  44               C、缺點
  45                   無法抑制那種周期性的干擾
  46                   平滑度差
  47                   1、限副濾波
  48             A值可根據(jù)實際情況調整
  49               value為有效值，new_value為當前采樣值
  50               濾波程序返回有效的實際值
  51          #define A 10
  52          
  53          char value;
  54          
  55          char filter()
C51 COMPILER V7.06   AUTO                                                                  07/20/2008 16:34:55 PAGE 2   

  56          {
  57              char   new_value;
  58              new_value = get_ad();
  59              if ( ( new_value - value > A ) || ( value - new_value > A )
  60                 return value;
  61              return new_value;
  62          
  63          }  */
  64          /* 2、中位值濾波法
  65               A、方法：
  66                   連續(xù)采樣N次（N取奇數(shù)）
  67                   把N次采樣值按大小排列
  68                   取中間值為本次有效值
  69               B、優(yōu)點：
  70                   能有效克服因偶然因素引起的波動干擾
  71                   對溫度、液位的變化緩慢的被測參數(shù)有良好的濾波效果
  72               C、缺點：對流量、速度等快速變化的參數(shù)不宜
  73               2、中位值濾波法
  74            N值可根據(jù)實際情況調整
  75               排序采用冒泡法
  76          #define N   11
  77          
  78          char filter()
  79          {
  80              char value_buf[N];
  81              char count,i,j,temp;
  82              for ( count=0;count<N;count++)
  83              {
  84                 value_buf[count] = get_ad();
  85                 delay();
  86              }
  87              for (j=0;j<N-1;j++)
  88              {
  89                 for (i=0;i<N-j;i++)
  90                 {
  91                    if ( value_buf[i]>value_buf[i+1] )
  92                    {
  93                       temp = value_buf[i];
  94                       value_buf[i] = value_buf[i+1];
  95                        value_buf[i+1] = temp;
  96                    }
  97                 }
  98              }
  99              return value_buf[(N-1)/2];
 100          } */
 101          /*3、算術平均濾波法
 102               A、方法：
 103                   連續(xù)取N個采樣值進行算術平均運算
 104                   N值較大時：信號平滑度較高，但靈敏度較低
 105                   N值較小時：信號平滑度較低，但靈敏度較高
 106                   N值的選取：一般流量，N=12；壓力：N=4
 107               B、優(yōu)點：
 108                   適用于對一般具有隨機干擾的信號進行濾波
 109                   這樣信號的特點是有一個平均值，信號在某一數(shù)值范圍附近上下波動
 110               C、缺點：
 111                   對于測量速度較慢或要求數(shù)據(jù)計算速度較快的實時控制不適用
 112                   比較浪費RAM
 113                   3、算術平均濾波法
 114          
 115          #define N 12
 116          
 117          char filter()
C51 COMPILER V7.06   AUTO                                                                  07/20/2008 16:34:55 PAGE 3   

 118          {
 119              int   sum = 0;
 120              for ( count=0;count<N;count++)
 121              {
 122                 sum + = get_ad();
 123                 delay();
 124              }
 125              return (char)(sum/N);
 126          }*/
 127          
 128          
 129           /*4、遞推平均濾波法（又稱滑動平均濾波法）
 130               A、方法：
 131                   把連續(xù)取N個采樣值看成一個隊列
 132                   隊列的長度固定為N
 133                   每次采樣到一個新數(shù)據(jù)放入隊尾,并扔掉原來隊首的一次數(shù)據(jù).(先進先出原則)
 134                   把隊列中的N個數(shù)據(jù)進行算術平均運算,就可獲得新的濾波結果
 135                   N值的選取：流量，N=12；壓力：N=4；液面，N=4~12；溫度，N=1~4
 136               B、優(yōu)點：
 137                   對周期性干擾有良好的抑制作用，平滑度高
 138                   適用于高頻振蕩的系統(tǒng)
 139               C、缺點：
 140                   靈敏度低
 141                   對偶然出現(xiàn)的脈沖性干擾的抑制作用較差
 142                   不易消除由于脈沖干擾所引起的采樣值偏差
 143                   不適用于脈沖干擾比較嚴重的場合
 144                   比較浪費RAM
 145                   4、遞推平均濾波法（又稱滑動平均濾波法）*/
 146          #define N 4
 147          unsigned  int value_buf[N];
 148          unsigned char i=0;
 149          long   float   filter(unsigned char port)
 150          {
 151   1      
 152   1      
 153   1          long float   sum=0;
 154   1          unsigned char count;
 155   1          value_buf[i++] =AD_Conver(port);
 156   1          if ( i == N )    i = 0;
 157   1          for ( count=0;count<N;count++)
 158   1             sum += value_buf[count];
 159   1             sum=sum*0.25;
 160   1             //sum*=10000;
 161   1             //sum=sum*9.999;
 162   1              sum=sum*5000/4095;
 163   1          return (long float)sum;
 164   1      }
 165          
 166          /*5、中位值平均濾波法（又稱防脈沖干擾平均濾波法）
 167               A、方法：
 168                   相當于“中位值濾波法”+“算術平均濾波法”
 169                   連續(xù)采樣N個數(shù)據(jù)，去掉一個最大值和一個最小值
 170                   然后計算N-2個數(shù)據(jù)的算術平均值
 171                   N值的選取：3~14
 172               B、優(yōu)點：
 173                   融合了兩種濾波法的優(yōu)點
 174                   對于偶然出現(xiàn)的脈沖性干擾，可消除由于脈沖干擾所引起的采樣值偏差
 175               C、缺點：
 176                   測量速度較慢，和算術平均濾波法一樣
 177                   比較浪費RAM #define N 12
 178          char filter()
 179          {
C51 COMPILER V7.06   AUTO                                                                  07/20/2008 16:34:55 PAGE 4   

 180              char count,i,j;
 181              char value_buf[N];
 182              int   sum=0;
 183              for   (count=0;count<N;count++)
 184              {
 185                 value_buf[count] = get_ad();
 186                 delay();
 187              }