% test dwtr
close all
clear all
%-------------------------
data = [ 1 0 -3 2 1 0 1 2 ];
H_haar = [1/sqrt(2) 1/sqrt(2)];
bar(data)
title('Original Data')
pause

aa = dwtr( data, 1,H_haar )
bar(aa)
title('Haar Transform: One Step')
pause
ab = dwtr( data, 2,H_haar )
bar(ab)
title('Haar Transform: Two Steps')
pause
ac = dwtr( data,3, H_haar )
bar(ac)
title('Haar Transform: Three Steps')
pause
clear data; 
clear aa; 
t=linspace(0,1,512);
data = sqrt(t.*(1-t)).*sin((2*pi*1.05) ./(t+.05));
filt = [ -0.07576571478934  -0.02963552764595  ...
          0.49761866763246   0.80373875180522  ...
          0.29785779560554  -0.09921954357694  ...
         -0.01260396726226   0.03222310060407]; 
aa = dwtr( data, 3, filt);
subplot(2,1,1)
plot(data)
axis([1 512 -2 2])
title('Doppler')
subplot(2,1,2)
plot(aa, '-');
axis([1 512 -2 2])
title('Wavelet Transform of Doppler by Symmlet 8')