IIR数字滤波器设计及实现.docVIP

下载本文档

24
0
约5.4千字
约 10页
2020-06-28 发布于湖北
举报
版权申诉

IIR数字滤波器设计及实现.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

实验三 IIR数字滤波器设计及实现一、实验目的（1）熟悉用脉冲响应不变法和双线性变换法设计IIR数字滤波器的原理与方法；（2）学会调用MATLAB信号处理工具箱中滤波器设计函数设计IIR数字滤波器，学会根据滤波需求确定滤波器指标参数。二、实验原理设计IIR数字滤波器一般采用脉冲响应不变法和双线性变换法。脉冲响应不变法：根据设计指标求出滤波器确定最小阶数N和截止频率Wc；计算相应的模拟滤波器系统函数；将模拟滤波器系统函数转换成数字滤波器系统函数双线性变换法：根据数字低通技术指标得到滤波器的阶数N；取合适的T值，几遍校正计算相应模低通的技术指标；根据阶数N查表的到归一化低通原型系统函数，将代入，去归一化得到实际的；用双线性变换法将转换成数字滤波器三、实验内容及步骤 1、用脉冲响应不变法设计（1）根据设计指标求出滤波器确定最小阶数N和截止频率Wc clear;close all;clc; % 开始准备 fp=3400;fs=5000;Fs=22050;Rp=2;Rs=20;T=1/Fs; % T=1s的模拟滤波器设计指标 W1p=fp/Fs*2; W1s=fs/Fs*2; % 求归一化频率 [N, Wn] = buttord(W1p, W1s, Rp, Rs, s; % 确定butterworth的最小阶数N和频率参数Wn 得到结果为： N = 7 Wn = 0.3266 即：该设计指标下的模拟滤波器最小阶数为N=7，其截至频率为Wn =0.3266；（2）计算相应的模拟滤波器系统函数 clear;close all;clc; % 开始准备 fp=3400;fs=5000;Fs=22050;Rp=2;Rs=20;T=1/Fs; % T=1s的模拟滤波器设计指标 W1p=fp/Fs*2; W1s=fs/Fs*2; % 求归一化频率 [N, Wn] = buttord(W1p, W1s, Rp, Rs, s; % 确定butterworth的最小阶数N和频率参数Wn [B,A]=butter(N,1,s %计算相应的模拟滤波器系统函数得到结果为： B = 1.0e-003 * 0 0 0 0 0 0 0 0.3966 A = 1.0000 1.4678 1.0773 0.5084 0.1661 0.0375 0.0055 0.0004 （3）将模拟滤波器系统函数转换成数字滤波器系统函数 clear;close all;clc; % 开始准备 fp=3400;fs=5000;Fs=22050;Rp=2;Rs=20;T=1/Fs; % T=1s的模拟滤波器设计指标 W1p=fp/Fs*2; W1s=fs/Fs*2; % 求归一化频率 [N, Wn] = buttord(W1p, W1s, Rp, Rs, s; % 确定butterworth的最小阶数N和频率参数Wn [B,A]=butter(N,1,s； %计算相应的模拟滤波器系统函数 [Bz,Az]=impinvar(B,A %用脉冲相应不变法将模拟滤波器转换成数字滤波器 sys=tf(Bz,Az,T； %得到传输函数 Bz = 1.0e-004 * -0.0000 0.0045 0.2045 0.8747 0.7094 0.1090 0.0016 0 Az = 1.0000 -5.5415 13.2850 -17.8428 14.4878 -7.1069 1.9491 -0.2304 即：由Bz和Az可以写出数字滤波器系统函数为： Transfer function: -9.992e-015 z^7 + 4.454e-007 z^6 + 2.045e-005 z^5 + 8.747e-005 z^4 + 7.094e-005 z^3 + 1.09e-005 z^2 + 1.561e-007 z --------------------------------------------------------------------------------------------------------- z^7 - 5.541 z^6 + 13.28 z^5 - 17.84 z^4 + 14.49 z^3 - 7.107 z^2 + 1.949 z - 0.2304 Sampling time: 4.5351e-005 （4）绘图 clear;close all;clc; % 开始准备 fp=3400;fs=5000;Fs=22050;Rp=2;Rs=20;T=1/Fs; % T=1s的模拟滤波器设计指标 W1p=fp/Fs*2; W1s=fs/Fs*2; % 求归一化频率 [N, Wn] = buttord(W1p, W1s, Rp, Rs,