一种用于空时自适应处理雷达的检测前跟踪方法
【专利摘要】该发明提供了一种用于空时自适应处理雷达的检测前跟踪方法，属于雷达目标跟踪【技术领域】，特别涉及了空时自适应处理(STAP)雷达微弱目标检测跟踪【技术领域】。首先利用全部帧全部距离单元计算统计STAP协方差矩阵；然后计算每帧每个距离单元在计算统计STAP协方差矩阵时所占的权重；进而利用每一帧中小于该帧权重均值的权重对应的距离单元重新估计统计STAP协方差矩阵；再计算每帧每个距离单元对应的检测统计量，并形成新数据平面；最后对新数据平面新距离单元进行值函数更新，最终得到目标航迹。从而具有减少估计性能损失，从而达到减少计算量的同时提高检测跟踪准确性、实时性的效果。
【专利说明】一种用于空时自适应处理雷达的检测前跟踪方法
【技术领域】
[0001] 本发明属于雷达目标跟踪【技术领域】，特别涉及了空时自适应处理（STAP)雷达微 弱目标检测跟踪【技术领域】。
【背景技术】
[0002] 微弱目标跟踪由于目标能量微弱、信噪比比较低，通常采用检测前跟踪方法。相对 于传统的目标跟踪方法，检测前跟踪方法并不在单帧内宣布目标检测结果，而是利用目标 在帧与帧之间的信息将多帧数据联合处理，充分积累目标在各帧之间的信息，提高信噪比， 并在宣布检测结果的同时获得目标航迹。动态规划检测前跟踪算法是检测前跟踪方法中的 一种。它是穷尽搜索法的一种等效实现算法，但是它的计算效率远远高于穷尽搜索。它的 主要思想是将n维优化问题转化为n个1维优化问题，对优化问题的分级处理大大降低了 它的计算量。
[0003] Danilo Orlando在2010年首次将动态规划检测前跟踪算法应用到STAP雷达中， 见文献 "Orlando, D.，et al. Track-before-detect strategies for STAP radars. IEEE Transactions on Signal Processing, 58, 2 (Feb. 2010), 933一938"。该算法首先利用回波 数据估计统计STAP协方差矩阵；然后对每个距离单元计算基于随机检验（ad hoc test)的 检测统计量，形成新的回波数据平面；对新的数据平面利用动态规划检测前跟踪算法进行 目标检测和跟踪。该算法针对时变场景（统计STAP协方差矩阵固定）和稳定场景（统计 STAP协方差矩阵不固定）采用两种不同的统计STAP协方差矩阵方法，结果显示场景信息对 算法性能具有重要的影响，若场景信息与真实场景不匹配，算法性能严重下降。另外，该算 法在时变场景下所采用的统计STAP协方差矩阵方法所需估计与求逆的统计STAP协方差矩 阵个数为距离单元数与积累帧数的乘积，而矩阵求逆计算量大、耗时，所以该算法不利于工 程实现。
【发明内容】
[0004] 针对【背景技术】存在的缺陷，本发明提出了一种用于空时自适应处理雷达的动态规 划检测前跟踪方法，该方法采用新的统计STAP协方差矩阵估计方法，减少估计性能损失， 从而达到减少计算量的同时提高检测跟踪准确性、实时性的目的。
[0005] 本发明的技术方案为一种用于空时自适应处理雷达的检测前跟踪方法，该方法包 括：
[0006] 步骤1 :系统参数初始化，初始化参数包括阵元（雷达传感器）个数Na，脉冲数目 Np，距离单元个数L，动态规划检测前跟踪算法处理帧数M，自然数r ;其中对于一个指定的 方位（方位角为fc)，雷达所有Na个阵元接收到的Np次脉冲回波数据称为一帧数据，将一 帧数据中各阵元接收到的每次脉冲回波数据平均分为若干段，并对各阵元接收到的每次 脉冲回波数据的分段顺序编号，将一帧数据中相同编号的数据段排列在一起称为一个距离 单元；
[0007] 步骤2 :估计统计STAP协方差矩阵
[0008] 步骤2. 1 :根据接收到的M帧回波数据，计算统计STAP协方差矩阵
【权利要求】
1. 一种用于空时自适应处理雷达的检测前跟踪方法，该方法包括： 步骤1 :系统参数初始化，初始化参数包括阵元（雷达传感器）个数Na，脉冲数目Np， 距离单元个数L，动态规划检测前跟踪算法处理帧数M，自然数r ;其中对于一个指定的方位 (方位角为氣），雷达所有Na个阵元接收到的Np次脉冲回波数据称为一帧数据，将一帧数 据中各阵元接收到的每次脉冲回波数据平均分为若干段，并对各阵元接收到的每次脉冲回 波数据的分段顺序编号，将一帧数据中相同编号的数据段排列在一起称为一个距离单元； 步骤2 :估计统计STAP协方差矩阵 步骤2. 1 :根据接收到的M帧回波数据，计算统计STAP协方差矩阵
其中zf表示第m帧第1个距离单元，其结构如下：
其中< gi)表示第m帧，第1个距离单元，第i个阵元，第j个脉冲对应的回波数据； 步骤2. 2 :计算第m帧第1个距离单元数据在计算协方差矩阵时所占的权重
步骤2. 3 :利用每一帧中小于该帧权重均值的权重对应的距离单元重新估计统计STAP 协方差矩阵，具体操作如下：对m = 1，2, ...，M，假设if, 1 = 1，2,…，L中小于权重均值的 权重为
，则利用权重
对应的距离单元重新估计统 计STAP协方差矩阵，即
步骤3 :计算每帧中每个距离单元基于ad hoc检验的检测统计量，并针对每一帧形成 新的LXl维回波数据平面； 步骤4 :利用动态规划检测前跟踪方法处理新的数据平面 步骤4. 1 :值函数初始化 对于第1帧新的数据平面，即k = 1，每个新距离单元的值函数为该新距离单元记录的 数据； 步骤4. 2 :值函数更新 对于第2帧的新数据平面的每个新距离单元，如果该新距离单元包含目标，计算目标 在第1帧新数据平面中可能所在的所有新距离单元，寻找出这些新距离单元中最大值函数 所在的位置，并将此最大值函数与第2帧当前新距离单元中记录的数据进行叠加，将叠加 之后的值赋值给第2帧当前新距离单元的值函数； 步骤4. 3 :采用与步骤4. 2相同的方法计算出第3帧到第K帧的值函数； 步骤4. 4 :将第K帧的最大值函数与门限Vt进行比较，若最大值函数高于门限，则认为 目标存在，其中门限Vt通过蒙特卡洛仿真得到； 步骤4. 5 :将步骤4. 4中认定的目标，利用步骤4. 1-4. 3中记录的目标在前K帧的目标 位置信息，恢复出目标航迹； 步骤4. 6、若目标存在，则宣布目标存在，并且输出目标航迹。
2.如权利要求1所述的一种用于空时自适应处理雷达的检测前跟踪方法，其特征在于 步骤3的具体步骤为： 步骤3. 1 :对于每帧中每个距离单元计算基于ad hoc检验的检测统计量
表不时间导向 矢量；
表示空间导向矢量；?为Kronecker乘积符 号
，d表不雷达传感器之间的间隔，λ表不波长；
，T表不脉冲 重复周期，f。表示载波频率，V表示目标的径向速度，c表示电磁波在介质中传播的速度； u (vm，vsm)H表示向量u (vm，vsm)的共轭转置；IT1表示步骤2中估计的统计STAP协方差矩阵 R的逆矩阵；max表示求最大值；I I表示取模运算
表示不模 糊多普勒频域，其中r为自然数（通常选择r时应使得rNp为2的指数次幂）； 步骤3. 2针对每一帧回波数据形成新数据平面 每帧新数据平面大小为LX 1维，并将第1个单元称之为新距离单元，第1个新距离单 元记录的数据为步骤3. 1中计算的该帧第1个距离单元基于ad hoc检验的检测统计量。
【文档编号】G01S13/66GK104391293SQ201410628005
【公开日】2015年3月4日 申请日期:2014年11月7日 优先权日:2014年11月7日
【发明者】易伟, 姜海超, 李小龙, 姬亚龙, 孔令讲, 崔国龙, 杨晓波, 杨益川 申请人:电子科技大学