发明名称 水下远场目标被动定位过程中的近场强干扰源抑制方法
摘要 本发明提供的是一种水下远场目标被动定位过程中的近场强干扰源抑制方法。声纳设备对远场目标探测前,利用近场声图扫描的方法提取声纳设备近场区域中干扰源的个数N及位置信息;建立水下远场目标探测时存在近场强干扰源的模型,计算出第i个干扰源所在位置坐标的导向矢量;确定波束零陷权和方位谱的形式;判断接收基阵接收到的信号是否属于宽带信号;将接收到的宽带信号在频域上划分成若干窄带信号,最后综合所有频段的输出得到宽带信号的估计结果。本发明将本地干扰视为近场,针对其所在位置坐标形成聚焦到点的波束指向性零点,给出了具体的波束零陷权和方位谱的形式,克服了以往波束零陷方法带来的“盲区”影响;同时又将真实目标视为远场进行测量。
申请公布号 CN102809744B 申请公布日期 2015.02.25
申请号 CN201210149999.7 申请日期 2012.05.15
申请人 哈尔滨工程大学 发明人 张亮;梅继丹;王逸林;蔡平
分类号 G01S7/537(2006.01)I;G01S15/06(2006.01)I 主分类号 G01S7/537(2006.01)I
代理机构 代理人
主权项 一种水下远场目标被动定位过程中的近场强干扰源抑制方法,其特征是:(1)声纳设备对远场目标探测前,利用近场声图扫描的方法提取声纳设备近场区域中干扰源的个数N及位置信息r<sub>i</sub>、θ<sub>i</sub>,i=1,…,N;其中r<sub>i</sub>为第i个干扰源与接收基阵中心的距离,θ<sub>i</sub>为第i个干扰源与接收基阵中心连线和基阵左侧延长线的夹角;(2)建立水下远场目标探测时存在近场强干扰源的模型,明确近场干扰源、远场信号与接收基阵的几何关系;近场干扰源按球面波扩展,远场信号按平面波传播;以接收基阵最左端的基元作为参考基元,按照干扰源的传播规律和干扰源与接收基阵的几何关系计算出第i个干扰源所在位置坐标的导向矢量α(r<sub>i</sub>,θ<sub>i</sub>);(3)根据步骤(1)和(2)中获取的近场强干扰源的个数信息、与声纳设备的几何关系,确定波束零陷权和方位谱的形式<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>W</mi><mi>opt</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>B&alpha;</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow><msqrt><msup><mi>&alpha;</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>B&alpha;</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></msqrt></mfrac></mrow>]]></math><img file="FDA0000560627360000011.GIF" wi="519" he="155" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>P</mi><mi>opt</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mi>B</mi><mi>H</mi></msup><msup><mi>&alpha;</mi><mi>H</mi></msup><mi>RB&alpha;</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow><msqrt><msup><mi>&alpha;</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>B&alpha;</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></msqrt></mfrac></mrow>]]></math><img file="FDA0000560627360000012.GIF" wi="554" he="161" /></maths>其中,B=I‑H(H<sup>H</sup>H)<sup>‑1</sup>H<sup>H</sup>,H=[α(r<sub>1</sub>,θ<sub>1</sub>),α(r<sub>2</sub>,θ<sub>2</sub>),…,α(r<sub>N</sub>,θ<sub>N</sub>)],I为M×M的单位矩阵,M为接收基阵阵元的个数,θ为目标方位,扫描范围[0,180°],R=E{XX<sup>H</sup>},X为接收基阵的观测数据矢量;(4)判断接收基阵接收到的信号是否属于宽带信号;(5)基阵接收到的信号为宽带时,将接收到的宽带信号在频域上划分成若干窄带信号,每个子带上波束零陷权和方位谱的形式分别为<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>W</mi><mi>opt</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>[</mo><mi>I</mi><mo>-</mo><mi>H</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><msup><mrow><mo>(</mo><msup><mi>H</mi><mi>H</mi></msup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>H</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>H</mi><mi>H</mi></msup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>]</mo><mi>&alpha;</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow><msqrt><mi>&alpha;</mi><msup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>H</mi></msup><mo>[</mo><mi>I</mi><mo>-</mo><mi>H</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><msup><mrow><mo>(</mo><msup><mi>H</mi><mi>H</mi></msup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>H</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>H</mi><mi>H</mi></msup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>]</mo><mi>&alpha;</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow></msqrt></mfrac></mrow>]]></math><img file="FDA0000560627360000013.GIF" wi="1337" he="165" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>P</mi><mi>opt</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>B</mi><msup><mi>&alpha;</mi><mi>H</mi></msup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>RB&alpha;</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow><mrow><msup><mi>&alpha;</mi><mi>H</mi></msup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>B&alpha;</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000560627360000014.GIF" wi="701" he="149" /></maths>其中,H=[α(f<sub>i</sub>,r<sub>1</sub>,θ<sub>1</sub>),α(f<sub>i</sub>,r<sub>2</sub>,θ<sub>2</sub>),…,α(f<sub>i</sub>,r<sub>N</sub>,θ<sub>N</sub>)],信号第i个频段的方向矢量为<img file="FDA0000560627360000015.GIF" wi="759" he="79" />f<sub>i</sub>为第i个频段的中心频率,τ=dsin(θ)/C,C为声速,d为接收基阵相邻阵元间的距离;最后综合所有频段的输出得到宽带信号的估计结果<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><msub><mi>&eta;</mi><mi>i</mi></msub><msub><mi>P</mi><mi>opt</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000560627360000021.GIF" wi="452" he="140" /></maths>η<sub>i</sub>为各子带能量加权系数,取η<sub>i</sub>=1。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室
您可能感兴趣的专利