mdwt.c
上传用户:speoil
上传日期:2022-06-23
资源大小:224k
文件大小:3k
源码类别:

波变换

开发平台:

Matlab

  1. /*
  2. File Name: mdwt.c
  3. Last Modification Date: 06/14/95 12:56:43
  4. Current Version: mdwt.c 1.5
  5. File Creation Date: Wed Oct 12 08:44:43 1994
  6. Author: Markus Lang  <lang@jazz.rice.edu>
  7. Copyright: All software, documentation, and related files in this distribution
  8.            are Copyright (c) 1994 Rice University
  9. Permission is granted for use and non-profit distribution providing that this
  10. notice be clearly maintained. The right to distribute any portion for profit
  11. or as part of any commercial product is specifically reserved for the author.
  12. Change History: Fixed code such that the result has the same dimension as the 
  13.                 input for 1D problems. Also, added some standard error checking.
  14. Jan Erik Odegard <odegard@ece.rice.edu> Wed Jun 14 1995
  15. MATLAB gateway for MDWT.c, discrete wavelet transform
  16. */
  17. #include <math.h>
  18. //#include <malloc.h>
  19. #include <stdio.h>
  20. #include "mex.h"
  21. #include "matrix.h"
  22. #define max(A,B) (A > B ? A : B)
  23. #define min(A,B) (A < B ? A : B)
  24. #define even(x)  ((x & 1) ? 0 : 1)
  25. #define isint(x) ((x - floor(x)) > 0.0 ? 0 : 1)
  26. void mexFunction(int nlhs,mxArray *plhs[],int nrhs,const mxArray *prhs[])
  27. {
  28.   double *x, *h,  *y, *Lf, *Lr;
  29.   int m, n, h_col, h_row, lh, L, i, po2, j;
  30.   double mtest, ntest;
  31.   /* check for correct # of input variables */
  32.   if (nrhs>3){
  33.     mexErrMsgTxt("There are at most 3 input parameters allowed!");
  34.     return;
  35.   }
  36.   if (nrhs<2){
  37.     mexErrMsgTxt("There are at least 2 input parameters required!");
  38.     return;
  39.   }
  40.   x = mxGetPr(prhs[0]);
  41.   n = mxGetN(prhs[0]); 
  42.   m = mxGetM(prhs[0]); 
  43.   h = mxGetPr(prhs[1]);
  44.   h_col = mxGetN(prhs[1]); 
  45.   h_row = mxGetM(prhs[1]); 
  46.   if (h_col>h_row)
  47.     lh = h_col;
  48.   else  
  49.     lh = h_row;
  50.   if (nrhs == 3){
  51.     L = (int) *mxGetPr(prhs[2]);
  52.     if (L < 0)
  53.       mexErrMsgTxt("The number of levels, L, must be a non-negative integer");
  54.   }
  55.   else /* Estimate L */ {
  56.     i=n;j=0;
  57.     while (even(i)){
  58.       i=(i>>1);
  59.       j++;
  60.     }
  61.     L=m;i=0;
  62.     while (even(L)){
  63.       L=(L>>1);
  64.       i++;
  65.     }
  66.     if(min(m,n) == 1)
  67.       L = max(i,j);
  68.     else
  69.       L = min(i,j);
  70.     if (L==0){
  71.       mexErrMsgTxt("Maximum number of levels is zero; no decomposition can be performed!");
  72.       return;
  73.     }
  74.   }
  75.   /* Check the ROW dimension of input */
  76.   if(m > 1){
  77.     mtest = (double) m/pow(2.0, (double) L);
  78.     if (!isint(mtest))
  79.       mexErrMsgTxt("The matrix row dimension must be of size m*2^(L)");
  80.   }
  81.   /* Check the COLUMN dimension of input */
  82.   if(n > 1){
  83.     ntest = (double) n/pow(2.0, (double) L);
  84.     if (!isint(ntest))
  85.       mexErrMsgTxt("The matrix column dimension must be of size n*2^(L)");
  86.   }
  87.   plhs[0] = mxCreateDoubleMatrix(m,n,mxREAL);
  88.   y = mxGetPr(plhs[0]);
  89.   plhs[1] = mxCreateDoubleMatrix(1,1,mxREAL);
  90.   Lr = mxGetPr(plhs[1]);
  91.   *Lr = L;
  92.   MDWT(x, m, n, h, lh, L, y);
  93. }