数字信号处理报告实验8----实验报告8

《数字信号处理报告实验8----实验报告8》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验报告实验名称______________________课程名称院系部:专业班级：学生姓名：学号:同组人:实验台号:指导教师：成绩：实验日期:华北电力大学（北京）一：实验目的及要求加深理解FIR数字滤波器的时域特性和频域特性，掌握FIR数字滤波器的设计原理与设计方法，以及FIR数字滤波器的应用。二：实验仪器仪器名称规格/型号数量备注计算机一台装有Matlab三：实验原理FIR数字滤波器可以设计成具有线性相位，在数据通信、图像处理、语音信号处理等实际应用领域得到广泛应用。M阶FIR数字滤波器的系统函数为：FI

2、R数字滤波器的单位脉冲响应h[k]是长度为M+1的有限长因果序列。当满足对称条件时，该FIR数字滤波器具有线性相位。FIR数字滤波器设计方法主要有窗口法、频率取样法及优化设计法。MATLAB中提供的常用FIR数字滤波器设计函数有：fir1窗函数法设计FIR数字滤波器（低通、高通、带通、带阻、多频带滤波器）fir2频率取样法设计FIR数字滤波器：任意频率响应firlsFIR数字滤波器设计：指定频率响应firrcos升余弦型FIR数字滤波器设计intfilt内插FIR数字滤波器设计kaiserord凯塞(Kaise

3、r)窗函数设计法的阶数估计firpmParks-McClellan算法实现FIR数字滤波器优化设计firpmordParks-McClellan数字滤波器的阶数选择cremez复系数非线性相位FIR等波纹滤波器设计1.窗口法设计FIR数字滤波器fir1函数可以很容易地实现FIR数字滤波器窗口法设计。可设计低通、高通、带通、带阻滤波器、多频带滤波器。b=fir1(M,Wc)b=fir1(M,Wc,'ftype')b=fir1(M,Wc,window)b=fir1(M,Wc,'ftype',window)b=fir

4、1(M,Wc,‘ftype',window)输出参数：b为FIR数字滤波器的M+1个系数构成的矩阵（即系统的单位脉冲响应）输入参数：M为FIR数字滤波器的阶数。Wc为3dB截频：0

5、不指定，默认为哈明窗。2.频率取样法设计FIR滤波器fir2函数可以实现FIR数字滤波器的频率取样法设计。可设计任意形状频率响应的滤波器。格式如下：b=fir2(M,f,m)b=fir2(M,f,m,window)输出参数：b为FIR数字滤波器的M+1个系数构成的矩阵。输入参数：M为滤波器的阶数。f指定归一化的各频带边界频率，从0到1递增，1对应于fsam/2，即数字频率W＝p。m指定各频带边界频率处的幅度响应，因此f和m的长度相等，即length(f)=length(a)。window指定窗函数，若不指定，默

6、认为哈明窗。四：实验步骤第一：实验内容1.分别使用矩形窗、汉明窗、汉宁窗设计一个阶数M=9的FIR数字低通滤波器，截频为（1）画出各种方法设计的数字滤波器的单位脉冲响应。（2）画出它们的幅频响应，并比较各滤波器的通带纹波和阻带纹波，有何结论?（3）若输入为计算各滤波器的输出并画出其波形。2.利用频率抽样方法设计FIR数字低通滤波器，并绘出衰耗特性。已知阶数M=15，给定指标为：改变Ad[4]的值，观察该FIR低通数字滤波器的衰耗特性的变化。3.利用频率抽样方法设计FIR数字带通滤波器，并绘出衰耗特性。已知阶数M

7、=15，给定指标为：改变Ad[2]或Ad[6]的值，观察该FIR带通数字滤波器的衰耗特性的变化。4.设计一窄带通FIR数字滤波器，通带中心频率带宽不大于。（1）利用fir1函数和kaiser窗设计该滤波器。（2）利用fir3函数设计该滤波器，达到fir1函数的设计效果。（3）分别画出上述两个滤波器的实现结构，并比较其经济性。第二：实验程序，波形：实验内容（1）%Rectangularwindowb1=fir1(9,1/pi,boxcar(10));[H1,w]=freqz(b1,1,512);H1_db=20*

8、log10(abs(H1));%hammingwindowb2=fir1(9,1/pi,hamming(10));[H2,w]=freqz(b2,1,512);H2_db=20*log10(abs(H2));%hanningwindowb3=fir1(9,1/pi,hanning(10));[H3,w]=freqz(b3,1,512);H3_db=20*log10(abs(H3));s