用双线性变换法设计iir数字左右滤波器实验报告.docxVIP

用双线性变换法设计iir数字滤波器实验报告 　　数学与软件科学学院实验报告　　学期：13至14第_2学期XX年5月3日　　课程名称：数字信号处理专业:信息与计算机科学实验编号：7实验项目：7　　指导教师姓名：　　学号：实验成绩：　　实验七用双线性变换法设IIR滤波器　　一、实验目的　　(1)熟悉用双线性变换法设计IIR滤波器的原理和方法；　　(2)掌握数字滤波器的计算机仿真算法；　　(3)通过观察对实际心电图信号的滤波作用。　　二、实验内容　　(1)用双线性变换法设计一个巴特沃斯低通IIR数字滤波器，设计的指标参数为：在通内频率低于?时，最大衰减度小于1dB,在阻带内[?,?]频率区间上，最小衰减度大于15dB；　　(2)以?为采样间隔，打印出数字滤波器在频率区间[0，?/2]上的幅频响应特性曲线；　　(3)用所设计的滤波器对实际心电图信号采样序列进行仿真滤波处理，并分别打印出滤波前后的心电图信号波形图，观察总结滤波作用和效果。　　输入序列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的计算机系统。　　四、实验步骤及结果　　一.用双线性变换法设计一个巴特沃斯低通IIR数字滤波器，设计的指标参数为：　　在通内频率低于?时，最大衰减度小于1dB,在阻带内[?,?]频率区间　　上，最小衰减度大于15dB；　　实验代码如下：　　wp=*pi;ws=*pi;rp=1;as=15;　　dispf=20;dispa=-80;　　gridon;　　axis([0,dispf,dispa,5]);　　[n,wc]=buttord(wp,ws,rp,as,s);　　[B,A]=butter(n,wc,s);　　k=0:511;fk=0:1400/512:1400;wk=2*pi*fk;　　Hk=freqs(B,A,wk);　　plot(fk/1000,20*log10(abs(Hk)));　　xlabel(频率,FontSize,16,FontWeight,bold);　　ylabel(幅度,FontSize,16,FontWeight,bold);　　legend(巴特沃斯滤波,Location,East);　　grid;　　程序运行后的结果如下：　　三.用所设计的滤波器对实际心电图信号采样序列进行仿真滤波处理，并分　　别打印出滤波前后的心电图信号波形图，观察总结滤波作用和效果。输入序列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]　　实验代码如下：　　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];　　k=1;　　closeall;　　figure(1);　　subplot(2,2,1);　　n=0:55;　　stem(n,x,.);　　axis([0056-10050]);　　holdon;　　n=0:60;　　m=zeros(61);　　plot(n,m);　　xlabel(n,color,g);　　ylabel(x(n),color,g);　　title(心电图信号采样序列x(n),color,r);　　B=[2*];　　A=[-];　　A1=[-];　　A2=[-];　　while(k;color,r);　　title(滤波器的幅频响应曲线,color,g);　　程序运行结果为：　　《数字信号处理》　　实践报告　　题目：实验三用双线性变换法设计IIR数字滤波器　　实验三：用双线性变换法设计IIR数字滤波器　　1.实验目的　　(1)熟悉用双线性变换法设计IIR数字滤波器的原理与方法。　　(2)掌握数字滤波器