THESE/codes/wave/allcode/sfb.m

   1 function y = sfb(lo, hi, sf)
   2
   3 % Synthesis filter bank
   4 %
   5 % USAGE:
   6 %    y = sfb(lo, hi, sf)
   7 % INPUT:
   8 %    lo - low frqeuency input
   9 %    hi - high frequency input
  10 %    sf - synthesis filters
  11 %    sf(:, 1) - lowpass filter (even length)
  12 %    sf(:, 2) - highpass filter (even length)
  13 % OUTPUT:
  14 %    y - output signal
  15 % See also afb
  16 %
  17 % WAVELET SOFTWARE AT POLYTECHNIC UNIVERSITY, BROOKLYN, NY
  18 % http://taco.poly.edu/WaveletSoftware/
  19
  20 N = 2*length(lo);
  21 L = length(sf);
  22 lo = upfirdn(lo, sf(:,1), 2, 1);
  23 hi = upfirdn(hi, sf(:,2), 2, 1);
  24 y = lo + hi;
  25 y(1:L-2) = y(1:L-2) + y(N+[1:L-2]);
  26 y = y(1:N);
  27 y = cshift(y, 1-L/2);
  28