在你上面的评论中,当你说“这行不通”时,你是什么意思?根据文档的给定链接直接修改代码可以按预期工作:
a=[-8/3 0 0;0-10 10;0 28-1];
Y=[35-10-7]';
H=0.01;
P=动画线(Y(1),Y(2),Y(3),“marker”,“.”,“markersize”,2);
轴([0 50-25 25-25 25])
坚持
而1
a(1,3)=y(2);
A(3,1)=-Y(2);
YDOT= a*y;
Y= Y+H*YDOT;
加点(P,Y(1),Y(2),Y(3)
刷新屏幕
结束
< /代码>

按预期缠绕工作:
A = [ -8/3 0 0; 0 -10 10; 0 28 -1 ];
y = [35 -10 -7]';
h = 0.01;
p = animatedline(y(1),y(2),y(3),'Marker','.','MarkerSize',2);
axis([0 50 -25 25 -25 25])
hold on
while 1
A(1,3) = y(2);
A(3,1) = -y(2);
ydot = A*y;
y = y + h*ydot;
addpoints(p,y(1),y(2),y(3))
drawnow
end
