[发明专利]一种判断航行舰船地震波信号线谱特征的方法

说明书

技术领域

本发明属于航行舰船地震波探测领域，涉及到一种判断航行舰船地震波信号线谱特征的方法，特别涉及一种利用线谱频率、线谱窄度和线谱信噪比进行综合判断的方法。

背景技术

随着新技术的不断应用，舰船辐射噪声得到很大的抑制，所以需要采用多种物理场进行综合探测。航行舰船地震波作为舰船辐射噪声在海底的二次传播，在某些特殊的海洋环境条件下，可以作为舰船声探测的重要补充手段。然而由于海水运动和海底微地震等干扰的存在，容易将目标信号混淆，使得数据的重要目标信息难以有效提取。

目前判断航行舰船地震波信号特征的主要方法是通过对时域信号进行频谱分析，在频谱图中观察所测数据的线谱特征，将这些线谱与其它稳定背景频点幅值进行比较，并设置一定的门限，从而得到相应的线谱特征。虽然传统的线谱特征提取方法容易操作，但也容易出现误判断，将海洋环境干扰线谱也纳入到信号特征中来，引起虚警，因此需要对航行舰船地震波信号线谱特征提取方法进行规范，减小虚警概率。

发明内容

本发明的目的是提供一种判断航行舰船地震波信号线谱特征的方法，解决利用地震波线谱信号准确检测航行舰船目标的问题。

本发明的技术方案是：表征线谱特征的三个基本要素为线谱频率、线谱窄度和线谱信噪比；首先对测试得到的时域数据，进行快速傅立叶变换，得到时域数据的幅频曲线；将幅频曲线换算成地震波级形式，采用双向曲线拟合方法计算地震波级的一次趋势项和二次趋势项；将地震波级频谱与地震波级频谱的二次趋势项相减，并设置门限，得到疑似线谱频率以及所对应的信噪比值；设置线谱窄度，从地震波级频谱一次趋势项中确定满足窄度要求的线谱频率；依据疑似线谱频率、设定的线谱窄度和地震波级频谱一次趋势项，确定航行舰船地震波真实的线谱频率和窄度，并由此从线谱频率参考曲线中确定线谱对应的信噪比值。

传统的判断航行舰船地震波线谱信号的方法主要是基于经验和简单分析，给出疑似线谱频率，并不能准确的判断航行舰船地震波信号的线谱特征；结合线谱频率、线谱窄度和线谱信噪比综合给出航行舰船地震波线谱信号特征，可以确定是否真正检测到信号而不是环境地震波干扰。

附图说明

图1航行舰船地震波级频谱曲线以及趋势项。

图中：1地震波级频谱曲线；2地震波级频谱一次趋势项；3地震波级频谱二次趋势项。

图2航行舰船地震波疑似线谱以及确定的线谱频率。

图中：4疑似线谱；5线谱频率参考曲线；6确定线谱；2地震波级频谱一次趋势项。

具体实施方式

第一步：幅频曲线提取

在采集得到的航行舰船地震波数据中截取长度为30s的数据，进行快速傅立叶变换的数据点数为整个数据的长度，快速傅立叶变换计算公式为：

式中，X为快速傅立叶变换值，k为正整数并在1和n之间取值，n为作快速傅立叶变换的点数，x为采集的原始数据，e为自然对数底，π为圆周率，i为虚数单位。计算之后得到快速傅立叶变换值X，对X取绝对值并利用公式 X/(n/2)计算得到二次变换值Fd，在Fd中取第一点到第n/2点的数值作为最终提取的频谱幅值Fa。频谱曲线对应的频点如下式所示：

Ff(r)＝r·fs/n(r＝0,1,...,n-1)

式中，Ff为快速傅立叶变换值对应频点，r为整数并在0和n-1之间取值， fs为数据采样频率，在Ff中取第一点到第n/2点的数值作为频谱幅值对应频点 Fw。得到了频谱幅值Fa和对应频点值Fw便可以确定航行舰船地震波信号的幅频曲线。

第二步：一次趋势项提取

将频谱幅值Fa换算成地震波级的形式，计算公式如下：

Ga＝20·log₁₀(Fa)+60

选择趋势项拟合的频谱宽度参数N_f为0.2·(fs/n)并取整数部分，分别从幅频曲线的两端拟合各自的趋势项。从最低频点拟合的趋势项G₁′的计算公式如下：

从上式中可以得到第1个至第N-N_f个点的地震波级G₁′，为了保证G₁′和Ga的数据长度一致，对于G₁′的第N-N_f+1个点至第N个点的数值直接由Ga对应的第 N-N_f+1个点至第N个点的数值进行赋值。

从最高频点拟合的趋势项G₂′的计算公式如下：