发明名称 一种二维振速梯度水听器的被动目标相关检测方法
摘要 本发明公开了一种二维振速梯度水听器的被动目标相关检测方法。该方法利用振速水听器的接收数据,先计算出中心点处的振速和振速梯度,后对二者进行互相关处理,利用互相关处理器的输出来检测是否存在被动目标,并根据不同角度上检测概率的分布粗略估计目标的方位。有益效果是:本发明充分利用了振速和振速梯度的指向性,以及噪声及目标信号振速与振速梯度信息相关性的差异性特点,降低了被动目标的可检测信噪比,从而有效的提高被动目标的探测距离。本发明方法对水下弱目标的远程探测与参数估计均有重要意义,可以在浮标、潜标、鱼探仪、水下航行器、水下定位与导航等产品中应用。
申请公布号 CN106066472A 申请公布日期 2016.11.02
申请号 CN201510606683.X 申请日期 2015.09.22
申请人 青岛理工大学 发明人 王绪虎;张群飞;吉爱国;吴则举
分类号 G01S5/18(2006.01)I 主分类号 G01S5/18(2006.01)I
代理机构 代理人
主权项 一种二维振速梯度水听器的被动目标相关检测方法,包括如下步骤:步骤一:将二维振速梯度水听器八个通道的接收数据v<sub>xi</sub>、v<sub>yi</sub>(i=1,2,3,4),通过希尔伯特变换器转换成复数信号<img file="FDA0000807964860000011.GIF" wi="183" he="81" />(i=1,2,3,4);步骤二:对四个通道的x轴复振速信号<img file="FDA0000807964860000012.GIF" wi="66" he="71" />(i=1,2,3,4)求平均,求取振速梯度水听器中心点处x轴方向振速<img file="FDA0000807964860000013.GIF" wi="94" he="77" />步骤三:对四个通道的y轴复振速信号<img file="FDA0000807964860000014.GIF" wi="69" he="79" />(i=1,2,3,4)求平均,求取振速梯度水听器中心点处y轴方向振速<img file="FDA0000807964860000015.GIF" wi="96" he="81" />步骤四:利用八通道振速复信号<img file="FDA0000807964860000016.GIF" wi="183" he="82" />(i=1,2,3,4)求振速梯度水听器中心点处横向振速梯度<img file="FDA0000807964860000017.GIF" wi="366" he="87" />和纵向振速梯度<img file="FDA0000807964860000018.GIF" wi="391" he="86" />步骤五:利用引导方向α,对振速梯度水听器中心点处振速分量<img file="FDA0000807964860000019.GIF" wi="186" he="80" />进行加权处理,得到合成振速<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>v</mi><mo>~</mo></mover><mi>c</mi></msub><mrow><mo>(</mo><mi>&alpha;</mi><mo>)</mo></mrow><mo>=</mo><msub><mover><mi>v</mi><mo>~</mo></mover><mrow><mi>o</mi><mi>x</mi></mrow></msub><mi>c</mi><mi>o</mi><mi>s</mi><mi>&alpha;</mi><mo>+</mo><msub><mover><mi>v</mi><mo>~</mo></mover><mrow><mi>o</mi><mi>y</mi></mrow></msub><mi>s</mi><mi>i</mi><mi>n</mi><mi>&alpha;</mi><mo>;</mo></mrow>]]></math><img file="FDA00008079648600000110.GIF" wi="587" he="90" /></maths>步骤六:利用引导方向α,对振速梯度水听器中心点处横向振速梯度<img file="FDA00008079648600000111.GIF" wi="190" he="86" /><img file="FDA00008079648600000112.GIF" wi="163" he="90" />和纵向振速梯度<img file="FDA00008079648600000113.GIF" wi="394" he="87" />进行加权处理,得到合成振速梯度<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>d</mi><mo>~</mo></mover><mi>c</mi></msub><mrow><mo>(</mo><mi>&alpha;</mi><mo>)</mo></mrow><mo>=</mo><mo>&part;</mo><msub><mover><mi>v</mi><mo>~</mo></mover><mi>x</mi></msub><mo>/</mo><mo>&part;</mo><msup><mi>xcos</mi><mn>2</mn></msup><mi>&alpha;</mi><mo>+</mo><mo>&part;</mo><msub><mover><mi>v</mi><mo>~</mo></mover><mi>y</mi></msub><mo>/</mo><mo>&part;</mo><mi>y</mi><mi> </mi><msup><mi>sin</mi><mn>2</mn></msup><mi>&alpha;</mi><mo>+</mo><mrow><mo>(</mo><mo>&part;</mo><msub><mover><mi>v</mi><mo>~</mo></mover><mi>x</mi></msub><mo>/</mo><mo>&part;</mo><mi>y</mi><mo>+</mo><mo>&part;</mo><msub><mover><mi>v</mi><mo>~</mo></mover><mi>y</mi></msub><mo>/</mo><mo>&part;</mo><mi>x</mi><mo>)</mo></mrow><mi>c</mi><mi>o</mi><mi>s</mi><mi>&alpha;</mi><mi>s</mi><mi>i</mi><mi>n</mi><mi>&alpha;</mi><mo>;</mo></mrow>]]></math><img file="FDA00008079648600000114.GIF" wi="1433" he="98" /></maths>步骤七:对M个采样点的合成振速<img file="FDA00008079648600000115.GIF" wi="138" he="88" />与合成振速梯度<img file="FDA00008079648600000116.GIF" wi="152" he="88" />进行相关处理,得到被动目标的检测统计量<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>D</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>,</mo><mi>M</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mo>|</mo><mo>|</mo><msub><mover><mi>v</mi><mo>~</mo></mover><mi>c</mi></msub><mrow><mo>(</mo><mi>&alpha;</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mover><mi>d</mi><mo>~</mo></mover><mi>c</mi></msub><mrow><mo>(</mo><mi>&alpha;</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>;</mo></mrow>]]></math><img file="FDA00008079648600000117.GIF" wi="725" he="144" /></maths>步骤八:根据系统要求的最大虚警概率P<sub>f</sub>和上述检验统计量D(α,M),在噪声环境中计算出不同引导方位上的检测门限D<sub>th</sub>(α,M);步骤九:实时计算D(α,M),将D(α,M)与D<sub>th</sub>(α,M)进行比较,如果D(α,M)>D<sub>th</sub>(α,M),则目标存在,反之目标不存在,并在不同的α上进行扫描,根据成功检测概率的分布粗略估计目标的方位。
地址 266520 山东省青岛市黄岛区嘉陵江路777号