一種改進(jìn)的小波閾值信號(hào)去噪方法

摘要 為改進(jìn)濾波效果，提高去噪質(zhì)量。通過(guò)分析軟硬閾值去噪的原理和方法，為小波闞值信號(hào)處理提出了一種改進(jìn)的去噪方法。該方法綜合了軟硬閾值的特點(diǎn)，對(duì)其參數(shù)進(jìn)行優(yōu)化設(shè)計(jì)，通過(guò)調(diào)節(jié)參數(shù)值以更好地獲得閾值估計(jì)。針對(duì)改進(jìn)后的去噪算法，通過(guò)Matlab仿真比較了傳統(tǒng)的小波軟硬閾值算法與該算法的消噪效果，結(jié)果表明，提出的方法有更好的消噪效果和穩(wěn)定性。
關(guān)鍵詞 小波閾值去噪；閾值函數(shù)；信噪比；均方誤差

小波變換與傅里葉變換、窗口傅里葉變換相比，它是一個(gè)時(shí)間和頻率的局部變換，因而能有效地從信號(hào)中提取信息，通過(guò)伸縮和平移等運(yùn)算功能對(duì)函數(shù)或信號(hào)進(jìn)行多尺度細(xì)化分析，解決了傅里葉變換不能解決的許多問(wèn)題。
小波變換被譽(yù)為“數(shù)學(xué)顯微鏡”，是調(diào)和分析發(fā)展史上里程碑式的進(jìn)展。小波理論被認(rèn)為是對(duì)傅里葉分析的重大突破，是近年來(lái)信號(hào)處理領(lǐng)域的研究熱點(diǎn)，許多學(xué)者將小波在理論上的研究成果應(yīng)用到諸如圖像壓縮、特征提取、信號(hào)濾波和數(shù)據(jù)融合等方面。小波之所以在信號(hào)處理領(lǐng)域具有很大的優(yōu)勢(shì)，在于小波變換可以獲得信號(hào)的多分辨率描述，同時(shí)，小波變換具有豐富的小波基可適應(yīng)不同特征的信號(hào)。隨著小波理論的發(fā)展，Mallat提出了模極大值重構(gòu)濾波，Xu提出了空域相關(guān)濾波和Donoho提出了小波域閾值濾波來(lái)消除噪聲。一般地，不同性質(zhì)的噪聲需要采用不同的消噪方法處理。
上述算法都存在參數(shù)的選取問(wèn)題，不同的參數(shù)選取對(duì)濾波的效果會(huì)有一些差異。Donoho的傳統(tǒng)小波閾值去噪方法的實(shí)現(xiàn)最簡(jiǎn)單、計(jì)算量最小。但其在理論上找到的最優(yōu)通用閾值，實(shí)際應(yīng)用中效果并不理想。因此，文中在Donoho的傳統(tǒng)小波閾值去噪方法的基礎(chǔ)上，改變了小波閾值函數(shù)的部分參數(shù)得到了一種新的小波閾值函數(shù)，既避免了硬閾值函數(shù)的不連續(xù)性，又可以自適應(yīng)去噪，從而有效地保存了信號(hào)的邊緣信息。仿真結(jié)果表明，改進(jìn)后方法有更好的去噪性能。

1 小波閾值的去噪原理
小波變換具有一種“集中”的能力。信號(hào)經(jīng)小波變換后，可以認(rèn)為由信號(hào)產(chǎn)生的小波系數(shù)包括有信號(hào)的重要信息，其幅值較大，但數(shù)目較小，而噪聲對(duì)應(yīng)的小波系數(shù)幅值小。通過(guò)在不同尺度上選取一合適的閾值，并將小于該閾值的小波系數(shù)置零，而保留大于閾值的小波系數(shù)，從而使信號(hào)中的噪聲得到有效的抑制，最后進(jìn)行小波逆變換，得到濾波后的重構(gòu)信號(hào)。
小波去噪的基本思路如圖1所示。信號(hào)先經(jīng)過(guò)預(yù)處理，然后利用小波變換講信號(hào)分解到多尺度上，再對(duì)每一層小波系數(shù)進(jìn)行閾值處理，最后對(duì)處理后的小波系數(shù)進(jìn)行信號(hào)重構(gòu)。

設(shè)一個(gè)含噪聲的一維信號(hào)模型可以表示為f(i)=s(i)+n(i)，其中s(i)為原始信號(hào)，n(i)為方差σ2的高斯白噪聲，服從N(0，σ2)。
1．1 小波硬閾值去噪的步驟
(1)對(duì)信號(hào)求小波變換。
(2)除了最粗尺度信號(hào)外，將各細(xì)節(jié)信號(hào)作閾值處理，閾值t取為，當(dāng)某位置小波變換值大于閾值時(shí)，保留原值，否則置零，用公式表示為

(3)利用小波變換重構(gòu)，求出信號(hào)的濾波值。
1．2 小波軟閾值去噪的步驟
(1)對(duì)信號(hào)求小波變換。
(2)除了最粗尺度信號(hào)外，將各細(xì)節(jié)信號(hào)作閾值處理，閾值t取為，當(dāng)某位置小波變換大于閾值時(shí)，向著減小系數(shù)幅值的方向作一個(gè)收縮t，否則置零，用公式表示為

其中，sgn(x)為符號(hào)函數(shù)。