IIR和FIR数字滤波器过滤信号的实现及比较

信号、系统与信号处理实验Ⅱ实验报告姓名：杨海清学号：14081331班级：14083413上课时间：2016.12.1实验名称：IIR和FIR数字滤波器过滤信号的实现及比较一、实验目的（1）掌握数字滤波器的计算机仿真方法；（2）通过观察对实际心电图信号的滤波作用，获得数字滤波器的感性知识。二、实验原理（1）编写FIR数字滤波器实现程序，计算其对心电图信号取样序列x(n)的响应序列yFIR(n)。（2）编写IIR滤波器实现程序，计算其对心电图信号取样序列x(n)的响应序列yIIR(n)。（3）在通过计算机上运行上述两个滤波器实现程序，在屏幕上打印出x(n)、yFIR(n)和yIIR(n)，并进行比较。三、仿真结果分析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];subplot(3,1,1);plot([0:length(x)-1],x,'blu');title('输入')N=21;wn=hamming(N+1);wp=0.2*pi;h=fir1(N,wp/pi,wn);y=conv(x,h);subplot(3,1,2);plot([0:length(y)-1],y,'red');title('经FIR滤波')A=0.2318;a=A*[1,2,1];b1=[1,0.4164,-0.6006];b2=[1,0.3093,-0.1888];b3=[1,0.2692,-0.0349];y1=filter(a,b1,x);y2=filter(a,b2,y1);y3=filter(a,b3,y2);subplot(3,1,3);plot([0:length(y3)-1],y3,'gre');title('经IIR滤波')四、实验总结对于FIR滤波器主要是使用了前面的fir1函数，求得系统的冲激响应函数，再利用conv函数求的输出，而对于IIR滤波器书上给出了级联型系统函数的各个系数，利用filter函数可以直接求得输出。