% test multicomponent signal

x = v(6301:11500); % free fall portion, 5200 pts long
x = hilbert(x);    % 

x1 = freqshift(x(1:2048),.06,.25);
x2 = x(1576:1576+2047);
x2 = freqshift(x2,.25);
x3 = x(length(x)-2047:length(x));  % last 2048 pts
x3 = freqshift(x3,-.25,-.5);

subplot(221)
%contour(abs(tf(x1,50,50)))
%contour(abs(tf(x2,50,50)))
%contour(abs(tf(x3,50,50)))
contour(abs(tf(x1,80,80)))
grid
contour(abs(tf(x2,80,80)))
grid
contour(abs(tf(x3,80,80)))
grid

s = x1 + x2 + x3;
contour(abs(tf(s,80,80)))
grid

fix(clock)
S = clt(s,.5,300);
fix(clock)

save multicomponentclt

