求matlab大神编程！急用！ |
| 2013-06-14 51自学网 |
|
哥们，李雅普诺夫指数的算法有很多种，不知你是需要哪种算法的呢？
比如Nicolis方法、Benettin方法、Wolf方法、Jacobia方法等等
我把我以前计算伊侬映射李指数的程序给你，你参考一下：
%---------伊侬吸引子最大Laypunov指数的计算----------%
clear all;clc;
a=0.9:0.001:1.4;k=length(a);
b=0.3;p=600;
for n=1:k
for m=2:p
x(1,n)=0.4;y(1,n)=0.6;
x(m,n)=1+b*y(m-1,n)-a(n)*x(m-1,n)^2;
y(m,n)=x(m-1,n);
end
end
for r=1:k %计算雅克比矩阵
for h=2:p
A{1,r}=[-2*a(r)*x(1,r),b;1,0];
A{h,r}=[-2*a(r)*x(h,r),b;1,0]*A{h-1,r}; %注意元胞数组相乘顺序
end
end
for t=1:k %计算最大李指数
vv(:,t)=eig(A{p,t});v=max(abs(vv));
LE1=1/p*log(v);
end
plot(a,LE1,'k');hold on;
plot(a,0,'k:');
axis([a(1),a(k),-1 1]);
xlabel('a');ylabel('LE1');title('最大李指数'); |
|
|
|
