lab3iirbilinear.m

基于双线性变换法的IIR滤波器范例 该滤波器为低通原型 可根据需要自行转换成高通 带通 以及重叠累加和重叠保留的例子

close all;
clear all;

%[Y,FS,NBITS]=wavread('Adiemus.wav')
[N, Wn] = buttord(1, 1.6974098, 0.1, 70,'s');
[B, A] = butter(N, Wn ,'s');

[Bs, As] = lp2bp(B, A, 0.608434, 0.2170171);
[Bz, Az] = bilinear(Bs, As, 0.5);

[H,W] = freqz(Bz, Az, 512);
plot(W/pi, 20*log10(abs(H)));
grid;