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