?? henon_lle.m
字號:
clear
close all
n=2; % 維數
%e=ones(n,1);
x(1)=0.6;
y(1)=0.4;
V=diag(ones(1,n));
S=V;
b=0;
k=300;
for i=1:k
J=[-2.8*x(i) 1;0.3 0];
B=J*V*S;
[V,S,U]=svd(B);
a_max=max(diag(S));
S=(1/a_max)*S;
b=b+log(a_max);
x(i+1)=1+y(i)-(1.4).*x(i).*x(i);
y(i+1)=0.3*x(i);
end
for i=1:n
Lyapunov(i)=(log(S(i,i))+b)/k;
fprintf('%f\n',Lyapunov(i));
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -