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

权利要求书

1.一种广义高阶同步挤压地震信号时频分解与重构方法，其特征在于：包括以下步骤：

(1)获取一地震信号x(t)；

(2)预设一阶数K和时窗函数g(t)，计算在窗函数t^kg(t)下的短时傅里叶变换值其中K≤6，k＝0,1,…,max{2K-2,1}；

(3)预设一阈值ε，根据下式构建K阶矩阵L_K(t,ω)，并查找矩阵L_K(t,ω)中所有顺序主子式det(L_m(t,ω))绝对值大于ε的最大值m；

m＝max{j:|det(L_j(t,ω))|＞ε,j＝1,2,…,K}

式中det表示求解矩阵行列式，|·|表示取绝对值；

(4)构建m阶矩阵L_m(t,ω)和M_m(t,ω)，其中，

(5)根据下式计算二维时频谱上的K阶瞬时频率估计值

式中，Re表示取复数的实部；

(6)根据下式对短时傅里叶变换值进行挤压，得到高阶同步挤压变换值

(7)采用下式对进行逆变换，重构地震信号x(t)；

2.根据权利要求1所述的一种广义高阶同步挤压地震信号时频分解与重构方法，其特征在于：步骤(6)得到高阶同步挤压变换值后，再对取模操作得到最终的高阶同步挤压变换时频谱

3.根据权利要求1所述的一种广义高阶同步挤压地震信号时频分解与重构方法，其特征在于：所述步骤(2)中，采用下式计算

式中，t为时间，ω为频率，σ为时窗函数的调节因子。

4.根据权利要求1所述的一种广义高阶同步挤压地震信号时频分解与重构方法，其特征在于：所述阈值ε，用于削弱步骤(5)对瞬时频率估计所产生的误差，使步骤(6)获得稳定的高度聚焦的时频表征结果。