【文章內(nèi)容簡(jiǎn)介】 所有分量 (fs 為采樣頻率 )，以防止混疊干擾。 ② 抑制 50Hz 的電源工頻干擾。這樣，預(yù)濾波器必須是一個(gè)帶通濾波器，設(shè)其上、下截止顏率分別是 fH 和 fL，則對(duì)于絕人多數(shù)語(yǔ)音編譯碼器， fH=3400Hz、 fL＝ 60~100Hz、采樣率為 fs＝ 8kHz；而對(duì)丁語(yǔ)音識(shí)別而言，當(dāng)用于電話用戶時(shí)，指標(biāo)與語(yǔ)音編譯碼器相同。當(dāng)使用要求較高或很高的場(chǎng)合時(shí) fH＝ 4500Hz 或 8000Hz、 fL＝ 60Hz、 fs＝ 10kHz 或 20kHz。 基于 Matlab 的聲信號(hào)采集與譜分析設(shè)計(jì) 作者： 尹青山 第 8 頁(yè) 共 22 頁(yè) 為了將原始模擬語(yǔ)音信號(hào)變?yōu)閿?shù)字信號(hào)，必須經(jīng)過(guò)采樣和量化兩個(gè)步驟，從而得到時(shí)間和幅度上均為離散的數(shù)字語(yǔ)音信號(hào)。 采樣也稱抽樣，是信號(hào)在時(shí)間上的離散化，即按照一定時(shí)間間隔△ t 在模擬信號(hào) x(t)上逐點(diǎn)采取其瞬時(shí)值。 采樣時(shí)必須要注意滿足 奈奎斯特定理 ，即采樣頻率 fs 必須 以高于受測(cè)信號(hào)的最高頻率兩倍以上的速度進(jìn)行取樣，才能正確地重建波 它是通過(guò)采樣脈沖和模擬信號(hào)相乘來(lái)實(shí)現(xiàn)的 。下圖（ 圖 3 3）是 一段語(yǔ)音信號(hào)在采樣頻率 情況下的頻譜圖。 圖 3 3 原始信號(hào) 由圖可知，這段語(yǔ)音信號(hào)的頻率主要集中在 1KHz 左右，當(dāng)采樣頻率為 時(shí)，由于采樣頻率比較大，所以 采樣點(diǎn)數(shù)就越密，所得離散信號(hào)就越逼近于原信號(hào) ，頻譜也沒(méi)有發(fā)生混疊 （ 圖 3 4） 。 基于 Matlab 的聲信號(hào)采集與譜分析設(shè)計(jì) 作者： 尹青山 第 9 頁(yè) 共 22 頁(yè) 圖 3 4抽樣后的信號(hào) 對(duì)上述信號(hào)進(jìn)行 1/80 采樣頻率抽取，即采樣頻率變?yōu)閷⒔?500Hz 時(shí)，由于采樣頻率比較小，所以采樣點(diǎn)數(shù)就稀疏 ，所得離散信號(hào)就越 偏離 于原信號(hào) ，頻譜也發(fā)生了混疊 [7]。 在采樣的過(guò)程中應(yīng)注意 采樣間隔的選擇和信號(hào)混淆：對(duì)模擬信號(hào)采樣首先要確定采樣間隔。如何合理選擇△ t 涉及到許多需要考慮的技術(shù)因素。一般而言，采樣頻率越高，采樣點(diǎn)數(shù)就越密，所得離散信號(hào)就越逼近于原信號(hào)。但過(guò)高的采樣頻率并不可取，對(duì)固定長(zhǎng)度（ T）的信號(hào)，采集到過(guò)大的數(shù)據(jù)量（ N=T/△ t），給計(jì)算機(jī)增加不必要的計(jì)算工作量和存儲(chǔ)空間；若數(shù)據(jù)量（ N）限定，則采 樣時(shí)間過(guò)短，會(huì)導(dǎo)致一些數(shù)據(jù)信息被排斥在外。采樣頻率過(guò)低，采樣點(diǎn)間隔過(guò)遠(yuǎn)，則離散信號(hào)不足以反映原有信號(hào)波形特征，無(wú)法使信號(hào)復(fù)原，造成信號(hào)混淆。 根據(jù)采樣定理，當(dāng)采樣頻率大于信號(hào)的兩倍帶寬時(shí)，采樣過(guò)程不會(huì)丟失信息，利用理想濾波器可從采樣信號(hào)中不失真地重構(gòu)原始信號(hào)波形。量化 是對(duì)幅值進(jìn)行離散化，即將振動(dòng)幅值用二進(jìn)制量化電平來(lái)表示。量化電平按級(jí)數(shù)變化，實(shí)際的振動(dòng)值是連續(xù)的物理量。具體振值用舍入法歸到靠近的量化電平上。 語(yǔ)音信號(hào)經(jīng)過(guò)預(yù)濾波和采樣后，由 A／ D 變換器變換為二址制數(shù)字碼。這種防混疊濾波通常與模數(shù)轉(zhuǎn)換器做在一個(gè) 集成塊內(nèi)，因此目前來(lái)說(shuō)，語(yǔ)音信號(hào)基于 Matlab 的聲信號(hào)采集與譜分析設(shè)計(jì) 作者： 尹青山 第 10 頁(yè) 共 22 頁(yè) 的數(shù)字化的質(zhì)量還是有保證的。市面上購(gòu)買到的普通聲卡在這方面做的都很好，語(yǔ)音聲波通過(guò)話筒輸入到聲卡后直接獲得的是經(jīng)過(guò)防混疊濾波、 A/D 變換、量化處理的離散的數(shù)字信號(hào) [8]。 在實(shí)際工作中，可以利用 windows 自帶的錄音機(jī)錄制語(yǔ)音文件， 圖 3 5 是基于 PC 機(jī)的語(yǔ)音信號(hào)采集過(guò)程，聲卡可以完成語(yǔ)音波形的 A/D 轉(zhuǎn)換，獲得WAVE 文件，為后續(xù)的處理儲(chǔ)備原材料。調(diào)節(jié)錄音機(jī)保存界面的“更改”選項(xiàng)，可以存儲(chǔ) 各種格式的 WAVE 文件。 圖 3 5 基于 PC 機(jī)的語(yǔ)音信號(hào)采集過(guò)程 在 MATLAB 中， [y,fs,bits]=wavread(39。Blip39。,[N1 N2])。用于讀取語(yǔ)音，采樣值放在向量 y 中， fs 表示采樣頻率 (Hz)， bits 表示采樣位數(shù)。 [N1 N2]表示讀取從 N1 點(diǎn)到 N2 點(diǎn)的值（若只有一個(gè) N 的點(diǎn)則表示讀取前 N 點(diǎn)的采樣值）。 sound(x,fs,bits)。 用于對(duì)聲音的回放。向量 y 則就代 表了一個(gè)信號(hào)（也即一個(gè)復(fù)雜的“函數(shù)表達(dá)式”）也就是說(shuō)可以像處理一個(gè)信號(hào)表達(dá)式一樣處理這個(gè)聲音信號(hào)。 Windows 自帶的錄音機(jī) 聲音 麥克風(fēng) 聲卡 濾波 采樣 A/D 轉(zhuǎn)換 Wav 基于 Matlab 的聲信號(hào)采集與譜分析設(shè)計(jì) 作者： 尹青山 第 11 頁(yè) 共 22 頁(yè) 第 4章 語(yǔ)音信號(hào)的分析 語(yǔ)音信號(hào)的時(shí)域分析 語(yǔ)音信號(hào)的時(shí)域分析就是分析和提取語(yǔ)音信號(hào)的時(shí)域參數(shù)。進(jìn)行語(yǔ)音分析時(shí)，最先接觸到并且也是最直觀的是它的時(shí)域波形。語(yǔ)音信號(hào)本身就是時(shí)域信號(hào)，因而時(shí)域分析是最早使用，也是應(yīng)用最廣泛的一種分析方法，這種方法直接利用語(yǔ)音信號(hào)的時(shí)域波形。時(shí)域分析通常用于最基本的參數(shù)分析及應(yīng)用，如語(yǔ)音的分割、預(yù)處理、大分類等。這種分析方法的特點(diǎn)是： ① 表示 語(yǔ)音信號(hào)比較直觀、物理意義明確。 ② 實(shí)現(xiàn)起來(lái)比較簡(jiǎn)單、運(yùn)算且少 。 ③ 可以得到語(yǔ)音的一些重要的參數(shù)。④只使用示波器等通用設(shè)備，使用較為簡(jiǎn)單等 [9]。 語(yǔ)音信號(hào)的時(shí)域參數(shù)有短時(shí)能量、短時(shí)過(guò)零率、短時(shí)白相關(guān)函數(shù)和短時(shí)平均幅度差函數(shù)等，這是語(yǔ)音信號(hào)的一組最基本的短時(shí)參數(shù)，在各種語(yǔ)音信號(hào)數(shù)字處理技術(shù)中都要應(yīng)用 [10]。在計(jì) 算這些參數(shù)時(shí)使用的一 般是方 窗 或漢明 窗 。 短時(shí)能量及短時(shí)平均幅度分析 設(shè)語(yǔ)音波形時(shí)域信號(hào)為 x(l)、加密分幀處理后得到的第 n 幀語(yǔ)音信號(hào)為 Xn(m)，則 Xn(m)滿足下式： ( ) ( ) ( ) ( )nx m w m x n m m?? 01mN? ? ? (41) 1 0 ~ ( 1 )( ) { mNmwm ???? ， 0， 其 他 值 (42) 其中， n＝ 0， 1T， 2T，?，并且 N 為幀長(zhǎng)， T 為幀移長(zhǎng)度。 設(shè) 第 n 幀語(yǔ)音信號(hào) Xn(m)的短時(shí)能量用 En 表示，則其計(jì) 算公式如下： 1 20 ()NnnmE x m???? (43) En 是一個(gè)度量語(yǔ)音信號(hào)幅度值變化的函數(shù)，但它有一 個(gè)缺陷，即它對(duì)高電平非常敏感 (因?yàn)樗?計(jì)算時(shí)用的是信號(hào)的平方 )。為此可采用另一個(gè)度量語(yǔ)音信號(hào)幅度值變化的函數(shù)．即短時(shí)平均幅度函數(shù) Mn， 它定義為： 10 ()NMnnm xm???? (44) 基于 Matlab 的聲信號(hào)采集與譜分析設(shè)計(jì) 作者： 尹青山 第 12 頁(yè) 共 22 頁(yè) Mn 也是一帕語(yǔ)音信號(hào)能量 大小的表征，它與 En 的區(qū)別在于計(jì)算時(shí)小取樣值和大取樣值不會(huì)因取平方而造成較大差異，在某些應(yīng)用領(lǐng)域呻會(huì)帶來(lái)一些好處。 短時(shí)能量和短時(shí)平均幅度 函數(shù) 的主要用途有：①可以區(qū)分濁音段與清音段，因?yàn)闈嵋魰r(shí) En 值比清音時(shí)大的多。 ② 可以用來(lái)區(qū)分聲母與韻母的分界，無(wú)聲與有盧的分界，連字 (指字之間無(wú)間隙 )的分界等。③作為一種超音段信息，用 于 語(yǔ)音識(shí)別 中 短時(shí)過(guò)零率分析 短時(shí)過(guò) 零 率表示一幀語(yǔ)音中語(yǔ)音信號(hào)波形穿過(guò)橫軸 （零 電平 ） 的次數(shù)。過(guò)零 分析是語(yǔ)音時(shí)域分析中最簡(jiǎn)單的一種 。 對(duì)于連續(xù)語(yǔ)音信號(hào)，過(guò)零即意味著 時(shí)域波形通過(guò)時(shí)間軸；而對(duì)于離散信號(hào)，如果相鄰的取樣值改變符號(hào)則稱為過(guò)零。過(guò)零率就是樣本改變符號(hào)的次數(shù) [11]。定義語(yǔ)音信號(hào) Xn（ m） 的短時(shí)過(guò)零 率Zn 為： ? ? ? ?101 s g n ( ) s g n ( 1 )2 Nn n nmZ x m x m??? ? ?? （ 45） 式中， sgn[ ]是符號(hào)函 數(shù)，即： ? ? ? 1 ,( 0 )1 ,( 0 )sg n x xx ???? （ 46） 在實(shí)際中求過(guò)零率參數(shù)時(shí)，需要十分注 意的一 個(gè)問(wèn)題是如果輸入信號(hào)中包含有 50Hz 的工頻干擾或者 A/D 變換器的工作點(diǎn)有偏移 (這等效于輸入信號(hào)有直流偏移 )，往往會(huì)使計(jì)算的過(guò)零率參數(shù)很不 準(zhǔn) 確 。 為了解決前一個(gè)問(wèn)題， A/D 變換器前的防混疊帶通濾波器的低端截頻應(yīng)高于 50Hz，以有效地抑制電源干擾。對(duì)于后一個(gè)問(wèn)題除了可以采用低直流漂移器件外，也 可以在軟件上加以解決，這就是算出每一幀的直流分量并予以濾除 。 對(duì)語(yǔ)音信號(hào)進(jìn)行分析，發(fā)現(xiàn)發(fā)濁音時(shí)，盡管聲道有若干個(gè)共振峰，但由于聲門波引起譜的 高頻跌落，所以其話音能量約集中在 3kHz 以下。而發(fā)清音時(shí)，多數(shù)能量出現(xiàn)在較高頻 率上。高頻就意味著高的平均過(guò)零率，低頻意味著低的平均過(guò)零率，所以可以認(rèn)為濁音時(shí)具有較低的過(guò)零率，而清音時(shí)具有較高的過(guò)零率。當(dāng)然，這種高低僅是相對(duì)而言，并沒(méi)方精確的數(shù)值關(guān)系。 基于 Matlab 的聲信號(hào)采集與譜分析設(shè)計(jì) 作者： 尹青山 第 13 頁(yè) 共 22 頁(yè) 利用短時(shí)平均過(guò)零率還可以從背景噪聲中找出語(yǔ)音信號(hào)，可用于判斷寂靜無(wú)聲段和有聲 段的起點(diǎn)和終點(diǎn)位置。在孤立詞的語(yǔ)音識(shí)別中，必須要在一 連串連續(xù)的語(yǔ)音信號(hào)中進(jìn)行適當(dāng)分割，用以確定一個(gè)一個(gè)單詞的語(yǔ)音信號(hào)，即找出每一個(gè)單詞的開(kāi)始和終止位置，這在語(yǔ)音處理中是一個(gè)基本問(wèn)題。此時(shí)，在背景噪聲較小時(shí)用平均能量識(shí)別較為有效，而在背景噪聲較大時(shí)用平均過(guò)零率識(shí)別較為有效。但是研究表明，在以某些音為開(kāi)始或結(jié)尾時(shí)．如當(dāng)弱摩擦音 (如[f]、 [h]等音素 )、弱燃破音 (如 [p]、 [t]、 [k]等 音素 )為語(yǔ)音的開(kāi)頭或結(jié)尾；以鼻音(如 [n]、 [m]等音素 )為語(yǔ)音的結(jié)尾時(shí)．只用其中一個(gè)參量來(lái)判別語(yǔ)音的起點(diǎn)和終點(diǎn)是有困難