基于稀疏表示的多传感器航迹融合方法
【专利摘要】本发明提供了一种基于稀疏表示的多传感器航迹融合方法。该多传感器航迹融合方法包括：根据各个局部传感器航迹的当前时刻状态信息构建无噪声训练样本集T和含噪声训练样本集R；基于两者建立字典训练模型；对字典训练模型进行训练，学习出了一对无噪声、含噪声字典D、C，使得无噪声、含噪声训练样本T、R分别在D、C中能够被稀疏的表示，并且具有相同的表示系数G；每一时刻把所有局部传感器得到的航迹状态向量表示成一个向量，找出该向量在含噪声字典中的稀疏表示系数，通过无噪声字典D和稀疏表示系数得到融合后的全局航迹状态向量。本发明能够得到了比线性最小均方误差(LMMSE)融合方法更好的融合效果；同时，能够融合各种各样的信息。
【专利说明】
基于稀疏表示的多传感器航迹融合方法
技术领域
[0001 ]本发明设及微波行业目标追踪技术领域，尤其设及一种基于稀疏表示的多传感器 航迹融合方法。
【背景技术】
[0002] 依据多传感器的信息流通形式及对信息的综合处理层次，多传感器数据融合目标 跟踪系统中多源信息的融合处理结构主要有Ξ种类型:集中式、分布式和混合式。分布式融 合系统结构具有较强的生存能力和稳定性，对于数据传输的信道带宽和融合中屯、的计算能 力要求都相对较低，所W分布式结构在实际应用中成为一种重要的选择。
[0003] 分布式结构中，各传感器首先在自身处理器下完成单传感器目标跟踪，然后将局 部传感器航迹送往融合中屯、，融合中屯、使用航迹融合算法把多个局部传感器航迹融合为全 局航迹，而航迹融合则是指把表示同一目标的多个局部航迹融合为一个全局航迹。目前，人 们已经提出了一些多传感器航迹融合方法，主要有简单凸组合航迹融合算法、互协方差组 合航迹融合算法、协方差交叉航迹融合算法、伪量测融合算法、线性最小均方差估计 (LMMSE)融合方法等。LMMSE融合方法具有简单、线性、无偏等优点，因而被广泛用于分布式 多传感器目标跟踪中。
[0004] 近年来，稀疏表示问题已成为人们的研究热点。信号的稀疏表示就是将信号投影 到正交变换基时，绝大部分变换系数的绝对值很小，所得到的变换向量是稀疏或者近似稀 疏的。传统的信号表示理论基于正交线性变换，但许多信号是各种自然现象的混合体，运些 混合信号在单一的正交基变换中不能非常有效地表现出来。最近几年，研究人员在改变传 统信号表示方面取得了很大的进展。新的信号表示理论的基本思想是:基函数用称之为字 典的超完备的冗余函数系取代，字典的选择尽可能好的符合被逼近信号的结构，其构成可 W没有任何限制，字典中的元素被称为原子。从字典中找到m个原子的最佳线性组合来表示 一个信号，称作信号的稀疏逼近或高度非线性逼近。
[0005] 然而，在实现本发明的过程中，
【申请人】发现现有的LMMSE融合方法的融合效果并不 理想。业内亟需一种具有更好融合效果的多传感器航迹融合方法。
【发明内容】
[0006] (一)要解决的技术问题
[0007] 鉴于上述技术问题，本发明提供了一种基于稀疏表示的多传感器航迹融合方法， W提高航迹融合效果。
[000引（二)技术方案
[0009] 本发明基于稀疏表示的多传感器航迹融合方法包括：
[0010] 步骤A:根据各个局部传感器航迹的当前时刻状态信息构建无噪声训练样本集T;
[0011] 步骤B:根据传感器产生的局部传感器航迹构建含噪声训练样本集R;
[0012] 步骤C:基于无噪声训练样本集T和含噪声训练样本集R，建立字典训练模型；
[0013] 步骤D:对字典训练模型进行训练，学习出了一对无噪声、含噪声字典D、C，使得无 噪声、含噪声训练样本T、R分别在无噪声、含噪声字典D、C中能够被稀疏的表示，并且具有相 同的表示系数G;W及
[0014] 步骤E:每一时刻把所有局部传感器得到的航迹状态向量表示成一个向量，找出该 向量在含噪声字典中的稀疏表示系数，通过无噪声字典D和稀疏表示系数得到融合后的全 局航迹状态向量。
[001引巧益效果
[0016] 从上述技术方案可W看出，本发明基于稀疏表示的多传感器航迹融合方法具有W 下有益效果：
[0017] (1)基于多传感器航迹融合和信号的稀疏表示，提出了一种基于稀疏表示的多传 感器航迹融合方法，能够得到了比线性最小均方误差(LMMSE)融合方法更好的融合效果；
[0018] (2)把信号的稀疏表示应用在多传感器航迹融合中，给出了一个多源信息融合的 框架，能够融合各种各样的信息。
【附图说明】
[0019] 图1为根据本发明实施例基于稀疏表示的多传感器航迹融合方法的流程图；
[0020] 图2给出了仿真中的目标真实轨迹和传感器位置分布；
[0021 ]图3给出了仿真中目标1滤波轨迹X向均方根误差；
[0022] 图4给出了仿真中目标1滤波轨迹y向均方根误差；
[0023] 图5给出了仿真中目标2滤波轨迹X向均方根误差；
[0024] 图6给出了仿真中目标2滤波轨迹y向均方根误差。
【具体实施方式】
[0025] 本发明首先基于给定的样本航迹构建训练样本集;然后通过字典学习给出一对能 够紧密表示训练样本的含噪声字典和无噪声字典;最后在每一跟踪时刻把所有传感器得到 的局部航迹状态向量按照一定的格式表示成一个向量，在含噪声字典中找出该向量的稀疏 表示系数，通过无噪声字典和稀疏表示系数得到融合后的全局航迹状态向量。
[0026] 为使本发明的目的、技术方案和优点更加清楚明白，W下结合具体实施例，并参照 附图，对本发明进一步详细说明。需要说明的是，在附图或说明书描述中，相似或相同的部 分都使用相同的图号。并且，附图中未绘示或描述的实现方式，为所属技术领域中普通技术 人员所知的形式。
[0027] 在本发明的一个示例性实施例中，提供了一种基于稀疏表示的多传感器航迹融合 方法。图1为根据本发明实施例基于稀疏表示的多传感器航迹融合方法的流程图。如图1所 示，本实施例基于稀疏表示的多传感器航迹融合方法包括：
[0028] 步骤A:根据各个局部传感器航迹的当前时刻状态信息构建无噪声训练样本集T; [00巧]该步骤A具体包括：
[0030]子步骤A1:在每隔Δ ( Δ为整数，Δ > 1)个跟踪时刻，根据各个局部传感器航迹的 当前时刻状态信息，限定一个位置范围和速度范围，在运个范围内随机生成η条样本航迹 (ti，t2，···，tn}，其中，η为样本航迹数目，ti表示第i条样本航迹，（i = l，2,…，η);
[0031] 该η条样本航迹中包括匀速直线运动目标航迹、机动目标航迹等。并且，样本航迹 数目η需要基于实际的计算性能和需求提前设定，跟踪时刻的个数、位置范围和速度范围由 用户根据自身需要设定。
[0032] 其中，对于样本航迹ti，对应着一个由各个时刻状态向量组成的矩阵：
[003；3] Ti=[xi(l),xi(2),...,xi(ki)]，i = l,2,...,n (1)
[0034] 其中，k功样本航迹ti的跟踪步数，Xi化)表示样本航迹ti在第k个跟踪时刻的状态 向量。
[0035] 子步骤A2:n条样本航迹中每条样本航迹的由各个时刻状态向量组成的矩阵构成 无噪声训练样本集Τ= [Τι，T2，…，Τη]。
[0036] 步骤Β:根据传感器产生的局部传感器航迹构建含噪声训练样本集R;
[0037] 该步骤Β具体包括：
[0038] 子步骤Β1:假定传感器数目为Ns，每个传感器产生的局部传感器航迹称为含噪声 航迹，对于样本航迹ti (i = 1，2，…，η)，传感器Si (1 = 1，2，…，Ns)有一个含噪声航迹与其对 应，该含噪声航迹对应着一个由各个时刻含噪声状态向量组成的矩阵化i:
[0039] Uii=[uii(l),山i(2),山i(ki)]，1 = 1,2,...,化，i = l ,2,... ,η (2)
[0040] 子步骤Β2:为了更好地学习含噪声航迹的特征，在化1中提取航迹的方向信息，与 uii化)对应，定义rii化)如下：
[0041]
[0044] 与化i相比，Rii不仅含有含噪声航迹每个时刻的状态向量，而且还包含每个时刻的 万向f曰息。
[0045] 子步骤B3:构建传感器Si(l = l，2,…，Ns)对应的含噪声训练样本集Ri=[Rn， Ri2 , ,Rln];
[0046] 子步骤B4:由所有传感器含噪声训练样本集构建含噪声训练样本集R:
[0047]
(5；
[004引此时，我们得到了最终用于学习无噪声字典和含噪声字典的无噪声训练样本集T 和含噪声训练样本集R。
[0049] 步骤C:基于无噪声训练样本集T和含噪声训练样本集R，建立字典训练模型；
[0050] 该步骤C具体包括：
[0051] 子步骤C1:基于无噪声训练样本集T和含噪声训练样本集R，无噪声字典D和含噪声 字典c训练模型如下：
[00 对
（6)
[0053] 其中，G表示系数矩阵，在求解该训练模型可同时得出，C为矩阵R的行数，d为矩阵T 的行数，λ为正则化参数(可取值为0.2)。表达式中的第一项和第二项是为了使得含噪声、无 噪声训练样本分别在含噪声、无噪声字典中具有相同的表示系数，第Ξ项是为了使得表示 系数更具稀疏性。
[0054] 子步骤C2:构建上述训练模型的等价形式：
[0057]相对于第一个模型要训练两个字典而言，该等价的模型可W通过单一字典的训练 达到同样效果。
[005引步骤D:对字典训练模型进行训练，学习出了一对无噪声、含噪声字典D、C，使得无 噪声、含噪声训练样本T、R分别在D、C中能够被稀疏的表示，并且具有相同的表示系数G;
[0059] 在本实施例中，基于上述等价的单一字典训练模型，利用K-SVD算法，依据无噪声、 含噪声训练样本T、R分别在D、C中能够被最稀疏表示的准则，迭代求解出一对最优的无噪 声、含噪声字典D、C，使得训练样本T、R分别在D、C中能够被稀疏表示，并且具有相同的表示 系数G。
[0060] 该步骤D具体包括：
[0061 ] Step 1:用随机生成方法初始化字典D和C;
[0062] Step 2:基于单一字典训练模型，利用K-SVD算法迭代求解出一对能够分别稀疏表 示T、R的字典D、C。
[0063] 步骤E:每一时刻把所有局部传感器得到的航迹状态向量表示成一个向量，找出该 向量在含噪声字典中的稀疏表示系数，通过无噪声字典D和稀疏表示系数得到融合后的全 局航迹状态向量。
[0064] 假设在某一跟踪时刻k(k> 3)，Ns个传感器给出的同一目标在该时刻的航迹状态 向量分别为_1.'1(。..，':(。、~，.，'、.斯，并已知该目标在时亥化-2^-1的状态向量分别为7化-2)、 y化-1)，令：
[00 化]
[0066]设
[0067]
(9)
[006引求解如下模型，得到向量z在含噪声字典C中的表示系数向量g:
[0069]
CIO)
[0070] 融合后的状态向量y化)为：
[0071] y(k)=Dg,k> 3 (11)
[0072] 其中，由航迹方向信息的定义，当k=l，2时，因为不能提取航迹的方向信息，所W 要单独计算。
[0073] 当跟踪时刻k = l，2时，令：
[0074]
(12)
[0075] 在含噪声字典C中选取与各个传感器航迹状态向量对应的行构成新的含噪声字典 C，求解如下模型：
[007引
（13)
[0077]得到向量2在新含噪声字典C中的表示系数向量f，侧融合后的状态向量y化)为： [007引.1作）=巧含，& = 1，2 (14)
[0079] 下面通过仿真的方法对本实施例进行验证。在仿真中，我们使用两个传感器对两 个目标进行跟踪。W融合中屯、为坐标原点，各传感器的位置分别为(-30000m，55000m)、（- 80000m,35000m)。两个目标均作匀速直线运动，目标的状态向量由目标的X和y轴位置、X和y 轴速度构成，两个目标的初始状态向量分别为
[0080] (-40000m,40000m,-400m/s,90m/s)
[0081 ] (-40000m，45000m，-450m/s，-20m/s)
[0082] 各个传感器的性能参数相同，参数设置如下:扫描周期为2s，探测概率为l，x轴和y 轴测量误差标准差均为100m。仿真的跟踪步数为40,过程噪声标准差为10m。图2给出了两个 目标的真实轨迹和传感器分布，"Δ"表示传感器所在位置。
[0083] 在每隔20个跟踪时刻，根据各个局部传感器航迹的位置信息和速度信息，在一定 的位置和速度范围内随机生成200条无噪声样本航迹，其中包含匀速直线运动目标航迹和 机动目标航迹。给无噪声航迹加入适当的过程噪声得到含噪声样本航迹，过程噪声标准差 为10m。通过字典学习得到含噪声字典和无噪声字典，字典原子个数为2048。
[0084] 针对上述给定的仿真场景和传感器参数设置，使用分布式多传感器多假设跟踪方 法进行数据关联和航迹关联，分别使用LMMSE融合方法和本实施例进行多传感器航迹融合。 图3和图4分别给出了目标1滤波轨迹在各个时刻的X向和y向均方根误差。图5和图6分别给 出了目标2滤波轨迹在各个时刻的X向和y向均方根误差。蒙特卡洛仿真次数为50。
[0085] 从图3~图6可W看出，当跟踪稳定后，与LMMSE融合方法相比，本实施例得到的滤 波轨迹在各个时刻有更低的均方根误差，也就是说，本实施例融合效果优于LMMSE融合方 法。
[0086] 至此，已经结合附图对本实施例进行了详细描述。依据W上描述，本领域技术人员 应当对本发明基于稀疏表示的多传感器航迹融合方法有了清楚的认识。
[0087] 此外，上述对各元件和方法的定义并不仅限于实施例中提到的各种具体结构、形 状或方式，本领域普通技术人员可对其进行简单地更改或替换，例如：
[008引（l)rii化)定义还可W改成如下形式：
[0089]
(14)
[0090] (2化-SVD方法可W用PCA方法来代替；
[0091] 另外，虽然本文可提供包含特定值的参数的示范，但应了解，参数无需确切等于相 应的值，而是可在可接受的误差容限或设计约束内近似于相应的值。实施例中提到的方向 用语，例如"上"、"下"、"前后"、"左"、"右"等，仅是参考附图的方向。因此，使用的方向用 语是用来说明并非用来限制本发明的保护范围。
[0092] 综上所述，本发明基于多传感器航迹融合和信号的稀疏表示，提出了一种基于稀 疏表示的多传感器航迹融合方法，能够得到了比线性最小均方误差(LMMSE)融合方法更好 的融合效果。同时，把信号的稀疏表示应用在多传感器航迹融合中，给出了一个多源信息融 合的框架，能够融合各种各样的信息。
[0093] W上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详 细说明，所应理解的是，W上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡 在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保 护范围之内。
【主权项】
1. 一种基于稀疏表示的多传感器航迹融合方法，其特征在于，包括： 步骤A:根据各个局部传感器航迹的当前时刻状态信息构建无噪声训练样本集T; 步骤B:根据传感器产生的局部传感器航迹构建含噪声训练样本集R; 步骤C:基于无噪声训练样本集T和含噪声训练样本集R，建立字典训练模型； 步骤D:对字典训练模型进行训练，学习出了一对无噪声、含噪声字典D、C，使得无噪声、 含噪声训练样本T、R分别在无噪声、含噪声字典D、C中能够被稀疏的表示，并且具有相同的 表示系数G;以及 步骤E:每一时刻把所有局部传感器得到的航迹状态向量表示成一个向量，找出该向量 在含噪声字典中的稀疏表示系数，通过无噪声字典D和稀疏表示系数得到融合后的全局航 迹状态向量。2. 根据权利要求1所述的多传感器航迹融合方法，其特征在于，所述步骤A具体包括： 子步骤A1:在每隔△个跟踪时刻，根据各个局部传感器航迹的当前时刻状态信息，限定 一个位置范围和速度范围，在这个范围内随机生成η条样本航迹{?ι，?2，···，?η}，其中，η为样 本航迹数目，ti表示第i条样本航迹，（i = l，2,…，η); 其中，对于样本航迹^，对应着一个由各个时刻状态向量组成的矩阵： Ti = [xi(l)，xi(2)，···，xi(ki)]，i = l，2,…，η 其中，ki为样本航迹ti的跟踪步数，xi (k)表示样本航迹ti在第k个跟踪时刻的状态向量； 子步骤A2:由η条样本航迹中每条样本航迹的由各个时刻状态向量组成的矩阵构成无 噪声训练样本集Τ = [ Ti，Τ2，…，Τη]。3. 根据权利要求2所述的多传感器航迹融合方法，其特征在于，所述步骤B具体包括： 子步骤B1:假定传感器数目为Ns，每个传感器产生的局部传感器航迹称为含噪声航迹， 对于样本航迹ti (i = l，2,…，η)，传感器Si (1 = 1，2，…，Ns)有一个含噪声航迹与其对应，该 含噪声航迹对应着一个由各个时刻含噪声状态向量组成的矩阵Uh: Uii = [uii(l)，uii(2)，…，uii(ki)]，1 = 1，2,…，NS，i = 1，2,…，n (2) 子步骤B2:为了更好地学习含噪声航迹的特征，在Uu中提取航迹的方向信息，与im(k) 对应，定义rii(k)如下：子步骤B3:构建传感器Si(l = l，2,…，Ns)对应的含噪声训练样本集Ri=[Rn，Ri2，···， Rln]； 子步骤B4:由所有传感器含噪声训练样本集构建含噪声训练样本集R:4. 根据权利要求3所述的多传感器航迹融合方法，其特征在于，所述步骤C具体包括： 子步骤C1:基于无噪声训练样本集T和含噪声训练样本集R，无噪声字典D和含噪声字典 C训练模型如下：其中，G表示系数矩阵，在求解该训练模型可同时得出，c为矩阵R的行数，d为矩阵T的行 数，λ为正则化参数； 子步骤C2:构建上述训练模型的等价形式：5. 根据权利要求4所述的多传感器航迹融合方法，其特征在于，所述步骤D具体包括:基 于所述字典训练模型，利用K-SVD算法，依据无噪声、含噪声训练样本T、R分别在D、C中能够 被最稀疏表示的准则，迭代求解出一对最优的无噪声、含噪声字典D、C，使得训练样本T、R分 别在D、C中能够被稀疏表示，并且具有相同的表示系数G。6. 根据权利要求5所述的多传感器航迹融合方法，其特征在于，所述步骤E具体包括： 假设在某一跟踪时刻k(k2 3)，NS个传感器给出的同一目标在该时刻的航迹状态向量分 别为艽(办艿(幻(幻，并已知该目标在时刻k-2、k_l的状态向量分别为y(k_2)、y(k_ 1)，令：求解如下模型，得到向量z在含噪声字典C中的表示系数向量g: 融合后的状态向量y(k)为：y(k)=Dg，k> 3。7. 根据权利要求6所述的多传感器航迹融合方法，其特征在于，所述步骤E还包括： 当跟踪时刻k=l，2时，令：在含噪声字典C中选取与各个传感器航迹状态向量对应的行构成新的含噪声字典C,求 解如下模型：得到向量f在新含噪声字典€中的表示系数向量#，则融合后的状态向量y(k)为： ,y(fc) = 2：f，k=l，2〇
【文档编号】G01S7/02GK105825239SQ201610196780
【公开日】2016年8月3日
【申请日】2016年3月31日
【发明人】张冰尘, 王欢, 孙进平, 洪文, 蒋成龙, 毕辉
【申请人】中国科学院电子学研究所