CN111211760A - 一种基于分布式扩散策略的反馈粒子滤波方法 - Google Patents
一种基于分布式扩散策略的反馈粒子滤波方法 Download PDFInfo
- Publication number
- CN111211760A CN111211760A CN202010042906.5A CN202010042906A CN111211760A CN 111211760 A CN111211760 A CN 111211760A CN 202010042906 A CN202010042906 A CN 202010042906A CN 111211760 A CN111211760 A CN 111211760A
- Authority
- CN
- China
- Prior art keywords
- particle
- sensor
- distributed
- feedback
- particle filtering
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
Landscapes
- Multi Processors (AREA)
- Complex Calculations (AREA)
Abstract
本发明属于信号处理领域,涉及信号处理领域的目标跟踪问题,具体为一种基于分布式扩散策略的反馈粒子滤波方法。本发明基于分布式网络框架,网络中的每个节点通过与邻居节点交换信息完成对目标的跟踪,降低了网络运算负担,且具有更强的稳健性;同时,本发明通过引入由目标状态的真实后验分布与粒子后验分布之间的Kullback‑Leibler(K‑L)散度定义代价函数,使每个粒子在基于其自身状态和经验分布特征的反馈控制下演化,从而不需要再构造建议分布和进行重采样过程;另外,本发明与经典分布式粒子滤波方法相比,能够实现更高的跟踪精度,更好的跟踪性能,且对粒子数变化敏感度低,稳定性更好。
Description
技术领域
本发明属于信号处理领域,涉及信号处理领域的目标跟踪问题,特别是涉及到分布式架构下的反馈粒子滤波目标跟踪问题,具体为一种基于分布式扩散策略的反馈粒子滤波方法。
背景技术
粒子滤波是一类基于递推贝叶斯法则,利用样本统计量对系统状态进行估计的方法,在处理非高斯非线性系统的参数/状态滤波问题方面具有独特的优势,已被广泛应用于目标跟踪、通信信号处理以及计算机视觉等领域。
一般针对粒子滤波方法研究主要基于文献“Gordon N J,Salmond D J,Smith A FM.Novel Approach to Nonlinear/Non-Gaussian Bayesian State Estimation[J].Radarand Signal Processing,IEE Proceedings F,1993,140(2):107-113.”中提出的自举滤波(Bootstrap Filter)框架,即序列重要性采样(Sequential Importance Sampling,SIS)或建议分布采样-重要性再采样(Sampling-Importance Resampling,SIR),如果建议分布产生的样本不能足量覆盖后验分布样本,则容易导致粒子贫化。针对贫化问题已有大量改进算法,主要针对重采样改进或(和)建议分布改进;重采样改进,如马尔科夫蒙特卡洛(MCMC)采样、辅助粒子滤波(APF)、正则化粒子滤波(RPF)等,目的在于从建议样本中更准确的抽取后验样本,并未从根本上克服贫化问题;建议分布改进,如无迹卡尔曼-粒子滤波(UKF-PF)、容积卡尔曼-粒子滤波(CKF-PF)及混合粒子滤波(MPF)等,目的是构造更接近后验概率密度的建议分布,这类改进方法对于不同系统模型的性能差异较大,难以达到贝叶斯意义下的最优估计。
近年来,文献“Yang T,Laugesen R S,Mehta P G,et al.Multivariablefeedback particle filter[C]2016:10-23.”中提出了一种反馈粒子滤波(FeedbackParticle Filter,FPF)方法,反馈粒子滤波是一种基于平均场博弈论(mean-field gametheory)的近似非线性滤波新方法;通过对先验样本构造反馈,采用Kullback-LeiblerDivergence(KLD)来度量反馈后概率分布与真实后验分布之间的差异,将最小化KLD等效为欧拉-拉格朗日边值问题求解反馈增益,使得由先验样本能够直接得出后验样本,每个粒子在基于其自身状态和经验分布特征的反馈控制下演化,从而不需要再构造建议分布和进行重采样过程。
目前,反馈粒子滤波方法都基于集中式架构,即将网络中各接收机接收的信号都传递到融合中心,在该特定的接收机上进行运算,即集中式的处理方式;然而,集中式处理方式扩展性较差,存在多跳通信问题,对网络的通信带宽要求较高,以及因为存在融合中心接收机和参考信号导致鲁棒性较差;同时,因为所有的位置估计运算都在融合中心接收机进行,所以融合中心接收机的运算负担和能量消耗都很大。
基于此,本发明提出一种基于分布式扩散策略的反馈粒子滤波方法。
发明内容
本发明的目的在于针对上述集中式反馈粒子滤波方法所存在的问题,提出一种基于分布式扩散策略的反馈粒子滤波方法。为实现该目的,本发明采用的技术方案如下。
一种基于分布式扩散策略的反馈粒子滤波方法,其特征在于,包括以下步骤:
其中,
其中,
步骤5:计算k时刻传感器l的中间状态估计值αl,k:
本发明的有益效果在于:
本发明提出的一种基于分布式扩散策略的反馈粒子滤波方法,具有如下优点:
1、本发明具有更强的稳健性;对集中式而言,当融合处理中心出现问题时,整个系统将会瘫痪,而分布式算法能有效避免融合中心出问题时整个系统崩溃的风险;
2、在大型的传感器网络中,带宽限制了邻居传感器之间的信息传输数量,每个节点处只需同其邻点交换信息,本发明提出的基于分布式网络框架的方法不需要网络所有节点的信息送至融合中心进行处理,允许网络中的每个节点通过与邻居节点交换信息完成对目标的跟踪,降低了网络运算负担;
3、本发明提出的方法是基于最优控制和平均场博弈论的概念的新型非线性分布式滤波粒子滤波方法,引入由目标状态的真实后验分布与粒子后验分布之间的Kullback-Leibler(K-L)散度定义的代价函数,最小化该代价函数,得到最优反馈增益函数进而使每个粒子在反馈控制下不断演化,从而不需要再构造建议分布和进行重采样过程;
4、本发明提出的反馈粒子滤波方法与经典分布式粒子滤波方法相比,能够实现更高的跟踪精度,更好的跟踪性能;同时,受粒子数影响更小,稳定性更好。
附图说明
图1为本发明基于分布式扩散策略的反馈粒子滤波方法的流程图。
图2为本发明实施例中传感器网络示意图。
图3为本发明实施例中使用的网络拓扑结构图。
图4~图7为本发明实施例中仿真结果图。
具体实施方式
下面结合附图对本发明作进一步说明:
本发明提供一种基于分布式扩散策略的反馈粒子滤波方法,考虑动态系统的状态空间模型为:
xk+1=a(xk)+ξk,
zk=h(xk)+ηk,
其中,k时刻:(表示xk是实数域内d维列向量、d为目标状态向量的维度)为目标状态向量,为观测向量;映射函数映射函数N为观测向量的维度;分别为系统噪声和观测噪声,均服从互不相关的零均值高斯分布,其协方差矩阵表示为:
其中,(·)T表示矩阵/向量的转置,δki为克罗内克尔记号;
本实施例中,传感器网络中有L个传感器,如图2所示,将传感器l及其邻居传感器(包括传感器l)的集合记为记表示集合的基数、即nl表示传感器l的邻居传感器数目(包含传感器l);即传感器l的邻居传感器集合由传感器2、3、4、l及j组成;基于分布式扩散策略的反馈粒子滤波方法的流程如图2所示,具体包括以下步骤:
其中,
仿真测试:
下面构建一个典型的基于距离的系统模型,将本发明方法(简称分布式反馈粒子滤波)与文献“Mohammadi A,Asif A.Diffusive Particle Filtering for DistributedMultisensor Estimation[J].Acoustics,Speech,and Signal Processing,1988.ICASSP-88.1988International Conference on,2016.”中提出的一种经典的分布式粒子滤波方法(简称分布式粒子滤波)对比来说明本发明方法的可行性、优越性。
仿真条件:
(1)信号模型:采用包含L个节点的分布式网络,网络拓扑结构如图2所示,假定目标在x-y平面移动,假设目标做匀速直线运动,目标状态为目标运动过程中必然受到风力、摩擦力等因素的干扰,这个干扰噪声可以认为随机加速度,采样间隔为T=1,因此状态方程的具体形式如下:
xk+1=Φxk+Γξk,
zl,k=hl(xk)+ηl,k,
(2)参数设置:网络中共有13个传感器,粒子初始估计位置[60,3,60,3];目标初始真实位置[-100,10,100,10];观测噪声方差为100;最大邻居节点数目4;邻域最大半径2000;分布式粒子滤波方法与分布式反馈粒子滤波方法都使用静态的扩散组合系数(Uniform)进行仿真,蒙特卡罗实验50次。
位置NRMSE(Network Root Mean Square Error)定义为:
每个传感器上的粒子数设为50时,分布式粒子滤波方法与分布式反馈粒子滤波方法的位置误差仿真结果如图3所示,二者速度误差的仿真结果如图4所示。
每个传感器上的粒子数设为200时,分布式粒子滤波方法与分布式反馈粒子滤波方法的位置误差仿真结果如图5所示,二者速度误差的仿真结果如图6所示。
由图4、5、6、7表明,本发明方法(即分布式反馈粒子滤波方法)在基于距离的观测模型下,相对于文献“Diffusive Particle Filtering for Distributed MultisensorEstimation”中的分布式粒子滤波方法,位置和速度的NRMSE结果具有更快的收敛速度以及更低的稳态,这意味着本发明方法(即分布式反馈粒子滤波方法)跟踪目标的准确度和速度方面明显优于文献“Diffusive Particle Filtering for Distributed MultisensorEstimation”中的分布式粒子滤波方法。
进一步对比图4、5与图6、7可见,在粒子数设置为50和200时,相对于文献“Diffusive Particle Filtering for Distributed Multisensor Estimation”中的分布式粒子滤波方法,本发明方法(即分布式反馈粒子滤波方法)位置和速度的NRMSE结果差异不大,因此,本发明方法对粒子数变化不敏感,稳定性更强。
以上所述,仅为本发明的具体实施方式,本说明书中所公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换;所公开的所有特征、或所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以任何方式组合。
Claims (1)
1.一种基于分布式扩散策略的反馈粒子滤波方法,其特征在于,包括以下步骤:
其中,
其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010042906.5A CN111211760B (zh) | 2020-01-15 | 2020-01-15 | 一种基于分布式扩散策略的反馈粒子滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010042906.5A CN111211760B (zh) | 2020-01-15 | 2020-01-15 | 一种基于分布式扩散策略的反馈粒子滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111211760A true CN111211760A (zh) | 2020-05-29 |
CN111211760B CN111211760B (zh) | 2023-04-11 |
Family
ID=70790126
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010042906.5A Active CN111211760B (zh) | 2020-01-15 | 2020-01-15 | 一种基于分布式扩散策略的反馈粒子滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111211760B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101162482A (zh) * | 2007-11-20 | 2008-04-16 | 南京邮电大学 | 基于节点协作的高斯和半粒子滤波方法 |
US20080215298A1 (en) * | 2006-10-10 | 2008-09-04 | Haney Philip J | Parameterization of non-linear/non-gaussian data distributions for efficient information sharing in distributed sensor networks |
CN102427341A (zh) * | 2011-11-22 | 2012-04-25 | 上海大学 | 基于卡尔曼滤波的远程迭代学习控制系统传输噪声抑制方法 |
US20130132454A1 (en) * | 2011-11-21 | 2013-05-23 | The Board Of Trustees Of The University Of Illinois | Feedback-Based Particle Filtering |
CN103888100A (zh) * | 2014-03-29 | 2014-06-25 | 北京航空航天大学 | 一种基于负熵的非高斯线性随机系统滤波方法 |
CN107367710A (zh) * | 2017-07-18 | 2017-11-21 | 电子科技大学 | 一种基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法 |
CN109710978A (zh) * | 2018-11-30 | 2019-05-03 | 电子科技大学 | 一种分布式异构自适应粒子滤波直接跟踪方法 |
CN110516198A (zh) * | 2019-07-17 | 2019-11-29 | 电子科技大学 | 一种分布式非线性卡尔曼滤波方法 |
CN110649911A (zh) * | 2019-07-17 | 2020-01-03 | 电子科技大学 | 一种基于α散度的分布式非线性卡尔曼滤波方法 |
-
2020
- 2020-01-15 CN CN202010042906.5A patent/CN111211760B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080215298A1 (en) * | 2006-10-10 | 2008-09-04 | Haney Philip J | Parameterization of non-linear/non-gaussian data distributions for efficient information sharing in distributed sensor networks |
CN101162482A (zh) * | 2007-11-20 | 2008-04-16 | 南京邮电大学 | 基于节点协作的高斯和半粒子滤波方法 |
US20130132454A1 (en) * | 2011-11-21 | 2013-05-23 | The Board Of Trustees Of The University Of Illinois | Feedback-Based Particle Filtering |
CN102427341A (zh) * | 2011-11-22 | 2012-04-25 | 上海大学 | 基于卡尔曼滤波的远程迭代学习控制系统传输噪声抑制方法 |
CN103888100A (zh) * | 2014-03-29 | 2014-06-25 | 北京航空航天大学 | 一种基于负熵的非高斯线性随机系统滤波方法 |
CN107367710A (zh) * | 2017-07-18 | 2017-11-21 | 电子科技大学 | 一种基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法 |
CN109710978A (zh) * | 2018-11-30 | 2019-05-03 | 电子科技大学 | 一种分布式异构自适应粒子滤波直接跟踪方法 |
CN110516198A (zh) * | 2019-07-17 | 2019-11-29 | 电子科技大学 | 一种分布式非线性卡尔曼滤波方法 |
CN110649911A (zh) * | 2019-07-17 | 2020-01-03 | 电子科技大学 | 一种基于α散度的分布式非线性卡尔曼滤波方法 |
Non-Patent Citations (2)
Title |
---|
XIAOLIN NING;ZHUO LI;WEIREN WU;YUQING YANG;JIANCHENG FANG;GANG LIU;: "Recursive adaptive filter using current innovation for celestial navigation during the Mars approach phase" * |
陈明方;刘金鑫;张永霞;葛天佑;臧家秀;陈浩;: "基于S7-1200AD采样的高效数字滤波算法的设计与实践" * |
Also Published As
Publication number | Publication date |
---|---|
CN111211760B (zh) | 2023-04-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yuan et al. | TOA-based passive localization constructed over factor graphs: A unified framework | |
Mohammadi et al. | Distributed particle filter implementation with intermittent/irregular consensus convergence | |
CN111985093A (zh) | 一种带噪声估计器的自适应无迹卡尔曼滤波状态估计方法 | |
CN109710978B (zh) | 一种分布式异构自适应粒子滤波直接跟踪方法 | |
CN109827579B (zh) | 一种组合定位中滤波模型实时校正的方法和系统 | |
CN111257865B (zh) | 一种基于线性伪量测模型的机动目标多帧检测跟踪方法 | |
CN112881979A (zh) | 一种基于ekf滤波的初始状态自适应融合定位方法 | |
CN109671100A (zh) | 一种分布式可变扩散组合系数粒子滤波直接跟踪方法 | |
CN110426671B (zh) | Wsn中基于模型概率实时修正的imm目标跟踪方法及装置 | |
Zhang et al. | High‐degree cubature Kalman filter for nonlinear state estimation with missing measurements | |
Nannuru et al. | Hybrid multi-Bernoulli and CPHD filters for superpositional sensors | |
Ge et al. | Adaptive cubature strong tracking information filter using variational Bayesian method | |
CN113543026A (zh) | 一种基于径向基函数网络的多楼层室内定位方法 | |
CN111211760B (zh) | 一种基于分布式扩散策略的反馈粒子滤波方法 | |
Zhao et al. | Factor graph based multi-source data fusion for wireless localization | |
Zhang et al. | Distributed adaptive cubature information filtering for bounded noise system in wireless sensor networks | |
Ji et al. | Data-driven distributed information-weighted consensus filtering in discrete-time sensor networks with switching topologies | |
Li et al. | Multi-sensor measurement fusion based on minimum mixture error entropy with non-Gaussian measurement noise | |
Ge et al. | Credible Gaussian sum cubature Kalman filter based on non-Gaussian characteristic analysis | |
Fu et al. | Distributed maximum correntropy Kalman filter with state equality constraints in a sensor network with packet drops | |
Fu et al. | Distributed state estimation with state equality constraints in the presence of packet dropping | |
CN115114985A (zh) | 一种基于集合理论的传感器系统分布式融合方法 | |
CN109581280A (zh) | 终端自适应在线定位方法、系统及装置 | |
Li et al. | Joint sequential target estimation and clock synchronization in wireless sensor networks | |
Huang et al. | VariFi: Variational Inference for Indoor Pedestrian Localization and Tracking Using IMU and WiFi RSS |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |