THESE/codes/wave/allcode/idualtree.m

   1 function y = idualtree(w, J, Fsf, sf)
   2
   3 % Inverse Dual-tree Complex DWT
   4 %
   5 % USAGE:
   6 %    y = idualtree(w, J, Fsf, sf)
   7 % INPUT:
   8 %    w - DWT coefficients
   9 %    J - number of stages
  10 %    Fsf - synthesis filters for the last stage
  11 %    sf - synthesis filters for preceeding stages
  12 % OUTUT:
  13 %    y - output signal
  14 % See dualtree
  15 %
  16 % WAVELET SOFTWARE AT POLYTECHNIC UNIVERSITY, BROOKLYN, NY
  17 % http://taco.poly.edu/WaveletSoftware/
  18
  19 % Tree 1
  20 y1 = w{J+1}{1};
  21 for j = J:-1:2
  22    y1 = sfb(y1, w{j}{1}, sf{1});
  23 end
  24 y1 = sfb(y1, w{1}{1}, Fsf{1});
  25
  26 % Tree 2
  27 y2 = w{J+1}{2};
  28 for j = J:-1:2
  29    y2 = sfb(y2, w{j}{2}, sf{2});
  30 end
  31 y2 = sfb(y2, w{1}{2}, Fsf{2});
  32
  33 % normalization
  34 y = (y1 + y2)/sqrt(2);
  35