基于matlab的IIR数字滤波器设计.docx

基于matlab的IIR数字滤波器设计IIR数字滤波器介绍1.IIR数字滤波器的基本原理所谓数字滤波器是指输入输出均为数字信号通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的硬件。实质上就是一个由有限精度算法实现的线性时不变离散系统。它的基本工作原理是利用离散系统的特性对系统输入信号进行加工和变换改变输入序列的频谱或信号波形让有用的频率分量通过抑制无用的信号分量输出因此数字滤波与模拟滤波的概念相同根据其频率特性同样可以分为低通高通带通带阻只是信号的形式和实现滤波方式有所不同。如果要处理的信号是模拟信号就可以通过A/D或者D/A转换在信号形式上进行匹配转换同样可以使用数字滤波器对模拟信号进行滤波。数字滤波器滤波的数学表达式：y（n）=x(n)*h(n);如果滤波器的输入输出信号都是离散信号那么该滤波器的脉冲响应也一定是离散信号这样的滤波器就成为了数字滤波器。上面的系统为时域离散系统时其频域特性为：其中分别是数字滤波器的输出序列和输入序列的频域响应是数字滤波器的频域响应。可以看见按照输入信号的频谱特点和处理信号的目的适当选择滤波器的频域响应使得滤波后的输出信号满足设计性能要求就是滤波器的滤波原理。2．IIR数字滤波器传输特性IIR数字滤波器的系统函数可以表示为：H(Z)=式中H(Z)称为N阶IIR滤波器函数。3．.数字滤波器的技术要求.我们通常设计的数字滤波器一般属于选频滤波器。我们的目的是要设计一个因果可实现的滤波器另外买也要考虑到成本和复杂性问题因此实用中通带和阻带都允许一定的误差容限即通带不一定是完全水平的阻带也不可能完全衰减到零。而且通带和阻带之间还要设置一定带宽的过渡带。如下图表示低通滤波器的技术要求：图中分别表示通带截止频率和阻带截止频率通带频率范围为0≤w≤通带中要求（1-δ1）≤|H≤1阻带截止频率范围≤w≤Π再阻带中要求≤δ2从到称为过渡带在这个频带内幅度响应从通带平滑的下落到阻带。IIR数字滤波器的设计方法1．常用的IIR滤波器设计方法（1）.以模拟滤波器函数为基础的变换法；即先设计一满足指定条件的模拟滤波器H(s)再将该模拟滤波器转化为数字滤波器H(z)。（2）.直接设计法：在z平面内根据零、极点对系统特性的影响调整零极点位置得H(z)。（3）.最优化设计法(计算机辅助设计)在某种最小化误差准则下建立差分方程系数ak、bi对理想特性的逼近方程使用迭代方法解方程组得到最佳逼近系统。由于此方法计算量大需要借助于计算机进行设计。下面分别介绍函数设计法和信号处理图形界面来设计IIR数字滤波器得出最优化的设计方法。2.基于matlab的函数设计IIR数字滤波器IIR数字滤波器设计的一般方法是先设计低通模拟滤波器进行频率变换将其转换为相应的（高通带通等）模拟滤波器在转换为高通带通或带阻数字滤波器由模拟滤波器设计数字滤波器的方法。这是因为模拟滤波器设计方法已经很成熟它不仅有完整的设计公式还有完善的图表供查阅另外还有一些典型的滤波器类型可供我们使用。对设计的全过程的各个步骤MATLAB都提供了了相应的工具箱函数使IIR数字滤波器设计变得非常简单。2.1.基于matlab函数的IIR滤波器设计2.1.1.设计方法选择：程序设计法是基于MATLAB中相应的工具箱函数来实现的IIR数字滤波器的设计步骤可由图1.2所示的流程图来表示。这个图也清晰的表示了5类20个信号处理工具箱函数的作用在MATLAB中模拟滤波器的系统函数H(S)=数字滤波器的系统函数H(Z)=在实际工程中需要的设计结果是系数向量B和A用B和A来综合滤波器的硬件实现结构或软件运算结构为了直观的看出设计结果本文的实例均以滤波器幅频响应曲线作为设计结果输出。如果需要滤波器系数在运行程序后只要在MATLAB命令窗口键入系数向量名则相应的系数就显示出来了。2.1.2.程序设计实例分析（a）设计高通和带通Butterworth数字滤波器我们给出四阶归一化Butterworth模拟滤波器的系统函数用双线性变换法从Ha（s）设计四阶带通butterworth数字滤波器并图示||设计采样周期T=1s指标如下现在我们分步进行：■建模由于本例主要涉及三个问题：由数字滤波器指标求相应的模拟滤波器指标；模拟滤波器频率变换（因为已给定阶数和模拟滤波器的归一化低通原型）；由相应的模拟滤波器到数字滤波器（双线性变换法）。由于调用bilinear函数将模拟滤波器转换成数字滤波器非常容易并且有效抑制频率失真的问题本例给定了数字滤波器指标所以首先要设计处与该指标相应的四阶Butterworth模拟滤波器然后调用bilinear函数将其转