中值濾波器超高效實現

Question

我正在尋找快速/高效中值濾波器的Ansi C實現。任何指針？

到目前爲止，我找到了following的實現，這很好，但我對更快的實現感到好奇。我只需要1維。

Answer 1

我需要從非常嘈雜的CPU消耗數據中提取信號。這裏是Jeff McClintock中值濾波器...

將平均值和中值初始化爲零，然後對於每個樣本的「英寸」，向輸入樣本的中值增加一個小增量。最終，它將在大約50％的投入樣本更大，50％比中位數少的點上解決。 增量的大小應該與實際的中位數成正比。由於我們不知道實際的中位數，我用平均數作爲粗略估計。步長計算爲估計值的0.01倍。越小的步長越準確，但需要更長時間才能解決。

float median = 0.0f; 
float average = 0.0f; 

// for each sample 
{ 
    average += (abs(sample) - average) * 0.1f; // rough running average magnitude. 
    median += _copysign(average * 0.01, sample - median); 
}