課程設(shè)計(jì)--語音信號濾波去噪

1、<p><b>　　課程設(shè)計(jì)報(bào)告</b></p><p>　　題 目： 語音信號濾波去噪 </p><p><b>　　語音信號濾波去噪</b></p><p>　　1 課程設(shè)計(jì)的任務(wù)與要求</p><p>　　首先錄制一段自己的語音信號，并對錄制的信號進(jìn)行采樣

2、，修改其參數(shù)設(shè)置方便實(shí)驗(yàn)，再畫出采樣后語音信號的時(shí)域波形和頻譜圖。再通過程序編寫在程序調(diào)試中加入噪聲，再畫出加入噪聲后的時(shí)域波形和頻譜圖。接著自行給定濾波器的性能指標(biāo)，采用IIR濾波器，并畫出濾波器的頻率響應(yīng)。用自己設(shè)計(jì)的濾波器對采集的信號進(jìn)行濾波，畫出濾波后信號的時(shí)域波形和頻譜，并對濾波前后的信號進(jìn)行對比，分析信號的變化并回放語音信號；通過利用各種不同的開發(fā)工具實(shí)現(xiàn)語音信號的濾波去噪，掌握數(shù)字信號的分析方法和處理方法.</p&g

3、t;<p><b>　　2 設(shè)計(jì)原理</b></p><p><b>　　2.1設(shè)計(jì)流圖</b></p><p>　　MATLAB 語言是一種數(shù)據(jù)分析和處理功能十分強(qiáng)大的計(jì)算機(jī)應(yīng)用軟件 ,它可以將聲音文件變換為離散的數(shù)據(jù)文件 , 然后利用其強(qiáng)大的矩陣運(yùn)算能力處理數(shù)據(jù) ,如數(shù)字濾波 、傅里葉變換 、時(shí)域和頻域分析 、聲音回放以及各種圖

4、的呈現(xiàn)等 , 信號處理是MATLAB 重要應(yīng)用的領(lǐng)域之一 。</p><p>　　本課程設(shè)計(jì)介紹了基于Matlab的對語音信號采集、處理及濾波器的設(shè)計(jì)，并使之實(shí)現(xiàn)的過程。理解與掌握課程中的基本概念、基本原理、基本分析方法，用Matlab進(jìn)行數(shù)字語音信號處理，并闡述了課程設(shè)計(jì)的具體方法、步驟和內(nèi)容。綜合運(yùn)用本課程的理論知識進(jìn)行頻譜分析以及濾波器設(shè)計(jì)，通過理論推導(dǎo)得出相應(yīng)結(jié)論，并利用MATLAB作為工具進(jìn)行實(shí)現(xiàn)，從而

5、復(fù)習(xí)鞏固課堂所學(xué)的理論知識，提高對所學(xué)知識的綜合應(yīng)用能力，并從實(shí)踐上初步實(shí)現(xiàn)對數(shù)字信號的處理。</p><p>　　采樣器的作用是把連續(xù)信號變?yōu)槊}沖或數(shù)字序列。圖中示出了一個(gè)連續(xù)信號f(t)經(jīng)采樣器采樣后變?yōu)殡x散信號的過程【1】。</p><p>　　圖1 連續(xù)信號f(t)經(jīng)采樣器采樣后變?yōu)殡x散信號</p><p>　　圖中f(t)為被采樣的連續(xù)信號，S（t）為周期

6、性窄脈沖信號，fs(t)為采樣后的離散信號，它用下式來表征：fs(t)=f(t)s(t) </p><p>　　采樣信號的頻率特性為：</p><p>　　如果｜F*(jω)｜中各個(gè)波形不重復(fù)搭接，相互間有一定的距離(頻率)即若</p><p>　　即采樣定理可敘述如下：如果采樣周期滿足下列條件，即：</p><p>　　式中ωmax為連續(xù)信

7、號f(t)的最高次諧波的角頻率。則采樣信號f*(t)</p><p>　　就可以無失真地再恢復(fù)為原連續(xù)信號f(t)。需要指出的是，采樣定理只是在理論上給出了信號準(zhǔn)確復(fù)現(xiàn)的條件。</p><p>　　我們可以利用windows自帶的錄音機(jī)錄制語音文件，進(jìn)行數(shù)字信號的采集。（開始—程序—附件—娛樂—錄音機(jī)，文件—屬性—立即轉(zhuǎn)換—8000KHz，8位，單聲道）或其他軟件，將話筒輸入計(jì)算機(jī)的語音

8、輸入插口上,啟動錄音機(jī)，錄制一段自己的話音。然后在MATLAB軟件平臺下，利用函數(shù)wavread對語音信號進(jìn)行采樣，記住采樣頻率和采樣點(diǎn)數(shù)。通過wavread函數(shù)的使用進(jìn)一步理解采樣頻率、采樣位數(shù)等概念。圖2是基于PC機(jī)的語音信號采集過程:</p><p><b>　　聲音</b></p><p><b>　　圖2 語音采集過程</b></

9、p><p>　　3 濾波去噪的仿真與實(shí)現(xiàn)</p><p><b>　　3.1錄制語音信號</b></p><p>　　此次錄制語音信號，需用到電腦附件中的錄音機(jī)功能，在機(jī)房的電腦系統(tǒng)為Windows XP系統(tǒng)，在筆記本的win7系統(tǒng)中，錄制的語音信號格式為.wma格式文件，需用格式轉(zhuǎn)換器進(jìn)行轉(zhuǎn)換。在機(jī)房電腦中，點(diǎn)擊電腦桌面上的“開始”菜單，再選擇“

10、程序”中的“附件”，在“附件”的菜單欄中選擇“娛樂”，最后點(diǎn)擊“錄音機(jī)”，就可以得到如圖3.2的圖。然后點(diǎn)擊開始錄制一段“數(shù)字通信原理課程設(shè)計(jì)”的語音信號，時(shí)間大約在2~3秒之間。錄制好語音信號之后點(diǎn)擊另存為，在另存為的對話框中右下角有“更改”選項(xiàng)，進(jìn)行更改音頻文件的參數(shù)，其參數(shù)設(shè)置的值如圖3.3所示。最后將這段語音信號保存在MATLAB下的Work文檔中以便以后調(diào)用。并命名為“ywt【4】”。</p><p>

11、　　圖3 XP系統(tǒng)下的錄音機(jī)</p><p>　　圖4 語音信號參數(shù)設(shè)置</p><p>　　3.2語音信號的讀取</p><p>　　在matlab軟件平臺下，利用wavread函數(shù)對語音信號進(jìn)行采樣，記住采樣頻率和采樣點(diǎn)數(shù)。通過使用wavread函數(shù)，理解采樣頻率、采樣位數(shù)等概念。</p><p>　　wavread函數(shù)調(diào)用格式<

12、;/p><p>　　y=wavread(file)，讀取file所規(guī)定的wav文件，返回采樣值放在向量y中。</p><p>　　[y,fs,nbits]=wavread(file)，采樣值放在向量y中，fs表示采樣頻率（hz），nbits表示采樣位數(shù)。</p><p>　　y=wavread（file，N），讀取前N點(diǎn)的采樣值放在向量y中。</p>&l

13、t;p>　　y=wavread（file，[N1,N2]），讀取從N1到N2點(diǎn)的采樣值放在向量y中。</p><p>　　對語音信號“ywt.wav”進(jìn)行采樣其程序如下：</p><p>　　[y,fs,nbits]=wavread('ywt.wav')</p><p>　　語音信號是一種非平穩(wěn)的時(shí)變信號，它攜帶著各種信息。在語音編碼、語音合成

14、、語音識別和語音增強(qiáng)等語音處理中無一例外需要提取語音中包含的各種信息。語音信號分析的目的就在與方便有效的提取并表示語音信號所攜帶的信息。語音信號分析可以分為時(shí)域和變換域等處理方法，其中時(shí)域分析是最簡單的方法，直接對語音信號的時(shí)域波形進(jìn)行分析，提取的特征參數(shù)主要有語音的短時(shí)能量，短時(shí)平均過零率，短時(shí)自相關(guān)函數(shù)等。</p><p><b>　　圖5 原始信號波形</b></p>&

15、lt;p>　　3.3 噪聲信號的構(gòu)建</p><p>　　在MATLAB軟件平臺下，給原始的語音信號疊加上噪聲，噪聲類型分為如下幾種：白噪聲、單頻噪聲（正弦干擾）、多頻噪聲（多正弦干擾）、其他干擾，可設(shè)置為低頻、高頻、帶限噪聲、Chirp干擾、沖激干擾。</p><p>　　繪制出疊加噪聲之后的語音信號時(shí)域圖形及頻域圖形，在視覺上與原始信號圖形對比，也可以通過Windows播放軟件或

16、者是利用MATLAB中的sound()函數(shù)從聽覺上進(jìn)行對比，分析并體會含噪語音信號頻譜和時(shí)域波形的改變。</p><p>　　利用MATLAB中的隨機(jī)函數(shù)(randn)產(chǎn)生噪聲加入到語音信號中，模仿語音信號被污染，并對其進(jìn)行頻譜分析。程序如下：</p><p>　　N=length(y); %求出語音信號的長度</p><p>　　noise

17、=rand(N,2)/20; %噪聲信號的函數(shù)</p><p>　　z=fft(noise); %快速傅里葉變換</p><p>　　subplot(2 ,1 ,1),plot(noise);title('噪聲信號波形');</p><p>　　subplot(2 ,1 ,2),plot(abs(z));title(&

18、#39;噪聲信號頻譜');</p><p>　　axis([0,250000,0,100]);</p><p><b>　　圖6 噪聲信號</b></p><p>　　3.4 濾波器的設(shè)計(jì)原理</p><p>　　利用模擬濾波器來設(shè)計(jì)數(shù)字濾波器，也就是使數(shù)字濾波器能模仿模擬濾波器的特性，這種模仿可以從不同的角度出

19、發(fā)。脈沖響應(yīng)不變法是從濾波器的脈沖響應(yīng)出發(fā)，使數(shù)字濾波器的單位脈沖響應(yīng)序列h(n)模仿模擬濾波器的沖激響應(yīng)ha(t)，即將ha(t)進(jìn)行等間隔采樣，使h(n)正好等于ha(t)的采樣值，滿足</p><p>　　h(n)=ha(nT)</p><p>　　式中,T是采樣周期。</p><p>　　如果令Ha(s)是ha(t)的拉普拉斯變換，H(z)為h(n)的Z變換

20、，利用采樣序列的Z變換與模擬信號的拉普拉斯變換的關(guān)系得</p><p><b>　　(1-1)</b></p><p>　　則可看出，脈沖響應(yīng)不變法將模擬濾波器的S平面變換成數(shù)字濾波器的Z平面，這個(gè)從s到z的變換z=esT是從S平面變換到Z平面的標(biāo)準(zhǔn)變換關(guān)系式【2】。</p><p>　　圖7 脈沖響應(yīng)不變法的映射關(guān)系</p>

21、<p>　　由（1-1）式，數(shù)字濾波器的頻率響應(yīng)和模擬濾波器的頻率響應(yīng)間的關(guān)系為</p><p><b>　　(1-2)</b></p><p>　　這就是說，數(shù)字濾波器的頻率響應(yīng)是模擬濾波器頻率響應(yīng)的周期延拓。正如采樣定理所討論的，只有當(dāng)模擬濾波器的頻率響應(yīng)是限帶的，且?guī)抻谡郫B頻率以內(nèi)時(shí)，即</p><p><b>　

22、　(1-3)</b></p><p>　　才能使數(shù)字濾波器的頻率響應(yīng)在折疊頻率以內(nèi)重現(xiàn)模擬濾波器的頻率響應(yīng)，而不產(chǎn)生混疊失真，即</p><p>　　|ω|<π (1-4)</p><p>　　但是，任何一個(gè)實(shí)際的模擬濾波器頻率響應(yīng)都不是嚴(yán)格限帶的，變換后就會產(chǎn)生周期延拓分量的頻譜交疊，即產(chǎn)生頻率響應(yīng)的混疊失真，如圖7-

23、4所示。這時(shí)數(shù)字濾波器的頻響就不同于原模擬濾波器的頻響，而帶有一定的失真。當(dāng)模擬濾波器的頻率響應(yīng)在折疊頻率以上處衰減越大、越快時(shí)，變換后頻率響應(yīng)混疊失真就越小。這時(shí)，采用脈沖響應(yīng)不變法設(shè)計(jì)的數(shù)字濾波器才能得到良好的效果。</p><p>　　圖8 脈沖響應(yīng)不變法中的頻響混疊現(xiàn)象</p><p>　　對某一模擬濾波器的單位沖激響應(yīng)ha(t)進(jìn)行采樣，采樣頻率為fs，若使fs增加，即令采樣時(shí)

24、間間隔（T=1/fs）減小，則系統(tǒng)頻率響應(yīng)各周期延拓分量之間相距更遠(yuǎn)，因而可減小頻率響應(yīng)的混疊效應(yīng)[3]。</p><p>　　從以上討論可以看出，脈沖響應(yīng)不變法使得數(shù)字濾波器的單位脈沖響應(yīng)完全模仿模擬濾波器的單位沖激響應(yīng)，也就是時(shí)域逼近良好，而且模擬頻率Ω和數(shù)字頻率ω之間呈線性關(guān)系ω=ΩT。因而，一個(gè)線性相位的模擬濾波器（例如貝塞爾濾波器）通過脈沖響應(yīng)不變法得到的仍然是一個(gè)線性相位的數(shù)字濾波器。</p&g

25、t;<p>　　脈沖響應(yīng)不變法的最大缺點(diǎn)是有頻率響應(yīng)的混疊效應(yīng)。所以，脈沖響應(yīng)不變法只適用于限帶的模擬濾波器(例如，衰減特性很好的低通或帶通濾波器)，而且高頻衰減越快，混疊效應(yīng)越小。至于高通和帶阻濾波器，由于它們在高頻部分不衰減，因此將完全混淆在低頻響應(yīng)中。如果要對高通和帶阻濾波器采用脈沖響應(yīng)不變法，就必須先對高通和帶阻濾波器加一保護(hù)濾波器，濾掉高于折疊頻率以上的頻率，然后再使用脈沖響應(yīng)不變法轉(zhuǎn)換為數(shù)字濾波器。當(dāng)然這樣會進(jìn)

26、一步增加設(shè)計(jì)復(fù)雜性和濾波器的階數(shù)。</p><p>　　脈沖響應(yīng)不變法的主要缺點(diǎn)是產(chǎn)生頻率響應(yīng)的混疊失真。這是因?yàn)閺腟平面到Ｚ平面是多值的映射關(guān)系所造成的。為了克服這一缺點(diǎn)，可以采用非線性頻率壓縮方法，將整個(gè)頻率軸上的頻率范圍壓縮到-π/T～π/T之間，再用z=esT轉(zhuǎn)換到Z平面上。也就是說，第一步先將整個(gè)S平面壓縮映射到S1平面的-π/T～π/T一條橫帶里；第二步再通過標(biāo)準(zhǔn)變換關(guān)系z=es1T將此橫帶變換到整個(gè)

27、Z平面上去。這樣就使S平面與Z平面建立了一一對應(yīng)的單值關(guān)系，消除了多值變換性，也就消除了頻譜混疊現(xiàn)象</p><p>　　IIR數(shù)字濾波器是一種離散時(shí)間系統(tǒng)，其系統(tǒng)函數(shù)為</p><p><b>　　（1-5）</b></p><p>　　假設(shè)M≤N，當(dāng)M＞N時(shí),系統(tǒng)函數(shù)可以看作一個(gè)IIR的子系統(tǒng)和一個(gè)(M-N)的FIR子系統(tǒng)的級聯(lián)。IIR數(shù)字

28、濾波器的設(shè)計(jì)實(shí)際上是求解濾波器的系數(shù)和 ，它是數(shù)學(xué)上的一種逼近問題，即在規(guī)定意義上（通常采用最小均方誤差準(zhǔn)則）去逼近系統(tǒng)的特性。如果在S平面上去逼近，就得到模擬濾波器；如果在z平面上去逼近，就得到數(shù)字濾波器【2】。</p><p><b>　　設(shè)計(jì)流程如圖所示：</b></p><p>　　頻域指標(biāo) 雙線性變換法</p

29、><p><b>　　脈沖響應(yīng)不變法</b></p><p>　　圖9 IIR濾波器的間接法設(shè)計(jì)</p><p><b>　　IIR低通濾波器：</b></p><p><b>　　Ft=8000;</b></p><p><b>　　Fp=100

30、0;</b></p><p><b>　　Fs=1200;</b></p><p>　　wp=2*pi*Fp/Ft;</p><p>　　ws=2*pi*Fs/Ft;</p><p>　　fp=2*Ft*tan(wp/2);</p><p>　　fs=2*Ft*tan(ws/2);&l

31、t;/p><p>　　[N,wc]=ellipord(wp,ws,1,50,'s');</p><p>　　[b,a]=ellip(N,wc,'s');</p><p>　　[bz,az]=bilinear(b,a,0.5);</p><p>　　[h,w]=freqz(bz,az);</p><

32、;p>　　plot(w*8000*0.5/pi,abs(h));</p><p>　　Legend('用橢圓設(shè)計(jì)');</p><p>　　title('IIR低通濾波器');</p><p><b>　　grid;</b></p><p>　　圖10 IIR低通濾波器</p

33、><p>　　3.5濾波器對信號濾波</p><p>　　用設(shè)計(jì)好的數(shù)字濾波器對含噪語音信號進(jìn)行濾波,在Matlab 中IIR 濾波器利用函數(shù)fftfilt 對信號進(jìn)行濾波, IIR濾波器利用函數(shù)filter 對信號進(jìn)行濾波，在一個(gè)窗口同時(shí)畫出濾波前后的波形及頻譜.</p><p>　　matlab程序如下： </p><p>　　y=wavre

34、ad('ywt.wav');</p><p>　　N=length(y);</p><p>　　noise=rand(N,2)/20;</p><p>　　s=y+noise; %噪聲信號的疊加</p><p>　　Subplot(2,1,1);</p><p>　　plot(s);titl

35、e ('濾波前的時(shí)域波形');</p><p><b>　　S=fft(s);</b></p><p>　　Subplot(2,1,2);</p><p>　　plot(abs(S));title ('濾波前的頻域波形');</p><p>　　axis([0,45000,0,200]);

36、</p><p>　　圖11 濾波前后的信號波形</p><p>　　從圖看出，加噪聲后的語音信號經(jīng)過FIR濾波器濾波后，時(shí)域和頻域圖與原始信號幾乎一樣，說明噪聲幾乎都被濾掉了，也說明FIR濾波器設(shè)計(jì)很理想，能滿足所需要求。</p><p><b>　　3.6 結(jié)果分析</b></p><p>　　要確定本課程設(shè)計(jì)是否

37、成功就得看原始信號的頻域圖和時(shí)域圖與經(jīng)過濾波器后的語音</p><p>　　信號的頻域圖和時(shí)域圖是否一樣，若一樣則表示該設(shè)計(jì)是成功的，否則是不成功的。</p><p>　　在第一個(gè)圖中：第一幅圖和第二幅圖是原始語音信號的時(shí)域圖和頻域圖,第二個(gè)圖中：第一幅圖和第二幅圖是加入隨機(jī)噪聲時(shí)域圖和頻域圖。在第三個(gè)圖中：為濾波前后的時(shí)域及頻域，從圖中可以看出，第一圖和第三圖相比因?yàn)榧尤朐肼暤木壒仕缘?/p>

38、三圖y軸的幅度要比第一圖要大，但其形狀還是基本沒有改變。</p><p><b>　　4 心得與體會</b></p><p>　　通過這次的課程設(shè)計(jì)我們可以學(xué)的到很多的東西，不僅可以鞏固以前所學(xué)過的知識，還可以學(xué)到很多在書本上所沒有學(xué)到過的知識。進(jìn)一步加深了對數(shù)字信號處理的了解，讓我對它有了更加濃厚的興趣。因?yàn)橐郧皩W(xué)過信號與系統(tǒng)和數(shù)字信號處理兩門課，但這只是理論知識，

39、通過實(shí)驗(yàn)我們才能真正理解其意義。語音信號處理在以前的學(xué)習(xí)中并沒有怎么接觸，但是通過網(wǎng)上查找資料，我很快的理解了其基本原理。</p><p>　　在設(shè)計(jì)的過程中我還是遇到不少的問題的，比如剛開始，要處理的是WAV格式的語音信號，可是我發(fā)現(xiàn)并不是所有的WAV文件都可以用WAVREAD函數(shù)讀取的，這才了解到WAV文件之間也是有區(qū)別的，并且如果格式不對的話，我們要把格式轉(zhuǎn)換成WAV格式。剛開始由于對濾波器的濾波原理并不是

40、很了解，于是我又翻出學(xué)過的數(shù)字信號處理課本，認(rèn)真研究起各種濾波器了，這才使我明白了大多數(shù)濾波器是如何工作地，不再單單只是懂理論，理論與實(shí)際相結(jié)合是很重要的，只有理論知識是遠(yuǎn)遠(yuǎn)不夠的，只有把所學(xué)的理論知識與實(shí)踐相結(jié)合起來，從理論中得出結(jié)論。</p><p>　　語音信號是一種非平穩(wěn)的時(shí)變信號，它攜帶著各種信息。語音信號分析的目的就在與方便有效的提取并表示語音信號所攜帶的信息。語音信號分析可以分為時(shí)域和變換域等處理方

41、法，其中時(shí)域分析是最簡單的方法，直接對語音信號的時(shí)域波形進(jìn)行分析，提取的特征參數(shù)。語音信號經(jīng)過FFT變換之后，就可以得到信號的頻域圖形，另外，頻域圖形使信號的某些特性變得更明顯，信號的傅立葉表示在信號的分析與處理中起著重要的作用。因?yàn)閷τ诰€性系統(tǒng)來說，可以很方便地確定其對正弦或復(fù)指數(shù)和的響應(yīng)，所以傅立葉分析方法能完善地解決許多信號分析和處理問題。</p><p>　　應(yīng)用MATLAB 根據(jù)給定的技術(shù)指標(biāo)進(jìn)行FIR

42、 數(shù)字濾波器設(shè)計(jì)時(shí)，可根據(jù)設(shè)計(jì)要求隨時(shí)改變參數(shù)，以使濾波器達(dá)到最優(yōu)化。 使用凱瑟窗設(shè)計(jì)低通、高通濾波器，并對疊加噪聲后的語音信號進(jìn)行濾波。比較濾波前后的時(shí)域波形和頻譜并進(jìn)行分析。在頻域波形中，我們可以明顯的看到設(shè)計(jì)的濾波器對語音信號進(jìn)行了濾波處理，將噪聲進(jìn)行了濾除。此次濾波基本達(dá)到了要求，完成了設(shè)計(jì)指標(biāo)。</p><p>　　總之，加噪聲的后的語音信號經(jīng)過IIR濾波器的濾波處理，時(shí)域和頻域圖與原始信號都

43、幾乎一樣，這說明噪聲幾乎全被濾掉了，同時(shí)也說明IIR濾波器設(shè)計(jì)還算理想，能滿足所需要求。也可看出濾波之后仍有一小段噪聲沒有濾除，后面經(jīng)過調(diào)試濾波器的參數(shù)發(fā)現(xiàn)當(dāng)fp、fs參數(shù)的數(shù)值相差增大時(shí)，剩余噪聲便會減少，于是通過進(jìn)一步調(diào)試得到了較好的效果。 </p><p>　　低通濾波后，聲音稍微有些發(fā)悶、低沉，原因是高頻分量被低通濾波器衰減。但是很接近原來的聲音。高通濾波器濾波以后，此時(shí)只有少許雜音，原因是低頻

44、分量被高通濾波器衰減，而人聲部分正好是低頻部分，所以只剩下雜音，或者發(fā)出高頻雜音但人的耳朵聽不到。</p><p>　　本課程設(shè)計(jì)的主要目的是在MATLAB環(huán)境下，使用窗函數(shù)設(shè)計(jì)法設(shè)計(jì)一個(gè)FIR濾波器，并對疊加了噪聲后的語音信號進(jìn)行濾波去噪。此次濾波基本達(dá)到了要求，完成了設(shè)計(jì)指標(biāo)。 </p><p>　　總的來說，通過這次的課程設(shè)計(jì)我對語音信號有了全面的認(rèn)識，對數(shù)字信號處理的知識

45、又有了深刻的理解，讓我感受到只有在充分理解課本知識的前提下，才能更好的應(yīng)用這個(gè)工具；并且熟練的應(yīng)用MATLAB也可以很好的加深我對課程的理解，方便我的思維。這次設(shè)計(jì)使我了解了MATLAB的使用方法，學(xué)會分析濾波器的優(yōu)劣和性能，提高了分析和動手實(shí)踐能力。</p><p><b>　　參考文獻(xiàn)</b></p><p>　　[1]陳懷琛.MATLAB及在電子信息課程中的應(yīng)用

46、[M].北京：電子工業(yè)出版社 2006年3月。</p><p>　　[2]張文.基于MATLAB的語音信號的濾波與實(shí)現(xiàn)[J].山西電子技   2002年2月。</p><p>　　[3]蔡蓮紅，黃德智．現(xiàn)代語音技術(shù)基礎(chǔ)與應(yīng)用    [M] 清華大學(xué)出版社   2003年11月。  &l