my_display_wavelet_ti.m

This document is intended to serve as an introduction to Wavelet processing through a set of Matlab

function my_display_wavelet_ti(f,T)

f = f(:);
n = size(f,1);
q = 5;
options.ti = 1;
fwti = perform_wavelet_transform(f,log2(n)-q,1,options);

if nargin==2
    fwti = fwti .* (abs(fwti)>1.1*T);
    f = perform_wavelet_transform(fwti,log2(n)-q,-1,options);
end

clf;
subplot(q+1,1,1);
plot(f); axis tight; ylabel('Original');
a = 1;
for i=[1 q:-1:2]
    a = a+1;
    subplot(q+1,1,a);
    plot(fwti(:,i)); axis tight; 
    if a>2
        ylabel(['scale ' num2str(a-2)] );
    else
        ylabel('coarse scale');
    end
end