plotstatauto.m

％ Atomizer Main Directory, Version .802 里面信号含有分解去噪合成过程的代码 ％---------------------------------------

function PlotStatAuto(StatWT)
% PlotStatTable -- Plot Stationary Wavelet Transform, with each scale
%			scaled by itself
%  Usage
%    	PlotStatTable(StatWT)
%  Inputs
%    	StatWT   1-d stationary wavelet transform
%  Side Effects
%    	A depiction of the stationary wavelet transform,
%    	much like multi-resolution decomposition
%    	of signal
%  Description
%	For use in bpfig41 only
%  See Also
%	PlotCoeff
%

	[n,D1] = size(StatWT); 
	D = D1-1;
	J = log2(n);
	L = J-D;
%
	x = StatWT(:,1)';
	LockAxes([0 1 -(J) (-L+3)]);
	t = (.5:(n-.5))/n;
	scal = .5 ./ max(max(abs(x)));
	plot(t,(-(L-1)) + scal.*x);
%
	for j=(J-1):-1:L
	   x = StatWT(:,J-j+1);
	   scal = .5 ./ max(max(abs(x)));
	   plot(t, -(j) + scal.*x);
	end
%
%
	xlabel('t'); ylabel('Dyad');
	title('Stationary Wavelet Transform');
%
	UnlockAxes;