数字信号处理上机实验报告.doc

实践 创新 努力 奋斗 PAGE PAGE 1 实验报告 数字信号处理实验 实 验 报 告 专业班级： 姓名： 学号： 20题 、调用MATLAB工具箱函数firl设计线性相位低通FIR滤波器，要求希望逼近的理想低通滤波器通带截止频率w=pi/4rad,滤波器长度N=21.分别选用矩形窗、Hanning窗、Hamming窗和Blackman窗进行设计，绘制用每种窗函数设计的单位脉冲响应h(n)及其损耗函数曲线，并进行比较，观察各种窗函数的设计性能。 实验代码如下： N=input(输入窗函数长度N=?(输入0=退出)); %注意加分号与不加分号的区别 while(N~=0) wc=input(输入希望逼近的理想低通滤波器的截止频率Wc=?); %注意截止频率pi/4的输入，matlab中已经默认定义了pi n=0:(N-1); alpha=(N-1)/2; m=n-alpha+eps; hd=sin(wc*m)./(pi*m); %得到理想低通滤波器 k=input(请选择窗函数类型(1=矩形；2=汉宁；3=海明；4=布莱克曼):); if k==1 B=boxcar(N); %产生矩形窗 string=[Boxcar, N=,num2str(N)]; %text函数使用字符串string在图中标明所用窗的类型及长度 elseif k==2 B=hamming(N); string=[Hamming, N=,num2str(N)]; elseif k==3 B=hanning(N); string=[Hanning, N=,num2str(N)]; else B=blackman(N) string=[Blackman, N=,num2str(N)]; end h=hd.*(B); %得到FIR数字滤波器h(n)=hd(n)w(n), 注意*是矩阵相乘，.*是矩阵的对应元素相乘 [H,w]=freqz(h,[1],1024); %求滤波器h(n)的频率响应；对FIR而言, H(z)分子分母多项式的系数向量b=[1], a=h；返回向量H的点数N =1024 db=20*log10(abs(H)+eps); %得到幅值 pha=angle(H); %得到相位 %绘制单位脉冲响应h(n)、幅频衰减特性20lg︱H(ejw)︱)、相频特性和幅频特性︱H(ejw)︱的波形 figure; subplot(2,1,1) stem(n,h,.r); %绘制单位脉冲响应h(n)波形 axis([0,N-1,-0.1,0.3]); xlabel(n); ylabel(h(n)); title(实际低通滤波器的h(n)); text((0.28*N),0.275,string); %创建text图形句柄对象函数，text(x,y,string)在图形中指定的位置上(x,y)上显示字符串string subplot(2,1,2) plot(w/pi,db,m); %绘制幅频衰减特性20lg︱H(ejw)︱)波形 axis([0,1,-100,0]); xlabel(w/pi); ylabel(dB); title(衰减特性(dB)); end 实验结果如图： 实验四 IIR数字滤波器设计及软件实现 一、实验目的 （1）熟悉用双线性变换法设计IIR数字滤波器的原理与方法； （2）学会调用MATLAB信号处理工具箱中滤波器设计函数（或滤波器设计分析工具fdatool）设计各种IIR数字滤波器，学会根据滤波需求确定滤波器指标参数。 （3）掌握IIR数字滤波器的MATLAB实现方法。 （4）通过观察滤波器输入输出信号的时域波形及其频谱，建立数字滤波的概念。 二、实验原理 设计IIR数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法），应用最广泛的是双线性变换法。基本设计过程是：①先将给定的数字滤波器的指标转换成过渡模拟滤波器的指标； ②设计过渡模拟滤波器；③将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。MATLAB信号处理工具箱中的各种IIR数字滤波器设计函数都是采用双线性变换法。第六章介绍的滤波器设计函数butter、cheby1 、cheby2 和ellip可以分别被调