[发明专利]人工鱼群粒子滤波方法有效
申请号: | 201110308010.8 | 申请日: | 2011-10-11 |
公开(公告)号: | CN102360451A | 公开(公告)日: | 2012-02-22 |
发明(设计)人: | 朱志宇;李阳;李冀;张冰;刘维亭;魏海峰;赵强;袁文华 | 申请(专利权)人: | 江苏科技大学 |
主分类号: | G06N3/00 | 分类号: | G06N3/00 |
代理公司: | 南京经纬专利商标代理有限公司 32200 | 代理人: | 楼高潮 |
地址: | 212003*** | 国省代码: | 江苏;32 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明人工鱼群粒子滤波方法,将人工鱼群算法引入粒子滤波,其基本思路是在选取目标函数即量测函数和调整重要性权值时,引入人工鱼群算法的聚群行为和觅食行为,引导先验粒子不断向高似然域移动,从而改善粒子分布,增加粒子的多样性,提高算法的滤波精度。此外,在人工鱼群算法中,当寻优的域较大或处于变化平坦的区域时,一部分人工鱼将处于无目的的随机移动中,影响了寻优效率,本发明进一步提出了一种自适应步长方法对人工鱼视野选取的随机性进行改进,不仅减轻了算法的计算负担,而且保证了算法的收敛性。 | ||
搜索关键词: | 人工 鱼群 粒子 滤波 方法 | ||
【主权项】:
1.一种人工鱼群粒子滤波方法,其特征在于包括如下步骤:Step1:取得量测值Z k : [ ( 2 π ) σ v 2 ] - 1 / 2 e - 1 2 σ v 2 [ ( z k - z ^ k / k - 1 i ) 2 ] ]]> 其中Zk为最新量测值,
为预测量测值,k表示时刻;Step2:初始化在k=0时刻,从重要性函数采样取N个粒子,抽样出的粒子用
表示,重要性密度函数取转移先验:X k i : q ( X k i | X k - 1 i , Z k ) = p ( X k i | X k - 1 i ) ]]>X k i = ( X ^ k / k - 1 1 , X ^ k / k - 1 2 , · · · , X ^ k / k - 1 n ) ]]> Step3:重要性权值计算w k i = w k - 1 i p ( Z k | X k - 1 i ) = w k - 1 i p ( Z k | X k i ) p ( X k i | X k - 1 i ) q ( X k i | X k - 1 i , Z k ) ]]>= w k - 1 i p ( Z k | X k i ) = w k - 1 i [ ( 2 π ) σ v 2 ] - 1 / 2 e - 1 2 σ v 2 [ ( z k - z ^ k / k - 1 i ) 2 ] ]]> 令目标函数Y = [ ( 2 π ) σ v 2 ] - 1 / 2 e - 1 2 σ v 2 [ ( z k - z ^ k / k - 1 i ) 2 ] ; ]]> (1)觅食行为当Yi<Yj时,粒子通过比较目标函数不断更新自己向更真实的状态靠近:X ^ k / k - 1 i m = X ^ k / k - 1 i m - 1 + r · s X ^ k / k - 1 j - X ^ k / k - 1 i m - 1 | | X ^ k / k - 1 j - X ^ k / k - 1 i m - 1 | | · | 1 - Y j Y i | ]]> (求解极小值)其中
为对
第m∈{1,2,…,n}次后的迭代值,否则X ^ k / k - 1 i m = X ^ k / k - 1 i m - 1 + r · s · X ^ k / k - 1 j - X ^ k / k - 1 i m - 1 | | X ^ k / k - 1 j - X ^ k / k - 1 i m - 1 | | · | 1 - Y i Y j | ]]> (求解极大值)(2)聚群行为设
在范围v内的其他预测观测值集合为Av,|Av|为集合Av的势,则其中心为Z ^ k / k - 1 ic = 1 | A v | Σ j ≠ i | A v | Z ^ k / k - 1 j , ]]> 若Y ic | A v | > δ · Y i , ]]> 则X ^ k / k - 1 i m = X ^ k / k - 1 i m - 1 + r · s X ^ k / k - 1 ic - X ^ k / k - 1 i m - 1 | | X ^ k / k - 1 ic - X ^ k / k - 1 i m - 1 | | , ]]> 否则转(1);Step4:权值更新并归一化权值根据w k i = w k / k - 1 i p ( Z k | X k i ) ]]> 更新权值;根据
对取出的N1个粒子的权值归一化;Step5:输出状态估计:X ^ k = Σ i = 1 N w k i X k i ; ]]> 方差估计:P k = Σ i = 1 N w k i ( X k i - X ^ k ) ( X k i - X ^ k ) T ; ]]> Step6:判断是否结束,当达到设置的迭代次数n或者目标函数前后偏差小于设定的阈值则退出本算法;否则返回step2。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于江苏科技大学,未经江苏科技大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201110308010.8/,转载请声明来源钻瓜专利网。