void dwt(double g[],double h[],int wlen,double c[],double d[],int m,int sca[]);//小波分解
void idwtc3(double h[],int wlen,double c3[],int m,int sca[]);///小波單支重構
void idwtd3(double g[],int wlen,double d3[],int m,int sca[]);///
void idwtd2(double g[],int wlen,double d2[],int m,int sca[]);///
void idwtd1(double g[],int wlen,double d1[],int m,int sca[]);////

void trmul(double a[],double b[],int m,int n,int k,double c[]);  //矩陣相乘
int ssgj(double a[],int n);           //對稱正定矩陣求逆

void netc3();
//void netd3();
//void netd2();
//void netd1();