本发明涉及数据处理技术领域，具体涉及一种基于指数平滑预测的dp-tbd算法跟踪方法及设备。
背景技术：
现有技术中，很多场景需要进行检测跟踪，例如飞行器航迹检测跟踪。检测前跟踪算法是区别于传统检测后跟踪算法的一种新的目标检测跟踪思路，是目前低信噪比下检测目标的一种有效的方法。要实现这种方法通常有两类方式：一类是基于给定模型预测的算法，如卡尔曼滤波、粒子滤波等一系列算法；另一类是无运动模型的优化算法，如动态规划等。
对于当前机动能力强的飞行器来说，第一类算法应用因为受到运动模型匹配性的限制失去了其原有的优越性，第二类算法能较好的适应需求。但是传统dp-tbd算法仅使用当前帧的观测数据对指标函数进行更新积累，容易受到尖峰噪声的干扰，造成航迹关联错误。
技术实现要素：
本发明的目的在于针对现有技术中传统dp-tbd算法仅使用当前帧的观测数据对指标函数进行更新积累，容易受到尖峰噪声的干扰，造成航迹关联错误的不足之处，提供一种基于指数平滑预测的dp-tbd算法跟踪方法及设备，以解决上述不足之处。
为实现上述目的，本发明提供一种基于指数平滑预测的dp-tbd算法跟踪方法，所述基于指数平滑预测的dp-tbd算法跟踪方法包括以下步骤：
对第k帧每一个设定目标状态，定义积累指标函数的初始值，定义回溯函数的初始值，并定义平滑初始值，其中k为大于等于1的整数；
确定第2帧的关联区域，根据所述关联区域，更新所述积累函数与所述回溯函数，根据更新后的积累函数和更新后的回溯函数，求解第2帧的平滑值；
确定第k帧的关联区域，根据所述第k帧的关联区域，计算第k帧的状态转移密度函数，根据所述第k帧的状态转移密度函数，更新第k帧的积累函数与回溯函数，根据更新的第k帧的积累函数与回溯函数，求解第k帧平滑值；
根据递归积累的结果，进行门限检测；
根据门限检测的结果，进行航迹回溯。
优选地，所述对第k帧每一个设定目标状态，定义积累指标函数的初始值，定义回溯函数的初始值，并定义平滑初始值，包括：
对于第k帧每一个假设目标状态xk，对第1帧各个假设目标状态所对应的分辨单元内的观测数据z1(x1)的λ1加权值v1(x1)＝λ1z1(x1)，获得积累指标函数的初始值v1(x1)；
定义回溯函数为ψk(xk)，用以记录与当前帧当前状态相关性最大的前一帧的目标状态，并定义ψ1(x1)初始值为(0,0)；
定义平滑初始值为
优选地，所述确定第2帧的关联区域，根据所述关联区域，更新所述积累函数与所述回溯函数，根据更新后的积累函数和更新后的回溯函数，求解第2帧的平滑值，包括：
第2帧所有分辨单元内的目标状态x2状态转移集合为
其中，t为雷达扫描周期，
根据所述关联区域，按照第一公式和第二公式更新所述积累指标函数和所述回溯函数；
所述第一公式为
所述第二公式为
其中，p(x2|x1)是状态转移集合t2(x2)上的均匀分布；
根据第三公式、第四公式和第五公式对第2帧的状态值进行指数平滑预测，获得第2帧的平滑值；
所述第三公式为
所述第四公式为
所述第五公式为
优选地，所述确定第k帧的关联区域，根据所述第k帧的关联区域，计算第k帧的状态转移密度函数，根据所述第k帧的状态转移密度函数，更新第k帧的积累函数与回溯函数，根据更新的第k帧的积累函数与回溯函数，求解第k帧平滑值，包括：
当3≤k≤k时，确定第k帧中所有分辨单元内的目标状态转移集合
根据已获得的第k-1帧平滑值，利用预测模型得到第k帧的预测
按第六公式和第七公式更新积累指标函数和回溯函数，所述第六公式
利用第八公式、第九公式和第十公式，计算获得第k帧各假设目标状态的平滑值；
所述第八公式为
所述第九公式为
所述第十公式为
优选地，所述根据递归积累的结果，进行门限检测，包括：
当k＝k时，对vk(xk)做门限检测；
当vk(xk)≥vt时，认定xk置于集合c；
当vk(xk)＜vt时，认定目标不存在；
其中，vt为设定的门限值，c代表所有候选航迹终止位置的集合。
优选地，所述根据门限检测的结果，进行航迹回溯，包括：
对所有的
所述第十一公式为
此外，为实现上述目的，本发明还提出一种基于指数平滑预测的dp-tbd算法跟踪装置，所述基于指数平滑预测的dp-tbd算法跟踪装置包括：
初始模块，用于对第k帧每一个设定目标状态，定义积累指标函数的初始值，定义回溯函数的初始值，并定义平滑初始值，其中k为大于等于1的整数；
起始积累模块，用于确定第2帧的关联区域，根据所述关联区域，更新所述积累函数与所述回溯函数，根据更新后的积累函数和更新后的回溯函数，求解第2帧的平滑值；
递归积累模块，用于确定第k帧的关联区域，根据所述第k帧的关联区域，计算第k帧的状态转移密度函数，根据所述第k帧的状态转移密度函数，更新第k帧的积累函数与回溯函数，根据更新的第k帧的积累函数与回溯函数，求解第k帧平滑值；
检测模块，用于根据递归积累的结果，进行门限检测；
航迹回溯模块，用于根据门限检测的结果，进行航迹回溯。
此外，为实现上述目的，本发明还提出一种基于指数平滑预测的dp-tbd算法跟踪设备，所述基于指数平滑预测的dp-tbd算法跟踪设备包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的基于指数平滑预测的dp-tbd算法跟踪程序，所述基于指数平滑预测的dp-tbd算法跟踪程序被所述处理器执行时实现如上文所述的基于指数平滑预测的dp-tbd算法跟踪方法的步骤。
此外，为实现上述目的，本发明还提出一种存储介质，所述存储介质上存储有基于指数平滑预测的dp-tbd算法跟踪程序，所述基于指数平滑预测的dp-tbd算法跟踪程序被处理器执行时实现如上文所述的基于指数平滑预测的dp-tbd算法跟踪方法的步骤。
本发明中，对第k帧每一个设定目标状态，定义积累指标函数的初始值，定义回溯函数的初始值，并定义平滑初始值，确定第2帧的关联区域，根据所述关联区域，更新所述积累函数与所述回溯函数，根据更新后的积累函数和更新后的回溯函数，求解第2帧的平滑值，确定第k帧的关联区域，根据所述第k帧的关联区域，计算第k帧的状态转移密度函数，根据所述第k帧的状态转移密度函数，更新第k帧的积累函数与回溯函数，根据更新的第k帧的积累函数与回溯函数，求解第k帧平滑值，根据递归积累的结果，进行门限检测，根据门限检测的结果，进行航迹回溯。本发明提出的基于测量值所在状态与目标真实状态的偏差角度的转移概率模型，由于在动态规划的积累过程中，当前帧的目标真实状态是未知的，因此以指数平滑预测状态来代替当前帧的真实状态进行计算。采用的指数平滑预测法利用了所有历史帧的数据信息，既增强了目标的帧间关联性，又具有较高的预测精度，能较为准确地计算状态转移概率从而提高tbd算法的性能，从而提高航迹回溯的准确性。
附图说明
图1是本发明实施例方案涉及的硬件运行环境的基于指数平滑预测的dp-tbd算法跟踪设备的结构示意图；
图2为本发明基于指数平滑预测的dp-tbd算法跟踪方法实施例的流程示意图；
图3为本发明基于指数平滑预测的dp-tbd算法跟踪装置实施例的结构框图。
本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
参照图1，图1为本发明实施例方案涉及的硬件运行环境的基于指数平滑预测的dp-tbd算法跟踪设备结构示意图。
如图1所示，该基于指数平滑预测的dp-tbd算法跟踪设备可以包括：处理器1001，例如中央处理器(centralprocessingunit，cpu)，通信总线1002、用户接口1003，网络接口1004，存储器1005。其中，通信总线1002用于实现这些组件之间的连接通信。用户接口1003可以包括显示屏(display)，可选用户接口1003还可以包括标准的有线接口、无线接口，对于用户接口1003的有线接口在本发明中可为usb接口。网络接口1004可选的可以包括标准的有线接口、无线接口(如无线保真(wireless-fidelity，wi-fi)接口)。存储器1005可以是高速的随机存取存储器(randomaccessmemory，ram)存储器，也可以是稳定的存储器(non-volatilememory，nvm)，例如磁盘存储器。存储器1005可选的还可以是独立于前述处理器1001的存储装置。
本领域技术人员可以理解，图1中示出的结构并不构成对基于指数平滑预测的dp-tbd算法跟踪设备的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
如图1所示，作为一种计算机存储介质的存储器1005中可以包括操作系统、网络通信模块、用户接口模块以及基于指数平滑预测的dp-tbd算法跟踪程序。
在图1所示的基于指数平滑预测的dp-tbd算法跟踪设备中，网络接口1004主要用于连接后台服务器，与所述后台服务器进行数据通信；用户接口1003主要用于连接用户设备；所述基于指数平滑预测的dp-tbd算法跟踪设备通过处理器1001调用存储器1005中存储的基于指数平滑预测的dp-tbd算法跟踪程序，并执行本发明实施例提供的基于指数平滑预测的dp-tbd算法跟踪方法。
基于上述硬件结构，提出本发明基于指数平滑预测的dp-tbd算法跟踪方法的实施例。
参照图2，提出本发明基于指数平滑预测的dp-tbd算法跟踪方法实施例。
在本实施例中，所述基于指数平滑预测的dp-tbd算法跟踪方法包括以下步骤：
步骤s10：对第k帧每一个设定目标状态，定义积累指标函数的初始值，定义回溯函数的初始值，并定义平滑初始值，其中k为大于等于1的整数。
应理解的是，本实施例的执行主体是所述基于指数平滑预测的dp-tbd算法跟踪设备，其中，所述基于指数平滑预测的dp-tbd算法跟踪设备可为个人电脑或服务器等电子设备，本实施例对此不加以限制。
从二阶markov递归积累模型的推导中可以看出，目标的状态转移概率直接影响积累指标函数的取值。因此，如何准确对状态转移概率进行建模对tbd的性能十分重要。假设目标运动状态服从二阶markov模型，用来对状态转移概率进行近似，但其关联帧数较少，当目标在帧间的关联性较强时，该模型无法准确描述目标运动状态的变化特性。而且该模型对目标进行了在短时间内做直线运动假设限制，当目标机动性较强时，转移概率的近似精度下降，导致算法帧间关联错误。
一些改进算法引入了多帧数据和状态转移概率，但是关联帧数较少，不能充分描述观测数据的帧间相关性，而且转移概率模型受目标直线运动状态的限制，在目标机动性较强的情况下，算法的检测和跟踪性能下降。本实施例对tonissen提出的基于幅度积累的基于基于指数平滑预测的dp-tbd算法跟踪(dp-tbd)算法算法进行了相应的改进。提出的方法利用所有历史帧的数据信息，采用三次指数平滑预测目标的当前状态，并以此为基础建立了余弦偏差角状态转移概率模型。由该模型计算出的转移概率精确度更高，而且不受目标运动形式的限制，适用于多种运动模型，仿真结果证明，本算法具有良好的检测跟踪效果。
可理解的是，基于指数平滑预测的动态规划tbd可总结为如下的步骤：步骤一：初始化。步骤二：起始积累。步骤三：递归积累。步骤四：检测。步骤五：航迹回溯。
在具体实现中，初始化包括，对于第k帧每一个假设目标状态xk，定义积累指标函数的初始值v1(x1)是对第1帧各个假设目标状态所对应的分辨单元内的观测数据z1(x1)的λ1加权值，即v1(x1)＝λ1z1(x1)。
定义回溯函数为ψk(xk)，用以记录与当前帧当前状态相关性最大的前一帧的目标状态，并定义ψ1(x1)初始值为(0,0)，即ψ1(x1)＝(0,0)。
定义平滑初始值为
步骤s20：确定第2帧的关联区域，根据所述关联区域，更新所述积累函数与所述回溯函数，根据更新后的积累函数和更新后的回溯函数，求解第2帧的平滑值。
需要说明的是，起始积累包括：
1)确定第2帧的关联区域。
第2帧所有分辨单元内的目标状态x2状态转移集合：
t为雷达扫描周期，
2)更新积累函数与回溯函数。
按照第一公式和第二公式更新积累指标函数和回溯函数，所述第一公式为
由于在航迹起始的第一帧和第二帧，目标有可能转向任何方向，因此p(x2|x1)是状态转移集合t2(x2)上的均匀分布。
3)求解第2帧的平滑值。
根据第三公式、第四公式和第五公式对第2帧的状态值进行指数平滑预测，获得第2帧的平滑值。
所述第三公式为
所述第四公式为
所述第五公式为
步骤s30：确定第k帧的关联区域，根据所述第k帧的关联区域，计算第k帧的状态转移密度函数，根据所述第k帧的状态转移密度函数，更新第k帧的积累函数与回溯函数，根据更新的第k帧的积累函数与回溯函数，求解第k帧平滑值。
应理解的是，递归积累包括：
1)确定第k帧的关联区域。
当3≤k≤k时，确定第k帧中所有分辨单元内的目标状态转移集合
2)计算第k帧的状态转移密度函数。
根据已获得的第k-1帧平滑值，并利用预测模型得到第k帧的预测
3)更新第k帧的积累函数与回溯函数。
按第六公式和第七公式更新积累指标函数和回溯函数，所述第六公式
4)求解第k帧平滑值。
利用第八公式、第九公式和第十公式，计算获得第k帧各假设目标状态的平滑值。
所述第八公式为
所述第九公式为
所述第十公式为
步骤s40：根据递归积累的结果，进行门限检测。
可理解的是，进行检测的步骤包括：
当k＝k时，对vk(xk)做门限检测,当vk(xk)≥vt时，认定xk置于集合c；当vk(xk)＜vt时，认定目标不存在，其中，vt为设定的门限值，c代表所有候选航迹终止位置的集合。
步骤s50：根据门限检测的结果，进行航迹回溯。
需要说明的是，对所有的
输出航迹序列的估计值
在本实施例中，传统dp-tbd在低信噪比条件下，尖峰噪声容易造成帧间关联错误，而引入了目标状态转移概率模型作为指标函数的二阶markov多帧关联dp-tbd算法却受到了关联帧间直线运动条件的制约。本实施例为了消除这一约束条件，引入了指数平滑预测来构造状态转移概率模型。对两种目标运动形式的仿真结果表明，此算法在各种运动形式下都能取得良好的检测跟踪结果。利用了所有历史帧的数据信息，采用指数平滑预测推算当前帧的预测状态，构建目标真实状态偏差角度的概率转移模型，既增强了目标的帧间关联性，又不受限于目标直线运动假设，提高了机动目标的预测精度。
此外，本发明实施例还提出一种存储介质，所述存储介质上存储有基于指数平滑预测的dp-tbd算法跟踪程序，所述基于指数平滑预测的dp-tbd算法跟踪程序被处理器执行时实现如上文所述的基于指数平滑预测的dp-tbd算法跟踪方法的步骤。
此外，参照图3，本发明实施例还提出一种基于指数平滑预测的dp-tbd算法跟踪装置，所述基于指数平滑预测的dp-tbd算法跟踪装置包括：
初始模块10，用于对第k帧每一个设定目标状态，定义积累指标函数的初始值，定义回溯函数的初始值，并定义平滑初始值，其中k为大于等于1的整数；
起始积累模块20，用于确定第2帧的关联区域，根据所述关联区域，更新所述积累函数与所述回溯函数，根据更新后的积累函数和更新后的回溯函数，求解第2帧的平滑值；
递归积累模块30，用于确定第k帧的关联区域，根据所述第k帧的关联区域，计算第k帧的状态转移密度函数，根据所述第k帧的状态转移密度函数，更新第k帧的积累函数与回溯函数，根据更新的第k帧的积累函数与回溯函数，求解第k帧平滑值；
检测模块40，用于根据递归积累的结果，进行门限检测；
航迹回溯模块50，用于根据门限检测的结果，进行航迹回溯。
本发明所述基于指数平滑预测的dp-tbd算法跟踪装置的其他实施例或具体实现方式可参照上述各方法实施例，此处不再赘述。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者系统不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者系统所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者系统中还存在另外的相同要素。
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。在列举了若干装置的单元权利要求中，这些装置中的若干个可以是通过同一个硬件项来具体体现。词语第一、第二、以及第三等的使用不表示任何顺序，可将这些词语解释为标识。
以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。
1.一种基于指数平滑预测的dp-tbd算法跟踪方法，其特征在于，所述基于指数平滑预测的dp-tbd算法跟踪方法包括以下步骤：
对第k帧每一个设定目标状态，定义积累指标函数的初始值，定义回溯函数的初始值，并定义平滑初始值，其中k为大于等于1的整数；
确定第2帧的关联区域，根据所述关联区域，更新所述积累函数与所述回溯函数，根据更新后的积累函数和更新后的回溯函数，求解第2帧的平滑值；
确定第k帧的关联区域，根据所述第k帧的关联区域，计算第k帧的状态转移密度函数，根据所述第k帧的状态转移密度函数，更新第k帧的积累函数与回溯函数，根据更新的第k帧的积累函数与回溯函数，求解第k帧平滑值；
根据递归积累的结果，进行门限检测；
根据门限检测的结果，进行航迹回溯。
2.如权利要求1所述的基于指数平滑预测的dp-tbd算法跟踪方法，其特征在于，所述对第k帧每一个设定目标状态，定义积累指标函数的初始值，定义回溯函数的初始值，并定义平滑初始值，包括：
对于第k帧每一个假设目标状态xk，对第1帧各个假设目标状态所对应的分辨单元内的观测数据z1(x1)的λ1加权值v1(x1)＝λ1z1(x1)，获得积累指标函数的初始值v1(x1)；
定义回溯函数为ψk(xk)，用以记录与当前帧当前状态相关性最大的前一帧的目标状态，并定义ψ1(x1)初始值为(0,0)；
定义平滑初始值为
3.如权利要求2所述的基于指数平滑预测的dp-tbd算法跟踪方法，其特征在于，所述确定第2帧的关联区域，根据所述关联区域，更新所述积累函数与所述回溯函数，根据更新后的积累函数和更新后的回溯函数，求解第2帧的平滑值，包括：
第2帧所有分辨单元内的目标状态x2状态转移集合为
其中，t为雷达扫描周期，
根据所述关联区域，按照第一公式和第二公式更新所述积累指标函数和所述回溯函数；
所述第一公式为
所述第二公式为
其中，p(x2|x1)是状态转移集合t2(x2)上的均匀分布；
根据第三公式、第四公式和第五公式对第2帧的状态值进行指数平滑预测，获得第2帧的平滑值；
所述第三公式为
所述第四公式为
所述第五公式为
4.如权利要求3所述的基于指数平滑预测的dp-tbd算法跟踪方法，其特征在于，所述确定第k帧的关联区域，根据所述第k帧的关联区域，计算第k帧的状态转移密度函数，根据所述第k帧的状态转移密度函数，更新第k帧的积累函数与回溯函数，根据更新的第k帧的积累函数与回溯函数，求解第k帧平滑值，包括：
当3≤k≤k时，确定第k帧中所有分辨单元内的目标状态转移集合
根据已获得的第k-1帧平滑值，利用预测模型得到第k帧的预测
按第六公式和第七公式更新积累指标函数和回溯函数，所述第六公式
利用第八公式、第九公式和第十公式，计算获得第k帧各假设目标状态的平滑值；
所述第八公式为
所述第九公式为
所述第十公式为
5.如权利要求4所述的基于指数平滑预测的dp-tbd算法跟踪方法，其特征在于，所述根据递归积累的结果，进行门限检测，包括：
当k＝k时，对vk(xk)做门限检测；
当vk(xk)≥vt时，认定xk置于集合c；
当vk(xk)＜vt时，认定目标不存在；
其中，vt为设定的门限值，c代表所有候选航迹终止位置的集合。
6.如权利要求5所述的基于指数平滑预测的dp-tbd算法跟踪方法，其特征在于，所述根据门限检测的结果，进行航迹回溯，包括：
对所有的
所述第十一公式为
7.一种基于指数平滑预测的dp-tbd算法跟踪设备，其特征在于，所述基于指数平滑预测的dp-tbd算法跟踪设备包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的基于指数平滑预测的dp-tbd算法跟踪程序，所述基于指数平滑预测的dp-tbd算法跟踪程序被所述处理器执行时实现如权利要求1至6中任一项所述的基于指数平滑预测的dp-tbd算法跟踪方法的步骤。
技术总结