toon0443.m

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

% toon0443 -- Best-m Compression Numbers in DCT & DWT of Four Objects
%
%  Compare the accuracy of best-m reconstructions of four
%  signals with the Wavelet and DCT transforms.
%
	n = 2048;
	Bumps     = MakeSignal('Bumps',n);
	Doppler   = MakeSignal('Doppler',n);
	HeaviSine = MakeSignal('HeaviSine',n);
	Blocks    = MakeSignal('Blocks',n);
	t         = (1:n)./n;
%
	S8QMF       = MakeONFilter('Symmlet',8);
	wcBumps     = FWT_PO(Bumps,4,S8QMF);
	wcBlocks    = FWT_PO(Blocks,4,S8QMF);
	wcDoppler   = FWT_PO(Doppler,4,S8QMF);
	wcHeaviSine = FWT_PO(HeaviSine,4,S8QMF);
%
	dctBumps     = dct_iv(Bumps);
	dctBlocks    = dct_iv(Blocks);
	dctDoppler   = dct_iv(Doppler);
	dctHeaviSine = dct_iv(HeaviSine);
%
	subplot(221); semilogy(reverse(cumsum(sort(wcBlocks.^2))),'y--');
	hold on; semilogy(reverse(cumsum(sort(dctBlocks.^2))),'r-');
	title('Blocks; Best-m')
	xlabel('m')
	ylabel('d_m')
	axis([0 512 1 10^5])
	h = text(350,10^4,'DCT -'); set(h,'Color','Red');
	h = text(350,10^(3.3),'DWT --'); set(h,'Color','Yellow');
%
	subplot(222); semilogy(reverse(cumsum(sort(wcBumps.^2))),'y--');
	hold on; semilogy(reverse(cumsum(sort(dctBumps.^2))),'r-');
	title('Bumps; Best-m')
	xlabel('m')
	ylabel('d_m')
	axis([0 512 1 10^5])
	h = text(350,10^4,'DCT -'); set(h,'Color','Red');
	h = text(350,10^(3.3),'DWT --'); set(h,'Color','Yellow');
%
	subplot(223); semilogy(reverse(cumsum(sort(wcHeaviSine.^2))),'y--');
	hold on; semilogy(reverse(cumsum(sort(dctHeaviSine.^2))),'r-');
	title('HeaviSine; Best-m')
	xlabel('m')
	ylabel('d_m')
	axis([0 512 1 10^5])
	h = text(350,10^4,'DCT -'); set(h,'Color','Red');
	h = text(350,10^(3.3),'DWT --'); set(h,'Color','Yellow');
%
	subplot(224); semilogy(reverse(cumsum(sort(wcDoppler.^2))),'y--');
	hold on; semilogy(reverse(cumsum(sort(dctDoppler.^2))),'r-');
	title('Doppler; Best-m')
	xlabel('m')
	ylabel('d_m')
	axis([0 512 1 10^5])
	h = text(350,10^4,'DCT -'); set(h,'Color','Red');
	h = text(350,10^(3.3),'DWT --'); set(h,'Color','Yellow');
    
    
    
 
 
%
%  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]