簡易な2次BandPassIIRフィルタの性能評価を行う。
利用するscilabライブラリは、
https://sites.google.com/site/k3ijogadget/digital-filter
に記載した通り。
今回評価するIIRフィルタの伝達関数H(z) を下記のように定義する。
サンプリング周波数 1000hzを前提とした際の上式のインパルス応答と周波数特性の解析
scilabコマンド
-->sfrq = 1000;tmax = 50;
-->b1 = 0.5; b2 = -0.5; h0 = (1+b2)/2;
-->hz = h0*(1-p^2)/(1-b1*p-b2*p^2);
-->ztspec(hz,1)
インパルス応答と、周波数特性(利得、位相)のグラフ表示
低周波数、高周波数領域になるにつれ、利得が下がっており、バンドパスフィルタ特性を持つフィルタだと言える。