wt04fig11.m

斯坦福大学的小波的工具箱 比matlab自带的工具箱更完善

%CAPTION
fprintf('\n');
disp('Figure 4.11')
disp('Window 1: Normalized scalogram of the signal shown in Figure 4.3,')
disp('calculated from the modulus of an analytic wavelet transform.')
disp('Dark points indicate large amplitude coefficients. ')
disp('Window 2: Complex phase of the analytic wavelet transform,')
disp('where the scalogram is non-zero.')

close all;
Name = 'Chirps';
N = 1024;
	f = MakeChirps(Name,N);
 	awt = AWT(f);
% axes handles	
	delta = 1/15;
	unit = (1-3*delta)/3;
	h3 = [1.5*delta delta 1-2.5*delta 2*unit];
	h2 = [delta 2*(unit+delta) 1-2*delta 2*unit];
	h1 = [delta 4*unit+3*delta 1-2*delta unit];
	h0 = [delta delta 1-2*delta 1-2*delta];
		sz     = size(awt);
	logawt = awt';
	linawt = log2linChap4(flipud(logawt),12);
	
figure(1);clf
        set(gcf, 'NumberTitle','off', 'Name','Window 1')
	axes('position',h3);
	ImageAWTnew(linawt','Overall','gray','lin');

figure(2);clf
        set(gcf,  'NumberTitle','off','Name','Window 2')
	axes('position',h3);
	phase = ourphase(linawt',1.e-1);
	pmax = max(max(phase));
	pmin = min(min(phase));
	phase = (phase-pmin)/(pmax-pmin);
	ImageAWTnew(phase,'Overall','gray','lin');


% Written by Maureen Clerc and Jerome Kalifa, 1997
% [e-mail], [e-mail]
    
    
%   
% Part of WaveLab Version 802
% Built Sunday, October 3, 1999 8:52:27 AM
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%