完整word版课程设计窗函数法设计FIR数字滤波器.docx
- 文档编号:8660219
- 上传时间:2023-02-01
- 格式:DOCX
- 页数:12
- 大小:269.13KB
完整word版课程设计窗函数法设计FIR数字滤波器.docx
《完整word版课程设计窗函数法设计FIR数字滤波器.docx》由会员分享,可在线阅读,更多相关《完整word版课程设计窗函数法设计FIR数字滤波器.docx(12页珍藏版)》请在冰豆网上搜索。
完整word版课程设计窗函数法设计FIR数字滤波器
课程设计(论文)任务书
院(系):
基层教学单位:
学号
学生姓名
专业(班级)
设计题目
利用窗函数法设计FIR数字滤波器
设
计
技
术
参
数
窗函数,通带截频,阻带截频,通带最大衰减,阻带最大衰减
幅度响应,相位响应等
设
计
要
求
(1)给定所要求的频率响应
(
);
(2)由
(
)的傅立叶反变换求出
(n);
(3)根据过渡带宽及阻带最小衰减的要求,选择窗的形状并估计窗口长度N;
(4)计算滤波器的单位脉冲响应h(n)=
(n)ω(n);
(5)用matlab编程实现,显示其幅频特性曲线,分析并得出结论。
工
作
量
共二周,每天工作8小时
工
作
计
划
第一到三天:
查资料并选定题目,完成设计方案
第三到七天:
根据题目编程并调试
第八到十四天:
完成论文
参
考
资
料
陈后金数字信号处理(2008年11月第二版)高等教育出版社
陈后金信后与系统(2007年12月1日第1版)
高会生MATLAB原理与工程应用(2006年1月第1版)电子工业出版社
丁玉美数字信号处理(2000年12月第二版)西安电子科技大学出版社
指导教师签字
基层教学单位主任签字
说明:
此表一式四份,学生、指导教师、基层教学单位、系部各一份。
年月日
课程设计评审意见表
指导教师评语:
成绩:
指导教师:
年月日
答辩小组评语:
成绩:
评阅人:
年月日
课程设计总成绩:
答辩小组成员签字:
年月日
利用窗函数法设计FIR数字滤波器
DesignFIRdigitalfilterswiththewindowfunction
摘要:
在数字信号处理中,数字滤波器是一种被广泛使用的信号处理部件。
分析FIR(有限冲激响应)数字滤波器的结构特征,得到了满足系统要求的数字滤波器设计法,结合实际工程所要求的数字滤波器指标,利用MATLAB对FIR数字滤波器进行了设计和仿真,最后通过讨论,给出窗函数法FIR滤波器的设计原则。
关键词:
FIR数字滤波器;窗函数;Matlab程序;设计
Abstract:
Thedigitalfilterisawidelyusedsignalprocessingsectionindigitalsignalprocessing,.AnalysisofthestructuralcharacteristicsoftheFIR(finiteimpulseresponse)digitalfilters,digitalfilterdesignmethodhasbeentomeetthesystemrequirements,combinedwithdigitalfilterindicatorsofpracticalengineeringrequired,usingMATLABtodesignandsimulationFIRdigitalfilter,andfinallygivenwindowfunctionFIRfilterdesignprinciplesthroughdiscussion.
Keywords:
FIRdigitalfilter;windowfunction;Matlabprogram;design
一.设计目的
1.掌握FIR数字滤波器的设计方法和步骤,深入了解其幅度和相位特性,并能在实际中学会选择与应用。
2.学习Matlab语言工具,并应用Matlab实现数字信号处理的各种方式,理解和熟悉数字信号处理的实际应用。
3.了解几种窗函数的性质,熟悉其波形图,分析加窗截断对波形的影响。
二.设计原理
1.FIR数字滤波器结构特征
FIR数字滤波器的冲激响应只能延续一定的时间,N阶因果有限冲激响应滤波器的差分方程表达式为:
其卷积和是一个有限和,可以直接计算
所以其涉及的基本操作就是简单的乘和加,其计算仅需要所给出的初始样本值及所要求的相关样本值。
通过精确设计,有限冲激响应滤波器在整个频率范围内均能提供精确的线性相位,由于系统的单位脉冲序列为有限
序列,当输入有限时,输出也必然为有限,这样其稳定性总可以独立于滤波器系数之外。
因此,在很多情况下,有限冲激响应滤波器成为首选,只要确定能满足要求的转移序列或者脉冲响应的常数,就可以准确地设计出满足要的FIR数字滤波器。
通常所采用的设计方法主要有窗函数、频率采样法和等波纹最佳逼近法等,其中窗函数法是从时域进行设计的,其算法简单、物理意义清晰,因此得到了较为广泛的应用。
2.算法研究
为了建立一个具有线性相位和稳定的非递归特性的有限脉冲响应滤波器即FIR滤波器,要考虑两个方面:
一是使用有限长的单位取样响应来逼近理想低通,二是单位取样响应对(N-1)/2对称,从而保证线性相位。
根据FIR数字滤波器输入输出关系的差分方程描述,有
上式反映了有限列长的单位取样响应。
滤波网络传输函数为:
利用傅立叶变换得单位取样响应
为:
是截止频率,
从而保证线性相位。
在窗函数法设计中使用有限列长的
逼近
,用Hamming窗函数
将
截断,进行加权处理后可得
。
此时频率响应函数为:
Hamming窗函数
]
,Blackman窗函数
。
其中
是幅度为1,长度为N的矩形序列。
3.窗函数法的原理
可以从时域或频域出发来设计FIR滤波器,从频域出发的方法称为频率采样设计法,从时域出发的方法称为窗函数法,也称为傅立叶级数法,这种方法应用较广泛。
窗函数法的设计思想是按照所要求的理想滤波器频率响应
设计一个FIR滤波器,使之频率响应
来逼近
。
因为设计是在时域中进行的,需先由
的傅立叶反变换导出序列
,即
由于
是矩形频率特性,有
是一无限长的序列,且是非因果的,而要设计的FIR滤波器的冲激响应序列是有限长的,所以要用有限长的序列
来逼近无限长的序列
,最有效的方法是截断
,或者说用一个有限长度的窗口函数
序列来截取
,即
按照复卷积公式,在时域中的乘积关系可表示成在频域中的周期性卷积关系,即可得所设计的FIR滤波器的频率响应
其中,
为截断窗函数的频率特性。
由此可见,实际的FIR数字滤波器的频率响应
逼近理想滤波器频率响应
的好坏,完全取决于窗函数的频率特性
。
4.窗函数法的设计步骤
(1)给定所要求的频率响应
;
(2)由
的傅立叶反变换求出
;
(3)根据过渡带宽及阻带最小衰减的要求,选择窗的形状并估计窗口长度N;
(4)计算滤波器的单位脉冲响应
(5)验算技术指标是否满足要求。
窗函数设计法设计时要满足以下两个条件:
(1)窗谱主辨尽可能窄,以获得较陡的过渡带;
(2)尽量减少窗谱的最大旁瓣的相对幅度,使能量尽量集中于主辨,进而增加阻带的衰减。
给定的滤波器指标参数一般为通带截止频率
、阻截止频率
、实际通带波动
和最小阻带衰减
。
窗函数设计的经验公式为
归一化过渡带
滤波器阶数
当
≥50时,
;当21≤
≤50时,
在实际工程中常用的函数有六种,海明窗(Hamming)、矩形窗(Boxcar)、三角窗(Bartlett)、汉宁窗(Hanning)、布拉克曼窗(Blackman)和凯泽窗(Kaiser)等。
三.设计内容
1.低通滤波器的设计指标。
通带截止频率
=0.2,阻带截止频率
=0.3,最小阻带衰减
=50dB。
2.设计分析
由以上参数可得过度带宽tr_width=ws-wp,wc=(ws+wp)/2。
考虑选择Hamming窗函数,则列长
N=8*pi/tr-width,
hd(n)=sin(wc*(n-(N-1)/2))/(pi*(n-(N-1)/2)),
h(n)=hd(n)*wh(n)。
利用MATLAB编程时,先计算N、tr-width、hd(n)和h(n),并利用MATLAB提供的相关工具计算该低通滤波器的幅度特性值,最终显示Hamming窗频谱、实际单位取样响应频谱、理想单位取样响应频谱和该FIR滤波器幅度特性曲线。
3.编写Matlab程序
(1)实现FIR滤波器的主程序:
%dfilter.m文件
%clear
wp=0.2*pi;ws=0.3*pi;%给出通带和阻带的截止频率
tr-width=ws-wp%过度带宽度
N=ceil(8*pi/tr_width)+1%计算单位取样响应列长,使其偶对称
n=[0:
1:
N-1];
wc=(ws+wp)/2%计算低通截止频率
hd=ideal_lp(wc,N);%计算理想低通单位取样响应
w_ham=(hamming(N))’;
h=hd.*w_ham;%应用窗函数加权后的有限长单位取样响应序列
[db,mag,pha,grd,w]=freqz_m(h,[1]);
%计算幅值响应、相位响应和群延迟响应
delta_w=2*pi/1000;
Rp=-(min(db(1:
1:
wp/delta_w+1)))
As=-(round(max(db(ws/delta_w+1:
1:
501))))
figure
(1);
subplot(2,2,2);stem(n,w_ham)
title(’HammingWindow’)
axis([0N01.1]);ylabel(’w(n)’)%显示Hamming窗频谱
subplot(2,2,3);stem(n,h)
title(’ActualImpulseResponse’)
axis([0N-0.10.3]);ylabel(’h(n)’)%显示实际单位取样响应频谱
subplot(2,2,4);plot(w/pi,db)
title(’MagnitudeResponseindb’);grid
set(gca,’XTickMode’,’manual’,’XTick’,[0,0.2,0.3,1])
axis([0,1,-100,20]);ylabel(’Decibels’)%显示FIR滤波器幅度特性曲线
subplot(2,2,1);stem(n,hd)
title(’IdealImpulseResponse’)
axis([0N-0.10.3]);ylabel(’hd(n)’)%显示理想单位取样响应频谱
(2)计算理想低通取样响应的程序:
%ideal_lp.m文件
functionhd=ideal_lp(wc,N);
%IdealLowpassfiltercomputation
%----------------------
%[hd]=ideal_lp(wc,M)
%hd=idealimpulseresponsebetween0toM-1
%wc=cutofffrequencyinredians
%M=lengthoftheidealfilter
%
alpha=(N-1)/2;
n=[0:
1:
N-1];
m=(n-alpha+eps);
hd=sin(wc*m)./(pi*m);
end
(3)计算幅值响应、相位响应和群延迟响应
%freqz_m.m文件
function[db,mag,pha,grd,w]=freqz_m(b,a);
%Modifiedversionoffreqzsubroutine
%--------------------------
%[db,mag,pha,grd,w]=freqz_m(b,a)
%db=RelativemagnitudeindBcomputedover0topiradians
%mag=absolutemagnitudecomputedover0topiradians
%pha=Phaseresponseinradiansover0topiradians
%grd=Groupdelayover0topiradians
%w=501frequencysamplesbetween0andpiradians
%b=nuberratorpolynomialofH(z)(forFIR:
b=h)
%a=denominatorpolynomialofH(z)(fora=[1]
%
[H,w]=freqz(b,a,1000,’whole’);
H=(H(1:
501));w=(w(1:
501))’;
mag=abs(H);
db=20*log10(mag+eps);
pha=angle(H);
grd=grpdelay(b,a,w);
end
运行后,有关参数:
N=81,tr_width=0.3142,wc=0.7854,As=53,相应幅频特性曲线为图1。
当改变参数N=201时,其幅频特性曲线为图2。
使用Blackman窗函数设计该滤波器,可将主程序中的w-ham=(hamming(N))′;改为w_ham=(Blackman(N))’;,N=81,相应幅频特性曲线为图3(图1~4中的子图依次为理想取样响应、窗函数、实际取样响应频谱和低通幅度特性曲线)。
图
(1)HammingN=81
四.结果分析:
在前面的FIR滤波器设计中,一是利用Hamming窗函数并且使用了不同的有限列长值(N),如图1和图2,二是利用不同的窗函数(Hamming窗和Blackman窗)使用同样的N值(N=81),如图1和图3。
通过分析图1和图2,注意到相同的窗函数下,N值越大,主瓣越窄,过度带越窄(在图中就是
=0.2处的下降沿的宽度),阻带的波动频率加快,但阻带的最小衰减没有改变(均接近-60dB),即最大的副瓣相对于主瓣的值没有改变。
虽然提高N值没有加大阻带衰减,但却使过度带变窄了。
再来分析图1和图3。
它们一个用的是Hamming窗,一个用的是Blackman窗,N值都是81。
使用Blackman窗函数时,过度带要比使用Hamming窗函数宽,但其阻带衰减却明显加大了(接近-80dB),其好处是减小了吉布斯效应,使阻带的震荡较小,通带也较平稳。
通过以上分析,可以得到利用窗函数设计FIR滤波器的具有实际指导意义的原则:
一是选择主瓣较宽的窗函数,从而加大阻带的衰减,保证通带的平稳,二是在保证阻带最小衰减指标的情况下,适当增加列长N值,使过度带窄一些。
但N值的变化同时影响过度带和主瓣的变宽或变窄,另外选择主瓣较宽的窗函数时(同等N值下),其过度带就宽一些了,可见这两个设计原则是无法同时满足的,这一结论可从上面的设计示例中得到验证。
因此设计FIR滤波器时,应根据技术指标,通过多次试验,找到合适的窗函数和N值。
与其他高级语言的程序设计相比,在MATLAB环境下利用程序法调用窗函数,可以更
方便、快捷地设计出具有严格线性相位的FIR滤波器,节省时间,提高编程效率,且参数的修改方便,还可以进一步进行优化设计。
随着信息技术的不断发展,利用MAT2LAB进行数字滤波器设计将发挥更大的作用。
同时,用MATLAB设计有关数字滤波器的参数,将会为硬件实现(比如说在DSP器件上的实现)提供一条简单而准确的途径和依据。
六.参考文献:
1.陈后金数字信号处理(2008年11月第二版)高等教育出版社
2.陈后金信后与系统(2007年12月1日第1版)
3.高会生MATLAB原理与工程应用(2006年1月第1版)电子工业出版社
4.丁玉美数字信号处理(2000年12月第二版)西安电子科技大学出版社
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 完整 word 课程设计 函数 设计 FIR 数字滤波器