% Thursday 10/7
 %------------Example 3 -------------------------------------

 %           This is famous MishMash =====>
 %
 %           input: n  
 %           t = (1:n) ./n;
 %           MishMash = QuadChirp + LinChirp + Sine
 %           sig = sin( (pi/3) .* t .* (n .* t.^2)) ; %<- QuadChirp
 %           sig = sig +  sin( pi * (n * .6902) .* t); %<- Sine
 %           sig = sig +  sin(pi .* t .* (n .* .125 .* t)); <- LinChirp
 %
        figure(1);
                set(0, 'DefaultAxesFontSize', 12);
                set(0, 'DefaultAxesFontName', 'Times')
        clf;
        n=1024; t=(0:(n-1))./n;
        depth=floor(log2(n));
        mm = MakeSignal('MishMash',n);
        plot(t, mm)
        xlabel('t'); ylabel('MishMash signal')
        print(gcf, '-deps', 'hand30.eps')
 %      -------------------------------
        figure(2);
                set(0, 'DefaultAxesFontSize', 12);
                set(0, 'DefaultAxesFontName', 'Times')
        clf;
        qmf=MakeONFilter('Daubechies',4);
        WPTour('P',mm,depth,qmf,'MishMash');
 %   
 %      WPTour: plot of Wavelet packet Table
 %       WPTour(Format,x,D,qmf,title)
 %      Format: 'P' for four-panel display ('S' for sequential display)
 %      x: 1-d signal to analyze; D-fineness of frequency splitting
 %      qmf: quadrature mirror filter; title-name of signal
 %
 %      tree plot, best basis
 %      phase plot, best basis
 %      printout, entropy, best basis

        subplot(2,2,3)
        axis([0 500 -10 10])
        subplot(2,2,4)
        brighten(-.8);
        print(gcf, '-deps', 'hand31.eps')