認(rèn)識(shí) Yn=a* Xn+（1-a）*Yn-1 請(qǐng)進(jìn)

jojolin

低通濾波算法如下：

   Yn=a* Xn+（1-a）*Yn-1

式中  Xn——本次采樣值

      Yn-1——上次的濾波輸出值；

      a——濾波系數(shù)，其值通常遠(yuǎn)小于1；

      Yn——本次濾波的輸出值。


我想問(wèn)下前輩這個(gè)a濾波系數(shù)如何確定呢？它是一個(gè)常數(shù)還是什么？？？



上學(xué)的時(shí)候光幸福了，今天打電話才知道我們?cè)瓉?lái)也開(kāi)過(guò)“信號(hào)與系統(tǒng)”，可是我認(rèn)識(shí)的同學(xué)沒(méi)有會(huì)的，哎......交友不慎阿~！��！

system888net

Yn=a* Xn+（1-a）*Yn-1是為了模擬了具有較大慣性的低通濾波器的功能,
因此a是一個(gè)遠(yuǎn)小于1的常數(shù),否則Yn則受Xn影響較大,根據(jù)實(shí)際情況和經(jīng)驗(yàn)設(shè)頂a的值(Xn對(duì)此慣性系統(tǒng)的影響力).

JohnBull

google一下“指數(shù)平滑”

system888net

原帖由 JohnBull 于 2008-4-17 16:14 發(fā)表
google一下“指數(shù)平滑”

一個(gè)經(jīng)典數(shù)學(xué)模型會(huì)在各個(gè)領(lǐng)域都會(huì)有應(yīng)用.

jojolin

十分感謝兩位，這個(gè)是我寫的關(guān)于此濾波算法的試驗(yàn)程序 。

#include <stdio.h> #define NUM_sig 110 const float come[NUM_sig]={1,2,124,3,-80,5,6,7,8,9,-60,11,12,13,254,15,16,17,78,100,20,21,22,23, 24,25,26,27,28,-100,30,31,32,33,134,35,36,37,38,123,40,41,42,43,44,45,46,47,48,49,50,51,52, 53,54,55,56,57,123,59,60,61,62,63,64,65,66,67,68,69,70,123,72,73,74,75,76,77,10,79,80,81, 82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,12,104,105,106,107,108,109 }; float DT_LB(float Yn_1,float a,float get_AD) { &nbsp;&nbsp;&nbsp;&nbsp;float Yn ; &nbsp;&nbsp;&nbsp;&nbsp;Yn=a*get_AD+(1-a)*Yn_1; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;return Yn ; } void main(int argc, char* argv[]) { &nbsp;&nbsp;&nbsp;&nbsp;float Y=1 ; &nbsp;&nbsp;&nbsp;&nbsp;float LB_xs=0.01 ; &nbsp;&nbsp;&nbsp;&nbsp;int i; &nbsp;&nbsp;&nbsp;&nbsp;for(i=0;i<NUM_sig;i++) &nbsp;&nbsp;&nbsp;&nbsp;{ &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Y=DT_LB(Y,LB_xs,come[i]) ;   //濾波處理 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;printf("%0.0f,",Y); &nbsp;&nbsp;&nbsp;&nbsp;} }

假設(shè)定義的數(shù)流為

const float come[NUM_sig]={1,2,124,3,-80,5,6,7,8,9,-60,11,12,13,254,15,16,17,78,100,20,21,22,23,
24,25,26,27,28,-100,30,31,32,33,134,35,36,37,38,123,40,41,42,43,44,45,46,47,48,49,50,51,52,
53,54,55,56,57,123,59,60,61,62,63,64,65,66,67,68,69,70,123,72,73,74,75,76,77,10,79,80,81,
82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,12,104,105,106,107,108,109
}

系數(shù)是0.01 處理后的結(jié)果為

1,1,2,2,1,1,2,2,2,2,1,1,1,1,4,4,4,4,5,6,6,6,6,7,7,7,7,7,8,6,7,7,7,7,9,9,9,10,10,
11,11,12,12,12,12,13,13,13,14,14,15,15,15,16,16,16,17,17,18,19,19,20,20,20,21,21
,22,22,23,23,24,25,25,25,26,26,27,27,27,28,28,29,29,30,30,31,32,32,33,33,34,34,3
5,36,36,37,37,38,38,39,40,40,41,41,41,42,43,43,44,44


可以看到紅色的值都濾掉了，可是這樣的算法滯后特別嚴(yán)重，110個(gè)數(shù)據(jù)到最后只有44，我現(xiàn)在定義的濾波系數(shù)是0.01，這是否就得根據(jù)實(shí)際情況找個(gè)合適的值了？？？

系數(shù)是0.1的結(jié)果：

1,1,13,12,3,3,4,4,4,5,-2,-0,1,2,27,26,25,24,30,37,35,34,32,31,31,30,30,29,29,16,
18,19,20,22,33,33,33,34,34,43,43,43,43,43,43,43,43,44,44,45,45,46,46,47,48,48,49
,50,57,57,58,58,58,59,59,60,61,61,62,63,63,69,70,70,70,71,71,72,66,67,68,70,71,7
2,73,74,76,77,78,79,80,81,82,83,84,85,87,88,89,90,91,92,93,85,87,88,90,92,93,95,

system888net

原帖由 jojolin 于 2008-4-17 18:10 發(fā)表
十分感謝兩位，這個(gè)是我寫的關(guān)于此濾波算法的試驗(yàn)程序 。


#include

#define NUM_sig 110

const float come[NUM_sig]={1,2,124,3,-80,5,6,7,8,9,-60,11,12,13,254,15,16,17,78,100,20,21,22,23,
24 ...

不錯(cuò),多取幾個(gè)a值試試.

system888net

如果你的機(jī)器夠快,寫個(gè)循環(huán)取a自動(dòng)檢查.

jojolin

多謝 ，
也看到了 您在另一個(gè)帖子中的闡述，呵呵，提醒我得改為定點(diǎn)運(yùn)算，我用的2410不支持浮點(diǎn)的。

system888net

原帖由 jojolin 于 2008-4-17 18:26 發(fā)表
多謝 ，
也看到了 您在另一個(gè)帖子中的闡述，呵呵，提醒我得改為定點(diǎn)運(yùn)算，我用的2410不支持浮點(diǎn)的。

[ 本帖最后由 system888net 于 2008-4-21 17:28 編輯 ]