wt08fig04.m

来自「老外写的小波变换的工具箱」· M 代码 · 共 41 行

% CAPTION :
fprintf('\n');
disp('Figure 8.4')
disp('Walsh wavelet packets at the depth j=3 of the Walsh tree,') 
disp('calculated with a Haar filter. They are ordered from')
disp('left to right and up to down, from low to high frequencies.')

close all;
n = 2048;
J = log2(n);
j = 7;
p = 10;
k = 2^(J-j-1);
figure(1);clf
for b = 0:7,
	w = MakeWaveletPacket(j,b,k,'Haar',p,n)*2^((j-3)/2);
	x = (abs(w)>1.e-6);
    x=double(x);
    length = x*x';
	i = (((1:n)-n/2)./2^(j-3)); 
	subplot(3,4,b+1)
	plot(i,w)
	axis([-1 9 min(w)-.05 max(w)+.05])
end;

% Written by Maureen Clerc and Jerome Kalifa, 1997
% [e-mail], [e-mail]

%Revision History
% 10/1/05        AM        logical vector x is converted to a double vector

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