原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
实验6无限冲激响应数字滤波器设计
实验目的
掌握双线性变换法及脉冲响应不变法设计IIR数字滤波器的具体设计方法
熟悉用双线性变换法及脉冲响应不变法，设计低通、高通和带通IIR数字滤波器的
计算机编程
实验原理
在MATLAB中，可以用下列函数辅助设计IIR数字滤波器：
1)利用buttord和cheb1ord可以确定低通原型巴特沃斯和切比雪夫滤波器的阶数和截止频率；
2)[num,den]=butter（N,Wn）(巴特沃斯)和[num,den]=cheby1（N,Wn）,[num,den]=cheby2（N,Wn）(切比雪夫1型和2型)可以进行滤波器的设计；
3）lp2hp，lp2bp，lp2bs可以完成低通滤波器到高通、带通、带阻滤波器的转换；
4)使用bilinear可以对模拟滤波器进行双线性变换，求得数字滤波器的传输函数系数；
5)利用impinvar可以完成脉冲响应不变法的模拟滤波器到数字滤波器的转换。
实验内容
利用MATLAB编程，用脉冲响应不变法和双线性变换法设计一个数字带通滤波器，指标要求如下：
通带边缘频率：，，通带峰值起伏：；
阻带边缘频率：，，最小阻带衰减：。
代码实现：
双线性变换法（巴特沃斯原型）
functiontest2sb
ws1=2*8000*tan(0.3*pi/2);求边界频率
ws2=2*8000*tan(0.8*pi/2);
wp1=2*8000*tan(0.45*pi/2);
wp2=2*8000*tan(0.65*pi/2);
ws=[ws1ws2];wp=[wp1wp2];
Rp=1;Rs=40;
[N,Wn]=buttord(wp,ws,Rp,Rs,s);
[num,den]=butter(N,Wn,s);
%传输函数的分子分母系数
[B,A]=bilinear(num,den,8000);
%模拟滤波器到数字滤波器的变化
[h,w]=freqz(B,A);
f=w/pi*4000;
subplot(2,1,1);
plot(f,20*log10(abs(h)));
axis([0,4000,-60,10]);
grid;xlabel(频率/Hz);ylabel(幅度/dB);
subplot(2,1,2);
plot(f,angle(h));
grid;xlabel(频率/Hz);ylabel(相位);
实验截图：
(切比雪夫原型)
functiontest2sq
ws1=2*8000*tan(0.3*pi/2);
ws2=2*8000*tan(0.8*pi/2);
wp1=2*8000*tan(0.45*pi/2);
wp2=2*8000*tan(0.65*pi/2);
ws=[ws1ws2];wp=[wp1wp2];
Rp=1;Rs=40;
[N,Wn]=cheb1ord(wp,ws,Rp,Rs,s);
[num,den]=cheby1(N,Wn,s);%s为模拟频率要乘以pi
[B,A]=bilinear(num,den,8000);
[h,w]=freqz(B,A);
f=w/pi*4000;
subplot(2,1,1);
plot(f,20*log10(abs(h)));
axis([0,4000,-60,10]);
grid;xlabel(频率/Hz);ylabel(幅度/dB);
subplot(2,1,2);
plot(f,angle(h));
grid;xlabel(频率/Hz);ylabel(相位);
实验结果截图：
脉冲响应不变法（巴特沃斯原型）
functiontestmb
fs=8000;
ws1=0.3*pi*fs;ws2=0.8*pi*fs;
wp1=0.45*pi*fs;wp2=0.65*pi*fs;
ws=[ws1ws2];wp=[wp1wp2];
Rp=1;Rs=40;
[N,Wn]=buttord(wp,ws,Rp,Rs,s);
[num,den]=butter(N,Wn,s);
[B,A]=impinvar(num,den,8000);
%脉冲响应不变法模拟滤波到数字滤波的转变
[h,w]=freqz(B,A);
f=w/pi*4000;
subplot(2,1,1);
plot(f,20*log10(abs(h)));
axis([0,4000,-80,10]);
grid;xlabel(频率/Hz);ylabel(幅度/dB);
subplot(2,1,2);
plot(f,angle(h));
grid;xlabel(频率/Hz);ylabel(相位);
结果截图：
（切比雪夫原型）
functiontest2mq
fs=
文档评论（0）