cao_m.m

本文是通过文中方法来求最小嵌入维程序.子程序,相重构程序.

function [embed_m]=cao_m(data,min_m,max_m,tau)
%该程序计算时间序列的嵌入维数
% clc
% clear
% load stock
% data=x; % data为原始数据,列向量，n行1列
% min_m=1;
% max_m=10;% min_m,max_m分别为最小和最大嵌入维数
% tau=2;    % tau为时间延迟
% 作者:Adu,武汉大学,[e-mail]

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


disp('--------cao法求m---------')
e=0.1;embed_m=0;
h=1:length(E1)-1;
delt=abs(E1(h)-E1(h+1));
num=find(delt<e);
num1=find(delt==max(delt(num(1):length(delt))));
e=mean(delt(num1(1):length(E1)-1));
for kk=num1(1):length(E1)
    if kk+2<=length(E1)
        delt1=abs(E1(kk)-E1(kk-1));
        delt2=abs(E1(kk+1)-E1(kk));
        delt3=abs(E1(kk+2)-E1(kk+1));
        
        if (delt1>delt2)&(delt2>delt3)&(delt2<e)
            embed_m=kk;
            break
        end
    end
end