模拟 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
No comments:
Post a Comment