clear;
% install Wavelab http://www-stat.stanford.edu/~wavelab/Wavelab_850/download.html
%and put the path

cd '/home/leonid/lectures/labs/Wavelab850'%path to wavelab
WavePath

N_signal=1024;
% generating two-sin signal
for (k=1:1:N_signal)
    if k<(N_signal/2)
      garm(k)=5*sin(2*pi/20*(k-1));
    else
      garm(k)=5*sin(2*pi/50*(k-1));
    end;
      chirp(k)=5*sin(2*pi/5000*(k-1)^2);

end;

plot(garm);

% ARMA process generating
ar(1)=0.5*randn(1);
ar(2)=-0.2*ar(1)+0.5*randn(1);

for (i=3:1:N_signal)
    ar(i)=0.9*ar(i-1)-0.7*ar(i-2)+0.5*randn(1);
end;

plot(ar);

signal=garm+ar+chirp;

plot(signal);

plot(chirp);
cd 'Orthogonal';

QMF_Filter  = MakeONFilter('Coiflet',2);
wcNaive      = FWT_PO(signal, 5, QMF_Filter);

%PlotWaveCoeff(wcNaive,1,0.5)
DisplayMultiRes(wcNaive,1,0.5,QMF_Filter,'Mesh')
%--------------------------------------------------------


specgm = WindowFT(signal,100,20,'Hamming','wft');
%--------------------------------------------------------

rwt = RWT(signal,7,'Morlet');

c = CWT_Wavelab(signal,5,'Morlet');
colors='20';
ImageCWT(c,'Individual',colors,'lin')