一种广义高阶同步挤压地震信号时频分解与重构方法技术

本发明专利技术公开了一种广义高阶同步挤压地震信号时频分解与重构方法，包括输入一维信地震号；选择阶数和时窗函数，计算在窗函数下的短时傅里叶变换值，设置阈值，计算K阶矩阵L

【技术实现步骤摘要】
一种广义高阶同步挤压地震信号时频分解与重构方法
本专利技术涉及一种地震信号处理方法，尤其涉及一种广义高阶同步挤压地震信号时频分解与重构方法。
技术介绍
作为非平稳信号处理领域的一个重要分支，时频分析一直是现代信号处理的研究热点之一。常用的时频分析方法有短时傅里叶变换(STFT)、小波变换(CWT)、S变换(ST)及广义S变换(GST)等。虽然这些时频分析方法的分辨率都有所提高，但由于他们都是借助某一窗函数(不断改进的窗函数)并在该窗函数的区间内对信号的频率成分进行估计，使得其时频能量总是分布在一个以信号真实的瞬时频率为中心的范围内。因此，这些时频变换得到的时频谱的能量聚焦性不够集中。同时受Heisenberg-Gabor不定问题影响，其时间分辨率和频率分辨率不能同时达到最优。同步挤压变换(SST)是Daubechies等在小波变换的基础上提出的一种新的时频变换方法。它通过严格的数学推导，把小波变换结果在一定频率范围内的时频能量“挤压”到信号真实的瞬时频率附近，极大地提高了时频分辨率。而在同步挤压变换处理过程中，对信号真实的瞬时频率估计的本文档来自技高网...

【技术保护点】
1.一种广义高阶同步挤压地震信号时频分解与重构方法，其特征在于：包括以下步骤：/n(1)获取一地震信号x(t)；/n(2)预设一阶数K和时窗函数g(t)，计算在窗函数t

【技术特征摘要】
1.一种广义高阶同步挤压地震信号时频分解与重构方法，其特征在于：包括以下步骤：
(1)获取一地震信号x(t)；
(2)预设一阶数K和时窗函数g(t)，计算在窗函数tkg(t)下的短时傅里叶变换值其中K≤6，k＝0,1,…,max{2K-2,1}；
(3)预设一阈值ε，根据下式构建K阶矩阵LK(t,ω)，并查找矩阵LK(t,ω)中所有顺序主子式det(Lm(t,ω))绝对值大于ε的最大值m；



m＝max{j:|det(Lj(t,ω))|＞ε,j＝1,2,…,K}
式中det表示求解矩阵行列式，|·|表示取绝对值；
(4)构建m阶矩阵Lm(t,ω)和Mm(t,ω)，其中，



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



式中，Re表示取复数的实部；
(6)根据下式对短时傅里叶变换值进行挤压，得到高阶同步挤压变换值
<...

【专利技术属性】
技术研发人员：陈辉，陈旭平，胡英，周怀来，朱冰雪，杨超，陈扬，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：四川;51