二、小波分析算法的C語言實(shí)現(xiàn) - 小波分析算法的公式與C語言實(shí)現(xiàn)

　　二、小波分析算法的C語言實(shí)現(xiàn)

　　小波變換程序：

　　voidDB4DWT（doubleData［］，intn）

　　{

　　if（n》=4）

　　}

　　inti，j;

　　intbalf=n》》1;

　　double*tmp=newdouble［n］;

　　i=0;

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

　　{

　　tmp［j］=Data［（2*j）%n］*h0+

　　Data（2*j+l）%n］*h1+

　　Data［（2*j+2）%n］*h2+

　　Data［（2*j+3）%n］*h3;

　　tmp［j+half］=Datal（2*j）%n］*g0+

　　Data［（2*j+I）%n］*g1+

　　Datal（2*j+2）%n］*g2+

　　Data［（2*j+3）%n］*g3;

　　}

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

　　{

　　Data［i］=tmp［i］;

　　}}}

　　提升算法的程序

　　小波的分解算法的提升過程如下：

　　使用提升算法的程序:

　　voidDB4LiftDWT（doubleData［］，intn）

　　應(yīng)用研究

　　inti，half=n》》1;

　　double*pS=newdouble［half］;

　　//臨時(shí)變量存放平滑系數(shù)

　　double*pD=new

　　double［half］;

　　//臨時(shí)變量存放細(xì)節(jié)系數(shù)

　　for（i=0;i《half;it+）/賦值

　　{

　　pS［j］=Data［2*i］：//even，

　　pD［i］=Data［2*i+1］：//odd

　　}

　　for（i=0;i《half;i++）11DB4變換Update1

　　{

　　pS［i］=pS［i］+pD［i］*sqrt_3;

　　}

　　for（i=0;i《half;i++）//DB4的predict

　　{

　　pD［i］=pD［i］-sqrt_3*p$［i14

　?。╯qt_3-2）*pS［（-1）》=0？（i-1）：（halfti-1）14;

　　//邊界是采用周期延拓

　　}

　　for（i=0;i《half;i++）//DB4的update2

　　{

　　pS［i］=pS［i］-pD［（i+1）%half］;

　　/1邊界采用周期延拓

　　}

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

　　{

　　p$［i］=（sqt_3-1）*pS［iV（sqnt_2）;/比例系數(shù)

　　pD［i］=（sqrt_3+1）*pD［i］（sqrt_2）;

　　}

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

　　{

　　Data［i］=pS［i］;//將平滑系數(shù)放回原數(shù)組

　　}

　　for（i=half;i《length;i++）

　　{

　　Data［i］=pD［i-half］;

　　/將細(xì)節(jié)系數(shù)放回原數(shù)組

　　}

　　deleteOpD;

　　delete［DpS;

　　在目述程序中用到了兩個(gè)臨時(shí)變量數(shù)組，雖然程

　　序結(jié)構(gòu)清楚，容易閱讀但消耗較大的內(nèi)存空間，這里給

　　出另一種方法不用輔助數(shù)組，帶來節(jié)約內(nèi)存，但耗時(shí)。

　　voidsplit（doubleData［］，intn）

　　{intstart=1;

　　intend=n-1;

　　while（start《end）

　　{for（inti=start;i《end;i=i+2）

　　{doubletmp=Data［i］;

　　Data［i］=Data［i+1］;

　　Data［i+l］=tmp;}

　　start=start+1;

　　end=end-1;}}

　　對原數(shù)組調(diào)用split函數(shù)進(jìn)行處理后，則Data數(shù)組的前半部分是偶數(shù)索引值，后半部分是奇數(shù)索引值。再對提升算法的程序進(jìn)行一些修正就可以起到節(jié)約內(nèi)存的目的。具體的實(shí)現(xiàn)限于篇幅不多介紹。利用提升算法的小波程序結(jié)構(gòu)清楚，與Mallat算法相比運(yùn)算量也?。?］。同時(shí)它的反變換也很易實(shí)現(xiàn)，這里由于篇幅限制，不對反變換作多的介紹。

　　本文給出了一維小波變換的C的實(shí)現(xiàn)，可在基礎(chǔ)上實(shí)現(xiàn)：二維的小波變換，希望對要自己實(shí)現(xiàn)小波變換的讀者有一定的幫助。