CN107367710B - 一种基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法 - Google Patents
一种基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法 Download PDFInfo
- Publication number
- CN107367710B CN107367710B CN201710584073.3A CN201710584073A CN107367710B CN 107367710 B CN107367710 B CN 107367710B CN 201710584073 A CN201710584073 A CN 201710584073A CN 107367710 B CN107367710 B CN 107367710B
- Authority
- CN
- China
- Prior art keywords
- receiver
- local
- neighbor
- positioning method
- time
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/06—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
Description
技术领域
本发明属于信号处理领域,特别是涉及基于时差和频差的无源定位系统中分布式自适应跟踪定位方法,具体为一种基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法。
背景技术
目前,基于时差的无源定位技术根据是否需要显式计算时差值而分为两大类:经典的两步定位方法和新兴的直接定位方法;两步定位方法首先第一步通过接收信号估计出时差值或频差,然后第二步利用估计出的时差值进行位置解算;然而由于两步算法在第一步估计时差或频差时没有考虑所有单目标定位中的约束问题,所以不能确保定位的准确性;对于移动目标来说,两步算法就更难准确的定位。
而直接定位算法则不需要显式地计算出时差值,而是直接利用接收信号估计目标的位置;在接收信号信噪比较低的情况下,直接定位方法的定位精度更高,直接定位方法又可分为批处理方法和自适应方法两类。对移动的目标需要进行跟踪定位,批处理方法需要对定位区域进行四维的网格式搜索,计算量很大,实时性差;后来提出了用代数方法来估计目标的位置和速度,但是这种方法只有在噪声较小的情况下定位准确;进一步用估计出来的目标的位置和速度作为卡尔曼滤波(KF)的输入,构造动态模型,提升了定位的精确度;后来人们又提出了extended KF(EKF),直接用时差或频差作为输入即可;但是由于KF和EKF主要解决线性问题,而在实际情况中,目标的运动轨迹大部分都是非线性的,现在人们又提出了粒子滤波算法,用来解决非线性的跟踪定位问题。
而粒子滤波算法中又存在粒子退化问题,为了解决这个问题,又提出了重采样(SIR)粒子滤波算法和高斯粒子滤波算法;这些粒子滤波算法又分为集中式的和分布式的。但是,目前基于时差和多普勒频移的自适应直接定位方法都是集中式的架构,即各接收机都将其接收的信号传递到定位运算接收机,在该接收机上进行定位运算,即集中式的处理方式。如文献《Delay and Doppler Induced Direct Tracking by Particle Filter》(Sidi,A.Y.,Weiss,A.J.) 就是采取的该方法;然而,这种集中式处理方式却存在着很大的缺陷;首先,由于所有的运算都集中在一台定位运算接收机上,所以当定位系统中接收机的个数较多时,定位运算接收机将承受很大的计算压力,其能量消耗也将很大,且随着定位系统中接收机个数的增加这种问题将越来越严重;其次,要将所有的信号都传输给定位运算接收机必然会存在多跳传输及传输路径选择等诸多问题,并且在多跳传输的过程中也必然存在着噪声干扰等情况;另外,靠近定位运算接收机的各接收机由于承担着很大的转发接收信号给定位接收机的任务,能量消耗也将变的很大;最后,将所有的定位运算都集中到一个接收机上,对于资源的利用也不是太合理,并且如果在定位的过程中定位运算接收机出现故障,那么整个的定位系统也就瘫痪了,即系统的稳健性不高;且在集中式方法中,会选取参考接收机的信号作为参考信号,定位的精度受该参考接收机信号的影响很大,导致其稳定性很差。
发明内容
本发明的目的在于提供一种基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法,用于进一步提高分布式自适应粒子滤波跟踪定位方法的性能。
为实现上述目的,本发明采用的技术方案为:
一种基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法,其特征在于,包括以下步骤:
步骤1:在接收机个数为N的无线传感器网络当中,局部GPF(Gaussian ParticleFilteri ng)从局部重要函数中随机采M个随机预测粒子,k时刻具体为:
步骤:2:第一次数据交换,接收机l接收信号传输给自己的邻居接收机,同时接收邻居接收机传来的接收信号;
步骤3:接收机l更新每个预测粒子所对应的权重值:
rl,k=hk(xk)+ul,k=bl,kAl,kFl,ksk+ul,k
其中,ul,k为零均值高斯白噪声,bl,k表示信号衰减,Aj,k和Fj,k分别表示多普勒频移和时延;
其中,C′和C″分别是两个常数,且:
不失一般性,||sk||2=1;
因为||rj,k||2和剩余其他参量相互独立,转换为求最大值问题:
其中,定义Ql,k=Vl,k(Vl,k)H、Ql,k为N×N的赫尔米特矩阵(Hermitian matrix),Vl,k=[(Fj,k)H(Aj,k)Hrj,k],j∈Nl、Vl,k为N×L维矩阵、L为集合中元素数量;
上述求最大值的问题等价于求Ql,k的最大特征值λmax(Ql,k),然而Ql,k的维数随着采样数据的增大而增大,这样会造成很大的计算负担,但是Ql,k的最大特征值和L×L矩阵的的大小一样,当L<<N时,这样的替代将会减小计算负担;因此最终将求最大似然估计的问题转化为求最大特征值的问题:
最后对粒子权重进行归一化,因此求得最大特征值即可计算得到每个粒子的归一化权重:
步骤4:接收机l更新它的局部中间状态估计值ψl,k以及相应的协方差Pl,k:
步骤5:第二次数据交换,接收机l将步骤4 计算得到的局部中间状态估计值ψl,k传输给自己的邻居接收机,同时接收邻居接收机传来的结果;
条件概率密度函数为:xk+1=φxk+nk,φ为转移矩阵,nk为零均值高斯白噪声;
本发明有益效果在于:
在本发明中分布式粒子滤波算法的推导过程为,提出新的局部似然函数,推广得到分布式粒子滤波算法。具体推导如下:
根据Diffusion LMS相关理论的启发,我们提出了一种基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法。在分布式算法中,每个接收机都将参与定位运算,在接收机 l(l=1,2,3,…,N)上做如下迭代运算:
从(1)(2)式中我们发现,每一个接收机的运算只需要其邻居接收机的信号,这样就避免了集中式算法中的多跳传输。
公式(1)中的计算细节描述如下:
本发明根据信号模型提出新的局部似然函数为:
其中,C′和C″分别是两个常数,且:
不失一般性,||sk||2=1;
因为||rj,k||2和剩余其他参量相互独立,转换为求最大值问题:
其中,定义Ql,k=Vl,k(Vl,k)H、Ql,k为N×N的赫尔米特矩阵(Hermitian matrix),Vl,k=[(Fj,k)H(Aj,k)Hrj,k],j∈Nl、Vl,k为N×L维矩阵、L为集合中元素数量;
上述求最大值的问题等价于求Ql,k的最大特征值λmax(Ql,k),然而Ql,k的维数随着采样数据的增大而增大,这样会造成很大的计算负担,但是Ql,k的最大特征值和L×L矩阵的的大小一样,当L<<N时,这样的替代将会减小计算负担;因此最终将求最大似然估计的问题转化为求最大特征值的问题:
利用该局部似然函数求出每个接收机上的预测粒子权重:
最后对粒子权重进行归一化,因此求得最大特征值即可计算得到每个粒子的归一化权重:
本发明分布式自适应直接定位方法与集中式自适应直接定位方法的区别及优势如下
(1)本发明分布式定位方法中,每个接收机都参与运算;避免了集中式定位方法中,只有一个接收机进行定位运算;并且每个接收机的运算量相对于集中式方法中的定位运算接收机的运算量要小很多,且在分布式算法中,每个接收机的运算量并不会随着接收机个数的增加而有明显的增加,可扩展性更强。在集中式方法当中,频移Aj,k和时延Fj,k是恒定不变的,但是在分布式方法当中,频移Aj,k和时延Fj,k的形式会发生变化的;
(2)本发明分布式定位方法中,每个接收机都只与自身直接相连的接收机(即邻居接收机)进行数据交互,即只存在单跳传输;避免了集中式定位方法中信号的多跳传输信号,信号在传递过程中受到的干扰将会大大的减小。
附图说明
图1为本发明工作流程示意图。
图2为本发明实施例中接收机网络拓扑结构的示例图。
图3为本发明实施例中发射机的运动轨迹。
图4为本发明实施例与现有集中式粒子滤波算法在信噪比为0的条件下的仿真对比图。
具体实施方式
下面结合附图与实施例对本发明作进一步详细的说明。
本实施例提供基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法,其流程如图1所示,具体包括以下步骤:
3.第一次数据交换:各接收机将自身接收的信号传给邻居接收机,同时接收邻居接收机传过来的相应信息;
4.更新粒子权重:各接收机根据其邻域内接收机上的信息求解局部似然函数,然后更新预测粒子对应的权重;
6.第二次数据交换:各接收机将上一步骤计算得到的局部中间估计值ψl,k传输给邻居接收机,同时接收邻居接收机传过来的相应信息;
对上述方法进行仿真,结果如图4所示,图中将本发明分布式算法与集中式算法进行比较,在相同接收机个数的仿真环境中,集中式的定位效果要比分布式的好,但是与集中式相比,本专利提出的分布式方法的鲁棒性更好,不存在多跳传输,且每个接收机的运算量大大降低。本发明通过引入时差和多普勒构造的分布式定位算法提升稳态性能,从图可以看出,分布式定位的稳态性能基本接近集中式定位稳态性能。
以上所述,仅为本发明的具体实施方式,本说明书中所公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换;所公开的所有特征、或所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以任何方式组合。
Claims (1)
1.一种基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法,其特征在于,包括以下步骤:
步骤1:在接收机个数为N的无线传感器网络当中,局部GPF从局部重要函数中随机采M个随机预测粒子,k时刻具体为:
步骤:2:第一次数据交换,接收机l将接收信号传输给自己的邻居接收机,同时接收邻居接收机传来的接收信号;
步骤3:接收机l更新每个预测粒子所对应的权重值:每个预测粒子所对应的权重值计算为求解如下局部似然函数:
rl,k=bl,kAl,kFl,ksk+ul,k
其中,sk为k时刻的发射信号,ul,k为零均值高斯白噪声,bl,k表示信号衰减,Aj,k和Fj,k分别表示多普勒频移和时延;
其中,C′和C″分别是两个常数,且:
将上式转换为求最大值问题:
最后对粒子权重进行归一化,计算得到每个粒子的归一化权重:
步骤4:接收机l更新它的局部中间状态估计值ψl,k以及相应的协方差Pl,k:
步骤5:第二次数据交换,接收机l将步骤4 计算得到的局部中间状态估计值ψl,k传输给自己的邻居接收机,同时接收邻居接收机传来的结果;
其中,ajl为预设非负加权系数,满足条件:
条件概率密度函数为:xk+1=φxk+nk,φ为转移矩阵,nk为零均值高斯白噪声;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710584073.3A CN107367710B (zh) | 2017-07-18 | 2017-07-18 | 一种基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710584073.3A CN107367710B (zh) | 2017-07-18 | 2017-07-18 | 一种基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107367710A CN107367710A (zh) | 2017-11-21 |
CN107367710B true CN107367710B (zh) | 2020-08-11 |
Family
ID=60306786
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710584073.3A Active CN107367710B (zh) | 2017-07-18 | 2017-07-18 | 一种基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107367710B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108445444B (zh) * | 2018-03-09 | 2020-05-05 | 中国人民解放军战略支援部队信息工程大学 | 修正容积卡尔曼滤波直接定位方法 |
CN109710978B (zh) * | 2018-11-30 | 2022-05-03 | 电子科技大学 | 一种分布式异构自适应粒子滤波直接跟踪方法 |
CN109671100B (zh) * | 2018-11-30 | 2020-09-25 | 电子科技大学 | 一种分布式可变扩散组合系数粒子滤波直接跟踪方法 |
CN110988796B (zh) * | 2019-12-20 | 2022-03-25 | 武汉理工大学 | 一种基于频偏机会信号的定位方法及系统 |
CN111211760B (zh) * | 2020-01-15 | 2023-04-11 | 电子科技大学 | 一种基于分布式扩散策略的反馈粒子滤波方法 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6882959B2 (en) * | 2003-05-02 | 2005-04-19 | Microsoft Corporation | System and process for tracking an object state using a particle filter sensor fusion technique |
CN101222240B (zh) * | 2008-01-22 | 2013-01-16 | 上海师范大学 | 异步ds-cdma盲多用户检测的粒子滤波方法 |
KR100951321B1 (ko) * | 2008-02-27 | 2010-04-05 | 아주대학교산학협력단 | 파티클 필터 기반의 음향 센서를 이용한 3차원 공간에서의객체 추적 방법 |
US9197996B2 (en) * | 2010-11-19 | 2015-11-24 | Commonwealth Scientific And Industrial Research Organisation | Tracking location of mobile device in a wireless network |
CN102256381A (zh) * | 2011-04-07 | 2011-11-23 | 江苏科技大学 | 分布式自适应粒子滤波的无线传感器网络目标跟踪方法 |
EP2831614B1 (en) * | 2012-03-29 | 2019-09-11 | SULU Networks Ltd. | Localization, synchronization and navigation using passive sensor networks |
CN103237348B (zh) * | 2013-05-10 | 2015-09-16 | 重庆大学 | 一种基于wsn改进粒子滤波的运动目标定位方法 |
CN104537257B (zh) * | 2015-01-12 | 2017-09-26 | 电子科技大学 | 一种基于时差的分布式自适应直接定位方法 |
CN104977562B (zh) * | 2015-01-16 | 2018-10-23 | 上海大学 | 完全分布式无线传感网中鲁棒性多声源定位方法 |
CN105704071A (zh) * | 2015-07-07 | 2016-06-22 | 大连大学 | 基于信息序列的自适应渐消扩展卡尔曼粒子滤波的多普勒频移估计方法 |
CN105068550B (zh) * | 2015-08-21 | 2017-10-20 | 燕山大学 | 一种基于拍卖模式的水下机器人多目标选择方法 |
CN105929389A (zh) * | 2015-12-05 | 2016-09-07 | 中国人民解放军信息工程大学 | 一种基于外辐射源时延和多普勒频率的直接定位方法 |
CN106683122B (zh) * | 2016-12-16 | 2020-02-04 | 华南理工大学 | 一种基于高斯混合模型和变分贝叶斯的粒子滤波方法 |
-
2017
- 2017-07-18 CN CN201710584073.3A patent/CN107367710B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107367710A (zh) | 2017-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107367710B (zh) | 一种基于时延和多普勒的分布式自适应粒子滤波直接跟踪定位方法 | |
CN111246491B (zh) | 一种智能反射表面辅助的太赫兹通信系统设计方法 | |
Bai et al. | A novel heavy-tailed mixture distribution based robust Kalman filter for cooperative localization | |
CN109710978B (zh) | 一种分布式异构自适应粒子滤波直接跟踪方法 | |
CN109671100B (zh) | 一种分布式可变扩散组合系数粒子滤波直接跟踪方法 | |
KR20010099929A (ko) | 적응배열안테나시스템의 최적 웨이트벡터 계산을 위한신호처리방법 및 장치 | |
Zhang et al. | An efficient machine learning approach for indoor localization | |
CN107592671B (zh) | 一种网络化多智能体主动变拓扑的自主协同定位方法 | |
Xu et al. | Distributed conjugate gradient strategies for distributed estimation over sensor networks | |
CN109379780B (zh) | 基于自适应差分进化算法的无线传感器网络定位方法 | |
Wang et al. | Two stage particle filter for nonlinear Bayesian estimation | |
Forsling et al. | Consistent distributed track fusion under communication constraints | |
Liu et al. | Channel tracking for RIS-aided mmWave communications under high mobility scenarios | |
CN106792982B (zh) | 一种基于自适应分簇策略的多目标直接定位方法 | |
CN110300380B (zh) | 移动wsn中平衡系统能耗与追踪精度的目标跟踪方法 | |
CN115499916B (zh) | 基于改进鲸鱼优化算法的无线传感器网络节点定位方法 | |
CN115494450B (zh) | 一种高精度的超宽带室内定位跟踪与控制方法及装置 | |
CN110940999A (zh) | 一种基于误差模型下的自适应无迹卡尔曼滤波方法 | |
CN107247250A (zh) | 一种基于粒子滤波的相干分布源波达方向跟踪方法 | |
Miao et al. | Localization accuracy of farmland wireless sensor network localization algorithm based on received signal strength indicator. | |
Chakkravarthy et al. | A machine learning based accurate localization technique for 5G networks | |
Kazemi et al. | Improved weighted RSS positioning algorithm for cognitive radio | |
Sepahvand et al. | Target tracking with unknown maneuvers using adaptive parameter estimation in wireless sensor networks | |
Luo et al. | An Improved Particle-Filter-Based Hybrid Optimization Algorithm for IoT Positioning in Uncertain WSNs | |
Gao et al. | An improved DV-Hop algorithm based on average hop distance and estimated coordinates |
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 |