matlab加速度积分计算速度和位移

1、内容简介
略

2、内容说明
略

3、仿真分析clc
close all
clear
load('data.txt')
% 起始时间：2017年10月22日16:55:21
% 时间(s)    x轴加速度    y轴加速度    z轴加速度    温度
t=data(:,1);
acc_x=data(:,2);
figure
plot(t,acc_x)
xlabel 时间/s
ylabel 加速度
f = 100;
% t_new = t(1):t_step:t(end);
% acc_x_new = interp1(t,acc_x,t_new);
vel_x = zeros(length(t),1);
dis_x = zeros(length(t),1);
for i = 2:length(t)
    vel_x(i) = (acc_x(i-1)+acc_x(i))*0.5/f;
end
figure
plot(t,vel_x)
xlabel 时间/s
ylabel 速度
for i = 2:length(t)
    dis_x(i) = (vel_x(i-1)+vel_x(i))*0.5/f;
end
figure
plot(t,dis_x)
xlabel 时间/s
ylabel 位移

4、参考论文
略