* 超てきとーな 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 になってるはず.