alian713 发表于 2007-1-10 18:49

大家好，求加速度问题！！着急！

大家好啊！我用ode45求微分方程，想叫画出加速度的曲线，不知道咋办？我用下面的办法输出的结果位移、速度都对， 但是加速度不对，不知道咋办！麻烦饿！
q_=; %各广义坐标的位移
dq_=; %各广义坐标的速度
D2q_=;=size(D2q_);D2q_(n+1,:)=D2q_(n,:);%加速度
figure;plot(t,dq_(:,1));
figure;plot(t,dq_(:,2));
figure;plot(t,dq_(:,3));
   
figure; plot(t,D2q_(:,1));
figure; plot(t,D2q_(:,2));
figure; plot(t,D2q_(:,3));

skywm 发表于 2007-1-10 21:48

可能是由于用ODE解微分方程,步长总是在随时调整,建议用从位移解速度的方法,从速度求加速度,而不是直接用diff().

alian713 发表于 2007-1-10 22:04

谢谢，问题已经解决
页: [1]
查看完整版本: 大家好，求加速度问题！！着急！