资源描述:
《基于神经网络补偿的多传感器航迹融合(2009研究生数学建模c题》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、基于神经网络补偿的多传感器航迹融合摘 要:针对多传感器环境的条件提出了一种基于神经网络补偿的航迹融合方法。各传感器的测量值用线性卡尔曼滤波器进行处理并将获得的局部航迹传送到融合中心。首先对局部航迹进行融合,然后引入神经网络来减少因共同过程噪声而导致的融合估计误差,其中神经网络采用DanSi2mon提出的网络结构,并对神经网络权值的优化采用无痕卡尔曼滤波(UKF)。仿真结果表明,这种融合方法对跟踪具有过程噪声的目标非常有效,而且过程噪声发生变化时该方法仍是有效的,从而使得它在很多实际应用中具有潜在的价值。关键词:航迹融合;多传感器;径向基函数神经网络;无痕卡尔曼滤波 近年来随着多传感器设
2、备的广泛使用,多传感器融合或数据融合已经受到越来越多的关注。目前,多传感器数据融合主要有两种方法:测量融合和航迹融合。测量融合是将传感器的测量值进行组合并获得目标状态向量的最优估计,这种方法的主要缺陷在于它需要传感器的测量值并且集中式方法缺乏鲁棒性。另外,传感器的测量值可能属于不同类型(比如红外,雷达等),而要同时处理这些不同类型的测量值则是非常困难的。因此,在很多的实际应用场合均使用航迹融合方法,每一个传感器都应用一个估计器来获取目标的状态向量及其协方差矩阵,并将其通过数据链传送到融合中心,在融合中心进行航迹关联和航迹融合从而获得更精确的目标状态向量。在航迹融合方法里,Singer等首
3、先提出了简单融合算法(SF)。它假定来自不同传感器的同一目标的航迹估计误差是不相关的,而实际上,这种假设:不成立,因为对各个估计器的滤波动态来说,由目标机动导致的过程噪声是共同的。为了消除不同航迹间的关联,Bar-Shalom等提出了加权方差融合方法(WCF),结果表明在协方差为正定的情况下,这种融合方法相对于简单融合方法的性能有很大幅度的提高,但是在非正定的情况下,加权方差融合方法的性能却下降了,甚至比简单融合方法还要差,另外,对交互协方差的计算也非常复杂。为了减少计算负担,Celine提出了一种自适应航迹融合方法。这种方法的要点在于它有一个决策逻辑,即通过距离矩阵和门限的比较来选择简
4、单融合或加权方差融合,以获得全局估计。但由于过程噪声的存在,大多数情况下仍然需要使用加权方差融合方法,无法避免计算交互协方差。稳定状态下交互式协方差的解析解可通过离散李亚普诺夫方程来求解,Saha通过引入双线性变换转化交互式协方差矩阵方程并通过积分求解而大大降低了计算负担。但是这些方法对于二维状态向量的融合是适合的,对于更复杂或更高维的航迹融合却不适合。本文针对2条航迹融合的问题提出了一种基于神经网络补偿的航迹融合方法。该方法不直接计算交互式协方差,而是补偿因交互式协方差而导致的融合误差。并通过跟踪实例,对所提出方法的可行性进行了验证。数学模型所跟踪目标的动态模型为式中X(k)=[x(k
5、)x(k)y(k)y(k)]T为状态向量Φ,G分别为状态矩阵和转移矩阵;输入噪声w(k)为零均值高斯白噪声,其方差。观测值为距离r和角度θ,它们的标准方差分别表示为σr和σθ,将其转化到笛卡尔坐标系下,x=rcosθ和y=rsinθ,这时测量方程为式中,vm(k)为笛卡尔坐标系下的测量噪声,它与过程噪声w(k)是相互独立的,假定其为零均值高斯白噪声且其方差矩阵为:融合系统结构如图1所示。2个传感器跟踪同一目标,融合中心的目的在于将2个跟踪器的输出进行融合而得到更精确的估计值,本文不考虑配准和航迹关联问题,并且假定整个系统是可观测的。目标状态估计可从以下的卡尔曼滤波方程得出:假定从2个传感
6、器来的航迹是同步的,否则在融合之前,这些航迹都必须进行时间配准。不考虑交互式协方差的简单融合方法为其方差矩阵可表示为由于共同过程噪声的存在,各个航迹之间其实是相关的,因此上述方法是次优的。本文将引入神经网络来减少由过程噪声所导致的估计误差。基于神经网络补偿的航迹融合本文提出一种神经网络补偿的方法来提高航迹融合的性能,从而提高跟踪的精确度。该方法的基本结构如图2所示,在融合中心首先对各个传感器的输出进行简单融合,然后引入径向基函数(RBF)神经网络辅助简单融合来减少因交互式协方差或模型变化所引起的估计误差。这里神经网络的训练采用4个输入信号2个为局部航迹的残余误差(εm(k),m=1,2)
7、,2个为传感器的当前状态估计值与当前状态预测值之差(^Xm(k)-Φ^Xm(k-1),m=1,2)。对神经网络输入的选择主要基于以下的考虑:从式(12)可知,交互式协方差主要是同卡尔曼滤波器的增益Km(k)相关,而由式(4)知Km(k)可由滤波器的残差矩阵及^Xm(k)和Φ^Xm(k-1)之间的差值获得,另外神经网络的输入Km(k)和^Xm(k)-Φ^Xm(k-1)都是和过程噪声相关的。因为就仿真而言目标的真实状态是已知的,因此,神