movielorez.m

lorenz吸引子的画法

% lorenz movie function
%this function can make the plot3 rotate
[t,y]=ode45('lorenz',[0 30],[12,9,7]);
axis([-20,20, -25, 25 ,10,50 ]);
view([20,42])
h=plot3(y(:,1),y(:,2),y(:,3)),grid on;
for j=1:100
        rotate(h,[0 1 0],1.8)
         M(j) = getframe;
end
movie(M,5)