matlab語音信號采集課程設(shè)計-wenkub.com
2025-01-13 05:28 本頁面
【正文】 以后不斷學(xué)習(xí)，不斷完善自己。所以根據(jù)各個濾波器的選頻范圍及濾波后信號的音效，可大致判斷原信號的頻率多集中在1000Hz以下。濾波器的作用就是濾除掉阻帶內(nèi)有頻率分量，保留通帶內(nèi)的頻率分量。通過老師所給的材料，結(jié)合自己所學(xué)的數(shù)字信號處理教程，以及通過網(wǎng)上查找的一些資料，在兩周時間內(nèi)的學(xué)習(xí)，使得自己能夠在本次課程設(shè)計中較好的完成老師分配的任務(wù)。信號經(jīng)過IIR帶通濾波器（頻域）39。)。 subplot(211)。)。 xlabel(39。) axis([0 1 80 5])。gain39。subplot(211)。%設(shè)計模擬原型濾波器 [numi3,deni3]= lp2bp(B3i,A3i,W0,B)。)。 %歸一化，求低通原型的參數(shù)。 W3i = Wdp1i*Wdp2i/Wds1i。Wdp2i=tan(wdp2i/2)。wdp2i=2*pi*fdp2i/fs。fds2i=3200。%設(shè)置頻譜的間隔，分辨率 ，這里保證了x軸的點數(shù)必須和y軸點數(shù)一致fs=Fs。 IIR帶通濾波器設(shè)計及信號濾波[x1,Fs,bits]=wavread(39。plot([Fs/2:derta_Fs: Fs/2derta_Fs],abs(fftshift(fft(x1_highi))))。title(39。%對信號進(jìn)行高通濾波 ,Y = filter(B,A,X) ，輸入X為濾波前序列，Y為濾波結(jié)果序列，B/A 提供濾波器系數(shù)，B為分子， A為分母 sound(x1_highi,Fs,bits)。 title(39。)。s39。Ws2i=tan(Ws2i2/2)。 Rs2i=100。%設(shè)置頻譜的間隔，分辨率 ，這里保證了x軸的點數(shù)必須和y軸點數(shù)一致% IIR高通按鈕消息響應(yīng)函數(shù).fs=Fs。 IIR高通濾波器設(shè)計及信號濾波[x1,Fs,bits]=wavread(39。plot([Fs/2:derta_Fs: Fs/2derta_Fs],abs(fftshift(fft(x1_lowi))))。title(39。%對信號進(jìn)行低通濾波 ,Y = filter(B,A,X) ，輸入X為濾波前序列，Y為濾波結(jié)果序列，B/A 提供濾波器系數(shù)，B為分子， A為分母 sound(x1_lowi,Fs,bits)。 title(39。)。s39。 Ts=1/fs。wp1i=2*pi*fp1i/fs。 derta_Fs = Fs/length(x1)。)。 subplot(212)。plot(x1_daitong)。x1_daitong = filter(b3,1,x1)。freqz(b3,1,512)。 %求凱澤窗窗函數(shù)b3=fir1(M3,wc3/pi,39。M3=ceil(()/(*BF3))+1。wp3=2*pi*fp3/fs。 FIR帶通濾波器設(shè)計及信號濾波fs=Fs。plot([Fs/2:derta_Fs: Fs/2derta_Fs],abs(fftshift(fft(x1_high))))。title(39。%對信號進(jìn)行高通濾波sound(x1_high,Fs,bits)。%數(shù)字濾波器頻率響應(yīng) title(39。high39。%按凱澤窗計算濾波器階數(shù)N2=M2+1。 ws2=2*pi*fs2/fs。 As2=100。39。 title(39。信號經(jīng)過FIR低通濾波器(時域)39。 figure(3)。FIR低通濾波器的頻率響應(yīng)39。 freqz(b1,1,512)。beta1=*()。 %BF1=ws1wp1。 fs1=1200。)。圖2 原始信號頻譜圖39。)。%用傅里葉變換求頻域subplot(2,1,1)。)。五．語音信號的采集及頻譜分析用WINDOWS下的錄音機(jī)，用單聲道錄制一段音樂或聲音，時間在5S內(nèi)。由于這種頻率之間的非線性變換關(guān)系，就產(chǎn)生了新的問題。S平面整個jΩ軸單值地對應(yīng)于Z平面單位圓一周，即頻率軸是單值變換關(guān)系。也就是說，S平面的左半平面映射到Z平面的單位圓內(nèi)，S平面的右半平面映射到Z平面的單位圓外，S平面的虛軸映射到Z平面的單位圓上。當(dāng)Ω1由π/T經(jīng)過0變化到π/T時，Ω由∞經(jīng)過0變化到+∞，也即映射了整個jΩ軸。為了克服這一缺點，可以采用非線性頻率壓縮方法，將整個頻率軸上的頻率范圍壓縮到π/T～π/T之間，再用z=esT轉(zhuǎn)換到Z平面上。 例如， 要設(shè)計線性相位低通特性可選擇h(n)=h(N1n)一類，而不能選h(n)=h(N1n)一類。； 給定濾波器的性能指標(biāo)如下： (1)低通濾波器的性能指標(biāo)：fb=1000Hz,fc=1200Hz,As=100dB,Ap=1dB. (2)高通濾波器的性能指標(biāo)：fb=2700Hz, fc=3000Hz,As=100dB,Ap=1dB. (3)帶通濾波器的性能指標(biāo)： fb1=1200Hz, fb2=3000Hz, As=100dB,Ap=1dB. fc1=1000Hz, fc2=3200Hz,As=100dB, Ap=1dB. 采用窗函數(shù)法和雙線性變換法設(shè)計上面要求的3種濾波器，并畫出濾波器的頻率響應(yīng)； 然后用自己設(shè)計的濾波器對采集到的信號進(jìn)行濾波，畫出濾波后信號的時域波形及頻譜，并對濾波前后的信號進(jìn)行對比，分析信號的變化； ，分析濾波前后的語音變化； 四．濾波器設(shè)計原理及設(shè)計方案根據(jù)過渡帶寬及阻帶衰減要求，選擇窗函數(shù)的類型并估計窗口長度N（或階數(shù)M=N1），窗函數(shù)類型可根據(jù)最小阻帶衰減As獨立選擇，因為窗口長度N對最小阻帶衰減As沒有影響，在確
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1