《FIR带通滤波器的设计》由会员分享，可在线阅读，更多相关《FIR带通滤波器的设计（18页珍藏版）》请在金锄头文库上搜索。
1、目录1技术要求 12 基本原理 12.1 FIR 带通滤波器简介11.2 窗函数法原理 33 建立模型描述 33.1 MATLAB常用函数33.1.1 窗函数33.1.2 fir1 函数 43.1.3 freqz 函数43.14 ceil 函数53.1.5 其他函数与命令53.2 程序流程图54 源程序代码（含注释）74.1 矩形窗74.2 凯泽窗74.3 布拉克曼窗84.4 海明窗95 调试过程及结论105.1 程序运行结果105.2 实验结果分析126 心得体会137 思考题138 参考文献14FIR 带通滤波器的设计1 技术要求用窗函数法设计FIR带通滤波器。要求低端阻带截止频率3 1s
2、=0.2 n，低端通带截止频率 3 , =0.35 n ,高端通带截止频率3 =0.65 n ,高端阻带截止频率3 =0.8 n。绘出h（n）及其 丄Pu pu p幅频响应特性曲线。2 基本原理2.1FIR 带通滤波器简介带通滤波器是从滤波器的特性上划分的，带通滤波器是指能通过某一频率范围内的频 率分量、但将其他范围的频率分量衰减到极低水平的滤波器，与带阻滤波器的概念相对。 从实现的网络结构或者从单位脉冲响应长度分类，可以分为无限长单位脉冲响应（IIR） 滤波器和有限长单位脉冲响应（FIR）滤波器。IIR数字滤波器设计方法是利用模拟滤波器 成熟的理论及设计图表进行设计的，因而保留了一些经典模拟
3、滤波器优良的幅度特性。但 设计中只考虑了幅度特性，没考虑相位特性，所设计的滤波器一般是某种确定的非线性相 位特性。为了得到线性相位特性，对IIR滤波器必须另外增加相位相校正网络，是滤波器 设计变得复杂，成本也高，又难以得到严格的线性相位特性。FIR滤波器在保证幅度特性 满足技术要求的同时，很容易做到有严格的线性相位特性。两者各有优点，择其而取之。 后面的FIR滤波器的设计中，为获得有限长单位取样响应，需要用窗函数截断无限长单位 取样响应序列。一个理想的滤波器应该有一个完全平坦的通带，例如在通带内没有增益或者衰减，并 且在通带之外所有频率都被完全衰减掉，另外，通带外的转换在极小的频率范围完成。实
4、 际上，并不存在理想的带通滤波器。滤波器并不能够将期望频率范围外的所有频率完全衰 减掉，尤其是在所要的通带外还有一个被衰减但是没有被隔离的范围。这通常称为滤波器 的滚降现象，并且使用每十倍频的衰减幅度dB来表示。通常，滤波器的设计尽量保证滚 降范围越窄越好，这样滤波器的性能就与设计更加接近。然而，随着滚降范围越来越小，通带就变得不再平坦开始出现“波纹”。这种现象在通带的边缘处尤其明显，这种效应称为吉布斯现象。下图1为理想模拟带通滤波器幅频特性：FIR滤波器的单位脉冲响应h (n)是有限长的(OWnWN-1),其z变换为z -1的(N-1) 阶多项式：H (z)二 Y (Z)= 乂 h (n)
5、z -nX ( Z ) = 0可得FIR滤波器的系统差分方程为：y (n) = b(0) x(n) + b(1) x(n 一 1) + b( N 一 1) x(n 一 N +1)=七 b(m)x(n 一 m) - b(n) x(n)m=0因此，FIR滤波器又称为卷积滤波器。FIR滤波器的频率响应表达式为：H (e j ) = 1 h (n) e - j nn-0信号通过FIR滤波器不失真条件是在通带内具有恒定的幅频特性和线性相位特性。理 论上可以证明：当FIR滤波器的系数满足下列中心对称条件：h (n) h (N 1 n)h (n) h (N 1 n)或者时，滤波器设计在逼近平直幅频特性的同时
6、，还能获得严格的线性相位特性。线性相位FIR 滤波器的相位滞后和群延迟在整个频带上是相等且不变的。对于一个N阶的线性相位FIR 滤波器，群延迟为常数，即滤波后的信号简单地延迟常数个时间步长。这一特性使通带频 率内信号通过滤波器后仍保持原有波形形状而无相位失真。1.2 窗函数法原理设计FIR数字滤波器的方法通常有三种：窗函数法，频率抽样法，等纹波逼近法。本次课程设计讨论的是第一种窗函数法。这种方法也叫傅里叶级数法。一般是先给定所要求的理想滤波器频率响应H (ejw)，导出h (n)，我们知道理想滤波 dd器的冲击响应h (n)是无限长的非因果序列，而我们要设计的是h(n)是有限长的FIR滤波 d
7、器，所以要用有限长序列h(n)来逼近无限长序列h (n)，设：dh (n) = J H(ejw)e/wdw d2 兀-Ka常用的方法是使用有限长的窗函数w(n)来截取h (n)即：dh(n)=w(n) h (n)d根据在时域是相乘关系，在频域则是卷积关系：1H (ejw) =JK H (ejw)W ej(w-0)d62兀-k dR其中W (ejw)为矩形窗谱，H (ejw)是FIR滤波器频率响应。R在设计过程中，将无限长序列变为有限长序列是通过时域加矩形窗乘积来实现的截断。常见的窗函数有：矩形窗，汉宁窗，海明窗，布拉克曼窗，凯泽窗等。 3 建立模型描述3.1MATLAB 常用函数3.1.1 窗
8、函数几种常见的窗函数如表 1 所示：表 1 几种常见的窗函数的函数表示定义式窗函数矩形窗(Boxcar)w(n) = RN(n)= 1 0 n N -10，其他三角窗(Triang)w(k)= 2k., n+1,1 k n +122(n 一 k +1) n +1, k nn +12海明窗(Hamming)w(k) = 0.54 一 0.46cos 2n0 n M-l汉宁窗(Hamming)( kw(k) = 0.5 1 一 cos 2nI In +1丿丿0 n M-l巴特利特窗（Bartlett）w(k)= 4,1 k n1n 一 1222(k 1) n +1 k 2 一, k nn 12名称
9、旁瓣峰值/dB近似过渡带宽精确过渡带宽最小阻带衰减/dB矩形窗-134n/N1.8n/N-21海明窗-418n/N6.6n/N-53布拉克曼窗-5712n/Nlln/N-74凯泽窗（卩=7.865）-5710n/N-80表 2 几种常见的窗函数的基本参数3.1.2 fir1 函数设计标准响应FIR滤波器可使用firl函数。fir1函数以经典方法实现加窗线性相位FIR 滤波器设计，它可以设计出标准的低通，带通，高通和带阻滤波器。形式为：b=fir1 (n,Wc,ftype,Window)各个参数的含义如下：b滤波器系数。对于一个n阶的FIR滤波器，其n+1个滤波器系数可表示为:b(z)= b(1
10、)+ b(2)z-1 + b (n + 1)z 一 nn滤波器阶数；Wc截止频率，0Wc1, Wc=1对应于采样频率的一半。当设计带通滤波器时，Wc=Wc1Wc2,Wc1ro clear;%清除工作空间close all;%关闭所有打开的窗口wls=0.2*pi;wlp=0.35*pi;%参数设置whp=0.65*pi;whs=0.8*pi;delta_w=min(wlp-wls),(whs-whp);%求两个过渡带的较小者wc1=(wls+wlp)/2;wc2=(whp+whs)/2;%截止频率取通带阻带边界频率的均值%矩形窗N1=ceil(1.8*pi/delta_w);%根据矩形窗精确过
11、渡带宽i.8n/N计算窗宽subplot(2,1,1);%把窗口分割成 2 行1 列n=0:N1-1;stem(n,hn1,.);%绘制矩形窗的单位脉冲响应axis(0,N1-1,-0.4,0.4);%设置显示范围xlabel(n);ylabel(h(n);grid on;%确定 x，y 轴坐标名称，加网格title（矩形窗单位冲击响应h（n）;%添加图形的标题hn1=fir1(N1-1,wc1,wc2/pi,boxcar(N1); %检验设计的滤波器单位脉冲响应 h1,w1=freqz(hn1,1);figure(1)%建立图形窗口subplot(2,1,2);plot(w1/pi,20*log10(abs(h1);%绘制矩形窗的幅频特性曲线axis(0,1,-150,5);%设置显示范围xlabel（归一化角频率）;%确定x坐标ylabel（幅度（单位:分贝）;%确定 y 坐标grid on;title（矩形窗幅频
