由振动加速度求振动烈度的问题，急急急！

wind__rain

大家好，我用的是振动加速度传感器和数据采集卡采集振动信号，采集到的是时域离散的加速度值，采集卡的采样频率是10k，我采集了2048个点（a0,a1............a2047),请问我如何根据这些信息求得振动烈度值啊？我看到一个公式,这个算法和角频率有关，但我不知道角频率是如何求取的？

VibrationMaster

烈度应该是速度的方差，如果记录的是加速度信号，数值积分之后再求方差

wind__rain

能说的详细些么？拜托!比如如何积分?初速度v0如何选定?

VibrationMaster

1dat=dat-mean(dat)%去掉加速度直流分量.
2.初速度取0即可，最简单的是cumsum(dat)*dt;