有关“离散求加速度”问题！

alian713

大家好，请问离散求加速度“咋求啊？我是用0de45求微分方程组，结果可输出位移、速度，但是在输出加速度的时候错误，估计求的时候错了！
   q_=[y(:,8),y(:,10),y(:,12)]; %各广义坐标的位移
    dq_=[y(:,9),y(:,11),y(:,13)]; %各广义坐标的速度
  D2q_=[diff(y(:,9)),diff(y(:,11)),diff(y(:,13))];[n,m]=size(D2q_);D2q_(n+1,:)=D2q_(n,:);%加速度！
这样能出结果，位移、速度图象还估计正确，但是加速度不对！请问那错了！多谢了啊！