加速度から変位の変換について

下記のようにMatlabにて積分して加速度時刻歴波形を速度、変位に変換しようとしています。
実測も変位データと照らし合わせてもに全然合っていなく、変位になっていないので困っています。
間違っているのでしょうか？完全に一致するとは思いませんが。。。。。
どなたかご教授下されば幸いです。


%t:時間[sec]
%acc：加速度波形[Gal]
acc=wave(7,:);
vel=cumtrapz(acc,t);
dis=cumtrapz(vel,t);
figure
subplot(3,1,1)
　plot(t,acc,'r'),grid on
　xlabel('time[sec]'),ylabel('ACC[cm/s~2]')
subplot(3,1,2)
　plot(t,vel,'r'),grid on
　xlabel('time[sec]'),ylabel('VEL[cm/s]')
subplot(3,1,3)
　plot(t,dis,'r'),grid on
　xlabel('time[sec]'),ylabel('DISP[cm]')

No.1 回答者： m0r1_2006 回答日時： 2013/07/06 14:41

cumtrapz を用いているので，

速度も位置も初期値 0 で計算しているので，あうわけがない．

変位データがあるなら，
速度の初期値を
v0 = (変位(t=dt) - 変位(t=0))/dt;
とかで計算して，
vel = v0 + vel;

dis = dis + 変位(t=0);

などにすれば，少しはましになるかも