gp_baozi9.m

来自「混沌的G-P算法」· M 代码 · 共 24 行

clc;
clear; 
data1 = load('e:\z.txt'); % z是我自己的实验数据
data=data1(1:200);
t=Mutual_Information_main(data) % 调用互信息法求时延

tic
R=[]; C=[];
r=[0.01:0.01:2];  %  尺度rexp([-1.99:0.01:1]);%[0.0002:0.0002:0.06];%

%%%%%%%%%
 for m=1:12; % 嵌入维
        Q=Phacespace(data,t,m);
        N=length(Q);
        D=zeros(N);
        [C,D]=Distance(Q,D,m,N,r);% 同时计算距离矩阵D和关联积分C
%%%%%%%%%%%%%%%%%%%%%

R=[log(r)'];

 plot(R,C)
end
toc