如何对采集到的信号做两次积分

[rose007]

各位大侠：
我是用采集到的加速度信号来做两次积分得到位移信号，由于采集到的加速度信号有直流偏移量，所以用了detrend命令做去趋势处理。为了检验代码是否正确，我用了一个标准的正弦函数来验证。代码如下:
t=(0:1:10240-1)/2000;
a=sin(5*t); %标准的正弦函数，准备对其进行两次积分处理
fs=2000;
N=10240;
za=detrend(a); %对采集到的加速度信号去趋势处理(因为采集到的加速度信号有直流偏移量)
v=cumtrapz(t,za); %对加速度信号积分得速度信号
zv=detrend(v); %对速度信号去趋势处理
l=cumtrapz(t,zv); %对速度信号积分得到位移信号
zl=detrend(l); %对位移信号去趋势处理
figure(1)
subplot(3,1,1)
plot(t,za)
title('加速度信号');
subplot(3,1,2)
plot(t,zv)
title('速度信号');
subplot(3,1,3)
plot(t,zl)
title('位移信号');



这个图算出来的正弦信号两次积分不对啊，这个正弦两次积分也应该是标准的正弦啊，请各位大侠帮我看看是哪的问题
先谢谢了。