CN101438184B - 一种跟踪移动电子设备的状态的方法 - Google Patents
一种跟踪移动电子设备的状态的方法 Download PDFInfo
- Publication number
- CN101438184B CN101438184B CN200680054493.0A CN200680054493A CN101438184B CN 101438184 B CN101438184 B CN 101438184B CN 200680054493 A CN200680054493 A CN 200680054493A CN 101438184 B CN101438184 B CN 101438184B
- Authority
- CN
- China
- Prior art keywords
- grid
- community
- state
- mobile electronic
- probable value
- 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
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/0294—Trajectory determination or predictive filtering, e.g. target tracking or Kalman filtering
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/029—Location-based management or tracking services
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Mobile Radio Communication Systems (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及跟踪移动电子设备的状态的方法以及包括被配置来执行所述方法的处理装置的移动电子设备。一种跟踪移动电子设备状态的方法,所述方法包括重复执行步骤(i)使用包括多个小区的网格表示移动电子设备的状态,每个小区表示由一个或多个状态变量定义的状态空间中的区域,并具有所述移动电子设备的状态在状态空间中的该区域内的概率值,所述网格被限制为仅包括具有高于预先确定的阈值的概率值的小区;(ii)获取指示一个或多个状态变量的值的测量信号;(iii)基于所述测量信号更新所述网格的所述概率值并重新限制所述网格。
Description
技术领域
本发明涉及跟踪移动电子设备状态的方法,以及包括被配置为执行该方法的处理装置的移动电子设备。
背景技术
个人定位的挑战之一是在仅有很少的并且可能具有带有异常分布状态的较大误差的可用测量资源的环境中,尤其是在需要由这些不同的测量源汇合形成非线性过滤方程的有效的数值解决方案的室内或城市区域,提供准确的位置信息。
诸如GPS的基于卫星的系统的状态最多是在室内以高灵敏度模式使用时是不可预测的。局域无线网络,例如蜂窝网络、WLAN或蓝牙提供了一定的定位能力,但是比起GPS具有较低的准确度。移动电子设备的其他可能的部件是单板(on-board)传感器,例如过载传感器(accelerometer)、气压计或数字罗盘(compasses)。
由于不同的误差特征、不可预测的失真、测量中的系统误差、强的非线性、复杂的时间相依性以及丢失的数据,因而合并各种测量源是困难的。以通常的方式对所有的情况进行建模并非易事,更不用说精确地解释这些模型。即使拥有正确的模型,常用的Kalman滤波器和其非线性扩展也可能没有警告地失效。
发明内容
根据本发明的第一方面,提供了跟踪移动电子设备的状态的方法,所述方法包括重复执行步骤:
(i)使用包括多个小区的网格表示所述移动电子设备的状态,每个小区表示由一个或多个状态变量定义的状态空间中的区域,并具有所述移动电子设备的状态在状态空间中的该区域内的概率值,所述网格被限制为仅包括具有高于预先确定的阈值的概率值的小区;
(ii)获得指示一个或多个状态变量的值的测量信号;
(iii)基于所述测量信号更新所述网格的所述概率值并重新限制所述网格。
不论是单独还是结合,所述移动电子设备的状态可以包括但不限于位置、速度、加速度以及时钟误差状态变量。
所述状态空间的维数对应于状态变量的总维数,每个小区包括与该状态空间的维数相等的维数。例如,状态空间和小区可以是六维的,以三维表示位置,以三维表示速度。
本发明提供了手段由此从不产生唯一的定位解决方案或产生多个解决方案的测量几何学保留所有信息,并由此在后分布的形状被保留而不漏掉其任何峰值。
本发明比粒子滤波器运行的快,比位置指纹方法更普遍,传播完整的位置分布从而能够精确地表示不确定的或多个解决方案的系统。
优选地,所述网格是均匀分隔的平行六面体网格,以缓解关于所述方法的计算负载的问题。
步骤(i)就一个或多个重复而言可以包括获取测量信号以及基于所述测量信号形成网格。位置估计可以从包括任何例如距离、距离差别或平面测量的测量信号获得。距离测量可以例如以时间延迟、往返行程(round-trip)或信号强度测量的形式,从蜂窝基站、WLAN或蓝牙发射机和/或声学传感器获得。指示位置的其他类型的测量信号包括但不限于到达角度的测量、最大-最小距离窗口、基本扇区(basesector)信息、以及单板气压计和数字罗盘。可以使用伪距(deltarange)和航向(heading)测量来测量速度。可以使用单板过载传感器来测量加速度。
应该理解测量信号的细节对于本发明并不重要,假设可以使用测量信号,不论单独使用还是与其他测量信号结合,来获得或估计所提及的状态变量的值。
步骤(i)就一个或多个第二个和后续重复而言,可以包括使用在先重复的步骤(iii)的所述更新网格。
所述一个或多个重复可以包括步骤(iv),基于所述更新网格和用于所述移动电子设备的运动模型,预测所述后续重复的所述网格的所述概率值以获得预测的网格。
步骤(iv)可以包括重新限制所述网格和/或在所述后续重复的步骤(i)中使用所述预测的网格。
重新限制所述网格可以包括移动所述网格的边界以排除概率值低于所述预先确定的阈值的小区,和/或移动边界以包括概率值高于所述预先确定的阈值的小区。在一些情况中,小区存在于所述边界之外但具有的概率值为零。如果,遵循任何处理,这样的小区的概率值增长了,则边界被移动以包括任何这种具有高于预先确定的阈值的概率值的小区。在一些情况中,在所述边界之外不存在小区。随着任何导致概率值的改变的处理,本发明可以包括定义边界之外的临时小区,为这些小区计算概率值,以及移动边界以包括任何这种具有高于预先确定的阈值的概率值的小区的步骤。
优选地,所述运动模型是线性的,以便缓解关于所述方法的计算负载的问题。
一个或多个重复可以包括步骤(v)基于所述更新网格计算用于所述移动电子设备的所述状态的期望的值和变化。
根据本发明的第二方面,提供了包括被配置为执行本发明的第一方面的方法的处理装置的移动电子设备。
本发明还包括计算机程序,被配置为执行所述方法以及其中运行了本发明的移动电子设备的系统。本发明涵盖了一个或多个方面以及各种组合的实施例,无论在该结合中是否明确地提及(或请求保护)。
附图说明
为了可以更容易地理解本发明,现在参考附图仅以实例的方式给出说明,其中:
图1至图5示例性说明本发明实现的方法的简单实例;
图6示出了根据本发明的设备。
具体实施方式
在下文的描述中,相同的附图标记代表相同的特征,不论这些特征属于哪个实施例。
图6是根据本发明的移动电子设备的局部概略图,移动电子设备包括与存储器102和多个收发器模块104a-d通信的处理器100。如下面将要描述的,存储器102包括软件代码部分用于执行本发明的方法。收发器模块104a-d分别能够接收来自GPS系统、蜂窝基站、WLAN发射机和蓝牙发射机的测量信号。为了方便收发器模块104a-d被如此称呼;事实上,一些收发器模块104a-d可能没有传输能力,而其他的可能被要求传输以便接收测量信号。
现在将使用下列名称描述本发明的方法。
下标k标记瞬时时间
k|k-1指在时间步骤k使用仅到时间k-1的数据做出的推论
k|k指在时间步骤k使用到时间k的数据做出的推论
x移动电子设备的状态
i小区的向量下标(index)
Gk(i)第i个小区
ck(i)第i个小区的中心
Ek方阵(squarematrix),其列定义小区的边缘
πk|k-1(i)第i个小区的在先概率值
πk|k(i)第i个小区的在后的概率值
φk(x|z)假设先前状态为z,状态x的概率密度
Lk(x)测量似然函数
τk(i-j)从第j个小区移动到第i个小区的概率(基于运动模型)
步骤(i)
在该步骤中,使用包括多个小区(cell)的网格(grid)表示移动电子设备的状态xk。每个小区表示状态空间中的区域,并且具有移动电子设备的状态xk在该区域内的概率值。
例如,状态xk可以包括三维状态变量位置rk和速度vk, 。那么每个小区是六维的用以表示状态空间的六维区域。然而,为了示例性目的,图1至5示出简单实例,其中状态xk包括二维位置rk。
图1示出在先(prior)分布10,对于移动电子设备的位置具有95%、66%和50%的置信(confidence)区域。可以从几个源之一得到在先分布10:(i)本方法在先迭代的在后分布;(ii)通过使用运动模型改变在后分布来获得预测的在先分布,所述运动模型用以预测移动电子设备的运动;以及(iii)通过任何收发器模块104a-d获得的测量信号。在上述源中的几个为可用的情况下,可以使用合并。
使用在先网格12来模拟(approximate)在先分布10,所述网格12由许多统一尺寸和形状的二维小区14组成。每个小区14代表地球表面上区域,并具有移动电子设备位于该区域内的在先概率值。
应当理解图1仅为示例性目的示出以网格12上的置信区域表示的在先分布10。尽管在图1中示出的在先分布10是连续的,但是它将仅由在先网格12的概率值以模拟的形式存在于移动电子设备的存储器102中。图1中小区14的相关阴影表示概率值,即,较深阴影的小区14表明对于小区14的较高的概率值。
在先网格12包括边界16,其中所有的小区14都具有高于预先确定的阈值的在先概率值。因此,在先分布的网格模拟被边界16截成表示重要的范围S,其是(简单连接)Rd中的区域,其中在先概率值是不可忽视的。
在图1的实例中,设置了阈值使得99%的在先分布由在先网格12表示。然而,应当理解阈值的级别是设计选择的问题。较低的阈值将导致更多的小区14被包括在边界16内,因此虽然以增加的计算负载为代价但是获得了更精确的模拟。相反,较高的阈值将导致较少的小区14被包括在边界内,因此以较少的计算负载获得了不太精确的模拟。
只有边界16内的小区14被定义在移动电子设备的存储器102中。在变型中,存储器102定义边界线外的小区,但将这些小区的在先概率值设置为零。在任一情况中,在先网格12对在先分布10的模拟由边界16截短以有助于计算。
在先方格12中的小区14的数量是设计选择问题,以便在计算负载和精确度两者之间找到平衡。一种极端的选择是生成大量的小小区。在这种情况下,如同已知的仅使用与根据本发明的小区14中心的密度值相当的密度值的点群(point-mass)滤波器的情况一样,即使概率值不是最理想的,模拟也是渐近精确的。另一个极端的选择是使用少数量的大小区14。其优点是在先概率值被尽可能精确地计算。当用大小区14模拟时,大部分在先分布10的结构将丢失。最优地,小区14不应当比在先分布10的最细微的(finest)特征小太多。
在该步骤中,时间下标k被设置为k=1。
步骤(ii)
在该步骤中,通过指示一个或多个状态变量值的收发器模块104a-d获得测量信号。
图2示出了基于通过收发器模块104a-d获得的测量信号表示移动电子设备的估计位置的测量似然函数18。如图所示,测量似然函数包括95%、66%和50%置信区域。
图3示出使用在先网格12以与在先分布10相似的方式模拟测量似然函数18。在图3中,概率值再次表示仅基于测量似然函数18,移动电子设备位于由各个小区14表示的区域内的概率。
尽管在图3中,使用与用于在先分布10相同的在先网格12模拟测量似然函数18,应当理解测量似然函数可以使用不同尺寸、形状和/或方向的网格来表示,但是这将需要更多的计算。
测量信号可以包括例如任何距离(range)、距离差或平面的测量,以便获得位置估计。
假设真实位置r,在位置s处的对站点的距离测量可以写为h(r)=||s-r||。相关的测量误差v不需要是标准的,并且由与实际情形相匹配的经验确定的分布表示。
从GPS系统获得的偏差的距离测量被视为距离差。一个站点被选为参考站点,所有差别都相对于它产生。如果参考站点位于s0,距离差测量为h(r)=||s-r||-||s0-r||。
最后,平面测量为h(r)=uTr,其中u为单位向量。
在该步骤中,所有可用的测量都被迭加(stack)到向量yk中并且相应的测量方程被迭加到向量函数hk(x)。
作为实例考虑具有nd距离差测量,nr距离测量,以及np平面测量的情况。则测量向量为y=[d1...dndr1...rnrα1...αnp]T,测量模型为
如果所有测量误差都具有标准的分布,则测量似然函数18为 其中∑为测量误差v的协方差矩阵。
步骤(iii)
在该步骤中,基于测量信号更新在先网格12的在先概率值以生成在后网格12′,在后网格12′被重新限界。
图4示出从在先分布10和测量似然函数18得到的在后分布20。在后分布20由具有小区14′的在后网格12′来表示,每个小区14′具有在后概率值以模拟在后分布20。
通过将每个小区14中的在先概率值乘以在小区14上对似然函数积分建立的小区14中的总的似然值,建立在后网格12′的在后概率值:
在后网格12′被用新的边界16′来重新限界,以便仅包括具有高于预先确定的阈值的概率值的小区14′。如果先前边界16内的小区14具有低于预先确定的阈值的在后概率值,则安置新的边界16′以排除该小区14。
步骤(iv)
在该可选的步骤中,基于当前迭代的在后网格12′和用于移动电子设备的动作模型,得到具有小区14″的预测网格12″。预测网格12″形成后续迭代的在先网格12。
图5示出遵循运动模型的应用的预测在先分布10′,该运动模型可以是任何适合的模型,例如恒定速度法向分布加速度(constant-velocity-normally-distributed-acceleration)模型或singer模型。预测的在先分布10′将构成本方法的后续迭代的在先分布10。
根据下列方程,通过对由转换到该小区14″的概率加权的当前迭代中的所有概率值来产生小区14″的预测在先概率值。在计算了所有的概率值之后预测在先分布10′被标准化。
用αk=|detEk|表示每个小区14″的容量(volume),通过在小区14″上对推测性的pdf积分来计算预测的在先概率值πk|k-1(i):
用其网格模拟当量(yeild)替代Pk-1|k-1
其中Γk(i|j)为从第(k-1)个网格的第j个小区到第k个网格的第i个小区的转换概率(transitionprobability)。
为了计算效率,使用诸如f(x)≡Tx的线性运动模型。通过将运动模型应用到在后网格12′形成预测网格12″。旧的网格中的第i个小区和新的网格的第j个小区两者之间的转换概率仅取决于差i-i,我们可以写为Γk(i|i)=τk(i-j)。预测的在先概率值成为 ,其作为d维离散线性卷积计算的较快。
现在转换概率τk(i-j)为
其中
其可简化为
其中是修改的处理噪声概率密度函数的简略的表达方式。具体地,如果,wk~N(O,Q),则
wk,i-j~N(j-i-Ek -1(ek-Tek-1),(Ek -1)TQEk -1)
那么积分仅是超盒(hyper-box)中的多重正态(multinormal)概率,并且能够在数值上计算。
如果wk是非高斯的,则可以使用累积分布计算转换概率Tk(i-j)。
遵照运动模型的应用,边界16′被移动,以形成边界16″,其包括先前边界16′内部和外部所有具有高于预先确定的阈值的预测概率值的小区14″。
在小区14″位于先前边界16′之外并且具有为零的在前和在后概率值的情况下,新的边界16″被放置以包括任何那些具有高于预先确定的阈值的预测概率值的小区14″。
在其中没有小区位于先前边界16′外的变型中,添加具有与在先小区14′对应的尺寸、形状和方位的新小区14″,放置新的边界16″以包括新小区14″。应当理解新小区14″的尺寸、形状和/或方位不需要与那些在先预测网格12″一致。通过定义先前边界16′外的若干临时小区14″并为临时小区14″计算预测的概率值,可以容易地确定这些不存在的小区是否将具有高于预先确定的阈值的概率值。这样的临时小区14″被不断地定义并且他们的概率值被计算,直到达到概率值低于阈值的多个临时小区14″的外形(shell)被创建的点。随后新的边界16″被放置以便包括所有具有高于阈值的概率值的临时小区14″。
步骤(v)
在该步骤中,根据下列方程,基于预测的网格12″计算用于移动电子设备状态的期期望的值和变化。
(10)
最后,增加时间下标k并从步骤(i)重复该方法。
可以理解的是本发明使用了网格块(grid-mass)的方式来执行跟踪移动电子设备的状态的方法。
Claims (15)
1.一种跟踪移动电子设备的状态的方法,所述方法包括重复执行步骤:
(i)使用包括多个小区的网格表示所述移动电子设备的状态,每个小区表示由一个或多个状态变量定义的状态空间中的区域,并具有所述移动电子设备的状态在状态空间中的该区域内的第一概率值,所述网格被限制为仅包括具有高于预先确定的阈值的概率值的小区;
(ii)获取指示一个或多个状态变量的值的测量信号;
(iii)更新所述网格的所述第一概率值并重新限制所述网格以获得包括具有第二概率值的小区的更新的网格,所述更新的网格的小区的所述第二概率值通过用测量似然函数在所述网格的对应小区上的积分乘以所述第一概率值而建立;并且
(iv)基于所述更新的网格和用于所述移动电子设备的运动模型,预测所述网格的第三概率值,以获得用于后续重复步骤(i)中的预测的网格,所述预测的网格的小区的所述第三概率值通过对由转换到所述预测的网格所述小区的概率加权的当前重复中的所有第二概率值求和而建立。
2.根据权利要求1的方法,其中所述网格是均匀分隔的平行六面体网格。
3.根据权利要求1的方法,其中步骤(i)就一个或多个重复而言包括获取测量信号以及基于所述测量信号形成所述网格。
4.根据权利要求1的方法,其中步骤(iv)包括重新限制所述网格。
5.根据权利要求1的方法,其中所述运动模型是线性的。
6.根据权利要求1的方法,其中所述方法就一个或多个重复而言包括步骤:
(v)基于所述预测网格计算用于所述移动电子设备的所述状态的期望的值和变化。
7.一种被配置为跟踪移动电子设备的状态的设备,所述设备包括:
装置,用于使用包括多个小区的网格表示移动电子设备的所述状态,每个小区表示由一个或多个状态变量定义的状态空间中的区域,并具有所述移动电子设备的状态在状态空间中的该区域内的第一概率值,所述网格被限制为仅包括具有高于预先确定的阈值的概率值的小区;
用于获取指示一个或多个状态变量的值的测量信号的装置;
用于更新所述网格的所述第一概率值并重新限制所述网格以获得包括具有第二概率值的小区的更新的网格的装置,所述更新的网格的小区的所述第二概率值通过用测量似然函数在所述网格的对应小区上的积分乘以所述第一概率值而建立;并且
基于所述更新的网格和用于所述移动电子设备的运动模型,预测所述网格的第三概率值,以获得用于后续重复步骤(i)中的预测的网格,所述预测的网格的小区的所述第三概率值通过对由转换到所述预测的网格所述小区的概率加权的当前重复中的所有第二概率值求和而建立。
8.根据权利要求7的设备,其中所述设备被配置为使用是均匀分隔的平行六面体网格的网格表示移动电子设备的所述状态。
9.根据权利要求7的设备,其中所述设备被配置为在一个或多个重复中通过获取测量信号以及基于所述测量信号形成所述网格表示移动电子设备的所述状态。
10.根据权利要求7的设备,其中所述设备被配置为在预测所述第三概率值后重新限制所述网格。
11.根据权利要求7的设备,其中所述设备被配置为使用线性运动模型表示移动电子设备的所述状态。
12.根据权利要求7的设备,其中所述设备被配置为,在一个或多个重复中,基于所述预测的网格计算用于所述移动电子设备的所述状态的期望的值和变化。
13.根据权利要求7的设备,其中所述设备被包含在移动电子设备中。
14.根据权利要求7的设备,其中所述设备被包含在系统中。
15.根据权利要求7的设备,所述设备包括:多个收发器模块,所述多个收发器模块被配置为获得所述测量信号。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/EP2006/002272 WO2007101453A1 (en) | 2006-03-07 | 2006-03-07 | A method of tracking a state of a mobile electronic device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101438184A CN101438184A (zh) | 2009-05-20 |
CN101438184B true CN101438184B (zh) | 2016-02-17 |
Family
ID=37101900
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200680054493.0A Active CN101438184B (zh) | 2006-03-07 | 2006-03-07 | 一种跟踪移动电子设备的状态的方法 |
Country Status (7)
Country | Link |
---|---|
US (1) | US8213956B2 (zh) |
EP (1) | EP1991882B1 (zh) |
JP (1) | JP5039063B2 (zh) |
KR (1) | KR101106276B1 (zh) |
CN (1) | CN101438184B (zh) |
BR (1) | BRPI0621604B1 (zh) |
WO (1) | WO2007101453A1 (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5094589B2 (ja) * | 2008-06-25 | 2012-12-12 | ヤフー株式会社 | 現在位置推定装置、方法及びシステム |
FR2944603B1 (fr) * | 2009-04-17 | 2014-07-11 | Univ Troyes Technologie | Systeme et procede de localisation de cible par un reseau d'emetteurs recepteurs |
US20110274094A1 (en) * | 2010-05-06 | 2011-11-10 | Qualcomm Incorporated | Methods and apparatus for estimating a geolocation of a wireless communications device |
US9062978B2 (en) * | 2011-05-31 | 2015-06-23 | Massachusetts Institute Of Technology | Tracking a body by nonlinear and non-Gaussian parametric filtering |
FR2979900B1 (fr) | 2011-09-12 | 2013-08-30 | Eurocopter France | Aeronef rapide a grande distance franchissable |
CN103916954B (zh) * | 2013-01-07 | 2017-11-03 | 华为技术有限公司 | 基于wlan的概率定位方法及定位装置 |
US20150181372A1 (en) * | 2013-12-20 | 2015-06-25 | Apple Inc. | Location Accuracy Prediction |
US10488505B2 (en) * | 2014-05-30 | 2019-11-26 | The Boeing Company | Positioning in indoor locations and other GPS-denied environments |
US9929817B2 (en) * | 2014-11-19 | 2018-03-27 | Oath Inc. | System and method for 3D tracking for ad-hoc cross-device interaction |
CN105699937B (zh) * | 2014-11-27 | 2017-08-25 | 中国科学院沈阳自动化研究所 | 一种基于css定位技术的变电站高精度混合定位方法 |
US10444269B2 (en) | 2017-05-26 | 2019-10-15 | Honeywell International Inc. | Apparatus and method for performing grid adaption in numerical solution of recursive bayesian estimators |
US11378403B2 (en) | 2019-07-26 | 2022-07-05 | Honeywell International Inc. | Apparatus and method for terrain aided navigation using inertial position |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1722897A (zh) * | 2004-07-13 | 2006-01-18 | 中兴通讯股份有限公司 | 一种移动通讯中移动终端的定位方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4796191A (en) | 1984-06-07 | 1989-01-03 | Etak, Inc. | Vehicle navigational system and method |
EP0514887B1 (en) | 1991-05-21 | 1997-04-16 | Matsushita Electric Industrial Co., Ltd. | Vehicle position detecting apparatus |
US6889053B1 (en) * | 1999-07-26 | 2005-05-03 | Lucent Technologies Inc. | Likelihood-based geolocation prediction algorithms for CDMA systems using pilot strength measurements |
JP3943339B2 (ja) * | 2000-02-25 | 2007-07-11 | 株式会社エヌ・ティ・ティ・ドコモ | 移動通信システムにおける移動機の位置推定方法及びシステム |
SE0004096D0 (sv) * | 2000-11-08 | 2000-11-08 | Nira Automotive Ab | Positioning system |
FI111901B (fi) * | 2000-12-29 | 2003-09-30 | Ekahau Oy | Sijainnin arviointi langattomissa tietoliikenneverkoissa |
US6542116B1 (en) * | 2001-06-22 | 2003-04-01 | Enuvis, Inc. | Determining the spatio-temporal and kinematic parameters of a signal receiver and its clock by information fusion |
FI113092B (fi) * | 2002-05-31 | 2004-02-27 | Ekahau Oy | Paikannusepävarmuuden mittauksia ja niiden sovelluksia |
-
2006
- 2006-03-07 CN CN200680054493.0A patent/CN101438184B/zh active Active
- 2006-03-07 WO PCT/EP2006/002272 patent/WO2007101453A1/en active Application Filing
- 2006-03-07 BR BRPI0621604-8A patent/BRPI0621604B1/pt not_active IP Right Cessation
- 2006-03-07 JP JP2008557597A patent/JP5039063B2/ja active Active
- 2006-03-07 EP EP06723379.1A patent/EP1991882B1/en active Active
- 2006-03-07 US US12/224,767 patent/US8213956B2/en active Active
- 2006-03-07 KR KR1020087024373A patent/KR101106276B1/ko active IP Right Grant
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1722897A (zh) * | 2004-07-13 | 2006-01-18 | 中兴通讯股份有限公司 | 一种移动通讯中移动终端的定位方法 |
Also Published As
Publication number | Publication date |
---|---|
BRPI0621604B1 (pt) | 2019-02-19 |
US20090233619A1 (en) | 2009-09-17 |
US8213956B2 (en) | 2012-07-03 |
EP1991882B1 (en) | 2016-11-02 |
BRPI0621604A2 (pt) | 2013-03-12 |
KR20080108521A (ko) | 2008-12-15 |
CN101438184A (zh) | 2009-05-20 |
WO2007101453A1 (en) | 2007-09-13 |
JP2009529128A (ja) | 2009-08-13 |
KR101106276B1 (ko) | 2012-01-18 |
EP1991882A1 (en) | 2008-11-19 |
JP5039063B2 (ja) | 2012-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101438184B (zh) | 一种跟踪移动电子设备的状态的方法 | |
St-Pierre et al. | Comparison between the unscented Kalman filter and the extended Kalman filter for the position estimation module of an integrated navigation information system | |
Schon et al. | Marginalized particle filters for mixed linear/nonlinear state-space models | |
Gustafsson | Particle filter theory and practice with positioning applications | |
Fernández-Prades et al. | Bayesian nonlinear filtering using quadrature and cubature rules applied to sensor data fusion for positioning | |
CN105424036A (zh) | 一种低成本水下潜器地形辅助惯性组合导航定位方法 | |
CN104408744A (zh) | 一种用于目标跟踪的强跟踪容积卡尔曼滤波方法 | |
CN111190211B (zh) | 一种gps失效位置预测定位方法 | |
CN105203110A (zh) | 一种基于大气阻力模型补偿的低轨卫星轨道预报方法 | |
CN103776416A (zh) | 一种基于应变传感器的天线结构变形的间接测量方法 | |
CN103152826A (zh) | 一种基于nlos状态检测补偿的移动目标追踪方法 | |
CN104266650A (zh) | 一种基于采样点继承策略的火星着陆器大气进入段导航方法 | |
CN116449374B (zh) | 一种水下声呐差分定位方法 | |
CN108562290A (zh) | 导航数据的滤波方法、装置、计算机设备及存储介质 | |
CN108508463B (zh) | 基于Fourier-Hermite正交多项式扩展椭球集员滤波方法 | |
Du et al. | A hybrid fusion strategy for the land vehicle navigation using MEMS INS, odometer and GNSS | |
CN106597428A (zh) | 一种海面目标航向航速估算方法 | |
CN105180955A (zh) | 机动车实时精准定位方法及装置 | |
CN103616024B (zh) | 一种行星探测进入段自主导航系统可观测度确定方法 | |
CN110231619B (zh) | 基于恩克法的雷达交接时刻预报方法及装置 | |
CN104050686A (zh) | 一种新型的密集空间目标跟踪方法 | |
CN106323280A (zh) | 用于bds和sins导航定位系统的滤波器和滤波方法 | |
CN106885567A (zh) | 一种惯导协作定位方法及定位设备 | |
CN105549003A (zh) | 一种汽车雷达目标跟踪方法 | |
CN111526487B (zh) | 基于gps和车载测距信息融合的协作车辆定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C41 | Transfer of patent application or patent right or utility model | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20151223 Address after: Espoo, Finland Applicant after: Technology Co., Ltd. of Nokia Address before: Espoo, Finland Applicant before: Nokia Oyj |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |