6-24.m

来自「《MATLAB 7.0编程基础》第6章 （源码实例）主要讲解MATLAB图形基础」· M 代码 · 共 20 行

A = [ -8/3 0 0; 0 -10 10; 0 28 -1 ];
y = [35 -10 -7]';
h = 0.01;
p = plot3(y(1),y(2),y(3),'.', ...
    'EraseMode','none','MarkerSize',5); 
% 设置擦除模式为none
axis([0 50 -25 25 -25 25])
hold on
for i=1:4000
    A(1,3) = y(2);
    A(3,1) = -y(2);
    ydot = A*y;
    y = y + h*ydot;
    % 改变坐标轴
    set(p,'XData',y(1),'YData',y(2),'ZData',y(3)) 
    drawnow 
    % 完成绘制事件
    i=i+1;
end