主权项 |
一种双基地多输入多输出雷达多目标跟踪方法,包括如下步骤:(1)建立信号模型:利用空间中雷达阵列天线的方向矢量,建立双基地多输入多输出雷达多目标信号模型,得到发射阵列天线和接收阵列天线的阵列响应矩阵以及雷达天线接收端的接收数据矩阵;(2)初始化:(2a)将目标的波离方向和波达方向的先验信息,作为待跟踪目标的位置预估计值;(2b)将待跟踪目标的位置预估计值,代入发射阵列天线和接收阵列天线的阵列响应矩阵,计算得到发射阵列天线和接收阵列天线的阵列响应矩阵初始值;(3)匹配滤波:(3a)脉冲信号发生器发射探测脉冲信号,将探测脉冲信号的样本数目作为采样数目,对雷达天线接收端接收到的信号进行采样,得到接收数据;(3b)对接收数据进行匹配滤波,得到匹配滤波后的接收数据矩阵;(4)更新阵列响应矩阵:(4a)利用匹配滤波后的接收数据矩阵,按照下式,计算阵列响应矩阵的估计误差:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>J</mi><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msup><mi>β</mi><mrow><mi>p</mi><mo>-</mo><mi>n</mi></mrow></msup><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>R</mi><mi>n</mi></msub><mo>-</mo><msub><mi>BΓ</mi><mi>n</mi></msub><msup><mi>A</mi><mi>H</mi></msup><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000686075960000011.GIF" wi="732" he="155" /></maths>其中,J表示阵列响应矩阵的估计误差,Σ表示作求和操作,p表示当前脉冲信号发生器发射的探测脉冲信号个数,n表示脉冲信号发生器发射的每一个探测脉冲信号个数,n的取值范围为n=1,...,p,β表示遗忘因子,β的取值范围为0<β≤1,R<sub>n</sub>表示第n个脉冲信号发生器发射的探测脉冲信号期间匹配滤波后的接收数据矩阵,B表示接收阵列天线阵列响应矩阵,Γ<sub>n</sub>表示第n个脉冲信号发生器发射的探测脉冲信号期间由待跟踪目标的雷达截面积和多普勒频移所确定的对角矩阵,A表示发射阵列天线阵列响应矩阵,||·||<sup>2</sup>表示取范数的平方操作,(·)<sup>H</sup>表示复共轭操作;(4b)按照下式,计算前一个脉冲信号发生器发射的探测脉冲信号所对应的虚拟阵列响应矩阵:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>X</mi><mrow><mo>(</mo><mi>p</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msup><mi>A</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>p</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>⊗</mo><mi>B</mi><mrow><mo>(</mo><mi>p</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000686075960000021.GIF" wi="693" he="82" /></maths>其中,X(p‑1)表示第p‑1个脉冲信号发生器发射的探测脉冲信号所对应的虚拟阵列响应矩阵,A(p‑1)表示第p‑1个脉冲信号发生器发射的探测脉冲信号期间发射阵列天线阵列响应矩阵,B(p‑1)表示第p‑1个脉冲信号发生器发射的探测脉冲信号时接收阵列天线阵列响应矩阵,<img file="FDA0000686075960000025.GIF" wi="169" he="74" />表示作克罗内克乘积运算,(·)<sup>*</sup>表示共轭操作;(4c)将匹配滤波后的数据矩阵列向量化,得到匹配滤波后的接收数据矩阵列向量化后的列向量;(4d)按照下式,计算当前脉冲信号发生器发射的探测脉冲信号期间:<img file="FDA0000686075960000022.GIF" wi="374" he="82" />其中,d<sub>p</sub>表示第p个脉冲信号发生器发射的探测脉冲信号期间,由待跟踪目标的雷达截面积和多普勒频移所确定的对角矩阵的对角元素形成的列向量,X(p‑1)表示第p‑1个脉冲信号发生器发射的探测脉冲信号所对应的虚拟阵列响应矩阵,r<sub>p</sub>表示第p个脉冲信号发生器发射的探测脉冲信号期间匹配滤波后的接收数据矩阵列向量化后的列向量,<img file="FDA0000686075960000023.GIF" wi="81" he="79" />表示作矩阵的伪逆操作;(4e)按照下式,计算当前脉冲信号发生器发射的探测脉冲信号期间:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>D</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mi>βD</mi><mrow><mo>(</mo><mi>p</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>d</mi><mi>p</mi></msub><msubsup><mi>d</mi><mi>p</mi><mi>H</mi></msubsup></mrow>]]></math><img file="FDA0000686075960000024.GIF" wi="561" he="77" /></maths>其中,D(p)表示第p个脉冲信号发生器发射的探测脉冲信号期间,由待跟踪目标的雷达截面积和多普勒频移所确定的对角矩阵的对角元素,形成的列向量的自相关矩阵,β表示遗忘因子,β的取值范围为0<β≤1,D(p‑1)表示第p‑1个脉冲信号发生器发射的探测脉冲信号期间,由对角矩阵的对角元素构成的列向量的自相关矩阵,d<sub>p</sub>表示由第p个脉冲信号发生器发射的探测脉冲信号期间,由对角矩阵的对角元素形成的列向量,(·)<sup>H</sup>表示复共轭操作;(4f)按照下式,计算当前脉冲信号发生器发射的探测脉冲信号期间:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>Q</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mi>βQ</mi><mrow><mo>(</mo><mi>p</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>r</mi><mi>p</mi></msub><msubsup><mi>d</mi><mi>p</mi><mi>H</mi></msubsup></mrow>]]></math><img file="FDA0000686075960000031.GIF" wi="548" he="84" /></maths>其中,Q(p)表示第p个脉冲信号发生器发射的探测脉冲信号期间,匹配滤波后的接收数据矩阵列向量化后的列向量,与由待跟踪目标的雷达截面积和多普勒频移所确定的对角矩阵的对角元素形成的列向量,二者的互相关矩阵,β表示遗忘因子,β的取值范围为0<β≤1,Q(p‑1)表示由第p‑1个脉冲信号发生器发射的探测脉冲信号期间,匹配滤波后的接收数据矩阵列向量化后的列向量,与由对角矩阵的对角元素形成的列向量的互相关矩阵,r<sub>p</sub>表示第p个脉冲信号发生器发射的探测脉冲信号期间,匹配滤波后的接收数据矩阵列向量化后的列向量,d<sub>p</sub>表示由第p个脉冲信号发生器发射的探测脉冲信号期间,由对角矩阵的对角元素形成的列向量,(·)<sup>H</sup>表示复共轭操作;(4g)按照下式,计算当前脉冲信号发生器发射的探测脉冲信号所对应的虚拟阵列响应矩阵:<img file="FDA0000686075960000032.GIF" wi="396" he="76" />其中,X(p)表示第p个脉冲信号发生器发射的探测脉冲信号所对应的虚拟阵列响应矩阵,Q(p)表示第p个脉冲信号发生器发射的探测脉冲信号期间,匹配滤波后的接收数据矩阵列向量化后的列向量,与由待跟踪目标的雷达截面积和多普勒频移所确定的对角矩阵的对角元素形成的列向量,二者的互相关矩阵,D(p)表示第p个脉冲信号发生器发射的探测脉冲信号期间,由待跟踪目标的雷达截面积和多普勒频移所确定的对角矩阵的对角元素,形成的列向量的自相关矩阵,<img file="FDA0000686075960000033.GIF" wi="68" he="74" />表示矩阵的伪逆操作;(4h)将当前脉冲信号发生器发射的探测脉冲信号所对应的虚拟阵列响应矩阵的每一列的元素,转换成M×N维的矩阵,M和N分别表示发射阵列和接收阵列的阵元数目;(4i)按照下式,计算当前脉冲信号发生器发射的探测脉冲信号期间估计得到的接收阵列响应矩阵的列向量:<img file="FDA0000686075960000041.GIF" wi="463" he="149" />其中,<img file="FDA0000686075960000042.GIF" wi="144" he="76" />表示第p个脉冲信号发生器发射的探测脉冲信号期间,待跟踪的第i个目标相对于接收阵列天线的方向向量,<img file="FDA0000686075960000043.GIF" wi="58" he="61" />表示待跟踪的第i个目标相对于接收阵列参考阵元的波达方向,i表示待跟踪目标标号,G<sub>i</sub>表示由第p个脉冲信号发生器发射的探测脉冲信号所对应的虚拟阵列响应矩阵的第i列的元素,转换成的M×N维的矩阵,a<sub>p‑1</sub>(<sub>i</sub>)表示第p‑1个脉冲信号发生器发射的探测脉冲信号期间,待跟踪的第i个目标相对于发射阵列天线的方向向量,θ<sub>i</sub>表示待跟踪的第i个目标相对于发射阵列参考阵元的波离方向,||·||<sup>2</sup>表示作取范数的平方操作;(4j)按照下式,计算当前脉冲信号发生器发射的探测脉冲信号期间估计得到的发射阵列响应矩阵的列向量:<img file="FDA0000686075960000044.GIF" wi="466" he="153" />其中,a<sub>p</sub>(θ<sub>i</sub>)表示第p个脉冲信号发生器发射的探测脉冲信号期间,待跟踪的第i个目标相对于发射阵列天线的方向向量,θ<sub>i</sub>表示待跟踪的第i个目标相对于发射阵列参考阵元的波离方向,i表示待跟踪目标标号,G<sub>i</sub>表示由第p个脉冲信号发生器发射的探测脉冲信号所对应的虚拟阵列响应矩阵的第i列的元素,转换成的M×N维的矩阵,<img file="FDA0000686075960000045.GIF" wi="174" he="77" />表示第p‑1个脉冲信号发生器发射的探测脉冲信号期间,待跟踪的第i个目标相对于接收阵列天线的方向向量,<img file="FDA0000686075960000046.GIF" wi="53" he="56" />表示待跟踪的第i个目标相对于接收阵列参考阵元的波达方向,||·||<sup>2</sup>表示作取范数的平方操作;(5)确定目标当前位置:(5a)采用逆映射方法,使用当前脉冲信号发生器发射的探测脉冲信号对应的发射阵列天线的阵列响应矩阵,得到待跟踪目标波离方向;(5b)采用逆映射方法,使用当前脉冲信号发生器发射的探测脉冲信号对应的接收阵列天线的阵列响应矩阵,得到待跟踪目标波达方向;(5c)将待跟踪目标的波离方向和波达方向作为待跟踪目标相对于天线阵列的当前位置;(6)判断是否继续跟踪目标,若是,执行步骤(3),否则,执行步骤(7);(7)目标跟踪结束。 |