日記/2010-05-02/超てきとーな IIR ディジタルフィルタ，の評価

超てきとーな IIR ディジタルフィルタ，の評価 †

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

てきとーな IIR ディジタルフィルタ †

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

インパルス応答を計算する †

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

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

ですね．

ということで，結果を直接代入します．

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

FFT にかけます．

octave:3> y = fft(h);

プロットしてみます．

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

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