C程序小波分析，使用DB4小波

【共享】贴一个自己写的小波变换的C源代码
#define N0 128
#include "stdio.h"
#include "stdlib.h"
#include "math.h"
#include "string.h"
void db4(double *h,double *g,double *hh,double *gg);
void wd(int N,double *h,double *g,double *c0,double *c,double *d);
void wr(int N,double *h,double *g,double *c, double *d,double *cd);
void main()
{
      double fk[N0],c0[N0],c[N0],d[N0];
      double h[8],g[8],hh[8],gg[8];
      float fk0[N0];
      FILE *fp;
      int i,k,j,n,l,N;
      fp=fopen("wdata.dat","rt");
      fscanf(fp,"%d",&N);
      for(k=0;k<N;k++)     fscanf(fp,"%f",&fk0[k]);
      fclose(fp);
      db4(h,g,hh,gg);
      for(k=0;k<N;k++)  {
           c0[k]=fk0[k];
           c[k]=0;
           d[k]=0;
       }
       wd(N,hh,gg,c0,c,d);
       wr(N,hh,gg,c,d,c0);
       for(k=0;k<N;k++)  printf("k=%d  c0=%f  c=%f\n",k,fk0[k],c0[k]);
       return;