一种自适应随机共振的地震波特征提取方法

摘要

本发明提供了一种自适应随机共振的地震波特征提取方法，采用最大似然估计估计噪声强度，建立二阶非线性双稳态随机共振系统，通过优化输出信噪比增益测度以及同步克莱姆逃逸率实现噪声匹配与频率匹配；按照设定的步长改变信号频率，选取信噪比最大值作为最优匹配值，并提取其对应的信号频率即为目标特征线谱频率。本发明能够大幅提升目标特征线谱处的信号能量及局部信噪比，对完全淹没的噪声中的特征信号亦能得到明显的增强输出。

权利要求

1.一种自适应随机共振的地震波特征提取方法，其特征在于包括下述步骤：步骤1，采用最大似然估计估计噪声强度，计算噪声方差 其中，N为信号点数，T(x)为检验统计量；

步骤2，建立二阶非线性双稳态随机共振系统 其中x

为系统输出，γ为二阶系统阻尼参数， 表示接收到的特征线谱信号，A0为输入信号幅值，f0为输入信号频率， 为初始相位；n(t)为噪声强度为D的背景噪声；非线性双稳态势函数 a＞0,b＞0，a、b为双稳态随机共振系统的势垒参数，系统势函数有两个稳态点在 处，被在零点处的势垒ΔV＝a2/4b分开；

步骤3，建立参数匹配的随机共振系统，通过优化输出信噪比增益测度以及同步克莱姆逃逸率实现噪声匹配与频率匹配；匹配参数 其中e为自然对数，信号频率f0通过设定信号频率搜索范围[fstart,fend]以及搜索步长fstep选取，D为噪声强度；

步骤4，将二阶非线性双稳态随机共振系统的二阶方程化简为两个一阶方程，然后用四阶龙戈库塔算法求解方程组，其中h为龙戈库塔步长参数；

步骤5，设定时间尺度搜索范围[hstart,hend]以及搜索步长hstep，根据设定的时间尺度搜索范围及搜索步长自适应寻找信号频率处最大的输出功率谱PMAX，并记录fstart所对应最优时间尺度hopt和PMAX值；

步骤6，计算局部信噪比测度 其中Pd表示特征频率d处的功率谱值，ΔB为信号频率附近的频带带宽， 表示局部噪声的能量；

步骤7，按照设定的步长fstep改变信号频率，重复步骤3～步骤6，选取信噪比最大值作为最优匹配值，并提取其对应的信号频率即为目标特征线谱频率f0。

2.根据权利要求1所述的自适应随机共振的地震波特征提取方法，其特征在于：所述的信号频率f0在海底地震波探测范围10Hz～100Hz中以搜索步长1Hz进行选取。

说明书

一种自适应随机共振的地震波特征提取方法

技术领域

[0001] 本发明涉及水下地震波信号提取和微弱信号检测领域。

背景技术

[0002] 随着水下大型航行器消磁降噪等措施的不断完善，其本身的自防护能力不断提高，尤其是安静型水下航行器无人潜航器等水下移动目标，其声学特性被不断削弱，传统声学探测方法难度增大。水中航行的舰艇是很大的能量载体，可以通过船体对水体的扰动、船体的固有频率振动、船上的机械振动辐射噪声、螺旋桨噪声以及水动力噪声等方式将能量释放到水体中向外传播。理论和实验已证明运动舰艇不平衡旋转部件、螺旋桨周期击水以及叶片共振会产生较高强度和稳定度的甚低频线谱特征，其中一部分能量被宽带空化噪声调制后向外辐射，另一部分能量经由水体耦合至海底，引起海底岩土层振动，进而产生舰艇地震波，可直接关联目标发动机、螺旋桨等物理特征。其作为水下大型航行器探测信号，携带水下航行目标的关键特征，可以不受水文条件的影响而通过海底传播，具有衰减慢、作用距离远且很难被抑制的特点。因此利用甚低频地震波信号进行水下目标探测，是实现对水下大型航行器远程探测的有效技术手段，是水下被动探测系统远距离感知目标的重点发展方向之一，在复杂的强海洋环境噪声中提取地震波信号特征显得尤为重要。

[0003] 近年来，伴随着非线性动力学和统计物理理论而突飞猛进，基于非线性随机共振的微弱周期信号检测理论方法得到快速发展，目前已在众多科学领域取得了丰硕的成果，其对于微弱周期信号的增强作用为低信噪比条件下的水中目标特征线谱检测与提取方法提供了一条新的途径。而随机共振技术则利用了噪声来增强目标信号，以提高微弱信号的检测性能。具体来说，随机共振技术通过建立合适的系统，将部分噪声能量转化为信号能量，能够大幅度提高输出信噪比。

[0004] 传统水中目标线谱特征检测与提取方法，如滤波、相关检测、时频分析等，能够处理一定信噪比的噪声，但是在强背景噪声干扰下其提取性能难度很大，无法在强背景噪声中提取微弱信号。

发明内容

[0005] 为了克服现有技术的不足，本发明提供一种自适应随机共振增强的微弱特征提取方法，通过建立自适应匹配随机共振模型，以最大局部信噪比为表征测度，在复杂的强海洋环境背景噪声中提取微弱的地震波特征信号。

[0006] 本发明解决其技术问题所采用的技术方案包括以下步骤：

[0007] 步骤1，采用最大似然估计估计噪声强度，计算噪声方差其中，N为信号点数，T(x)为检验统计量；

[0008] 步骤2，建立二阶非线性双稳态随机共振系统其中x为系统输出，γ为二阶系统阻尼参数， 表示接收到的特征
线谱信号，A0为输入信号幅值，f0为输入信号频率， 为初始相位；n(t)为噪声强度为D的背景噪声；非线性双稳态势函数 a＞0,b＞0，a、b为双稳态随机共振系统的
势垒参数，系统势函数有两个稳态点在 处，被在零点处的势垒ΔV＝a2/4b分开；

[0009] 步骤3，建立参数匹配的随机共振系统，通过优化输出信噪比增益测度以及同步克莱姆逃逸率实现噪声匹配与频率匹配；所述的匹配参数 其中e为自然对数，信号频率f0通过设定信号频率搜索范围[fstart,fend]以及搜索步长fstep选取，D为噪声强度；

[0010] 为了尽可能覆盖海底地震波的所有频段，信号频率f0在海底地震波探测范围10Hz～100Hz中以搜索步长1Hz进行选取。

[0011] 步骤4，将二阶非线性双稳态随机共振系统的二阶方程化简为两个一阶方程，然后用四阶龙戈库塔算法求解方程组，

[0012]

[0013] 其中h为龙戈库塔步长参数；

[0014] 步骤5，设定时间尺度搜索范围[hstart,hend]以及搜索步长hstep，根据设定的时间尺度搜索范围及搜索步长自适应寻找信号频率处最大的输出功率谱PMAX，并记录fstart所对应最优时间尺度hopt和PMAX值；

[0015] 步骤6，计算局部信噪比测度 其中Pd表示特征频率处的功率谱值，ΔB为信号频率附近的频带带宽， 表示局部噪声的能量；

[0016] 步骤7，按照设定的步长fstep改变信号频率，重复步骤3～步骤6，选取信噪比最大值作为最优匹配值，并提取其对应的信号频率即为目标特征线谱频率f0。

[0017] 本发明的有益效果是：在前期参数匹配随机共振的框架下进一步提出了一种自适应参数匹配随机共振微弱特征信号提取方法，仿真分析了其非线性滤波性能以及特征提取性能，并利用实测数据验证了本方法的实用性能。结果显示利用自适应随机共振可以大幅提升目标特征线谱处的信号能量及局部信噪比，对完全淹没的噪声中的特征信号亦能得到明显的增强输出，能够验证本方法对弱小目标关键特征提取的有效性。研究成果可为新型水中兵器引信、编队反潜作战系统、预警防护系统、我国海域安全监测网络构建等提供技术支撑，相关理论方法同样适用于其他各类水中舰艇目标，可用于提升水中兵器、水面及水下舰艇在复杂海洋环境噪声下对敌方弱、小目标的被动探测能力。

附图说明

[0018] 图1是本发明的自适应随机共振地震波特征提取方法原理框图。

[0019] 图2是本发明的仿真输入信号时域与频域图。

[0020] 图3是本发明的自适应随机共振系统输出信号时域与频域图。

[0021] 图4是本发明对实测水下航行器地震波数据的时频分析图。

具体实施方式

[0022] 下面结合附图和实施例对本发明进一步说明，本发明包括但不仅限于下述实施例。

[0023] 本发明建立自适应步长的随机共振信号检测系统，同时针对特征频率未知的问题，建立频率匹配的搜索提取方法。

[0024] 具体包括以下步骤：

[0025] 步骤1：噪声强度估计

[0026] 采用最大似然估计的方法，估计噪声强度，即只有噪声输入时的方差，计算公式如下：

[0027]

[0028] 其中， 为噪声的方差，N为信号点数，T(x)为检验统计量；

[0029] 步骤2：随机共振系统构建设计

[0030] 建立二阶非线性双稳态随机共振系统：

[0031]

[0032] 其中γ为二阶系统阻尼参数， 表示接收到的特征线谱信号，A0为输入信号幅值，f0为输入信号频率， 为初始相位；n(t)为噪声强度为D的背景噪声；V(x)为非线性双稳态势函数；

[0033]

[0034] 其中a、b为双稳态随机共振系统的势垒参数，系统势函数有两个稳态点在处，被在零点处的势垒ΔV＝a2/4b分开。

[0035] 步骤3：选取起始信号及匹配系统参数

[0036] 快速建立参数匹配的随机共振系统，通过优化输出信噪比增益测度以及同步克莱姆逃逸率实现噪声匹配与频率匹配，其匹配参数选取由式(4)计算可得：

[0037]

[0038] 其中e为自然对数，信号频率f0通过设定信号频率搜索范围[fstart，fend]以及搜索步长fstep中选取，噪声强度D可以通过经典统计信号处理理论估计得到，由于实际接收信号的噪声强度一般都比较大，因此可以简单设定参数a＝1以便在任意环境噪声下快速建立参数匹配的随机共振系统。

[0039] 步骤4：系统数值求解

[0040] 首先，将二阶方程化简为两个一阶方程，

[0041]

[0042] 然后用如下四阶龙戈库塔算法(RK4)求解方程组：

[0043]

[0044] 其中h为龙戈库塔步长参数。

[0045] 步骤5：搜索最优步长参数

[0046] 设定时间尺度搜索范围[hstart,hend]以及搜索步长hstep，根据设定的时间尺度搜索范围及搜索步长自适应寻找信号频率处最大的输出功率谱PMAX，并记录fstart所对应最优时间尺度hopt和PMAX值。

[0047] 步骤6：系统输出局部信噪比与存储

[0048] 以局部信噪比测度指标，局部信噪比测度的表示如下：

[0049]

[0050] 其中Pd表示特征频率处的功率谱值，ΔB为信号频率附近的频带带宽，表示局部噪声的能量。

[0051] 将信噪比及其对应的信号频率和输出信噪比。

[0052] 步骤7：改变信号频率，执行步骤3到步骤6。

[0053] 步骤8：在存储的信噪比中选取最大值作为最优匹配值，并提取其对应的信号频率即为目标特征线谱频率f0。

[0054] 步骤9：输出信号及目标特征线谱频率f0。