QUOTE(libintr @ Feb 15 2008, 13:20)

2 band for now...but i need many band implementation also....
For the 2-channel case:
'lp' and 'hp' are the QMF low-/highpass filter fair (impulse response). Usually only the lowpass filter coeffs are given. You can derive the highpass filter coeffs by alternating the signs a la hp = lp .* ((-1).^(1:length(lp)));
'x' = original signal
"analysis":
xlo = upfirdn(x,lp,1,2);
xhi = upfirdn(x,hp,1,2);
"snythesis":
x2 = 2 * (upfirdn(xlo,lp,2,1) + upfirdn(xhi,-hp,2,1));
Note1: '-hp' = hp(length(hp):-1:1) = synthesis filter in the special case for 2-channel QMF filters.
Note2: You may need to remove some elements of 'x2' at the start and the end so it matches x.
HTH,
SG