function  tau=Mutual_Information_main(data)
%互信息法求tau
%data;     % 時間序列，列向量
max_t = 20;  % 本程序默認最大時延
%Part = 128;     % 本程序默認box大小

[entropy]=mutual(data,max_t);
for i = 1:length(entropy)-1           
    if (entropy(i)<=entropy(i+1))
        tau = i;            % 第一個局部極小值位置
        break;
    end
end
tau = tau -1 ;              
plot(0:length(entropy)-1,entropy)
xlabel('Lag');
title('互信息法求時延');