%?互信息法求時延
%?使用平臺?-?Matlab6.5?/?Matlab7.0
%?作者：陸振波，海軍工程大學
%?歡迎同行來信交流與合作，更多文章與程序下載請訪問我的個人主頁
%?電子郵件：luzhenbo@yahoo.com.cn
%?個人主頁：http://luzhenbo.88uu.com.cn

clc
clear
close?all

%----------------------------------------------------------------------%
%?產(chǎn)生混沌時間序列

sigma?=?10;?????????????%?Lorenz?方程參數(shù)?a
b?=?8/3;????????????????%?????????????????b
r?=?34;?????????????????%?????????????????c????????????

y?=?[-101];???????????%?起始點?（1?x?3?的行向量）
h?=?0.01;???????????????%?積分時間步長

k1?=?6000;??????????????%?前面的迭代點數(shù)
k2?=?3000;??????????????%?后面的迭代點數(shù)

z?=?LorenzData（yhk1+k2sigmarb）;
z?=?z（k1+1:end:）;

%----------------------------------------------------------------------%
%?調(diào)用?mex?函數(shù)

X?=?z（:1）;?????????????%?時間序列，列向量
maxLags?=?100;??????????%?最大時延
Part?=?128;?????????????%?每一座標劃分的份數(shù)
r?=?Mutual_Information（XmaxLagsPart）;

%----------------------------------------------------------------------%
%?尋找第一個局部極小點

for?i?=?1:length（r）-1???????????
????if?（r（i）<=r（i+1））
????????tau?=?i;????????%?第一個局部極小值位置
????????break;
????end
end
optimal_tau?=?tau?-1????%?r?的第一個值對應(yīng)?tau?=?0所以要減?1

%----------------------------------------------------------------------%
%?圖形顯示

plot（0:length（r）-1r）
xlabel（‘Lag‘）;
title（‘互信息法求時延‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1397??2007-03-15?15:13??Main_Mutual_Information\Main_Mutual_Information.m
?????文件????????1618??2019-04-17?13:35??Main_Mutual_Information\mutual_information.m
?????目錄???????????0??2019-04-17?14:03??Main_Mutual_Information\