2022年用窗函数法设计FIR数字滤波器 .pdf

《2022年用窗函数法设计FIR数字滤波器 .pdf》由会员分享，可在线阅读，更多相关《2022年用窗函数法设计FIR数字滤波器 .pdf（5页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、用窗函数法设计FIR 数字滤波器实验六用窗函数法设计FIR 数字滤波器实验目的掌握用窗函数法设计FIR 数字滤波器的原理与方法. 熟悉线性相位FIR 数字滤波器的特性 . 了解各种窗函数对滤波特性的影响. 实验内容和要求复习用窗函数法设计FIR 数字滤波器一节内容,阅读本实验原理 ,掌握设计步骤 . 用升余弦窗设计一线性相位低通FIR 数字滤波器 ,截止频率 .窗口长度 N =15,33.要求在两种窗口长度情况下,分别求出 ,打印出相应的幅频特性和相频特性曲线,观察 3dB带宽和 20dB 带宽 .总结窗口长度N 对滤波器特性的影响 . 设计低通 FIR 数字滤波器时 ,一般以理想低通滤波特性

2、为逼近函数,即其中,用四种窗函数设计线性相位低通滤波器,绘制相应的幅频特性曲线,观察 3dB 带宽和20dB 带宽以及阻带最小衰减,比较四种窗函数对滤波器特性的影响. 实验方法 ,步骤及结果测试如果所希望的滤波器的理想的频率响应函数为,则其对应的单位脉冲响应为(4.1) 窗函数设计法的基本原理是用有限长单位脉冲响应序列逼近.由于往往是无限长序列,而且是非因果的 ,所以用窗函数将截断,并进行加权处理 ,得到: (4.2) 就作为实际设计的FIR 数字滤波器的单位脉冲响应序列,其频率响应函数为(4.3) 式中,N 为所选窗函数的长度. 我们知道 ,用窗函数法设计的滤波器性能取决于窗函数的类型及窗口

3、长度N 的取值 .设计过程中 ,要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度 N .各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见表4.1. 表 1 各种窗函数的基本参数窗函数旁瓣峰值幅度 /dB 过渡带宽阻带最小衰减 /dB 矩形窗-13 4/N-12 三角形窗-25 8/N-25 汉宁窗-31 8/N-44 哈明窗名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 5 页 - - - - - - - - - -41 8/N-53 不莱克曼窗-5

4、7 12/N-74 凯塞窗 (=7.865)-57 10/N-80 这样选定窗函数类型和长度N 之后,求出单位脉冲响应 ,并按照式 (4.3)求出.是否满足要求,要进行演算 .一般在尾部加零使长度满足2 的整数次幂 ,以便用 FFT 计算.如果要观察细节 ,补零点数增多即可 .如果不满足要求 ,则要重新选择窗函数类型和长度N ,再次验算 ,直至满足要求 . 如果要求线性相位特性,则还必须满足根据上式中的正 ,负号和长度N 的奇偶性又将线性相位FIR 滤波器分成四类 .要根据所设计的滤波特性正确选择其中一类,例如,要设计线性相位低通特性,可以选择这一类,而不能选择这一类 . 主程序框图如图4.1

5、 所示.其中幅度特性要求用dB 表示 . 设画图时 ,用打印幅度特性.第 k 点对应的频率.为使曲线包络更接近的幅度特性曲线,DFT 变换区间要选大些 .例如窗口长度N=33 时,可通过在末尾补零的方法,使长度变为 64,再进行 64 点 DFT,则可以得到更精确的幅度衰减特性曲线. 下面给出 MATLAB 主程序 : %实验四 ,用窗函数法设计FIR 数字滤波器b=1; close all; i=0; while(b); temp=menu(选择窗函数长度N,N=10,N=15,N=20,N=25,N=30,N=33,N=35,N=40,N=45,N=50,N=55,N=60,N=64);

6、menu1=10,15,20,25,30,33,35,40,45,50,55,60,64; N=menu1(temp); temp=menu(选择逼近理想低通滤波器截止频率Wc,Wc=pi/4,Wc=pi/2,Wc=3*pi/4,Wc=pi,Wc=0.5,Wc=1.0,Wc=1.5,Wc=2.0,Wc=2.5,Wc=3.0); menu2=pi/4,pi/2,3*pi/4,pi,0.5,1,1.5,2,2.5,3; w=menu2(temp); n=0:(N-1); hd=ideal(w,N); %得到理想低通滤波器k=menu(请选择窗口类型 :,boxcar,hamming,hanning

7、,blackman); if k=1 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 5 页 - - - - - - - - - B=boxcar(N); string=Boxcar,N=,num2str(N); else if k=2 B=hamming(N); string=Hamming,N=,num2str(N); else if k=3 B=hanning(N); string=Hanning,N=,num2str(N); else if k=4 B=black

8、man(N); string=Blackman,N=,num2str(N); end end end end h=hd.*(B); % 得到 FIR 数字滤波器H,m=freqz(h,1,1024,whole); % 求其频率响应mag=abs(H); %得到幅值db=20*log10(mag+eps)/max(mag); pha=angle(H); %得到相位i=i+1; figure(i) subplot(2,2,1); n=0:N-1; stem(n,h,.); axis(0,N-1,-0.1,0.3); hold on; n=0:N-1; x=zeros(N); plot(n,x,-)

9、; xlabel(n); ylabel(h(n); title(实际低通滤波器的h(n); text(0.3*N),0.27,string); hold off; subplot(2,2,2); plot(m/pi,db); axis(0,1,-100,0); xlabel(w/pi); ylabel(dB); title(衰减特性 (dB); grid; subplot(2,2,3); plot(m,pha); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 5 页 -

10、 - - - - - - - - hold on; n=0:7; x=zeros(8); plot(n,x,-); title(相频特性 ); xlabel(频率(rad); ylabel(相位(rad); axis(0,3.15,-4,4); subplot(2,2,4); plot(m,mag); title(频率特性 ); xlabel(频率 W(rad); ylabel(幅值); axis(0,3.15,0,1.5); text(0.9,1.2,string); b=menu(Do You want To Continue ,Yes,No); if b=2 b=0; end end t

11、emp=menu(Close All Figure ,Yes,No); if temp=1 close all end 程序运行结果 : 运行程序 ,根据实验内容要求和程序提示选择你要进行的实验参数.三个实验参数选定后,程序运行输出用所选窗函数设计的实际FIR 低通数字滤波器的单位脉冲响应h(n),幅 频 衰 减 特 性 (20lgH(ejw), 相 频 特性 及 幅 频 特 性 H(ejw) 的 波 形 ,h(n)和H(ejw) 图中标出了所选窗函数类型及其长度N 值.对四种窗函数 (N=15 和 N=33)的程序运行结果如图4-2 到图 4-9 所示,由图可以看出用各种窗函数设计的FIR

12、滤波器的阻带最小衰减及过渡带均与教材中一致.在通带内均为严格相位特性. 实验报告要求简述实验原理及目的. 按照实验步骤以及要求,比较各种情况下的滤波性能,说明窗口长度N 和窗函数类型对滤波特性的影响 . 总结用窗函数法设计FIR 滤波器的主要特点 . 简要回答思考题 . 思考题如果给定通带截止频率和阻带截止频率以及阻带最小衰减,如何用窗函数法设计线性相位低通滤波器 ,写出设计步骤 . 如果要求用窗函数法设计带通滤波器,而且给定上 ,下边带截止频率为和,试求理想带通的单位脉冲响应 %实验中的子函数 :产生理想低通滤波器单位脉冲响应hd(n) function hd=ideal(w,N); alp

13、ha=(N-1)/2; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 5 页 - - - - - - - - - n=0:(N-1); m=n-alpha+eps; hd=sin(w*m)./(pi*m); 相关练习 : 例 1 用凯塞窗设计一FIR 低通滤波器 ,低通边界频率,阻带边界频率,阻带衰减不小于 50dB. 解 首先由过渡带宽和阻带衰减来决定凯塞窗的N 和图 4.1给出了以上设计的频率特性,(a) 为N=30 直接截取的频率特性(b)为凯塞窗设计的频率特性

14、.凯塞窗设计对应的MATLAB程序为: wn=kaiser(30,4.55); nn=0:1:29; alfa=(30-1)/2; hd=sin(0.4*pi*(nn-alfa)./(pi*(nn-alfa); h=hd.*wn; h1,w1=freqz(h,1); plot(w1/pi,20*log10(abs(h1); axis(0,1,-80,10); grid; xlabel(归一化频率 /) ylabel(幅度/dB) 例 2 利用雷米兹交替算法 ,设计一个线性相位低通FIR 数字滤波器 ,其指标为 :通带边界频率fc=800Hz,阻带边界fr=1000Hz,通带波动阻带最小衰减At

15、=40dB,采样频率fs=4000Hz. 解在 MATLAB 中可以用 remezord 和 remez两个函数设计 ,其结果如图 4.2,MATLAB 程序如下 : fedge=800 1000; mval=1 0; dev=0.0559 0.01; fs=4000; N,fpts,mag,wt=remezord(fedge,mval,dev,fs); b=remez(N,fpts,mag,wt); h,w=freqz(b,1,256); plot(w*2000/pi,20*log10(abs(h); grid; xlabel(频率/Hz) ylabel(幅度/dB) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 5 页 - - - - - - - - -