pangjialai.m

混沌映射庞加莱界面Matlab程序

clear
ff=0.222;
options=odeset('RelTol',1e-7);
tt=2*pi/1.2
[t,x]=ode45(@dafen,[0:tt/100:80*tt],[0,0],options,[],ff);
plot(x(2000:end,1),x(2000:end,2),'-')
%pojialai
hold on
i=2000:100:3000
plot(x(i,1),x(i,2),'*')