數(shù)字信號處理實驗 濾波器在語音信號中的應(yīng)用 (1)

數(shù)字信號處理實驗姓名：彭金金 王春香學(xué)號： 濾波器在語音信號中的應(yīng)用一、實驗?zāi)康模?1、了解語音信號在MATLAB中的輸入和輸出； 2、通過對語音信號的處理進(jìn)一步了解濾波器的設(shè)計及應(yīng)用。二、實驗原理： 1、濾波器設(shè)計調(diào)用fir1（）函數(shù)生成有限長數(shù)字濾波器。2、采集語音信號在MATLAB軟件平臺下利用函數(shù)wavread對語音信號進(jìn)行采樣，得到聲音數(shù)據(jù)變量x，同時把x的采樣頻率fs和數(shù)據(jù)位n放進(jìn)MATLAB工作空間。x,fs,n=wavread('D:matas.wav');語音信號輸出指令為：sound（x，fs）；然后對語音信號進(jìn)行頻譜分析：X=fft（x）；3、加噪聲x1=x+randn(size(x);三、實驗內(nèi)容：程序如下：x,fs,n=wavread('D:matas.wav');t=0:0.0005:2;x1=x+randn(size(x);N=1:length(x);n1=0:400;M=64;subplot(3,1,1);plot(n1,x(3256+n1)title('原始信號x')subplot(3,1,2)plot(n1,x1(3256+n1)title('在原始信號上加噪聲')Y=abs(fft(x1);G=abs(fft(x);df=fs/length(x);wn=175 285 346 445/1000;b=fir1(M,wn);z=filter(b,length(x),x1);zk=abs(fft(z);subplot(3,1,3)plot(n1,z(3256+n1)title('濾波后的語音信號z')figure(2)n2=0:length(G)-1;n3=n2*df;subplot(3,1,1)plot(n3,G);title('原始信號頻譜')subplot(3,1,2)plot(n3,Y)title('加入噪聲后信號的頻譜')subplot(3,1,3)plot(n3,zk)title('濾波后信號的頻譜')sound(z,fs);輸出波形如下：四、實驗結(jié)論此次實驗中使用的FIR濾波器的濾波效果不是很好，可以根據(jù)原信號頻譜的范圍通過改變?yōu)V波器的參數(shù)改善濾波效果；axis(0,2500,0,5000)例如：當(dāng)wn=120 200 316 422/10000;時各信號頻譜為：當(dāng)wn=10 500 1000 1500/10000;時濾波后信號的頻譜如下：