語(yǔ)音處理課程設(shè)計(jì)

1、mainb=menu('請(qǐng)選擇選項(xiàng):','IIR濾波器','FIR濾波器','退出');if b=3 b=0;endif b=1 temp=menu('請(qǐng)選擇選項(xiàng):','IIR低通濾波器','IIR高通濾波器','IIR帶通濾波器','返回'); if temp=1 IIR_LP; main; elseif temp=2 IIR_HP; main; elseif temp=3 IIR_BP; main; else main; endelseif b=2

2、 temp=menu('請(qǐng)選擇選項(xiàng):','FIR低通濾波器','FIR高通濾波器','FIR帶通濾波器','返回'); if temp=1 FIR_LP; main; elseif temp=2 FIR_HP; main; elseif temp=3 FIR_BP; main; else main; endendIIRLPfs=22050;x1,fs,bits=wavread('d:徐亮.wav'); sound(x1,fs);x2=wavread('d:XP.wav');Ts=1/

3、fs;R1=10;wp=2*pi*1000/fs;ws=2*pi*1200/fs;Rp=1;Rl=100;wp1=2/Ts*tan(wp/2); %將模擬指標(biāo)轉(zhuǎn)換成數(shù)字指標(biāo)ws1=2/Ts*tan(ws/2); N,Wn=buttord(wp1,ws1,Rp,R1,'s'); %選擇濾波器的最小階數(shù)Z,P,K=buttap(N); %創(chuàng)建butterworth模擬濾波器Bap,Aap=zp2tf(Z,P,K);b,a=lp2lp(Bap,Aap,Wn); bz,az=bilinear(b,a,fs); %用雙線性變換法實(shí)現(xiàn)模擬濾波器到數(shù)字濾波器的轉(zhuǎn)換H,W=freqz(bz,a

4、z); %繪制頻率響應(yīng)曲線figure(1)plot(W*fs/(2*pi),abs(H)axis(0 3000 0 1.4)gridxlabel('頻率Hz')ylabel('頻率響應(yīng)幅度')title('IIR低通濾波器')f1=filter(bz,az,x2);figure(2)subplot(2,2,1)plot(x2) %畫出濾波前的時(shí)域圖title('IIR低通濾波器濾波前的時(shí)域波形');subplot(2,2,3)plot(f1); %畫出濾波后的時(shí)域圖title('IIR低通濾波器濾波后的時(shí)域波形'

5、);sound(f1,44100); %播放濾波后的信號(hào)F0=fft(f1,1024);f=fs*(0:511)/1024;y2=fft(x2,1024);subplot(2,2,2);plot(f,abs(y2(1:512); %畫出濾波前的頻譜圖axis(0 10000 0 150)title('IIR低通濾波器濾波前的頻譜')xlabel('頻率/Hz');ylabel('幅值')subplot(2,2,4)F1=plot(f,abs(F0(1:512); %畫出濾波后的頻譜圖axis(0 5000 0 150)title('IIR

6、低通濾波器濾波后的頻譜')xlabel('頻率/Hz');ylabel('幅值');IIRHPclear;close allx1,fs,bits=wavread('d:徐亮.wav'); sound(x1,fs)fs=22050;x1=wavread('d:XP.wav');Fs=22050;Ts=1/Fs;R1=50;Wp=2*pi*5000/fs;Ws=2*pi*4800/fs;Rp=1;Rl=100;Wp1=2/Ts*tan(Wp/2); %將模擬指標(biāo)轉(zhuǎn)換成數(shù)字指標(biāo)Ws1=2/Ts*tan(Ws/2); N,Wn=c

7、heb2ord(Wp1,Ws1,Rp,Rl,'s'); %選擇濾波器的最小階數(shù)Z,P,K=cheb2ap(N,Rl); %創(chuàng)建切比雪夫模擬濾波器Bap,Aap=zp2tf(Z,P,K);b,a=lp2hp(Bap,Aap,Wn); bz,az=bilinear(b,a,Fs); %用雙線性變換法實(shí)現(xiàn)模擬濾波器到數(shù)字濾波器的轉(zhuǎn)換H,W=freqz(bz,az); %繪制頻率響應(yīng)曲線figure(1)plot(W*fs/(2*pi),abs(H)gridxlabel('頻率Hz')ylabel('頻率響應(yīng)幅度')title('IIR高通濾波器

8、')f1=filter(bz,az,x1);figure(2)subplot(2,2,1)plot(x1) %畫出濾波前的時(shí)域圖title('IIR高通濾波器濾波前的時(shí)域波形');subplot(2,2,3)plot(f1); %畫出濾波后的時(shí)域圖title('IIR高通濾波器濾波后的時(shí)域波形');sound(f1,44100); %播放濾波后的信號(hào)F0=fft(f1,1024);f=fs*(0:511)/1024;y2=fft(x1,1024);subplot(2,2,2);plot(f,abs(y2(1:512); %畫出濾波前的頻譜圖axis(0

9、5000 0 150)title('IIR高通濾波器濾波前的頻譜')xlabel('頻率/Hz');ylabel('幅值');subplot(2,2,4)plot(f,abs(F0(1:512); %畫出濾波后的頻譜圖axis(3000 10000 0 3)title('IIR高通濾波器濾波后的頻譜')xlabel('頻率/Hz');ylabel('幅值');IIRBPfs=22050;x1,fs,bits=wavread('d:徐亮.wav'); sound(x1,fs);x1=w

10、avread('d:XP.wav');Fs=22050;Ts=1/Fs;R1=30;fb1=1200;fb2=3000;fc1=1000;fc2=3200;fs=22050;W1=2*fb1*pi/fs;W2=2*fc1*pi/fs;W3=2*fb2*pi/fs;W4=2*fc2*pi/fs;Wp=W1,W3;Ws=W2,W4;Rp=1;Rl=100;Wp1=2/Ts*tan(Wp/2); %將模擬指標(biāo)轉(zhuǎn)換成數(shù)字指標(biāo)Ws1=2/Ts*tan(Ws/2); N,Wn=cheb2ord(Wp1,Ws1,Rp,R1,'s'); %選擇濾波器的最小階數(shù)Z,P,K=che

11、b2ap(N,Rl); %創(chuàng)建切比雪夫模擬濾波器Bap,Aap=zp2tf(Z,P,K);b,a=lp2bp(Bap,Aap,2100*2*pi,1800*2*pi); bz,az=bilinear(b,a,Fs); %用雙線性變換法實(shí)現(xiàn)模擬濾波器到數(shù)字濾波器的轉(zhuǎn)換H,W=freqz(bz,az); %繪制頻率響應(yīng)曲線figure(1)plot(W*fs/(2*pi),abs(H)axis(500 4000 0 1.4)gridxlabel('頻率Hz')ylabel('頻率響應(yīng)幅度')title('IIR帶通濾波器')f1=filter(bz,

12、az,x1);figure(2)subplot(2,2,1)plot(x1) %畫出濾波前的時(shí)域圖title('IIR帶通濾波器濾波前的時(shí)域波形');subplot(2,2,3)plot(f1); %畫出濾波后的時(shí)域圖title('IIR帶通濾波器濾波后的時(shí)域波形');sound(f1,44100); %播放濾波后的信號(hào)F0=fft(f1,1024);f=fs*(0:511)/1024;y2=fft(x1,1024);subplot(2,2,2);plot(f,abs(y2(1:512);%畫出濾波前的頻譜圖axis(0 5000 0 120)title(

13、9;IIR帶通濾波器濾波前的頻譜')xlabel('頻率/Hz');ylabel('幅值');subplot(2,2,4)plot(f,abs(F0(1:512); %畫出濾波后的頻譜圖axis(0 5000 0 12)title('IIR帶通濾波器濾波后的頻譜')xlabel('頻率/Hz');ylabel('幅值');FIRLPx1,fs,bits=wavread('d:徐亮.wav'); sound(x1,fs);fs=10000;x1=wavread('d:XP.wav

14、9;);wp=2*pi*1000/fs;ws=2*pi*1200/fs;Rp=1;Rs=100;wdelta=ws-wp;N=ceil(8*pi/wdelta); %取整wn=(wp+ws)/2;b,a=fir1(N,wn/pi,kaiser(N+1); %選擇窗函數(shù)，并歸一化截止頻率figure(1)freqz(b,a,512);title('FIR低通濾波器');f2=filter(b,a,x1);figure(2)subplot(2,2,1)plot(x1)title('FIR低通濾波器濾波前的時(shí)域波形');subplot(2,2,3)plot(f2);t

15、itle('FIR低通濾波器濾波后的時(shí)域波形');sound(f2,44100); %播放濾波后的語(yǔ)音信號(hào)F0=fft(f2,1024);f=fs*(0:511)/1024;y2=fft(x1,1024);subplot(2,2,2);plot(f,abs(y2(1:512);%axis(0 2000 0 150)title('FIR低通濾波器濾波前的頻譜')xlabel('頻率/Hz');ylabel('幅值');subplot(2,2,4)F2=plot(f,abs(F0(1:512);%axis(0 1500 0 150)t

16、itle('FIR低通濾波器濾波后的頻譜')xlabel('頻率/Hz');ylabel('幅值');FIRHPx1,fs,bits=wavread('d:徐亮.wav'); sound(x1,fs);fs=22050;x1=wavread('d:XP.wav');wp=2*pi*5000/fs;ws=2*pi*4800/fs;Rp=1;Rs=100;wdelta=wp-ws;N=ceil(8*pi/wdelta); %取整wn=(wp+ws)/2;b,a=fir1(N,wn/pi,'high');

17、 figure(1)freqz(b,a,512);title('FIR高通濾波器');f2=filter(b,a,x1);figure(2)subplot(2,2,1)plot(x1)title('FIR高通濾波器濾波前的時(shí)域波形');subplot(2,2,3)plot(f2);title('FIR高通濾波器濾波后的時(shí)域波形');sound(f2,44100); %播放濾波后的語(yǔ)音信號(hào)F0=fft(f2,1024);f=fs*(0:511)/1024;y2=fft(x1,1024);subplot(2,2,2);plot(f,abs(y2(1:

18、512);axis(0 4000 0 150)title('FIR高通濾波器濾波前的頻譜')xlabel('頻率/Hz');ylabel('幅值');subplot(2,2,4)F2=plot(f,abs(F0(1:512);axis(4000 10000 0 1.5)title('FIR高通濾波器濾波后的頻譜')xlabel('頻率/Hz');ylabel('幅值');FIRBPx1,fs,bits=wavread('d:徐亮.wav'); sound(x1,fs);fs=22050;x1=wavread('d:徐亮.wav');wp1=2*pi*1200/fs;wp2=2*pi*3000/fs;ws1=2*pi*1000/fs;ws2=2*pi*3200/fs;Rp=1;Rs=100;wp=(wp1+ws1)/2;ws=(wp2+ws2)