主权项 |
一种水下远场目标被动定位过程中的近场强干扰源抑制方法,其特征是:(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>θ</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>Bα</mi><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow></mrow><msqrt><msup><mi>α</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><mi>Bα</mi><mrow><mo>(</mo><mi>θ</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>θ</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mi>B</mi><mi>H</mi></msup><msup><mi>α</mi><mi>H</mi></msup><mi>RBα</mi><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow></mrow><msqrt><msup><mi>α</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><mi>Bα</mi><mrow><mo>(</mo><mi>θ</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>θ</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>α</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>θ</mi><mo>)</mo></mrow></mrow><msqrt><mi>α</mi><msup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>θ</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>α</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>θ</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>θ</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>B</mi><msup><mi>α</mi><mi>H</mi></msup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>θ</mi><mo>)</mo></mrow><mi>RBα</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>θ</mi><mo>)</mo></mrow></mrow><mrow><msup><mi>α</mi><mi>H</mi></msup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>θ</mi><mo>)</mo></mrow><mi>Bα</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mi>θ</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>θ</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><msub><mi>η</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>θ</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000560627360000021.GIF" wi="452" he="140" /></maths>η<sub>i</sub>为各子带能量加权系数,取η<sub>i</sub>=1。 |