[发明专利]一种广义高阶同步挤压地震信号时频分解与重构方法

摘要

本发明公开了一种广义高阶同步挤压地震信号时频分解与重构方法，包括输入一维信地震号；选择阶数和时窗函数，计算在窗函数下的短时傅里叶变换值，设置阈值，计算K阶矩阵LK(t,ω)并根据阈值得到m值；构建m阶矩阵Lm(t,ω)和Mm(t,ω)；计算在二维时频谱上的K阶瞬时频率估计值再用K阶瞬时频率估计值对短时傅里叶变换值进行挤压。本发明根据信号的特征灵活地选择合适的时窗函数，得到聚集的短时傅里叶变换时频谱，且能够直接求解任意阶估计值，有效提高了计算效率，提升了瞬时频率估计值的精确度，使得高阶同步挤压变换的时频谱能量更为集中，极大地提高了信号的时间和频率分辨率，同时在挤压过程中也保持了该变换的可逆性，能够高精度的重构原信号。