数值积分

摘要：本次主要介绍利用MATLAB对加速度信号进行两次数值积分到位移的算法。

1. 假如加速度信号是一个正弦信号。

A = 1;
f = 40;
w = 2*pi*f;
t = 0:0.001:0.4;
fs = 1000;
acc = A*sin(w*t);

2. 对加速度信号进行两次积分。

spd = cumtrapz(t,acc);
spd = spd * 1000;
pos = cumtrapz(t,spd);
pos = pos * 1000;

3. 画出加速度、速度、位移与时间的曲线图。

可以看到，位移信号存在类似递增的趋势。

4) 用高通滤波对位移信号去趋势。去趋势后的结果如下：

可以看到，前两个周期内的位移信号存在一定程度的信号失真。但是后面的积分后的位移信号是可以接受的结果。