toon0421.m

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

% toon0421 -- Compare Wavelet & Haar Partial Reconstructions
%
%  Try to reconstruct the Doppler signal using the best 50
%  coefficients from both the Haar and S8 wavelet transforms.
%  The S8 wavelet obviously does a better job of reconstruction.
%
	n       = 2048;
	Doppler = MakeSignal('Doppler',n);
	HaarQMF = MakeONFilter('Haar');
	S8QMF   = MakeONFilter('Symmlet',8);
	wcHaar  = FWT_PO(Doppler,1,HaarQMF);
	wcS8    = FWT_PO(Doppler,4,S8QMF);
	t       = (1:n)./n;
%
	BigwcHaar = KeepBiggest(wcHaar,50);
	subplot(221); PlotWaveCoeff(BigwcHaar,1,1);
	title('50 largest Haar[Doppler]')
%
	BigwcS8   = KeepBiggest(wcS8,50);
	subplot(222); PlotWaveCoeff(BigwcS8  ,4,1);
	title('50 largest in S8 Wavelet[Doppler]')
%
	subplot(223); PartialHaar = IWT_PO(BigwcHaar,1,HaarQMF);
	plot(t,PartialHaar)
	title('Haar Reconstruction from 50 largest')
	axis([ 0 1 -.5 .5])
%
	subplot(224); PartialS8 = IWT_PO(BigwcS8,4,S8QMF);
	plot(t,PartialS8)
	title('S8 Wavelet Reconstruction from 50 largest')
	axis([ 0 1 -.5 .5])
    
    
    
 
 
%
%  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]