plot3c.m

一种新的时频分析方法的matlab源程序。

function  nt= plot3c(t,c,i)

% The function PLOT3C plots 3 subplots starting from the I-th component
% of the data c(n,k), where n is the number of data points and
% k is the number of IMF components.
% The number of subplots in vertical direction is equal to 3.
% The image on the screen closes all existing figures.
%
% Calling sequence-
% nt=plot3c(t,c,i)
%
% Input-
%	t		    - vector t(n) representing the 
%				time scale values
%	c		    - component data array c(n,k)
%	i		    - index of the first plotted component 
% Output-
%	nt		    - the handle of the figure
 
% Z. Shen (Caltech)		Initial

close;

tn=length(t);
t1=min(t(1),t(tn));
t2=max(t(1),t(tn));
for j=i:i+2
   subplot(3,1,j+1-i),plot(t,c(:,j));
   axis([t1 t2 min(c(:,j)) max(c(:,j))]);

%   if j==1 
%     s2='X';
%   elseif j==2
%     s2='Y';
%   else
%     s2='Z';
%   end
	
   s1='c';
   s2=[s1 num2str(j)];
   ylabel(s2);
end
nt=1;