[发明专利]高精度地震模拟振动台的波形再现在线迭代方法

摘要

高精度地震模拟振动台的波形再现在线迭代方法，对迭代修正的传递函数矩阵进行在线更新提高振动台波形再现精度。在线迭代控制包括：输入原始加速度信号，对命令数据进行数据前加零值和添加白噪声处理后，求自相关函数和后修正傅里叶变换，得到自功率谱Gxx，命令信号傅式谱乘以传递函数矩阵的逆得到加速度信号傅式谱，然后反傅氏变换后得到输入加速度命令，振动台内部控制器执行加速度命令，收集振动台输出的加速度信号，求二者互相关函数，后修正傅里叶变换后得互功率谱Gxy，互功率谱除以自功率谱得到传递函数矩阵，通过判断条件采样步数是否大于傅里叶变换一帧的数据个数，是，更新传递函数矩阵；否，保持传递函数矩阵不变，从而完成在线迭代控制。

主权项

1.高精度地震模拟振动台的波形再现线迭代方法，主要由输入10原始加速度命令，将10原加速度命令前增加31‑傅里叶变换一帧的数据个数m的1.2倍个零值，得到11新的加速度命令，输入或生成白噪声加速度命令，将白噪声加速度命令与新的加速度命令对应数据相加得到12输入加速度命令，13传递12输入加速度命令acx、acy、acz后，求14自相关函数Rxx，再进行15后修快速傅里叶变换FFT得到其16自功率谱Gxx和其17加速度谱Acx、Acy、Acz，用18传递函数矩阵进行变换，得到19加速度信号Acx’、Acy’、Acz’ ，之后进行20反傅氏变换IFFT，得到21输入给振动台的加速度信号acx’、acy’、acz’，3输入加速度信号acx’、acy’、acz’输入给6振动台，由传感器收集到的9振动台输出的加速度信号amx、amy、amz，求22互相关函数Rxx，进行15后修正快速傅里叶变换FFT，得到23互功率谱Gxx，求24 Gxy/Gxx得到25传递函数矩阵FIO6x6，最后进行判断26 n ≥ m是否成立，27当n ≥ m成立时，更新传递函数矩阵FIO6x6= FIO6x6K+ FIO6x6 (1‑K)，28 当n ≥ m不成立时，缓存传递函数矩阵FIO6x6;29进行下一轮迭代，等步骤构成，从而使传递函数矩阵能够实时更新，保证传递函数矩阵实时性，从而提高了模拟地震波动再现精度。