baker.m

来自「baker（面包师变换）的混沌图像」· M 代码 · 共 15 行

x(1)=0.1,y(1)=0.0001;a=1/4;
for n=1:100
    if(x(n)>=0)&(x(n)<=1/2)
      x(n+1)=2*x(n);
    elseif(x(n)>1/2)&(x(n)<=1)
      x(n+1)=2*x(n)-1;
    end;
   if(y(n)>=0)&(y(n)<=1/2)
      y(n+1)=a*y(n);
   elseif(y(n)>1/2)&(y(n)<=1)
      y(n+1)=a*x(n)+1/2;
   end;
      plot(x,y);
       xlabel('x     x(1)=0.1,y(1)=0.0001;a=1/4'),ylabel('y'),title('baker');
end;