模拟 Lorenz 系统的 Poincare 映像

Z=[]; [T,Y]=ode45('Lorenz3',[0,6000],[0.5;0.5;2;.9;.1;2;0]); for k=1:length(Y) if (abs(Y(k,2)-0.2)<1e-2) Z=[Z Y(k,1)+i*Y(k,3)]; end end plot(Z,'.','markersize',2) title('模拟 Lorenz 系统的 Poincare 映像 y=0.2') xlabel('x'),ylabel('z')

相关程序 Lorenz3.m