[0077] 以下结合附图对本发明作进一步说明。
[0078] 本发明提出了一种针对有限集跟踪滤波器的均匀密集杂波稀疏方法,其具体实施方式如下:
[0079] 步骤1.系统建模;
[0080] 假设目标和杂波产生的的量测服从下面的二项式混合分布:
[0081]
[0082] 其中,f(·|·)是条件分布函数,zk,i是目标或杂波的量测, 是目标状态分布的权重, 其中上标t表示目标,c表示杂波,假定在k时刻有 个杂波状态和 个目标状态, 是杂波在k时刻的状态集, 是目标在k时刻
的状态集,杂波和目标的量测的混合分布 分别具有如下形式:
[0083]
[0084]
[0085] 其中, 和 分别是杂波和目标的第 和第 个元素的组成权重。这里假设杂波的量测 服从泊松随机有限集分布,且泊松强度为 同样的,目标的量测 服从伯努利随机有限集分布。
[0086] 当杂波的量测服从均匀分布时,公式(2)中的杂波混合分布变成则杂波强度混合分布具有如下形式:
[0087]
[0088] 其中Sc是杂波面积或体积。
[0089] 步骤2杂波稀疏过程
[0090] 步骤2.1杂波分布
[0091] 杂波密度在公式(2)和(4)中被描述为混合分布,因此,杂波的随机有限集Θk具有如下分布:Θk~U(zk,i|Sc) (5)
[0092] 在公式(2)和(3)中分别给出了杂波和目标的量测的混合分布,即这里, 和 分别是杂波和目标的量测集中的某个元
素。如步骤1中所述,杂波的量测 服从泊松随机有限集分布,且泊松强度为 同样的,目标的量测 服从多伯努利随机有限集分布。
[0093] 步骤2.2假设检验原理在随机有限集量测中的应用
[0094] 假设z是某一个时刻的量测,为了简单起见,省略了时间下标,二元假设检验问题可以通过杂波和目标的量测表示成下面的形式:
[0095]
[0096] 因此,似然函数的条件假设为
[0097]
[0098]
[0099] 似然比检验为
[0100]
[0101] 这里的阈值η可以由下面的虚警概率导出:
[0102]
[0103] 其中 是积分区间。
[0104] 对于服从简单高斯分布的单目标来说,它就变成了普通的二元假设,其阈值η可以通过下面的公式获得:
[0105]
[0106] 积分区间可以为[zη,∞)或者(-∞,zη],为了不失一般性,采用[zη,∞)作为积分区间,且通过r(z)>η将其导出;
[0107] 步骤2.3基本命题
[0108] 首先通过一个引理给出向量切比雪夫不等式
[0109] 引理1:假设随机变量z的维数是n,E(z)是均值,Σ是协方差并且γ>0,那么下面的不等式成立:
[0110]
[0111] 令 并且假设目标的量测z服从下面的高斯混合分布:
[0112]
[0113] 目标的量测集定义为:
[0114]
[0115] 这里的p(·)代表混合目标的第 个组成 的概率分布,当目标量测概率小于公式(9)中的假报警概率PF时,Γ将是一个充分大的阈值。
[0116] 杂波的量测集定义为:
[0117]
[0118] 这里的p(·)代表混合目标的第 个组成 的概率分布,ζ>0,β是杂波量测概率的阈值。
[0119] 当 时,下面的命题成立:
[0120] 命题1:如果满足公式(13)和(14),那么下面的不等式成立
[0121] p(z∈Tk|H1)≥Γ (15)
[0122] p(z∈Ck|H0)≥β (16)
[0123] 这里的p(z∈Tk|H1)是(12)式中给出的混合分布,定义上面的公式证明如下:
[0124]
[0125] 同理可证不等式(16)。
[0126] 命题2:假设对于任意的z∈Sk,有Sk=Tk∪Ck,那么下面的等式成立
[0127]
[0128]
[0129] 这个命题给出了似然比的范围。很明显r(z∈Sk)≥0,因此,仅考虑r(z∈Sk)≤rk,max的情况。
[0130] 证明如下:
[0131]
[0132] 因此
[0133] 即证。
[0134] 当混合分布p(z|H1)为高斯分布时,有下面的推论
[0135] 推论1:当目标量测满足(12)式中的高斯混合分布时,似然比可以简化成如下形式:
[0136]
[0137] 命题3:假设目标量测满足(12)式中的高斯混合分布,并且检测概率的阈值是pD,min,那么似然比的最小阈值可通过下面的式子求得
[0138]
[0139]
[0140] 命题4:定义量测集为Rk={z:rk(z)≥rk,min},则当Fk=Rk∩Ck时,假报警概率具有如下形式:
[0141] 证明如下:
[0142]
[0143] 命题5:假设Sk=Tk∪Ck,并且杂波服从均匀分布,则则稀疏过程的杂波强度为:λk,s(Sk)=λk,cVol(Tk),这里的λk,s(Sk)是区域Sk上的杂波强度,Vol(Tk)代表目标量测集的体积。
[0144] 图1给出了杂波和目标的概率密度函数和似然比函数,假设目标量测的混合分布为:f(y)=0.5N(y;-30,4)+0.5N(y;40,2),杂波的量测的分布为c(y)=U([-100,100]),图1中上面的图是目标量测和杂波的量测的概率密度函数,下面的图是似然比函数,从图中可以看出,似然比函数与目标的概率密度函数具有相同的波形。
[0145] 图2给出了检测概率PD和假报警概率PF与似然比的关系图,从图中可以得到,如果假报警概率等于10%,相应的检测概率将会大于90%。
[0146] 为了更好地阐释说明本发明,在本发明实验中,我们选择高斯混合势概率假设密度和高斯混合多伯努利滤波器并通过假设检验来验证稀疏过程。目标的检测概率PD(xk)=0.98。监控区域S=[-1000,1000]×[-1000,1000]m2。根据转移密度pk|k-1(xk|xk-1)=N(xk,Fk-1,Qk-1),知单目标状态空间模型是线性高斯分布的,似然比函数为gk(zk|xk)=N(zk;Hkxk,
Rk),其中 Rk
=diag([100,100])m2,这里的 是Kronecker积,Δ=1s是采样周期。过程噪声协方差为Wk=diag([25,25])m2。新生目标出现在四个不同的位置 其中
wγ(1)=0.1,wγ(2)=wγ(3)=wγ(4)=0.05μγ(1)=[0,0,0,0]T,μγ(2)=[400,0,-600,0]T,μγ(3)=[-800,0,-200,0]T,μγ(4)=[-200,0,800,0]T,pγ=diag([10,10,10,10])。在仿真实验中,杂波密度的范围为1.25×10-5m-2到7.5×10-5m-2,相应的杂波个数的范围是50到300。
[0147] 图3给出了通过100次蒙特卡洛试验在GM-CPHD滤波器下稀疏算法与传统算法的关于最优次模式分配(OSPA)距离与杂波个数的关系的对比图。
[0148] 图4给出了通过100次蒙特卡洛试验在GM-CPHD滤波器下稀疏算法与传统算法的关于势距离和位置距离与杂波个数的关系的对比图。
[0149] 图5给出了通过100次蒙特卡洛试验在GM-CPHD滤波器下稀疏算法与传统算法的关于计算时间与杂波个数的关系的对比图。从图中可以看到,传统算法的计算时间随着杂波密度的增加呈指数增长,而稀疏算法的计算时间却呈线性增长。比如当杂波个数等于300时,稀疏算法100次蒙特卡洛试验需要30s,而传统算法则需要500s。
[0150] 图6给出了通过100次蒙特卡洛试验在GMmulti-Bernoulli滤波器下稀疏算法与传统算法的关于最优次模式分配(OSPA)距离与杂波密度的关系的对比图。
[0151] 图7给出了通过100次蒙特卡洛试验在GMmulti-Bernoulli滤波器下稀疏算法与传统算法的关于势距离和位置距离与杂波密度的关系的对比图。
[0152] 图8给出了通过100次蒙特卡洛试验在GMmulti-Bernoulli滤波器下稀疏算法与传统算法的关于计算时间与杂波密度的关系的对比图。
[0153] 通过图3~图8的对比分析,我们可以得出以下结论:
[0154] (1)GM-CPHD滤波器在目标数目估计方面比GMmulti-Bernoulli滤波器具有更好的鲁棒性。
[0155] (2)不管是在GM-CPHD滤波器下还是在GMmulti-Bernoulli滤波器下,稀疏算法都可以大大减少计算时间,从而体现了本发明方案的优越性。
[0156] 最后说明,以上描述仅用以说明本发明的技术方案而非限制其所包含范围,即对本发明的技术方案进行修改或者等同替换,而并未脱离其目的和范围的,均应涵盖于本发明的权利要求范围当中。