- 1、本文档共26页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
北京交通大学(数字信号处理研究生课程)chIIR数字滤波器的基本结构概要
IIR数字滤波器的基本结构 滤波器结构及有限字长 近代数字信号处理(Advanced Digital Signal Processing) 电子信息工程学院 信号与图像处理研究室 IIR数字滤波器设计的基本思想 模拟低通滤波器设计 模拟域频率变换 脉冲响应不变法 双线性变换法 IIR数字滤波器的基本结构 利用MATLAB设计IIR DF IIR数字滤波器的基本结构 直接型结构 级联型结构 并联型结构 IIR数字滤波器的直接型结构 直接 I 型结构 设M=N IIR数字滤波器的直接型结构 直接 I 型结构 IIR数字滤波器的直接型结构 将H2(z)系统的延时器与H1(z)系统的延时器共用 交换H1(z) H2(z)两级联子系统的级联顺序 直接 II 型结构 IIR数字滤波器的直接型结构 转置直接 II 型结构 IIR数字滤波器的直接型结构 IIR数字滤波器的直接型结构优缺点 优点:简单直观 缺点: 1. 改变某一个{ak }将影响所有的极点。 2. 改变某一个{bk }将影响所有的零点。 3. 对有限字长效应太敏感,容易出现不稳定现象。 对于三阶以上的IIR滤波器,几乎都不采用直接型结构,而是采用级联型、并联型等其他形式的结构。 IIR数字滤波器的直接型结构 将滤波器系统函数H(z)的分子和分母分解为一阶和二阶实系数因子之积的形式 画出各二阶基本节的直接型结构,再将它们级联。 二阶基本节 IIR数字滤波器的级联型结构 IIR数字滤波器的级联型结构 基于直接II型的级联型结构 级联型结构信号流图 IIR数字滤波器的级联型结构 基于转置直接II型的级联型结构 级联型结构信号流图 IIR数字滤波器的级联型结构优点 硬件实现时,可以用一个二阶节进行时分复用。 每一个基本节系数变化只影响该子系统的零极点。 对系数变化的敏感度小,受有限字长的影响比直接型低。 IIR数字滤波器的级联型结构 将滤波器系统函数H(z)展开成部分分式之和,并将一阶系统仍采用二阶基本节表示 画出各二阶基本节的直接型结构,再将它们并联。 IIR数字滤波器的并联型结构 IIR数字滤波器的并联型结构 并联型结构信号流图 基于直接II型的 并联型结构 基于转置直接II型 的并联型结构 IIR数字滤波器的并联型结构 并联型结构信号流图 IIR数字滤波器的并联型结构优缺点 优点: 1. 运算速度快。 2. 各基本节的误差互不影响。 3. 可以单独调整极点的位置。 缺点: 不能向级联型那样直接调整零点。 IIR数字滤波器的并联型结构 例:已知某三阶数字滤波器的系统函数为 试画出其直接型、 级联型和并联型结构。 将系统函数H(z)表达为 解: 直接型 例:已知某三阶数字滤波器的系统函数为 试画出其直接型、 级联型和并联型结构。 解: 级联型 将系统函数H(z)表达为一阶、二阶实系数分式之积 例:已知某三阶数字滤波器的系统函数为 试画出其直接型、 级联型和并联型结构。 解: 并联型 将系统函数H(z)表达为 部分分式之和的形式 例题结论:由系统函数画直接型结构的规律 (1) 画出n个级联的延时器。 (2) 将各延时器的输出反馈连接到输入端的加法器形成反馈回路,这些反馈回路的系统函数分别为 (3) 将输入端加法器的输出和各延时器的输出正向连接到输出端的加法器构成前向通路,各条前向通路的系统函数分别为 例题结论:由系统函数画级联型和并联型结构 级联型结构: (1) 将系统函数H(z)表达为一阶、二阶实系数分式之积 (2) 分别画出每个一阶、二阶子系统的直接型结构,再将它们级联。 并联型结构: (1) 将系统函数H(z)表达为一阶、二阶实系数分式之和 (2) 分别画出每个一阶、二阶子系统的直接型结构,再将它们并联。 H(z)=1+a z-N H0(z)=1/H(z)=1/(1+a z-N) 信号处理(回声消除) 信号处理(回声消除) 信号处理(回声消除) load mtlb; x=mtlb; M=length(x); a=0.9; N=3000; x0=[x,zeros(1,N)]; xN=[zeros(1,N),x]; y=x0+a.*xN; wavwrite(y,8000,32,H:\2013dsp\Matlab_programs\echo11.wav) %sound(y); yy=xcorr(y,y); L=length(yy); %wavwrite(y,a.wav); subplot(321);stem(x0,.); title(Signal x[k]));grid on subplot(323);
文档评论(0)