chaos_parabola.m

来自「分形计算的几个程序原码」· M 代码 · 共 21 行

function chaos_parabola; 
% 抛物线-k*x^2+(k+1)*x混沌分岔行为 
% Author's email:[e-mail] 
% 
fff=inline('-k*x^2+(k+1)*x','x','k'); 
x0=0.2; 
kk=0:.01:3; 
xc=[];kc=[]; 
for k=kk; 
   x=x0; 
   for p=1:10; 
       x=fff(x,k); 
   end 
   for p=1:80; 
       x=fff(x,k); 
       xc=[xc,x]; 
       kc=[kc,k]; 
   end 
end 
plot(kc,xc,'k.','markersize',2)