原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验2数字信号处理
wp=2*1000*tan(2*pi*300/(2*1000));
ws=2*1000*tan(2*pi*200/(2*1000));
[N,wn]=cheb1ord(wp,ws,0.8,20,s);
[B,A]=cheby1(N,0.8,wn,s);
[bn1,an1]=bilinear(B,A,2000);
[h,w]=freqz(bn1,an1);
f=w/pi*500;plot(f,20*log10(abs(h)));
axis([0,500,-30,10]);
grid on;
xlabel(频率);
ylabel(幅度响应);
wp=2*1000*tan(2*pi*200/(2*1000));
ws=2*1000*tan(2*pi*300/(2*1000));
Ap=1;As=25;
[n,Wn]=buttord(Wp,Ws,Ap,As,s);
[b,a]=butter(n,Wn,s);
[bn1,an1]=bilinear(b,a,1000);
[H1,W]=freqz(bn1,an1);
[bn2,an2]=impinvar(b,a,1000);
[H2,W]=freqz(bn2,an2);
f=W/(2*pi)*1000;plot(f,abs(H1),-.,f,abs(H2),-);
grid on
xlabel(频率);ylabel(幅度响应);
legend(双线性变法,脉冲响应不变法);
wp=2*8000*tan(2*pi*1200/(2*8000));
ws=2*8000*tan(2*pi*2000/(2*8000));
[N1,wn1]=buttord(wp,ws,0.5,40,s);
[B1,A1]=butter(N1,wn1,s);
[bn1,an1]=bilinear(B1,A1,8000);
[h1,w1]=freqz(bn1,an1);
f1=w1/(2*pi)*8000;
subplot(311);plot(f1,20*log10(abs(h1)));
axis([0,4000,-80,10]);
grid on;
xlabel(频率1);ylabel(幅度响应1);title(butterworth滤波器);
[N2,wn2]=cheb1ord(wp,ws,0.5,40,s);
[B2,A2]=cheby1(N2,0.5,wn2,s);
[bn2,an2]=bilinear(B2,A2,8000);
[h2,w2]=freqz(bn2,an2);
f2=w2/(2*pi)*8000;
subplot(312);plot(f2,20*log10(abs(h2)));
axis([0,4000,-80,10]);
grid on;
xlabel(频率2);ylabel(幅度响应2);title(chebyshev滤波器);
[N3,wn3]=ellipord(wp,ws,0.5,40,s);
[B3,A3]=ellip(N3,0.5,40,wn3,s);
[bn3,an3]=bilinear(B3,A3,8000);
[h3,w3]=freqz(bn3,an3);
f3=w3/(2*pi)*8000;
subplot(313);plot(f3,20*log10(abs(h3)));
axis([0,4000,-80,10]);
grid on;
xlabel(频率3);ylabel(幅度响应3);title(椭圆滤波器);
w1=2*30000*tan(2*pi*2000/(2*30000));
w2=2*30000*tan(2*pi*3000/(2*30000));
wr1=2*30000*tan(2*pi*1500/(2*30000));
wr2=2*30000*tan(2*pi*6000/(2*30000));
[N,wn]=buttord([w1 w2],[wr1 wr2],3,20,s);
[B,A]=butter(N,wn,s);
[bn1,an1]=bilinear(B,A,30000);
[h,w]=freqz(bn1,an1);
f=w/(2*pi)*30000;plot(f,20*log10(abs(h)));
axis([0,15000,-80,10]);
grid on;
xlabel(频率3);ylabel(幅度响应3);
N=45，计算并画出矩形窗、汉明窗、布莱克曼窗的归一化的幅度谱，并比较各自的主要
特点。
N=45;
wn1=kaiser(N,0);
wn2=hamming(N);
wn3=blackman(N);
[h1,w1]=freqz(wn1,N);
[h2,w2]=freqz(wn2,N);
[h3,w3]=freqz(wn3,N);
plot(
文档评论（0）