matlab矩形窗设计FIR高通,基于窗函数法FIR数字滤波器的设计最终版（全文完整版）...

《基于窗函数法FIR数字滤波器的设计.doc》由会员分享，可免费在线阅读全文，更多与《基于窗函数法FIR数字滤波器的设计(最终版)》相关文档资源请在帮帮文库(www.woc88.com)数亿文档库存里搜索。

1、的窗谱为：TTTTTTW)sin()sin(sin)((式)由此式可以看出，汉宁窗可以看作是个矩形时间窗的频谱之和，或者说是个sin(t)型函数之和，而括号中的两项相对于第一个谱窗向左、右各移动了πT，从而使旁瓣互相抵消，消去高频干扰和漏能。可以看出，汉宁窗主瓣加宽并降低，旁瓣则显著减小，从减小泄漏观点出发，汉宁窗优于矩形窗。但汉宁窗主瓣加宽，相当于分析带宽加宽，频率分辨率下降。()海明(Hamming)窗海明窗也是余弦窗的一种，又称改进的升余弦窗，其时间函数表达式为：TtTtTtTt,),cos()((式)其窗谱为：TTTTTTW)sin()sin(sin)((式)海明窗与汉宁窗都是余弦窗，只是加权系数不同。海明窗加权的系数能使旁瓣达到更小。分析表明，海明窗的第一旁瓣衰减为dB。海明窗的频谱也是由个矩形窗的频谱合成，但其旁瓣衰减速度为dB／(oct)，这比汉宁窗衰减速度慢。海明窗与汉宁窗都是很有用的窗函数。()高斯窗高斯窗是一种指数窗。其时域函数为：TtTteTtat,,)((式)式中a为常数，决定了函数曲线衰减的快慢。a值如果选取适当，可以使截断点(T为有限值)处的函数值比较小，则截断造成的影响就比较小。高斯窗谱无负的旁瓣，第一旁瓣衰减达一dB。高斯窗的主瓣较宽，故而频率分辨率低。高斯窗函数常被用来截断一些非周期信号，如指数衰减信号等。不同的窗函数对信号频谱的影响是不一样的，这主要是因为不同的窗函数，产生泄漏的大小不一样，频率分辨能力也不一样。信号的截断产生了能量泄漏，而用FFT算法计算频谱又产生了栅栏效应，从原理上讲这两种误差都是不能消除的，但是我们可以通过选择不同的窗函数对它们的影响进行抑制。图是几种常用的窗函数的时域和频域波形，其中矩形窗主瓣窄，旁瓣大，频率识别精度最高，幅值识别精度最低；布莱克曼窗主瓣宽，旁瓣小，频率识别精度最低，但幅值识别精度最高。图几种常用的窗函数的时域和频域波形对于窗函数，还有一些要求：)。

2、满足设计要求。取N=,偶对称,得:过渡带宽Δω:π,第一通带波纹:dB,第一阻带最小衰减:dB。我们设计的低通FIR数字滤波器的理想性能指标是：通带截止频率W=π,阻带截止频率Ws=π，阻带衰减At不小于dB，通带衰减不大于dB,fs=。Matlab程序w=*i；ws=*i；wdelta=wsw；N=ceil(*iwdelta)；Wn=(+)*i；b=fir(N,Wni,hanning(N+))；freqz(b,,)N=ceil()；n=:N；window=hanning(N)；[h,w]=freqz(window,)；figure()；stem(window)；axis([])；grid；xlabel('n')；title('Hanning窗函数')；figure()；lot(wi,*log(abs(h)abs(h())))；axis([])；grid；xlabel('wi')；ylabel('幅度(db)')；title('Hanning窗函数的频谱')；hn=fir(N,wc,hanning(N))；[h,w]=freqz(hn,,)；figure()；stem(n,hn)；axis([])；grid；xlabel('n')；ylabel('h(n)')；title('Hanning窗函数的单位脉冲响应')；figure()；lot(wi,*log(abs(h)abs(h())))；grid；xlabel('wi')；ylabel('幅度(db)')；gtgt运行结果：nHanning窗函数wi幅度(db)Hanning窗函数的频谱nh(n)Hanning窗函数的单位脉冲响应滤波器，其它滤波器都不能设计；第四种情况，奇对称、偶数点，可设计高、带通滤波器，不能设计低通和带阻。FIR数字滤波器的设计原理一个截止频率为c(rads)的理想数字低通滤波器，其传递函数的表达式是：ccjjdeeH,,)((式)由式可以看出，这个滤波器在物理上是不可实现的，因为冲激响应具有无限性和因果性。为了产生有限长度的冲激响应函数，我们取样。

3、是频域抽样内插函数，其幅度函数)(RW在N之内为一个主瓣，两侧形成许多衰减振荡的旁瓣，如果将理想频率响应也写成)()()(NjdjdeHeH(式)则其幅度函数为ccdH,,)((式)窗函数介绍实际应用的窗函数，可分为以下主要类型：、幂窗采用时间变量某种幂次的函数，如矩形、三角形、梯形或其它时间(t)的高次幂；、三角函数窗应用三角函数，即正弦或余弦函数等组合成复合函数，例如汉宁窗、海明窗等；、指数窗采用指数时间函数，如ste形式，例如高斯窗等。下面介绍几种常用窗函数的性质和特点。()矩形窗矩形窗属于时间变量的零次幂窗，函数形式为:TtTtTt,,)((式)相应的窗谱为：TTWsin)((式)矩形窗使用最多，习惯上不加窗就是使信号通过了矩形窗。这种窗的优点是主瓣比较集中，缺点是旁瓣较高，并有负旁瓣，导致变换中带进了高频干扰和泄漏，甚至出现负谱现象。图矩形窗的时域及频域波形()三角窗亦称费杰(Fejer)窗，是幂窗的一次方形式，其函数形式是：TtTtTtTt,),()((式)三角窗与矩形窗比较，主瓣宽约等于矩形窗的两倍，但旁瓣小，而且无负旁瓣，如图所示。图三角窗的时域及频域波形()汉宁(Hanning)窗汉宁窗又称升余弦窗，其时域表达式为：TtTtTtTt,),cos()((式)相应延时响应。图各种理想数字滤波器的幅度频率响应窗函数设计法窗函数设计原理分析设数字滤波器的传输函数为)(jeH,)(nhd是与其对应的单位脉冲响应,)(zH为系统函数。低通高通带通带阻全通)(jdeH)(jdeH)(jdeH)(jdeH)(jdeH)()(NnnjjenheH(式)deeHnhnjjdd)()((式))()(NnnznhzH(式)一般说来,)(nhd是无限长的,需要求对)(jdeH的一个逼近。采用窗函数设计法时,可通过对理想滤波器的单位采样响应加窗设计滤波器)()()(nhnnhd(式)其中,。

4、相当于一个全通滤波器减去一个低通滤波器。、一个带通滤波器相当于两个低通滤波器相减。、一个带阻滤波器相当于一个低通滤波器加上一个高通滤波器。这里的相加相减都是相当于并联结构。图中所示的各种数字滤波器理想频率响应只表示了正频率部分，这样的理想频率响应是不可能实现的，原因是频带之间幅度响应是突变的，因而其单位抽样响应是非因果的。因此要给出实际逼近容限。数字滤波器的系统函数)(zH，它在z平面单位圆上的值为滤波器频率响应)(jeH，表征数字滤波器频率响应特征的三个参量是幅度平方响应、相位响应和群。利用()式可得)()(sin)(nndeenhcccnjjdcc(式))(nhd是中心点在的偶对称无限长非因果序列，要得到有限长的)(nh，一种最简单的方法就是取矩形窗)(nRN，即)()(nRnN但是按照线形相位滤波器的约束，)(nh必须是偶对称的，对称中心应为长度的一半(N)，因而必须=(N)，所以有N,),()()()(为其他nNnnhnnhnhdd(式)将(式)代入(式)，可得,,)()(sin)(NnnNnNnnhccc为其他值(式)此时，一定满足)()(nNhnh这一线性相位的条件。下面求)(nh的傅里叶变换，也就是找出待求FIR滤波器的频率特性，以便能看出加窗处理后究竟对频率响应有何影响。按照复卷积公式，在时域是相乘、频域上是周期性卷积关系，即deeHeHjjdj)()()((式)因而)(jeH逼近)(jdeH的好坏，完全取决于窗函数的频率特性)(jeW。窗函数)(n的频率特性)(jeW为)()(NnnjjeneW(式)对矩形窗)(nRN，则有)sin()sin()(NNeeeWNjNnnjjR(式)也可表示成幅度函数与相位函数)()()(NjRjNeWeW(式)其中)sin()sin()(NNWR(式))(jReW就。

5、dB带宽B，它是主瓣归一化的幅度下降到dB时的带宽。当数据长度为N时，矩形窗主瓣两个过零点之间的宽度为πN。)最大边瓣峰值A(dB)。)边瓣谱峰渐进衰减速度D(dBcot)。所以，理想的窗函数应当具有最小的B和A，和最大的D。窗函数法设计步骤、首先是给定所要求的频率响应函数)(jdeH；、其次，求单位冲激响应deeHnhnjjdd)()(；、再次，有过渡带宽及阻带最小衰减的要求，查表选定窗函数及N的大小，一般N的大小要通过几次试探而后确定；、求得所设计的FIR滤波器的单位冲激响应；)()()(nhnnhd，n=，，,N；、求)()(NnnjjenheH，检验是否满足设计要求，如不满足，则需要重新设计。设计实例线性相位FIR低通滤波器的设计(用窗函数法)。指标要求:通带截止频率:π,阻带起始频率:π,阻带最小衰减:dB。()设)(jeH为理想线性相位滤波器其他,)(,cjjeeHwi幅度(db)mKaiser窗函数wi幅度(db)Kaiser窗函数的频谱NormalizedFrequency(radsamle)Phase(degrees)NormalizedFrequency(radsamle)Magnitude(dB)由所需低通滤波器的过渡带求出理想低通滤波器的截止数字频率ω=π,得出:。为线性相位所需的移位,,,)],(sin[)()(Nnnnndedeenhccnjnjjdcc()由阻带衰减确定窗函数,由过渡带宽确定N值。阻带最小衰减dB,比对种窗函数基本参数选定窗函数为海明窗。所要求的过渡带宽:Δω=ππ=πN=ππ=,τ=(N)=()由海明窗函数确定FIR滤波器的h(n)。)()](sin[)()()]cos([nNnNnnhnRNncdN)(得出：)()]cos([)()(sin[)()()(nRNnnnnhnhNd()仿真检验各项指标,得出结论。

6、响应为)(nh，长度为N，其系数函数为)(zH：)()(NnnznhzH(式)用)(nh表示截取)(nhd后冲激响应，即)()()(nhnnhd，式子中)(n为窗函数，长度为N。当τ=(N)时，截取的一段)(nh对(N)对称，可保证所设计的滤波器具有线性相位。一般来说，FIR数字滤波器输出)(ny的Z变换形式Y(z)与输入)(nx的Z变换形式之间的关系如下：)())()()(()()()(zXznhzhhzXzHzYn(式)从上面的Z变换和结构图可以很容易得出FIR滤波器的差分方程表示形式。对式进行反Z变换，可得：)()()()()()()(xnhnxhnxhny(式)图卷积型滤波器式()为FIR数字滤波器的时域表示方法，其中)(nx是在时间n的滤波器的输入抽样值。根据式()即可对滤波器进行设计。从上面的公式我们可以看出，在对滤波器实际设计时，整个过程的运算量很大。设计完成后对已设计的滤波器的频率响应进行校核，运算量也很大。并且在数字滤波器设计的过程中，要根据设计要求和滤波效果不断地调整，以达到设计的最优化。在这种情况下，要进行大量复杂的运算，单纯靠公式计算和编制简单的程序很难在短时间内完成。而利用MATLAB工具进行计算机辅助设计，则可以快速有效地设计数字滤波器，大大的减少了计算量。数字滤波器的性能指标我们在进行滤波器设计时，需要确定其性能指标。一般来说，滤波器的性能要求往往以频率响应的幅度特性的允许误差来表征。以低通滤波器特性为例，频率响应有通带、过渡带及阻带三个范围。在通带内：AP≤)(jeH≤c≤c在阻带中：)(jeH≤stAst≤≤c其中c为通带截止频率,st为阻带截止频率，A为通带误差,stA为阻带误差。图低通滤波器的幅度特性与模拟滤波器类似，数字滤波器按频率特性划分为低通、高通、带通、带阻、全通等类型，由于数字滤波器的频率响应是周期性的，周期为π。由于频率响应的周期性，频率变量以数字频率来表示，所以数字滤波器设计中必须给出抽样频率。图为各种数字滤波器理想幅度，可以看出：、一个高通滤波。