matlab 求加速度,【求助】位移转加速度(谱转换法)

【求助】位移转加速度(谱转换法)

那啥,小弟初学matlab,比较菜

有这样一组数据(样本见附件)

0.3987 1.3757 -5.0000 -0.8385

0.3984 1.3759 -5.0000 -0.8220

0.3989 1.3765 -5.0000 -0.7935

0.3986 1.3760 -5.0000 -0.7515

0.3984 1.3759 -5.0000 -0.7213

0.3976 1.3756 -5.0000 -0.7024

0.3973 1.3753 -5.0000 -0.6691

0.3978 1.3756 -5.0000 -0.6296

0.3969 1.3754 -5.0000 -0.5919

... .... ... ...

说明:

1.采样速率20K，k=1000

2.文件长度22K

3.第1，2 列激光位移，第3列是干扰,(忽略)第4列加速度

现在要求用频谱转换法,将加速度转为位移

算法是这个样子的:

如果在时间T内采集N个数据,

x(n)经离散傅里叶变换后得到的X(k)是一个长度为N的复数序列(频谱),

X(k)=DFT[x(n)]=[(a0，jb0)，(a1，jb1 )，...,(an-1 ，jbn-1 )]

x(n)中各谐波分量的幅值、圆频率及初相角可由式(1)~(3)求出

Ak=sqrt(ak^2