CN107462882A - 一种适用于闪烁噪声的多机动目标跟踪方法及系统 - Google Patents

Abstract

本发明适用于多传感器信息融合领域，提供了一种适用于闪烁噪声的多机动目标跟踪方法与系统，通过预测步骤、更新步骤、融合步骤、生成步骤以及提取步骤，用t分布建模闪烁噪声，应用变分贝叶斯方法近似地求出不同模型下观测噪声的协方差、目标状态以及自由度的联合概率密度，有效解决了闪烁噪声下多机动目标的跟踪问题，提高了多机动目标的跟踪精度。

Description

一种适用于闪烁噪声的多机动目标跟踪方法及系统

技术领域

本发明属于多传感器信息融合技术领域，尤其涉及一种适用于闪烁噪声的多机动目标跟踪方法及系统。

背景技术

闪烁噪声下的概率假设密度滤波器是适用于闪烁噪声多目标跟踪的有效方法，该滤波器能够对闪烁噪声下的非机动目标进行有效的跟踪，但该滤波器不能对在不同模型之间转换的机动目标进行跟踪。

发明内容

本发明提供了一种适用于闪烁噪声的多机动目标跟踪方法及系统，旨在解决闪烁噪声下多机动目标的有效跟踪问题。

为解决上述技术问题，本发明是这样实现的，本发明提供了一种适用于闪烁噪声的多机动目标跟踪方法，所述方法包括：

步骤1、利用前一时刻各个目标的运动模型、以及高斯伽玛混合分布和存在概率产生伽玛分布的形状参数和尺度参数；根据前一时刻各个目标的运动模型、高斯伽玛混合分布和存在概率、当前时刻与前一时刻的时间差、模型间的马尔可夫转移概率、以及伽玛分布的形状参数和尺度参数，得到当前时刻各个目标在不同运动模型下的预测高斯伽玛混合分布和预测存在概率；

步骤2、根据当前时刻各个目标在不同运动模型下的预测高斯伽玛混合分布和预测存在概率，利用变分贝叶斯方法对当前时刻的测量数据进行序贯处理，得到当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布和更新存在概率；

步骤3、对当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布和更新存在概率分别进行融合，得到当前时刻已存在目标的运动模型、高斯伽玛混合分布和存在概率；

步骤4、利用当前时刻的测量数据生成新生目标的高斯伽玛混合分布，为所述新生目标指定存在概率，并选取匀速运动模型作为新生目标的运动模型，将所述新生目标的运动模型、高斯伽玛混合分布及存在概率分别与所述当前时刻已存在目标的运动模型、高斯伽玛混合分布及存在概率进行合并，得到当前时刻各个目标的运动模型、高斯伽玛混合分布和存在概率；

步骤5、从所述当前时刻各个目标中裁减掉存在概率小于第一阈值的目标，并将裁减后余下的目标的运动模型、高斯伽玛混合分布和存在概率作为滤波器下一次递归的输入；

从所述裁减后余下的目标中提取存在概率大于第二阈值的目标，并将提取出的目标的高斯伽玛混合分布作为当前时刻的输出，将输出的高斯伽玛混合分布中，高斯分布的均值作为当前时刻目标的状态估计。

进一步地，所述步骤1包括：

以k-1表示前一时刻，k表示当前时刻，t_k-1表示前一时刻的时间，t_k表示当前时刻的时间，r_i,k-1表示前一时刻目标i的运动模型标识，且1≤r_i,k-1≤M_r，r_i,k表示当前时刻目标i的模型标识，且r_i,k＝1,2,…,M_r，其中，M_r表示运动模型的总数目；当前时刻的测量噪声服从ξ维的学生t分布，以S(y_j,k；H_kx_i,k,R_i,k,ν_i,k)表示当前时刻测量的概率密度函数，其中，y_j,k表示当前时刻的一个测量，H_k表示测量矩阵，x_i,k表示目标的状态，R_i,k表示精度矩阵，ν_i,k表示t分布的自由度，且为精度矩阵R_i,k主对角线上的各个元素；

已知k-1时刻目标i的运动模型为r_i,k-1，则在该模型r_i,k-1下的高斯伽玛混合分布为存在概率为ρ_i,k-1(r_i,k-1)；其中，N表示高斯分布，g表示伽玛分布，x_i,k-1表示k-1时刻目标i的状态向量，m_i,k-1(r_i,k-1)和P_i,k-1(r_i,k-1)分别表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中高斯分布的均值和协方差，表示k-1时刻精度矩阵的对角线元素，ν_i,k-1表示k-1时刻t分布的自由度，和γ_i,k-1(r_i,k-1)表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中伽玛分布的形状参数，和表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中伽玛分布的尺度参数，ξ表示状态向量的维数，i＝1,2,...,N_k-1，N_k-1表示前一时刻目标的总数；

根据k-1时刻目标i的运动模型、高斯伽玛混合分布和存在概率、当前时刻与前一时刻的时间差、以及模型间的马尔可夫转移概率，得到当前时刻目标i在模型r_i,k下的预测高斯伽玛混合分布为预测存在概率为ρ_i,k|k-1(r_i,k)；其中，m_i,k|k-1(r_i,k)＝F_k-1(r_i,k)m_i,k-1(r_i,k-1)， γ_i,k|k-1(r_i,k)＝ρ_γγ_i,k-1(r_i,k-1)，η_i,k|k-1(r_i,k)＝ρ_ηη_i,k-1(r_i,k-1)，ρ_i,k|k-1(r_i,k)＝p_s,k(t_k-t_k-1)t_k|k-1(r_i,k|r_i,k-1)ρ_i,k-1(r_i,k-1)，F_k-1(r_i,k)表示k-1时刻目标i在模型r_i,k下的状态转移矩阵，Q_k-1(r_i,k)表示k-1时刻目标i在模型r_i,k下的过程噪声协方差矩阵，上标T表示矩阵或向量的转置，ρ_α，ρ_β，ρ_γ，ρ_η表示传递因子，取值范围为(0,1]，p_s,k(t_k-t_k-1)表示目标的幸存概率，且t_k|k-1(r_i,k|r_i,k-1)表示模型间的马尔可夫转移概率，T表示采样周期，δ为预定的常数，l＝1,...,ξ，ξ表示状态向量的维数。

进一步地，设当前时刻的测量集为y_k＝(y_1,k,...,y_M,k)，其中，M为当前时刻的测量数据的总数，则所述步骤2包括：

以当前时刻各个目标在M_r个运动模型下的预测高斯伽玛混合分布和预测存在概率，作为当前时刻各个目标在相应运动模型下的初始高斯伽玛混合分布和初始存在概率，则目标i在模型r_i,k下的初始高斯伽玛混合分布取为

初始存在概率取为其中i＝1,2,...,N_k-1，r_i,k＝1,2,…,M_r，

利用变分贝叶斯方法对第1个至第M个当前时刻的测量数据进行序贯处理，具体包括以下：

设第j个测量数据处理前，目标i在模型r_i,k下的高斯伽玛混合分布及存在概率分别为和其中，i＝1,2,...,N_k-1,r_i,k＝1,2,…,M_r,j＝1,2,…,M；

根据变分贝叶斯方法，计算得到第j个测量更新时，目标i在模型r_i,k下的高斯伽玛混合分布为存在概率分别为其中，

；

表示伽玛函数，tr表示矩阵的迹，表示均值向量，表示协方差矩阵， 表示滤波器增益，伽玛分布的形状参数为和伽玛分布的尺度参数为和精度矩阵为 其中， H_k为模型r_i,k的观测矩阵，P_D,k为目标的检测概率，λ_c,k为杂波密度，y_j,k为当前时刻接收到的第j个测量数据，上标T表示矩阵或向量的转置，ξ为状态向量的总维数；

若则第j个测量数据处理后，目标i在模型r_i,k下的的高斯伽玛混合分布为

目标i的存在概率为其中

若则第j个测量数据处理后，目标i在模型r_i,k下的高斯伽玛混合分布为目标i的存在概率为其中，

在第M个测量数据处理后，目标i在模型r_i,k下的高斯伽玛混合分布和存在概率分别为和其中，i＝1,2,...,N_k-1，r_i,k＝1,2,…,M_r；

将第M个测量数据处理后，各个目标在不同模型下的高斯伽玛混合分布及存在概率分别作为当前时刻各个目标的更新高斯伽玛混合分布和存在概率，则所述当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布为

更新存在概率为其中i＝1,…,N_k-1，r_i,k＝1,2,…,M_r，

进一步地，所述步骤3包括：

对当前时刻目标i在M_r个运动模型下的更新的高斯伽玛混合分布和存在概率分别进行融合，得到当前时刻已存在目标q的运动模型为r_q,k、高斯伽玛混合分布为和存在概率为ρ_q,k(r_q,k)，其中，q为当前时刻已存在目标的索引号，且q＝1,…,N_k-1， 伽玛分布的形状参数为和伽玛分布的尺度参数为和

进一步地，所述步骤4包括：

选取匀速运动模型作为新生目标的运动模型，设此时新生目标的运动模型为r_j,k＝1，其中，j＝1,…,M；

利用M个测量数据生成当前时刻新生目标的高斯伽玛混合分布为其中，j＝1,…,M，为第j个新生目标的协方差，由第j个测量数据y_j，k＝[x_j，k y_j，k]^T产生，且和为第j个新生目标分布函数中伽玛分布的形状参数，和为第j个新生目标分布函数中伽玛分布的尺度参数；

指定当前时刻各新生目标的存在概率为其中，j＝1,…,M，ρ_υ为所指定的存在概率；

将所述当前时刻已存在目标的运动模型与所述当前时刻新生目标的运动模型合并，得到当前时刻各个目标的运动模型为其中，N_k＝N_k-1+M；

将所述当前时刻已存在目标的高斯伽玛混合分布与所述当前时刻新生目标的高斯伽玛混合分布进行合并，得到当前时刻各个目标的高斯伽玛混合分布为

将所述当前时刻各个目标的存在概率与所述当前时刻新生目标的存在概率进行合并，得到所述当前时刻各个目标的存在概率为

为解决上述技术问题，本发明还提供了一种适用于闪烁噪声的多机动目标跟踪系统，所述系统包括：

预测模块，用于利用前一时刻各个目标的运动模型、以及高斯伽玛混合分布和存在概率产生伽玛分布的形状参数和尺度参数；根据前一时刻各个目标的运动模型、高斯伽玛混合分布和存在概率、当前时刻与前一时刻的时间差、模型间的马尔可夫转移概率、以及伽玛分布的形状参数和尺度参数，得到当前时刻各个目标在不同运动模型下的预测高斯伽玛混合分布和预测存在概率；

更新模块，用于根据当前时刻各个目标在不同运动模型下的预测高斯伽玛混合分布和预测存在概率，利用变分贝叶斯方法对当前时刻的测量数据进行序贯处理，得到当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布和更新存在概率；

融合模块，用于对当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布和更新存在概率分别进行融合，得到当前时刻已存在目标的运动模型、高斯伽玛混合分布和存在概率；

生成模块，用于利用当前时刻的测量数据生成新生目标的高斯伽玛混合分布，为所述新生目标指定存在概率，并选取匀速运动模型作为新生目标的运动模型，将所述新生目标的运动模型、高斯伽玛混合分布及存在概率分别与所述当前时刻已存在目标的运动模型、高斯伽玛混合分布及存在概率进行合并，得到当前时刻各个目标的运动模型、高斯伽玛混合分布和存在概率；

提取模块，用于从所述当前时刻各个目标中裁减掉存在概率小于第一阈值的目标，并将裁减后余下的目标的运动模型、高斯伽玛混合分布和存在概率作为滤波器下一次递归的输入；从所述裁减后余下的目标中提取存在概率大于第二阈值的目标，并将提取出的目标的高斯伽玛混合分布作为当前时刻的输出，将输出的高斯伽玛混合分布中，高斯分布的均值作为当前时刻目标的状态估计。

本发明与现有技术相比，有益效果在于：

本发明所提供的方法通过用t分布建模闪烁噪声，应用变分贝叶斯方法近似地求出不同模型下观测噪声的协方差、目标状态以及自由度的联合概率密度，有效解决了闪烁噪声下多机动目标的跟踪问题，提高了多机动目标的跟踪精度。

附图说明

图1是本发明第一个实施例提供的一种适用于闪烁噪声的多机动目标跟踪方法示意图；

图2是本发明实施例提供的传感器经过70个扫描周期得到的测量数据示意图；

图3是在闪烁噪声下采用JMS-PHD滤波器的多目标跟踪方法处理得到的结果示意图；

图4是在闪烁噪声下采用本发明所提供的多机动目标跟踪方法处理得到的结果示意图；

图5是本发明实施例提供的按照本发明所提供的多目标跟踪方法与JMS-PHD滤波方法在经过100次实验得到的平均OSPA距离对比示意图；

图6是本发明第二个实施例提供一种适用于闪烁噪声的多目标跟踪系统的结构示意图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

作为本发明的第一个实施例，如图1所示，本发明提供的一种适用于闪烁噪声的多机动目标跟踪方法，该方法包括以下步骤：

步骤S101：在马尔科夫链控制的不同系统模型中，利用前一时刻各个目标的运动模型、以及高斯伽玛混合分布和存在概率产生伽玛分布的形状参数和尺度参数；根据前一时刻各个目标的运动模型、高斯伽玛混合分布和存在概率、当前时刻与前一时刻的时间差、模型间的马尔可夫转移概率、以及伽玛分布的形状参数和尺度参数，得到当前时刻各个目标在不同运动模型下的预测高斯伽玛混合分布和预测存在概率。步骤S101具体包括以下步骤：

以k-1表示前一时刻，k表示当前时刻，t_k-1表示前一时刻的时间，t_k表示当前时刻的时间，r_i,k-1表示前一时刻目标i的运动模型标识，且1≤r_i,k-1≤M_r，r_i,k表示当前时刻目标i的模型标识，且r_i,k＝1,2,…,M_r，其中，M_r表示运动模型的总数目；当前时刻的测量噪声服从ξ维的学生t分布，以S(y_j,k；H_kx_i,k,R_i,k,ν_i,k)表示当前时刻测量的概率密度函数，其中，y_j,k表示当前时刻的一个测量，H_k表示测量矩阵，x_i,k表示目标的状态，R_i,k表示精度矩阵，ν_i,k表示t分布的自由度，且为精度矩阵R_i,k主对角线上的各个元素；

已知k-1时刻目标i的运动模型为r_i,k-1，则在该模型r_i,k-1下的高斯伽玛混合分布为存在概率为ρ_i,k-1(r_i,k-1)；其中，N表示高斯分布，g表示伽玛分布，x_i,k-1表示k-1时刻目标i的状态向量，m_i,k-1(r_i,k-1)和P_i,k-1(r_i,k-1)分别表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中高斯分布的均值和协方差，表示k-1时刻精度矩阵的对角线元素，ν_i,k-1表示k-1时刻t分布的自由度，和γ_i,k-1(r_i,k-1)表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中伽玛分布的形状参数，和η_i,k-1(r_i,k-1)表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中伽玛分布的尺度参数，ξ表示状态向量的维数，i＝1,2,...,N_k-1，N_k-1表示前一时刻目标的总数；

根据k-1时刻目标i的运动模型、高斯伽玛混合分布和存在概率、当前时刻与前一时刻的时间差、以及模型间的马尔可夫转移概率，得到当前时刻目标i在模型r_i,k下的预测高斯伽玛混合分布为预测存在概率为ρ_i,k|k-1(r_i,k)；其中，m_i,k|k-1(r_i,k)＝F_k-1(r_i,k)m_i,k-1(r_i,k-1)， γ_i,k|k-1(r_i,k)＝ρ_γγ_i,k-1(r_i,k-1)，η_i,k|k-1(r_i,k)＝ρ_ηη_i,k-1(r_i,k-1)，ρ_i,k|k-1(r_i,k)＝p_s,k(t_k-t_k-1)t_k|k-1(r_i,k|r_i,k-1)ρ_i,k-1(r_i,k-1)，F_k-1(r_i,k)表示k-1时刻目标i在模型r_i,k下的状态转移矩阵，Q_k-1(r_i,k)表示k-1时刻目标i在模型r_i,k下的过程噪声协方差矩阵，上标T表示矩阵或向量的转置，ρ_α，ρ_β，ρ_γ，ρ_η表示传递因子，取值范围为(0,1]，p_s,k(t_k-t_k-1)表示目标的幸存概率，且t_k|k-1(r_i,k|r_i,k-1)表示模型间的马尔可夫转移概率，T表示采样周期，δ为预定的常数，l＝1,...,ξ，ξ表示状态向量的维数。

步骤S102：根据当前时刻各个目标在不同运动模型下的预测高斯伽玛混合分布和预测存在概率，利用变分贝叶斯方法对当前时刻的测量数据进行序贯处理，得到当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布和更新存在概率。设当前时刻的测量集为y_k＝(y_1,k,...,y_M,k)，其中，M为当前时刻的测量数据的总数，则此时步骤S102具体包括以下步骤：

以当前时刻各个目标在M_r个运动模型下的预测高斯伽玛混合分布和预测存在概率，作为当前时刻各个目标在相应运动模型下的初始高斯伽玛混合分布和初始存在概率，则目标i在模型r_i,k下的初始高斯伽玛混合分布取为

初始存在概率取为其中i＝1,2,...,N_k-1，r_i,k＝1,2,…,M_r，

利用变分贝叶斯方法对第1个至第M个当前时刻的测量数据进行序贯处理，具体包括以下：

设第j个测量数据处理前，目标i在模型r_i,k下的高斯伽玛混合分布及存在概率分别为和其中，i＝1,2,...,N_k-1,r_i,k＝1,2,…,M_r,j＝1,2,…,M；

根据变分贝叶斯方法，计算得到第j个测量更新时，目标i在模型r_i,k下的高斯伽玛混合分布为存在概率分别为其中，

；

表示伽玛函数，tr表示矩阵的迹，表示均值向量，表示协方差矩阵，表示滤波器增益，伽玛分布的形状参数为和伽玛分布的尺度参数为和精度矩阵为其中， H_k为模型r_i,k的观测矩阵，P_D,k为目标的检测概率，λ_c,k为杂波密度，y_j,k为当前时刻接收到的第j个测量数据，上标T表示矩阵或向量的转置，ξ为状态向量的总维数；

若则第j个测量数据处理后，目标i在模型r_i,k下的的高斯伽玛混合分布为

目标i的存在概率为其中

若则第j个测量数据处理后，目标i在模型r_i,k下的高斯伽玛混合分布为

目标i的存在概率为其中，

在第M个测量数据处理后，目标i在模型r_i,k下的高斯伽玛混合分布和存在概率分别为和其中，i＝1,2,...,N_k-1，r_i,k＝1,2,…,M_r；

将第M个测量数据处理后，各个目标在不同模型下的高斯伽玛混合分布及存在概率分别作为当前时刻各个目标的更新高斯伽玛混合分布和存在概率，则所述当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布为

更新存在概率为其中i＝1,…,N_k-1，r_i,k＝1,2,…,M_r，

步骤S103：对当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布和更新存在概率分别进行融合，得到当前时刻已存在目标的运动模型、高斯伽玛混合分布和存在概率。步骤S103具体包括以下步骤：

对当前时刻目标i在M_r个运动模型下的更新的高斯伽玛混合分布和存在概率分别进行融合，得到当前时刻已存在目标q的运动模型为r_q,k、高斯伽玛混合分布为和存在概率为ρ_q,k(r_q,k)，其中，q为当前时刻已存在目标的索引号，且q＝1,…,N_k-1， 伽玛分布的形状参数为和伽玛分布的尺度参数为和

步骤S104：利用当前时刻的测量数据生成新生目标的高斯伽玛混合分布，为所述新生目标指定存在概率，并选取匀速运动模型作为新生目标的运动模型，将所述新生目标的运动模型、高斯伽玛混合分布及存在概率分别与所述当前时刻已存在目标的运动模型、高斯伽玛混合分布及存在概率进行合并，得到当前时刻各个目标的运动模型、高斯伽玛混合分布和存在概率。步骤S104具体包括以下步骤：

选取匀速运动模型作为新生目标的运动模型，设此时新生目标的运动模型为r_j,k＝1，其中，j＝1,…,M。该匀速运动模型为多种模型中的其中一种类型(可以定义为模型1)。

利用M个测量数据生成当前时刻新生目标的高斯伽玛混合分布为其中，j＝1,…,M，为第j个新生目标的协方差，由第j个测量数据y_j，k＝[x_j，k y_j，k]^T产生，且和为第j个新生目标分布函数中伽玛分布的形状参数，和为第j个新生目标分布函数中伽玛分布的尺度参数；

指定当前时刻各新生目标的存在概率为其中，j＝1,…,M，ρ_υ为所指定的存在概率；

将所述当前时刻已存在目标的运动模型与所述当前时刻新生目标的运动模型合并，得到当前时刻各个目标的运动模型为其中，N_k＝N_k-1+M；

将所述当前时刻已存在目标的高斯伽玛混合分布与所述当前时刻新生目标的高斯伽玛混合分布进行合并，得到当前时刻各个目标的高斯伽玛混合分布为

将所述当前时刻各个目标的存在概率与所述当前时刻新生目标的存在概率进行合并，得到所述当前时刻各个目标的存在概率为

步骤S105：从所述当前时刻各个目标中裁减掉存在概率小于第一阈值的目标，并将裁减后余下的目标的运动模型、高斯伽玛混合分布和存在概率作为滤波器下一次递归的输入；从所述裁减后余下的目标中提取存在概率大于第二阈值的目标，并将提取出的目标的高斯伽玛混合分布作为当前时刻的输出，将输出的高斯伽玛混合分布中，高斯分布的均值作为当前时刻目标的状态估计。在本实施例中，经过多次的仿真实验得出，将第一阈值设定为10^-3，将第二阈值设定为0.5，可以达到较好的效果。

下面结合图2至图5对本实施例进行进一步地解释：

在本实施例中，考虑二维空间[-400(m),400(m)]×[-400(m),400(m)]中运动的机动目标。目标的状态由位置和速度构成，表示为其中x和y分别表示位置分量，和分别表示速度分量，上标T表示向量的转置，状态转移矩阵为过程噪声方差矩阵为Δt_k＝t_k-t_k-1为当前时刻与前一时刻的时间差，σ_v为过程噪声标准差；目标的运动模型包括r_i,k＝1，r_i,k＝2，r_i,k＝3，模型1是一个匀速运动模型，该模型1的状态转移矩阵和协方差矩阵分别为F_k-1(r_i,k＝1)＝F(ω＝0s^-1)和Q_k-1(r_i,k＝1)＝Q(σ_v＝1ms^-2)；模型2是一个匀转弯模型，该模型2的状态转移矩阵和协方差矩阵分别为F_k-1(r_i,k＝2)＝F(ω＝5°s^-1)和Q_k-1(r_i,k＝2)＝Q(σ_v＝3ms^-2)；模型3是一个匀转弯模型，该模型3的状态转移矩阵和协方差矩阵分别为F_k-1(r_i,k＝3)＝F(ω＝-5°s^-1)和Q_k-1(r_i,k＝3)＝Q(σ_v＝3ms^-2)；不同运动模型之间的马尔科夫转移矩阵为观测矩阵观测噪声假设服从r_k＝10的t分布。

为了产生仿真数据，取幸存概率p_S,k＝1.0、检测概率p_D,k＝0.9、杂波密度λ_c,k＝4.6875×10^-6m^-2。一次实验中传感器经过70个扫描周期得到的仿真观测数据如图2所示。本实施例在闪烁噪声下，分别采用JMS-PHD(跳变马尔科夫系统模型概率假设密度)滤波器的多目标跟踪方法和本发明所提供的适用于闪烁噪声的多机动目标跟踪方法对图2中的仿真数据分别进行了处理，在采用JMS-PHD滤波器进行处理时，相关参数设置为p_S,k＝1.0、p_D,k＝0.9、λ_c,k＝4.6875×10^-6m^-2、第一阈值为10^-3、第二阈值为0.5，传递因子ρ_α＝ρ_β＝ρ_γ＝ρ_η＝0.98，伽玛分布形状参数初始值伽玛分布尺度参数初始值在闪烁噪声下，采用JMS-PHD滤波器时的新生目标的权重为w_υ＝0.1，采用本发明实施例所提供的方法时的新生目标的存在概率为ρ_υ＝0.1、新生目标的协方差图3和图4分别为JMS-PHD滤波器方法所产生的结果和本发明实施例提供的多机动目标跟踪方法所产生的结果，通过两图的对比结果可以直观的得出，采用本发明所提供的目标跟踪方法所得到的运动轨迹更加连续、完整、精确。通过本发明实施例所提供的方法与现有的闪烁噪声下的JMS-PHD滤波器对图2的仿真数据分别进行处理，经过100次Monte Carlo实验，得到平均OSPA(Optimal Subpattern Assignment,最优亚模式分配)距离，如图5所示。通过图中对比可以看出，本发明的多目标跟踪方法在闪烁噪声下对于关联不确定、检测不确定的多机动目标跟踪可获得更为精确和可靠的目标状态估计、其滤波后的OSPA距离比闪烁噪声下的JMS-PHD滤波器滤波后的OSPA距离要小。

综上所述，本发明第一个实施例所提供的方法，根据前一时刻各个目标的运动模型、高斯伽玛混合分布和存在概率，以及模型间的马尔可夫转移概率，得到当前时刻各个目标在不同运动模型下的预测高斯伽玛混合分布及预测存在概率，根据各个目标在不同模型中预测的高斯伽玛混合分布和存在概率，利用变分贝叶斯方法序贯处理当前时刻接收到的各测量得到各目标在不同模型下的更新高斯伽玛混合分布和存在概率，对当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布和存在概率进行融合，得到当前时刻各个已存在目标的运动模型、高斯伽玛混合分布和存在概率，利用当前时刻的测量数据生成新目标的高斯伽玛混合分布，为所述新目标指定存在概率，并选取匀速运动模型作为新目标的运动模型，再将新生目标的运动模型、高斯伽玛混合分布及存在概率分别与当前时刻已存在目标的运动模型、高斯伽玛混合分布及存在概率进行合并，生成当前时刻各个目标的运动模型、高斯伽玛混合分布和存在概率，从而很好地解决了闪烁噪声下多机动目标的跟踪问题，可以达到有效跟踪的目的。

作为本发明的第二个实施例，如图6所示，本发明提供了一种适用于闪烁噪声的多机动目标跟踪系统，所述系统包括：

预测模块101，用于利用前一时刻各个目标的运动模型、以及高斯伽玛混合分布和存在概率产生伽玛分布的形状参数和尺度参数；根据前一时刻各个目标的运动模型、高斯伽玛混合分布和存在概率、当前时刻与前一时刻的时间差、模型间的马尔可夫转移概率、以及伽玛分布的形状参数和尺度参数，得到当前时刻各个目标在不同运动模型下的预测高斯伽玛混合分布和预测存在概率。预测模块101具体用于：

以k-1表示前一时刻，k表示当前时刻，t_k-1表示前一时刻的时间，t_k表示当前时刻的时间，r_i,k-1表示前一时刻目标i的运动模型标识，且1≤r_i,k-1≤M_r，r_i,k表示当前时刻目标i的模型标识，且r_i,k＝1,2,…,M_r，其中，M_r表示运动模型的总数目；当前时刻的测量噪声服从ξ维的学生t分布，以S(y_j,k；H_kx_i,k,R_i,k,ν_i,k)表示当前时刻测量的概率密度函数，其中，y_j,k表示当前时刻的一个测量，H_k表示测量矩阵，x_i,k表示目标的状态，R_i,k表示精度矩阵，ν_i,k表示t分布的自由度，且为精度矩阵R_i,k主对角线上的各个元素；

已知k-1时刻目标i的运动模型为r_i,k-1，则在该模型r_i,k-1下的高斯伽玛混合分布为存在概率为ρ_i,k-1(r_i,k-1)；其中，N表示高斯分布，g表示伽玛分布，x_i,k-1表示k-1时刻目标i的状态向量，m_i,k-1(r_i,k-1)和P_i,k-1(r_i,k-1)分别表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中高斯分布的均值和协方差，表示k-1时刻精度矩阵的对角线元素，ν_i,k-1表示k-1时刻t分布的自由度，和γ_i,k-1(r_i,k-1)表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中伽玛分布的形状参数，和表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中伽玛分布的尺度参数，ξ表示状态向量的维数，i＝1,2,...,N_k-1，N_k-1表示前一时刻目标的总数；

根据k-1时刻目标i的运动模型、高斯伽玛混合分布和存在概率、当前时刻与前一时刻的时间差、以及模型间的马尔可夫转移概率，得到当前时刻目标i在模型r_i,k下的预测高斯伽玛混合分布为预测存在概率为ρ_i,k|k-1(r_i,k)；其中，m_i,k|k-1(r_i,k)＝F_k-1(r_i,k)m_i,k-1(r_i,k-1)， γ_i,k|k-1(r_i,k)＝ρ_γγ_i,k-1(r_i,k-1)，η_i,k|k-1(r_i,k)＝ρ_ηη_i,k-1(r_i,k-1)，ρ_i,k|k-1(r_i,k)＝p_s,k(t_k-t_k-1)t_k|k-1(r_i,k|r_i,k-1)ρ_i,k-1(r_i,k-1)，F_k-1(r_i,k)表示k-1时刻目标i在模型r_i,k下的状态转移矩阵，Q_k-1(r_i,k)表示k-1时刻目标i在模型r_i,k下的过程噪声协方差矩阵，上标T表示矩阵或向量的转置，ρ_α，ρ_β，ρ_γ，ρ_η表示传递因子，取值范围为(0,1]，p_s,k(t_k-t_k-1)表示目标的幸存概率，且t_k|k-1(r_i,k|r_i,k-1)表示模型间的马尔可夫转移概率，T表示采样周期，δ为预定的常数，l＝1,...,ξ，ξ表示状态向量的维数。

更新模块102，用于根据当前时刻各个目标在不同运动模型下的预测高斯伽玛混合分布和预测存在概率，利用变分贝叶斯方法对当前时刻的测量数据进行序贯处理，得到当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布和更新存在概率。设当前时刻的测量集为y_k＝(y_1,k,...,y_M,k)，其中，M为当前时刻测量的总数，则更新模块102具体用于：

以当前时刻各个目标在M_r个运动模型下的预测高斯伽玛混合分布和预测存在概率，作为当前时刻各个目标在相应运动模型下的初始高斯伽玛混合分布和初始存在概率，则目标i在模型r_i,k下的初始高斯伽玛混合分布取为

初始存在概率取为其中i＝1,2,...,N_k-1，r_i,k＝1,2,…,M_r，

利用变分贝叶斯方法对第1个至第M个当前时刻的测量数据进行序贯处理，具体包括以下：

设第j个测量数据处理前，目标i在模型r_i,k下的高斯伽玛混合分布及存在概率分别为和其中，i＝1,2,...,N_k-1,r_i,k＝1,2,…,M_r,j＝1,2,…,M；

根据变分贝叶斯方法，计算得到第j个测量更新时，目标i在模型r_i,k下的高斯伽玛混合分布为存在概率分别为其中，

；

表示伽玛函数，tr表示矩阵的迹，表示均值向量，表示协方差矩阵， 表示滤波器增益，伽玛分布的形状参数为和伽玛分布的尺度参数为和精度矩阵为 其中， H_k为模型r_i,k的观测矩阵，P_D,k为目标的检测概率，λ_c,k为杂波密度，y_j,k为当前时刻接收到的第j个测量数据，上标T表示矩阵或向量的转置，ξ为状态向量的总维数；

若则第j个测量数据处理后，目标i在模型r_i,k下的的高斯伽玛混合分布为

目标i的存在概率为其中

若则第j个测量数据处理后，目标i在模型r_i,k下的高斯伽玛混合分布为

目标i的存在概率为其中，

在第M个测量数据处理后，目标i在模型r_i,k下的高斯伽玛混合分布和存在概率分别为和其中，i＝1,2,...,N_k-1，r_i,k＝1,2,…,M_r；

将第M个测量数据处理后，各个目标在不同模型下的高斯伽玛混合分布及存在概率分别作为当前时刻各个目标的更新高斯伽玛混合分布和存在概率，则所述当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布为

更新存在概率为其中i＝1,…,N_k-1，r_i,k＝1,2,…,M_r，

融合模块103，用于对当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布和更新存在概率分别进行融合，得到当前时刻已存在目标的运动模型、高斯伽玛混合分布和存在概率。融合模块103具体用于：

对当前时刻目标i在M_r个运动模型下的更新的高斯伽玛混合分布和存在概率分别进行融合，得到当前时刻已存在目标q的运动模型为r_q,k、高斯伽玛混合分布为和存在概率为ρ_q,k(r_q,k)，其中，q为当前时刻已存在目标的索引号，且q＝1,…,N_k-1， 伽玛分布的形状参数为和伽玛分布的尺度参数为和

生成模块104，用于利用当前时刻的测量数据生成新生目标的高斯伽玛混合分布，为所述新生目标指定存在概率，并选取匀速运动模型作为新生目标的运动模型，将所述新生目标的运动模型、高斯伽玛混合分布及存在概率分别与所述当前时刻已存在目标的运动模型、高斯伽玛混合分布及存在概率进行合并，得到当前时刻各个目标的运动模型、高斯伽玛混合分布和存在概率。生成模块104具体用于：

选取匀速运动模型作为新生目标的运动模型，设此时新生目标的运动模型为r_j,k＝1，其中，j＝1,…,M；

利用M个测量数据生成当前时刻新生目标的高斯伽玛混合分布为其中，j＝1,…,M，为第j个新生目标的协方差，由第j个测量数据y_j，k＝[x_j，k y_j，k]^T产生，且和为第j个新生目标分布函数中伽玛分布的形状参数，和为第j个新生目标分布函数中伽玛分布的尺度参数；

指定当前时刻各新生目标的存在概率为其中，j＝1,…,M，ρ_υ为所指定的存在概率；

将所述当前时刻已存在目标的运动模型与所述当前时刻新生目标的运动模型合并，得到当前时刻各个目标的运动模型为其中，N_k＝N_k-1+M；

将所述当前时刻已存在目标的高斯伽玛混合分布与所述当前时刻新生目标的高斯伽玛混合分布进行合并，得到当前时刻各个目标的高斯伽玛混合分布为

将所述当前时刻各个目标的存在概率与所述当前时刻新生目标的存在概率进行合并，得到所述当前时刻各个目标的存在概率为

提取模块105，用于从所述当前时刻各个目标中裁减掉存在概率小于第一阈值的目标，并将裁减后余下的目标的运动模型、高斯伽玛混合分布和存在概率作为滤波器下一次递归的输入；从所述裁减后余下的目标中提取存在概率大于第二阈值的目标，并将提取出的目标的高斯伽玛混合分布作为当前时刻的输出，将输出的高斯伽玛混合分布中，高斯分布的均值作为当前时刻目标的状态估计。

综上所述，本发明第二个实施例所提供的系统，通过预测模块、更新模块、融合模块、生成模块以及提取模块，用t分布建模闪烁噪声，应用变分贝叶斯方法近似地求出不同模型下观测噪声的协方差、目标状态以及自由度的联合概率密度，有效解决了闪烁噪声下多机动目标的跟踪问题，提高了多机动目标的跟踪精度。

以上所述仅为本发明的较佳实施例而已，并不用以限制发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种适用于闪烁噪声的多机动目标跟踪方法，其特征在于，所述方法包括：

步骤1、利用前一时刻各个目标的运动模型、以及高斯伽玛混合分布和存在概率产生伽玛分布的形状参数和尺度参数；根据前一时刻各个目标的运动模型、高斯伽玛混合分布和存在概率、当前时刻与前一时刻的时间差、模型间的马尔可夫转移概率、以及伽玛分布的形状参数和尺度参数，得到当前时刻各个目标在不同运动模型下的预测高斯伽玛混合分布和预测存在概率；

步骤2、根据当前时刻各个目标在不同运动模型下的预测高斯伽玛混合分布和预测存在概率，利用变分贝叶斯方法对当前时刻的测量数据进行序贯处理，得到当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布和更新存在概率；

步骤3、对当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布和更新存在概率分别进行融合，得到当前时刻已存在目标的运动模型、高斯伽玛混合分布和存在概率；

步骤4、利用当前时刻的测量数据生成新生目标的高斯伽玛混合分布，为所述新生目标指定存在概率，并选取匀速运动模型作为新生目标的运动模型，将所述新生目标的运动模型、高斯伽玛混合分布及存在概率分别与所述当前时刻已存在目标的运动模型、高斯伽玛混合分布及存在概率进行合并，得到当前时刻各个目标的运动模型、高斯伽玛混合分布和存在概率；

步骤5、从所述当前时刻各个目标中裁减掉存在概率小于第一阈值的目标，并将裁减后余下的目标的运动模型、高斯伽玛混合分布和存在概率作为滤波器下一次递归的输入；

从所述裁减后余下的目标中提取存在概率大于第二阈值的目标，并将提取出的目标的高斯伽玛混合分布作为当前时刻的输出，将输出的高斯伽玛混合分布中，高斯分布的均值作为当前时刻目标的状态估计。

2.如权利要求1所述的方法，其特征在于，所述步骤1包括：

以k-1表示前一时刻，k表示当前时刻，t_k-1表示前一时刻的时间，t_k表示当前时刻的时间，r_i,k-1表示前一时刻目标i的运动模型标识，且1≤r_i,k-1≤M_r，r_i,k表示当前时刻目标i的模型标识，且r_i,k＝1,2,…,M_r，其中，M_r表示运动模型的总数目；当前时刻的测量噪声服从ξ维的学生t分布，以S(y_j,k；H_kx_i,k,R_i,k,ν_i,k)表示当前时刻测量的概率密度函数，其中，y_j,k表示当前时刻的一个测量，H_k表示测量矩阵，x_i,k表示目标的状态，R_i,k表示精度矩阵，ν_i,k表示t分布的自由度，且 为精度矩阵R_i,k主对角线上的各个元素；

已知k-1时刻目标i的运动模型为r_i,k-1，则在该模型r_i,k-1下的高斯伽玛混合分布为存在概率为ρ_i,k-1(r_i,k-1)；其中，N表示高斯分布，g表示伽玛分布，x_i,k-1表示k-1时刻目标i的状态向量，m_i,k-1(r_i,k-1)和P_i,k-1(r_i,k-1)分别表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中高斯分布的均值和协方差，表示k-1时刻精度矩阵的对角线元素，ν_i,k-1表示k-1时刻t分布的自由度，和γ_i,k-1(r_i,k-1)表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中伽玛分布的形状参数，和η_i,k-1(r_i,k-1)表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中伽玛分布的尺度参数，ξ表示状态向量的维数，i＝1,2,...,N_k-1，N_k-1表示前一时刻目标的总数；

根据k-1时刻目标i的运动模型、高斯伽玛混合分布和存在概率、当前时刻与前一时刻的时间差、以及模型间的马尔可夫转移概率，得到当前时刻目标i在模型r_i,k下的预测高斯伽玛混合分布为预测存在概率为ρ_i,k|k-1(r_i,k)；其中，m_i,k|k-1(r_i,k)＝F_k-1(r_i,k)m_i,k-1(r_i,k-1)， γ_i,k|k-1(r_i,k)＝ρ_γγ_i,k-1(r_i,k-1)，η_i,k|k-1(r_i,k)＝ρ_ηη_i,k-1(r_i,k-1)，ρ_i,k|k-1(r_i,k)＝p_s,k(t_k-t_k-1)t_k|k-1(r_i,k|r_i,k-1)ρ_i,k-1(r_i,k-1)，F_k-1(r_i,k)表示k-1时刻目标i在模型r_i,k下的状态转移矩阵，Q_k-1(r_i,k)表示k-1时刻目标i在模型r_i,k下的过程噪声协方差矩阵，上标T表示矩阵或向量的转置，ρ_α，ρ_β，ρ_γ，ρ_η表示传递因子，取值范围为(0,1]，p_s,k(t_k-t_k-1)表示目标的幸存概率，且t_k|k-1(r_i,k|r_i,k-1)表示模型间的马尔可夫转移概率，T表示采样周期，δ为预定的常数，l＝1,...,ξ，ξ表示状态向量的维数。

3.如权利要求2所述的方法，其特征在于，设当前时刻的测量集为y_k＝(y_1,k,...,y_M,k)，其中，M为当前时刻的测量数据的总数，则所述步骤2包括：

以当前时刻各个目标在M_r个运动模型下的预测高斯伽玛混合分布和预测存在概率，作为当前时刻各个目标在相应运动模型下的初始高斯伽玛混合分布和初始存在概率，则目标i在模型r_i,k下的初始高斯伽玛混合分布取为 初始存在概率取为其中i＝1,2,...,N_k-1，r_i,k＝1,2,…,M_r，

利用变分贝叶斯方法对第1个至第M个当前时刻的测量数据进行序贯处理，具体包括以下：

设第j个测量数据处理前，目标i在模型r_i,k下的高斯伽玛混合分布及存在概率分别为和其中，i＝1,2,...,N_k-1,r_i,k＝1,2,…,M_r,j＝1,2,…,M；

根据变分贝叶斯方法，计算得到第j个测量更新时，目标i在模型r_i,k下的高斯伽玛混合分布为存在概率分别为其中，

<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>q</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mo>{</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>ln</mi> <mo>|</mo> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>t</mi> <mi>r</mi> <mrow> <mo>{</mo> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mo>-</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>t</mi> <mi>r</mi> <mrow> <mo>{</mo> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>m</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msubsup> <mi>m</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>m</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msubsup> <mi>m</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>T</mi> </msup> </mrow> <mo>}</mo> </mrow> <mo>+</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;xi;</mi> </munderover> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;beta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>ln</mi> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>+</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>ln</mi> <mo>|</mo> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <mo>-</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;xi;</mi> </munderover> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;beta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>ln</mi> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>+</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;eta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>ln</mi> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;eta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>ln</mi> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mi>n</mi> <mo>+</mo> <mn>2</mn> </mrow> <mn>2</mn> </mfrac> <mo>}</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> ；

表示伽玛函数，tr表示矩阵的迹，表示均值向量，表示协方差矩阵， 表示滤波器增益，伽玛分布的形状参数为和伽玛分布的尺度参数为 和 精度矩阵为其中， H_k为模型r_i,k的观测矩阵，P_D,k为目标的检测概率，λ_c,k为杂波密度，y_j,k为当前时刻接收到的第j个测量数据，上标T表示矩阵或向量的转置，ξ为状态向量的总维数；

若则第j个测量数据处理后，目标i在模型r_i,k下的的高斯伽玛混合分布为 目标i的存在概率为其中

若则第j个测量数据处理后，目标i在模型r_i,k下的高斯伽玛混合分布为 目标i的存在概率为其中，

在第M个测量数据处理后，目标i在模型r_i,k下的高斯伽玛混合分布和存在概率分别为和其中，i＝1,2,...,N_k-1，r_i,k＝1,2,…,M_r；

将第M个测量数据处理后，各个目标在不同模型下的高斯伽玛混合分布及存在概率分别作为当前时刻各个目标的更新高斯伽玛混合分布和存在概率，则所述当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布为 更新存在概率为其中i＝1,…,N_k-1，r_i,k＝1,2,…,M_r，

4.根据权利要求3所述的方法，其特征在于，所述步骤3包括：

对当前时刻目标i在M_r个运动模型下的更新的高斯伽玛混合分布和存在概率分别进行融合，得到当前时刻已存在目标q的运动模型为r_q,k、高斯伽玛混合分布为和存在概率为ρ_q,k(r_q,k)，其中，q为当前时刻已存在目标的索引号，且q＝1,…,N_k-1， 伽玛分布的形状参数为和伽玛分布的尺度参数为和

5.如权利要求4所述的方法，其特征在于，所述步骤4包括：

选取匀速运动模型作为新生目标的运动模型，设此时新生目标的运动模型为r_j,k＝1，其中，j＝1,…,M；

利用M个测量数据生成当前时刻新生目标的高斯伽玛混合分布为

其中，j＝1,…,M，为第j个新生目标的协方差，由第j个测量数据y_j，k＝[x_j，k y_j，k]^T产生，且 和为第j个新生目标分布函数中伽玛分布的形状参数，和为第j个新生目标分布函数中伽玛分布的尺度参数；

指定当前时刻各新生目标的存在概率为其中，j＝1,…,M，ρ_υ为所指定的存在概率；

将所述当前时刻已存在目标的运动模型与所述当前时刻新生目标的运动模型合并，得到当前时刻各个目标的运动模型为其中，N_k＝N_k-1+M；

将所述当前时刻已存在目标的高斯伽玛混合分布与所述当前时刻新生目标的高斯伽玛混合分布进行合并，得到当前时刻各个目标的高斯伽玛混合分布为

<mrow> <mtable> <mtr> <mtd> <mrow> <msubsup> <mrow> <mo>{</mo> <mrow> <mi>N</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>;</mo> <msub> <mi>m</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <munderover> <mi>&amp;Pi;</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;xi;</mi> </munderover> <mi>g</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>e</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mo>;</mo> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msubsup> <mi>&amp;beta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>g</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>v</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>;</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msub> <mi>&amp;eta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>k</mi> </msub> </msubsup> <mo>=</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mrow> <mo>{</mo> <mrow> <mi>N</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>x</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>;</mo> <msub> <mi>m</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msub> <mi>P</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <munderover> <mi>&amp;Pi;</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;xi;</mi> </munderover> <mi>g</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>e</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mo>;</mo> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msubsup> <mi>&amp;beta;</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>g</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>v</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>;</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msub> <mi>&amp;eta;</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mrow> <mi>q</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </msubsup> <mo>&amp;cup;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msubsup> <mrow> <mo>{</mo> <mrow> <mi>N</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>x</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>;</mo> <msubsup> <mi>m</mi> <mi>&amp;upsi;</mi> <mi>j</mi> </msubsup> <mrow> <mo>(</mo> <mrow> <msub> <mi>r</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>=</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <mo>,</mo> <msubsup> <mi>P</mi> <mi>&amp;upsi;</mi> <mi>j</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <munderover> <mi>&amp;Pi;</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;xi;</mi> </munderover> <mi>g</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>e</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mo>;</mo> <msubsup> <mi>&amp;alpha;</mi> <mi>&amp;upsi;</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msubsup> <mi>&amp;beta;</mi> <mi>&amp;upsi;</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>g</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>v</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>;</mo> <msubsup> <mi>&amp;gamma;</mi> <mi>&amp;upsi;</mi> <mi>j</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msubsup> <mi>&amp;eta;</mi> <mi>&amp;upsi;</mi> <mi>j</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </msubsup> </mtd> </mtr> </mtable> <mo>;</mo> </mrow>

将所述当前时刻各个目标的存在概率与所述当前时刻新生目标的存在概率进行合并，得到所述当前时刻各个目标的存在概率为

6.一种适用于闪烁噪声的多机动目标跟踪系统，其特征在于，所述系统包括：

预测模块，用于利用前一时刻各个目标的运动模型、以及高斯伽玛混合分布和存在概率产生伽玛分布的形状参数和尺度参数；根据前一时刻各个目标的运动模型、高斯伽玛混合分布和存在概率、当前时刻与前一时刻的时间差、模型间的马尔可夫转移概率、以及伽玛分布的形状参数和尺度参数，得到当前时刻各个目标在不同运动模型下的预测高斯伽玛混合分布和预测存在概率；

更新模块，用于根据当前时刻各个目标在不同运动模型下的预测高斯伽玛混合分布和预测存在概率，利用变分贝叶斯方法对当前时刻的测量数据进行序贯处理，得到当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布和更新存在概率；

融合模块，用于对当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布和更新存在概率分别进行融合，得到当前时刻已存在目标的运动模型、高斯伽玛混合分布和存在概率；

生成模块，用于利用当前时刻的测量数据生成新生目标的高斯伽玛混合分布，为所述新生目标指定存在概率，并选取匀速运动模型作为新生目标的运动模型，将所述新生目标的运动模型、高斯伽玛混合分布及存在概率分别与所述当前时刻已存在目标的运动模型、高斯伽玛混合分布及存在概率进行合并，得到当前时刻各个目标的运动模型、高斯伽玛混合分布和存在概率；

提取模块，用于从所述当前时刻各个目标中裁减掉存在概率小于第一阈值的目标，并将裁减后余下的目标的运动模型、高斯伽玛混合分布和存在概率作为滤波器下一次递归的输入；从所述裁减后余下的目标中提取存在概率大于第二阈值的目标，并将提取出的目标的高斯伽玛混合分布作为当前时刻的输出，将输出的高斯伽玛混合分布中，高斯分布的均值作为当前时刻目标的状态估计。

7.如权利要求6所述的系统，其特征在于，所述预测模块具体用于：

以k-1表示前一时刻，k表示当前时刻，t_k-1表示前一时刻的时间，t_k表示当前时刻的时间，r_i,k-1表示前一时刻目标i的运动模型标识，且1≤r_i,k-1≤M_r，r_i,k表示当前时刻目标i的模型标识，且r_i,k＝1,2,…,M_r，其中，M_r表示运动模型的总数目；当前时刻的测量噪声服从ξ维的学生t分布，以S(y_j,k；H_kx_i,k,R_i,k,ν_i,k)表示当前时刻测量的概率密度函数，其中，y_j,k表示当前时刻的一个测量，H_k表示测量矩阵，x_i,k表示目标的状态，R_i,k表示精度矩阵，ν_i,k表示t分布的自由度，且 为精度矩阵R_i,k主对角线上的各个元素；

已知k-1时刻目标i的运动模型为r_i,k-1，则在该模型r_i,k-1下的高斯伽玛混合分布为存在概率为ρ_i,k-1(r_i,k-1)；其中，N表示高斯分布，g表示伽玛分布，x_i,k-1表示k-1时刻目标i的状态向量，m_i,k-1(r_i,k-1)和P_i,k-1(r_i,k-1)分别表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中高斯分布的均值和协方差，表示k-1时刻精度矩阵的对角线元素，ν_i,k-1表示k-1时刻t分布的自由度，和γ_i,k-1(r_i,k-1)表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中伽玛分布的形状参数，和η_i,k-1(r_i,k-1)表示k-1时刻目标i在模型r_i,k-1下的高斯伽玛混合分布中伽玛分布的尺度参数，ξ表示状态向量的维数，i＝1,2,...,N_k-1，N_k-1表示前一时刻目标的总数；

根据k-1时刻目标i的运动模型、高斯伽玛混合分布和存在概率、当前时刻与前一时刻的时间差、以及模型间的马尔可夫转移概率，得到当前时刻目标i在模型r_i,k下的预测高斯伽玛混合分布为预测存在概率为ρ_i,k|k-1(r_i,k)；其中，m_i,k|k-1(r_i,k)＝F_k-1(r_i,k)m_i,k-1(r_i,k-1)， γ_i,k|k-1(r_i,k)＝ρ_γγ_i,k-1(r_i,k-1)，η_i,k|k-1(r_i,k)＝ρ_ηη_i,k-1(r_i,k-1)，ρ_i,k|k-1(r_i,k)＝p_s,k(t_k-t_k-1)t_k|k-1(r_i,k|r_i,k-1)ρ_i,k-1(r_i,k-1)，F_k-1(r_i,k)表示k-1时刻目标i在模型r_i,k下的状态转移矩阵，Q_k-1(r_i,k)表示k-1时刻目标i在模型r_i,k下的过程噪声协方差矩阵，上标T表示矩阵或向量的转置，ρ_α，ρ_β，ρ_γ，ρ_η表示传递因子，取值范围为(0,1]，p_s,k(t_k-t_k-1)表示目标的幸存概率，且t_k|k-1(r_i,k|r_i,k-1)表示模型间的马尔可夫转移概率，T表示采样周期，δ为预定的常数，l＝1,...,ξ，ξ表示状态向量的维数。

8.如权利要求7所述的系统，其特征在于，设当前时刻的测量集为y_k＝(y_1,k,...,y_M,k)，其中，M为当前时刻测量的总数，则所述更新模块具体用于：

以当前时刻各个目标在M_r个运动模型下的预测高斯伽玛混合分布和预测存在概率，作为当前时刻各个目标在相应运动模型下的初始高斯伽玛混合分布和初始存在概率，则目标i在模型r_i,k下的初始高斯伽玛混合分布取为 初始存在概率取为其中i＝1,2,...,N_k-1，r_i,k＝1,2,…,M_r，

利用变分贝叶斯方法对第1个至第M个当前时刻的测量数据进行序贯处理，具体包括以下：

设第j个测量数据处理前，目标i在模型r_i,k下的高斯伽玛混合分布及存在概率分别为和其中，i＝1,2,...,N_k-1,r_i,k＝1,2,…,M_r,j＝1,2,…,M；

根据变分贝叶斯方法，计算得到第j个测量更新时，目标i在模型r_i,k下的高斯伽玛混合分布为存在概率分别为其中，

<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>q</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mo>{</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>ln</mi> <mo>|</mo> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>t</mi> <mi>r</mi> <mrow> <mo>{</mo> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mo>-</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>t</mi> <mi>r</mi> <mrow> <mo>{</mo> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>m</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msubsup> <mi>m</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>m</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msubsup> <mi>m</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>T</mi> </msup> </mrow> <mo>}</mo> </mrow> <mo>+</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;xi;</mi> </munderover> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;beta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>ln</mi> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>+</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>ln</mi> <mo>|</mo> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <mo>-</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;xi;</mi> </munderover> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;beta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>ln</mi> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>+</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;eta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>ln</mi> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;eta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>ln</mi> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <msubsup> <mi>ln&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mrow> <mi>a</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mi>n</mi> <mo>+</mo> <mn>2</mn> </mrow> <mn>2</mn> </mfrac> <mo>}</mo> </mrow> </mtd> </mtr> </mtable> </mfenced>

表示伽玛函数，tr表示矩阵的迹，表示均值向量，表示协方差矩阵， 表示滤波器增益，伽玛分布的形状参数为和伽玛分布的尺度参数为和精度矩阵为 其中， H_k为模型r_i,k的观测矩阵，P_D,k为目标的检测概率，λ_c,k为杂波密度，y_j,k为当前时刻接收到的第j个测量数据，上标T表示矩阵或向量的转置，ξ为状态向量的总维数；

若则第j个测量数据处理后，目标i在模型r_i,k下的的高斯伽玛混合分布为 目标i的存在概率为其中

若则第j个测量数据处理后，目标i在模型r_i,k下的高斯伽玛混合分布为 目标i的存在概率为其中，

在第M个测量数据处理后，目标i在模型r_i,k下的高斯伽玛混合分布和存在概率分别为和其中，i＝1,2,...,N_k-1，r_i,k＝1,2,…,M_r；

将第M个测量数据处理后，各个目标在不同模型下的高斯伽玛混合分布及存在概率分别作为当前时刻各个目标的更新高斯伽玛混合分布和存在概率，则所述当前时刻各个目标在不同运动模型下的更新高斯伽玛混合分布为

更新存在概率为其中i＝1,…,N_k-1，r_i,k＝1,2,…,M_r，

9.如权利要求8所述的系统，其特征在于，所述融合模块具体用于：

对当前时刻目标i在M_r个运动模型下的更新的高斯伽玛混合分布和存在概率分别进行融合，得到当前时刻已存在目标q的运动模型为r_q,k、高斯伽玛混合分布为和存在概率为ρ_q,k(r_q,k)，其中，q为当前时刻已存在目标的索引号，且q＝1,…,N_k-1， 伽玛分布的形状参数为和伽玛分布的尺度参数为和

10.如权利要求9所述的系统，其特征在于，所述生成模块具体用于：

选取匀速运动模型作为新生目标的运动模型，设此时新生目标的运动模型为r_j,k＝1，其中，j＝1,…,M；

利用M个测量数据生成当前时刻新生目标的高斯伽玛混合分布为

其中，j＝1,…,M，为第j个新生目标的协方差，由第j个测量数据y_j，k＝[x_j，k y_j，k]^T产生，且 和为第j个新生目标分布函数中伽玛分布的形状参数，和为第j个新生目标分布函数中伽玛分布的尺度参数；

指定当前时刻各新生目标的存在概率为其中，j＝1,…,M，ρ_υ为所指定的存在概率；

将所述当前时刻已存在目标的运动模型与所述当前时刻新生目标的运动模型合并，得到当前时刻各个目标的运动模型为其中，N_k＝N_k-1+M；

将所述当前时刻已存在目标的高斯伽玛混合分布与所述当前时刻新生目标的高斯伽玛混合分布进行合并，得到当前时刻各个目标的高斯伽玛混合分布为

<mrow> <mtable> <mtr> <mtd> <mrow> <msubsup> <mrow> <mo>{</mo> <mrow> <mi>N</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>;</mo> <msub> <mi>m</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <munderover> <mi>&amp;Pi;</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;xi;</mi> </munderover> <mi>g</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>e</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mo>;</mo> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msubsup> <mi>&amp;beta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>g</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>v</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>;</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msub> <mi>&amp;eta;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>k</mi> </msub> </msubsup> <mo>=</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mrow> <mo>{</mo> <mrow> <mi>N</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>x</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>;</mo> <msub> <mi>m</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msub> <mi>P</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <munderover> <mi>&amp;Pi;</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;xi;</mi> </munderover> <mi>g</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>e</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mo>;</mo> <msubsup> <mi>&amp;alpha;</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msubsup> <mi>&amp;beta;</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>g</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>v</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>;</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msub> <mi>&amp;eta;</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mrow> <mi>q</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </msubsup> <mo>&amp;cup;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msubsup> <mrow> <mo>{</mo> <mrow> <mi>N</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>x</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>;</mo> <msubsup> <mi>m</mi> <mi>&amp;upsi;</mi> <mi>j</mi> </msubsup> <mrow> <mo>(</mo> <mrow> <msub> <mi>r</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>=</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <mo>,</mo> <msubsup> <mi>P</mi> <mi>&amp;upsi;</mi> <mi>j</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <munderover> <mi>&amp;Pi;</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&amp;xi;</mi> </munderover> <mi>g</mi> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>e</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> <mi>l</mi> </msubsup> <mo>;</mo> <msubsup> <mi>&amp;alpha;</mi> <mi>&amp;upsi;</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msubsup> <mi>&amp;beta;</mi> <mi>&amp;upsi;</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>l</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>g</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>v</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>;</mo> <msubsup> <mi>&amp;gamma;</mi> <mi>&amp;upsi;</mi> <mi>j</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <msubsup> <mi>&amp;eta;</mi> <mi>&amp;upsi;</mi> <mi>j</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </msubsup> </mtd> </mtr> </mtable> <mo>;</mo> </mrow>

将所述当前时刻各个目标的存在概率与所述当前时刻新生目标的存在概率进行合并，得到所述当前时刻各个目标的存在概率为