数字滤波器设计

1、实验六 数字滤波器设计实验室名称： 信息学院2204 实验时间：2015年11月26日姓 名： 专业： 通信工程 指导教师：陶大鹏成绩教师签名： 年 月 日一、实验目的1. 会设计满足某个给定幅度或增益响应的无限冲激响应或有限冲激响应的数字滤波器。2. 能设计四种类型滤波器并分析各自的特点。3. 利用matlab绘制滤波器的频率特征曲线并分析滤波器特点。利用matlab设计出符合条件的滤波器。4. 利用matlab绘制滤波器的增益曲线和频率特征曲线。二、实验内容q7.1 用matlab确定一个数字无限冲激响应低通滤波器所有四种类型的最低阶数。指标如下：40khz的抽样率，4khz的通带边界频率

2、，8khz的阻带边界频率，0.5db的带通波纹，40db的最小阻带衰减。评论你的结果。q7.2 用matlab确定一个数字无限冲激响应高通滤波器所有四种类型的最低阶数。指标如下：3500hz的抽样率，1050hz的通带边界频率，600hz的阻带边界频率，1db的带通波纹，50db的最小阻带衰减。评论你的结果。q7.5通过运行程序p7.1来设计巴特沃兹带阻滤波器。写出所产生的传输函数的准确表达式。滤波器的指标是什么？你的设计符合指标吗？使用matlab，计算并绘制滤波器的未畸变的相位响应及群延迟响应。q7.6修改程序p7.1来设计符合习题q7.1所给指标的切比雪夫1型低通滤波器。写出所产生的传输

3、函数的准确表达式。你的设计符合指标吗？使用matlab，计算并绘制滤波器的未畸变的相位响应及群延迟响应。q7.20 使用函数firl，设计一个线性相位有限冲激响应低通滤波器，使其满足习题q7.23给出的指标，并画出其增益和相位响应。使用习题q7.13中用凯泽公式估计出的阶数。用表格形式显示滤波器的系数。你的设计满足指标吗？若不满足，调整滤波器阶数直到设计满足指标。满足指标的滤波器阶数是多少？q7.23 用凯泽窗设计一个有限冲激响应低通滤波器。滤波器的指标是：wp=0.31，ws=0.41，as=50db。注意，函数kaiser需要参数及阶数n的值，它们必须先用式（7.36）和式（7.37）分别

4、算出。你的设计满足指标吗？q7.25 用fir2设计一个95阶有限冲激响应滤波器，它具有三个不同的常数幅度级：在频率范围0到0.25中为0.4，在频率范围0.3到0.45中为1.0，在频率范围0.5到1.0中为0.8.画出所设计的滤波器的幅度响应。你的设计满足指标吗？q7.27 用remez设计具有如下指标的有限冲激响应带通滤波器：通带边界为1.8khz,阻带边界为1.5khz和4.2khz，通带波纹=0.1，阻带波纹=0.02，抽样频率为12khz。用kaiserord估计滤波器的阶数。你的设计是一个最优有限冲激响应滤波器吗？你的设计满足指标吗？若不满足，增加滤波器阶数在满足指标方面有用吗？

5、指标由一个较低阶数的滤波器来满足而不是由kaiserord得到的来满足吗？在不等过渡带的情形下，用remez设计的滤波器可能在较大的过渡带宽中以增益响应表现不满意的行为。改进该行为的一种方法是：通过移动阻带边界减少过渡带宽，直到使设计在过渡带中以平滑的下降来满足指标。在通带边界保持固定的情况下，尝试这种方法并确定新的指标，它在过渡带中提供平滑的下降。三、实验器材及软件1. 微型计算机1台2. matlab 7.0软件四、实验原理1、可以通过几种设计方式，有滤波器指标来直接估计最小滤波器长度n。2、设计无限冲激响应滤波器的方法是基于s平面到z平面的双线性变换。3、对理想频率响应hd（ej）实行离

6、散时间傅里叶变换，可得到其理想无限长冲激响应hdn。4、对于带通和带阻滤波器的设计，使用适当滤波器命令得到的传输函数的实际阶次为2n。五、实验步骤1、 打开计算机中的matlab2、 在matlab中新建一个m-file，在其中输入实验要求的相关代码3、 输入代码后调试代码，代码无误后可以得到相应的实验结果4、 观察实验结果回答相关问题5、记录实验结果，将图形等截图粘贴，完成实验报告六、实验记录（数据、图表、波形、程序等）q 7.1 用matlab确定一个数字无限冲激响应低通滤波器所有四种类型的最低阶数。指标如下：40khz的抽样率，4khz的通带边界频率，8khz的阻带边界频率，0.5db的

7、带通波纹，40db的最小阻带衰减。评论你的结果。答：由题意可知其wp=0.1;ws=0.2巴特沃兹滤波器:n,wn=buttord(0.1,0.2,0.5,40)n = 8 wn = 0.1151切比雪夫1型滤波器：n,wn=cheb1ord(0.1,0.2,0.5,40)n = 5 wn = 0.1000切比雪夫2型滤波器：n,wn=cheb2ord(0.1,0.2,0.5,40)n = 5 wn = 0.2000椭圆滤波器：n,wn=ellipord(0.1,0.2,0.5,40)n = 4 wn =0.1000q7.2用matlab确定一个数字无限冲激响应高通滤波器所有四种类型的最低阶数

8、。指标如下：3500hz的抽样率，1050hz的通带边界频率，600hz的阻带边界频率，1db的带通波纹，50db的最小阻带衰减。评论你的结果。答：由题意可知wp=0.3;ws=0.17巴特沃兹滤波器:n,wn=buttord(0.3,0.17,1,50)n =11 wn = 0.2754切比雪夫1型滤波器：n,wn=cheb1ord(0.3,0.17,1,50)n =6 wn =0.3000切比雪夫2型滤波器：n,wn=cheb2ord(0.3,0.17,1,50)n = 6 wn =0.1700椭圆滤波器：n,wn=ellipord(0.3,0.17,1,50)n =5 wn =0.300

9、0q7.5通过运行程序p7.1来设计巴特沃兹带阻滤波器。写出所产生的传输函数的准确表达式。滤波器的指标是什么？你的设计符合指标吗？使用matlab，计算并绘制滤波器的未畸变的相位响应及群延迟响应。ws = 0.4 0.6; wp = 0.2 0.8; rp = 0.4; rs = 50;n1, wn1 = buttord(wp, ws, rp, rs); num,den = butter(n1,wn1,stop); disp(numerator coefficients are );disp(num);disp(denominator coefficients are );disp(den);

10、g, w = gain(num,den); h = freqz(num, den, w);gd,w1=grpdelay(num,den);plot(w/pi,g);grid axis(0 1 -60 5);xlabel(omega /pi); ylabel(gain in db);title(gain response of a butterworth bandstop filter);pauseplot(w/pi,unwrap(angle(h);gridtitle(phase spectrum argh(ejomega)xlabel(omega /pi);ylabel(phase in ra

11、dians);pauseplot(w1/pi,unwrap(gd);gridtitle(delay characteristic)xlabel(omega /pi);ylabel(group delay);numerator coefficients are columns 1 through 10 0.0493 0.0000 0.2465 0.0000 0.4930 0.0000 0.4930 0.0000 0.2465 0.0000 column 11 0.0493denominator coefficients are columns 1 through 10 1.0000 0.0000

12、 -0.0850 0.0000 0.6360 0.0000 -0.0288 0.0000 0.0561 0.0000 column 11 -0.0008q7.6修改程序p7.1来设计符合习题q7.1所给指标的切比雪夫1型低通滤波器。写出所产生的传输函数的准确表达式。你的设计符合指标吗？使用matlab，计算并绘制滤波器的未畸变的相位响应及群延迟响应。% design of a butterworth bandstop digital filterws = 0.4 0.6; wp = 0.2 0.8; rp = 0.4; rs = 50;% estimate the filter ordern1

13、, wn1 = cheb1ord(wp, ws, rp, rs); % design the filternum,den = cheby1(n1,rp,wn1,stop); % display the transfer functiondisp(numerator coefficients are );disp(num);disp(denominator coefficients are );disp(den);g, w = gain(num,den); h = freqz(num, den, w);gd,w1=grpdelay(num,den);plot(w/pi,g);grid axis(

14、0 1 -60 5);xlabel(omega /pi); ylabel(gain in db);title(gain response of a butterworth bandstop filter);pauseplot(w/pi,unwrap(angle(h);gridtitle(phase spectrum argh(ejomega)xlabel(omega /pi);ylabel(phase in radians);pauseplot(w1/pi,unwrap(gd);gridtitle(delay characteristic)xlabel(omega /pi);ylabel(gr

15、oup delay);numerator coefficients are 0.0326 -0.0000 0.1304 -0.0000 0.1956 -0.0000 0.1304 -0.0000 0.0326denominator coefficients are 1.0000 -0.0000 -1.3141 0.0000 1.3919 -0.0000 -0.7390 0.0000 0.2076q7.20 使用函数firl，设计一个线性相位有限冲激响应低通滤波器，使其满足习题q7.23给出的指标，并画出其增益和相位响应。使用习题q7.13中用凯泽公式估计出的阶数。用表格形式显示滤波器的系数。你

16、的设计满足指标吗？若不满足，调整滤波器阶数直到设计满足指标。满足指标的滤波器阶数是多少？ft=10000; % 采样频率fp=2000; % 通带截止频率fs=2500; % 阻带截止频率dev = 0.005 0.005;fedge = fp fs;aval = 1 0;n,wn,beta,ftype = kaiserord(fedge,aval,dev,ft);b = fir1(n,wn);g, w = gain(b,1); h = freqz(b,1,w);plot(w/pi,g);grid xlabel(omega /pi); ylabel(gain in db);title(gain

17、);pauseplot(w/pi,unwrap(angle(h);gridtitle(phase spectrum argh(ejomega)xlabel(omega /pi);ylabel(phase in radians);分析：是移除跳变后的相位。q7.23 用凯泽窗设计一个有限冲激响应低通滤波器。滤波器的指标是：wp=0.31，ws=0.41，as=50db。注意，函数kaiser需要参数及阶数n的值，它们必须先用式（7.36）和式（7.37）分别算出。你的设计满足指标吗？根据式（7.36）和式（7.37）计算得出n=184， =4.6123。程序如下：w1=kaiser(184,4.

18、6);b=fir1(183,0.31,w1);g, w = gain(b,1); h = freqz(b,1,w);plot(w/pi,g);grid xlabel(omega /pi); ylabel(gain in db);title(gain);pauseplot(w/pi,unwrap(angle(h);gridtitle(phase spectrum argh(ejomega)xlabel(omega /pi);ylabel(phase in radians); 分析：由图可知，该滤波器的设计满足指标。q7.25 用fir2设计一个95阶有限冲激响应滤波器，它具有三个不同的常数幅度级

19、：在频率范围0到0.25中为0.4，在频率范围0.3到0.45中为1.0，在频率范围0.5到1.0中为0.8.画出所设计的滤波器的幅度响应。你的设计满足指标吗？ fpts=0 0.25 0.3 0.45 0.5 1; mval=0.4 0.4 1 1 0.8 0.8; b=fir2(95,fpts,mval); h = freqz(b,1,w); plot(w/pi,abs(h);grid title(amplitude response) xlabel(omega /pi); ylabel(amplitude);q7.27 用remez设计具有如下指标的有限冲激响应带通滤波器：通带边界为1.

20、8khz,阻带边界为1.5khz和4.2khz，通带波纹=0.1，阻带波纹=0.02，抽样频率为12khz。用kaiserord估计滤波器的阶数。你的设计是一个最优有限冲激响应滤波器吗？你的设计满足指标吗？若不满足，增加滤波器阶数在满足指标方面有用吗？指标由一个较低阶数的滤波器来满足而不是由kaiserord得到的来满足吗？在不等过渡带的情形下，用remez设计的滤波器可能在较大的过渡带宽中以增益响应表现不满意的行为。改进该行为的一种方法是：通过移动阻带边界减少过渡带宽，直到使设计在过渡带中以平滑的下降来满足指标。在通带边界保持固定的情况下，尝试这种方法并确定新的指标，它在过渡带中提供平滑的下

21、降。fedge = 1500 1800 3000 4200;aval = 0 1 0;dev = 0.1 0.02 0.1;ft = 12000;n wn beta ftype = kaiserord(fedge,aval,dev,ft);fpts=fedge/ft; % 求解归一化通、阻带的截止频率fpts=0 fpts 1;mval= 0 1 0 0 1 0;b = remez(n,fpts,mval); b = fir1(n,wn);g, w = gain(b,1); h = freqz(b,1,w);plot(w/pi,g);grid xlabel(omega /pi); ylabel(gain in db);title(gain); 分析：观察上图可知，其阻带波纹过大，不符合题意。右因为remez函数采用的是park-mcclellan算法，所以得到最优有限冲激响应滤波器。将滤波器阶数更改滤波器阶数为30，得到结果：下图左为将滤波器阶数该成100时的结果：实验证明：通过增加滤