我对 Matlab 很陌生。我有一个脚本,在 Matlab 上使用ode45
并显示摆动弹簧在通过 3-D 空间时的运动。arrow.m
该程序几乎正在做我想做的事。现在,钻石的密度正在显示有效的弹簧速度(除了当ode45
将其作为个人最喜欢的数字样本),并且速度几乎可以通过函数的步长(至少以我的计算机运行代码的速度)准确计算。我想要做的是让我在代码中注释掉的位置矢量只显示在质量的瞬时位置,即曲线的端点,而不是钻石出现的每个点。我四处寻找帮助,但我尝试的一切似乎都会导致错误。如果有人能指出我正确的方向,那将不胜感激。请尝试运行程序,你会明白我的意思,也可以玩弄函数的参数。
function elasticPendulum(totalTime)
time=1;
totalTime=20
X=1
Vx=0
Y=0
Vy=0
Z=1.1
Vz=0
w=3;
g=9;
l=1;
while (time<totalTime)
tspan=[0,time];
x0=[X,Vx,Y,Vy,Z,Vz];
[t,x]=ode45(@F,tspan,x0);
plot3(x(:,1),x(:,3),x(:,5),'dk'), axis([-2 2 -2 2 -10 2]);
grid on;
o=[0, 0, 0];
Xeq=[0, 0, -(g/(w^2)+l)];
arrow(o,Xeq,'Length',5);
%{
Xt=[x(:,1) x(:,3) x(:,5)]
arrow(o,Xt,'Length',5);
%}
time=time+.1*(((x(2))^2+(x(4))^2+(x(6))^2)^(1/2))
pause(.1);
end
function xprime=F(t,x)
r=sqrt(x(1)^2+x(3)^2+x(5)^2);
xprime=[x(2);-w*(r-l)/r*x(1);x(4);-w*(r-l)/r*x(3);x(6);-w*(r-l)/r*x(5)-g];
end
end