DSP实验七用双线性变换法设IIR滤波器.docxVIP

数学与软件科学学院 实验报告 学期：_2016_至_2017_ 第_一_ 学期 2016年10月26日 课程名称:_数字信号处理_ 专业:_信息与计算科学_实验编号： 07 实验项目：用双线性变换法设IIR滤波器 实验成绩：＿＿＿＿＿实验七 用双线性变换法设IIR滤波器一、实验目的及要求(1) 熟悉用双线性变换法设计IIR滤波器的原理和方法；(2) 掌握数字滤波器的计算机仿真算法；(3) 通过观察对实际心电图信号的滤波作用。二、实验内容(1) 用双线性变换法设计一个巴特沃斯低通IIR数字滤波器，设计的指标参数为：在通内频率低于0.2时，最大衰减度小于1dB,在阻带内[0.3,]频率区间上，最小衰减度大于15dB；(2) 以0.02为采样间隔，打印出数字滤波器在频率区间[0，/2]上的幅频响应特性曲线；(3) 用所设计的滤波器对实际心电图信号采样序列进行仿真滤波处理，并分别打印出滤波前后的心电图信号波形图，观察总结滤波作用和效果。输入序列为： =[-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,00,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0]三、实验准备安装MATLAB的计算机系统。四、实验步骤(该部分不够填写.请填写附页) T=1;Fs=1/T; wpz=0.2;wsz=0.3; wp=2*tan(wpz*pi/2);ws=2*tan(wsz*pi/2);rp=1;rs=15; [N,wc]=buttord(wp,ws,rp,rs,s); [B,A]=butter(N,wc,s); fk=0:1/512:1;wk=2*pi*fk; Hk=freqs(B,A,wk); subplot(2,2,1); plot(fk,20*log10(abs(Hk)));grid on; xlabel(\omega/\pi);ylabel(幅度(dB)); axis([0,1,-100,5]);title((b)); [N,wc]=buttord(wpz,wsz,rp,rs); [Bz,Az]=butter(N,wc); wk=0:pi/512:pi; Hz=freqz(Bz,Az,wk); subplot(2,2,4); plot(wk/pi,20*log10(abs(Hz)));grid on; xlabel(\omega/\pi);ylabel(幅度(dB)); axis([0,1,-100,5]);title((b)); wp =0.2*pi; %通带边缘频率 ws =0.3*pi; %阻带边缘频率 Rp =1; %通带最大衰减 As=15; %阻带最小衰减 T=1; %采样周期 Fs=1/T; Omgp=(2/T)*tan(wp/2); %原型通带频率预修正 Omgs=(2/T)*tan(ws/2); %原型阻带频率预修正 [m,Omgc]=buttord(Omgp,Omgs,Rp,As,s); %计算阶数和截止频率 [ba1,aa1]=butter(m,Omgc,s); %求模拟滤波器系数 [bd,ad]=bilinear(ba1,aa1,Fs); %用双线性变换法求数字滤波器系数b,a [H,w]=freqz(bd,ad); %求数字系统的频率特性 x=[-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,... -16,-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,... 6,4,0,0,0,0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0]; figure(1); %绘制滤波器的频率响应 subplot(211); plot(w/pi,abs(H)); title(幅频响应); axis([0,1,0,1]); subplot(212); plot(w/pi,angle(H)/pi); title(相位响应); axis([0,1,-1,1]); figure(2); %绘制滤波前后的心电图时域特性 subplot(211); stem(x,.); title(心电图信号采样序列x(n)); xlabel(n); ylabel(x(n)); axis([0,56,-100,50]); y=filter(bd,ad,x); subplot(212); stem(y,.); title(滤波后的心电图信号); xlabel(n);