bz2 =[0.0014 -0.0000 -0.0042 -0.0000 0.0042 -0.0000 -0.0014] az2 =[1.0000 -4.8071 10.2473 -12.2838 8.7245 -3.4849 0.6176] 因此双线性变换法的系统函数为:

H(z)2

0.0014 0.0042z

1 4.8071z

1

2

0.0042z

3

4

0.0014z

4

6

5

10.2473z

2

12.2838z 8.7245z 3.4849z 0.6176z

6

（5）利用双线性变换法设计满足下列指标的切比雪夫型数字带阻滤波器，并作图验证设计结果：当1kHz f 2kHz时，At 18dB，当f 500Hz以及f 3kHz时， 3dB；采样频率fs 10kHz。 解: 程序:

clear;

fs=10000;fc=[500 3000]; fr=[1000 2000];rp=3;rs=20; wp=2*fs*tan(2*pi*fc/fs/2); ws=2*fs*tan(2*pi*fr/fs/2); [N,wn]=cheb1ord(wp,ws,rp,rs,'s'); [b,a]=cheby1(N,rp,wp,‘stop’,'s');

[bz,az]=bilinear(b,a,fs);

[h,w]=freqz(bz,az); f=w/(2*pi)*fs; figure;plot(f,20*log10(abs(h))); axis([0,2000,-200,10]);

grid; xlabel('频率/Hz'); ylabel('幅度/dB'); title('切比雪夫带阻滤波器');

切比雪夫带通滤波器

-20

-40

幅度/dB

-60

-80

-100

-120

0500100015002000

25003000频率/Hz

3500400045005000