CN111024207B - 一种矢量水听器线谱自动检测与判决方法 - Google Patents

Abstract

本发明公开了一种矢量水听器线谱自动检测与判决方法，涉及水声信号处理与信号检测技术领域，是能够通过幅度与方位的变化判断运动目标的存在，有效发挥水下目标线谱检测的应有效果，其检测技术可靠、虚警率低。本发明的技术方案包括如下步骤：获取矢量水听器接收的每一帧矢量数据，求解每一帧数量数据的平均声强谱，采用斜率门限与峰高门限从平均声强谱中提取线谱，线谱由线谱频率组成。从第2帧矢量数据开始，将每一帧矢量数据的线谱频率与之前帧矢量数据的线谱频率相比较，其中满足多普勒容限的线谱频率归为同一目标的线谱。对属于同一目标的线谱频率进行比较判断，若存在正向幅度变化以及单调方位变化，则判断目标存在，发出报警信号。

Description

一种矢量水听器线谱自动检测与判决方法

技术领域

本发明涉及水声信号处理与信号检测技术领域，具体涉及一种矢量水听器线谱自动检测与判决方法。

背景技术

水声监测小平台如噪声监测潜标一般利用水下目标的宽带信号检测目标。随着消声降噪技术的发展，水下目标辐射的宽带能量越来越低，给传统监测小平台采用的宽带检波积分检测方法带来很大困难。而水下目标低频线谱特征很难消除，所以低频线谱检测方法受到重视。现有的小平台线谱检测技术存在的主要问题是：线谱检测方法通过在某一窄带在连续一段时间内存在稳定的能量来判断线谱存在，容易把浅海环境下常见的稳定干扰判决为目标，导致虚警率较高，严重影响监测效果。

因此目前缺少一种可靠的、虚警率较低的针对小平台线谱检测技术。

发明内容

有鉴于此，本发明提供了一种矢量水听器线谱自动检测与判决方法，是能够通过幅度与方位的变化判断运动目标的存在，有效发挥水下目标线谱检测的应有效果，其检测技术可靠、虚警率低。

为达到上述目的，本发明的技术方案包括如下步骤：

获取矢量水听器接收的每一帧矢量数据，求解每一帧数量数据的平均声强谱，采用斜率门限与峰高门限从平均声强谱中提取线谱，线谱由线谱频率组成。

从第2帧矢量数据开始，将每一帧矢量数据的线谱频率与之前帧矢量数据的线谱频率相比较，其中满足多普勒容限的线谱频率归为同一目标的线谱。

对属于同一目标的线谱频率进行比较判断，若存在正向幅度变化以及单调方位变化，则判断目标存在，发出报警信号。

进一步地，获取矢量水听器接收的每一帧矢量数据，求解每一帧数量数据的平均声强谱，采用斜率门限与峰高门限从平均声强谱中提取线谱，线谱由线谱频率组成，具体为：

S0、对矢量水听器接收的一帧矢量数据包括声压p(n)、振速v_x(n)、振速v_y(n)，其中n＝0,…,N-1，n为采样点，N为一帧数据的长度。

矢量水听器的振速轴包括x轴y轴和z轴，其中x轴和y轴为水平面上的轴、z轴为垂直方向的轴。

顺次取每一帧矢量数据作为当前处理帧，执行如下步骤S1～S3：

S1、把当前处理帧的N点数据分为L段，每段M点，段与段之间有重叠，重叠部分占每段设定比例如50％，则分段数为

式中符号

表示向下取整；

当前处理帧的第l段(l＝0,1,…,L-1)数据的互功率谱为

式中P_l(k)是p(n)数据第l段(l＝0,1,…,L-1)的快速傅里叶变换(FFT)结果， V_xl(k)是v_x(n)数据第l段的FFT结果，V_yl(k)是v_y(n)数据第l段的FFT结果。

分别对L段C_xl(k)与C_yl(k)做平均，得到x轴与y轴方向的平均互功率谱为 C_x(k)与C_y(k)，k为频率点；分别取C_x(k)与C_y(k)实部，得到I_x(k)与I_y(k)。

S2、对当前处理帧求平均声强谱I₀。

S3、求平均声功率谱I₀(k)的连续谱背景B(k),k＝0,1,…,M-1；求平均声功率谱I₀(k)的差分得到Y(k),k＝0,1,…,M-2。对Y(k),k＝0,1,…,M-2逐点与设定的斜率门限进行比较，记录所有Y(k-1)>C₀且Y(k)<-C₀的点的序号k_i和相应的声功率谱幅度I₀(k_i)，其中C₀为斜率门限，为正数。

取满足条件I₀(k_i)/B(k_i)≥C₁的k_i，得到共K个满足条件的k_i，记为序列 k’₀,k’₁,…,k'_K-1。

其中C₁为设定的峰高门限。

k’₀,k’₁,…,k'_K-1为当前处理帧中线谱的数字频率，转换为实际频率为：

其中f_s是矢量水听器所处的水声监测系统对接收信号采样频率。

f_j为线谱频率，j＝0,1,…,K-1，则K个组成当前处理帧的线谱。

进一步地，从第2帧矢量数据开始，将每一帧矢量数据的线谱频率与之前帧矢量数据的线谱频率相比较，其中满足多普勒容限的线谱频率归为同一目标的线谱，具体为：

S4、取矢量水听器接收的第2帧矢量数据，第2帧数量数据的线谱中频率数量为K₂，将第2帧矢量数据得到的每个线谱频率均与第1帧的K₁个频率比较，在设定多普勒容限范围内的保留并归为同一目标的线谱，剩下的作为新出现的线谱保留；设定多普勒容限为

其中v为目标运动速度，c是水下声速， f是根据水听器探测目标的类型估计的线谱频率。

S5、取矢量水听器接收的第3帧矢量数据，第3帧数量数据的线谱中频率数量为K₃，将第3帧矢量数据得到的每个线谱频率均与第2帧的K₂个频率比较，在设定多普勒容限范围内的保留并归为同一目标的线谱，剩下的作为新出现的线谱保留。

S6、取矢量水听器接收的第R帧矢量数据，R＝4,5,…，长度为N点，第R 帧矢量数据的线谱中频率数量为K_R，把线谱频率依次与第R-1～R-p帧数据的线谱频率比较，在多普勒容限范围内的保留并归为同一目标的线谱，剩下的作为新出现的线谱保留。

p为设定比较帧数量门限，当R-p-1≥1时，对R-p-1帧数据的线谱，如果在后续R-p,…,R-2,R-1,R帧未再出现，则舍弃。

对于保留的线谱，记录其幅度I₀(k⁰)，同时根据下式计算方位

进一步地，对属于同一目标的线谱频率进行比较判断，若存在正向幅度变化以及单调方位变化，则判断目标存在，发出报警信号，具体为：

S7、针对同一目标A的线谱，从获得第q帧线谱结果开始，在每一帧中均有目标A的线谱频率时，如果第q-8、q-7以及q-6帧中对应该目标A的幅度均值为I₁，方位均值为φ₁，第q-5、q-4以及q-3帧中对应该目标A幅度均值为I₂，方位均值为φ₂，第q-2、q-1以及q帧中对应目标A幅度均值为I₃，方位均值为φ₃。

q至少取9。

且满足如下条件：I₃>I₂>I₁，且φ₃>φ₂>φ₁或者φ₃<φ₂<φ₁，则可判决目标A存在，可发出报警信号。

进一步地，S7之后还包括：

S8、若所有目标均不存在，则删掉第q-8帧的线谱结果，后续帧的线谱结果的帧序号依次前移。

继续获取矢量水听器的下一帧数据，按照S1～S4的方法提取下一帧数据的线谱，并将该帧的线谱结果作为第q帧线谱结果。

返回S7，直至能判决到存在目标为止。

有益效果：

本发明所述矢量水听器线谱自动检测与判决方法利用单个矢量水听器，采用线谱自动提取与相邻数据帧中线谱的关联性实现线谱检测，根据目标接近正横时线谱频率的幅度正向变化与方位单调变化特征，实现目标存在与否的判决，实现简单，且能有效防止传统线谱检测方法不考虑幅度与方位变化时虚警较高的缺点。

附图说明

图1为本发明提供的矢量水听器线谱自动检测与判决方法流程示意图；

图2为自动提取的若干帧数据线谱时间历程；

图3为目标方位变化历程；

图4为目标幅度变化历程；

图5为根据线谱方位与幅度变化给出判决启动信号。

具体实施方式

下面结合附图并举实施例，对本发明进行详细描述。

本发明提供了一种矢量水听器线谱自动检测与判决方法，其流程如图1所示，包括如下步骤：

步骤一、获取矢量水听器接收的每一帧矢量数据，求解每一帧数量数据的平均声强谱，采用斜率门限与峰高门限从平均声强谱中提取线谱，线谱由线谱频率组成。

本发明实施例中，该步骤包括如下具体的步骤：

S0、对矢量水听器接收的一帧矢量数据包括声压p(n)、振速v_x(n)、振速v_y(n)，其中n＝0,…,N-1，n为采样点，N为一帧数据的长度。

矢量水听器的振速轴包括x轴y轴和z轴，其中x轴和y轴为水平面上的轴、z轴为垂直方向的轴。

顺次取每一帧矢量数据作为当前处理帧，执行如下步骤S1～S3：

S1、把当前处理帧的N点数据分为L段，每段M点，段与段之间有重叠，重叠部分占每段设定比例，本发明实施例中，设置重叠50％，则分段数为

式中符号

表示向下取整。

当前处理帧的第l段(l＝0,1,…,L-1)数据的互功率谱为

式中P_l(k)是p(n)数据第l段(l＝0,1,…,L-1)的快速傅里叶变换(FFT)结果， V_xl(k)是v_x(n)数据第l段的FFT结果，V_yl(k)是v_y(n)数据第l段的FFT结果。

分别对L段C_xl(k)与C_yl(k)做平均，得到x轴与y轴方向的平均互功率谱为 C_x(k)与C_y(k)，k为频率点；分别取C_x(k)与C_y(k)实部，得到I_x(k)与I_y(k)。

S2、对当前处理帧求平均声强谱I₀。

S3、求平均声功率谱I₀(k)的连续谱背景B(k),k＝0,1,…,M-1；平均声功率谱由线谱叠加连续谱组成。采用“双通分离窗算法(Two-Pass Split-Window Algorithm)”得到平均声功率谱的连续谱背景B(k),k＝0,1,…,M-1。

求平均声功率谱I₀(k)的差分得到Y(k),k＝0,1,…,M-2。对Y(k),k＝0,1,…,M-2逐点与设定的斜率门限进行比较，记录所有Y(k-1)>C₀且Y(k)<-C₀的点的序号k_i和相应的声功率谱幅度I₀(k_i)，其中C₀为斜率门限，为正数，本发明实施例中取0.4，根据经验进行取值。

取满足条件I₀(k_i)/B(k_i)≥C₁的k_i，得到共K个满足条件的k_i，记为序列 k’₀,k’₁,…,k'_K-1。

其中C₁为设定的峰高门限，本发明实施例中C₁取2，经验值即线谱峰值比背景高6dB。

k’₀,k’₁,…,k'_K-1为当前处理帧中线谱的数字频率，转换为实际频率为：

其中f_s是矢量水听器所处的水声监测系统对接收信号采样频率。

f_j为线谱频率，j＝0,1,…,K-1，则K个组成当前处理帧的线谱。

步骤二、从第2帧矢量数据开始，将每一帧矢量数据的线谱频率与之前帧矢量数据的线谱频率相比较，其中满足多普勒容限的线谱频率归为同一目标的线谱。

S4、取矢量水听器接收的第2帧矢量数据，第2帧数量数据的线谱中频率数量为K₂，将第2帧矢量数据得到的每个线谱频率均与第1帧的K₁个频率比较，在设定多普勒容限范围内的保留并归为同一目标的线谱，剩下的作为新出现的线谱保留；设定多普勒容限为

其中v为目标运动速度，c是水下声速， f是根据水听器探测目标的类型估计的线谱频率。v对水下目标一般最大取8kn，水面目标取30kn。

由于目标由远到近通过观测点时，多普勒频偏由正到负变化，对于匀速直航目标，线谱频率总是变小的，所以多普勒容限负向取-Δf，正向可取小一些，如Δf/2。

S5、取矢量水听器接收的第3帧矢量数据，第3帧数量数据的线谱中频率数量为K₃，将第3帧矢量数据得到的每个线谱频率均与第2帧的K₂个频率比较，在设定多普勒容限范围内的保留并归为同一目标的线谱，剩下的作为新出现的线谱保留。

S6、取矢量水听器接收的第R帧矢量数据，R＝4,5,…，长度为N点，第R 帧矢量数据的线谱中频率数量为K_R，把线谱频率依次与第R-1～R-p帧数据的线谱频率比较，在多普勒容限范围内的保留并归为同一目标的线谱，剩下的作为新出现的线谱保留。

p为设定比较帧数量门限，当R-p-1≥1时，对R-p-1帧数据的线谱，如果在后续R-p,…,R-2,R-1,R帧未再出现，则舍弃。

对于保留的线谱，记录其幅度I₀(k⁰)，同时根据下式计算方位

步骤三、对属于同一目标的线谱频率进行比较判断，若存在正向幅度变化以及单调方位变化，则判断目标存在，发出报警信号。

S7、针对同一目标A的线谱，从获得第q帧线谱结果开始，在每一帧中均有目标A的线谱频率时，如果第q-8、q-7以及q-6帧中对应该目标A的幅度均值为I₁，方位均值为φ₁，第q-5、q-4以及q-3帧中对应该目标A幅度均值为I₂，方位均值为φ₂，第q-2、q-1以及q帧中对应目标A幅度均值为I₃，方位均值为φ₃，

q至少取9；

且满足如下条件：I₃>I₂>I₁，且φ₃>φ₂>φ₁或者φ₃<φ₂<φ₁，则可判决目标A存在，可发出报警信号。

之后还包括如下步骤：

S8、若所有目标均不存在，则删掉第q-8帧的线谱结果，后续帧的线谱结果的帧序号依次前移；

继续获取矢量水听器的下一帧数据，按照S1～S4的方法提取下一帧数据的线谱，并将该帧的线谱结果作为第q帧线谱结果；

返回S7，直至能判决到存在目标为止。

采用上述线谱的自动提取方法，得到线谱时间历程图如图2所示，对比线谱频率可知，本发明自动提取线谱的方法正确提取了平均声强谱中的线谱序列，线谱频率约250Hz。图3是目标运动过程中的方位变化历程，图4是目标运动过程中的线谱幅度变化历程。

图5是根据图3中250Hz频率附近线谱方位变化与图4中线谱幅度变化给出的目标报警启动信号，在280s时满足方位与幅度变化条件，给出高电平报警信号，说明本发明正确检测了目标线谱并在一定距离上启动了报警。

综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.一种矢量水听器线谱自动检测与判决方法，其特征在于，该方法利用单个矢量水听器，采用线谱自动提取与相邻数据帧中线谱的关联性实现线谱检测，包括如下步骤：

获取所述矢量水听器接收的每一帧矢量数据，求解每一帧数量数据的平均声强谱，采用斜率门限与峰高门限从所述平均声强谱中提取线谱，线谱由线谱频率组成；包括如下具体步骤：

S0、对所述矢量水听器接收的一帧矢量数据包括声压p(n)、振速v_x(n)、振速v_y(n)，其中n＝0,…,N-1，n为采样点，N为一帧数据的长度；

所述矢量水听器的振速轴包括x轴y轴和z轴，其中x轴和y轴为水平面上的轴、z轴为垂直方向的轴；

顺次取每一帧矢量数据作为当前处理帧，执行如下步骤S1～S3：

S1、把当前处理帧的N点数据分为L段，每段M点，段与段之间有重叠，重叠部分占每段设定比例，则分段数为

式中符号

表示向下取整；

当前处理帧的第l段数据的互功率谱如下：其中l＝0,1,…,L-1；

C_xl(k)＝P_l ^*(k)V_xl(k)

C_yl(k)＝P_l ^*(k)V_yl(k),k＝0,1,…,M-1

式中P_l(k)是p(n)数据第l段的快速傅里叶变换FFT结果，V_xl(k)是v_x(n)数据第l段的FFT结果，V_yl(k)是v_y(n)数据第l段的FFT结果；

分别对L段C_xl(k)与C_yl(k)做平均，得到x轴与y轴方向的平均互功率谱为C_x(k)与C_y(k)，k为频率点；分别取C_x(k)与C_y(k)实部，得到I_x(k)与I_y(k)；

S2、对当前处理帧求平均声强谱I₀；

S3、求平均声功率谱I₀(k)的连续谱背景B(k),k＝0,1,…,M-1；求平均声功率谱I₀(k)的差分得到Y(k),k＝0,1,…,M-2；对Y(k),k＝0,1,…,M-2逐点与设定的斜率门限进行比较，记录所有Y(k-1)>C₀且Y(k)<-C₀的点的序号k_i和相应的声功率谱幅度I₀(k_i)，其中C₀为斜率门限，为正数；

取满足条件I₀(k_i)/B(k_i)≥C₁的k_i，得到共K个满足条件的k_i，记为序列k’₀,k’₁,…,k'_K-1；

其中C₁为设定的峰高门限；

k’₀,k’₁,…,k'_K-1为当前处理帧中线谱的数字频率，转换为实际频率为：

其中f_s是矢量水听器所处的水声监测系统对接收信号采样频率；

f_j为线谱频率，j＝0,1,…,K-1，则K个组成当前处理帧的线谱；

从第2帧矢量数据开始，将每一帧矢量数据的线谱频率与之前帧矢量数据的线谱频率相比较，其中满足多普勒容限的线谱频率归为同一目标的线谱；包括如下具体步骤：

S4、取所述矢量水听器接收的第2帧矢量数据，第2帧数量数据的线谱中频率数量为K₂，将第2帧矢量数据得到的每个线谱频率均与第1帧的K₁个频率比较，在设定多普勒容限范围内的保留并归为同一目标的线谱，剩下的作为新出现的线谱保留；所述设定多普勒容限为

其中v为目标运动速度，c是水下声速，f是根据水听器探测目标的类型估计的线谱频率；

S5、取所述矢量水听器接收的第3帧矢量数据，第3帧数量数据的线谱中频率数量为K₃，将第3帧矢量数据得到的每个线谱频率均与第2帧的K₂个频率比较，在所述设定多普勒容限范围内的保留并归为同一目标的线谱，剩下的作为新出现的线谱保留；

S6、取所述矢量水听器接收的第R帧矢量数据，R＝4,5,…，长度为N点，第R帧矢量数据的线谱中频率数量为K_R，把线谱频率依次与第R-1～R-p帧数据的线谱频率比较，在多普勒容限范围内的保留并归为同一目标的线谱，剩下的作为新出现的线谱保留；

p为设定比较帧数量门限，当R-p-1≥1时，对R-p-1帧数据的线谱，如果在后续R-p,…,R-2,R-1,R帧未再出现，则舍弃；

对于保留的线谱，记录其幅度I₀(k⁰)，同时根据下式计算方位

对属于同一目标的线谱频率进行比较判断，若存在正向幅度变化以及单调方位变化，则判断目标存在，发出报警信号；包括如下具体步骤：

S7、针对同一目标A的线谱，从获得第q帧线谱结果开始，在每一帧中均有目标A的线谱频率时，如果第q-8、q-7以及q-6帧中对应该目标A的幅度均值为I₁，方位均值为φ₁，第q-5、q-4以及q-3帧中对应该目标A幅度均值为I₂，方位均值为φ₂，第q-2、q-1以及q帧中对应目标A幅度均值为I₃，方位均值为φ₃，

q至少取9；

且满足如下条件：I₃>I₂>I₁，且φ₃>φ₂>φ₁或者φ₃<φ₂<φ₁，则可判决目标A存在，可发出报警信号；

所述S7之后还包括：

S8、若所有目标均不存在，则删掉第q-8帧的线谱结果，后续帧的线谱结果的帧序号依次前移；

继续获取矢量水听器的下一帧数据，按照S1～S4的方法提取所述下一帧数据的线谱，并将该帧的线谱结果作为第q帧线谱结果；

返回S7，直至能判决到存在目标为止。

Images