数字信号处理教程—matlab释义与实现第08章.ppt

数字信号处理教程—matlab释义与实现第08章.ppt

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数字信号处理教程—matlab释义与实现第08章

第八章 IIR滤波器设计 IIR滤波器设计 本章主要内容 8.1 关于模拟滤波器的预备知识 8.2 模拟滤波器原型的设计 8.3 由模拟滤波器变换为数字滤波器 8.4 模拟滤波器的频带变换 8.5 数字滤波器的频带变换 8.6 IIR数字滤波器的直接设计函数 8.7 IIR和FIR数字滤波器的比较 8.8 MATLAB中的滤波器设计工具 8.1 模拟滤波器的预备知识 模拟滤波器以是平方幅度响应作为指标。通常把通带最大幅特性归一化为一,故以相对指标的形式给出。对于低通滤波器,在通带和阻带内分别为: 关于模拟滤波器的预备知识 图中ε为通带波动系数,Ωp和Ωs是通带和阻带边缘频率。A为阻带衰减参数。参数ε和A化成dB值后可与第七章中的数字滤波器指标Rp和As建立关系。其关系式为: 关于模拟滤波器的预备知识 模拟滤波器技术指标以平方幅度响应的形式给出,所以它不包含相位信息。 为了从平方幅度响应计算s域系统函数,数学上就是要由复平面频率轴上的函数,推算出全复平面上的函数,那是没有唯一解的。从可实现性考虑,要作两点限制:一是规定函数的结构形式为s的实有理分式;二是此函数(也即其零极点)要关于虚轴对称,满足 关于模拟滤波器的预备知识 于是模拟滤波器设计问题可以表述如下。先找到一个满足幅特性要求的Ω2的有理真分式, 然后取出它的左半平面零点zj和极点pk,用这些零极点可以组成一个稳定的模拟滤波器。 关于模拟滤波器的预备知识 在这样的意义下,我们说模拟滤波器的传递函数Ha(s)可由平方幅度函数指标直接设计出来。 IIR数字滤波器的设计方法 IIR滤波器设计的基本技术就是利用复变函数映射把模拟滤波器转换成数字滤波器。这种技术的优点是:模拟滤波器设计方法已很成熟。而模拟到数字滤波器的复值映射也是成熟的。设计过程中的关键是进行两种不同的变换:(1)频带变换(低通到高通,带通,带阻等的变换);(2)模拟/数字变换;根据这两种变换的先后次序,引出两种设计方法: IIR数字滤波器的两种设计方法 8.2 模拟滤波器原型设计 巴特沃思低通滤波器 这种滤波器的特征是其通带和阻带都有平坦的幅度响应。N阶滤波器的平方幅度表达式为: 其中,N为滤波器的阶数,Ωc是边缘频率(单位:弧度/秒)。 模拟滤波器原型的设计 取N=1,2,10,得到的平方幅度响应见图8.2.1。 模拟滤波器原型的设计 为确定系统函数Ha(s),对J(Ω)进行零极点分解, 特征方程 的根就是系统的极点, 且此滤波器没有零点。 这2N个极点为 它们关于jΩ对称,并绝不会落在虚轴上。 模拟滤波器原型的设计 巴特沃思滤波器平方幅特性的2N个极点分布图。其左半平面的N个极点构成了它的系统函数。 模拟滤波器原型的设计 例 8.2.1 给定 试确定模拟 滤波器的系统函数。 解:根据给定的条件:N=3,Ωc=0.5。Ha(s)Ha(-s)的极点如上图左所示。取出左半平面的三个极点,可写出系统函数 模拟滤波器原型的设计 MATLAB提供了一个叫做buttap的函数,来设计N阶归一化(即Ωc=1)巴特沃思模拟原型滤波器,其调用格式为[z0,p0,k0]=buttap(N)。 去归一化的零极点及增益为z=z0Ωc、p=p0Ωc 及k=k0ΩcN 。 要求出系统函数分子分母的系数b和a,可用 b=poly(z); p=poly(p) 这使得高阶巴特沃斯滤波器的设计计算非常方便,不必用查表的老办法。 模拟滤波器原型的设计 计算本例的程序hc821应为: N=input(N= ); % 输入滤波器参数 OmegaC=input(OmegaC= ); [z0,p0,k0] = buttap(N) % 调用buttap函数, p = p0*OmegaC; % 非归一化极点 z = z0*OmegaC; % 非归一化零点 k = k0*OmegaC^N; % 非归一化增益 b = real(poly(z));b = k*b % 分子系数向量 a = real(poly(p)) % 分母系数向量 由此程序可求b,a,并得到系统函数的任何形式。 模拟滤波器原型的设计 设计过程中,需要由给定的滤波器指标参数 Ωp,Rp, Ωs和As,求得滤波器的阶数N和边缘频率Ωc,即要把滤波器公式变为它的设计方程。 令滤波器在通带和阻带边缘频率的幅特性满足指标Rp和As,建立两个联立方程。解出N和Ωc: 模拟原型的设计例8.2.2 例8.2.2 设通带边缘频率fp=5kHz,最大衰减Rp为1dB,阻带边缘频率fs=12kHz,最大衰减As

文档评论(0)

sandaolingcrh + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档