mpwork21.m

老外写的小波变换的工具箱

% MPWork21 -- Matching-Pursuit, Wavelet Packet Tour on Four Basic Signals
%
%     	HeaviSine
% 		Doppler
%		QuadChirp
%		MishMash
%
	n   = 1024;
	qmf = MakeONFilter('Coiflet',3);
%
	compressfig  = figure;
	equivratefig = figure;
	equivdimfig  = figure;
	figarray = zeros(4,1);
%
	for i=1:4,
		if i==1,
			SigName = 'LinChirp';
		elseif i==2,
			SigName = 'TwoChirp';
		elseif i==3,
			SigName = 'QuadChirp';
		elseif i==4,
			SigName = 'MishMash';
		end			
%
		hs = MakeSignal(SigName,n);
		figarray(i) = figure;
		[mpatoms,mpresid] = WPPursuitTour('P',hs,10,qmf,300,SigName);
		subplot(2,2,3)
		axis([0 500 0 10])
	
		energy = norm(hs).^2;
		mpcnums = [energy ; energy - cumsum(mpatoms(:,1).^2)];

		figure(compressfig);
		subplot(2,2,i)
		PlotSemiCompressNum(mpcnums,'MPWP', SigName);

		figure(equivratefig);
		subplot(2,2,i)
		PlotEquivRate(mpcnums,'WP', SigName)

		figure(equivdimfig);
		subplot(2,2,i)
		PlotEquivDimension(mpcnums,'WP', SigName)

	end
%
	%ifprint(compressfig,'mp_fig215.ps')
	%ifprint(equivratefig,'mp_fig216.ps')
	%ifprint(equivdimfig,'mp_fig227.ps')
	
	for i=1:4,
		ifprint(figarray(i),['mp_fig21' int2str(i) '.ps']);
	end

			
    
    
 
 
%
%  Part of Wavelab Version 850
%  Built Tue Jan  3 13:20:43 EST 2006
%  This is Copyrighted Material
%  For Copying permissions see COPYING.m
%  Comments? e-mail [e-mail]