CN105279481A - 基于级联原子库的海洋目标辨识探测方法 - Google Patents

Abstract

基于级联原子库的海洋目标辨识探测方法，涉及海洋目标辨识探测。通过水下声信号本身及其回波信号的相关性建立过完备原子库，并根据其呈现的多分量辐射源信号的形式，选择时间尺度原子和时频原子组成级联原子库，将待识别信号自适应分解在原子库上，得到信号的最稀疏表示；通过建立时间尺度原子库并结合radon变换来处理水下信号并进行时间尺度原子参数t₀估计；在估计出t₀值的基础上，建立级联原子库，根据t₀的估计可以抽取出这一确定的行，将接收的信号分解到这一确定的行上，找出最大匹配值，进而可以直接提取出时频尺度原子参数k；利用级联原子库准确地提取出水下目标信号的主要特征，为最终的信号识别奠定坚实的基础，提高可靠性。

Description

基于级联原子库的海洋目标辨识探测方法

技术领域

本发明涉及海洋目标辨识探测，尤其是涉及一种基于级联原子库的海洋目标辨识探测方法。

背景技术

海洋占地球表面积的71％，在海底及海洋中，蕴藏着极其丰富的矿产资源和生物资源。在大海中有潜艇、鱼雷、暗礁、鱼群及生物群等各式各样的目标，如何通过大海中传播的声音去识别目标，将目标分类，这在军用或民用上都具有很重要的意义。国内的一些单位，如：西北工业大学、东南大学、哈尔滨工程大学、760所、715所、中科院声学所等在水下目标特性的提取、识别和神经网络分类器的设计等方面进行了大量的研究，比如基于IMF能量谱的水声信号特征提取与选择方法、基于关联维度特征的改进水下目标模式识别方法、基于Walsh谱统计分布特征的水下入侵目标分类方法、基于时频分析的S变换的水下目标识别方法，这几种方法虽然识别率尚可，复杂度低，但实时性和鲁棒性差；而基于小波包和支持向量机的水下目标识别方法、基于小波包和概率神经网络的水下目标识别方法以及采用小波统计特征的方法对舰船目标进行分类，这三种方法虽然识别率较高，鲁棒性较好，但是实时性较差，复杂度较高；基于单目视觉的水下目标三维信息重建方法、基于不相似度评价的水下声目标分类与听觉特征辨识、基于听觉谱的声音特征提取方法以及基于伽马通滤波器和应用听觉慢特征的水下目标识别方法，这几种基于听觉模型的水下目标识别方法虽然识别率较高且改善了以往的特征提取方法普遍存在的鲁棒性差等问题，但存在复杂度较高、实时性差的缺点。

发明内容

本发明的目的在于提供可在声波信号复杂、密集、多变且交叠严重的海洋环境中有效识别目标声信号的一种基于级联原子库的海洋目标辨识探测方法。

本发明包括以下步骤：

步骤1，通过水下声信号本身及其回波信号的相关性来建立过完备原子库，并根据其呈现的多分量辐射源信号的形式，选择时间尺度原子和时频原子组成级联原子库，将待识别信号自适应分解在原子库上，得到信号的最稀疏表示；

步骤2，通过建立时间尺度原子库并结合radon变换来处理水下信号并进行时间尺度原子参数t₀估计；

步骤3，在估计出时间尺度原子参数t₀值的基础上，建立级联原子库，根据t₀的估计可以抽取出这一确定的行，将接收的信号分解到这一确定的行上，找出最大匹配值，进而可以直接提取出时频尺度原子参数k；

步骤4，利用级联原子库准确地提取出水下目标信号的主要特征，为最终的信号识别奠定坚实的基础，提高可靠性。

在步骤1)中，所述级联原子库是将待识别信号自适应分解在时间尺度原子和时频原子构成的原子库上，其中：

(1)时间尺度原子

设信号为f(t)，在尺度s、位置μ的时间尺度变换运算W[]可以定义为以下形式：

$Wf(\mathrm{\μ},s)={\∫}_{-\mathrm{\∞}}^{+\mathrm{\∞}}f\left(t\right)\frac{1}{\sqrt{s}}\mathrm{\ψ}\left(\frac{t-\mathrm{\μ}}{s}\right)dt---\left(1\right)$

其中，f＝(μ,s)是原子库的参数组，是过完备原子库D中的原子，在过完备原子库D中选择与信号f(t)最匹配的原子；

(2)时频原子

在过完备原子库D中选择与信号f(t)最匹配的原子，其选择原则为内积最大，设第一个最匹配原子为g₁，则信号f(t)可以分解为两部分：原子g₁上的分量和第一次分解的残差R¹f，即：

f＝<f,g₁>g₁+R¹f(2)

接着对每次匹配的残差进行分解，直到终止条件，若对信号分解了n次，则可以表示成：

$f=\underset{k=0}{\overset{n-1}{\&Sigma;}}<R^{k}f,g_{r_k}>g_{r_k}+R^{n}f---\left(3\right)$

其中，R^kf是在第k个原子上的分量，Rⁿf是第k次分解的残差。

在步骤2)中，所述时间尺度原子参数t₀的估计的具体方法如下：

第一：建立声信号的时间尺度原子库：

其中，s是尺度参数，t是时间参数，τ_k是第k条路径的时延参数，整个原子库的个数为n×m；

第二：将信号分解在过完备原子库上：

把接收到的信号分解在过完备原子库上，信号在过完备原子库的某一个分量上出现最优匹配，现提取最优匹配；最大匹配是一条直线，其轨迹表达式为：

$a=\frac{1}{s_0}-\frac{1}{t_0}\&CenterDot;(b-{\mathrm{\&tau;}}_0)---\left(5\right)$

其中，a为原子库中的尺度因子匹配，b为原子库中的时延因子匹配，s₀是原子库的尺度参数，t₀是原子库的时间参数，即时间尺度原子参数；

第三：通过radon变换检测直线：

把以上匹配后的结果的信号经过radon变换后，会出现一个峰值点，通过该峰值点，可以估计出t₀；

在理想背景下，峰值出现在式(6)所示直线在(a,b)平面上的对应点(a₀,b₀)处，可得：

$a_0=\frac{1}{\sqrt{1+t_0^2}}\&CenterDot;(\frac{t_0}{s_0}+{\mathrm{\&tau;}}_0)---\left(6\right)$

$t_0=-\frac{1}{20\tan b}---\left(7\right).$

在步骤3)中，所述时频尺度原子参数k的估计的具体方法如下：。

建立关于估计参数t₀，k的级联原子库：

其中，s是尺度参数，t是时间参数，τ_k是第k条路径的时延参数，整个原子库的个数为n×m；

因为已经估计出t₀的值，所以可以建立级联原子库，再根据t₀的估计可以抽取出这一确定的行，将接收的信号分解到这一确定的行上，找出最大匹配值，进而可以直接估计出k。

虽然单个的原子库已经是过完备的，但是其仍然是由一个函数离散化得出，具有单一结果，通常只适合分解一类信号，进行单一参数的估计。而不同原子库对具有不同特点的信号分解时各有优劣。因此，若分解多分量信号并进行多参数估计，单个原子库分解必然会有一些局限，导致分解后的稀疏度不高。而使用级联原子库进行稀疏分解则能很好的解决这个问题。水下目标信号本身具有复杂，密集、多变，且交叠严重这些特点，在实际中更是呈现为多分量辐射源信号的形式，使得单一原子库下原子分解效果明显下降。因此，本发明提出基于级联原子库的原子分解方法，用来分解多分量的水下目标源信号。

本发明与常规海洋声目标识别相比具有以下优点：

本发明率先将建立原子库估计声信号参数与稀疏分解这种新的信号处理方法应用于海洋目标辨识技术中。通过水下声信号本身及其回波信号的相关性来建立过完备原子库，可进行未知信号的参数估计。同时，稀疏分解是一种根据信号特点进行信号自适应的分解方法，即将信号自适应分解在一组高度冗余的过完备基(即原子库)上，在一定的误差下的得到信号的最稀疏表示。由于原子库是过完备的，因此可以得到该信号的多个表示，而最稀疏的表示只有一个。稀疏分解就是为了寻找信号在该原子库下的最稀疏表示，即用最少的原子来表示该信号。从而可以有效地解决水下目标识别系统的识别率低、鲁棒性差、可靠性不高而不能满足新形势下的水下目标识别要求的问题。

综上所述，本发明可以有效地解决水下目标识别系统的识别率低、鲁棒性差、可靠性不高而不能满足新形势下的水下目标识别要求的问题。

本发明提供了实时性高、可靠且识别效果理想的水下目标信号识别系统，有效解决了现有水下声目标信号本身具有复杂，密集、多变，且交叠严重的特点，在实际中呈现为多分量辐射源信号的形式，使得水下信号识别效率低、实时性差等问题。因此，本发明提出基于级联原子库的时频原子分解方法，用来分解多分量的水下目标源信号。由于不同的时频原子库能稀疏表示适合分解的那一分量信号，因此可以得到多分量水下目标信号更稀疏的表示，能提取出水下目标信号的主要特征，为信号尚未识别提高了可靠性，以分辨海洋声信号。

附图说明

图1为本发明实施例的过完备原子库的最优匹配示意图。

图2为本发明实施例的最优匹配值经radon变换示意图。

图3为本发明实施例的时频原子k的最佳原子示意图。

具体实施方式

以下结合实施例和附图对本发明作进一步说明。

构造级联原子库的时候，需要充分考虑海洋声信号的特点，由信号的特点来选择合适的原子库组成级联原子库。组成级联原子库的各个原子库之间通常要求有不同的结构，以便于分解不同分量的信号。本发明所选择的级联原子库是时间尺度原子和时频原子。

(1)时间尺度原子

设信号为f(t)，在尺度s，位置μ的时间尺度变换运算W[]可以定义为以下形式：

$Wf(\mathrm{\&mu;},s)={\&Integral;}_{-\mathrm{\&infin;}}^{+\mathrm{\&infin;}}f\left(t\right)\frac{1}{\sqrt{s}}\mathrm{\&psi;}\left(\frac{t-\mathrm{\&mu;}}{s}\right)dt---\left(1\right)$

其中，f＝(μ,s)是原子库的参数组。是过完备原子库D的原子，在过完备原子库D中选择与信号f(t)最匹配的原子。

(2)时频尺度原子

设信号为f(t)，其长度为N，D为过完备原子库，g为原子库中经归一化的原子||g||＝1，信号分解时，在过完备原子库D中选择与信号f(t)最匹配的原子，其选择原则为内积最大，设第一个最匹配原子为g₁，则

$\left|<f,g_1>\right|=\underset{g\&Element;D}{sup}\left|<f,g>\right|---\left(2\right)$

已知原子满足归一化条件，因此，信号f可以分解为两部分：原子g₁上的分量和残差R¹f，即

f＝<f,g₁>g₁+R¹f(3)

接着对每次匹配的残差进行分解，直到终止条件，若对信号分解了n次，则可以表示成

$f=\underset{k=0}{\overset{n-1}{\&Sigma;}}<R^{k}f,g_{r_k}>g_{r_k}+R^{n}f---\left(4\right)$

基于级联原子库的时频原子分解方法在级联的2个原子库中搜索每个原子库中的最优原子，然后对搜索到的2个库中的最优原子进行比较，选择出最优原子。具体算法流程描述如下：

步骤1：将信号残差初始化为原信号；

步骤2：时间尺度原子种群初始化。随机产生一定数量的DNA链；

步骤3：进行DNA进化算法搜索时间尺度中的最优原子。主要包括选择、交叉和变异、灾变等操作，具体如下：

步骤3.1：按照DNA码表进行解码(见表1),将DNA链转化为实数；每3个碱基对应一位0到19的实数，如碱基UCA对应3，AAU对应12。再由解码后的实数映射到时间尺度原子参数，由时间尺度原子参数计算出种群中的每个时间尺度原子；

表1DNA密码表

步骤3.2：对种群中的每个时间尺度原子进行适应度值评价，并保存最优原子信息；

步骤3.3：对适应度值进行排序，然后按轮盘赌规则进行选择；其后按照一定概率进行交叉和变异；

步骤3.4：判断是否进行灾变。计算种群中DNA链的多样性值，若其小于预定值则进行灾变。灾变策略为保留5％的最优个体和步骤5％的次优个体，剩下的90％个体重新初始化；

步骤3.5：重复步骤3.1直到满足终止条件；

步骤4：时频原子种群初始化。并按照步骤3进行最优时频原子搜索；

步骤5：对搜索到的2个原子进行比较，选择适应度值最大为最优原子；

步骤6：将当前信号残差更新为原残差减去最优原子乘上最优原子与原残差的内积；

步骤7：重复步骤2直到满足分解的终止条件。

(3)海洋声信号的级联原子库设计及参数估计

时间尺度原子参数t₀的估计

下面通过建立时间尺度原子库并结合radon变换来处理双曲调频信号并进行参数估计。

第一：建立基HFM信号的时间尺度原子库

设发射信号为f(t)，其表达式为：

$\mathrm{\&mu;}\left(t\right)=\frac{1}{\sqrt{T}}rect\left(\frac{t}{T}\right)*\exp \{-j\&lsqb;2\mathrm{\&pi;}kln(1-\frac{t}{t_0})\&rsqb;\}---\left(5\right)$

为了在干扰背景下进行HFM信号的参数估计，通过建立关于尺度因子和时延因子的过完备原子库，进而实现具有不同尺度和时延的副本信号与接收信号的相关。下面建立过完备原子库。假设τ₀为时延参数，s₀为尺度参数，则可以构造出信号关于尺度和时延的过完备原子库G_f{g_r}：

其中，m,n为τ₀，s₀的集合，整个原子库的个数为n×m。

第二：将信号分解在过完备原子库上。

把接收到的信号分解在过完备原子库上，如图1所示。

信号在过完备原子库的某一个分量上出现最优匹配，提取最优匹配。

最优匹配中可以提取出其最大匹配，已经证明其最大匹配是一条直线，其轨迹表达式为：

$a=\frac{1}{s_0}-\frac{1}{t_0}.(b-{\mathrm{\&tau;}}_0)---\left(7\right)$

其中a为原子库中的尺度因子匹配，b为原子库中的时延因子匹配。

第三：通过radon变换检测直线。

把上面匹配后的结果的信号经过radon变换后，会出现一个峰值点，通过这个峰值点，可以估计出时间尺度原子参数t₀，如图2示。

在理想背景下，峰值出现在式(7)所示直线在(a,b)平面上的对应点(a₀,b₀)处，可得：

$a_0=\frac{1}{\sqrt{1+t_0^2}}.(\frac{t_0}{s_0}+{\mathrm{\&tau;}}_0)---\left(10\right)$

$t_0=-\frac{1}{20\tan b}---\left(11\right)$

由上述原理，估计时频尺度原子参数k时，可首先建立关于估计时间尺度原子参数t₀，k的级联原子库H_f{g_d}：

因为已经估计出t₀的值，所以可以建立级联原子库，再根据t₀的估计可以抽取出这一确定的行，将接收的信号分解到这一确定的行上，找出最大匹配值，进而可以直接提取出k，其最佳原子如图3示。

本发明提出了一种具有低计算复杂度、快实时性、高鲁棒性与识别效果理想的水下目标信号辨识系统。由于水下声目标信号本身具有复杂，密集、多变以及交叠严重的特点，在海洋环境中呈现为多分量辐射源信号的形式，使得水下信号识别效率低、实时性差。因此，本发明通过水下声信号本身及其回波信号的相关性来建立过完备原子库，将待识别信号自适应分解在原子库上，得到信号的最稀疏表示，估计此信号的参数，最终准确地提取出目标信号的主要特征，从而实现海洋信号辨识。在有潜艇、鱼雷、暗礁、鱼群及生物群等各式各样目标的海洋中，基于级联原子库的探测技术可以通过海洋传播的声音去识别目标，可以有效地解决水下目标识别系统的识别率低、鲁棒性差、可靠性不高而不能满足新形势下的水下目标识别要求的问题，这在军用或民用上都具有很重要的意义。

Claims (4)

1.基于级联原子库的海洋目标辨识探测方法，其特征在于包括以下步骤：

步骤1，通过水下声信号本身及其回波信号的相关性来建立过完备原子库，并根据其呈现的多分量辐射源信号的形式，选择时间尺度原子和时频原子组成级联原子库，将待识别信号自适应分解在原子库上，得到信号的最稀疏表示；

步骤2，通过建立时间尺度原子库并结合radon变换来处理水下信号并进行时间尺度原子参数t₀估计；

步骤3，在估计出时间尺度原子参数t₀值的基础上，建立级联原子库，根据t₀的估计可以抽取出这一确定的行，将接收的信号分解到这一确定的行上，找出最大匹配值，进而可以直接提取出时频尺度原子参数k；

步骤4，利用级联原子库准确地提取出水下目标信号的主要特征，为最终的信号识别奠定坚实的基础，提高可靠性。

2.基于级联原子库的海洋目标辨识探测方法，其特征在于在步骤1中，所述级联原子库是将待识别信号自适应分解在时间尺度原子和时频原子构成的原子库上，其中：

(1)时间尺度原子

设信号为f(t)，在尺度s、位置μ的时间尺度变换运算W[]可以定义为以下形式：

$Wf(\mathrm{\&mu;},s)={\&Integral;}_{-\mathrm{\&infin;}}^{+\mathrm{\&infin;}}f\left(t\right)\frac{1}{\sqrt{s}}\mathrm{\&psi;}\left(\frac{t-\mathrm{\&mu;}}{s}\right)dt---\left(1\right)$

其中，f＝(μ,s)是原子库的参数组，是过完备原子库D中的原子，在过完备原子库D中选择与信号f(t)最匹配的原子；

(2)时频原子

在过完备原子库D中选择与信号f(t)最匹配的原子，其选择原则为内积最大，设第一个最匹配原子为g₁，则信号f(t)可以分解为两部分：原子g₁上的分量和第一次分解的残差R¹f，即：

f＝<f,g₁>g₁+R¹f(2)

接着对每次匹配的残差进行分解，直到终止条件，若对信号分解了n次，则表示成：

$f=\underset{k=0}{\overset{n-1}{\&Sigma;}}<R^{k}f,g_{r_k}>g_{r_k}+R^{n}f---\left(3\right)$

其中，R^kf是在第k个原子上的分量，Rⁿf是第k次分解的残差。

3.基于级联原子库的海洋目标辨识探测方法，其特征在于在步骤2中，所述时间尺度原子参数t₀的估计的具体方法如下：

第一：建立声信号的时间尺度原子库：

其中，s是尺度参数，t是时间参数，τ_k是第k条路径的时延参数，整个原子库的个数为n×m；

第二：将信号分解在过完备原子库上：

把接收到的信号分解在过完备原子库上，信号在过完备原子库的某一个分量上出现最优匹配，现提取最优匹配；最大匹配是一条直线，其轨迹表达式为：

$a=\frac{1}{s_0}-\frac{1}{t_0}.(b-{\mathrm{\&tau;}}_0)---\left(5\right)$

其中，a为原子库中的尺度因子匹配，b为原子库中的时延因子匹配，s₀是原子库的尺度参数，t₀是原子库的时间参数，即时间尺度原子参数；

第三：通过radon变换检测直线：

把以上匹配后的结果的信号经过radon变换后，会出现一个峰值点，通过该峰值点，可以估计出t₀；

在理想背景下，峰值出现在式(6)所示直线在(a,b)平面上的对应点(a₀,b₀)处，可得：

$a_0=\frac{1}{\sqrt{1+t_0^2}}.(\frac{t_0}{s_0}+{\mathrm{\&tau;}}_0)---\left(6\right)$

$t_0=-\frac{1}{20\tan b}---\left(7\right).$

4.基于级联原子库的海洋目标辨识探测方法，其特征在于在步骤3中，所述时频尺度原子参数k的估计的具体方法如下：

建立关于估计参数t₀，k的级联原子库：

其中，s是尺度参数，t是时间参数，τ_k是第k条路径的时延参数，整个原子库的个数为n×m；

因为已经估计出t₀的值，所以可以建立级联原子库，再根据t₀的估计可以抽取出这一确定的行，将接收的信号分解到这一确定的行上，找出最大匹配值，进而可以直接估计出k。