《FIR数字滤波器设计及MATLAB实现》由会员分享，可在线阅读，更多相关《FIR数字滤波器设计及MATLAB实现（15页珍藏版）》请在金锄头文库上搜索。
1、?数字信号处理?课程设计报告FIR数字滤波器设计及MATLAB实现专 业： 通信工程 班 级： 通信09-1BF 组 次： 第16组 姓 名： 李 涛 学 号： FIR数字滤波器设计及MATLAB实现摘要：本文采用窗函数法和等波纹最正确逼近法两种方法分别设计一个FIR数字带通滤波器，运用MATLAB 进行FIR数字滤波器的设计与仿真，并比拟了两种方法的特点。结果说明，在同样的设计指标下，利用等波纹切比雪夫逼近法那么的设计可以获得最正确的频率特性和衰耗特性，具有通带和阻带平坦，过渡带窄等优点。关键词：FIR滤波器；MATLAB；窗函数法；等波纹最正确逼近法1 FIR滤波器的根本概念FIR滤波器：
2、有限长单位冲激响应滤波器，是数字信号处理系统中最根本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。滤波器设计是根据给定滤波器的频率特性，求得满足该特性的传输函数。2 FIR滤波器的特点有限长单位冲激响应FIR滤波器有以下特点： 1 系统的单位冲激响应在有限个n值处不为零； 2 系统函数在处收敛，极点全部在z = 0处因果系统； 3 结构上主要是非递归结构，没有输出到输入的反应，但有些结构中例如频率抽样结构也包含有反应的递归局部。 设FIR滤波器的单位冲激响应
3、为一个N点序列，那么滤波器的系统函数为 2-1就是说，它有N1阶极点在z = 0处，有N1个零点位于有限z平面的任何位置。优点 ：1很容易获得严格的线性相位，防止被处理的信号产生相位失真，这一特点在宽频带信号处理、阵列信号处理、数据传输等系统中非常重要；2可得到多带幅频特性；3极点全部在原点永远稳定，无稳定性问题；4任何一个非因果的有限长序列，总可以通过一定的延时，转变为因果序列，所以因果性总是满足；5无反应运算，运算误差小。缺点：1因为无极点，要获得好的过渡带特性，需以较高的阶数为代价；2无法利用模拟滤波器的设计结果，一般无解析设计公式，要借助计算机辅助设计程序完成。3 FIR 数字滤波器设
4、计及matlab实现 FIR滤波器设计的任务是选择有限长度的，使传输函数满足一定的幅度特性和线性相位要求。由于FIR 滤波器很容易实现严格的线性相位，所以FIR 数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。设计过程一般包括以下三个根本问题： 1 根据实际要求确定数字滤波器性能指标；2 用一个因果稳定的系统函数去逼近这个理想性能指标；3 用一个有限精度的运算去实现这个传输函数。3.1窗函数法3.1.1窗函数法设计原理设计FIR数字滤波器的最简单的方法是窗函数法，通常也称之为傅立叶级数法。FIR数字滤波器的设计首先给出要求的理想滤波器的频率响应,设计一个FIR数字滤波器频率响
5、应，去逼近理想的滤波响应。然而，窗函数法设计FIR数字滤波器是在时域进行的，因而必须由理想的频率响应推导出对应的单位取样响应，再设计一个FIR数字滤波器的单位取样响应去逼近。设计过程如下： 3-1加窗的作用是通过把理想滤波器的无限长脉冲响应乘以窗函数来产生一个被截断的脉冲响应，即并且对频率响应进行平滑。窗函数主要用来减少序列因截断而产生的Gibbs效应。但当这个窗函数为矩形时，得到的FIR滤波器幅频响应会有明显的Gibbs效应，并且任意增加窗函数的长度即FIR滤波器的抽头数Gibbs效应也不能得到改善。为了克服这种现象,窗函数应该使设计的滤波器：1 频率特性的主瓣宽度应尽量窄,且尽可能将能量集
6、中在主瓣内；2 窗函数频率特性的旁瓣趋于 的过程中，其能量迅速减小为零。实际工程中常用窗的特性及MATLAB函数比拟如表1所示。表1 常用窗函数性能比拟窗类型最小阻带衰减主瓣宽度精确过渡带宽窗函数矩形窗21dB4/M1.8/Mboxcar三角窗25dB8/M6.1/Mbartlett汉宁窗44dB8/M6.2/Mhanning哈明窗53dB8/M6.6/Mhamming布莱克曼窗74dB12/M11/Mblackman取凯塞窗时用kaiserord函数来得到长度M和kaiser3.1.2窗函数法设计FIR数字滤波器的步骤(1) 选择窗函数类型和长度，写出窗函数w(n)表达式。根据阻带最小衰减选
7、择窗函数w(n)的类型，再根据过渡带宽度确定所选窗函数的长度N。(2) 构造希望逼近的频率响应函数Hd(ej)。(3) 计算： 或者可以直接写出低通单位脉冲响应(4) 加窗得到设计结果：h(n)=hd(n)w(n)5检验滤波器性能： 本文将针对一个含有5Hz，15Hz和30Hz的混合正弦波信号，用设计的带通滤波器对其进行滤波，检验滤波器性能。3.1.3 设计题目及程序的代码 用窗函数法设计第二类线性相位FIR带通滤波器。要求阻带下截止频率Ws1=0.2pi，通带下截止频率Wp1=0.35pi，通带上截止频率Wp2=0.65pi，阻带上截止频率Ws2=0.8pi，通带最大衰减Rp=1dB，阻带最
8、小衰减As=60dB。采样频率fs=150Hz根据窗函数最小阻带衰减的特性。以及关参照窗函数的根本参数表，选择布莱克曼窗可到达75dB最小阻带衰减，其过渡带为11pi/N。本文的带通滤波器设计及滤波程序如下：clear;close all;clc; %输入设计指标Wp1=0.35*pi;Wp2=0.65*pi;Ws1=0.2*pi;Ws2=0.8*pi;As=60;Rp=1;Fs=150;Bt=min(Wp1-Ws1),(Ws2-Wp2);%求过渡带宽M=ceil(11*pi/Bt)+1; %求滤波器长度if mod(M,2)=0; N=M, else N=M+1, end; %选用第二类滤波
9、器n=0:1:N-1;m=0:1:N;wc1=(Ws1+Wp1)/2/pi;%理想带通滤波器的下截止频率wc2=(Ws2+Wp2)/2/pi; %理想带通滤波器的上截止频率wc=wc1,wc2;%设置理想带通截止频率w_bla=(blackman(N);%布莱克曼窗hn=fir1(N,wc,blackman(N+1);%调用fir1计算带通FIR数字滤波器的hnt=(0:100)/Fs;s=sin(2*pi*t*5)+sin(2*pi*t*15)+sin(2*pi*t*30);%混和正弦波信号sf=filter(hn,1,s);%对信号s进行滤波%滤波器幅频特性和相频特性figure(1)fr
10、eqz(hn,1,512);grid on;% 绘制结果并加网络%布莱克曼窗w(n)figure(2)stem(n,w_bla);title(blackman window) axis(0 N-1 0 1.1);xlabel(n);ylabel(w(n);%绘制滤波前后波形figure(3)subplot(2,1,1) plot(t,s);grid on;subplot(2,1,2)plot(t,sf);grid on;%实际单位脉冲响应hnfigure(4)stem(m,hn);title(实际单位脉冲响应h(n)xlabel(m);ylabel(h(n)3.1.4 设计并用matlab仿真
11、后的结果N = 76图1 滤波器幅频特性和相频特性 图2 滤波效果图滤波前后的波形 图3 布莱克曼窗函数波形图 图4 实际单位脉冲响应31.5 结果分析经计算N为76，设计的是第二类滤波器，一般来说选取的N越大，误差越小，但同时h(n)长度也增加，使本钱和计算量增大，因此在满足指标的情况下N只取76。另外从频率特性图可以看出，阻带衰减刚好满足指标，通带衰减有较大的充裕。从相位特性图可以看出，在截止频率内，显现出严格的线性特性。3.2 等波纹最正确逼近法等波纹最正确逼近法是指采用最优化准那么来设计的方法。在 FIR DF的最优化设计中 ,最优化准那么有均方误差最小化准那么和等波纹切比雪夫逼近(也
12、称最大误差最小化)准那么两种。实际设计中 ,只有采用窗函数法中的矩形窗 才能满足前一种最优化准那么 ,但由于吉布斯 (Gibbs )效应的存在，使其根本不能满足设计的要求。为了满足设计的要求 ,可以采用其它的窗函数来消除吉布斯效应 ,但此时的设计已经不能满足该最优化准那么了。因此 ,要完成 FIR DF的最优化设计 ,只能采用后一种优化准那么来实现。 等波纹最正确逼近法设计原理通常线性相位滤波在不同的频带内逼近的最大容许误差要求不同。等波纹切比雪夫逼近准那么就是通过通带和阻带使用不同的加权函数,实现在不同频段(通常指的是通带和阻带) 的加权误差最大值相同,从而实现其最大误差在满足性能指标的条件
13、下到达最小值，即使得 和之间的最大绝对误差最小。 等波纹切比雪夫逼近是采用加权逼近误差，它可以表示为： (3-4)其中，为逼近误差加权函数在误差要求高的频段上，可以取较大的加权值，否那么，应当取较小的加权值。 尽管按照 FIR 数字滤波器单位取样响应 h(n)的对称性和 N的奇、偶性，FIR 数字滤波器可以分为 4 种类型，但滤波器的频率响应可以写成统一的形式： (3-5)其中，k0 ,1 , H ()为幅度函数，且是一个纯实数，表达式也可以写成统一的形式： (3-6)其中，为的固定函数，为M个余弦函数的线性组合。设计题目及程序的代码用最优设计法设计线性相位FIR带通滤波器。要求阻带下截止频率
14、Ws1=0.2pi，通带下截止频率Wp1=0.35pi，通带上截止频率Wp2=0.65pi，阻带上截止频率Ws2=0.8pi，通带最大衰减Rp=1dB，阻带最小衰减As=60dB。采样频率fs=150Hz本文将针对一个含有5Hz，15Hz和30Hz的混合正弦波信号，用设计的带通滤波器对其进行滤波，检验滤波器性能。本文的带通滤波器设计及滤波程序如下：clear all,close all,clc;%输入设计指标Wp1=0.35*pi;Wp2=0.65*pi;Ws1=0.2*pi;Ws2=0.8*pi;f=0.2 0.35 0.65 0.8;m=0 1 0;rp=1;rs=60;Fs=150;%求波纹幅度dat1、dat
