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

《數(shù)字信號處理實(shí)驗(yàn) 濾波器在語音信號中的應(yīng)用 (1)》由會員分享，可在線閱讀，更多相關(guān)《數(shù)字信號處理實(shí)驗(yàn) 濾波器在語音信號中的應(yīng)用 (1)（5頁珍藏版）》請在裝配圖網(wǎng)上搜索。

1、數(shù)字信號處理實(shí)驗(yàn)姓名：彭金金 王春香學(xué)號： 濾波器在語音信號中的應(yīng)用一、實(shí)驗(yàn)?zāi)康模?1、了解語音信號在MATLAB中的輸入和輸出； 2、通過對語音信號的處理進(jìn)一步了解濾波器的設(shè)計及應(yīng)用。二、實(shí)驗(yàn)原理： 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(si

2、ze(x);三、實(shí)驗(yàn)內(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)

3、,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);輸出波形如下：四、實(shí)驗(yàn)結(jié)論此次實(shí)驗(yàn)中使用的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;時濾波后信號的頻譜如下：