donghua2.m

来自「包括MATLAB7.0实用指南（1）书中所有的例子的源代码。」· M 代码 · 共 27 行

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);
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
hold on
p=plot3(y(1),y(2),y(3),'square','EraseMode','background','MarkerSize',15,'MarkerEdgeColor',[1 .7 .7],'MarkerFaceColor',[0 0 0]);
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
hold off