代码21-7.txt

matlab宝典书籍源码 全部代码 希望大家喜欢

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                                       % 马上画图
end