关于课程设计1.docVIP

下载本文档

8
0
约2.73千字
约 7页
2018-02-04 发布于重庆
举报
版权申诉

关于课程设计1.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文档。上传文档

关于课程设计1

巴特沃斯模拟地通数字滤波器 %声音信号 [x,fs,bits]=wavread(D:\wode\shangbuqi);%采集语音信号 sound(x,fs,bits); N=length(x); X=fft(x,N);%求原始信号频谱 n=0:N-1; subplot(421);plot(x);grid; title(原波形图); subplot(422);plot(abs(X));grid; title(原频谱图) %加噪音 y=x+0.01rand(N,1);%加信噪比为13的白噪声 sound(y,fs,bits); subplot(423);plot(y);grid; title(加噪声后的波形图); Y=fft(y,N);%对加噪后的信号求频谱 subplot(424);plot(abs(Y));grid; title(加噪声后频谱图); %滤波器 Fs=8000; %滤波器指标 fp=1200;wp=2*pi*fp/Fs; fs=1600;ws=2*pi*fs/Fs; wap=2*Fs*tan(wp/2); was=2*Fs*tan(ws/2); Rp=1;As=100; [M,wc]=buttord(wap,was,Rp,As,s);%求阶数 [z,p,k]=buttap(M);%求零极点 [bp,ap]=zp2tf(z,p,k);%求模拟滤波器系数 [bs,as]=lp2lp(bp,ap,wap);%去归一化 [bz,az]=bilinear(bs,as,Fs/2);%用双线性变换法转为数字系数 [h,w]=freqz(bz,az);%计算数字滤波器频响 subplot(425); plot(w,abs(h));grid; title(巴特沃斯数字低通滤波器特性); %滤波 l=filter(bz,az,y);%对原声信号进行过滤 sound(l,fs,bits); subplot(426);plot(l);grid; title(滤波后波形图); L=fft(l,N);%对过滤后的声音信号进行傅里叶变换求幅频响应 subplot(427);plot(abs(L));grid; title(滤波后图谱); 二．椭圆数字高通滤波器 %声音信号 [x,fs,bits]=wavread(D:\wode\shangbuqi);%采集原始语音信号 sound(x,fs,bits); N=length(x); X=fft(x,N);%对原始信号进行傅里叶变换求频谱 n=0:N-1; subplot(421); plot(x);grid; title(原信号波形图); subplot(422); plot(abs(X));grid; title(原信号频谱图) %加噪音 y=x+0.01*rand(N,1);%对原始语音信号加噪 sound(y,fs,bits); subplot(423); plot(y);grid; title(加噪声后的波形图); Y=X+0.01*(N,1); subplot(424); plot(abs(Y));grid; title(加噪声后频谱图); %加滤波器 fp=5000;Rp=1;%指标 fs=4800;As=100; Fs=16000; wp=2*fp/Fs; ws=2*fs/Fs; [n,wn]=ellipord(wp,ws,Rp,As);%求椭圆数字滤波器阶数和通带截止频率 subplot(425); [B,A]=ellip(n,Rp,As,wn,high);%求滤波器系统函数的系数 freqz(B,A);%幅频响应 %滤波 h=filter(B,A,y);%对原始信号滤波 sound(h,fs,bits); subplot(426); plot(h);grid; title(滤波后信号波形图); H=fft(h,N);%对滤波后的原始信号求频响 subplot(427); plot(abs(H));grid; title(滤波后频谱图); 声音信号三．布莱克曼窗函数滤波器 %声音信号 [x,fs,bits]=wavread(D:\wode\shangbuqi);%采集语音信号 sound(x,fs,bits); N=length(x); X=fft(x,N);%原信号求频谱 n=0:N-1; subplot(421); plot(x);grid; title(原信号波形图); subplot(422); plot(abs(X));grid; title(原信号频谱图) %加噪音 y=x+0.01rand(N,1);%原信号加噪 sound(y,fs,bits); subplot(423); plot(y);grid; title(加噪声后的波形图);