BlackFinを使ってみるなり(29) 間引きフィルタ

前回はFIR filterを実装してみました．

今回は、blackfinライブラリのfilter.hにある間引き機能（decimation）を動かしてみます．

まず間引かない普通のFIRのおさらいをします．

用意するもの．構造体、FIFO．

_fir_fr32_state firparam;

fract32 fir_input[64], fir_output[64];

fract32 fir_delay[2 * FILTER_TAP_NUM];

64の意味は、FIRに64 sampleのデータを一括入力するという意味です．出力も64 sampleが出てきます．

FIR構造体を設定します．

fir_init(firparam, filter_taps, fir_delay, FILTER_TAP_NUM, 1);

引数１　構造体

引数２　タップ配列

引数３　work area　タップ数の２倍

引数４　タップ数

引数５　間引かないので１

FIRを呼び出します．

fir_fr32 (

　fir_input,

　fir_output,

　64,

　&firparam　　);

ーーーー

次は、1/32に間引くやりかたです．

つまり、ADCはsampling freq. 32kHz ですが、これを1kHz sampleに変換するという意味です．

用意するもの．構造体、FIFO．

_fir_fr32_state firparam;

fract32 fir_input[64], fir_output[2];

fract32 fir_delay[2 * FILTER_TAP_NUM];

outputサイズが２になってます．理由は、1/32に間引くので64個が２個に減っちゃうからです．

FIR構造体を設定します．

fir_init(firparam, filter_taps, fir_delay, FILTER_TAP_NUM, 32);

引数５　1/32に間引くので３２

FIRのdecimationバージョンを呼び出します．

fir_decima_fr32 (

　fir_input,

　fir_output,

　64,

　&firparam　　);

データを間引いたメリットはこうです．

間引く前：

ADCのcallbackは32kHzの６４サンプル分＝2mSec周期です．2mSecで６４個のデータを処理しなくちゃいけません．

1/32に間引くと：

2mSecで２個のデータを処理するだけですのでとても軽くなりました．いろんなフィルタ処理をやりやすくなるでしょう．

次はFFTやってみる．

２８へ　　　　３０へ

かしこ