1、小波閾值去噪理論

小波閾值去噪就是對(duì)信號(hào)進(jìn)行分解，然后對(duì)分解后的系數(shù)進(jìn)行閾值處理，最后重構(gòu)得到去噪信號(hào)。該算法其主要理論依據(jù)是：小波變換具有很強(qiáng)的去數(shù)據(jù)相關(guān)性，它能夠使信號(hào)的能量在小波域集中在一些大的小波系數(shù)中;而噪聲的能量卻分布于整個(gè)小波域內(nèi)。因此，經(jīng)小波分解后，信號(hào)的小波系數(shù)幅值要大于噪聲的系數(shù)幅值。可以認(rèn)為，幅值比較大的小波系數(shù)一般以信號(hào)為主，而幅值比較小的系數(shù)在很大程度上是噪聲。于是，采用閾值的辦法可以把信號(hào)系數(shù)保留，而使大部分噪聲系數(shù)減小至零。小波閾值收縮法去噪的具體處理過(guò)程為：將含噪信號(hào)在各尺度上進(jìn)行小波分解，設(shè)定一個(gè)閾值，幅值低于該閾值的小波系數(shù)置為0，高于該閾值的小波系數(shù)或者完全保留，或者做相應(yīng)的“收縮（shrinkage）”處理。最后將處理后獲得的小波系數(shù)用逆小波變換進(jìn)行重構(gòu)，得到去噪后的信號(hào).

2、小波閾值去噪c語(yǔ)言程序

此程序是用于信號(hào)處理分析，突出奇異值的前段處理，對(duì)信號(hào)進(jìn)行小波包分解，用C語(yǔ)言實(shí)現(xiàn)的，僅供參考。

#include

#include

#include

#include

#defineLENGTH4096//信號(hào)長(zhǎng)度

#defineDB_LENGTH8//Daubechies小波基緊支集長(zhǎng)度

/******************************************************************

*　一維卷積函數(shù)

*　說(shuō)明：循環(huán)卷積，卷積結(jié)果的長(zhǎng)度與輸入信號(hào)的長(zhǎng)度相同

*　輸入參數(shù)：data［］，輸入信號(hào);core［］，卷積核;cov［］，卷積結(jié)果;

*n，輸入信號(hào)長(zhǎng)度;m，卷積核長(zhǎng)度。

******************************************************************/

/*voidCovlution（doubledata［］，doublecore［］，doublecov［］，intn，intm）

{

inti=0;

intj=0;

intk=0;

//將cov［］清零

for（i=0;i《n;i++）

{

cov［i］=0;

}

//前m/2+1行

i=0;

for（j=0;j《m/2;j++，i++）

{

for（k=m/2-j;k《m;k++）

{

cov［i］+=data［k-（m/2-j）］*core［k］;//k針對(duì)core［k］

}

for（k=n-m/2+j;k《n;k++）

{

cov［i］+=data［k］*core［k-（n-m/2+j）］;//k針對(duì)data［k］

}

}

//中間的n-m行

for（i=m/2;i《=（n-m）+m/2;i++）

{

for（j=0;j《m;j++）

{

cov［i］+=data［i-m/2+j］*core［j］;

}

}

//最后m/2-1行

i=（n-m）+m/2+1;

for（j=1;j《m/2;j++，i++）

{

for（k=0;k《j;k++）

{

cov［i］+=data［k］*core［m-j-k］;//k針對(duì)data［k］

}

for（k=0;k《m-j;k++）

{

cov［i］+=core［k］*data［n-（m-j）+k］;//k針對(duì)core［k］

}

}

}

*/

//定義一個(gè)線性卷積

voidCovlution（doubledata［］，doublecore［］，doublecov［］，intn，intm）

{

inti=0;

intj=0;

intt=0;

//將cov［］清零

for（j=0;j《n+m-1;j++）

{

cov［j］=0;

}

for（j=0;j《m+n-1;j++）

{

if（j《=m-1）//前面m行

{

for（i=0，t=j;t》=0;i++，t--）

cov［j］+=data［i］*core［t］;

}

elseif（j《=n-1）//中間n-m行

{

for（i=j-m+1，t=m-1;t》=0;i++，t--）

cov［j］+=data［i］*core［t］;

}

else//后面m行

{

for（i=j-m+1，t=m-1;i《n;i++，t--）

cov［j］+=data［i］*core［t］;

}

}

}