FIR数字滤波器设计666.docVIP

《数字信号处理》 课程设计报告 FIR数字滤波器设计及MATLAB实现 目录 一、设计目的 ………………………………………………………………3 二、设计任务…………………………………………………………………3 三、设计原理…………………………………………………………………3 3.1窗函数法 ................................................3 3.2频率采样法 ................................................4 3.3最优化设计.................................................5 3.3.1等波纹切比雪夫逼近准则................................5 3.3.2仿真函数..............................................6 四、设计过程…………………………………………………………………7 五、收获与体会 ……………………………………………………………13 参考文献………………………………………………………………………13 FIR数字滤波器设计及MATLAB实现 一、设计目的 FIR滤波器：有限长单位冲激响应滤波器，是数字信号处理系统中最基本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。FIR滤波器设计的任务是选择有限长度的，使传输函数满足一定的幅度特性和线性相位要求。由于FIR滤波器很容易实现严格的线性相位，所以FIR 数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。 设计过程一般包括以下三个基本问题： （1） 根据实际要求确定数字滤波器性能指标； （2） 用一个因果稳定的系统函数去逼近这个理想性能指标； （3） 用一个有限精度的运算去实现这个传输函数。 三、设计原理 FIR滤波器设计的任务是选择有限长度的，使传输函数满足一定的幅度特性和线性相位要求。由于FIR滤波器很容易实现严格的线性相位，所以FIR数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。 设计过程一般包括以下三个基本问题： （1） 根据实际要求确定数字滤波器性能指标； （2） 用一个因果稳定的系统函数去逼近这个理想性能指标； （3） 用一个有限精度的运算去实现这个传输函数。 3.1窗函数法 设计FIR数字滤波器的最简单的方法是窗函数法，通常也称之为傅立叶级数法。FIR数字滤波器的设计首先给出要求的理想滤波器的频率响应,设计一个FIR数字滤波器频率响应，去逼近理想的滤波响应。然而，窗函数法设计FIR数字滤波器是在时域进行的，因而必须由理想的频率响应推导出对应的单位取样响应，再设计一个FIR数字滤波器的单位取样响应去逼近。设计过程如下： （3-1） 加窗的作用是通过把理想滤波器的无限长脉冲响应乘以窗函数来产生一个被截断的脉冲响应，即并且对频率响应进行平滑。MATLAB工具箱提供的窗函数有：矩形窗(Rectangularwindow)、三角窗(Triangular window)、布拉克曼窗(Blackman window)、汉宁窗(Hanningwindow)、海明窗(Hamming window)、凯塞窗(Kaiser window)、切比雪夫窗(Chebyshev window)。 窗函数主要用来减少序列因截断而产生的Gibbs效应。但当这个窗函数为矩形时，得到的FIR滤波器幅频响应会有明显的Gibbs效应，并且任意增加窗函数的长度（即FIR滤波器的抽头数）Gibbs效应也不能得到改善。为了克服这种现象,窗函数应该使设计的滤波器： （1） 频率特性的主瓣宽度应尽量窄,且尽可能将能量集中在主瓣内； （2） 窗函数频率特性的旁瓣ω趋于π 的过程中，其能量迅速减小为零。 3.2频率采样法 频率采样法是从频域出发,根据频域采样定理，对给定的理想滤波器的频率响应加以等间隔的抽样 ，得到： k=0,1,…,N-1 （3-2） 再利用可求得FIR滤波器的系统函数及频率响应。 而在各采样点间的频率响应则是其的加权内插函数延伸叠加的结果。但对于一个 无限长的序列，用频率采样法必然有一定的逼近误差,误差的大小取决于理想频响曲线的形状, 理想频响特性变换越平缓, 则内插函数值越接近理想值,误差越小。为了提高逼近的质量，可以通过在频率相应的过渡带内插入比较连续的采样