等波纹低通滤波器的设计及与其他滤波器的比较

燕山大学课程设计说明书题目：等波纹低通滤波器的设计学院（系）：里仁学院年级专业：仪表10-2学号：学生姓名：指导教师：教师职称：燕山大学课程设计说明书-1-燕山大学课程设计（论文）任务书院（系）：电气工程学院基层教学单位：自动化仪表系学号学生姓名专业（班级）仪表10-2设计题目等波纹低通滤波器的设计设计技术参数等波纹（equiripple）滤波器指标：通带截止频率为4KHz,阻带截止频率为5KHz,阻带至少衰减40dB采样截止频率20kHz。设计要求仿真产生一个连续信号，包含低频，高频分量，对其进行采样，进行频谱分析。设计等波纹低通滤波器对信号进行滤波处理，观察滤波后的频谱。分析该类型滤波器与其他类型低通滤波器（如Butterworth）优势及特点。参考资料布置课程设计任务,查阅资料,确定方案；进行编程设计；完成课程设计报告书。周次前半周后半周应完成内容收集消化资料、学习MATLAB软件，进行相关参数计算。编写仿真程序、调试。指导教师签字基层教学单位主任签字说明：此表一式四份，学生、指导教师、基层教学单位、系部各一份。2013年7月5日燕山大学课程设计说明书-2-摘要等波纹最佳逼近法是一种优化设计法，它克服了窗函数设计法和频率采样法的缺点，使最大误差（即波纹的峰值）最小化，并在整个逼近频段上均匀分布。用等波纹最佳逼近法设计的FIR数字滤波器的幅频响应在通带和阻带都是等波纹的，而且可以分别控制通带和阻带波纹幅度。这就是等波纹的含义。最佳逼近是指在滤波器长度给定的条件下，使加权误差波纹幅度最小化。与窗函数设计法和频率采样法比较，由于这种设计法使滤波器的最大逼近误差均匀分布，所以设计的滤波器性能价格比最高。阶数相同时，这种设计法使滤波器的最大逼近误差最小，即通带最大衰减最小，阻带最小衰减最大；指标相同时，这种设计法使滤波器阶数最低。实现FIR数字滤波器的等波纹最佳逼近法的MATLAB信号处理工具函数为remez和remezord。Remez函数采用数值分析中的remez多重交换迭代算法求解等波纹最佳逼近问题，求的满足等波纹最佳逼近准则的FIR数字滤波器的单位脉冲响应h（n）。由于切比雪夫和雷米兹对解决该问题做出了贡献，所以又称之为切比雪夫逼近法和雷米兹逼近法。关键词：FIR数字滤波器MATLABremez函数remezord函数等波纹燕山大学课程设计说明书-3-目录摘要--------------------------------------------------------------------------------------------2关键字------------------------------------------------------------------------------------------2第一章第一章数字滤波器的基本概-------------------------------------------------41.1滤波的涵义----------------------------------------------------------------------41.2数字滤波器的概述-------------------------------------------------------------41.3数字滤波器的实现方法-------------------------------------------------------41.4.数字滤波器的可实现性------------------------------------------------------51.5数字滤波器的分类-------------------------------------------------------------51.6FIR滤波器简介及其优点----------------------------------------------------5-第二章等波纹最佳逼近法的原理-------------------------------------------------------52.1等波纹最佳逼近法概述-------------------------------------------------------92.2.等波纹最佳逼近法基本思想-------------------------------------------------92.3等波纹滤波器的技术指标及其描述参数介绍---------------------------102.3.1滤波器的描述参数-----------------------------------------------------102.3.2设计要求-----------------------------------------------------------------10第三章matlab程序------------------------------------------------------------------------11第四章该型滤波器较其他低通滤波器的优势及特点--------------------12第五章课程设计总结---------------------------------------------------------------------15参考文献资料-------------------------------------------------------------------------------15燕山大学课程设计说明书-4-第一章数字滤波器的基本概念介绍1.1滤波的涵义a)将输入信号的某些频率成分或某个频带进行压缩、放大;b)对信号进行检测;c)对参数估计;1.2数字滤波器的概述所谓数字滤波器，是指输入、输出均为数字信号，通过数值运算处理改变输入信号所含频率成分的相对比例，或者滤除某些频率成分的数字器件或程序。下图给出了一个具有模拟输入信号和输出信号的实时数字滤波器的简化方框图。这个带限模拟信号被周期地抽样，且转化成一系列数字X（n）（n=0,1,…）。数字处理器依据滤波器的计算算法，执行滤波运算，把输入系列X（n）映射到输出系列Y（n）。DAC把数字滤波后的输出转化成模拟值，这些模拟值接着被模拟滤波器平滑，并且消去不想要的高频分量。输入滤波器带抽样的ADC数字处理器DAC输出滤波器X(t)模拟输入X(n)Y(n)Y(t)模拟输出一个具有模拟输入和输出信号的实时数字滤波器的简化方框图数字滤波器在数字信号处理中具有非常重要的地位。在许多应用中（例如数据压缩，生物医学信号处理、语音处理、图象处理、数据传输、数字音频、电话回声对消，等等），数字滤波器和模拟滤波器比数字滤波器的有优势更加明显。1.3数字滤波器的实现方法a)用软件在计算机上实现燕山大学课程设计说明书-5-b)用专用的数字信号处理芯片c)用硬件1.4.数字滤波器的可实现性a)要求系统因果稳定，即所设计的系统极点全部集中在单位圆内。b)要求系统的差分方程的系数或者系统函数的系数为实数，即系统的零极点必须共轭成对出现，或者是实数。1.5数字滤波器的分类数字滤波器可以分为经典滤波器和现代滤波器。经典滤波器按照滤波特性可以分为数字高通、数字低通、数字带通、数字带阻等类型。数字滤波器从实现的网络结构或者从单位脉冲响应长度分类可以分为无限长单位脉冲（IIR）和有限长单位脉冲（FIR）。1.6FIR滤波器简介及其优点1.6.1FIR滤波器即有限长单位脉冲响应滤波器，是数字信号处理中最基本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位冲击响应是有限的，没有输入到输出的反馈，是稳定的系统。因此，FIR滤波器在通信、图像、语音、雷达等许多领域都有着十分广泛的应用。1.6.2FIR滤波器具有以下主要优点：a)FIR滤波器具有准确的线性相位b)FIR滤波器永远稳定c)FIR滤波器设计方法一般是线性的d)FIR滤波器在硬件上具有更高的运行效率e)FIR滤波器启动传输时间只需要有限时间1.6.3FIR滤波器的基本结构FIR滤波器的单位抽样响应为有限长度，一般采用非递归形式实现。通常的FIR数字滤波器有横截性和级联型两种。(a)FIR滤波器的横截型结构表示系统输入输出关系的差分方程可写作：燕山大学课程设计说明书-6-10()()()Nmynhmxnm(2-3)直接由差分方程得出的实现结构如图1所示：图1、横截型(直接型﹑卷积型)若h(n)呈现对称特性，即此FIR滤波器具有线性相位，则可以简化加横截型结构，下面分情况讨论：图2、N为奇数时线形相位FIR滤波器实现结构图3、N为偶数时线性相位FIR滤波器实现结构(b)FIR滤波器的级联型结构将H(z)分解成实系数二阶因子的乘积形式：[]121201201()()NNNkkkNkHzhnzbbzbz(2-4)这时FIR滤波器可用二阶节的级联结构来实现，每个二阶节用横截型结构实现。如下图所示燕山大学课程设计说明书-7-1.6.4线性相位FIR数字滤波器：对于长度为N的h(n)，频率响应函数为10)(NnjjeenhHH(ejω)=Hg（ω）ejθ(ω)式中，Hg（ω）称为幅度特性；θ（ω）称为相位特性。这里Hg（ω）为ω的实函数，可能取负值。线性相位FIR滤波器指θ（ω）是ω的线性函数。有两种情况：即θ（ω）=-τω和θ（ω）=θ0-τω。第一种是第一类线性相位，第二种为第二类线性，通常取θ0=-π/2。1.6.4.1线性相位FIR滤波器的时域约束条件线性相位FIR滤波器的时域约束条件是指满足线性相位时，对h(n)的约束条件。为使滤波器对实信号信号的处理结果仍为实信号，要求h(n)是实系列。线性相位FIR滤波器的时域约束条件可分为以下几种情况：第一类线性相位：h(n)关于n=(N-1)/2点偶对称对，即h(n)=h(N-n-1)。N确定时，其线性特性是一明确的线性函数：θ（ω）=-ω（N-1）/2。第二类线性相位：h(n)关于n=(N-1)/2点奇对称，即h(n)=-h(N-n-1)。N确定时，一般取其线性特性是一明确的线性函数：θ（ω）=-π/2-ω（N-1）/2。1.6.4.2线性相位FIR滤波器幅度特性Hg（ω）的特点:实质上，，幅度特性的特点就是线性相位FIR滤波器的频域约束条件。结合时域约束条件)1()(nNhnh，引入两个参数τ=（N-1）/2，2/)1(NM。分四种情况讨论幅度特性特点。情况1：)1()(nNhnh；N为奇数。此时：10)](cos[)(2)()(MnnnhhHg所以，Hg（ω）关于ω=0，π，2π对称。因此情况1可以实现低通、高通、带通、带阻等各种滤波器。燕山大学课程设计说明书-8-情况2：)1()(nNhnh；N为偶数。此时：MnnnhHg0)](cos[)(2)(所以Hg（π）=0，Hg（ω）关于ω=π奇对称，关于ω=0和2π偶对称。因此，情况2不能实现高通和带阻滤波器。情况3：)1()(nNhnh；N为奇数。此时：10)](sin[)(2)(MnnnhHg所以Hg（ω）关于ω=0，π,2π三点奇对称。因此，情况3只能实现带通滤波器。情况4：)1()(nNhnh；N为偶数。此时：MnnnhHg0)](sin[)(2)(所以Hg（ω）关于ω=0,2π两点奇对称。因此，情况4不能实现低通和带阻滤波器。线性相位FIR数字滤波器的零极点分布特点：将)1()(nNhnh代入-10)()(NnznhzHn得到)()(1)1(zzHzHN由上式可以得到，如果zzi是零点，则zi1，*zi，*)(1zi也必是其零点。FIR滤波器在z=0