发明名称 一种粒子滤波方法、装置及目标跟踪方法、装置
摘要 本发明公开了一种粒子滤波方法、装置及目标跟踪方法、装置,该粒子滤波方法包括利用高斯-厄米特积分构建上一目标观测时刻的多个积分点概率密度函数;根据多个积分点概率密度函数获取积分点的近似粒子集;根据目标的相关特性修正近似粒子集以获取预测粒子集;根据预测粒子集获取当前目标观测时刻的目标状态预测概率密度函数;根据当前目标观测时刻的目标状态预测概率密度函数获取当前目标观测时刻的目标状态后验概率密度函数。通过上述方式,本发明能够增强粒子的多样性和准确性,有效提高滤波精度以及目标状态的估计性能。
申请公布号 CN103902812B 申请公布日期 2016.05.04
申请号 CN201410079861.3 申请日期 2014.03.05
申请人 深圳大学 发明人 李良群;谢维信;刘宗香
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 深圳市威世博知识产权代理事务所(普通合伙) 44280 代理人 何青瓦
主权项 一种粒子滤波方法,其特征在于,包括:获取k‑1时刻的目标状态后验概率密度函数,具体如下式所示:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>|</mo><msub><mi>z</mi><mrow><mn>1</mn><mo>:</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>N</mi><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>^</mo></mover><msub><mi>x</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>,</mo></mrow></msub></msub><msub><mi>P</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000756382730000011.GIF" wi="1142" he="94" /></maths>其中,k‑1时刻表示上一目标观测时刻,<img file="FDA0000756382730000012.GIF" wi="141" he="78" />为均值,P<sub>k‑1|k‑1</sub>为协方差;估计k‑1时刻的目标状态后验概率密度函数对应的积分点x<sub>l,k‑1|k‑1</sub>,具体如下式所示:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>x</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msqrt><msub><mi>P</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></msqrt><msub><mi>&xi;</mi><mi>l</mi></msub><mo>+</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000756382730000013.GIF" wi="1134" he="84" /></maths>其中,l=1,2,...,m,ξ<sub>l</sub>表示高斯‑厄米特积分点;以k‑1时刻的目标状态后验概率密度函数对应的积分点x<sub>l,k‑1|k‑1</sub>为均值,P<sub>k‑1|k‑1</sub>为协方差,构建k‑1时刻的多个积分点概率密度函数<img file="FDA0000756382730000014.GIF" wi="350" he="87" />根据所述多个积分点概率密度函数获取积分点的近似粒子集;所述积分点x<sub>l,k‑1|k‑1</sub>的近似粒子集为<img file="FDA0000756382730000015.GIF" wi="261" he="103" />其中<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mrow><mo>{</mo><msubsup><mi>x</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>}</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>s</mi></msub></msubsup><mo>~</mo><mi>N</mi><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>l</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>,</mo></mrow></msub><msub><mi>P</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000756382730000016.GIF" wi="637" he="102" /></maths>i=1,2,...,N<sub>s</sub>;获取所述近似粒子集<img file="FDA0000756382730000017.GIF" wi="230" he="101" />中每个近似粒子的期望值<img file="FDA0000756382730000018.GIF" wi="198" he="79" />具体如下式所示:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>x</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>=</mo><msub><mi>f</mi><mi>k</mi></msub><mo>(</mo><mrow><msubsup><mi>x</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>,</mo><msub><mi>v</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow><mo>)</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000756382730000019.GIF" wi="1078" he="95" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>&alpha;</mi><mi>k</mi><mi>i</mi></msubsup><mo>=</mo><mfrac><mrow><msup><mi>T</mi><mn>2</mn></msup><mo>&CenterDot;</mo><msup><mi>v</mi><mn>2</mn></msup><mo>&CenterDot;</mo><msubsup><mi>&sigma;</mi><mi>v</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mi>&lambda;</mi><mo>&CenterDot;</mo><msubsup><mi>&sigma;</mi><mi>m</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>T</mi><mn>2</mn></msup><mo>&CenterDot;</mo><msup><mi>v</mi><mn>2</mn></msup><mo>&CenterDot;</mo><msubsup><mi>&sigma;</mi><mi>v</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00007563827300000110.GIF" wi="1126" he="159" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>&mu;</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi></mrow><mi>i</mi></msubsup><mrow><mo>(</mo><msub><mi>A</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>x</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>+</mo><msubsup><mi>&alpha;</mi><mi>k</mi><mi>i</mi></msubsup><mrow><mo>(</mo><msub><mi>z</mi><mi>k</mi></msub><mo>-</mo><msub><mi>h</mi><mi>k</mi></msub><mo>(</mo><msubsup><mi>x</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00007563827300000111.GIF" wi="1277" he="102" /></maths>其中,所述目标的相关特性包括目标观测、采样时间间隔以及目标速度,z<sub>k</sub>表示k时刻的目标观测,k时刻表示当前目标观测时刻,T表示采样时间间隔,v表示目标速度,λ为常数,σ<sub>m</sub>(k)表示观测误差标准差,σ<sub>v</sub>(k)表示观测新息标准差,f<sub>k</sub>表示系统状态转移函数,v<sub>k‑1</sub>表示状态的过程噪声,h<sub>k</sub>表示非线性的观测函数;获取预测粒子集<img file="FDA0000756382730000021.GIF" wi="215" he="103" />具体如下式所示:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>x</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>=</mo><msubsup><mi>&mu;</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi></mrow><mi>i</mi></msubsup><mrow><mo>(</mo><msub><mi>A</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000756382730000022.GIF" wi="1003" he="79" /></maths>根据所述预测粒子集获取当前目标观测时刻的目标状态预测概率密度函数;根据所述当前目标观测时刻的目标状态预测概率密度函数获取当前目标观测时刻的目标状态后验概率密度函数,完成粒子滤波过程。
地址 518060 广东省深圳市南海大道3688号