CN110649911B - 一种基于α散度的分布式非线性卡尔曼滤波方法 - Google Patents
一种基于α散度的分布式非线性卡尔曼滤波方法 Download PDFInfo
- Publication number
- CN110649911B CN110649911B CN201910645751.1A CN201910645751A CN110649911B CN 110649911 B CN110649911 B CN 110649911B CN 201910645751 A CN201910645751 A CN 201910645751A CN 110649911 B CN110649911 B CN 110649911B
- Authority
- CN
- China
- Prior art keywords
- node
- distributed
- divergence
- representing
- distribution function
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000001914 filtration Methods 0.000 title claims abstract description 39
- 239000002245 particle Substances 0.000 claims abstract description 31
- 238000009826 distribution Methods 0.000 claims abstract description 29
- 238000005315 distribution function Methods 0.000 claims abstract description 26
- 238000012545 processing Methods 0.000 claims abstract description 10
- 238000005070 sampling Methods 0.000 claims description 10
- 238000009792 diffusion process Methods 0.000 claims description 9
- 230000007704 transition Effects 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 description 11
- 238000011160 research Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000004927 fusion Effects 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 229940060587 alpha e Drugs 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0248—Filters characterised by a particular frequency response or filtering method
- H03H17/0255—Filters based on statistics
- H03H17/0257—KALMAN filters
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0223—Computation saving measures; Accelerating measures
- H03H17/0238—Measures concerning the arithmetic used
- H03H17/0241—Distributed arithmetic
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Physics (AREA)
- Probability & Statistics with Applications (AREA)
- Computer And Data Communications (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明本发明属于信号处理领域,具体提出一种基于α散度的分布式非线性卡尔曼滤波方法,能有效减小***自身噪声和外部噪声的影响,得到对真实值高精度的估计,可用于解决信号处理中的滤波问题和参数估计问题,尤其是目标跟踪问题。本发明中,首先,通过最小化每个节点的中间近似后验分布函数与真实后验分布函数之间的α散度以获取中间状态估计;然后,通过最小化每个节点的最终后验分布近似函数与其邻居节点的中间后验分布近似函数之间前向KL散度的凸组合来计算每个节点的最终状态估计结果。本发明不同于现有的基于最小方差准则的分布式粒子滤波和基于最小均方误差准则的分布式扩展卡尔曼滤波,本有效提高了非线性分布式卡尔曼滤波的稳态性能。
Description
技术领域
本发明属于信号处理领域,涉及信号处理领域的目标跟踪问题,尤其涉及到分布式无线 传感器网络上目标跟踪问题,具体为一种基于α散度的分布式非线性卡尔曼滤波方法。
背景技术
卡尔曼滤波算法不仅可以对平稳的一维随机过程进行估计,也可以对非平稳的多维随机 过程进行估计,而且卡尔曼滤波算法是递推的,具有存储量小,收敛速度快,实时处理速度快 等优点,因此卡尔曼滤波在复杂***中有更广泛的应用,如导航、目标跟踪、定位等。此外,卡尔曼滤波算法也被用于动态***的预测,如星体运动轨迹,商品交换价格的变化趋势。
目前,对于线性卡尔曼滤波的研究已经取得了巨大的进展,然而对非线性卡尔曼滤波的 研究仍有所欠缺。在卡尔曼滤波中常用后验分布表示其状态估计,由于非线性卡尔曼滤波的 高斯后验分布不具有封闭形式,无法直接求解其后验分布,这就为非线性卡尔曼的研究带来了巨大困难。
随着无线传感器技术的快速发展,近年来,分布式网络成为焦点,大量低成本的无限传 感器被广泛应用于分布式网络中。在分布式网络中,每个节点仅同其邻点进行信息交流,如 此可以大大提高网络的可扩展性和灵活性,能大量节省通信资源,在环境监测,救灾管理,参数估计等领域有着广泛应用。然而目前,分布式的扩散融合技术在分布式线性卡尔曼滤波 中的应用较为成熟,而对于非线性动力***中的研究尚不够充分。现有的分布式非线性卡尔 曼滤波算法有分布式扩展卡尔曼滤波和分布式粒子滤波等。扩展卡尔曼滤波算法对非线性动 力学***做了线性近似处理,再利用线性卡尔曼滤波***的方法进行处理,引入了线性化误差,具有收敛速度慢等局限性。
发明内容
本发明的目的在于提出一种基于α散度的分布式非线性卡尔曼滤波方法,采用一个指数 族的分布函数q(x),来近似真实的后验分布,并采用α散度衡量两个分布之间的差距,α散 度越小表示两个分布越接近;本发明通过直接优化α散度的方法来计算近似分布函数,从而得到对目标真实状态的估计结果。
为实现上述目的,本发明采用的技术方案如下:
一种基于α散度的分布式非线性卡尔曼滤波方法,包括以下步骤:
步骤1.针对节点r,根据其t-1的后验分布函数与状态转移方程计算t 时刻的先验分布函数/>
步骤2:设定t时刻的采样的建议分布函数以及初始化中间近似 分布函数/>
步骤3:节点r根据建议分布进行采样,使得每个粒子都独立同分布于建议分布:
其中,S为输入采样粒子数,表示节点r在t时刻采样的第s个粒子的状态向量;
步骤4:根据下式计算节点k的粒子权重,并进行归一化处理:
其中,表示节点r在t时刻第s个采样粒子的权重,Wr,t表示节点r在t时刻所有采样粒子 的权重之和,yu,t表示节点r的邻居节点u在t时刻的观测值;/>表示邻居节点u对/>观测值的预测,由观测方程计算得到;/>表示节点r对/>的预测,由状态转移方 程计算得到;
步骤5:根据粒子权重计算近似中间后验分布函数其服从高斯分布
步骤6:在分布式网络中,将每个节点得到的近似中间后验分布函数进行扩散, 计算得到节点r的最终状态估计qr(xt):
其中,au,r表示在扩散时邻点u对当前节点r的权重,满足
并计算qr(xt)的均值μr,t和协方差矩阵
μr,t即为节点r的状态估计向量。
本发明的有益效果在于:
本发明提出的一种基于α散度的分布式非线性卡尔曼滤波方法,具有如下优点:
1.本发明提出的方法可用于非线性动力***,相比经典的线性卡尔曼滤波算法,应用范 围更广;
2.本发明提出的是一种分布式扩散策略,每个节点处只需同其邻点交流信息,允许各个 节点同时进行处理,不需要将各个节点的信息送至融合中心进行处理,需要的通信能量少, 运算效率高;
3.本发明提出的分布式算法比相应的集中式算法的稳健性强。对于集中式而言,当处理 中心出现问题时,整个***将无法正常工作,分布式算法能有效避免融合中心出现问题时导 致整个***崩溃的风险;
4.本发明在计算每个节点的中间状态估计时,基于蒙特卡洛技术直接优化的无偏算法直 接优化近似分布与真实后验分布之间的α散度,不需要对非线性函数做线性近似;
5.本发明使用的α散度可以看作是KL散度的推广,其中涉及的参数α满足α∈(-∞,∞)。 当α=1时,α散度可退化为后向KL散度;当α=0时,α散度可退化为前向KL散度;当α=0.5, α散度可退化为海林格距离(Hellinger distance)。本发明可以根据具体的应用场景选择 合适的α的取值,通过调整α的大小控制散度范围,可能获得更优的估计结果(不同α值近 似分布结果如图1所示)。
附图说明
图1为本发明基于α散度的分布式非线性卡尔曼滤波方法中每个节点的流程示意图。
图2为不同α值的近似分布结果。
图3为实施例中采用的分布式网络拓扑结构(以σQ=0.05网络中有30个节点为例)。
图4为实施例中本发明方法与分布式扩展卡尔曼滤波,分布式粒子滤波和集中式α散度 方法在某一次蒙特卡罗实验中一个节点的跟踪结果图。
图5为实施例中过程噪声取σQ=0.05时,本发明提出的基于α散度的分布式非线性卡尔 曼滤波与分布式扩展卡尔曼滤波,分布式粒子滤波和集中式α散度方法的位置均方误差(MS E)的对比图。
图6为实施例中过程噪声取σQ=0.05时,本发明提出的基于α散度的分布式非线性卡尔 曼滤波与分布式扩展卡尔曼滤波,分布式粒子滤波和集中式α散度方法的速度MSE对比图。
表1为实施例中不同的过程噪声参数下的基于α散度的分布式非线性卡尔曼滤波,分布 式扩展卡尔曼滤波,分布式粒子滤波和集中式α散度方法的平均误差结果。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
本发明中,考虑有个节点的分布式网络结构,每个节点r都有一个传感器可对目标进 行跟踪,与节点r直接相连的节点称之为他的邻点(每个节点都与自身相连)表示为u,将r 的邻点网络记为Nr;本发明提出的基于α散度的分布式非线性卡尔曼滤波方法主要基于扩散 策略,可分为两步:在自适应阶段得到每个节点的中间状态估计,再通过组合阶段将中间状 态估计在每个节点的邻域内进行扩散;每个节点根据其所有邻点观测信息累积量得到的的真实后验分布为/>其中,xt表示t时刻的状态向量;/>表示节点k的所有邻居 节点在1~t时刻观测数据的集合:
每个节点的中间状态估计用中间近似后验分布函数来描述,每个节点的最终状态 估计用近似后验分布函数qr(xt)来描述,qr(xt)的均值即为对目标的状态估计值。
用α散度来衡量每个节点的真实后验分布与近似后验分布函数之间的差距,定义每个节 点的α散度表达式为:
其中,α为预设参数,如图2所示,α的取值范围为{-∞,+∞}/{0}。假设各邻节的观测信息 彼此独立,即之间彼此独立。
优化每个节点的α散度,得到在α散度意义下与真实后验分布最接近的中间近似后验分 布函数令/>使α散度最小即对/>与做矩匹配,得到:/>E表示期望;由于此等式两边均有/>不能直接求解,因此本发明通过一次迭代的方法得到/>
本实施例提供一种基于α散度的分布式非线性卡尔曼滤波方法,其流程如图1所示,具 体包括以下步骤:
步骤1.针对节点r,根据其t-1的后验分布函数与状态转移方程计算t 时刻的先验分布函数/>
步骤2:设定t时刻的采样的建议分布函数以及初始化中间近似 分布函数/>
步骤3:节点r根据建议分布进行采样,使得每个粒子都独立同分布于建议分布:
其中,S为输入采样粒子数,表示节点r在t时刻采样的第s个粒子的状态向量;
步骤4:根据下式计算节点k的粒子权重,并进行归一化处理:
其中,表示节点r在t时刻第s个采样粒子的权重,Wr,t表示节点r在t时刻所有采样粒子 的权重之和,yu,t表示节点r的邻居节点u在t时刻的观测值;/>表示邻居节点u对/>观测值的预测,由观测方程计算得到;/>表示节点r对/>的预测,由状态转移方 程计算得到;
步骤5:根据上式求得的粒子权重计算中间近似后验分布其服从高斯分布/>
步骤6:在分布式网络中,将每个节点得到的中间后验分布近似函数进行扩散, 最小化节点r的最终近似后验分布与其所有邻点的中间近似后验分布函数的凸组合之间的前 向KL散度,定义其KL散度为:
其中,au,r表示在扩散时邻点u对当前节点r的贡献,满足
令上述前向KL散度最小,得到:
因为上式中均服从高斯分布,因此qr(xt)也服从高斯分布,其均值和协方差矩阵计算 如下:
得到每个节点的后验分布近似函数qr(xt),其均值即为该节点的状态估计向量;
步骤7:进行时间迭代,得到下一时刻的状态估计结果。
仿真测试
将本发明提出的方法用于分布式网络的目标跟踪中,与分布式扩展卡尔曼滤波,分布式 粒子滤波和集中式α散度方法进行比较。整个分布式网络的拓扑结构共30个节点,拓扑结构 如图3所示,每个节点处一个传感器分别对目标进行跟踪,各节点的观测噪声功率一样。目 标运动的状态空间模型如下:
xt=Ftxt-1+wt,wt~N(0,Qt)
yt=h(xt)+vt,vt~N(0,Rt)
其中,Ft表示状态转移矩阵,h(xt)表示观测函数,wt是状态噪声,vt是观测噪声,wt与vt是均值为零,协方差矩阵分别为Qt和Rt的高斯白噪声:
其中s表示传感器的位置信息。本发明中采取恒定的测量速率,令Δt=1;观测噪声的协方差 矩阵为σR=20;状态向量是四维的:/>包括两个方向的位置信 息,第二维和第四维表示两个方向的速度信息。采样粒子数为1000,迭代次数500次,蒙特卡洛实验100次,α散度参数的取值为α∈{0.5,1},过程噪声的取值为 σQ∈{0.001,0.005,0.01,0.05,0.1,0.3}。
将本发明提出的基于α散度的分布式卡尔曼滤波方法,α分别取0.5和1(图中分别标注 为“分布式AKF,α=0.5”和“分布式AKF,α=1”)与“Nonlinear Kalman FilteringWith Divergence Minimization”所提出的集中式基于α散度卡尔曼滤波方法(图中标注为“集中式 AKF,α=0.5”),“Diffusive Particle Filtering For DistributedMultisensor Estimation”提出的分布式粒子滤波方法(图中标注为“分布式PF”)和“Distributed nonlinear Kalman filtering with applications to wirelesslocalization”提出的分布式扩展卡尔曼滤波(图中标注为“分 布式EKF”)仿真结果如图4、5、6和表1所示。
图4表明本发明提出的基于α散度分布式非线性卡尔曼滤波方法同其他几种方法都能有 效地对目标的位置进行跟踪。图5表明本发明提出的分布式方法与其他方法相比,本发明方 法的位置状态估计MSE曲线的收敛速度与其他方法相差不大,稳态性能略差于集中式的方 法,但比其他两种分布式方法的稳态误差都小,并且在α=0.5时比α=1时(退化为KL散度) 稳态误差更小。
图6表明在过程噪声取值为0.05时,本发明提出的方法与其他三种方法相比,他们的速 度状态估计MSE曲线的收敛速度大致都相同,但本发明提出的方法的稳态误差比其他方法都小,性能更好,甚至比集中式的方法都好,并且在α=0.5时效果更好。表1表明在过程噪声 较小时,取值在一定范围内,本发明提出的方法相比其他两种分布式方法的稳态误差更好, 但会略小于集中式的方法,α=0.5比α=1的平均误差更小。
表1
以上所述,仅为本发明的具体实施方式,本说明书中所公开的任一特征,除非特别叙述, 均可被其他等效或具有类似目的的替代特征加以替换;所公开的所有特征、或所有方法或过 程中的步骤,除了互相排斥的特征和/或步骤以外,均可以任何方式组合。
Claims (1)
1.一种基于α散度的分布式非线性卡尔曼滤波方法,该方法用于分布式网络的目标跟踪中,设置分布式网络有个节点,每个节点r都有一个传感器对目标进行跟踪;该方法包括以下步骤:
步骤1.针对节点r,根据其t-1时刻的后验分布函数与状态转移方程计算t时刻的先验分布函数/>
步骤2:设定t时刻的采样的建议分布函数以及初始化中间近似分布函数/>
步骤3:节点r根据建议分布进行采样,使得每个粒子都独立同分布于建议分布:
其中,S为输入采样粒子数,表示节点r在t时刻采样的第s个粒子的状态向量;
步骤4:根据下式计算节点r的粒子权重,并进行归一化处理:
其中,表示节点r在t时刻第s个采样粒子的权重,Wr,t表示节点r在t时刻所有采样粒子的权重之和,yu,t表示节点r的邻居节点u在t时刻的观测值;/>表示邻居节点u对观测值的预测,由观测方程计算得到;/>表示节点r对/>的预测,由状态转移方程计算得到;
步骤5:根据粒子权重计算近似中间后验分布函数其服从高斯分布
步骤6:在分布式网络中,将每个节点得到的近似中间后验分布函数进行扩散,计算得到节点r的最终状态估计qr(xt):
其中,au,r表示在扩散时邻点u对当前节点r的权重,满足
并计算qr(xt)的均值μr,t和协方差矩阵
μr,t即为节点r的状态估计向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910645751.1A CN110649911B (zh) | 2019-07-17 | 2019-07-17 | 一种基于α散度的分布式非线性卡尔曼滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910645751.1A CN110649911B (zh) | 2019-07-17 | 2019-07-17 | 一种基于α散度的分布式非线性卡尔曼滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110649911A CN110649911A (zh) | 2020-01-03 |
CN110649911B true CN110649911B (zh) | 2023-10-27 |
Family
ID=68989877
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910645751.1A Active CN110649911B (zh) | 2019-07-17 | 2019-07-17 | 一种基于α散度的分布式非线性卡尔曼滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110649911B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111211760B (zh) * | 2020-01-15 | 2023-04-11 | 电子科技大学 | 一种基于分布式扩散策略的反馈粒子滤波方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102082560A (zh) * | 2011-02-28 | 2011-06-01 | 哈尔滨工程大学 | 一种基于集合卡尔曼滤波的粒子滤波方法 |
CN102779333A (zh) * | 2012-07-10 | 2012-11-14 | 武汉大学 | 一种基于卡尔曼滤波的光学影像复原方法 |
CN105704071A (zh) * | 2015-07-07 | 2016-06-22 | 大连大学 | 基于信息序列的自适应渐消扩展卡尔曼粒子滤波的多普勒频移估计方法 |
WO2016155241A1 (zh) * | 2015-04-03 | 2016-10-06 | 百度在线网络技术(北京)有限公司 | 基于Kalman滤波器的容量预测方法、***和计算机设备 |
CN108566178A (zh) * | 2018-04-16 | 2018-09-21 | 武汉理工大学 | 一种非稳态随机机会网络特征值滤波方法 |
CN108599737A (zh) * | 2018-04-10 | 2018-09-28 | 西北工业大学 | 一种变分贝叶斯的非线性卡尔曼滤波器的设计方法 |
EP3435101A1 (de) * | 2017-07-26 | 2019-01-30 | INVENOX GmbH | Verfahren und vorrichtung zur überwachung eines stabilen konvergenzverhaltens eines kalman-filters |
CN109508445A (zh) * | 2019-01-14 | 2019-03-22 | 哈尔滨工程大学 | 一种带有色量测噪声和变分贝叶斯自适应卡尔曼滤波的目标跟踪方法 |
-
2019
- 2019-07-17 CN CN201910645751.1A patent/CN110649911B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102082560A (zh) * | 2011-02-28 | 2011-06-01 | 哈尔滨工程大学 | 一种基于集合卡尔曼滤波的粒子滤波方法 |
CN102779333A (zh) * | 2012-07-10 | 2012-11-14 | 武汉大学 | 一种基于卡尔曼滤波的光学影像复原方法 |
WO2016155241A1 (zh) * | 2015-04-03 | 2016-10-06 | 百度在线网络技术(北京)有限公司 | 基于Kalman滤波器的容量预测方法、***和计算机设备 |
CN105704071A (zh) * | 2015-07-07 | 2016-06-22 | 大连大学 | 基于信息序列的自适应渐消扩展卡尔曼粒子滤波的多普勒频移估计方法 |
EP3435101A1 (de) * | 2017-07-26 | 2019-01-30 | INVENOX GmbH | Verfahren und vorrichtung zur überwachung eines stabilen konvergenzverhaltens eines kalman-filters |
CN108599737A (zh) * | 2018-04-10 | 2018-09-28 | 西北工业大学 | 一种变分贝叶斯的非线性卡尔曼滤波器的设计方法 |
CN108566178A (zh) * | 2018-04-16 | 2018-09-21 | 武汉理工大学 | 一种非稳态随机机会网络特征值滤波方法 |
CN109508445A (zh) * | 2019-01-14 | 2019-03-22 | 哈尔滨工程大学 | 一种带有色量测噪声和变分贝叶斯自适应卡尔曼滤波的目标跟踪方法 |
Non-Patent Citations (2)
Title |
---|
State Estimation Using Non-uniform and Delayed Information: A Review;Jhon A.Isaza;Hector A.Botero;Hernan Alvarez;;International Journal of Automation and Computing(第02期);315-320 * |
带扩展卡尔曼滤波的柔性关节机器人虚拟分解控制;刘华山;李生;王保相;;控制理论与应用(第06期);212-217 * |
Also Published As
Publication number | Publication date |
---|---|
CN110649911A (zh) | 2020-01-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109508445B (zh) | 一种带有色量测噪声和变分贝叶斯自适应卡尔曼滤波的目标跟踪方法 | |
CN110244715B (zh) | 一种基于超宽带技术的多移动机器人高精度协同跟踪方法 | |
CN108318856B (zh) | 一种异构网络下快速精准的目标定位与跟踪方法 | |
Yan et al. | An improved NLOS identification and mitigation approach for target tracking in wireless sensor networks | |
CN110673089B (zh) | 未知视距和非视距分布情况下基于到达时间的定位方法 | |
CN108134640A (zh) | 一种基于节点运动状态约束的协作定位***和方法 | |
Xie et al. | An improved algorithm based on particle filter for 3D UAV target tracking | |
CN111711432B (zh) | 一种基于ukf和pf混合滤波的目标跟踪算法 | |
CN111798494A (zh) | 广义相关熵准则下的机动目标鲁棒跟踪方法 | |
CN110649911B (zh) | 一种基于α散度的分布式非线性卡尔曼滤波方法 | |
CN110516198B (zh) | 一种分布式非线性卡尔曼滤波方法 | |
CN112034713B (zh) | 非理想网络环境下移动目标的最优状态估计方法及*** | |
CN112904274A (zh) | 一种改进tdoa/fdoa算法的多运动目标定位方法 | |
Xu et al. | Improved Indoor Localization Based on Received Signal Strength Indicator and General Regression Neural Network. | |
CN115685128A (zh) | 一种机动目标场景下的雷达目标跟踪算法及电子设备 | |
CN112800889B (zh) | 一种基于分布式矩阵加权融合高斯滤波的目标跟踪方法 | |
CN115114985A (zh) | 一种基于集合理论的传感器***分布式融合方法 | |
CN105654517A (zh) | 基于分层空间的rb化粒子滤波算法 | |
CN112596383A (zh) | 一种提升低成本mems阵列式传感器环境鲁棒性的控制方法 | |
CN111695617A (zh) | 一种基于改进的协方差交叉算法的分布式火控融合方法 | |
Sepahvand et al. | Target tracking with unknown maneuvers using adaptive parameter estimation in wireless sensor networks | |
CN112491393B (zh) | 一种基于观测噪声协方差矩阵未知的线性滤波方法 | |
Bhavesh et al. | Smart Device Localization under Generic α-κ-µ Shadowed Fading IoT Environment using Signed Gradient Method | |
CN112748416B (zh) | 一阶传播的多节点分散式gm-phd融合方法 | |
Bian et al. | An indoor localization algorithm based on IMM in mixed LOS/NLOS environments |
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 |