一种双层圆柱阵水下被动目标检测方法
【技术领域】
[0001] 本发明属于阵列信号处理领域，尤其涉及一种双层圆柱阵水下被动目标检测方 法。
【背景技术】
[0002] 利用目标辐射的噪声进行探测的声纳被称为被动声纳，具有隐蔽性高和能进行目 标识别等优点，是对各类舰船和潜艇等水面和水下目标进行探测的主要手段。被动声呐探 测的首要任务是通过目标辐射的声信号发现目标，并确定目标的距离、方位、速度和类型等 信息，分别对应目标信号检测、参数估计和目标识别的技术。在声呐探测系统的信号处理流 程中，目标信号检测通常位于系统的前段，是目标分类和定位的前提和基础。
[0003] 目前水下被动目标信号检测的技术主要分为两种：一种是基于时域或频域的检测 技术，主要应用于复杂噪声环境下的特定声信号的检测，常常需要背景噪声和目标特性等 先验知识的帮助，才能获得较好的性能；另一种是基于空域的信号检测技术，主要利用信号 与噪声和干扰的空间分布的不同来进行目标信号的检测。基于空域信息的检测技术常常与 探测系统的接收前端密切结合，在实际工程中得到较为广泛的应用。
[0004] 水下目标的辐射噪声由于海水的吸收衰减，经过长距离传播后，到达接收系统时 的信号强度往往比较微弱。利用基阵增加目标辐射噪声信号的空间增益，以提高探测距离， 是水下目标探测普遍采用的方法。传统的采用线阵或者圆阵等二维阵进行水下目标探测比 较常见，但随着对水下目标探测技术指标的进一步提高，传统的方法已经很难满足需求，采 用更复杂的组合体积阵成为一种趋势。采用体积阵，可以在基阵上安放更多的阵元，在进一 步提高空间增益的同时，可以采用更加灵活多样的处理算法，有助于提高综合探测性能。
[0005] 现有的基于空域的检测技术，往往需要事先根据方位历程图估计噪声门限。然而 一方面：方位历程图在系统的处理流程中，往往是在检测到目标信号之后才进行的操作； 另一方面水下噪声的时变和不均匀分布导致信号检测的性能很难获得较好的效果。
【发明内容】
[0006] 本发明的目的是提供一种具有强稳健性的，双层圆柱阵水下被动目标检测方法。
[0007] -种双层圆柱阵水下被动目标检测方法，包括以下步骤，
[0008] 步骤一：接收基阵是由双层圆柱阵组成的体积阵，其柱面的内外圈半径分别为^ 和^，双层圆柱阵的内外圈圆柱分别分布N个线阵，每个线阵上等间隔分布着P个阵元，内 外圈圆柱阵的阵元输出信号为Xl]1(n)，i = 1，2分别代表内外圈，j = 1，2···Ν代表圆柱圆 周上阵元数目，1 = 1，2…Ρ代表每个线阵上的阵元数目；
[0009] 步骤二：将每个线阵的阵元输出相加，则圆柱阵等效成圆阵，得到2*Ν路输出 Υ?,?(η)：
[0011] 步骤三：根据双层圆阵输出的声信号yi] (η)，对内外层圆阵分别预成Μ个波束，获 得2Μ个波束：
[0013] 其中，分别用匕^,)和匕^,)表示内外圈的波束指向性函数，0k代表第k个预 成角度，9k=2Ji (k-l)/M，k= 1…M，波束与波束之间的夹角为Λ Θ =2jt/M，N代表圆周 上的阵元个数，fs是采样频率，
[0015] 其中，c = 1500m/s是水声声速，α = 2 π/N代表圆周上阵元之间的夹角；
[0016] 步骤四：将内外圈的相同编号的波束输出结果进行相关运算：
[0017] rk (m) = E {Ylk (n) Y2k (n+m)}, Mj < m < M 2
[0018] 其中，参数吣〈0〈]?2;
[0019] 步骤五：计算每个相关函数的峰值pk和所有波束相关函数的最大值p
[0020] pk= max {r k (m) | m < M 2}
[0021] pmax= max {p k | k = 1. . . M}；
[0022] 步骤六：进行目标判断；
[0023] (1)单目标情况，用所有峰值的平均值和最大值进行比较，根据阈值判断目标是否 存在，如果存在，目标方向与最大相关峰对应的波束方向最接近，
[0024] H。:目标存在，如果
[0025] H1:目标不存在，如果
[0026] 其中0〈 η '1是预设的阈值；
[0027] (2)多目标情况，如果目标数目为Μτ，且目标数目远小于波束个数，即Μ τ< < Μ，先 计算最大的Μτ个相关峰值，用所有峰值的平均值分别与这M 峰值进行比较，根据阈值判 断目标是否存在，如果存在，目标方向与相关峰值对应的波束方向接近，
[0028] H。:目标存在，如果
[0029] H1:目标不存在，如果：
[0030] 其中，Ρηι是所有相关峰值中最大的Μτ值，1彡m彡Μ τ，0〈 nM〈l是预设的阈值。
[0031] 有益效果：
[0032] 本发明的有益效果在于提供一种基于体积阵的水下被动目标空域检测方法，解决 了传统空域检测技术依赖方位历程图估计噪声阈值，受噪声起伏影响较大的缺陷。本方法 采用体积阵，具有比二维阵更高的探测增益。同时，利用本方法即可检测单目标也可检测多 目标，且检测结果可为后续的精确方向估计提供大致方位。
【附图说明】
[0033] 图1双层圆柱阵示意图；
[0034] 图2圆阵波束形成示意图；
[0035] 图3内外圈对应波束零点附近的互相关函数示意图；
[0036] 图4确定来波方位所在波束示意图；
[0037] 图5检测算法流程图。
【具体实施方式】
[0038] 下面将结合附图对本发明做进一步详细说明。
[0039] 本发明的目的是提供一种基于双层圆柱阵的水下被动目标信号空域检测技术。与 传统的空域信号检测技术相比：一是采用了更复杂的双层圆柱阵，具有更高的空间增益，可 探测更远的距离；二是所采用的检测算法利用有无信号时接收能量的空间分布差别进行检 测，与噪声起伏无关，不需要通过方位历程图进行预估噪声门限，因此具有更好的稳健性。
[0040] 本发明的目是这样实现的：
[0041] (1)将双层圆柱阵的内外圈圆柱分别等效成N个线阵，每个线阵上等间隔分布着P 个阵元，假设双层圆柱阵每个阵元的输出为XiuOi), i = 1，2分别代表内外圈，j = 1，2…N 代表圆柱圆周上阵元数目，1 = 1，2…P代表每个线阵上阵元数目。将每个线阵的阵元输出 相加，即每P路信号变成1路信号，则圆柱阵等效成圆阵，得到M*N路输出yi] (η):
[0043] (2)根据双层圆阵输出的声信号yi] (η)，对内外层圆阵分别预成Μ个波束，获得2Μ 个波束。分别用131(014)和132(014)表示内外圈的波束指向性函数，0 14代表第1^个预成角 度，Θ k= 2 31 (k-l)/M，k = 1···Μ，波束与波束之间的夹角为2 31 /Μ，波束形成的公式如（1) 所示：
[0045] 其中，Ν代表圆周上的阵元个数，fs是采样频率，i = 1，2分别代表内外圈。
[0047] 其中，c是水声声速，a = 23T/N代表圆周上阵元之间的夹角。
[0048] (3)将内外圈的相同编号的波束输出结果进行相关运算，理论上只需计算零点的 相关函数，为了降低干扰的影响，计算零点附近的一些相关函数值，计算公式如（4)所示：
[0049] rk (m) = E {Ylk (n) Y2k (n+m)}, Mj < m < M 2 (4)
[0050] 其中M'CKMz，即只需计算零点左右的几个相关值。
[0051] (4)计算每个相关函数的峰值pk和所有波束相关函数的最大值p
[0052] pk= max {r k (m) | m < M 2} (5)
[0053] pmax= max {p k I k = 1. · · M} (6)
[0054] (5)对于单目标情况，用所有峰值的平均值和最大值进行比较，超过预设的阈值就 可判断是否有目标存在，且目标方向与最大相关峰对应的波束方向最接近，可为后续的目 标精确方位估计奠定基础，即：
[0055] H。：目标存在，如果
[0056] H1:目标不存在，如果
[0057] 其中0〈 η '1是预设的阈值。
[0058] 对于多目标情况，如果目标数目为Μτ，且目标数目远小于波束个数，即Μτ< < Μ， 贝1J，可先计算最大的仏个相关峰值，用所有峰值的平均值分别与这^个峰值进行比较，超过 预设的阈值就可判断是否有目标存在，且目标方向与相关峰值对应的波束方向接近。
[0059] Η。：目标存在，如果
[0060] Η1:目标不存在，如果
[0061]