henon.m

来自「混沌与分岔」· M 代码 · 共 13 行

%x(n+1)=1-a*x(n)^2;y(n+1)=b*x(n)
N=1000;
x=0.3;y=0.5;
a=1.2;b=0.3;
for i=1:N
    x1=1-a*x^2+y;
    y1=b*x;
    plot(x1,y1);
    hold on;
    x=x1;
    y=y1;
end