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

*　一維小波變換函數(shù)

*　說明：一維小波變換，只變換一次

*　輸入?yún)?shù)：input［］，輸入信號;output［］，小波變換結(jié)果，包括尺度系數(shù)和

*　小波系數(shù)兩部分;temp［］，存放中間結(jié)果;h［］，Daubechies小波基低通濾波器系數(shù);

*　g［］，Daubechies小波基高通濾波器系數(shù);n，輸入信號長度;m，Daubechies小波基緊支集長度。

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

voidDWT1D_1（doubleinput［］，doubleoutput0［］，doubleoutput1［］，doubletemp［］，doubleh［］，

doubleg［］，intn，intm）

{

//doubletemp［LENGTH］={0};//？？？？？？？？？？？？？

inti=0;

/*

//尺度系數(shù)和小波系數(shù)放在一起

Covlution（input，h，temp，n，m）;

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

{

output［i］=temp［i］;

}

Covlution（input，g，temp，n，m）;

for（i=1;i《n;i+=2）

{

output［i］=temp［i］;

}

*/

//尺度系數(shù)和小波系數(shù)分開

Covlution（input，h，temp，n，m）;

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

{

output0［i/2］=temp［i］;//尺度系數(shù)，進(jìn)行了2抽值，即尺度空間，低頻概貌部分

}

Covlution（input，g，temp，n，m）;

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

{

//output［i］=temp［i］;

output1［（i-1）/2］=temp［i］;//小波系數(shù)，已經(jīng)進(jìn)行了2抽取，即高頻細(xì)節(jié)部分

}

}

voidDWT1D_2（doubleinput［］，doubleAA2［］，doubleDA2［］，doubleAD2［］，doubleDD2［］，doubletemp0［］，doubletemp1［］，doubletemp［］，doubleh［］，

doubleg［］，intn，intm）

{

DWT1D_1（input，temp0，temp1，temp，h，g，n，m）;

inta1=（m+n）/2;

DWT1D_1（temp0，AA2，DA2，temp，h，g，a1，m）;

intd1=（n+m-4）/2;

DWT1D_1（temp1，AD2，DD2，temp，h，g，d1，m）;

}

voidDWT1D_3（doubleinput［］，doubleAAA3［］，doubleDAA3［］，doubleADA3［］，doubleDDA3［］，doubleAAD3［］，doubleDAD3［］，doubleADD3［］，doubleDDD3［］，

doubletemp0［］，doubletemp1［］，doubletemp2［］，doubletemp3［］，doubletemp00［］，doubletemp11［］，doubletemp［］，doubleh［］，doubleg［］，intn，intm）

{

DWT1D_2（input，temp0，temp1，temp2，temp3，temp00，temp11，temp，h，g，n，m）;

intaa2=（m+n）/4;

DWT1D_1（temp0，AAA3，DAA3，temp，h，g，aa2，m）;

intda2=（n+3*m-8）/4;

DWT1D_1（temp1，ADA3，DDA3，temp，h，g，da2，m）;

intad2=（n+3*m-4）/4;

DWT1D_1（temp2，AAD3，DAD3，temp，h，g，ad2，m）;

intdd2=（n+3*m-12）/4;

DWT1D_1（temp3，ADD3，DDD3，temp，h，g，dd2，m）;

}

voidmain（）

{

doubledata［LENGTH］;//輸入信號

doubletemp0［（LENGTH+DB_LENGTH）/4］;//先定義了一個中間結(jié)果

doubletemp1［（LENGTH+DB_LENGTH*3-8）/4］;

doubletemp2［（LENGTH+DB_LENGTH*3-4）/4］;

doubletemp3［（LENGTH+DB_LENGTH*3-12）/4］;

doubletemp00［（LENGTH+DB_LENGTH）/2］;

doubletemp11［（LENGTH+DB_LENGTH-4）/2］;

doubletemp［LENGTH+DB_LENGTH-1］;

doubleAAA3［（LENGTH+DB_LENGTH*5）/8］;//一維小波變換后的結(jié)果

doubleDAA3［（LENGTH+DB_LENGTH*5-16）/8］;

doubleADA3［（LENGTH+DB_LENGTH*7-8）/8］;

doubleDDA3［（LENGTH+DB_LENGTH*7-24）/8］;

doubleAAD3［（LENGTH+DB_LENGTH*7-4）/8］;

doubleDAD3［（LENGTH+DB_LENGTH*7-20）/8］;

doubleADD3［（LENGTH+DB_LENGTH*7-12）/8］;

doubleDDD3［（LENGTH+DB_LENGTH*7-28）/8］;

intaaa3=（LENGTH+DB_LENGTH*5）/8;//一維小波變換后的結(jié)果數(shù)組的長度

intdaa3=（LENGTH+DB_LENGTH*5-16）/8;

intada3=（LENGTH+DB_LENGTH*7-8）/8;

intdda3=（LENGTH+DB_LENGTH*7-24）/8;

intaad3=（LENGTH+DB_LENGTH*7-4）/8;

intdad3=（LENGTH+DB_LENGTH*7-20）/8;

intadd3=（LENGTH+DB_LENGTH*7-12）/8;

intddd3=（LENGTH+DB_LENGTH*7-28）/8;

intn=0;//輸入信號長度

intm=8;//Daubechies正交小波基長度

inti=0;

chars［32］;//從txt文件中讀取一行數(shù)據(jù)

/*//DB3

staticdoubleh［］={.332670552950，.806891509311，.459877502118，-.135011020010，

-.085441273882，.035226291882};

staticdoubleg［］={.035226291882，.085441273882，-.135011020010，-.459877502118，

.806891509311，-.332670552950};

*/

//DB4

staticdoubleh［］={0.2303778133088964，0.7148465705529154，0.6308807679398587，-0.0279837694168599，-0.1870348117190931，0.0308413818355607，0.0328830116668852，-0.0105974017850690};//h［］，Daubechies小波基低通濾波器系數(shù);

staticdoubleg［］={-0.0105974017850690，-0.0328830116668852，0.0308413818355607，0.1870348117190931，-0.0279837694168599，-0.6308807679398587，0.7148465705529154，-0.2303778133088964};//g［］，Daubechies小波基高通濾波器系數(shù)