function Lyapunov(x)
if length(x)<400
helpdlg('数据量需要大于400，请重新选定时间','提示信息');pause(3),
end
dm=26;tao=5;set_error=0.001;
for i=0:(dm-1)
for j=1:(length(x)/1.2)
k=j+i*tao;
Y(i+1,j)=x(k);
end
end
l=1; i=1; p=1; % i为Y的一列，j为Y的一列
while i<length(Y)
for j=p:length(Y) %求两点之间距离
if j>i
e_normal=max(abs(Y(:,i)-Y(:,j)));
end
end
flag=e_normal-set_error;
if flag>0
d(l)=e_normal;l=l+1;
%i=i+1;
for j=1:length(Y) %找离点最小距离的点
l_normal(j)=max(abs(Y(:,i)-Y(:,j)));
end
p=find((l_normal-min(l_normal))==0); %找到最近点Y(:,j)
else
j=j+1;i=i+1; %范数小于set_error，Y(:,i)和Y(:,j)继续向前一步
end
end
LE1=0;
for i=1:l
f1=d(i);f2=d(1);
LE1=log10(f1/f2)+LE1;
end
LE1=1/(l*0.014)*LE1,
计算logistic映射李雅普诺夫指数的程序
4星 · 超过85%的资源 需积分: 46 118 浏览量
2008-10-18
19:25:07
上传
评论 7
收藏 2KB RAR 举报 
logistic映射李雅普诺夫指数的matlab程序.rar （6个子文件） 
计算lorenz李雅普洛夫指数 
function Lyapunov(x).txt 1KB
logistic.m 7B
liyapunuofu.rar 307B liyapunuofu.m 326B logistic.asv 169B www.pudn.com.txt 218B
评论4