资源描述：
《IIR数字滤波器在的MATLAB实现》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、第6章IIR数字滤波器的MATLAB实现6.2实验原理及实例分析6.2.1IIR数字滤波器的传递函数及特点设IIR滤波器的输入序列为x(n)，则IIR滤波器的输入序列x(n)与输出序列y(n)之间的关系可以用下面的方程式表示：其中，和是滤波器的系数，其中中至少有一个非零。与之相对应的差分方程为：由传递函数可以发现无限常单位冲激响应滤波器有如下特点：（1）单位冲激响应h(n)是无限长的。（2）系统传递函数H(z)在有限z平面上有极点存在。（3）结构上存在着输出到输入的反馈，也就是结构上是递归型的。6.2.2IIR数字滤波器的设计与实
2、现IIR数字滤波器的设计有多种方法，如频率变换法、数字域直接设计以及计算辅助设计等。下面只介绍频率变换设计法。首先考虑由模拟低通滤波器到数字低通滤波器的转换，其基本的设计过程如下:（1）将数字滤波器的技术指标转换为模拟滤波器的技术指标；（2）设计模拟滤波器G(S)；（3）将G(S)转换成数字滤波器H(Z)；在低通滤波器的设计基础上，可以得到数字高通、带通、带阻滤波器的设计流程如下：（1）给定数字滤波器的设计要求（高通、带阻、带通）；（2）转换为模拟（高通、带阻、带通）滤波器的技术指标；（3）转换为模拟低通滤波器的指标；（4）设计得
3、到满足第三步要求的低通滤波器传递函数；（5）通过频率转换得到模拟（高通、带阻、带通）滤波器；（6）变换为数字（高通、带阻、带通）滤波器。在matlab中设计IIR滤波器的方法及其它们所用到的函数如表6－1所示。表6－1matlab中设计IIR滤波器的方法列表方法描述函数模拟原型法采用经典低通滤波器作为连续域上的设计模型，通过频率变换得到IIR数字滤波器，最后进行离散化处理完整设计函数：Beself,butter,cheby1,cheby2,ellip滤波器的阶估计函数：Buttord,cheb1ord,cheb2ord,ellip
4、ord低通模拟滤波器原型函数：beselap,buttap,cheb1ap,cheb2ap,ellipap频域变换函数：Lp2bp,lp2bs,lp2hp,lp2lp其他函数：Bilinear,impinvar直接设计方法直接在离散时域上估计线性的幅度响应yulewalk通用butterworth设计方法使用butterworth设计低通数字滤波器Maxflat参数建模方法寻找接近于所需要设计的滤波器的通用原型时域上的建模函数：Lpc,prony,stmcb频域上的建模函数：Invfreqs,invfreqz表6－2频率转换函数列
5、表频率转换转换函数低通到低通[numt,dent]=lp2lp(num.den,w0)[At,Bt,Ct,Dt]=lp2lp(A,B,C,D,w0)低通到高通[numt,dent]=lp2hp(num.den,w0)[At,Bt,Ct,Dt]=lp2hp(A,B,C,D,w0)低通到带通[numt,dent]=lp2bp(num.den,w0)[At,Bt,Ct,Dt]=lp2bp(A,B,C,D,w0)低通到带阻[numt,dent]=lp2bs(num.den,w0)[At,Bt,Ct,Dt]=lp2bs(A,B,C,D,w0
6、)6.2.3标准数字滤波器设计函数Matlab提供了一组标准的数字滤波器设计函数，大大简化了滤波器的设计过程。1、butter功能：Butterworth模拟/数字滤波器设计格式：[b,a]=butter(n,wn,'ftype',’s’)[b,a]=butter(n,wn,'ftype')说明：l选项中加入‘S’用于设计各种模拟Butterworth滤波器；不加设计各种数字Butterworth滤波器lFtype为缺省，设计低通滤波器lFtype＝hign，设计高通滤波器lFtype＝stop，设计带阻滤波器【实例6-1】设计一
7、个5阶Butterworth数字高通滤波器，阻带截止频率为250Hz。设采样频率为1000Hz。解：源代码如下：[b,a]=butter(5,250/500,'high')[z,p,k]=butter(5,250/500,'high')freqz(b,a,512,1000)程序运行后，产生结果如下所示。b=0.0528-0.26390.5279-0.52790.2639-0.0528a=1.0000-0.00000.6334-0.00000.0557-0.0000z=11111p=0.0000+0.7265i0.0000-0.72
8、65i0.0000+0.3249i0.0000-0.3249i0.0000k=0.0528图6－15阶Butterworth数字高通滤波器1、Cheby1、Cheby2功能：chebyshevI、chebyshevII型模拟/数字滤波器设计格式：[
