一种用于异步多基地雷达系统的检测跟踪一体化方法与流程

本发明属于雷达信号处理技术领域，特别涉及一种用于异步多基地雷达系统的检测跟踪一体化方法，适用于实现异步多基地雷达系统对目标的检测跟踪联合处理，并能够提升对目标的跟踪精度。

背景技术：

目标跟踪在雷达及声呐系统的很多应用中扮演着重要的角色，在战场监控、空防、空中交通管制和火控等方面都起着重要作用；近年来，随着科学技术和武器装备的发展，单个雷达站已难以应对越来越复杂的作战环境，使用组网雷达系统可有效弥补单个雷达站对目标探测跟踪的不足；目前，在多基地雷达系统信息融合理论中，研究较多的是同步多基地雷达系统目标跟踪问题；然而，实际中经常遇到异步多基地雷达系统目标跟踪问题，因为各部雷达可能具有不同的采样频率，预处理时间和传输时延等。

同时，从目标跟踪的角度来看，传统的检测跟踪算法将检测和跟踪过程分离，并且使用的各个雷达检测器都采取单向、开环的信息处理方式，各个雷达检测器在完成跟踪波门内回波数据的检测后，将过门限的点迹信息传送给跟踪器进行数据关联与融合；但是，传统的检测跟踪算法会丢失掉跟踪器输出的先验信息。

技术实现要素：

针对上述现有技术存在的不足，本发明的目的在于提出一种用于异步多基地雷达系统的检测跟踪一体化方法，该种用于异步多基地雷达系统的检测跟踪一体化方法用于在异步多基地雷达系统的条件下，实现检测跟踪联合处理，能够增强异步多基地雷达系统的检测跟踪性能。

本发明的基本思路：首先建立异步观测模型和目标运动模型目标的运动模型以及多基地雷达系统，异步多基地雷达系统中每个雷达站都结合从融合中心得到的先验信息，指导其设置自身的检测门限，从而提升目标检测概率，进而提升异步多基地雷达系统的目标跟踪性能。

为达到上述技术目的，本发明采用如下技术方案予以实现。

一种用于异步多基地雷达系统的检测跟踪一体化方法，包括以下步骤：

步骤1，建立异步多基地雷达系统，所述异步多基地雷达系统包含n个雷达站和一个融合中心，n个雷达站的探测区域中存在目标；n个雷达站向其探测区域内的目标发射信号并接收回波数据；n为大于0的正整数；

初始化：令k表示第k时刻，k的初始值为1，k∈{1,2,…,k}，k为预先设定的最大跟踪时刻，k为大于0的正整数；

步骤2，n个雷达站将接收到的回波数据发送至融合中心进行融合处理，融合中心根据所述回波数据构建目标量测模型，确定第k时刻融合中心的目标状态，得到第k时刻异步多基地雷达系统中第l个雷达站的非线性测量函数；l∈{1,2,…,n}；

步骤3，根据第k时刻融合中心的目标状态，确定目标的运动方程；

步骤4，根据第k时刻异步多基地雷达系统中第l个雷达站的非线性测量函数和目标的运动方程，计算得到第k时刻异步多基地雷达系统中第l个雷达站的检测门限；

步骤5，根据第k时刻异步多基地雷达系统中第l个雷达站的检测门限，得到第时刻异步多基地雷达系统得到的所有有效量测值；

步骤6，根据第时刻异步多基地雷达系统得到的所有有效量测值，计算第k时刻融合中心得到的目标状态估计值xk|k和第k时刻融合中心得到的目标状态估计值的协方差ck|k；

步骤7，令加1，返回步骤2，直到得到第1时刻融合中心得到的目标状态估计值x1|1至第k时刻融合中心得到的目标状态估计值xk|k，以及第1时刻融合中心得到的目标状态估计值的协方差c1|1至第k时刻融合中心得到的目标状态估计值的协方差ck|k，并记为一种用于异步多基地雷达系统的目标检测跟踪一体化结果。

本发明与现有技术相比具有如下优点：

第一，由于本发明的公式是在多基地雷达系统异步观测条件下对目标进行的推导，所以可以用于异步组网雷达系统对目标的跟踪。

第二，由于本发明通过将跟踪器当前的跟踪状态反馈到检测中心，计算目标检测雷达的恒虚警检测门限值，克服现有技术无法根据目标跟踪的位置来自适应调整波门内各检测单元的门限，使得本发明设计的在波门内平均虚警率恒定的前提下，提升目标的平均检测概率和系统的跟踪性能。个雷达站的发射信号功率和发射信号带宽，使得本发明方法能够提高组网雷达系统对目标的跟踪精度。

附图说明

下面结合附图和具体实施方式对本发明作进一步详细说明。

图1为本发明的一种用于异步多基地雷达系统的检测跟踪一体化方法流程图；

图2为异步多基地雷达系统的采样示意图；

图3为多基地雷达系统中各个雷达站和目标的空间关系图；

图4为本发明方法与未加反馈的异步组网雷达系统对目标跟踪的性能对比图；

图5为本发明方法相对于未加反馈的异步组网雷达系统对目标跟踪的性能提升率示意图。

具体实施方式

参照图1，为本发明的一种用于异步多基地雷达系统的目标检测跟踪一体化方法流程图；其中所述用于异步多基地雷达系统的目标检测跟踪一体化方法，包括以下步骤：

步骤1，建立异步多基地雷达系统，所述异步多基地雷达系统包含n个雷达站和一个融合中心，n为大于1的正整数，本实施例中n取值为6。

异步多基地雷达系统中n个雷达站的探测区域内存在目标；以目标0时刻的位置为原点o、东西方向为轴、南北方向为轴，建立平面直角坐标系时刻目标位置为异步多基地雷达系统中第i个雷达站的坐标为i＝1,2,...,n，表示异步多基地雷达系统中第i个雷达站在轴方向的位置，表示异步多基地雷达系统中第i个雷达站在轴方向的位置。

n个雷达站向其探测区域内的目标发射信号并接收回波数据；初始化：将0时刻的目标状态向量预测协方差矩阵记为c0，c0为设定的维对角阵，为每个时刻的目标状态向量维数，为大于0的正整数，本实施例中取4；其中目标状态包括任意时刻目标在轴方向的坐标值、任意时刻目标在轴方向的速度、任意时刻目标在轴方向的坐标值和任意时刻目标在轴方向的速度；k表示第k时刻，k的初始值为1，k∈{1,2,…,k}，k为预先设定的最大跟踪时刻，k为大于0的正整数，本实施例中，该最大跟踪时刻k取值为60。

步骤2，设定异步多基地雷达系统的异步观测模型，参考图2，由于异步多基地雷达系统中每个雷达站采样周期不同，并且每个雷达站分别周期性地对目标进行测量，本实施例中6个雷达站的采样周期依次为2s、4s、6s、2s、4s和6s；同时为了同步各个雷达站的时间，将多基地雷达系统中n个雷达站的采样时刻对应投影到融合中心上，令k表示融合中心第k时刻，t(k)表示融合中心第k时刻对应的时间；由于各个雷达站采样周期不同，所以在第k时刻，多基地雷达系统中对目标进行观测的雷达站个数会有不同，用nk表示第k时刻对目标进行观测的雷达站个数，用uk表示第k时刻对目标观测的雷达站标号集合，uk中元素个数为例如如果n取值为6，k时刻对目标进行观测的雷达站为异步多基地雷达系统中第1个雷达站、第3个雷达站、第5个雷达站以及第6个雷达站，即共有4个雷达站进行观测，那么的取值为4，

n个雷达站将接收到的回波数据发送至融合中心进行融合处理，融合中心根据所述回波数据按照下式，构建目标量测模型：

其中，表示融合中心第时刻异步多基地雷达系统中第l个雷达站对目标测量得到过门限γl,k的第j个含噪声的有效量测值，设定γl,k为第k时刻异步多基地雷达系统中第l个雷达站的检测门限，l取值为中的元素，∈表示属于；设定第时刻异步多基地雷达系统中第l个雷达站得到过门限的含噪声的有效量测值个数j取值为hl(xk)表示第k时刻异步多基地雷达系统中第l个雷达站的非线性测量函数，wl,k表示第k时刻异步多基地雷达系统中第l个雷达站的测量误差，υl,k表示第k时刻异步多基地雷达系统中第l个雷达站测得的虚假量测值，具体为融合中心第时刻异步多基地雷达系统中第l个雷达站对目标测量得到的除了来源于目标量测值外的剩余所有量测值；xk表示第k时刻融合中心的目标状态，其定义表达式为：

其中，xt(k)表示t(k)时刻目标在轴方向的位置，表示t(k)时刻目标在轴方向的速度，yt(k)表示t(k)时刻目标在轴方向的位置，表示t(k)时刻目标在轴方向的速度，t(k)表示融合中心第k时刻对应的时间，上标t表示转置。

hl(xk)表示第k时刻异步多基地雷达系统中第l个雷达站的非线性测量函数，其表达式为：

其中，r^l(xk)表示第k时刻异步多基地雷达系统中第l个雷达站对目标径向距离的量测函数，θ^l(xk)表示第k时刻异步多基地雷达系统中第l个雷达站对目标俯仰角度的量测函数，xl′表示异步多基地雷达系统中第l个雷达站在轴方向的位置，yl′表示异步多基地雷达系统中第l个雷达站在轴方向的位置，arctan表示求反正切。

将第k时刻异步多基地雷达系统中第l个雷达站的测量误差记为wl,k，wl,k服从零均值的高斯分布；第k时刻异步多基地雷达系统中第l个雷达站的测量误差wl,k的协方差矩阵为wl,k，其表达式为：

其中，表示第k时刻第l个雷达站对目标径向距离测量误差的方差，表示第k时刻第l个雷达站对目标俯仰角度测量误差的方差。

第k时刻第l个雷达站对目标径向距离测量误差的方差和第k时刻第l个雷达站对目标俯仰角度测量误差的方差的大小分别与第k时刻异步多基地雷达系统中第l个雷达站从目标接收的回波数据信噪比(snr)μl,k有关，其关系为：

其中，∝表示正比于，μl,k表示第k时刻异步多基地雷达系统中第l个雷达站从目标接收的回波数据信噪比(snr)，βl表示异步多基地雷达系统中第l个雷达站向目标发射的波束带宽，表示异步多基地雷达系统中第l个雷达站的3db接收波束宽度，上标-1表示求逆。

步骤3，设定异步多基地雷达系统中的目标运动为匀速直线运动，并设定目标的运动方程。

具体地，设定异步多基地雷达系统中的目标运动为匀速直线运动，并设定目标在平面直角坐标系内做匀速运动，目标的运动方程可写为：

xk＝fk-1xk-1+uk-1

其中，xk表示第k时刻融合中心的目标状态，xk-1表示第k-1时刻融合中心的目标状态，uk-1表示第k-1时刻零均值、白色高斯过程噪声，fk-1表示第k-1时刻目标状态的转移矩阵，

其中，△tk-1表示融合中心在第k时刻和第k-1时刻之间的时间差，△tk-1＝t(k)-t(k-1)，t(k)表示融合中心第k时刻对应的时间，t(k-1)表示融合中心第k-1时刻对应的时间，表示矩阵直积的运算符号，i2表示2×2维单位矩阵；将第k-1时刻零均值、白色高斯过程噪声记为uk-1，用于衡量目标状态转移的不确定性，进而得到第k-1时刻零均值、白色高斯过程噪声uk-1的协方差矩阵为qk-1：

其中，i2表示2×2维单位矩阵，表示矩阵直积的运算符号，q0表示目标运动过程的噪声强度，q0越大，目标运动的波动越大；本实施例中，q0取值为10。

t(0)表示融合中心第0时刻对应的时间，u0表示第0时刻零均值、白色高斯过程噪声，x0表示第0时刻融合中心的目标状态。

步骤4，确定检测门限值。

(4a)计算平均检测概率：

将vl,k表示第k时刻异步多基地雷达系统中第l个雷达站跟踪波门的大小，其表达式为：π表示圆周率，γ表示第k时刻异步多基地雷达系统中第l个雷达站跟踪波门系数，本实施例中γ取值为16；dl,k|k-1表示第k时刻异步多基地雷达系统中第l个雷达站的预测量测协方差，ηl,k表示第k时刻异步多基地雷达系统中第l个雷达站的雅克比矩阵，表示对xk求偏导操作，ck|k-1表示目标在第k-1时刻的预测协方差，hl(xk)表示第k时刻异步多基地雷达系统中第l个雷达站的非线性测量函数。

在第k时刻异步多基地雷达系统中第l个雷达站跟踪波门内，以第l个雷达站测量分辨率为刻度划分分辨单元，对应得到个分辨单元；其中雷达站测量分辨率包括雷达站测量距离分辨率和雷达站测量角度分辨率；本实施例中，每个雷达站测量距离分辨率均为150m，每个雷达站测量角度分辨率均为1度；表示第k时刻异步多基地雷达系统中第l个雷达站包括的分辨单元个数。

对于每个分辨单元内有无目标做出两种假设，用h0表示分辨单元内没有目标，且接收的数据仅来自于噪声；用h1表示分辨单元内存在目标，且接收的数据将是目标的回波数据和噪声的叠加；在两种假设条件下第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元接收到的回波数据幅度的概率密度函数可表示为:

其中，表示在h0假设下第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元接收到的回波数据幅度的概率密度函数，表示在h1假设下第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元接收到的回波数据幅度的概率密度函数，表示第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元接收到的回波数据幅度，表示第k时刻异步多基地雷达系统中第l个雷达站包括的分辨单元个数；μl,k表示第k时刻异步多基地雷达系统中第l个雷达站从目标接收的回波数据信噪比(snr)。

在异步多基地雷达系统检测跟踪一体化算法中，第k时刻融合中心会将目标在第k-1时刻的预测状态xk|k-1以及目标在第k-1时刻的预测协方差ck|k-1发送给第k时刻对目标观测的雷达站标号集合uk。

根据目标的运动方程，计算目标在第k时刻的预测状态xk|k-1，其表达式为：

xk|k-1＝fk-1xk-1|k-1

其中，fk-1表示第k-1时刻目标状态的转移矩阵，xk-1|k-1表示第k-1时刻融合中心得到的目标状态估计值；当k时取值为1时，x0|0表示第0时刻融合中心得到的目标状态估计值，第0时刻融合中心得到的目标状态估计值x0|0为第0时刻融合中心的目标状态x0，第0时刻融合中心的目标状态x0为目标初始状态。

按照下式计算目标在第k-1时刻的预测协方差ck|k-1：

ck|k-1＝fk-1ck-1|k-1fk-1^t+qk-1

其中，fk-1表第k-1时刻目标状态的转移矩阵，ck-1|k-1表示第k-1时刻融合中心得到的目标状态估计值的协方差，当k时取值为1时，c0|0表示第0时刻融合中心得到的目标状态估计值的协方差，第0时刻融合中心得到的目标状态估计值的协方差c0|0为0时刻的目标状态向量预测协方差矩阵c0；qk-1表示第k-1时刻零均值、白色高斯过程噪声的协方差矩阵。

异步多基地雷达系统中每个雷达站根据目标在第k时刻的预测状态xk|k-1以及目标在第k-1时刻的预测协方差ck|k-1，可以得到在h0假设下第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元对目标的量测向量概率密度函数和在h1假设下第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元对目标的量测向量概率密度函数

其中，vl,k表示第k时刻异步多基地雷达系统中第l个雷达站跟踪波门的大小，π表示圆周率，γ表示第k时刻异步多基地雷达系统中第l个雷达站跟踪波门系数，本实施例中取值为16；dl,k|k-1表示第k时刻异步多基地雷达系统中第l个雷达站的预测量测协方差，ηl,k表示第k时刻异步多基地雷达系统中第l个雷达站的雅克比矩阵，表示对xk求偏导操作，ck|k-1表示目标在第k-1时刻的预测协方差，zl,k|k-1表示第k时刻异步多基地雷达系统中第l个雷达站对目标预测状态的量测值，zl,k|k-1＝hl(xk|k-1)，hl(xk)表示第k时刻异步多基地雷达系统中第l个雷达站的非线性测量函数，hl(xk|k-1)表示第k时刻异步多基地雷达系统中第l个雷达站对目标预测状态的非线性量测函数，xk|k-1表示目标在第k时刻的预测状态，xk|k-1表示目标在第k时刻的预测状态中轴方向的位置，表示目标在第k时刻的预测状态中轴方向的速度，yk|k-1表示目标在第k时刻的预测状态中轴方向的位置，表示目标在第k时刻的预测状态中轴方向的速度，r^l(xk|k-1)表示预测第k时刻异步多基地雷达系统中第l个雷达站对目标径向距离的量测函数，θ^l(xk|k-1)表示预测第k时刻异步多基地雷达系统中第l个雷达站对目标俯仰角度的量测函数，x′l表示异步多基地雷达系统中第l个雷达站在轴方向的位置，x′l表示异步多基地雷达系统中第l个雷达站在轴方向的位置，y′l表示异步多基地雷达系统中第l个雷达站在轴方向的位置，arctan表示求反正切，上标t表示转置，表示服从正态分布的概率密度函数，表示服从以zl,k|k-1为均值、dl,k|k-1为协方差的正态分布的概率密度函数，表示第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元的量测值。

根据贝叶斯最小风险准则可得到第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元的检测门限

其中，表示用于确保第k时刻异步多基地雷达系统中第l个雷达站在恒虚警率条件下进行检测的工作常数。

各分辨单元的检测概率可通过下式积分得到第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元的检测概率其表达式为：

其中，表示在h1假设下第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元接收到的回波数据幅度的概率密度函数；由于各分辨单元检测门限不同，所以给出第k时刻异步多基地雷达系统中第l个雷达站跟踪波门内的平均检测概率其表达式为：

其中，表示在h1假设下第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元对目标的量测向量概率密度函数，n表示设定的量测向量维度，本实施例中n取值为2。

(4b)计算平均虚警概率

按照下式计算得到第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元的虚警概率其表达式为：

其中，表示在h0假设下第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元接收到的回波数据幅度的概率密度函数。

以及计算第k时刻异步多基地雷达系统中第l个雷达站跟踪波门的平均虚警概率其表达式为：

其中，表示在h0假设下第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元对目标的量测向量概率密度函数，表示对进行微分，表示第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元接收到的回波数据幅度；nz表示设定的正整数，本实施例中nz取值为2。

(4c)目标检测中通常需要保证虚警率恒定，在本实施例中，取值为10^-4，根据所述平均虚警概率计算公式，可以求得为用来确保第k时刻异步多基地雷达系统中第l个雷达站在恒虚警率条件下进行检测的工作常数进而代入第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元的检测门限公式中，得到第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第r个分辨单元的检测门限

(4d)令r的值分别取1至重复执行(4c)，进而分别得到第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第1个分辨单元的检测门限至第k时刻异步多基地雷达系统中第l个雷达站在其跟踪波门内第个分辨单元的检测门限记为第k时刻异步多基地雷达系统中第l个雷达站的检测门限γl,k。

步骤5，确定有效量测值。

根据第k时刻异步多基地雷达系统中第l个雷达站跟踪波门的大小vl,k对目标进行检测的过程中，将第k时刻异步多基地雷达系统中第l个雷达站跟踪波门内超过检测门限γl,k的量测值作为有效量测值。

设定第时刻异步多基地雷达系统中第l个雷达站得到过门限含噪声的有效量测值个数为ml,k，进而将融合中心第时刻异步多基地雷达系统中第l个雷达站得到的有效量测值表示为zl,k，其表达式为：

其中，表示融合中心第时刻异步多基地雷达系统中第l个雷达站对目标测量得到过门限的第j个含噪声的量测值，l取值为中的元素，第时刻异步多基地雷达系统中第l个雷达站得到过门限含噪声的有效量测值个数为j取值为

融合中心第时刻异步多基地雷达系统得到的所有有效量测值表示为zk：

zk＝{zl,k},l∈uk

步骤6，估计第k时刻融合中心得到的目标状态估计值xk|k和第k时刻融合中心得到的目标状态估计值的协方差ck|k。

具体地，在异步多基地雷达系统中，在每次的检测过程中，每个雷达站都可能获得多个过门限的量测值，个数设定为ml,k，ml,k个量测值中每个量测值可能来源于目标，也可能来源于虚警；使用表示融合中心第时刻的联合事件j，j表示联合事件的标号，j是由个数字依次排列而成的数组，其中每个数字表示对应量测值来源于目标或来源于虚警；nk表示第k时刻对目标进行观测的雷达站个数，表示第时刻异步多基地雷达系统中第l个雷达站得到过门限的含噪声的有效量测值个数。

考虑到单个雷达站从一个目标之多获得一个量测的限制，此处使用表示融合中心第时刻异步多基地雷达系统中第l个雷达站中第j(l′)个量测值来源于目标这一事件，j(l′)取值范围为0,1,...,ml,k，l′表示融合中心第时刻异步多基地雷达系统中第l个雷达站在uk中的序号，取值为1,2,...,nk，j(l′)＝0；则表示融合中心第时刻异步多基地雷达系统中第l个雷达站没有量测值来源于目标，所有量测值均源于虚警；为详细解释以上符号的含义，现举例说明，例如n＝4，若j＝1021则表示以下事件：多基地雷达系统中第1个雷达站的第1个量测值来源于目标，第2个雷达站的所有量测值均源于虚警，第3个雷达站的第2个量测值来源于目标，第4个雷达站的第1个量测值来源于目标。

进而得到第k时刻异步多基地雷达系统中n个雷达站联合事件的数量nk：

其中，ml,k表示第时刻异步多基地雷达系统中第l个雷达站得到过门限的含噪声的有效量测值个数。

按照下式计算联合事件的概率。

将融合中心第时刻的联合事件的联合概率定义为：

其中，表示在已知量测zk的情况下联合事件发生的概率，表示定义。

因为不同雷达站的量测之间相互独立，所以计算得到在已知量测zk的情况下联合事件发生的概率其表达式为：

其中，表示第k时刻异步多基地雷达系统中第l个雷达站的第j(l′)个有效量测值来源于目标的概率，表示融合中心第时刻异步多基地雷达系统中第l个雷达站中第j(l′)个量测值来源于目标这一事件，j(l′)取值范围为0,1,...,ml,k，l′表示融合中心第时刻异步多基地雷达系统中第l个雷达站在uk中的序号，取值为1,2,...,nk，j(l′)＝0；通过对采用概率数据关联算法(pda)计算得到；zk表示融合中心第时刻异步多基地雷达系统得到的所有有效量测值。

按照下式，计算第k时刻融合中心得到的目标状态估计值xk|k和第k时刻融合中心得到的目标状态估计值的协方差ck|k，其表达式分别为：

其中，j表示联合事件的标号，xj,k|k表示对融合中心第时刻的联合事件通过多探测器-概率数据关联算法(md-pda)计算得到的目标更新后状态，表示融合中心时间轴的第时刻的联合事件的联合概率，通过对采用概率数据关联算法(pda)计算得到；cj,k|k表示对融合中心第时刻的联合事件通过多探测器-概率数据关联算法(md-pda)得到的状态估计协方差，md-pda算法为文献“amultiple-detectionjointprobabilisticdataassociationfilter”中所提算法。

步骤7，当未达到预先设定的最大跟踪时刻k时，令加1，返回步骤2；当达到预先设定的最大跟踪时刻k时迭代重复执行过程停止，此时得到第1时刻融合中心得到的目标状态估计值x1|1至第k时刻融合中心得到的目标状态估计值xk|k，以及第1时刻融合中心得到的目标状态估计值的协方差c1|1至第k时刻融合中心得到的目标状态估计值的协方差ck|k，并记为一种用于异步多基地雷达系统的目标检测跟踪一体化结果。

通过以下仿真实验对本发明效果作进一步验证说明。

(一)仿真条件：

本发明的仿真运行系统为intel(r)core(tm)i5-4590cpu@3.30ghz，64位windows7操作系统，仿真软件采用matlab(r2014b)。

(二)仿真内容与结果分析：

参照图3本发明的仿真实验设定异步多基地雷达系统的布阵情况，异步多基地雷达系统包含的雷达站总个数n＝6,目标的初始位置在(0,0)km，并以速度为(300,20)m/s做匀速运动；仿真序列数据为60帧，每个雷达站发射信号的有效带宽为2mhz，每个雷达站发射信号波长为λc＝1m；每个雷达站天线孔径d＝10λc；相关波门系数g＝4；目标运动过程的噪声强度q0＝100；每个雷达站采样周期分别为2s、4s、6s、2s、4s和6s，目标距离72km，反射系数为1时信噪比snr设置为14db，在仿真中，反射系数均设置为1，设置波门平均虚警概率为10^-4。

参考图4给出的两种跟踪算法目标跟踪的均方根误差rmse对比，算法1为传统无反馈信息的异步多基地雷达系统跟踪算法，算法2为本发明方法；可以看到，由于目标远离雷达站飞行，所以跟踪的误差也会变大；同时可以看到跟踪性能会出现锯齿状的波动，这是因为在每个时刻对目标跟踪的雷达个数不同，当雷达站个数多的时候，跟踪性能比雷达个数少的时候好，例如图4中的第40帧对目标进行跟踪的雷达站个数比前后两帧都多，所以跟踪性能也较好。

参照图5，算法2相对于算法1性能提升率示意图，算法2相比于算法1的提升率为γ，其表达式为：

γ＝(rmse¹-rmse²)/rmse¹

rmse¹表示算法1的rmse，rmse²表示算法2的rmse；可以看到，相比算法1，算法2的跟踪精度能够提升超过5％。

综上所述，仿真实验验证了本发明的正确性，有效性和可靠性。

显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围；这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。