对加速度二次积分求位移

n = 10 * 1024;                          % cai yang dian shu
fs = 1000;                              % cai yang ping lv
t = (0 : 1 : n) / fs;                   % shi jian
fre = 2;                                % xin hao ping lv
jiasudu = 1000 * sin(2 * pi * fre * t); % jia su du
sudu = cumtrapz(t, jiasudu);            % su du
sudu = sudu - mean(sudu);
weiyi = cumtrapz(t, sudu);              % wei yi

subplot(3, 1, 1)
plot(t, jiasudu);
subplot(3, 1, 2)
plot(t, sudu);
subplot(3, 1, 3)
plot(t, weiyi);

加速度一次积分求速度 / 2 * pi * f

速度一次积分求位移 / 2 * pi * f