c_113.m

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

%混沌初值敏感性分析
N=10000;
x=0.2;y=0.2;z=0.2;q=1;
x0=0.2001;y0=0.2;z0=0.2;c=zeros(3,N);
for i=1:N
      x1=(-0.3-x+x^2-x*y+y+y^2)*q;
      y1=(0.2-x+x*y)*q;
      z1=(-x*y-x*z-z*y)*q;
      x=x1;y=y1;z=z1;    
    
      x2=(-0.3-x0+x0^2-x0*y0+y0+y0^2)*q;
      y2=(0.2-x0+x0*y0)*q;
      z2=(-x0*y0-x0*z0-z0*y0)*q;
      x0=x2;y0=y2;z0=z2;     
    
      c(1,i)=x2-x1;
      c(2,i)=y2-y1;
      c(3,i)=z2-z1;
end
for j=1:3
    subplot(1,3,j);
    plot(1:10:N,c(j,1:10:N));
end