C语言实现97提升小波变换代码分享

#include <stdlib.h>    
#include <stdio.h>    
#include <math.h>    
#include <string.h>    
#include <time.h>    
   
#define F_EXTPAD 4    
#define D_EXTPAD 2    
   
/*  float 9/7  
you can get the coefficient by using [a,b,c,d]=WFILTERS('bior4.4') in matlab  
********************************/   
double h[5] ={ 0.85269867900889 ,  0.37740285561283 , -0.11062440441844,  -0.02384946501956 ,  0.03782845550726 };   
double g[4] ={ -0.78848561640558 ,  0.41809227322162 , 0.04068941760916,  -0.06453888262870  };   
double h_b[4] ={ 0.78848561640558 ,  0.41809227322162 , -0.04068941760916,  -0.06453888262870 };   
double g_b[5] ={-0.85269867900889 ,  0.37740285561283 , 0.11062440441844,  -0.02384946501956 ,  -0.03782845550726  };   
   
/*********  float 5/3 **********  
you can get the coefficient by using [a,b,c,d]=WFILTERS('bior2.2') in matlab  
**********************************  
double h[5] ={  1.06066017177982 ,  0.35355339059327  ,-0.17677669529664 ,0};  
double g[4] ={ -0.70710678118655 ,  0.35355339059327 , 0 };  
double h_b[4] ={ 0.70710678118655 ,  0.35355339059327, 0 };  
double g_b[5] ={-1.06066017177982 ,  0.35355339059327,   0.17677669529664,0 };  
*/   
static float   *x_alloc = NULL; // **    
   
unsigned char **Alloc_Char_Matrix(short r, short c)   
{   
   unsigned char  *x, **y;   

   short           n;   
   // calloc    
   x = (unsigned char *) calloc((long) r * c, sizeof(unsigned char));   
   y = (unsigned char **) calloc(r, sizeof(unsigned char *));   
   for (n = 0; n <= r - 1; n++)   
       y[n] = &x[(long) c * n];   
   return y;   
}   
   
short **Alloc_Short_Matrix(short r, short c)   
{   
   short          *x, **y;   
   short           n;   
   // calloc    
   x = (short *) calloc((long) r * c, sizeof(short));   
   y = (short **) calloc(r, sizeof(short *));   
   for (n = 0; n <= r - 1; n++)   
      y[n] = &x[(long) c * n];   
   return y;   
}   
   
float **Alloc_Float_Matrix(short r, short c)   
{   
   float          *x, **y;   
   short           n;   
   // calloc    
   x = (float *) calloc((long) r * c, sizeof(float));   
   y = (float **) calloc(r, sizeof(float *));   
   for (n = 0; n <= r - 1; n++)   
       y[n] = &x[(long) c * n];