振动台试验数据处理指点

振动台试验的数据量大，常规的数据处理的办法可能效率极低，拟采用Matlab编程的方法解决此问题，作为规划，也为了备忘，整理需完成的内容如下：1.起始零点的调整；按照台面加速度时程最大值与原始波形最大值重合的方法，已经完成。2.峰值调整；试验时最大加速度峰值可能与目标值略有偏差，按线性方法调整，已经完成。3.时程曲线基线调整考虑到部分时程曲线有趋势项，进行去趋势项操作，可挑战拟合阶数。已经完成。4.滤波部分时程曲线含有高频噪声，采用带通滤波的方法滤波。结合积分的滤波已经完成，单独滤波部分未完成。5.时程曲线积分对加速度时程曲线进行两次积分，采用频域积分的办法。已经完成。6.各通道之间的运算考虑进行各通道之间的运算（如求层间位移），未完成。7.应变数据时标的调整由于应变数据的采集与加速度时程的时标不一致，考虑一种合适的方法进行应变数据的时标调整。未完成。8.位移数据的时标调整与应变数据一样，考虑调整位移数据的时标调整。未完成。9.GUI界面考虑到易用性，采用MATLAB的GUI编程，降低使用难度。未完成。振动台试验非常关心结构的动力特性及其变化,常用的测试方法有:1.白噪声扫描;2.敲击激振法;经验证,采用地震波激振的方法同样可以得出其自振频率的分布,利用求输入信号(通常为台面加速度信号)与输出信号(如顶层加速度)之间的传递函数.MATLAB中可采用语句:[Txy,F]=tfestimate(x(:,3),x(:,12),1024,[],[],128);其中:x(:,3),x(:,12)分别为输入和输出信号列,128为采样频率.加速度信号由于测试方便,在振动台试验中有广泛的应用,加速度信号进行两次积分后就是位移的信号,但由于常数的存在,往往发生积分结果的漂移.采用加速度信号的频域积分可以解决此问题，根据王济的《MATLAB在振动信号处理中的应用》中的程序，编程后运行效果良好。按照王济的《MATLAB在振动信号处理中的应用》一书的原程序，我运行验证是可行的。运行结果（采样频率128Hz，多项式阶数：3）：