1.外辐射源雷达中一种基于随机有限集的多目标跟踪方法，其特征在于，包括以下步骤：
步骤1，利用外辐射源雷达接收数字电视信号，对外辐射源雷达接收的数字电视信号依次进行解调、重编码，重构出各个发射站的直达波信号；
步骤2，利用每个发射站的直达波信号对外辐射源雷达接收的数字电视信号进行自适应杂波相消处理，得到杂波相消后信号；根据所述杂波相消后信号、以及重构出的各个发射站的直达波信号，得到各个目标的观测量；
步骤3，采用概率假设密度粒子滤波方法对外辐射源雷达矩形观测区域内的每个目标实施定位跟踪，每个目标的定位跟踪方法为到达时间定位方法。
2.如权利要求1所述的外辐射源雷达中一种基于随机有限集的多目标跟踪方法，其特征在于，所述步骤2的具体子步骤为：
(2.1)利用每个发射站的直达波信号对外辐射源雷达接收的数字电视信号进行自适应杂波相消处理，得到杂波相消后信号；
(2.2)对所述杂波相消后信号、以及重构出的各个发射站的直达波信号进行距离—多普勒二维相关运算，得到各个目标的观测量；每个目标的观测量包括：对应目标的双基地距离和以及对应目标的速度。
3.如权利要求1所述的外辐射源雷达中一种基于随机有限集的多目标跟踪方法，其特征在于，所述步骤3具体包括以下子步骤：
(3.1)设置参数k，k＝0,1,2...；当k＝0时，在k时刻用L0个粒子表示目标状态的先验分布，所述L0个粒子表示为 表示0时刻所述L0
个粒子中第i0个粒子的权值， 表示所述L0个粒子中第i0个粒子在0时刻的粒子状态；
(3.2)当k≥1时，将k-1时刻的粒子集为 Lk-1为k-1时刻的粒子
数，ik-1取1至Lk-1， 为k-1时刻第ik-1个粒子的权值， 为k-1时刻第ik-1个粒子的粒子状态；在k-1时刻每个粒子处产生一组西格玛点，k-1时刻每个粒子处产生的一组西格玛点包括 至 以及 至 n为粒子状态的维数；
将k-1时刻每个粒子处产生的一组西格玛点带入目标状态转移方程中，得到k时刻的一步预测状态值 j取0至2n；对 至 进行加权求和，得到k
时刻粒子状态的初步一步预测值 计算出k时刻 的初步一步预测状态协方差矩阵
将 代入目标观测方程中，得到k时刻 的一步预测观测值 所述目
标观测方程为zk＝h(xk,νk)，其中，νk为观测噪声，h(xk,νk)表示目标的观测函数；对至 进行加权求和，得到k时刻观测量的初步一步预测值
在k时刻目标的观测集Zk中，选取出与 具有最近欧式距离的对应目标观测量zk，k时刻目标的观测集Zk为：由k时刻每个目标的观测量组成的集合；用zk对粒子状态的初步一步预测值 和 的初步一步预测状态协方差矩阵 做进一步修正，得到初步一步预测均值 及 的初步一步预测状态协方差矩阵 在均值为 协
方差为 的高斯分布上采样，得到k时刻第ik-1个粒子的粒子状态的一步预测值根据 的概率密度qk，并按照下式计算出 所对应的权值
其中，ps,k为k-1时刻的目标在k时刻仍然存在的概率, 为目标状
态转移概率密度函数；
(3.3)当k≥1时，得出k时刻每个发射站对应的每个目标的双基地双基地椭圆，在k时刻每个发射站对应的每个目标的双基地椭圆曲线上，每个点对应的双基地距离和为一固定值，每个点对应的双基地距离和指：每个点与对应发射站的距离、以及对应点与外辐射源雷达的距离之和；根据k时刻每个发射站对应的每个目标的双基地椭圆曲线与外辐射源雷达的矩形观测区域边缘的交点，分布k时刻新生粒子，并计算k时刻每个新生粒子的权值，k时刻新生粒子的个数表示为Jk，k时刻第inew个新生粒子的粒子状态表示为 k时刻第inew个新生粒子的权值表示为 inew取Lk-1+1至Lk-1+Jk；
(3.4)当k≥1时，根据子步骤(3.2)和子步骤(3.3)，得到预测的k时刻的粒子集其中， 表示预测的k时刻的粒子集中第i'k个粒子的粒子状态，
表示预测的k时刻的粒子集中第i'k个粒子的权值，i'k取1至Lk-1+Jk；利用k时刻目标的观测集Zk，对预测的k时刻的粒子集中第i'k个粒子的权值 进行修正，得出预测的k时刻的粒子集中第i'k个粒子的修正后权值
(3.5)按照预测的k时刻的粒子集中第i'k个粒子的修正后权值 对预测的k时刻的粒子集进行重采样，得出k时刻的粒子集；对k时刻的粒子集进行聚类分析，得到每个目标k时刻的状态。
4.如权利要求3所述的外辐射源雷达中一种基于随机有限集的多目标跟踪方法，其特征在于，在子步骤(3.3)中，计算k时刻每个发射站对应的每个目标的双基地椭圆曲线与外辐射源雷达的矩形观测区域边缘的交点的位置；根据k时刻每个发射站对应的每个目标的双基地椭圆曲线与外辐射源雷达的矩形观测区域边缘的交点的位置，在外辐射源雷达的矩形观测区域边缘中的上边缘、下边缘、左边缘或右边缘中，找出任一个满足以下条件的边缘：与k时刻每个发射站对应的每个目标的双基地椭圆曲线存在交点；将找出的满足以上条件的边缘记为：k时刻目标的新生粒子分布边缘；
在k时刻目标的新生粒子分布边缘、以及k时刻每个发射站对应的每个目标的双基地椭圆曲线的所有交点中，找出相距最近的l个交点，这l个交点与l个发射站一一对应；得出所述相距最近的l个交点的位置的均值，将所述相距最近的l个交点的位置的均值记为：
k时刻共同交点；根据目标观测方程，计算k时刻共同交点对应的目标速度；
在k时刻共同交点处采样Jk个点，得出k时刻Jk个新生粒子、k时刻Jk个新生粒子的粒子状态；按下式计算k时刻第inew个新生粒子的权值
其中，inew取Lk-1+1至Lk-1+Jk， 为k时刻第inew个新生粒子的位置在x方向上的概率密度函数， 为k时刻第inew个新生粒子的位置在y方向上的概率密度函数， 为k时刻第inew个新生粒子的速度在x方向上的概率密度函数， 为k时刻第inew个新生粒子的速度在y方向上的概率密度函数；当k时刻目标的新生粒子分布边缘为外辐射源雷达的矩形观测区域边缘中的上边缘或下边缘时，当k时刻目标的新生粒子分布边缘为外辐射源雷达的矩
形观测区域边缘中的左边缘或右边缘时，
5.如权利要求3所述的外辐射源雷达中一种基于随机有限集的多目标跟踪方法，其特征在于，在子步骤(3.4)中，首先由下式计算 处的检测概率
其中， 表示 处的信噪比，pFA为设定的虚警概率，Q[·]表示marcum Q
函数， I0表示零阶贝塞尔函数；
然后，按照下式计算k时刻杂波平均数λk和杂波概率密度ck：
其中，L1为外辐射源雷达接收信号的距离范围，L2为外辐射源雷达接收信号的多普勒频率范围，l1外辐射源雷达接收信号的距离分辨率，l2为外辐射源雷达接收信号的多普勒分辨率；
利用zk计算似然函数
其中，Rν是观测噪声νk的协方差矩阵，上标‐1表示矩阵的逆，上标T表示矩阵或向量的转置；
按照下式计算出预测的k时刻的粒子集中第i'k个粒子的修正后权值
其中，i'k取1至Lk-1+Jk。
6.如权利要求3所述的外辐射源雷达中一种基于随机有限集的多目标跟踪方法，其特征在于，在子步骤(3.5)中，首先根据预测的k时刻的粒子集中第i'k个粒子的修正后权值计算k时刻的估计目标数
针对每个目标固定采样m个粒子，对预测的k时刻的粒子集进行重采样，得到k时刻的粒子集 ik取1至Lk， 表示k时刻的粒子集中第ik个粒
子的粒子状态， 表示k时刻的粒子集中第ik个粒子权值。
7.如权利要求3所述的外辐射源雷达中一种基于随机有限集的多目标跟踪方法，其特征在于，在步骤3之后，针对每个目标0时刻至k时刻的状态，进行航迹关联，剔除虚假航迹，得出每个目标的航迹。