《數(shù)字信號(hào)處理實(shí)驗(yàn) 濾波器在語(yǔ)音信號(hào)中的應(yīng)用 (1)》由會(huì)員分享，可在線閱讀，更多相關(guān)《數(shù)字信號(hào)處理實(shí)驗(yàn) 濾波器在語(yǔ)音信號(hào)中的應(yīng)用 (1)（5頁(yè)珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。
1、數(shù)
字
信
號(hào)
處
理
實(shí)
驗(yàn)
姓名：彭金金 王春香
學(xué)號(hào)：
濾波器在語(yǔ)音信號(hào)中的應(yīng)用
一、實(shí)驗(yàn)?zāi)康模?
1、了解語(yǔ)音信號(hào)在MATLAB中的輸入和輸出；
2、通過(guò)對(duì)語(yǔ)音信號(hào)的處理進(jìn)一步了解濾波器的設(shè)計(jì)及應(yīng)用。
二、實(shí)驗(yàn)原理：
1、濾波器設(shè)計(jì)
調(diào)用fir1（）函數(shù)生成有限長(zhǎng)數(shù)字濾波器。
2、采集語(yǔ)音信號(hào)
在MATLAB軟件平臺(tái)下利用函數(shù)wavread對(duì)語(yǔ)音信號(hào)進(jìn)行采樣，得到聲音數(shù)據(jù)變量x，同時(shí)把x的采樣頻率fs和數(shù)據(jù)位n放進(jìn)MATLAB工作空間。
[x,fs,n]=wavread('D:\mat\as.wav')
2、;
語(yǔ)音信號(hào)輸出指令為：sound（x，fs）；
然后對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析：X=fft（x）；
3、加噪聲
x1=x+randn(size(x));
三、實(shí)驗(yàn)內(nèi)容：
程序如下：
[x,fs,n]=wavread('D:\mat\as.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('原始信號(hào)x')
subplot(3,1,2)
plot(n1,x1(3256+n1))
ti
3、tle('在原始信號(hào)上加噪聲')
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('濾波后的語(yǔ)音信號(hào)z')
figure(2)
n2=0:length(G)-1;
n3=n2*df;
subplot(3,1,1)
plot(n3,G);
title('原始信號(hào)頻譜')
4、subplot(3,1,2)
plot(n3,Y)
title('加入噪聲后信號(hào)的頻譜')
subplot(3,1,3)
plot(n3,zk)
title('濾波后信號(hào)的頻譜')
sound(z,fs);
輸出波形如下：
四、實(shí)驗(yàn)結(jié)論
此次實(shí)驗(yàn)中使用的FIR濾波器的濾波效果不是很好，可以根據(jù)原信號(hào)頻譜的范圍通過(guò)改變?yōu)V波器的參數(shù)改善濾波效果；
axis([0,2500,0,5000])
例如：當(dāng)wn=[120 200 316 422]/10000;時(shí)各信號(hào)頻譜為：
當(dāng)wn=[10 500 1000 1500]/10000;時(shí)濾波后信號(hào)的頻譜如下：
數(shù)字信號(hào)處理實(shí)驗(yàn) 濾波器在語(yǔ)音信號(hào)中的應(yīng)用 (1)
