設(shè)計(jì)數(shù)字低通IIR和FIR濾波器對語音信號(hào)進(jìn)行濾波處理
- 文件介紹：
- 該文件為 doc 格式，下載需要 0 積分
- 設(shè)計(jì)數(shù)字低通IIR和FIR濾波器對語音信號(hào)進(jìn)行濾波處理
一．實(shí)驗(yàn)?zāi)康?br />
1.鞏固所學(xué)的數(shù)字信號(hào)處理理論知識(shí)， 理解信號(hào)的采集、處理、傳輸、顯示和存儲(chǔ)過程；
2.綜合運(yùn)用專業(yè)及基礎(chǔ)知識(shí)，解決實(shí)際工程技術(shù)問題的能力；
3.學(xué)習(xí)資料的收集與整理，學(xué)會(huì)撰寫課程設(shè)計(jì)報(bào)告。
二．實(shí)驗(yàn)內(nèi)容
1.選擇一個(gè)語音信號(hào)作為分析的對象，對其進(jìn)行頻譜分析；
2.設(shè)計(jì)FIR和IIR數(shù)字濾波器， 并對加噪語音信號(hào)進(jìn)行濾波，分析濾波后信號(hào)的時(shí)域和頻域特征。
三．設(shè)計(jì)過程
1.原語音信號(hào)的時(shí)域，頻域圖
[y,fs,nbits]=wavread ('d:\1\liuwei.wav'); %IIR低通sound(y,fs,bits);
figure(1);
plot(y);%做原始語音信號(hào)的時(shí)域圖形
title('原始語音信號(hào)');
xlabel('時(shí)間 t');
ylabel('音量 n');
figure(2);
y1=fft(y);
y1=fftshift(y1); %平移，是頻率中心為0
derta_fs = fs/length(y); %設(shè)置頻譜的間隔，分辨率
plot([-fs/2:derta_fs: fs/2-derta_fs],abs(y1));%畫出原始語音信號(hào)的頻譜圖
title('原始語音信號(hào)的頻譜');
2.低通濾波器的設(shè)計(jì)
設(shè)計(jì)指標(biāo)：
fp=1000Hz，fs=1200Hz，As=100db ,Ap=1dB
（1）低通IIR濾波器
Ft=8000;
Fp=1000;
Fs=1200;
wp=2*pi*Fp/Ft;
ws=2*pi*Fs/Ft;
fp=2*Ft*tan(wp/2);
fs=2*Fs*tan(wp/2);
[n11,wn11]=buttord(wp,ws,1,50, 's');%求低通濾波器的階數(shù)和截止頻率
[b11,a11]=butter(n11,wn11, 's'); %求S域的頻率響應(yīng)的參數(shù)
[num11,den11]=bilinear(b11,a11,0.5); %利用雙線性變換實(shí)現(xiàn)頻率響應(yīng)S域到Z域的變換
[h,w]=freqz(num11,den11);
figure(3);
plot(w*8000*0.5/pi,abs(h));
legend('IIR低通濾波器','Location','NorthWest');
grid;
z11=filter(num11,den11,y);
sound(z11);
m11=fft(z11); %求濾波后的信號(hào)...
熱門關(guān)鍵詞： 
