imagephasevarious.m

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

function ImagePhaseVarious(wp,wp_stree,cp,cp_stree,titlestr)
% ImagePhaseVarious -- Compare four phase plane displays
%  Usage
%    ImagePhaseVarious(wp,wp_stree,cp,cp_stree,titlestr)
%  Inputs  
%    wp          Wavelet Packet table
%    wp_stree    Wavelet Packet Stat tree
%    cp          Cosine Packet table
%    cp_stree    Cosine Packet Stat tree
%    titlestr    title for plots
%
%  Side Effects
%    For Both Wavelet packets & Cosine Packets:
%       phase plot, best level
%       phase plot, best basis
%       printout, entropies best level, best basis
%
%  See Also
%    ImagePhasePlane, PlotPhaseVarious
%
	[n,L] = size(wp); D=L-1;
	disp(['Best Level Search in Wavelet Packet Table for ' titlestr]) 
	bestwplevel = CompareStdBases(wp_stree,D);
	header = sprintf('WP Best Level(=%g) ',bestwplevel);
	signaltitle = [header titlestr];
	subplot(221);
	ImageGaborPhase('WP',wp,bestwplevel,signaltitle);
%
	subplot(222);
	bwptree = BestBasis(wp_stree,D);
	signaltitle = ['WP Best Basis; ' titlestr];
	ImagePhasePlane('WP',bwptree,wp,signaltitle);
%
	[n,L] = size(cp); D=L-1;
	disp(['Best Level Search in Cosine Packet Table for ' titlestr]) 
	bestcplevel = CompareStdBases(cp_stree,D);
	header = sprintf('CP Best Level(=%g) ',bestcplevel);
	signaltitle = [header titlestr];
	subplot(223);
	ImageGaborPhase('CP',cp,bestcplevel,signaltitle);
%
	subplot(224);
	signaltitle = ['CP Best Basis; ' titlestr];
	bcptree = BestBasis(cp_stree,D);
	ImagePhasePlane('CP',bcptree,cp,signaltitle);

%
% Copyright (c) 1993. David L. Donoho
%     
    
    
 
 
%
%  Part of Wavelab Version 850
%  Built Tue Jan  3 13:20:40 EST 2006
%  This is Copyrighted Material
%  For Copying permissions see COPYING.m
%  Comments? e-mail [e-mail]