% 互信息法求時延% 使用平臺 - Matlab6.5 / Matlab7.0% 作者：陸振波，海軍工程大學(xué)% 歡迎同行來信交流與合作，更多文章與程序下載請訪問我的個人主頁% 電子郵件：luzhenbo@yahoo.com.cn% 個人主頁：http://luzhenbo.88uu.com.cnclcclearclose all%----------------------------------------------------------------------%% 產(chǎn)生混沌時間序列sigma = 10;             % Lorenz 方程參數(shù) ab = 8/3;                %                 br = 34;                 %                 c            y = [-1,0,1];           % 起始點(diǎn) (1 x 3 的行向量)h = 0.01;               % 積分時間步長k1 = 6000;              % 前面的迭代點(diǎn)數(shù)k2 = 3000;              % 后面的迭代點(diǎn)數(shù)z = LorenzData(y,h,k1+k2,sigma,r,b);z = z(k1+1:end,:);%----------------------------------------------------------------------%% 調(diào)用 mex 函數(shù)X = z(:,1);             % 時間序列，列向量maxLags = 100;          % 最大時延Part = 128;             % 每一座標(biāo)劃分的份數(shù)r = Mutual_Information(X,maxLags,Part);%----------------------------------------------------------------------%% 尋找第一個局部極小點(diǎn)for i = 1:length(r)-1               if (r(i)<=r(i+1))        tau = i;        % 第一個局部極小值位置        break;    endendoptimal_tau = tau -1    % r 的第一個值對應(yīng) tau = 0,所以要減 1%----------------------------------------------------------------------%% 圖形顯示plot(0:length(r)-1,r)xlabel('Lag');title('互信息法求時延');