mutual_information_main.m

来自「互信息法求混沌序列的最优时延 文件说明： 1、Mutual_Inform」· M 代码 · 共 36 行

% MEX 文件调用 Matlab 函数

clc
clear
close all
disp('------------- Input --------------');

sigma = 10;       % Lorenz 方程参数 a
b = 8/3;          %                 b
r = 34;           %                 c            

y = [-1,0,1];    % 起始点 (1 x 3 的行向量)
h = 0.01;        % 积分时间步长

k1 = 6000;       % 前面的迭代点数
k2 = 3000;       % 后面的迭代点数

z = LorenzData(y,h,k1+k2,sigma,r,b);
z = z(k1+1:end,:);

X = z(:,1);     % 时间序列，列向量
maxLags = 100;  % 最大时延
Part = 128;     % 
r = Mutual_Information(X,maxLags,Part);

for i = 1:length(r)-1           
    if (r(i)<=r(i+1))
        tau = i;            % 第一个局部极小值位置
        break;
    end
end
optimal_tau = tau -1               % r 的第一个值对应 tau = 0,所以要减 1

plot(0:length(r)-1,r)
xlabel('Lag');
title('互信息法求时延');