一种基于复伪谱的速度滤波检测前跟踪方法与流程
本发明涉及目标检测技术领域，尤其涉及一种基于复伪谱的速度滤波检测前跟踪方法。
背景技术：
现代战争要求警戒雷达、预警雷达在进行远程警戒、预警任务的时候，具备对弱目标的检测能力。然而，目标的多样化和环境的复杂化，使雷达的探测能力面临巨大的挑战。隐身技术的应用使飞行器等目标的雷达横截面积大幅度下降，目标反射能量大幅度减弱，雷达探测威力显著下降。在强杂波环境(如海杂波、地物杂波、电离层干扰等)中目标信杂比明显降低，这就需要雷达具备更强的微弱目标探测能力。此外，目标(如超音速飞机、高超声速导弹等)的运动速度大幅增加，使得雷达的预警时间急剧减少，这就需要雷达能够探测回波更弱的远距离目标。
传统的检测跟踪方法被称为检测后跟踪(detect-before-track,dbt)算法。在检测阶段，对每一帧的回波数据设置检测门限，进行门限判决，只有具有足够强度的点迹才会被保留，然后在跟踪阶段利用这些点迹估计出目标的航迹。dbt算法在信噪比较高的时候性能较好，但是对于弱目标的检测跟踪性能较差。这是因为在单帧检测阶段，反射回波较小的目标容易被丢失。尽管可以通过降低检测门限来提高弱目标的检测概率，然而目标信号被保留的同时，大量的杂波、噪声也会被保留下来，这会增加跟踪阶段数据关联的难度。
为了解决传统dbt算法在低信噪比下检测跟踪性能不佳的问题，利用检测前跟踪(track-before-detect,tbd)算法对单帧回波数据不进行门限检测，而是充分利用目标回波和杂波噪声在时间和空间上的差异性，通过对雷达数据进行多帧积累和联合处理，提高目标的信噪比。由于tbd算法摒弃了单帧门限检测，避免了单帧门限判决造成的目标信息丢失问题，因而可以有效的改善雷达系统对微弱目标的检测和跟踪性能。
现有技术中的弱目标检测前跟踪方法推导的检测策略是基于高斯噪声背景的，推导结果无法适用于瑞利等雷达的观测背景。该方法在多帧积累过程中，只考虑了观测的幅度信息。然而，由于正交双通道处理，雷达的观测通常是复数，即雷达观测不仅存在幅度信息，而且存在相位信息。而该方法没有考虑观测的相位信息，这种信息丢弃会导致算法性能的退化。
因此，针对以上不足，需要提供一种同时考虑雷达观测的幅度和相位信息的速度滤波检测前跟踪方法。
技术实现要素：
本发明要解决的技术问题在于现有的检测前跟踪方法仅考虑了雷达观测的幅度信息，针对现有技术中的缺陷，提供了一种基于复伪谱的速度滤波检测前跟踪方法。
为了解决上述技术问题，本发明提供了一种基于复伪谱的速度滤波检测前跟踪方法，包括以下步骤：
数据获取步骤、获取雷达回波的复值观测；
滤波器参数设定步骤、设置滤波器的假设速度，由匀速直线运动模型及目标在各帧出现的分辨单元位置预测得到批处理的最后一帧目标的预测位置；
复伪谱构造步骤、构造复伪谱，所述复伪谱以目标在各帧出现的分辨单元位置的复值观测为峰值，以对应预测位置为中心；
帧内相参积累步骤、将来自同一帧的复伪谱的复值相互累加，得到来自该帧的帧内积累结果；
帧间非相参积累步骤、将来自批处理中各帧的帧内积累结果的幅度值累加到最后一帧，得到帧间积累结果；
目标检测和参数估计步骤、判断是否有速度滤波器的所述帧间积累结果超过预设的检测门限，否则判断没有检测到目标，是则选择目标所在的速度滤波器的假设速度为估计的目标速度，该速度滤波器所输出的积累平面的最大峰值所在位置为估计的目标位置，再根据目标空间域和速度域波形，对目标位置和目标速度进行进一步的精确估计；
航迹回溯并输出步骤、根据估计出的目标位置和目标速度进行航迹回溯并输出。
优选地，所述滤波器参数设定步骤中通过以下公式得到批处理的最后一帧目标的预测位置：
其中，(ncr,ncθ)为目标在第k帧时出现的分辨单元位置，ncr＝1,2,…,nr，ncθ＝1,2,…,nθ，k＝1,…,k，nr和nθ分别表示观测平面每个方向的分辨单元个数，k表示一个处理批次的积累帧数，
优选地，所述复伪谱构造步骤中根据以下公式构造复伪谱：
fc(nr,nθ,zk(ncr,ncθ),npr,npθ)＝zk(ncr,ncθ)exp(-ηr(nr-npr)2-ηθ(nθ-npθ)2)
其中，zk表示雷达第k帧的复值观测，zk(ncr,ncθ)表示第k帧分辨单元位置(ncr,ncθ)的复值观测，fc表示构造的复伪谱，(ηr,ηθ)表示包络扩散程度。
优选地，所述帧内相参积累步骤中通过以下公式得到来自该帧的帧内积累结果：
其中，
优选地，所述帧间非相参积累步骤中通过以下公式得到帧间积累结果：
其中，
式中·表示取模运算，
优选地，所述目标检测和参数估计步骤中，通过以下公式得到检测门限：
其中
上式中vt表示检测门限的值，
优选地，所述航迹回溯并输出步骤中，根据估计的第k帧的目标位置和估计的目标速度，回溯出一个批处理中前k-1帧的目标位置。
优选地，多个滤波器通过cuda或者mpi编程并行处理。
实施本发明的基于复伪谱的速度滤波检测前跟踪方法，具有以下有益效果：本发明根据单元上雷达的复值观测构造复伪谱，保留了雷达观测的相位信息，用于能量积累，解决了现有速度滤波检测前跟踪算法在处理雷达观测时存在相位信息丢弃的问题；本发明将来自同一帧的复伪谱先相互累加，实现帧内的相参积累，这种帧内相参积累处理，提高了目标能量积累效率，改善了算法性能；本发明将每一帧的帧内相参积累结果的幅度值累加到最后一帧，实现帧间能量积累，可以有效地提高目标的信噪比，且积累效果随着积累帧数的增加而改善。本发明根据复高斯分布推导出恒虚警检测门限，对检测门限的选取提供理论支持。此外，本发明中多个速度滤波器通过cuda或者mpi编程并行处理，减少了运行时间，提高了系统效率。
附图说明
图1是本发明优选实施例的基于复伪谱的速度滤波检测前跟踪方法的流程图；
图2是目标回波扩散及其相位，其中，横线阴影为目标峰值位置，竖线阴影为目标能量扩散位置，φk表示第k帧目标回波的相位，
图3是不同输入信噪比下检测概率对比；
图4是不同积累帧数下检测概率对比；
图5是不同输入信噪比下位置估计性能对比；
图6是不同输入信噪比下速度估计性能对比；
图7是不同积累帧数下位置估计性能对比；
图8是不同积累帧数下速度估计性能对比；
其中，图3-图8中，cps-vf-tbd表示本发明提出的基于复伪谱的速度滤波检测前跟踪方法；aps-vf-tbd表示基于幅度伪谱的速度滤波检测前跟踪方法；am-vf-tbd表示基于幅度值的速度滤波检测前跟踪方法。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明实施例提供的一种基于复伪谱的速度滤波检测前跟踪方法，具体包括如下步骤：
s1、数据获取步骤：获取雷达回波的复值观测。
数据获取步骤中，设雷达回波数据为zk，表示雷达第k帧的复值观测矩阵，k＝1,…,k，k表示为一个处理批次的积累帧数。雷达的回波数据平面是一个nr*nθ的离散点迹平面，每一小格代表一个单位距离。
如图2所示，本发明中考虑的目标回波模型为高斯点扩散模型，即目标回波占据了多个分辨单元，公式可以表示为：
s(nr,nθ,lr,k,lθ,k)＝a·exp(jφk)exp(-ηr(nr-lr,k)2-ηθ(nθ-lθ,k)2)(1)
其中，a表示目标回波峰值强度，φk表示第k帧目标回波的相位，(ηr,ηθ)表示包络扩散程度，nr＝1,…,nr，nθ＝1,…,nθ，(lr,k,lθ,k)表示目标第k帧的真实位置。
s2、滤波器参数设定步骤：设置滤波器的假设速度，由匀速直线运动模型及目标在各帧出现的分辨单元位置预测得到批处理的最后一帧目标的预测位置。
优选地，本发明中假设目标的速度在积累的时间内是恒定不变的，并且空间分布是恒定的。本发明考虑目标可能出现在量测平面的任意位置。本发明猜想目标在第k帧时出现在分辨单元位置(ncr,ncθ)上，ncr＝1,2,…,nr，ncθ＝1,2,…,nθ，k＝1,…,k，nr和nθ分别表示观测平面每个方向的分辨单元个数，根据匀速直线运动的模型，将该分辨单元预测到批处理的最后一帧，得到预测位置：
其中，
优选地，本发明中的多个滤波器通过cuda或者mpi编程并行处理。本发明中多个速度滤波器通过cuda或者mpi编程并行处理，能够减少运行的时间，提高了系统效率。
s3、复伪谱构造步骤：构造复伪谱，所述复伪谱以目标在各帧出现的分辨单元位置的复值观测为峰值，以对应预测位置为中心。
优选地，本发明中以第k帧分辨单元位置(ncr,ncθ)的复值观测zk(ncr,ncθ)为峰值，以对应预测位置(npr,npθ)为中心，根据目标回波的点扩散函数，构造复伪谱，具体根据以下公式得到复伪谱：
fc(nr,nθ,zk(ncr,ncθ),npr,npθ)＝zk(ncr,ncθ)exp(-ηr(nr-npr)2-ηθ(nθ-npθ)2)(3)
其中，fc表示构造出的复伪谱，(ηr,ηθ)表示包络扩散程度。
传统的速度滤波检测前跟踪方法只利用观测的幅度信息进行多帧能量积累，然而，由于雷达的正交双通道处理，所得到的观测往往是复数的，即同时存在幅度和相位信息。由于能量积累过程中相位信息的丢弃，使得传统速度滤波检测前跟踪方法的积累效率降低，算法性能退化。本发明以单元上的复值观测为峰值，以预测位置为中心，构造复伪谱，用于能量积累，保留了观测的相位信息，能够解决现有速度滤波检测前跟踪算法在处理雷达观测时存在相位信息丢弃的问题。
s4、帧内相参积累步骤、将来自同一帧的复伪谱的复值相互累加，得到来自该帧的帧内积累结果。
优选地，本发明对于来自同一帧的复伪谱，首先将其复值相互累加，根据以下公式得到来自该帧的帧内积累结果：
其中，
本发明首先将来自同一帧的复伪谱的复值相互累加，实现帧内的相参积累，这种帧内相参积累处理，利用了雷达观测的相位信息，提高了速度滤波检测前跟踪算法多帧能量积累的效率，改善了算法的检测性能和估计精度。
s5、帧间非相参积累步骤、将来自批处理中各帧的帧内积累结果的幅度值累加到最后一帧，得到帧间积累结果。
由于目标回波的相位未知，且随着帧数的变化而变化，本发明通过将步骤s4中帧内积累结果的幅度值累加到批处理的最后一帧，实现帧间的能量积累。优选地，本发明中根据以下公式得到帧间积累结果：
其中，
式中·表示取模运算，
本发明将每一帧的帧内相参积累结果的幅度值累加到最后一帧，实现帧间能量积累，可以有效地提高目标的信噪比，且积累效果随着积累帧数的增加而改善。
s6、目标检测和参数估计步骤、判断是否有速度滤波器的所述帧间积累结果超过预设的检测门限，否则判断没有检测到目标，是则选择目标所在的速度滤波器的假设速度为估计的目标速度，该速度滤波器所输出的积累平面的最大峰值所在位置为估计的目标位置，再根据目标空间域和速度域波形，对目标位置和目标速度进行进一步的精确估计。
优选地，本发明假设噪声服从复高斯分布，根据上述能量积累过程，本方法的检测门限可以表示为：
其中
上式中vt表示检测门限的值，
本发明中根据复高斯分布推导出恒虚警检测门限，检测目标和估计目标参数，对检测门限的选取提供理论支持。
本发明的目标检测准则可以表示为：
其中h0和h1分别表示不存在目标和存在目标。检测到的目标所在的速度滤波器的假设速度即为该目标速度的估计，该速度滤波器所输出的积累平面的最大峰值所在位置为估计的目标位置。为了得到更精确的参数估计结果，本发明可以根据目标输出包络的形状，通过插值、拟合等手段实现对目标位置和目标速度的进一步精确估计。
s7、航迹回溯并输出步骤、根据估计出的目标位置和目标速度进行航迹回溯并输出。
优选地，本发明根据估计的第k帧的目标位置和估计的目标速度，回溯出一个批处理中前k-1帧的目标位置。
本发明还通过将所提出的方法与传统的am-vf-tbd方法及aps-vf-tbd方法进行对比实验，验证了本发明所提出的方法在检测概率、估计精度等方面的优越性。结合图3到图8可以看出，相比于传统的am-vf-tbd方法，aps-vf-tbd方法和所提出的cps-vf-tbd方法的性能更好，这是由于aps-vf-tbd和所提出的cps-vf-tbd的额外帧内积累使得扩散的目标回波能量得到充分积累，改善了输出信噪比。相比于传统的am-vf-tbd方法和aps-vf-tbd方法，所提出的cps-vf-tbd方法能够获得最好的检测概率和更好的估计精度，这是因为复值观测中相位信息的使用，使得cps-vf-tbd方法能够进行额外的帧内相参积累，相比于传统方法的非相参积累策略，所提出的cps-vf-tbd方法中的相参积累提高了能量积累的效率，改善了算法的性能。
综上所述，由于传统的速度滤波检测前跟踪算法只考虑雷达观测的幅度信息，在能量积累过程中，丢弃了雷达观测的相位信息，这会造成积累效率的退化和算法性能的下降。本发明将雷达复值观测的相位信息考虑在多帧能量积累过程中，提出一种两步积累策略的新方法，实现目标能量的充分积累，提高算法的性能。具体地，本发明根据雷达复值观测构造复伪谱，保留了观测的相位信息；本发明首先将来自同一帧的复伪谱相互累加，实现帧内相参积累，这种帧内相参积累处理，提高了目标能量积累效率，改善了算法性能；再将批处理中每一帧的帧内积累结果的幅度值累加到最后一帧，实现帧间能量积累，能够有效地提高目标的信噪比，且积累效果随着积累帧数的增加而改善。本发明根据复高斯分布推导出恒虚警检测门限，对检测门限的选取提供理论支持，此外，本发明中多个速度滤波器通过cuda或者mpi编程并行处理，减少了运行时间，提高了系统效率。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
- 还没有人留言评论。精彩留言会获得点赞！