?? henon_lle_nicolis2.m
字號:
% clear;
close all;
k=300;
for i=1:301
a=(i-1)*5/1000;
x(1)=0.6;
y(1)=0.4;
dx=[0.002;0.01];
d0=norm(dx);
e=1/d0*dx;
le=0;
for j=1:k
J=[-2*a*x(j) 1;0.3 0];
vec=J*e;
d=norm(vec);
e=1/d*vec;
le=le+1/k*log(d);
x(j+1)=1+y(j)-a.*x(j).*x(j);
y(j+1)=0.3*x(j);
end
lle(i)=le;
end
a=0:0.005:1.5;
plot(a,zeros(1,301),'k');
hold on
plot(a,lle,'b');
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -