模極大值與閾值決策融合的小波語音數(shù)據(jù)去噪方法

 收稿日期：2007-10-10；

修回日期：2008-03-07

基金項目：國家自然科學(xué)基金資助項目（70371008）



作者簡介：劉健（1975-），男，上海人，博士研究生，主要研究方向為信息資源管理與智能決策(sockl@163.com)；陶玉靜 (1977-)，女，講師，博士，主要研究方向為非線性系統(tǒng)仿真；張維明 (1962-)，男，教授 博導(dǎo)，博士，主要研究方向為信息資源管理與智能決策.

(國防科技大學(xué) a.信息系統(tǒng)與管理學(xué)院; b.武器裝備發(fā)展研究中心 長沙 410073)

摘要：小波在時域、頻域都具有表征信號局部特征的能力，對于語音數(shù)據(jù)的去噪，與其他方法比較具有優(yōu)勢，但是閾值的選擇缺乏明確的指標(biāo)。模極大值法保留幅度隨尺度增加而增大的點（對應(yīng)于有用信號的極值點）來達(dá)到去噪的目的，但將非極值處變換系數(shù)置零，造成信號損失。提出了一種變化模極大值與閾值決策相融合的去噪方法，有效地結(jié)合了兩者的優(yōu)點。經(jīng)實驗證明，該方法精度高、算法簡單、計算速度快。

關(guān)鍵詞：小波去噪； 模極大值； 閾值； 語音數(shù)據(jù)處理

中圖分類號：TP391

文獻(xiàn)標(biāo)志碼：A

文章編號：1001-3695(2008)10-3134-02

Modulus maximum and nonlinear threshold based wavelet-denoising method

LIU Jiana TAO Yu-jingb ZHANG Wei-minga

（a.College of Information System Management b.Institute of Equipment Development National University of Defense Technology Changsha 410073 China）

Abstract:This paper proposed a new wavelet-denoising method which was based on wavelet transform modulus maximum and nonlinear threshold denoising method on wavelet. This method absorbed the advantages and overcomed the disadvantages of traditional methods. Take an experiment to validate the improved method which shows the new method has the advantage of the high precision and the simple arithmetic and the high computing efficiency.

Key words：wavelet-denoising; maximum modulus; threshold; speech data processing

小波變換能有效完成信號的時間、空間局部化，具有多分辨率、多尺度的特點，可以由粗到精地逐步觀察信號；同時還具有品質(zhì)因數(shù)恒定，即相對帶寬(帶寬與中心頻率之比)恒定的特點，適當(dāng)?shù)剡x擇基本小波，便可以利用其在時、頻域都具有表征信號局部特征的能力，因此能有效區(qū)分信號中的突變部分和噪聲，從而實現(xiàn)對信號的去噪。利用小波分析去噪，即是在不同尺度下進(jìn)行小波變換，其實質(zhì)就是用不同中心頻率的帶通濾波器對信號進(jìn)行濾波，將那些主要反映噪聲性質(zhì)的成分去掉，從而得到質(zhì)量較好的觀測數(shù)據(jù)。

本文通過對現(xiàn)有小波去噪方法進(jìn)行改進(jìn)，提出一種基于小波變換模極大值與閾值決策相融合的去噪方法，并用經(jīng)典的Lorenz非線性動力系統(tǒng)驗證該方法的有效性，適合對復(fù)雜環(huán)境下錄制的語音數(shù)據(jù)進(jìn)行去噪處理，為后續(xù)分析提供干凈、真實的信號。

1小波去噪方法改進(jìn)

目前，常用的小波去噪方法有非線性小波變換閾值法、平移不變量法和小波變換模極大值法。

利用小波變換閾值法去噪，閾值的選取是一個關(guān)鍵問題[1，2]，其直接影響著去噪效果。閾值選取過高，會過多地將信號當(dāng)做噪聲去掉；閾值選取過低，則保留的噪聲信號過多，影響信號的進(jìn)一步分析。但閾值的確定往往缺乏明確的指標(biāo)，很大程度上依賴于經(jīng)驗，沒有可以掌握的定量標(biāo)準(zhǔn)。

小波變換模極大值去噪法是基于信號的模極大值隨尺度的增大而增加，噪聲的模極大值隨尺度的增大而減少，從而去除幅度隨尺度增加而減小的點（對應(yīng)噪聲的極值點），保留幅度隨尺度增加而增大的點（對應(yīng)于有用信號的極值點）來達(dá)到去噪的目的。經(jīng)模極大值法去噪后，小波變換系數(shù)僅剩下模極大值點處的值，而其余部分都被置為零。因此僅僅通過這有限個模極大值點直接重構(gòu)信號，誤差是很大的，必須找到一種有效的方法，可以根據(jù)這幾個模極大值點求出原始小波變換的一個近似。Mallat[3]提出的交替投影算法給出了一種逼近小波系數(shù)的方法，但由于其計算量很大、程序復(fù)雜、速度較慢，而且其計算過程可能不穩(wěn)定，不便于實際應(yīng)用。

基于上述原因，出現(xiàn)了很多小波去噪的改進(jìn)方法。文獻(xiàn)[4]基于混沌信號具有分形結(jié)構(gòu)的特點，提出了一種分?jǐn)?shù)維與小波變換相結(jié)合的帶觀測噪聲混沌信號的噪聲去除新方法，利用噪聲小波變換特性提取信號的真實分?jǐn)?shù)維，然后根據(jù)預(yù)處理混沌時間序列分?jǐn)?shù)維與信號真實分?jǐn)?shù)維之差控制小波包除噪的閾值及位置，達(dá)到噪聲與信號的自適應(yīng)分離。雖然去噪精度得到了提高，但需要反復(fù)計算分?jǐn)?shù)維和調(diào)整小波分解層數(shù)，使得計算過于復(fù)雜。文獻(xiàn)[5]根據(jù)不同尺度下小波信號的信噪比靈活地選擇閾值，但噪聲主要集中在哪些尺度上沒有定量的標(biāo)準(zhǔn)。

鑒于此，本文提出了一種基于小波變換模極大值與閾值決策相融合的去噪方法。具體算法如下：

a)對含噪信號進(jìn)行離散小波變換或小波包變換，得到不同尺度上的小波系數(shù)，即信號經(jīng)過小波分解后在不同尺度上的高通分量。

b)計算各個尺度上的小波變換模極大值。

c)計算各尺度上小波模極大值系數(shù)的功率。

Pj(x)=∑k∈Z｜W2j，k｜2 j=1，2，…

(1)

低尺度下，功率Pj(x)主要由噪聲控制，隨尺度的增加，噪聲變換模極大值逐漸減小，而信號變換模極大值逐漸增大，故功率因噪聲的影響而快速下降，首次達(dá)最小時對應(yīng)的尺度jm可作為尺度取舍的一個閾點。文獻(xiàn)[4]中令W2j，k=0(j=1，2，… jm)，此時不僅去除了噪聲，也去除了部分有用信號成分。本文融合閾值決策對此進(jìn)行了改進(jìn)，認(rèn)為信號的噪聲主要存在于小于jm的尺度上，這些尺度上噪聲能量較大，信噪比很低。以尺度jm作為分界點對小波系數(shù)分別進(jìn)行適當(dāng)?shù)拈撝氮┨幚怼＊オ?/p>

d)對于最大尺度J上的細(xì)節(jié)信號，其信噪比高，有用信號的能量較大，占主要部分，因此閾值的選取不應(yīng)太大，以免去除過多的有用信號。選取閾值為

tJ=σJ2 log N/J(2)

其中：tJ為相應(yīng)的閾值；σJ為信號在最大尺度上的方差，N為信號的長度。

e)對尺度j=jm+1，jm+2，…，J-1，信噪比較高，閾值應(yīng)當(dāng)適當(dāng)提高，用下式選取閾值：

tj=σ2 log N/ln (j+1)(3)

其中:σ為信號在最大尺度上的方差。隨著尺度j的增大，tj的值逐漸減小，使得其與噪聲在小波變換各個尺度上的傳播特性相一致。

f)對尺度j=1，2，…，jm，噪聲的能量較大，信噪比很低，因此用Donoho等人[1]提出的廣義閾值：

t^=σ^2 log N(4)

去除此部分的噪聲。其中：σ^=m/0.674 5；m為最高尺度J的絕對變差中值。

2小波去噪改進(jìn)方法仿真驗證

選用經(jīng)典的Lorenz系統(tǒng)來驗證上述方法，選取參數(shù)σ=10，r=28，b=8/3，用4階Runge-Kutta法進(jìn)行數(shù)值計算，步長為0.01，產(chǎn)生10 000個點的信號，用中間2 000個點疊加噪聲信號（均勻分布白噪聲和高斯分布白噪聲），原始信號及疊加噪聲后信號分別如圖1、2所示。

經(jīng)改進(jìn)算法步驟a)~c)處理后結(jié)果如圖3所示，可見其小波系數(shù)功率Pj(x)由于噪聲的影響隨尺度的增大而下降，在j為4時首次達(dá)最小。筆者認(rèn)為噪聲主要集中在前四個尺度中，以此為判據(jù)利用步驟d)~f)進(jìn)行去噪，結(jié)果如圖4所示。

圖5、6分別為利用模極大值法和閾值法去噪的結(jié)果。與圖4比較可以發(fā)現(xiàn)，利用模極大值法去噪后，由于非模極大值點的小波系數(shù)全置為零，使得有用信息有很大的損失；經(jīng)閾值法去噪后，由于閾值設(shè)置的隨意性，仍存在較多的噪聲；而利用本文提出的改進(jìn)方法去噪，既達(dá)到了良好的去噪效果，又保留了原始信號，從而驗證了改進(jìn)方法的有效性。

3語音數(shù)據(jù)的小波去噪

上文已經(jīng)驗證了基于小波變換模極大值與閾值決策相融合的小波去噪方法的有效性，下面采用此方法對語音數(shù)據(jù)進(jìn)行去噪處理，排除噪聲干擾，以提取數(shù)據(jù)中反映本質(zhì)特性的真實信號用于語音識別。

以圖7所示的帶有白噪聲的真實語音數(shù)據(jù)為對象，進(jìn)行去噪處理。

首先，采用sym6小波函數(shù)，將信號分解為六層；然后計算各尺度上小波系數(shù)模極大值系數(shù)的功率。結(jié)果如圖8所示。從圖中可以看出，j=4時功率達(dá)最小，即在前四個尺度上信號主要由噪聲控制；j＞4后功率增大，主要由真實信號控制。因此 j=1~4時，去噪閾值采用式(3)計算； j=5、6時，去噪閾值采用式(2)計算。去噪后的小波系數(shù)如圖9所示。可以看出，改進(jìn)方法不像模極大值法那樣將1~4層的小波系數(shù)全置為零，只有1、2層的小波系數(shù)經(jīng)閾值計算后最終全置為零，而3、4層的小波系數(shù)有一部分得到了保留，從而既達(dá)到了去噪的目的，又確保了相應(yīng)的有用信息。

圖10為實際語音圖，同時采用模極大值法去噪后的信號如圖11所示。最后，將各尺度上去噪后的小波系數(shù)進(jìn)行逆變換，重構(gòu)出去噪后的信號，如圖12所示。可以看出，利用本文提出的改進(jìn)方法達(dá)到了良好的去噪效果，保留了原信號的主要特征。

4結(jié)束語

本文提出的基于小波變換模極大值與閾值決策相融合的去噪方法，融合了小波變換模極大值去噪法和閾值去噪法的優(yōu)點，針對不同尺度下模極大值的情況，進(jìn)行合理的閾值設(shè)置，解決了單純利用模極大值法需要尋找各尺度上各個模極大值點的復(fù)雜點，也即尋找模極大值線的麻煩；克服了將非模極大值點的小波系數(shù)置為0從而舍棄部分有用信息的缺陷，并且經(jīng)實際應(yīng)用在數(shù)字校園語音平臺中，取得了理想的效果。

參考文獻(xiàn)：

[1]DONOHO D L JOHNSTONE. Ideal spatial adaptation by wavelet shrinkage[J].Biometrika 1994，81(3): 425-455.

[2]MEYER Y. Wavelets and operators[M]. Cambridge:-Cambridge-University Press，1992:133-156.

[3]MALLAT S. Characterization of signals from multiscale edges[J].IEEE Trans on PAMI 1992，14(7):710-732.

[4]馬麗萍 石炎福 余華瑞. 含噪聲混沌信號的小波去噪方法研究[J].信號處理，2002，18(1):83-87.

[5]韓敏 劉玉花 席劍輝，等. 基于小波變換閾值決策的混沌信號去噪研究[J].信息與控制，2005 34(5):33-37.

[6]楊福生. 小波變換的工程分析與應(yīng)用[M].北京：科學(xué)出版社，2000.

[7]顏龍，劉剛，郭軍. 基于小波分析的大詞匯漢語連續(xù)語音識別系統(tǒng)魯棒性的研究[J].中文信息學(xué)報，2006，20(2):60-65.