数字滤波器设计

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时的结果：试验证明：通过增加滤波器阶数可