Tweet


* 超てきとーな IIR ディジタルフィルタ,の評価 [#ja7a1712]

octave を使って,「超てきとーな IIR ディジタルフィルタ」の評価をしてみます.
てきとーなのはフィルタの方で,評価じゃないよ.
とはいっても,かっちりした評価でもないけれど.

** てきとーな IIR ディジタルフィルタ [#jc6542f0]

#ref(iir.png,center)

こんな感じの構成の IIR フィルタを考えます.
a0 = 0.9 とします.

** インパルス応答を計算する [#zceb9829]

シミュレーションするまでもなく,インパルス応答は

 h(n) = a0^(n+1)

ですね.

ということで,結果を直接代入します.

 octave:1> h = power(0.9,1:1024);
 octave:2> plot(h)

#ref(iir_impulse_resp.png,center)

FFT にかけます.

 octave:3> y = fft(h);

プロットしてみます.

 octave:5> plot ( 20*log10(abs(y)/max(abs(y))))

#ref(iir_freq_resp.png,center)

このグラフは横軸がリニアなので妙な形になってますが,対数軸にしてやれば見慣れた -6dB/oct の 1 次の LPF になってるはず.

トップ   編集 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS