[发明专利]一种分量数时变的水声信号瞬时频率估计方法

说明书

本发明涉及一种分量数时变的水声信号瞬时频率估计方法，包括：输入多分量水声信号s(t)，利用扩展的同步压缩变换ESST得到多分量水声信号s(t)的时频谱T(t，ω)；利用自动脊线提取AREA，追踪到时频谱T(t，ω)中所有可能的脊线，然后利用脊线融合标准，融合脊线；对得到的脊线用傅里叶级数展开，用Tikhonov正则化对脊线的傅里叶系数重构，得到重构的瞬时频率。本发明是对自动脊线提取，与传统重叠分量的瞬时频率估计相比，本发明能够处理分量数随时间变化的信号；本发明结合分数时频谱，采用扩展同步压缩变换提取高能量集中度的时频谱，能够克服噪声干扰影响水声信号瞬时频率估计的精度；本发明方法能够准确的估计水声信号的瞬时频率。

技术领域

本发明涉及水声信号处理技术领域，尤其是一种分量数时变的水声信号瞬时频率估计方法。

背景技术

在海洋实时监测、异常声波事件检测等应用中，通常需要重构特定分量的信息。瞬时频率是时频变换的一个重要属性，通过时频谱和瞬时频率，能重构水声信号所有分量的信息。因此，估计时频谱中的瞬时频率，能精确提取信号中特殊的目标信息。

对于多分量的信号，最为直接的方式是经验模式分解或变模式分解，将非平稳信号分解成一系列固有模式函数，所有模式参数可以用于重构分量的瞬时频率。然而，对于具有紧邻甚至频点重叠特性的水声信号，这些方法难以提取出所有分量的瞬时频率。一种替代方法是时频分析方法，如短时傅里叶变换、调频变换和二次时频变换方法，提取信号的时频谱。在时频平面，用时频变换的拟变换重构出分量信息。由于这些方法受噪声干扰，提取的时频分辨率低，估计的瞬时频率误差较大。为了提高时频分辨率，采用重分配方法改进时频谱的可读性，但是重分配方法破坏了时频变换的可逆性。为了解决这个问题，同步压缩变换采用频域重分配方法，在增强时频谱能量集中度的同时，保留时频变换的可逆性。不幸的是，同步压缩变换逆变换受分离条件限制，不能重构出紧邻甚至重叠分量的瞬时频率。固有模式调频分解、固有调频分量分解、核稀疏学习方法能估计具有重叠分量信号的瞬时频率，但是必须设定好信号分量数。

由于水声信号具有快时变、低信噪比、紧邻甚至频点重叠特性，至今为止，对于这种水声信号中分量数随时间变化的情况，还没有采用同步压缩变换和自动脊线提取相结合进行瞬时频率估计。

发明内容

本发明的目的在于提供一种能够处理分量数随时间变化的信号，克服噪声干扰影响水声信号瞬时频率估计的精度，准确地估计水声信号的瞬时频率的分量数时变的水声信号瞬时频率估计方法。

为实现上述目的，本发明采用了以下技术方案：一种分量数时变的水声信号瞬时频率估计方法，该方法包括下列顺序的步骤：

(1)输入多分量水声信号s(t)，利用扩展的同步压缩变换ESST得到多分量水声信号s(t)的时频谱T(t，ω)；

(2)利用自动脊线提取AREA，追踪到时频谱T(t，ω)中所有可能的脊线，然后利用脊线融合标准，融合脊线；

(3)对得到的脊线用傅里叶级数展开，用Tikhonov正则化对脊线的傅里叶系数重构，得到重构的瞬时频率。

在步骤(1)中，所述扩展的同步压缩变换ESST的公式如下：

式中，为公共实数集，h(0)为0时刻的高斯窗函数值，δ(·)为delta函数，x^[N](t，f)为第N阶同步压缩变换SST算子，S_MSTFFT(t，f)采用最大化短时间分数傅里叶变换的时频谱，定义为：