plotnc.m

这是一个关于hht变换很有用的工具箱

function  nt= plotnc(t,c,start,stop,fig_num)
% plotnc(t,c,start,stop,fig_num,stitle,xunits)
%		t is time axis
%		c is component array
%		plot begins from the component 'start'
%		plot ends at component stop
%		plot appears in Figure fig_num
%
% --- Steven R. Long at NASA GSFC / WFF --- Version 03/25/99 -- %
% --- Based on an Older Approach by Z. Shen ------------------- %

figure(fig_num);
clf;
tn=length(t);
t1=min(t(1),t(tn));
t2=max(t(1),t(tn));
for j=start:stop;
   subplot(stop-start+1,1,j+1-start),plot(t,c(:,j));
   axis([t1 t2 min(c(:,j)) max(c(:,j))]);	
   s1='c';
   s2=[s1 num2str(j)];
   ylabel(s2);
   if(j < stop);
      set(gca,'xticklabel','');
   end;
end;