demo.m

来自「时间序列工具箱」· M 代码 · 共 52 行

% DEMO program of the TSA-toolbox

%	Version 2.44
%	last revision 22.06.1998
%	Copyright (c) 1997-1998 by Alois Schloegl
%	e-mail: [e-mail]	

% This library is free software; you can redistribute it and/or
% modify it under the terms of the GNU Library General Public
% License as published by the Free Software Foundation; either
% Version 2 of the License, or (at your option) any later version.
%
% This library is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
% Library General Public License for more details.
%
% You should have received a copy of the GNU Library General Public
% License along with this library; if not, write to the
% Free Software Foundation, Inc., 59 Temple Place - Suite 330,
% Boston, MA  02111-1307, USA.

K=1;
while K
	K=menu('Select a demo ',...
		'Display demo Time series',...
	'1st investigation of a TS (invest1)',...
	'Compare performance',...
	'Inverse filtering',...
	'Bispectral analysis',...
	'Histogram and entropy of TS',...
	  'end');
subplot(111);
clear eeg8s;
if K==1	
	load eeg8s.mat;
	plot((1:length(eeg8s))/128,eeg8s);
	title('Demo time series: EEG sampled with 128Hz for 8 seconds');
elseif K==2	tsademo;
elseif K==3	demoperf;
elseif K==4	invfdemo;
elseif K==5	bisdemo;
elseif K==6	
	load eeg8s.mat;
	[H,X,E]=histo(eeg8s);
 	bar(X,H);
	fprintf(1,'The entropy of information is E=%1.3f\n',E);
elseif K==7
	K=0;
end;
end; %WHILE