IIR滤波器背景、现状分析及设计实现浅析.doc

IIR滤波器背景、现状分析及设计实现浅析 摘要：目前IIR滤波器的设计方法主要有两种，数模转换法和计算机辅助设计法。数模转化法就是借助模拟滤波器进行设计。先按照数字滤波器的设计指标得到模拟滤波器的系统函数，再通过一定的转换方式将其转换为数字滤波器。而计算机辅助设计法就是将滤波器系数作为要求解的未知数，根据设计指标设联立方程，最后求解得出滤波器系数。 在IIR滤波器的设计当中，使用PSO算法能简化步骤。种群中每个粒子的位置都是每一组滤波器系数解，而适应函数则是该系数所代表的滤波器幅频响应和理想滤波器之间的均方差之和。这样适应函数越小，滤波器越接近设计目标。经过多次迭代，我们求得的全局最佳值即是最好的一组滤波器系数。 关键词：IIR滤波器，双线性转换法，脉冲响应不变法，巴特沃兹滤波器，粒子群优化算法 第一章 绪论 1.1 研究数字滤波器的背景 1.1.1 数字滤波器 现代社会中，人们的日常生活越来越数字化，越来越多的数字信号需要处理。处理数字信号，不免要用到数字滤波器。滤波器可以用来滤除信号中的噪声部分，提取出我们需要的信号频段，因此在图像和语音处理，医学生物信号处理以及其他领域都有着广泛应用。 我们把一个能够接收数字信号序列，然后将接收到的信号序列经过某种具有步骤有穷性和精度有限性的算法，输出符合使用者要求的选择性传输特性的数字信号的数字信号变换系统叫做数字信号滤波器[1]。按照滤波器的处理的信号类型，可以将其分为数字信号滤波器和模拟信号滤波器两大类。按照滤波器的作用范围，又可分为带通，高通，低通等。数字滤波器由于其灵活性强，可随时调整滤波器频谱,滤波稳定性高，精度高，并且数字信号处理场景要远远多于模拟信号，所以应用范围一般比模拟滤波器要广。按照数字滤波器的脉冲响应时域特性，我们把其分为两种，无限脉冲响应(IIR，Infinite Impulse Response)滤波器和有限脉冲响应(FIR，Finite Impulse Response)滤波器[3]。 1.1.2 IIR滤波器 IIR滤波器是一种递归型的滤波器。它的单位脉冲响应为无限长，且具有反馈回路。而FIR滤波器正好与其相反，它的单位脉冲响应为有限长，且没有反馈回路。IIR滤波器的幅频特性较好，但输出信号相位非严格线性，即具有不同频率分量的原始信号在经过IIR系统之后，他们的时间差会发生变化。这就导致IIR滤波器通常只适用于语音通话，以及其他对信号相位要求不严格的场合。在对信号要求较高时，需要添加相位校准网络。 1.2 国内外研究现状 目前IIR滤波器的设计方法主要有两种，数模转换法和计算机辅助设计法。数模转化法就是借助模拟滤波器进行设计。先按照数字滤波器的设计指标得到模拟滤波器的系统函数，再通过一定的转换方式将其转换为数字滤波器。而计算机辅助设计法就是将滤波器系数作为要求解的未知数，根据设计指标设联立方程，最后求解得出滤波器系数。 国内专家学者研究的重点主要针对特定用途数字滤波器，从设计到硬件实现，主要还是借助模拟滤波器的设计方法进行设计。模拟滤波器设计理论非常成熟，而且有许多现成的滤波器可供选择，例如巴特沃兹滤波器，切比雪夫滤波器，椭圆滤波器等等。这些滤波器设计方案相对成熟，且幅频响应特性较好，所以数模转换法一般应用比较广泛。 而国外的专家则侧重于采用计算机辅助的方法，通过优化算法得到数字滤波器系统函数的参数。例如遗传算法GA，模拟退火算法SA，粒子群算法PSO等等。这些算法一般是根据最小均方误差准则或绝对误差准则进行计算，然后改变参数赋值再计算误差，如此迭代直至误差最小，从而得到滤波器的系数。 1.3 本文的研究内容 本文主要内容是基于MATLAB平台完成IIR滤波器的设计和实现。 IIR滤波器的设计主要有两类方法。第一类采用的是间接设计法，即通过数模转换的方式进行设计。主要有两种转换方式，双线性转换法和脉冲响应不变法。在本文第二章中，着重介绍了这两种方法的具体原理，给出了转换关系，并使用MATLAB编写程序实现了相关算法，设计出了符合参数要求的IIR滤波器，并给出其幅频响应特性曲线。 第二类方法是直接设计法，即采用计算机辅助的方式确定IIR滤波器的参数系数。本文主要采用粒子群算法PSO，第三章主要介绍了PSO算法的基本原理，并使用MATLAB编写程序实现PSO算法，设计出符合参数要求的IIR滤波器，并绘制出其幅频响应特性曲线。 第二章 采用数模转换法设计IIR滤波器 2.1 IIR滤波器的简介 2.1.1 IIR滤波器的特性 IIR滤波器的系统函数如下： 其中N为该IIR系统的阶数