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