版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、数字信号处理:已知通带截止频率fp=5kHz,通带最大衰减ap=2dB,阻带截止频率fs=2kHz,阻带最小衰减as=30dB,按照以上技术指标设计巴特沃斯低通滤波器:wp=2*pi*5000;ws=2*pi*12000;Rp=2;As=30;N,wc=buttord(wp,ws,Rp,As,s);B,A=butter(N,wc,s);k=0:511;fk=0:14000/512:14000;wk=2*pi*fk;Hk=freqs(B,A,wk);subplot(2,2,1);plot(fk/1000,20*log10(abs(Hk);grid onxlabel(频率(kHz));ylabel
2、(幅度(dB))axis(0,14,-40,5)切比雪夫1型低通滤波器:wp=2*pi*3000;ws=2*pi*12000;Rp=0.1;As=60;N1,wpl=cheb1ord(wp,ws,Rp,As,s); %cheb1ord,里面的是1,不是LB1,A1=cheby1(N1,Rp,wpl,s);subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;Hk=freqs(B1,A1,wk);plot(fk/1000,20*log10(abs(Hk);grid onxlabel(频率(kHz));ylabel(幅度(dB))axis(0,12,-70,
3、5)椭圆模拟低通滤波器:wp=2*pi*3000;ws=2*pi*12000;Rp=0.1;As=60;N,wpo=ellipord(wp,ws,Rp,As,s);B,A=ellip(N,Rp,As,wpo,s);subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;Hk=freqs(B1,A1,wk);plot(fk/1000,20*log10(abs(Hk);grid onxlabel(频率(kHz));ylabel(幅度(dB))axis(0,12,-70,5)p195-14wp=2*4/80;ws=2*20/80;rp=0.5;rs=45;N,w
4、c=buttord(wp,ws,rp,rs);B,A=butter(N,wc);hk,wk=freqz(B,A);fk=wk/pi*40;plot(fk,20*log10(abs(hk);axis(0,30,-100,0)xlabel(频率(kHZ);ylabel(幅度(db);grid onP195-16wp=2*325/2500;ws=2*225/2500;rp=1;rs=40;N,wc=ellipord(wp,ws,rp,rs);B,A=ellip(N,rp,rs,wc);hk,wk=freqz(B,A);fk=wk/pi*40;plot(fk,20*log10(abs(hk);axis
5、(0,30,-100,0)xlabel(频率(kHZ);ylabel(幅度(db);grid onP195-15wp=2*4/80;ws=2*20/80;rp=0.5;rs=45;N,wc=cheb1ord(wp,ws,rp,rs);B,A=cheby1(N,rp,wc);hk,wk=freqz(B,A);fk=wk/pi*40;plot(fk,20*log10(abs(hk);axis(0,30,-100,0)xlabel(频率(kHZ);ylabel(幅度(db);grid on切比雪夫低通滤波器wp=2*pi*3000;ws=2*pi*12000;rp=0.1;as=60;N1,wp1=
6、cheb1ord(wp,ws,rp,as,s);B1,A1=cheby1(N1,rp,wp1,s);subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;hk=freqs(B1,A1,wk);plot(fk/1000,20*log10(abs(hk);grid onxlabel(频率(kHZ);ylabel(幅度(db);axis(0,12,-70,5)双音频检测audiofile=test.wavin_audio,fs,bits=wavread(audiofile);b,a=cheby1(5,0.1,0.3);out_audio=filter(b,a,
7、in_audio);sound(out_audio,fs,bits);wavwrite(out_audio,fs,bits,test_out);xk1=fft(in_audio,512);xk2=fft(out_audio,512);subplot(2,1,1);stem(abs(xk1);subplot(2,1,2);stem(abs(xk2);巴特沃斯模拟高通滤波器。wp=2*pi*4000;ws=2*pi*1000;rp=0.1;as=40;n,wc=buttord(wp,ws,rp,as,s);bh,ah=butter(n,wc,high,s);subplot(2,2,1);fk=0:
8、12000/512:12000;wk=2*pi*fk;hk=freqs(bh,ah,wk);plot(fk/1000,20*log10(abs(hk);grid onxlabel(频率(khz);ylabel(幅度(db)axis(0,12,-70,5)巴特沃斯模拟低通滤波器。wp=2*pi*5000;ws=2*pi*12000;rp=2;as=30;n,wc=buttord(wp,ws,rp,as,s);b,a=butter(n,wc,s);k=0:511;fk=0:14000/512:14000;wk=2*pi*fk;hk=freqs(b,a,wk);subplot(2,2,1);plot
9、(fk/1000,20*log10(abs(hk);grid onxlabel(频率(khz);ylabel(幅度(db)axis(0,14,-40,5)巴特沃斯模拟带通滤波器wp=2*pi*4000,7000;ws=2*pi*2000,9000;n,wc=buttord(wp,ws,rp,as,s);bb,ab=butter(n,wc,s);subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;hk=freqs(bb,ab,wk);plot(fk/1000,20*log10(abs(hk);grid onxlabel(频率(khz);ylabel(幅度
10、(db)axis(0,12,-70,5)巴特沃斯带阻滤波器wp=2*pi*2000,9000;ws=2*pi*4000,7000;rp=1;as=20;nb,wc=buttord(wp,ws,rp,as,s);bsb,asb=butter(nb,wc,stop,s);subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;hk=freqs(bsb,asb,wk);plot(fk/1000,20*log10(abs(hk);grid onxlabel(频率(khz);ylabel(幅度(db)axis(0,12,-70,5)第14题。巴特沃斯低通滤波器(采样
11、频率80,同代便捷频率4,最大衰减0.5)wp=2*4/80;ws=2*20/80;rp=0.5;rs=45; N,wc=buttord(wp,ws,rp,rs);B,A=butter(N,wc);hk,wk=freqz(B,A);fk=wk/pi*40;plot(fk,20*log10(abs(hk);axis(0,30,-100,0);xlabel(频率(kHz);ylabel(幅度(dB)设计切比雪夫I型和切比雪夫II型模拟低通滤波器。 I型:wp=2*pi*3000;ws=2*pi*12000;rp=0.1;as=60;n1,wp1=cheb1ord(wp,ws,rp,as,s);b1
12、,a1=cheby1(n1,rp,wp1,s);subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;hk=freqs(b1,a1,wk);plot(fk/1000,20*log10(abs(hk);grid onxlabel(频率(khz);ylabel(幅度(db)axis(0,12,-70,5)II型:wp=2*pi*3000;ws=2*pi*12000;rp=0.1;as=60;n2,wp2=cheb2ord(wp,ws,rp,as,s);b2,a2=cheby2(n2,rp,wp2,s);subplot(2,2,1);fk=0:12000/51
13、2:12000;wk=2*pi*fk;hk=freqs(b2,a2,wk);plot(fk/1000,20*log10(abs(hk);grid onxlabel(频率(khz);ylabel(幅度(db)axis(0,12,-70,5)切比雪夫I型高通滤波器wp=2*pi*3000;ws=2*pi*12000;rp=0.1;as=60;n1,wp1=cheb1ord(wp,ws,rp,as,s);b1,a1=cheby1(n1,rp,wp1,HIGH,s);subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;hk=freqs(b1,a1,wk);pl
14、ot(fk/1000,20*log10(abs(hk);grid onxlabel(频率(khz);ylabel(幅度(db)axis(0,12,-70,5)设计椭圆模拟低通滤波器。wp=2*pi*3000;ws=2*pi*12000;rp=0.1;as=60;n,wpo=ellipord(wp,ws,rp,as,s);b,a=ellip(n,rp,as,wpo,s);subplot(2,2,1);fk=0:12000/512:12000;wk=2*pi*fk;hk=freqs(b,a,wk);plot(fk/1000,20*log10(abs(hk);grid onxlabel(频率(khz);ylabel(幅度(db)axis(0,12,-70,5)椭圆带阻滤波器wp=2*pi*2000,9000;ws=2*pi*4000,7000;r
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 数据湖建设项目解决方案
- (10)一元一次不等式和一元一次不等式组(B卷)-七年级数学冀教版暑假作业(含答案)
- 初中+物理电功率(第二课时)课件+人教版物理九年级全一册物理
- 大数据背景下企业财务风险防范
- 《1 例短肠综合征患者置入 PICC 导管并发血栓的护理》
- 大学美育 课件 第三篇 第三章 第三节 影视求真
- 物质的分类及转化 同步练习 高一上学期化学人教版(2019)必修第一册+
- Unit+1+A+New+Start+语法基础练 高中英语外研版(2019)必修第一册
- 高等数学(第五版)课件 5.1 定积分的概念与性质
- 教师诚信演讲稿范文
- 某工程室外箱式变电站安装施工方案完整
- 新加坡建筑业管理体制
- 社团活动记录手册(共14页)
- 取水头部及原水自流管专项施工方案(共56页)
- 常见材料耐腐蚀性表
- 推广普通话.写好规范字-主题班会_ppt课件
- 会计信息系统教学实用教案
- 直线的法向量和点法式方程演示教学
- 薪酬制度的发展演变
- 水文水井钻机安全操作规程示范文本
- 肺肿瘤生活质量调查表
评论
0/150
提交评论