关于时域积分问题

无量石

我的是随机激励，采集时域响应加速度信号，采样频率25600，为什么积分出来的速度好像是一直衰减的，怎么没有振荡呢（我要做的是统计能量分析方向，只要求得振动能量就可以，所以只要积分得到速度就可以，但是这个结果感觉不太正确）谢谢指点！

load zhangfeiwentai1.tim;     %调用data中的数据，第一通道时间，第十通道激励力，其他加速度响应
x=zeros(2278784,9);
x1=zeros(2278784,9);
yvs=zeros(2278784,9);
for i=1:9
x(:,i)=zhangfeiwentai1(:,i+1);   %取通道中的数据 x(9)是力信号
end
%% 先只对第一通道积分
x1(:,1)=x(:,1)-mean(x(:,1));
sf=25600;%采样频率
t1=1/sf;    %采样时间间隔
%%%%%辛普森（simpson）算法时域积分求速度
yvs(1,1)=t1*(x1(1,1)+x1(2,1))/2;
n=length(zhangfeiwentai1(:,1));
for k=2:n-1;
yvs(k,1)=yvs(k-1,1)+t1*(x1(k-1,1)+4*x1(k,1)+x1(k+1,1))/6;
end
yvs(n,1)=yvs(n-1,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot (2,1,1);
title('实际信号输入信号')
plot(zhangfeiwentai1(:,1),x(:,1));
xlabel('时间(s)');
ylabel('加速度(g)');
subplot (2,1,2);
plot(zhangfeiwentai1(:,1),yvs(:,1));
xlabel('时间(s)');
ylabel('速度(m/s)');
grid on;

fanokay

楼主的问题解决了吗？
我用时域积分时也碰到了相似的问题，但是对同一组数据用频域积分算法，却能得到正常的速度曲线，这是为什么呢？理论上来说，算法应该是没问题的。。。
走过路过的高人，不吝赐教啊！3Q

无量石

还没解决，我就用这个速度接着往下做，最终的结果误差十几个数量级，只好从新做了，用频域积分怎么做？谢谢！