clc;
clear all;
M=10000;N=10000;
D2=1;D3=0.45;D4=0;
L1=0;L2=0;
q=1;
for k=1:M;
x=zeros(1,N);
y=zeros(1,N);
x(1)=0;
y(1)=0;
for i=1:N-1;
x(i+1)=1-1.15.*x(i)^2+y(i);
y(i+1)=0.45*x(i);
end
if abs(x(end))<2;
D1=-2.3*x(end);
JT=[D1,D2;D3,D4];
[v,d]=eig(JT);
d=diag(d);
L1=L1+log(abs(d(1)));
L2=L2+log(abs(d(2)));
Xp(q)=x(end);
Yp(q)=y(end);
q=q+1;
end
end
L1=L1/(q-1);
L2=L2/(q-1);
评论0