对位移二次微分求加速度

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
weiyi = 10 * sin(2 * pi * fre * t);     % wei yi

sudu = diff(weiyi) * fs;                % su du
tsudu = (0 : 1 : n - 1) / fs;
jiasudu = diff(sudu) * fs;              % jia su du
tjiasudu = (0 : 1 : n - 2) / fs;

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

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

速度一次微分求加速度 * 2 * pi * f