% Thursday 10/7
 %------------Example 2 -------------------------------------
 %       Wavelet Piacket Analysis on Linear Chirp Signal

        n = 1024;
        t = (0:(n-1))./n;
        qmf = MakeONFilter('Daubechies',4);
        D   = log2(n);

        xi              = pi/3 .* t;
        LinChirp  = cos(xi .* (1:n));

        figure;
        plot(t, LinChirp)
        title('LinChirp');
        print(gcf,'-deps', 'hand20.eps')

%       Illustrate WP on Linear Chirp

        x   = LinChirp;
        figure;
        titlestr = 'LinChirp';

        wp  = WPAnalysis(x,D,qmf);
        PlotPacketTable(wp);
        signaltitle = [ 'Wavelet Packet Decomposition of ' titlestr];
        title(signaltitle); ylabel('Split Level'); xlabel('Frequency[Time]')
        print(gcf,'-deps', 'hand21.eps')

        stree = CalcStatTree(wp,'Entropy');
        [btree,vtree] = BestBasis(stree,D);
        signaltitle = [' WP Best Basis; ' titlestr];
        PlotBasisTree(btree,D,stree,signaltitle);
        xlabel('Frequency Domain Splits')
        ylabel('Entropy Gain')
        print(gcf,'-deps', 'hand22.eps')

        ImagePhasePlane('WP',btree,wp,signaltitle,256,qmf)
        print(gcf,'-deps', 'hand23.eps')