function y = wavelet( x ) % WAVELET Take Daubechies 4 tap wavelet transform [ m, n ] = size( x ); y = x; while ( n >= 4 ) y( 1:n ) = daubech4( y( 1:n ) ); n = n / 2; end