function  cao_m(data,min_m,max_m,tau)
%該程序用Cao氏法計(jì)算時(shí)間序列的嵌入維數(shù)
% clc
% clear
% load stock
% data=x; % data為原始數(shù)據(jù),列向量，n行1列
  min_m=1;
% max_m=10;% min_m,max_m分別為最小和最大嵌入維數(shù)
% tau=2;    % tau為時(shí)間延遲
% 作者:Adu,武漢大學(xué),adupopo@163.com

[E1,E2]=cao1(data,min_m,max_m,tau);
n=length(E1);
plot(1:n,E1,'-bs',1:n,E2,'-r*');xlabel('維數(shù)');ylabel('E1&E2');
title('Cao氏法求最小嵌入維數(shù)');
legend('E1','E2');
grid on