遗传重采样粒子滤波器
第33卷第8期
2007年8月
自动化
ACTAAUTOMATICASINICA
Vo1.33,No.8
August,2007
遗传重采样粒子滤波器
叶龙王京玲张勤
摘要粒子退化现象是影响粒子滤波器性能的一个重要因素.本文针
对粒子退化,将遗传机制应片j于粒子重采样,以进化
设计
领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计
解决退化问题.
分析并给出了平衡粒子集的有效性与多样性的手段以取得最佳性能的遗
传粒子滤波结构的方法.
关键词粒子滤波器,粒子退化,遗传机制,有效性粒子,多样性粒子
中图分类号TP3
GeneticResamplingParticleFilter
YELongWANGJing—LingZHANGQin
AbstractParticledegenerationisakeyissueintheperfor—
manceofaparticlefilter.Inthispaperweintroducegenetic mechanismintoparticleresamplingprocess.Itisshownthatthe newparticlefiltercaneffectivelyeliminateparticledegeneration andreduceitsdependencyontheparticlevalidity.Furthermore, thenewgeneticparticlefiltercanbeoptimizedbythreekeyge—
neticfactorsselection,crossoverandmutationprobabilities. KeywordsParticlefilter,particledegeneration,genetic mechanism,effectiveparticle,diversiformparticle 1引言
近年米,针对状态估计与运动跟踪,粒子滤波
器[h3_(Particlefilter)由于采用蒙特卡罗采样(MonteCarlo sampling)结构存非线性,非高斯系统状态跟踪上体现出 越来越大的优越性,并得到了广泛的应用[,.粒子退化[.】 是粒子滤波器中不可避免的,同时严重影响粒子滤波器性 能的现象,在粒子滤波器经过几次迭代之后,很多粒子只有 很小甚至接近于零的权值,这些权值存进行粒子的重要度更 新的时候虽然还要计算但是对整个系统的帮助很小,基奉 上属于无用的粒子,这样一方面浪费丫大量计算资源,同时 也容易造成跟踪的精度降低甚至日标丢失.粒子重采样是解 决粒子退化问题的一种重要方法,常用的重采样算法有累 积分布重采样[2](Binarysearch),系统重采样[.](Systematic
resampling),剩余重采样[](Residualresampling)等.这些 算法通过增加粒子的有效性解决了粒子退化问题,但是在实 际应用中可影响系统的鲁棒性.重采样完成后,重要度高的 粒子通过重采样被多次选取,这在一定程度上丢失了粒子的 多样性,由此造成的后果就是一旦目标丢失或跟踪精度不够, 系统自动收敛的可能性很小.
针对这一问题,本文将遗传机制【8_应用于粒子重采样, 利用进化思想解决粒子退化问题,即针对重采样中粒子有效 性与多样性的两个矛盾,提出了一种新型的在保证粒子有效 性的I时科学地增加粒子多样性的方法一遗传重采样算法. 1刊时通过单变量与多变量跟踪问题的实例,给出了针对具体 问题寻找最佳遗传系统的
方案
气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载
,通过引入置信区间分析并演 收稿日期2006—2—16收修改稿日期2007-4-19 ReceivedFebruary16,2006;inrevisedformApril19,2007
国家自然科学基金(60572041)资助
SupportedbyNationalNaturalScienceFoundationofChina
C60572041)
1.中国传媒大学信息工程学院北京100024 1.InformationEngineeringSchool,CommunicationUniversityof
China,Beijing100024
DoI:10.1360/aas-007-0885 示r多样性粒子与有效性粒子的比例对于跟踪准确性的影 响.
2遗传重采样算法与遗传粒子滤波器
假设一个动态状态空间中的目标状态方程与观测方程定 义如下
k=,(k一1,Vk)(1)
zk=h(xk,Uk)(2)
其中,zk与zk分别表示目标状态值与观测值,A:R× R,一R表示日标状态非线性转移函数,k为状态转移 噪声,n与礼分别表示目标状态矢量维数与状态噪声矢量 维数;hk:R×R,一R表示目标状态非线性观测函 数,uk为观测噪声,礼与礼分别表示目标观测矢量维数与 观测噪声矢量维数,跟踪的目的就是通过目标的观测状态zk 得到k的估计.
粒子滤波器是通过递归蒙特卡罗采样实现跟踪的一种统 计计算方法,其算法包括两个主要步骤:预测和更新.预测是 指根据Chapman—Kolmogorov方程得到目标运动的先验概 率密度函数(Probabilitydensityfunction,PDF);向更新是 指通过叭叶斯
公式
小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载
与目标状态似然度对先验PDF做出的修 正.在粒子滤波器中,PDF用带有权值的粒子组成的粒子集 近似表示,因此粒子滤波器是一种通过粒子状态与权值的预 测与更新而实现目标状态后验概率密度估计的方法.粒子退 化是粒子滤波器中一个不可避免的问题.在粒子状态更新后, 为解决粒子退化问题,需要对粒子集合进行重采样操作以去 除不重要的粒子,但是这种去除往往使粒子集丢失其多样性
(如引言中所述).因此,在本文中,我们将遗传操作引入了粒 子重采样,解决r这一问题.
将遗传机制应用于粒子的重采样.首先,作为一种进化 思想的理论基础,遗传机制为解决粒子退化问题提供J,重要 的指导思想;其次,遗传机制不仅仅通过选择算子【8_遴选 优良个体,还可以通过重组算子与变异算子操作产 生新的个体.凶此,适当地优化调整选择概普蔓,重组概率 Pc以及变异概率,可以在保证粒子有效性的时兼顾粒 子的多样性.遗传重采样过程是遗传机制的执行过程,即针 对每一个带有权值的粒子状态,首先对状态进行二进制编码, 然后按照设定的选择,交叉,变异概率对于粒子集依次进行 相虑的算子计算,得到的粒子集合在二进制解码后得到最终 的重采样粒子集.因此,遗传重采样下的粒子滤波器算法可 以描述为算法1.
算法1.遗传重采样粒子滤波器算法(GRPF) [{,ki,tN:1]=GRPF[{x~~一1,一1)v_1,zk] 步骤1.初始化粒子集忙&,6}:N fork=1:K
步骤2.粒子状态预测
P(,一1)i=1:?
步骤3.粒子状态权衡(更新)
o(一1P(zklx~)
:
?
i=1
步骤4.遗传粒子重采样
自动化33卷
dowhilei<N//选择操作
7"8:S一//表示说明:s表示原粒子,R表示重
采样粒子集,?表示粒子总数,此操作表示从?个粒子的粒子集中通过
选择算子得到新粒子集中的第J个粒子.(下同) 歹=歹+l
enddo
dowhilei<NPc/2//交义操作
:S=Rj:j+l
J=J+2;
enddo
dowhile<NPm//变异操作
:S__+Rf
J=J+l:
enddo
{似t}1,2.,N=百1
end
注.奉算法中,在进行粒子权值更新时假设重要度函数 取目标运动的前向概率密度函数.
3遗传粒子滤波器的性能分析与优化方法
奉文选用粒子滤波器性能评测时常用的两个戍用【】进 行遗传粒子滤波器的性能分析,即单变量的非平稳经济学变 化估计问题与多变量Bearings-only跟踪问题. 1)单变量经济学变化估计
该问题的运动模型与观测模型按式(3)与(4)得出. ()++8cos(1_2(t
:t)
(3)
其中,(t)与(t)为零均值高斯噪声,方差分别为10与1, 初始状态取0:0.1.本文中对于粒子状态进行[-2o201区 间内位数为9的二进制编码.跟踪输出采用最大后验概率密
度输出法(Maximumaposterior,MAP)[91.
我们参照文献[6】中的方法设定后验概率密度分布中 2.5%至97.5%的间隔范同为置信度95%的置信区间,高置信 水平的代价是宽的置信区间,有效性粒子与多样性粒子数目 分配的最佳即通过保证置信区间最有效的包含目标状态的真 实值而使目标估计具有最小的误差.
本文选择采用累积分布重采样,P.=0.6与P=0.8的 遗传重采样的三种形式的粒子滤波器对J—J一状态变化进行估 计,粒子个数取500,跟踪结果如图1.
由1可以看出,与累计分布重采样相比,只:0.8的
遗传重采样算法通过增加粒子集中多样性粒子数目扩大估 计的置信区间,存使置信区间更好地包含状态真实值的时 也使跟踪准确性有所增加.但当只=0.6时,多样性粒子数 目的增加虽然扩大丫估计的置信区间,但是导致准确性降低. 这说明,在有效性粒子与多样性粒子之间存在着一个使跟踪 性能最优的最佳分配问题.
为研究这一问题,我们针对10个选择概率取值分析,跟 踪误差取10次踪过程的误差平均值.如图2(a)所示,选择 概率取值0.8时跟踪性能最好,此时粒子集中新旧粒子数之 比为8:2.需要说明的一点是,0.8只是在奉仿真中特定条件 下的最佳值,存不的应用场合下,最佳选择概率会随预测 模型的准确性的不,预测模型与观测模型中信噪比的不同 而相应变化.
(a)累积分布重采样
(b)=0.8GRPF
(c)=O.6GRPF
图l跟踪效果的置信区间分析图
Fig.1Believingrangeanalysisoftracking
Chooserote
图2(a)遗传重采样粒子滤波器中跟踪误差一选择概率曲线 Fig.2(a)Trackingerror—P8curveofGRPF lnteration
图2(b)各种重采样算法的性能比较图
Fig.2(b)PerformancecomparisonbetweenGRPFandother
resamplingalgorithms
8期叶龙等:遗传重采样粒子滤波器887
21Bearings—only模型
Bearings—only运动模型与观测模型按式(5) x(t)=(t一1)+w(t)
与(6)得出
(5)
其中五=(X,,Y,Vy)t,妒为状态转移矩阵,(f)为
(0,0.001)高斯噪声.
(t)=tan一(Y/x)+(t)(6)
其中v(t)为f0,0.005)高斯噪声.状态初始值与转移矩阵按 文献『61中仿真方案取值.
与单变量估值问题不,对于高维变量使用遗传重采样 算法需要解决状态编码的问题,在Bearings—only模型中,V 与V的状态只与自身前一时刻的状态相关,而X与Y状态 除与自身状态相关外,还与V与V状态相关,变化复杂度 相对较高,并且观测值的确定也只与X和Y有关.因此在进 行遗传重采样的交叉与变异操作时,对于需要操作的父代个 体,我们只针对其状态中的X与Y进行操作.跟踪过程中,选 择概率=0.8,总粒子个数1000.跟踪结果如图3所示. 图3(a)遗传重采样粒子滤波器高维跟踪效果图
Fig.3(a)GRPFsperformanceinmulti—dimensionstrackingarea
4结论
图3(b)各状态变量跟踪曲线
Fig.3(b)Trackingcurvesofstateelement
本文针对粒子滤波器中的退化现象,通过改进重采样算
法,设计出遗传粒子滤波器结构,以进化设计解决_厂退化问
题.给出一种存保证粒子有效性的剧时科学地增加粒子多
样性的算法,解决J,粒子滤波器跟踪问题中粒子有效性与多
样性的矛盾问题,提高了跟踪的准确度并增加跟踪的鲁棒
性.针对粒子有效性与多样性两种矛盾的概念,奉文还通过
仿真实例给出r遗传粒子滤波器优化设计方法,使遗传粒子
滤波器能够达到最佳跟踪与最良好的系统鲁棒性.
References
1HuangAJ.AtutorialonBayesianestimationandtrack- ingtechniquesapplicabletonon—linearandnon—Gaussian
process【Online],available:;papers-05/05—0211/05—0211.pdf,February
11,2005
2DoucetA,GodsillS,ChistopheA.OnsequentialMonte CarlosamplingmethodsforBayesianfiltering.Statisticsand Computing,2000,10(3):197~208
3IsardM,BlakeA.Condensation—conditionaldensitypropa-
gationforvisualtracking.InternationalJournalofComputer Vision,1998,29(1):5^一28
4ChoJU,JinSH,PhamXD,JeonJW,ByunJE,Kang
H.Areal—timeobjecttrackingsystemusingaparticlefilter. In:ProceedingsofIEEE/RSJInternationalConferenceon IntelligentRobotsandSystems.IEEE,2006.2822~2827 5HaykinS,HuberK,ChenZ.Bayesiansequentialstatees- timationformimowirelesscommunications.Proceedingsof the?EE,2004,92(3):439~454
6GordonN,SalmondDJ,SmithAFM.Novelapproach
tononlinear/non—GaussianBayesianstateestimation.忸E
ProceedingsFRadarandSignalProcessing,1993,140(2): 107^一113
7LiuJS,ChenR,LogvinenkoT.Atheoreticalframeworkfor sequentialimportancesamplingandresampling.Sequential MonteCarloinPractice.NewYork:Springer-Verlag,2001. 225^一246
8ZhangWen—Xiu,LiangYi.MathematicalFoundationofGe- neticAlgorithms.Xian:XianJiaotongUniversityPress, 2000.15^一45
(张文修,梁怡.遗传算法的数学基础.两安:西安交通大学出版社
2000.15~45)
9DoucetA,GodsillSJ,WestM.MonteCarlofilteringand smoothingwithapplicationtotime-varyingspectralesti- marion.In:ProceedingsofIEEEInternationalConference onAcoustics,Speech,andSignalProcessing.IEEE,2000. 1701^一1704
叶龙助理研究员.主要研究方向为计算机视觉,统计计算和数学建模.
E-mail:yelong~cuc.edu.an
(YELongLecturer.HisresearchinterestcoveFScomputer vision,statisticalmodeling,andcomputing.) 王京玲教授.主要研究方向为计算机视觉,统计计算和通信.本文通信
作者.E-mail:wjl@CUC.edu.an
(WANGJing-LingProfessor.HerresearchinterestCOVETS computervision,statisticalmodeling,andtelecommunications. c0rresp0ndingauthorofthispaper.)
张勤教授.主要研究方向为多媒体信号处理和通信.
E—mail:zhangqin~cuc.edu.an
(ZHANGQinProfessor.HisresearchinterestCOVESmulti—
mediasignalprocessingandtelecommunications.)