数字信号处理用窗函数法设计FIR数字滤波器第三次实验.docVIP

数字信号实验报告 实验项目名称：用窗函数法设计FIR数字滤波器 所属课程名称： 数字信号处理 实 验 类 型 ： 综合型 指 导 教 师 ： 实 验 日 期 ： 2013.12.20 班 级 ： 学 号 ： 姓 名 ： 目录 一、实验目的································1 二、实验原理与方法······················1 三、实验内容与步骤·····················1 四、上级内容·······························2 1.程序 2.运行截图 3.注释 五、思考题··································8 PAGE 1用窗函数发设计 PAGE 1 PAGE PAGE 3 实验目的： 掌握用窗函数法设计FIR数字滤波器的原理和方法。 熟悉线性相位FIR数字滤波器特性。 了解各种窗函数对滤波特性的影响。 实验原理与方法： 如果所希望的滤波器的理想频率响应函数为 Hd(e jω)， 则其对应的单位脉冲响应为 用窗函数w(n)将hd(n)截断， 并进行加权处理， 得到: h(n)就作为实际设计的FIR数字滤波器的单位脉冲响应序列， 其频率响应函数H(e jω)为 如果要求线性相位特性， 则h(n)还必须满足： 根据上式中的正、 负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。 要根据所设计的滤波特性正确选择其中一类。 例如， 要设计线性相位低通特性， 可选择h(n)=h(N-1-n)一类， 而不能选h(n)=-h(N-1-n)一类。 三、实验内容与步骤 (1) 复习用窗函数法设计FIR数字滤波器一节内容， 阅读本实验原理， 掌握设计步骤。 (2) 编写程序。 ① 编写能产生矩型窗、 升余弦窗、 改进升余弦窗和二阶升余弦窗的窗函数子程序。 ② 编写主程序。 主程序框图如图 10.6.1 所示， 仅供参考。 其中幅度特性要求用dB表示。 用窗函数法设计滤波器主程序框图 设 四、上机实验内容 a、用升余弦窗设计一线性相位低通FIR数字滤波器，截止频率Wc= rad 。窗口N=15，33。要求在两种窗口长度情况下，分别求出h(n)，打印出相应的幅频特性和相频曲线。观察3db带宽和20db带宽，总结窗口N对滤波特性的影响。 设计低通FIR数字滤波器时，一般以理想低通滤波特性为逼近函数 即 其中 其中 close all; N=input(所有取值(10,15,20,25,30,33,35,40,45,50,55,60,64)\nN=); n=[0:(N-1)]; w=input(所有取值(pi/4,pi/2,3*pi/4,pi,0.5,1,1.5,2,2.5,3)\nw=); m=n-(N-1)/2+eps; hd=sin(w*m)./(pi*m); %得到理想低通滤波器 k=input(所有取值(boxcar=1,hamming=2,hanning=3,blackman=4)\nk=); if k==1 B=boxcar(N);%矩形窗 string=[Boxcar,N=,num2str(N)]; elseif k==2 B=hamming(N); string=[Hamming,N=,num2str(N)]; elseif k==3 B=hanning(N); string=[Hanning,N=,num2str(N)]; elseif k==4 B=blackman(N); string=[Blackman,N=,num2str(N)]; end h=hd.*(B); %得到FIR数字滤波器 [H,m]=freqz(h,[1],1024,whole); %求其频率响应 mag=abs(H); %得到幅值 db=20*log10((mag+eps)/max(mag)); %db计算公式 pha=angle(H); %得到相位 figure(1) subplot(2,2,1); n=0:N-1; stem(n,h,.); ax