數字信号处理实验报告实验三.docVIP

数字信号处理实验报告 实验三：IIR滤波器设计实验 班 级： 姓 名： 徐可 学 号：1308024121 实验三 IIR滤波器设计实验 一、实验目的 1.学习模拟－数字变换滤波器的设计方法 2.掌握双线性变换滤波器的设计方法 3.掌握实现数字滤波的具体方法。 二、实验要求 1. 用双线性变换法设计一个巴特沃斯低通IIR数字滤波器。设计指标参数为：在通带内频率低于0.2π时，最大衰减小于1dB；在阻带内[0.3π, π]频率区间上，最小衰减大于15dB. 2.0.02π为采样间隔，打印出数字滤波器在频率区间[0, π/2]上的频率响应特性曲线。 3. 用所设计的滤波器对实际心电图信号采样序列进行仿真滤波处理，观察总结滤波作用与效果 附：心电图采样序列x(n) 人体心电图信号在测量过程中往往受到工业高频干扰，所以必须经过低通滤波处理后，才能作为判断心脏功能的有用信息。下面给出一实际心电图信号采样序列样本x(n),其中存在高频干扰。在实验中以x(n)作为输入序列，滤除其中的干扰成分。 {x(n)}={-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} 三、实验设备 计算机、Matlab软件 四、实验报告要求 1.给出详细的滤波器设计说明书； 2.整理好经过运行并证明是正确的程序并且加上详细注释； 3.用所设计的滤波器对心电信号进行滤波，打印滤波后的波形； 五、实验结果及分析 IIR滤波器设计步骤： 将数字滤波器的技术指标转换为巴特沃斯模拟低通滤波器的技术指标。 算出巴特沃斯低通滤波器的阶数及 确定归一化模拟滤波器模型的零极点形式，确定。 用双线性变换法将转换为。 根据双线性法设计滤波器的步骤实现。 ap=1;as=15; %数字滤波器指标 T=0.02*pi; %采样频率与间隔 fsa=1/T; wp=0.2*pi;ws=0.3*pi; %转换为数字角频率 Wp=2/T*tan(wp/2);Ws=2/T*tan(ws/2); %由数字角频率转换为模拟角频率 [N,Wc]=buttord(Wp,Ws,ap,as,s); %获取模拟滤波器的阶数和3dB截止频率 [Z,P,K]=buttap(N); %归一化模拟滤波器模型的零极点形式参数 [B,A]=zp2tf(Z,P,K); %归一化模拟滤波器传递函数的系数 [Bl,Al]=lp2lp(B,A,Wc); %把模拟滤波器原型转换成截至频率为Wc的低通滤波器 [b,a]=bilinear(Bl,Al,fsa);%用双线性变换法实现模拟滤波器到数字滤波器的转换 [H,w]=freqz(b,a); %获取频率响应 L=length(w)/2+1; figure(1),plot(w(1:L)/pi,abs(H(1:L))),grid,xlabel(角频率（\pi）),ylabel(频率响应幅度); %plot(w*fsa/(2*pi),abs(H));grid; %绘制频率响应曲线 %xlabel(频率(Hz));ylabel(频率响应幅度); 2.消除心电图中干扰成分 ap=1;as=15; %数字滤波器指标 T=0.02*pi; %采样频率与间隔 fsa=1/T; wp=0.2*pi;ws=0.3*pi; %转换为数字角频率 Wp=2/T*tan(wp/2);Ws=2/T*tan(ws/2); %由数字角频率转换为模拟角频率 [N,Wc]=buttord(Wp,Ws,ap,as,s); %获取模拟滤波器的阶数和3dB截止频率 [Z,P,K]=buttap(N); %归一化模拟滤波器模型的零极点形式参数 [B,A]=zp2tf(Z,P,K); %归一化模拟滤波器传递函数的系数 [Bl,Al]=lp2lp(B,A,Wc); %把模拟滤波器原型转换成截至频率为Wc的低通滤波器 [b,a]=bilinear(Bl,Al,fsa);%用双线性变换法实现模拟滤波器到数字滤波器的转换 [H,w]=freqz(b,a); %获取频率响应 L=length(w)/2+1; figu