- 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滤波器直接型与级联型的设计Matlab代码
clear;clc;
%%
[y,fs]=wavread(music1-2);
qq=16;
%% 设置滤波器参数
Wp=2*900;Ws=2*1500;Rp=0.1;Rs=60;Fs=32000;
wp=Wp/Fs;ws=Ws/Fs;
%% 选择滤波器类型
%[N,Wpo]=cheb2ord(wp,ws,Rp,Rs);%求切比雪夫2型归一化模拟低通滤波器的阶数
%[b,a]=cheby2(N,Rs,Wpo);%设计滤波器
[N,Wpo]=cheb1ord(wp,ws,Rp,Rs);%求切比雪夫1型归一化模拟低通滤波器的阶数
[b,a]=cheby1(N,Rp,Wpo);%设计滤波器
%[N,Wpo]=buttord(wp,ws,Rp,Rs);%求巴特沃斯型归一化模拟低通滤波器的阶数
%[b,a]=butter(N,Wpo);%设计滤波器
%[N,Wpo]=ellipord(wp,ws,Rp,Rs);%求椭圆型归一化模拟低通滤波器的阶数
%[b,a]=ellip(N,Rp,Rs,Wpo);%设计椭圆滤波器
%% 滤波器频响
[h,w]=freqz(b,a,1024);%求滤波器的幅频响应
h=20*log10(abs(h)/max(abs(h)));%幅度归一
figure;plot(w/2/pi*Fs,h);
title(滤波器的幅频响应);grid on;
bqr=intbR(b,qq);aqr=intbR(a,qq);
[hr,wr]=freqz(bqr,aqr,1024);
hr=20*log10(abs(hr)/max(abs(hr)));
%figure;plot((0:1023)/1024,hr);
figure;plot(w/2/pi*Fs,hr);
title(滤波器的量化后的幅频响应);grid on;
%% 级联处理
[sos,g]=tf2sos(b,a);
b1=sos(1,1:3);a1=sos(1,4:6);
b2=sos(2,1:3);a2=sos(2,4:6);
b3=sos(3,1:3);a3=sos(3,4:6);
b4=sos(4,1:3);a4=sos(4,4:6);
b5=sos(5,1:3);a5=sos(5,4:6);
% b6=sos(6,1:3);a6=sos(6,4:6);
% b7=sos(7,1:3);a7=sos(7,4:6);
% b8=sos(8,1:3);a8=sos(8,4:6);
% b9=sos(9,1:3);a9=sos(9,4:6);
%% 级联量化
bqr1=intbR(b1,qq); aqr1=intbR(a1,qq);
[h1,w1]=freqz(bqr1,aqr1,1024);
bqr2=intbR(b2,qq); aqr2=intbR(a2,qq);
[h2,w1]=freqz(bqr2,aqr2,1024);
bqr3=intbR(b3,qq); aqr3=intbR(a3,qq);
[h3,w1]=freqz(bqr3,aqr3,1024);
bqr4=intbR(b4,qq); aqr4=intbR(a4,qq);
[h4,w1]=freqz(bqr4,aqr4,1024);
bqr5=intbR(b5,qq); aqr5=intbR(a5,qq);
[h5,w1]=freqz(bqr5,aqr5,1024);
% bqr6=intbR(b6,qq); aqr6=intbR(a6,qq);
% [h6,w1]=freqz(bqr6,aqr6,1024);
% bqr7=intbR(b7,qq); aqr7=intbR(a7,qq);
% [h7,w1]=freqz(bqr7,aqr7,1024);
% bqr8=intbR(b3,qq); aqr8=intbR(a8,qq);
% [h8,w1]=freqz(bqr8,aqr8,1024);
% bqr9=intbR(b3,qq); aqr9=intbR(a9,qq);
%[h9,w1]=freqz(bqr9,aqr9,1024);
hh=h1.*h2.*h3;%.*h4;%.*h5.*h6.*h7;%.*h8.*h9;
hh=20*log10(abs(hh)/max(abs(hh)));
sos1=intbR(sos,qq);
%% 信号过系统
%yy=filter(bqr,aqr,y);
yy=filter(b,a,y);
Y=sosfilt(sos1,y);
%% 画图
figure;
%specgram(Y);
%sound(yy,fs);
%plot(20*log10
您可能关注的文档
最近下载
- 三峡左岸电站#~#机组压力钢管的制作与安装(二)的论文.doc VIP
- 网络安全集成项目投标方案(安全设计-实施方案-售后服务方案).docx
- 《面向智能仓储的AGV物流小车结构设计》8000字(论文).doc
- 日本横河WT210230说明书.pdf VIP
- 中国武器装备大全课件PPT.pptx
- 3500词(带音标)——完整打印版 .pdf VIP
- 七年级上册语文第一单元作业设计 (优质案例41页) .pdf VIP
- 2025年度检验科生物安全培训上岗考核试题(含答案).docx VIP
- DB62T 3031-2018 通风与空调工程施工工艺规程.docx VIP
- 《华为薪酬体系》课件.ppt VIP
文档评论(0)