数字信号处理实验四-IIR及FIR滤波器设计.pdf

西安交通大学实验报告 成绩 课 程： 数字信号处理 第 1 页共 页 系 别： 生命学院 实 验 日 期： 年 月 日 专业班级： 医电 53 组别： 交 报告日 期： 年 月 日 姓 名：李竞捷 学号：2151500084 报 告 退 发： （订正、重做） 同 组 者： Null 教师审批签字： 实验名称： IIR 及 FIR 滤波器设计 一、 实验目的 学习用 MATLAB 完成 IIR 及 FIR 数字滤波器的设计 了解 IIR 及 FIR 数字滤波器的设计方法; 二、 实验过程 1. 用冲激响应不变法设计Butterworth IIR数字低通滤波器。其中通带边界频 率fp=1000Hz，阻带边界频率fst=4000Hz，通带最大波动Rp=3dB， 阻带最小衰 减As=15dB，抽样频率fs=20kHz，求H(z)及频率响应H(ejw)并画出幅频曲线。 2. 已知输入信号x=sin(2 π20t)+sin(2 π5000t), 抽样频率 fs=20kHz，请用步 骤1设计的滤波器过滤该信号，要求画出输入与输出信号的时域波形并求它们的 谱。 3. 用窗函数法设计FIR数字低通滤波器。设通带边界频率fp=50Hz，阻带边界频 率fst=100Hz，抽样频率fs=400Hz，阻带最小衰减As=40dB。 请选择合适的窗函 数及阶数n进行设计，求出h(n)与H(z)及频率响应H(j ω) 并画出幅频曲线。 实验一：代码 wp=2*pi*1000; ws=2*pi*4000; rp=3; rs=15; fs=20000; [n,wn]=buttord(wp,ws,rp,rs,'s') [z,p,k]=buttap(n) [b,a]=zp2tf(z,p,k); [bt,at]=lp2lp(b,a,wp) [bz,az]=impinvar(bt,at,fs) freqz(bz,az) 运行结果： 实验二： ts=0.2; t=0:1/fs:ts; ff=-0.5*fs:1/ts:0.5*fs; x=sin(2*pi*20*t)+sin(2*pi*5000*t); y=filter(bz,az,x); X=fftshift(fft(x)); Y=fftshift(fft(y)); figure subplot(1,2,1) plot(t,x); subplot(1,2,2)