数字信号处理实验指导书(四).docVIP

“数字信号处理”实验指导书（四） 一、实验课程编码：105003 二、实验课程名称：数字信号处理 三、实验项目名称：应用MATLAB设计IIR数字滤波器 四、实验目的 掌握应用MATLAB设计IIR数字滤波器的方法，即熟悉应用MATLAB设计IIR滤波器相关的函数。 五、主要设备 安装有MATLAB软件的电脑 六、实验内容 编写MATLAB程序，实现下列题目： 1．设计数字低通滤波器，要求在通带内频率低于0.2πrad时，允许幅度误差在1dB以内，在频率0.3πrad~πrad之间的阻带衰减大于15dB， 即ωp=0.2πrad，Rp＝1dB，ωs=0.3πrad，AS ＝15dB 用脉冲响应不变法设计数字滤波器，模拟滤波器采用切比雪夫I型滤波器原型，T＝1。画出所设计的滤波器的幅度响应。 2．设计低通数字滤波器，要求 ωp=0.2πrad , Rp =1dB, ωs=0.3πrad, AS=15dB 用双线性变换法设计数字滤波器，模拟滤波器采用巴特沃斯滤波器原型，T＝1。画出所设计的滤波器的幅度响应。并与上题结果进行比较。 3.设计一个巴特沃斯高通滤波器，要求 ωp=0.6πrad , Rp =1dB, ωs=0.4πrad, AS=15dB 用MATLAB提供的直接设计IIR数字滤波器的函数buttord及butter。T＝1。 4. 利用切比雪夫I型设计一个数字高通滤波器，使其满足： 采用数字域频率变换法、双线性变换法。T＝1。 5. 利用切比雪夫I型设计一个数字带通滤波器，使其满足： 采用数字域频率变换法、双线性变换法。T＝1。 七、实验步骤 1、熟悉与离散信号频谱分析相关的MATLAB函数（参考附录1）； 2、通过运行附录2中提供的例题，熟悉用MATLAB设计IIR数字滤波器的基本方法； 3、根据“六、实验内容”中各个题目的要求，编写MATLAB程序代码，调试程序，分析并保存结果。 八、实验结果 对实验练习题编写MATLAB程序并运行，在计算机上输出仿真结果。 附录1 主要的相关MATLAB函数 1.buttord.m 调用格式：〔N,Wn〕＝buttord（Wp，Ws， Rp， Rs） 〔N,Wn〕＝buttord（Wp，Ws， Rp， Rs，’s’） 2.buttap.m 调用格式：〔z,p,k〕＝buttap（N） 3.butter.m 调用格式：〔B,A〕＝butter（N,Wn,’TYPE’） 4.cheb1ord.m 5.cheb1ap.m 6.cheby1.m 7.cheb2ord.m 8.cheb2ap.m 9.cheby2.m 文件4～9的调用格式与文件1~3的调用格式类似。 10.bilinear.m 调用格式：〔Bz,Az〕＝bilinear（B，A，Fs） 11.impinvar.m 调用格式：〔Bz,Az〕＝impinvar（B，A，Fs） 附录2 例题 例1、设计一个低通模拟巴特沃斯滤波器，使其满足： 通带： 阻带： [MATLAB程序]： % design Butterworth lowpass filter. % Wp,Ws:rad/s,Rp,As:dB Wp=20000*pi; Ws=30000*pi; Rp=1; As=15; Ripple=10^(-Rp/20); Attn=10^(-As/20); % Analog filter design; [b,a]=afd_butt(Wp,Ws,Rp,As); %calculation of second-order sections: [C,B,A]=sdir2cas(b,a); % calculation of frequency response; [dB,mag,pha,w]=freqs_m(b,a,50000*pi); % calculation of impulse response: [ha,x,t]=impulse(b,a); % plot figures figure(1),clf; subplot(2,2,1); plot(w/pi,mag); title(Magnitude Response); xlabel(Analog frequency in pi units); ylabel(H); axis([0,50000,0,1.1]); set(gca,XTickMode,manual,XTick,[0,20000,30000,50000]); set(gca,YTickMode,manual,YTick,[0,Attn,Ripple,1]); grid; subplot(2,2,2); plot(w/pi,dB); title(Magnitu