1、 IIR 数字低通滤波器的 MATLAB 实现 一、设及计任务要求一、设及计任务要求： 设计任务：设计任务： 在 matlab 基础上，利用相关的命令实现语音信号的处理。 要要 求求： 1、查阅资料，写出命令。 2、 实现傅立叶变换 3、 设计数字低通滤波器。 4、 利用 Matlab GUI 工具做一个图形用户界面。 指导教师签名： 2012 年 6 月 日 二、指导教师评语二、指导教师评语： 指导教师签名： 2011 年 12 月 日 三、成绩三、成绩 验收盖章: 2011 年 12 月 日 目录目录错误错误!未找到引用源。未找到引用源。 一、引言 4 （一）设计目的 4 （二）设计思想
2、错误错误!未定义书签。未定义书签。 二、IIR 数字滤波器的 MATLAB 实现 三、实验步骤 错误错误!未定义书签。未定义书签。 四、实验结果及截图 错误错误!未定义书签。未定义书签。 结果 错误错误!未定义书签。未定义书签。 参考文献 错误错误!未定义书签。未定义书签。 附录 错误错误!未定义书签。未定义书签。 错误!未找到引用源。一、引言一、引言 设计目的设计目的 （1）熟悉 Matlab 在信号处理方面的应用 （2）深刻了解频域上的频谱特性。 （3）了解数字滤波器的设计。 （4）学会比较滤波前后语音信号的差异。 （5）利用 GUI 工具做一个简洁大方的图形用户界面。 二、二、IIRII
3、R 数字低通滤波器的数字低通滤波器的 MATLABMATLAB 实现实现 （一）（一）IIR 数字滤波器的传递函数及特点 设 IIR 滤波器的输入序列为 x(n)，则 IIR 滤波器的输入序列 x(n)与输出序列 y(n)之间的关系可以用 下面的方程式表示： )()()( 10 jnxainxbny N j j M i i 其中， j a和 i b是滤波器的系数，其中 j a中至少有一个非零。与之相对应的差分方程为： N N M M zaza zbzbb zX zY ZH 1 )( )( )( 1 1 1 10 由传递函数可以发现无限常单位冲激响应滤波器有如下特点： （1） 单位冲激响应 h(
4、n)是无限长的。 （2） 系统传递函数 H(z)在有限 z 平面上有极点存在。 （3） 结构上存在着输出到输入的反馈，也就是结构上是递归型的。 （二）（二）IIR 数字滤波器的设计与实现 IIR 数字滤波器的设计有多种方法，如频率变换法、数字域直接设计以及计算辅助设计等。下面只 介绍频率变换设计法。首先考虑由模拟低通滤波器到数字低通滤波器的转换，其基本的设计过程如下: （1） 将数字滤波器的技术指标转换为模拟滤波器的技术指标； （2） 设计模拟滤波器 G(S)； （3） 将 G(S)转换成数字滤波器 H(Z)； 在低通滤波器的设计基础上，可以得到数字高通、带通、带阻滤波器的设计流程如下： （1
5、） 给定数字滤波器的设计要求（高通、带阻、带通） ； （2） 转换为模拟（高通、带阻、带通）滤波器的技术指标； （3） 转换为模拟低通滤波器的指标； （4） 设计得到满足第三步要求的低通滤波器传递函数； （5） 通过频率转换得到模拟（高通、带阻、带通）滤波器； （6） 变换为数字（高通、带阻、带通）滤波器。 在 matlab 中设计 IIR 滤波器的方法及其它们所用到的函数如表 61 所示。 表 61 matlab 中设计 IIR 滤波器的方法列表 方法方法 描述描述 函数函数 模拟原型法模拟原型法 采用经典低通滤波器作为连续域上的设计模 型，通过频率变换得到 IIR 数字滤波器，最后 进行离
6、散化处理 完整设计函数： Beself,butter,cheby1,cheby2,ellip 滤波器的阶估计函数： Buttord,cheb1ord,cheb2ord,ellipord 低通模拟滤波器原型函数： beselap,buttap,cheb1ap,cheb2ap,ellipap 频域变换函数： Lp2bp,lp2bs,lp2hp,lp2lp 其他函数： Bilinear,impinvar 直接设计方直接设计方法法 直接在离散时域上估计线性的幅度响应 yulewalk 通用通用 butterworth 设计方法设计方法 使用 butterworth 设计低通数字滤波器 Maxflat 参数建模方法参数建模方法 寻找接近于所需要设计的滤波器的通用原型 时域上的建模函数： Lpc,prony,stmcb 频域上的建模函数： Invfreqs,invfreqz （三） 双线性变换法 与冲激响应不变法比较， 双线性变换的主要优点是靠频率的非线性关系得到 S 平面与 Z
