function?lambda_1=lyapunov_wolf（dataNmtauP）
%?該函數用來計算時間序列的最大Lyapunov?指數--Wolf?方法
%?m:?嵌入維數
%?tau:時間延遲
%?data:時間序列
%?N:時間序列長度
%?:時間序列的平均周期選擇演化相點距當前點的位置差，即若當前相點為I，則演化相點只能在|I－J|>的相點中搜尋
%?lambda_1:返回最大lyapunov指數值
min_point=1?;?%&&要求最少搜索到的點數
MAX_CISHU=5?;?%&&最大增加搜索范圍次數
%FLYINGHAWK
delt_t1?=?0.001;
t1?=?0:delt_t1:60;
[tt1y1]=ode45（‘D‘t1[-1010]）;
xx1?=?y1（:1）‘;
x1?=?spline（tt1?xx1?t1）;
data=?x1（20000:10:60000）;%采樣
N=length（data）;
m=4;
tau=11;
%*****************************************************
%?FFT計算平均周期
%**********************************************************
x=data;
xPower=abs（fft（x））.^2;
NN=length（xPower）;
xPower（1）=[];%去除直流分量
NN=floor（NN/2）;
xPower=xPower（1:NN）;
freq=（1:NN）/NN*0.5;
[mPindex]=max（xPower）;
P=index;
%?求最大、最小和平均相點距離
max_d?=?0;?%最大相點距離
min_d?=?1.0e+100;?%最小相點距離
avg_dd?=?0;
Y=reconstitution（dataNmtau）;?%相空間重構
M=N-（m-1）*tau;?%重構相空間中相點的個數
for?i?=?1?:?（M-1）
for?j?=?i+1?:?M
d?=?0;
for?k?=?1?:?m
d?=?d?+?（Y（ki）-Y（kj））*（Y（ki）-Y（kj））;
end
d?=?sqrt（d）;
if?max_d?max_d?=?d;
end
if?min_d?>?d
min_d?=?d;
end
avg_dd?=?avg_dd?+?d;
end
end
avg_d?=?2*avg_dd/（M*（M-1））;?%平均相點距離

dlt_eps?=?（avg_d?-?min_d）?*?0.02?;?%若在min_eps～max_eps中找不到演化相點時，對max_eps的放寬幅度
min_eps?=?min_d?+?dlt_eps?/?2?;?%演化相點與當前相點距離的最小限
max_eps?=?min_d?+?2?*?dlt_eps?;?%&&演化相點與當前相點距離的最大限

%?從P+1～M-1個相點中找與第一個相點最近的相點位置（Loc_DK）及其最短距離DK
DK?=?1.0e+100;?%第i個相點到其最近距離點的距離
Loc_DK?=?2;?%第i個相點對應的最近距離點的下標
for?i?=?（P+1）:（M-1）?%限制短暫分離，從點P+1開始搜索
d?=?0;
for?k?=?1?:?m
d?=?d?+?（Y（ki）-Y（k1））*（Y（ki）-Y（k1））;
end
d?=?sqrt（d）;
if?（d??min_eps）
DK?=?d;
Loc_DK?=?i;
end
end
%?以下計算各相點對應的李氏數保存到lmd（）數組中
%?i?為相點序號，從1到（M-1），也是i-1點的演化點；Loc_DK為相