c語言 算術平均濾波法 微控制器數字濾波的演算法

2021-10-14 11:31:25 字數 3350 閱讀 7887

微控制器主要作用是控制外圍的器件,並實現一定的通訊和資料處理。但在某些特定場合,不可避免地要用到數**算,儘管微控制器並不擅長實現演算法和進行複雜的運算。下面主要是介紹如何用微控制器實現數字濾波。

在微控制器進行資料採集時,會遇到資料的隨機誤差,隨機誤差是由隨機干擾引起的,其特點是在相同條件下測量同一量時,其大小和符號會現無規則的變化而無法**,但多次測量的結果符合統計規律。為克服隨機干擾引起的誤差,硬體上可採用濾波技術,軟體上可採用軟體演算法實現數字濾波。濾波演算法往往是系統測控演算法的乙個重要組成部分,實時性很強。

採用數字濾波演算法克服隨機干擾的誤差具有以下優點

l 數字濾波無需其他的硬體成本,只用乙個計算過程,可靠性高,不存在阻抗匹配問題。尤其是數字濾波可以對頻率很低的訊號進行濾波,這是模擬濾波器做不到的。

l 數字濾波使用軟體演算法實現,多輸入通道可共用乙個濾波程式,降低系統開支。

l 只要適當改變濾波器的濾波程式或運算,就能方便地改變其濾波特性,這對於濾除低頻干擾和隨機訊號會有較大的效果。

l 在微控制器系統中常用的濾波演算法有限幅濾波法、中值濾波法、算術平均濾波法、加權平均濾波法、滑動平均濾波等。

(1)限幅濾波演算法

該運算的過程中將兩次相鄰的取樣相減,求出其增量,然後將增量的絕對值,與兩次取樣允許的最大差值a進行比較。a的大小由被測物件的具體情況而定,如果小於或等於允許的最大差值,則本次取樣有效;否則取上次取樣值作為本次資料的樣本。

演算法的程式**如下:

#define a //允許的最大差值

char data; //上一次的資料

char filter()

char datanew; //新資料變數

datanew=get_data(); //獲得新資料變數

if( (datanew-data)>a||(data-datanew>a) )

return data;

else

return datanew;

說明:限幅濾波法主要用於處理變化較為緩慢的資料,如溫度、物體的位置等。使用時,關鍵要選取合適的門限制a。通常這可由經驗資料獲得,必要時可通過實驗得到。

(2)中值濾波演算法

該運算的過程是對某一引數連續取樣n次(n一般為奇數),然後把n次取樣的值按從小到大排列,再取中間值作為本次取樣值,整個過程實際上是乙個序列排序的過程。

演算法的程式**如下:

#define n 11 //定義獲得的資料個數

char filter()

char value_buff[n]; //定義儲存資料的陣列

char count,i,j,temp;

for(count=0;countvalue_buf[count]=get_data();

delay(); //如果採集資料比較慢,那麼就需要延時或中斷

for(j=0;jfor(value_buff[i]>value_buff[i+1]

temp=value_buff[i];

value_buff[i]=value_buff[i+1];

value_buff[i+1]=temp;

return value_buff[(n-1)/2];

說明:中值濾波比較適用於去掉由偶然因素引起的波動和取樣器不穩定而引起的脈動干擾。若被測量值變化比較慢,採用中值濾波法效果會比較好,但如果資料變化比較快,則不宜採用此方法。

(3)算術平均濾波演算法

該演算法的基本原理很簡單,就是連續取n次取樣值後進行算術平均。

演算法的程式**如下:

char filter()

int sum=0;

for (count=0;countsum+=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 code jq[n]=; //code陣列為加權係數表,存在程式儲存區

char code sum_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;countvalue_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);

c語言 算術平均濾波法 算術平均濾波C程式

檔名 arithmeticfilter.c 功能描述 算術平均濾波,乙個週期內取樣 n次資料,取n 次取樣值的算術平均值做為這個週期的取樣值 include define uchar unsigned char define uint unsigned int 取樣次數 define n 4 次取樣...

c語言 算術平均濾波法 常見的濾波演算法 C語言

說明 假設從8位ad中讀取資料 若是更高位的ad,可將資料型別定義為int 子程式為get ad 一 限幅濾波法 優點 能有效克服因偶然因素引起的脈衝干擾。缺點 無法抑制那種週期性的干擾,且平滑度差。1 限幅濾波 a值可根據實際情況調整 value為有效值,new value為當前取樣值 濾波程式返...

c語言 算術平均濾波法 C語言十大濾波演算法

限幅濾波法 函式名稱 amplitudelimiterfilter 限幅濾波法 優點 能有效克服因偶然因素引起的脈衝干擾 缺點 無法抑制那種週期性的干擾,且平滑度差 說明 1 呼叫函式 getad 該函式用來取得當前值 2 變數說明 value 最近一次有效取樣的值,該變數為全域性變數 newval...