plotpackettable.sci

小波分解源代码

function PlotPacketTable(pkt,scal)
// PlotPacketTable -- Display entries in wavelet, cosine packet tables
//  Usage
//    PlotPacketTable(pkt,scal)
//  Inputs
//    pkt     wavelet or cosine packet table
//    scal    optional scaling factor (0 => autoscale)
//
//  Side Effects
//    spike plot of coefficients in packet table
//
//  Copyright Aldo I Maalouf

  [lhs,rhs]=argn();
	[n,L] = size(pkt);
	D = L-1;
	if rhs < 2,
		scal = 0.;
	end

	LockAxes([0 1 (-D-1) (1)]); 

	t =( .5:(n-.5)) ./n;
	tt = zeros(1,3*n+2); 
	nxx = mtlb_length(tt);
	tt(1) = 0; tt(nxx)=1;
	tt(2:3:nxx-3) = t;
	tt(3:3:nxx-2) = t;
	tt(4:3:nxx-1) = t;
	xx = ones(size(tt));

	for d=0:D
		if(scal == 0),
			mult = .5 / (%eps + mtlb_max(abs(pkt(:,d+1))));
		else
			mult = scal;
		end

		xx = xx-1;
		xx(3:3:(nxx-2)) = -d + pkt(:,d+1) .* mult;
		plot(tt, xx); 

		nbox = 2^d;
		for b=0:(nbox-1)
			   mtlb_plot([ (b/nbox) (b+1)/nbox (b+1)/nbox (b/nbox)] , ...
			  [  -d-.5    -d-.5      -d+.5        -d+.5    ], ':');
		end
	end
	mtlb_hold off;//UnlockAxes;
	endfunction