% fbanl2.m % usage: Y = fbanl2(X,l,h,s) % Performs an s-stage logartihmic DWT transform on X. (analysis filter bank) % It is assumed that l and h are smaller than #rows(X)/2^s and #cols(X)/2^s % and that 2^s divides #rows(X) and #cols(X) % time reversals of the masks l and h are used function Y=fbanl2(X,l,h,s) [m,n] = size(X); Y = X; for r = 0:(s-1) mr = m/2^r nr = n/2^r; Y(1:mr,1:nr) = hiloanl2(Y(1:mr,1:nr),l,h,l,h); end;