基于matlab各滤波器源代码_第1页
基于matlab各滤波器源代码_第2页
基于matlab各滤波器源代码_第3页
基于matlab各滤波器源代码_第4页
基于matlab各滤波器源代码_第5页
已阅读5页，还剩14页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、汉明窗的FIR 低通滤波：Fs=22050;x,FS,bits=wavread('G:UsersDPDesktopSoundTest.wav');%G:UsersDPDesktop 好铃网-湖人掘金宣传片.wavx=x(:,1);figure(1);subplot(2,1,1);plot(x);%sound(x,FS,bits); % 回放语音title(' 语音信号时域波形图')y=fft(x,3260);f=(FS/3260)*1:1630;subplot(2,1,2);plot(f(1:1630),abs(y(1:1630);title(' 语音信

2、号频谱图');% 产生噪声信号并加到语音信号t=0:length(x)-1;zs0=0.05*cos(2*pi*100*t/22050);figure(2);subplot(2,1,1)plot(zs0)title(' 噪声信号波形');zs1=fft(zs0,1200);%sound(zs0,FS,bits);% 回放噪音f=(FS/1200)*1:600;subplot(2,1,2)plot(f(1:600),abs(zs1(1:600);title（' 噪声信号频谱'）;x1=x+zs0'%sound(x1,FS,bits); % 回放加入

3、噪声后的语音y1=fft(x1,1200);figure(3);subplot(2,1,1);plot(x1);f=(FS/1200)*1:600;subplot(2,1,2);plot(f(1:600),abs(y1(1:600);title(' 加入噪声后的信号波形');title(' 加入噪声后的信号频谱');% 加窗滤波wp=0.25*pi;ws=0.3*pi;wdelta=ws-wp;N=ceil(6.6*pi/wdelta); % 取整t=0:(size(x1)-1);wn=(0.2+0.3)*pi/2;b=fir1(N,wn/pi,hamming(

4、N+1); % 选择窗函数，并归一化截止频率f1=fftfilt(b,x1);figure(4)freqz(b,1,512)h1,w1=freqz(b,1);plot(w1*FS/(2*pi),20*log10(abs(h1);figure(5)subplot(2,1,1)plot(t,x1)title(' 滤波前的时域波形');subplot(2,1,3)plot(t,f1);title(' 滤波后的时域波形');sound(f1); % 播放滤波后的语音信号F0=fft(f1,1024);f=FS*(0:511)/1024;figure(6)y2=fft(x

5、1,1024);subplot(2,1,1);plot(f,abs(y2(1:512);% 画出滤波前的频谱图title(' 滤波前的频谱')xlabel('Hz');ylabel('fuzhi');subplot(2,1,2)F1=plot(f,abs(F0(1:512);% 画出滤波后的频谱图title(' 滤波后的频谱')xlabel('Hz');ylabel('fuzhi');巴特沃兹低通滤波器：Fs=22050;x,FS,bits=wavread('G:UsersDPDesktop

6、SoundTest.wav');%G:UsersDPDesktop 好铃网-湖人掘金宣传片.wavx=x(:,1);figure(1);subplot(2,1,1);plot(x);%sound(x,FS,bits); % 回放语音title(' 语音信号时域波形图')y=fft(x,3260);f=(FS/1630)*1:1630;subplot(2,1,2);plot(f(1:1630),abs(y(1:1630);title(' 语音信号频谱图');% 产生噪声信号并加到语音信号t=0:length(x)-1;zs=0.05*cos(2*pi*10

7、000*t/22050);zs0=0.05*cos(2*pi*10000*t/22050000);figure(2);subplot(2,1,1)plot(zs0)title(' 噪声信号波形');zs1=fft(zs,1200);%sound(zs,FS,bits); % 回放噪音subplot(2,1,2)plot(f(1:600),abs(zs1(1:600);title(' 噪声信号频谱');x1=x+zs'sound(x1,FS,bits); % 回放加入噪声后的语音y1=fft(x1,1200);figure(3);subplot(2,1,1

8、);plot(x1);title(' 加入噪声后的信号波形');subplot(2,1,2);plot(f(1:600),abs(y1(1:600);title(' 加入噪声后的信号频谱');% 低通滤波fp=3000;fs=3500;Fs=22050;rp=1;rs=10;wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;Fs1=1;wap=2*tan(wp/2);was=2*tan(ws/2);N,wc=buttord(wap,was,rp,rs,'s');B,A=butter(N,wc,'s');Bz,Az=bilin

9、ear(B,A,Fs1);figure(4);h,w=freqz(Bz,Az,512,Fs1*22050);plot(w,abs(h);title(' 巴特沃斯低通滤波器');xlabel(' 频率(HZ) ');ylabel(' 耗损(dB) ');grid on;yd=filter(Bz,Az,x1);figure(5);subplot(2,1,1);plot(yd); title(' 滤波后信号波形');ydd=fft(yd,1200);subplot(2,1,2);plot(f(1:600),abs(ydd(1:600);

10、title(' 滤波后信号频谱');%sound(yd,FS,bits);巴特沃兹高通滤波器：Fs=22050;x,FS,bits=wavread('G:UsersDPDesktopSoundTest.wav');x=x(:,1);figure(1);subplot(2,1,1);title(' 加入噪声后的信号频谱');plot(x);%sound(x,FS,bits); % 回放语音title(' 语音信号时域波形图')y=fft(x,3260);f=(FS/1630)*1:1630;subplot(2,1,2);plot(f

11、(1:1630),abs(y(1:1630);title(' 语音信号频谱图');% 产生噪声信号并加到语音信号t=0:length(x)-1;zs0=0.25*cos(2*pi*100*t/22050);figure(2);subplot(2,1,1)plot(zs0)title(' 噪声信号波形');zs1=fft(zs0,1200);%sound(zs0,FS,bits); % 回放噪音subplot(2,1,2)plot(f(1:600),abs(zs1(1:600);title(' 噪声信号频谱');x1=x+zs0'%soun

12、d(x1,FS,bits); % 回放加入噪声后的语音y1=fft(x1,1200);figure(3);subplot(2,1,1);plot(x1);title(' 加入噪声后的信号波形');subplot(2,1,2);plot(f(1:600),abs(y1(1:600);% 高通滤波fp=400;fs=300;Fs=22050;rp=1;rs=10;wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;T=1;Fs1=1;wap=2*tan(wp/2);was=2*tan(ws/2);N,wc=buttord(wap,was,rp,rs,'s');B

13、,A=butter(N,wc,'high','s');Bz,Az=bilinear(B,A,Fs1);figure(4);h,w=freqz(Bz,Az,512,Fs1*22050);plot(w,abs(h);title(' 巴特沃斯高通滤波器');xlabel(' 频率(HZ) ');ylabel(' 耗损(dB) ');grid on;yd=filter(Bz,Az,x1);figure(5);subplot(2,1,1);plot(yd); title(' 滤波后信号波形');ydd=fft

14、(yd,1200);subplot(2,1,2);plot(f(1:600),abs(ydd(1:600);title(' 滤波后信号频谱');%sound(yd,FS,bits)巴特沃兹带通滤波器：Fs=22050;x,FS,bits=wavread('G:UsersDPDesktopSoundTest.wav');x=x(:,1);figure(1);subplot(2,1,1);plot(x);%sound(x,FS,bits); % 回放语音title(' 语音信号时域波形图')y=fft(x,3260);f=(FS/1630)*1:16

15、30;subplot(2,1,2);plot(f(1:1630),abs(y(1:1630);title(' 语音信号频谱图');% 产生噪声信号并加到语音信号t=0:length(x)-1;zs0=0.05*cos(2*pi*100*t/22050);figure(2);subplot(2,1,1)plot(zs0)title(' 噪声信号波形');zs1=fft(zs0,1200);%sound(zs,FS,bits); % 回放噪音subplot(2,1,2)plot(f(1:600),abs(zs1(1:600);title(' 噪声信号频谱&#

16、39;);x1=x+zs0'%sound(x1,FS,bits); % 回放加入噪声后的语音y1=fft(x1,1200);figure(3);subplot(2,1,1);plot(x1);subplot(2,1,2);plot(f(1:600),abs(y1(1:600);title(' 加入噪声后的信号频谱');% 带通滤波fp=800,9000;fs=700,10000;Fs=22050;rp=1;rs=10;wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;T=1;Fs1=1;wap=2*tan(wp/2);was=2*tan(ws/2);N,wc=bu

17、ttord(wap,was,rp,rs,'s');B,A=butter(N,wc,'s');Bz,Az=bilinear(B,A,Fs1);figure(4);h,w=freqz(Bz,Az,512,Fs1*22050);plot(w,abs(h);title(' 巴特沃斯带通滤波器');xlabel(' 频率(HZ) ');ylabel(' 耗损( dB) ');grid on;yd=filter(Bz,Az,x1);figure(5);subplot(2,1,1);plot(yd);ydd=fft(yd,120

18、0);subplot(2,1,2);plot(f(1:600),abs(ydd(1:600);%sound(yd,FS,bits)Blackman 加窗高通滤波器：Fs=22050;x,FS,bits=wavread('G:UsersDPDesktopSoundTest.wav');x=x(:,1);figure(1);subplot(2,1,1);plot(x);%sound(x,FS,bits); % 回放语音title(' 语音信号时域波形图')y=fft(x,3260);f=(FS/1630)*1:1630;subplot(2,1,2);plot(f(1:1630),abs(y(1:1630);title(' 语音信号频谱图');% 产生噪声信号并加到语音信号% 产生噪声信号并加到语音信号t=0:length(x)-1;zs0=0.05*cos(2*pi*100*t/22050);figure(2);subplot(2,1,1)plot(zs0)title(' 噪声信号波形');zs1=fft(zs

人人文库> 全部分类> 行业资料 > 信息产业

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论