无限长单位脉冲响应滤波器设计剖析.docxVIP

实验四 无限长单位脉冲响应滤波器设计 一、实验目的 1.掌握双线性变换法及脉冲响应不变法设计 IIR 数字滤波器的具体设计方法及其原理，熟悉 用双线性变换法及脉冲响应不变法设计低通、高通和带通 IIR 数字滤波器的计算机编程。 2.观察双线性变换及脉冲响应不变法设计的滤波器的频率特性， 了解双线性变换法及脉冲响 应不变法的特点。 3.熟悉巴特沃思滤波器、切比雪夫滤波器和椭圆滤波器的频率特性。 二、实验原理 实验中有关变量的定义： fc 通带边界频率， fr 阻带边界频率， tao 通带波动， at 最小阻带衰减， fs 采样频率， t 采样周期。 设计一个数字滤波器一般包括以下两步： 按照任务要求，确定滤波器性能指标 用一个因果稳定的离散时不变系统的系统函数去逼近这一性能要求 数字滤波器的实现：对于 IIR 滤波器，其逼近问题就是寻找滤波器的各项系数，使其系统函数逼近一个所要求的特性。 先设计一个合适的模拟滤波器， 然后变换成满足约定指标的数字滤波器。 用双线形变换法设计 IIR 数字滤波器的过程： a. 将设计性能指标中的关键频率点进行“预畸” b. 利用“预畸”得到的频率点设计一个模拟滤波器。 c. 双线形变换，确定系统函数 三、实验内容 1、设计一切比雪夫高通滤波器，性能指标如下：通带边界频率 fc=0.4 kHz ，通带波动 δ=0.5 dB ，阻带边界频率 fr s =0.3 kHz ，阻带最小衰减 At=20dB，采样频率 f =1000Hz ，观察其 通带波动和阻带衰减是否满足要求。 （绘制对数幅度谱） 2、设计一巴特沃思低通滤波器，性能指标如下：通带边界频率 fc=0.4 kHz ，通带波动 δ=1dB，阻带边界频率 fr =0.6 kHz ，阻带最小衰减 At=40dB，采样频率 fs=2000Hz，分别用脉 冲响应不变法和双线性变换法进行设计，比较两种方法的优缺点。 （绘制线性幅度谱） 3、用双线性变换法设计巴特沃思、切比雪夫和椭圆低通滤波器，其性能指标如下：通 带边界频率 fc=1.8 kHz ，通带波动 δ≤ 1dB，阻带边界频率 fr=2.6 kHz ，阻带最小衰减 At ≥ 50dB，采样频率 fs=8kHz 。（绘制对数幅度谱） 4、设计一巴特沃思带通滤波器，性能指标如下：通带频率 3kHz≤ f≤ 4kHz ，通带波动 δ≤ 1dB；上阻带 f≥ 5kHz ，阻带最小衰减 At≥ 15dB；下阻带 f≤ 2kHz ，阻带最小衰减 At≥ 20dB； 采样频率 fs=20kHz ，分别用脉冲响应不变法和双线性变换法进行设计。 （绘制线性幅度谱） 5、用双线性变换法设计一椭圆带阻滤波器， 性能指标如下： 阻带频率 2kHz≤ f≤ 3kHz ， - 1 - 阻带最小衰减 At≥ 30dB ，通带频率 f ≥ 4kHz 和 f≤ 1kHz ，通带波动 δ≤ 3dB ，采样频率 fs=10kHz 。（绘制对数幅度谱） 注意：设计结果要求给出程序、数字滤波器的 系统函数 H (z) 和幅频响应 ， H ( z) 请用 word 自带的公式编辑器编辑（插入－对象－ Microsoft 公式 3.0 ）。 四、实验结果与分析 1.实验结果： 切比雪夫高通滤波器 5 0 -5 -10 d / 度 幅 -15 -20 -25 -30 0 50 100 150 200 250 300 350 400 450 500 频 率 /Hz 程序： fs=1000;fc=400;fr=300;rp=0.5;rs=20; ws=2*fs*tan(2*pi*fc/fs/2); wp=2*fs*tan(2*pi*fr/fs/2); [N,wn]=cheb1ord(wp,ws,rp,rs,s); [b,a]=cheby1(N,rp,wp,high,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,500,-30,5]); grid; xlabel( 频率 /Hz); ylabel( 幅度 /dB); title( 切比雪夫高通滤波器 ); 实验结果分析： 由图像观察可知，通带波动为 0.5dB，阻带最小衰减 20dB,满足要求。 2.实验结果： - 2 - 巴特沃思低通滤波器，线性幅度谱 1.4 脉冲响应不变法 双线性变换法 1.2 1 0.8 度 幅 0.6 0.4 0.2 0 0 100 200 300 400 500 600 700 800 900 1000 频 率 /Hz 程序： fs=2000;fc=400;fr=600;