THESE/codes/wave/allcode/dwt.m

   1 function w = dwt(x, J, af)
   2
   3 % Discrete 1-D Wavelet Transform
   4 %
   5 % USAGE:
   6 %    w = dwt(x, J, af)
   7 % INPUT:
   8 %    x - N-point vector, where
   9 %            1) N is divisible by 2^J
  10 %            2) N >= 2^(J-1)*length(af)
  11 %    J - number of stages
  12 %    af - analysis filters
  13 %    af(:, 1) - lowpass filter (even length)
  14 %    af(:, 2) - highpass filter (evenlength)
  15 % OUTPUT:
  16 %    w{j}, j = 1..J+1 - DWT coefficients
  17 % EXAMPLE:
  18 %    [af, sf] = farras;
  19 %    x = rand(1,64);
  20 %    w = dwt(x,3,af);
  21 %    y = idwt(w,3,sf);
  22 %    err = x - y;
  23 %    max(abs(err))
  24 %
  25 % WAVELET SOFTWARE AT POLYTECHNIC UNIVERSITY, BROOKLYN, NY
  26 % http://taco.poly.edu/WaveletSoftware/
  27
  28 for k = 1:J
  29     [x w{k}] = afb(x, af);
  30 end
  31 w{J+1} = x;
  32