发明名称 基于粒子群优化算法的无迹粒子滤波方法
摘要 本发明提供的是一种基于粒子群优化算法的无迹粒子滤波方法。第一步,初始时刻,由初始分布p(x0)中得到一组初始粒子,并设置其初始的均值和方差;第二步,重要性采样;第三步,权值更新;第四步,得到归一化的权值;第五步,再采样;第六步,状态更新;第七步,求解当前时刻的全局最优解G(t)。本发明通过粒子群优化过程,使得粒子集在权值更新之前更加趋向于高似然区域,从而在一定程度上解决了样本贫化问题。优化过程使得远离真实状态的粒子趋向真实状态出现概率较大的区域,提高了每个粒子的作用效果。与其他智能优化算法相比,粒子群优化算法具有容易实现且没有许多参数需要调整的优势,并且降低了精确估计所需要的粒子数,提高了滤波的计算效率。
申请公布号 CN101826852A 申请公布日期 2010.09.08
申请号 CN201010121571.2 申请日期 2010.03.11
申请人 哈尔滨工程大学 发明人 杨萌;高伟;郝燕玲
分类号 H03H17/00(2006.01)I;G06N3/12(2006.01)I 主分类号 H03H17/00(2006.01)I
代理机构 代理人
主权项 1.一种基于粒子群优化算法的无迹粒子滤波方法,其特征是主要包括如下步骤:第一步,初始时刻,由初始分布p(x<sub>0</sub>)中得到一组初始粒子,并设置其初始的均值和方差,设置粒子初速度,设置惯性因数w,学习因子c<sub>1</sub>和c<sub>2</sub>,求解初始时刻的全局最优解G(0);第二步,重要性采样(1)调整粒子的速度和位置V<sub>i</sub>(t)=w×V<sub>i</sub>(t-1)+c<sub>2</sub>×r<sub>2</sub>×[G(t-1)-X<sub>i</sub>(t-1)],其中V<sub>i</sub>(t)为第i个粒子在t时刻的速度,X<sub>i</sub>(t)为其位置;(2)根据无迹粒子滤波算法对粒子的状态进行更新,在该时刻的状态估计中引入粒子速度值,即X<sub>i</sub>(t)=X<sub>i</sub>(t-1)+V<sub>i</sub>(t),(3)求粒子集的均值<img file="FSA00000055063800011.GIF" wi="40" he="56" />和方差P<sub>k</sub><sup>i</sup>,(4)从重要性密度函数中抽取粒子;第三步,权值更新<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>w</mi><mi>k</mi><mi>i</mi></msubsup><mo>=</mo><msubsup><mi>w</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mfrac><mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>k</mi></msub><mo>|</mo><msubsup><mi>x</mi><mi>k</mi><mi>i</mi></msubsup><mo>)</mo></mrow><mi>p</mi><mrow><mo>(</mo><msubsup><mi>x</mi><mi>k</mi><mi>i</mi></msubsup><mo>|</mo><msubsup><mi>x</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>)</mo></mrow></mrow><mrow><mi>&pi;</mi><mrow><mo>(</mo><msubsup><mi>x</mi><mi>k</mi><mi>i</mi></msubsup><mo>|</mo><msubsup><mi>x</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>,</mo><msub><mi>y</mi><mrow><mn>1</mn><mo>:</mo><mi>k</mi></mrow></msub><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>上式中,w<sub>k</sub><sup>i</sup>为粒子对应的权值,p(·|·)为概率密度函数,重要性概率密度函数<img file="FSA00000055063800013.GIF" wi="560" he="74" /><img file="FSA00000055063800014.GIF" wi="198" he="74" />为由均值和方差得到的建议分布;第四步,得到归一化的权值<maths num="0002"><![CDATA[<math><mrow><msubsup><mover><mi>w</mi><mo>~</mo></mover><mi>k</mi><mi>i</mi></msubsup><mo>=</mo><mfrac><msub><mover><mi>w</mi><mi>i</mi></mover><mi>k</mi></msub><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>w</mi><mi>k</mi><mi>j</mi></msubsup></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>第五步,再采样定义阈值N<sub>eff</sub>来衡量有效粒子数量N<sub>thr</sub>为有效粒子数,如N<sub>eff</sub><N<sub>thr</sub>,即对粒子集<img file="FSA00000055063800016.GIF" wi="327" he="75" />重新采样,产生新的集合{x<sub>k</sub><sup>i</sup>,i=1,...,N},重新设定粒子的权值为:<img file="FSA00000055063800021.GIF" wi="216" he="56" />N为粒子数;第六步,状态更新:<maths num="0003"><![CDATA[<math><mrow><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>x</mi><mi>k</mi><mi>i</mi></msubsup><msubsup><mi>w</mi><mi>k</mi><mi>i</mi></msubsup><mo>;</mo></mrow></math>]]></maths>第七步,求解当前时刻的全局最优解G(t)。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公办公室