一种二维振速梯度水听器的被动目标相关检测方法
【专利摘要】本发明公开了一种二维振速梯度水听器的被动目标相关检测方法。该方法利用振速水听器的接收数据，先计算出中心点处的振速和振速梯度，后对二者进行互相关处理，利用互相关处理器的输出来检测是否存在被动目标，并根据不同角度上检测概率的分布粗略估计目标的方位。有益效果是：本发明充分利用了振速和振速梯度的指向性，以及噪声及目标信号振速与振速梯度信息相关性的差异性特点，降低了被动目标的可检测信噪比，从而有效的提高被动目标的探测距离。本发明方法对水下弱目标的远程探测与参数估计均有重要意义，可以在浮标、潜标、鱼探仪、水下航行器、水下定位与导航等产品中应用。
【专利说明】
-种二维振速梯度水听器的被动目标相关检测方法
技术领域
[0001] 本发明设及的是一种二维振速梯度水听器的被动目标相关检测方法，具体的说， 是一种利用二维振速梯度水听器中屯、点处的振速和振速梯度信息，通过互相关处理来检测 是否存在被动目标的方法。
【背景技术】
[0002] 随着科学的发展、技术的进步，军用、民用海洋装备对小孔径、远程、高精度的水下 探测设备需求增加，运使得水听器检测技术不断发展进步，也促使新的水下传感技术发展， 如近十几年来快速发展的矢量水听器技术。矢量水听器可同时共点的获取声场的声压（标 量）、振速（矢量）信息。振速具有与频率无关的偶极子指向性，且声压和振速信息联合处 理可W有效地抑制噪声，提高被动目标的探测能力。在低频段，矢量水听器相比于声压水听 器阵具有体积小，重量轻，便于布放等优点，广泛应用于军事探潜、海洋生物研究、地震监测 和环境噪声测量等系统中。
[0003] 为了得到指向性更加尖锐的传感器，很多学者对振速梯度水听器开展了系列研 究，已经通过理论分析和实场实验证明了振速梯度水听器可W获得更好的指向性。如何利 用振速梯度水听器实现被动目标的远程检测和有效参数估计，还有待开展进一步的研究。 各向同性噪声的振速和振速梯度是不相关的，而目标信号的振速与振速梯度是相干的。如 果有一种处理方法或技术，能够将振速与振速梯度信息联合处理，则可提高振速梯度水听 器的检测性能，有效增强其探测距离，运种技术必然会扩展振速梯度水听器的应用范围，特 别是在水下小孔径测量和监视系统中的应用。
【发明内容】
[0004] 本发明的目的是：利用目标信号和环境噪声信号振速与振速梯度相关信息的差异 性，提供一种简单实用、稳健可靠、可检测信噪比低的二维振速梯度水听器被动目标的相关 检测方法。 阳〇化]本发明的技术方案是：一种二维振速梯度水听器被动目标的相关检测方法，包括 如下的步骤：
[0006] 步骤一：将二维振速梯度水听器接收的八个通道实数信号数据V"、Vyi(i = 1，…，4)，通过希尔伯特变换器转换成复数信号咬,·、( ?
[0007] 步骤二：对四个通道的X轴复振速信号吃。=1:風永4)求平均，求取振速梯度水 听器中屯、点处X轴方向振速9。.,；
[000引步骤立：对四个通道的y轴复振速信号( /' = 1、2,3,4 )求平均，求取振速梯度水 听器中必点处y轴方向振速^
[0009] 步骤四：利用八通道振速复信号€:、(/' = 1，2,3,4 )求振速梯度水听器中屯、点 处横向振速梯度洗*,/奴、沸,./卸和纵向振速梯度孤,./郵、沉,/历Γ;:
[0010] 步骤五：利用引导方向α，对振速梯度水听器中屯、点处振速分量、^1,进行加权 处理，得到合成振速
W11] 步骤六：利用引导方向α，对振速梯度水听器中屯、点处横向振速梯度 挪1./如、36,/妙和纵向振速梯度W、./卸、进行加权处理，得到合成振速梯度
阳〇1引步骤屯：对Μ个采样点的合成振速(β)与合成振速梯度是(α)进行相关处理，得 到被动目标的检测统计量
[0013] 步骤八：根据系统要求的最大虚警概率和上述检验统计量D(a，Μ)，在噪声环境 中计算出不同引导方位上的检测口限Dth ( α，Μ);
[0014] 步骤九：实时计算D(a，Μ)，将D(a，Μ)与Dth(a，Μ)进行比较，如果D(a，Μ) > Dth(a，M)，则目标存在，反之目标不存在，并在不同的α上进行扫描，根据成功检测概率的 分布粗略估计目标的方位区域。
[0015] 本发明的有益效果是：本发明提供了一种稳健、低信噪比的振速梯度水听器被动 目标检测方法，根据水听器的接收数据，计算出中屯、点处组合振速和组合振速梯度，对二者 进行互相关处理，利用互相关处理器的输出来检测被动弱目标有无，并粗略估计目标的方 位区域，对水声探测与通信均有重要意义。
[0016] 本发明充分利用了振速和振速梯度的指向性，W及噪声和目标信号振速与振速梯 度相关信息的差异性特点，降低了被动目标的可检测信噪比，从而有效的提高被动目标的 探测距离。本发明方法对水下弱目标检测与参数估计具有重要意义，可W在浮标、潜标、鱼 探仪、水下航行器、水下定位与导航等产品中应用。
【附图说明】
[0017] 图1是本发明提出的二维振速梯度水听器被动目标检测的信号处理流程图。
[0018] 图2是二维振速梯度水听器的结构示意图。
[0019] 图3为实施例中目标信号中屯、频率fu= 1曲Ζ，系统采样频率为f ,= 8f。，噪声为 数值模拟的各向同性噪声，系统地虚警概率Pf= 0.05,目标在90°方位时，振速梯度水听 器振速与振速梯度互相关检测器检测成功概率统计结果图。
[0020] 图4为实施例中目标信号中屯、频率fu= 1曲Z，系统采样频率为f ,= 8f。，噪声为 数值模拟的各向同性噪声，系统地虚警概率Pf= 0. 05,目标在45°方位时，振速梯度水听 器振速与振速梯度互相关检测器检测成功概率统计结果图。
【具体实施方式】
[0021] 现结合实施例、附图对本发明做进一步描述：
[0022] 如图1所示，本发明利用四个对角振速传感器的平均和差分运算，计算出中屯、点 处的振速和振速梯度，并通过两者的互相关处理来检测是否存在目标。振速梯度水听器是 由对角布放的四个振速传感器组成，如图2所示。
[0023] 本实施例中振速梯度水听器通过振速与振速梯度相关处理进行被动目标检测的 方法，包括如下的步骤： W24] 步骤一：二维振速梯度水听器八个通道的接收信号用Vxl、Vyl、Vx2、Vy2、Vx3、Vy3、 Vw和Vy4来表示，通过希尔伯特变换器将接收的实信号变换成复信号f、l、^,.1、f、2、 巧d、％、起4和％4;
[0025] 步骤二：由四通道x轴复振速信号'^,、*、4求振速梯度水听器中屯、点 处X轴方向振廷
[0026] 步骤Ξ:由四通道y轴复振速信号P,.:、Cy3、f,4求振速梯度水听器中屯、点 处y轴方向振速
[0027] 步骤四：利用％、求解X轴横向振速梯度如-奴=j代，利用f，.3、巧Μ 求解y轴横向振速梯度
[00測步骤五：利用V求解纵向振速梯度滿,/觀=]於2-峰4巧，利用iW、起4求 解纵向振速梯适
[0029] 步骤六：利用引导方向α，对振速梯度水听器中屯、点处振速分量??、C。,进行加 权处理，得到合成振速
W30] 步骤屯：利用引导方向α，对振速梯度水听器中屯、点处横向振速梯度 、W,./妙和纵向振速梯度5C,./卸、进行加权处理，获得合成振速梯度
[0031] 步骤八：对Μ个采样点的合成振速V (α)与合成振速梯度进行相关处理，得 到被动目标的检测统计量
[0032] 步骤九：根据系统限定最大虚警概率Pf和上述检测统计量D(a，Μ)，在噪声环境 中统计计算出不同引导方向上的检测口限Dth(a，M);
[003引步骤十：在引导方向α上，实时计算D(a ,Μ)，将D(a ,Μ)与Dth(a ,Μ)进行比较， 如果D(a, Μ)〉Dth(a, Μ),则目标存在，反之目标不存在；在不同的α上进行扫描，根据检 测成功概率的大小可粗略估计目标存在的方位区域。
[0034] 图3是实施例振速梯度水听器半径为0. 01米，目标信号为中屯、频率为fe= 1曲Ζ 的CW脉冲信号，采样频率为8f。，噪声为数值模拟的各向同性噪声，系统地虚警概率Pf =0. 05,目标在90°方位时，单个振速梯度水听器振速与振速梯度互相关检测器检测成功 概率统计结果；图4是相同条件下，目标从45°方向入射时，单个振速梯度水听器振速与振 速梯度互相关检测器检测成功概率统计结果。为了比较，同时给出了声压能量检测器和矢 量水听器声能流检测器的成功检测概率随信噪比的变化曲线。从图3和图4可W看出，在 给定的仿真条件下，本发明方法（振速与振速梯度相关检测）的可检测信噪比低于矢量水 听器的声强检测法和声压水听器的能量检测法。运说明振速梯度水听器的振速与振速梯度 相关检测法对被动目标探测距离更远，在军用、民用水下探测设备中有广泛的应用价值。
[0035] 当然，本发明还可W有其他的实施例，熟悉本领域的技术人员可根据本发明做相 应的改变，但运些改变都应该属于本发明所附权利要求的保护范围。
【主权项】
1. 一种二维振速梯度水听器的被动目标相关检测方法，包括如下步骤： 步骤一：将二维振速梯度水听器八个通道的接收数据vxl、Vyl(i = 1，2,3,4)，通过希尔 伯特变换器转换成复数信号弋、％_ (i = l，2,3,4); 步骤二：对四个通道的X轴复振速信号^ (i = 1，2, 3, 4)求平均，求取振速梯度水听 器中心点处X轴方向振速t; 步骤三：对四个通道的y轴复振速信号匕（i = 1，2, 3, 4)求平均，求取振速梯度水听 器中心点处y轴方向振速^; 步骤四：利用八通道振速复信号t^ (i = 1，2, 3, 4)求振速梯度水听器中心点处 横向振速梯度、巧./#和纵向振速梯度％/%、; 步骤五：利用引导方向α，对振速梯度水听器中心点处振速分量I、&进行加权处 理，得到合成振速A (?) = t cos? +1\, sintt: 步骤六：利用引导方向α，对振速梯度水听器中心点处横向振速梯度、 A/办和纵向振速梯度化/w、进行加权处理，得到合成振速梯度 1 {a)-dvx /& c 齡2 α + r )，si it2 α 十0-ν' /以)cos α 沒} 步骤七：对Μ个采样点的合成振速^ (α)与合成振速梯度< (α):进行相关处理，得到被 动目标的检测统计量)=Σ||Μα，《)·0，《 )||: ?-1 步骤八：根据系统要求的最大虚警概率Pf和上述检验统计量D ( α，Μ)，在噪声环境中 计算出不同引导方位上的检测门限Dth(a，M); 步骤九：实时计算D( α，Μ)，将D( α，Μ)与Dth( α，Μ)进行比较，如果D( α，Μ) > Dth(a，M)，则目标存在，反之目标不存在，并在不同的α上进行扫描，根据成功检测概率的 分布粗略估计目标的方位。
【文档编号】G01S5/18GK106066472SQ201510606683
【公开日】2016年11月2日
【申请日】2015年9月22日 公开号201510606683.X, CN 106066472 A, CN 106066472A, CN 201510606683, CN-A-106066472, CN106066472 A, CN106066472A, CN201510606683, CN201510606683.X
【发明人】王绪虎, 张群飞, 吉爱国, 吴则举
【申请人】青岛理工大学