fir数字滤波器设计及matlab实现,fir数字滤波器设计及matlab使用.doc

fir数字滤波器设计及matlab使用.doc

数字信号处理课程设计数字信号处理课程设计报告FIR数字滤波器设计及MATLAB实现专业通信工程班级通信1101班组次第9组姓名及学号姓名及学号数字信号处理课程设计第1页组员承担任务负责编写基于MATLAB的时域采样理论研究及实现源程序，进行仿真实验，调试、运行，截取实验图形，对实验结果进行分析。深入理解课题含义及设计要求，注意材料收集与整理，对课程设计要求进行最后审核。负责上网查阅或者去图书馆查阅有关课题的资料，并且撰写及修改课程设计的报告，对实验结论的总结，制作课程设计答辩的PPT。指导教师评价意见数字信号处理课程设计第2页目录一、设计目的3二、设计任务3三、设计原理331窗函数法332频率采样法433最优化设计5331等波纹切比雪夫逼近准则5332仿真函数6四、设计过程7五、收获与体会13参考文献13数字信号处理课程设计第3页FIR数字滤波器设计及MATLAB实现一、设计目的FIR滤波器有限长单位冲激响应滤波器，是数字信号处理系统中最基本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。滤波器设计是根据给定滤波器的频率特性，求得满足该特性的传输函数。二、设计任务FIR滤波器设计的任务是选择有限长度的，使传输函数满足一HNJWHE定的幅度特性和线性相位要求。由于FIR滤波器很容易实现严格的线性相位，所以FIR数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。设计过程一般包括以下三个基本问题(1)根据实际要求确定数字滤波器性能指标；(2)用一个因果稳定的系统函数去逼近这个理想性能指标；(3)用一个有限精度的运算去实现这个传输函数。三、设计原理FIR滤波器设计的任务是选择有限长度的，使传输函数满足一HNJWHE定的幅度特性和线性相位要求。由于FIR滤波器很容易实现严格的线性相位，所以FIR数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。设计过程一般包括以下三个基本问题(1)根据实际要求确定数字滤波器性能指标；(2)用一个因果稳定的系统函数去逼近这个理想性能指标；(3)用一个有限精度的运算去实现这个传输函数。31窗函数法设计FIR数字滤波器的最简单的方法是窗函数法，通常也称之为傅立叶级数数字信号处理课程设计第4页法。FIR数字滤波器的设计首先给出要求的理想滤波器的频率响应,设计JWDHE一个FIR数字滤波器频率响应，去逼近理想的滤波响应。然而，JWHEJ窗函数法设计FIR数字滤波器是在时域进行的，因而必须由理想的频率响应推导出对应的单位取样响应，再设计一个FIR数字滤波器的单位取JWDHEDHN样响应去逼近。设计过程如下HNDN(3IDTFWNDTFJJDEHE1)加窗的作用是通过把理想滤波器的无限长脉冲响应乘以窗函数来DHNWN产生一个被截断的脉冲响应，即并且对频率响应进行平滑。DHNWMATLAB工具箱提供的窗函数有矩形窗RECTANGULARWINDOW、三角窗TRIANGULARWINDOW、布拉克曼窗BLACKMANWINDOW、汉宁窗HANNINGWINDOW、海明窗HAMMINGWINDOW、凯塞窗KAISERWINDOW、切比雪夫窗CHEBYSHEVWINDOW。窗函数主要用来减少序列因截断而产生的GIBBS效应。但当这个窗函数为矩形时，得到的FIR滤波器幅频响应会有明显的GIBBS效应，并且任意增加窗函数的长度(即FIR滤波器的抽头数)GIBBS效应也不能得到改善。为了克服这种现象,窗函数应该使设计的滤波器(1)频率特性的主瓣宽度应尽量窄,且尽可能将能量集中在主瓣内；(2)窗函数频率特性的旁瓣Ω趋于Π的过程中，其能量迅速减小为零。32频率采样法频率采样法是从频域出发,根据频域采样定理，对给定的理想滤波器的频率响应加以等间隔的抽样，得到JWHEDHKK0,1,,N12KNJWDDKHE(32)再利用可求得FIR滤波器的系统函数及频率响应。DHZJWE而在各采样点间的频率响应则是其的加权内插函数延伸叠加的结果。但数字信号处理课程设计第5页对于一个无限长的序列，用频率采样法必然有一定的逼近误差,误差的大小取决于理想频响曲线的形状,理想频响特性变换越平缓,则内插函数值越接近理想值,误差越小。为了提高逼近的质量，可以通过在频率相应的过渡带内插入比较连续的采样点，扩展过渡带使其比较连续，从而使得通带和阻带之间变换比较缓慢，以达到减少逼近误差的目的。选取W∈0,2Π内N个采样点的约束条件为(3HKKMM01N3)(1)增大阻带衰减三种方法1)加宽过渡带宽，以牺牲过渡带换取阻带衰减的增加。2)过渡带的优化设计利用线性最优化的方法确定过渡带采样点的值，得到要求的滤波器的最佳逼近(而不是盲目地设定一个过渡带值)。3)增大N。如果要进一步增加阻带衰减，但又不增加过渡带宽，可增加采样点数N。代价是滤波器阶数增加，运算量增加。直接从频域进行设计，物理概念清楚，直观方便；适合于窄带滤波器设计，这时频率响应只有少数几个非零值，但是截止频率难以控制。典型应用用一串窄带滤波器组成多卜勒雷达接收机，覆盖不同的频段，多卜勒频偏可反映被测目标的运动速度；33最优化设计最优化设计方法是指采用最优化准则来设计的方法。在FIRDF的最优化设计中,最优化准则有均方误差最小化准则和等波纹切比雪夫逼近也称最大误差最小化准则两种。实际设计中,只有采用窗函数法中的矩形窗才能满足前一种最优化准则,但由于吉布斯GIBBS效应的存在，使其根本不能满足设计的要求。为了满足设计的要求,可以采用其它的窗函数来消除吉布斯效应,但此时的设计已经不能满足该最优化准则了。因此,要完成FIRDF的最优化设计,只能采用后一种优化准则来实现。数字信号处理课程设计第6页331等波纹切比雪夫逼近准则在滤波器的设计中,通常情况下通带和阻带的误差要求是不一样的。等波纹切比雪夫逼近准则就是通过对通带和阻带使用不同的加权函数,实现在不同频段通常指的是通带和阻带的加权误差最大值相同,从而实现其最大误差在满足性能指标的条件下达到最小值。尽管窗函数法与频率采样法在FIR数滤波器的设计中有着广泛的应用,但两者不是最优化的设计。通常线性相位滤波在不同的频带内逼近的最大容许误差要求不同。等波纹切比雪夫逼近准则就是通过通带和阻带使用不同的加权函数,实现在不同频段通常指的是通带和阻带的加权误差最大值相同,从而实现其最大误差在满足性能指标的条件下达到最小值，即使得和之JWDHEJWE间的最大绝对误差最小。等波纹切比雪夫逼近是采用加权逼近误差，它可以表示为JWEE34JWJWJJDEEWEH其中，为逼近误差加权函数在误差要求高的频段上，可以取较大的J加权值，否则，应当取较小的加权值。尽管按照FIR数字滤波器单位取样响应HN的对称性和N的奇、偶性，FIR数字滤波器可以分为4种类型，但滤波器的频率响应可以写成统一的形式3512JNWJKJWHEEH其中，K∈{0,1},HΩ为幅度函数，且是一个纯实数，表达式也可以写成统一的形式36JWDEQP其中，为Ω的固定函数，为M个