CN111183335A - 表征由磁场测量值确定的航向的方法和设备 - Google Patents

表征由磁场测量值确定的航向的方法和设备 Download PDF

Info

Publication number
CN111183335A
CN111183335A CN201880063738.9A CN201880063738A CN111183335A CN 111183335 A CN111183335 A CN 111183335A CN 201880063738 A CN201880063738 A CN 201880063738A CN 111183335 A CN111183335 A CN 111183335A
Authority
CN
China
Prior art keywords
heading
magnetic
magnetic field
determined
magnetometer
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
Application number
CN201880063738.9A
Other languages
English (en)
Other versions
CN111183335B (zh
Inventor
大卫·维西埃
马蒂厄·希尔伊恩
查尔斯-伊凡·切斯诺
让-菲力浦·杜加尔德
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi Sinawei
Original Assignee
Xi Sinawei
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xi Sinawei filed Critical Xi Sinawei
Publication of CN111183335A publication Critical patent/CN111183335A/zh
Application granted granted Critical
Publication of CN111183335B publication Critical patent/CN111183335B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C17/00Compasses; Devices for ascertaining true or magnetic north for navigation or surveying purposes
    • G01C17/38Testing, calibrating, or compensating of compasses
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/04Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
    • G01C21/08Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means involving use of the magnetic field of the earth
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/165Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
    • G01C21/1654Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments with electromagnetic compass
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations
    • G01C21/206Instruments for performing navigational calculations specially adapted for indoor navigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Electromagnetism (AREA)
  • Measuring Magnetic Variables (AREA)

Abstract

本发明涉及一种用于确定航向的方法,其中磁场由至少一个磁力计测量,并且磁航向由所述测量值确定,并且其中计算装置针对磁航向的每次测量实施对由此确定的航向的表征值进行计算的处理,其特征在于,所述表征值是测量的磁场梯度的范数的单调函数,所述值用于确定航向。

Description

表征由磁场测量值确定的航向的方法和设备
一般技术领域和现有技术
总体而言,本发明涉及磁惯性技术。
更准确地说,本发明涉及通过磁力计来确定航向(heading)。
特别是,本发明有利地应用于在城市地区或“室内”,即在建筑物内部进行测量的情况。
传统上,磁力计用于计算嵌入式系统中的航向。
在这种情况下,假设传感器测量的磁场是地球磁场,其水平分量指向磁北。磁北和地理北的方向之间的差异(称为磁偏角)是已知的,并制成表格。在不失一般性的情况下,此后将因此认为磁北和地理北是合并的,并且因此磁力计指示了什么将被称为北。
通常,用于根据磁力计测量值计算航向的确定方法基于以下:
-将磁力计提供的测量值和与系统航向相关的信息联系起来的建模;
-这种关系对于所执行的测量的相关性的表征。
在常规方法中,所实现的计算可以是以磁场为测度的卡尔曼类型滤波,使得可以利用作为相关性表征的函数的测量噪声来修正包含在状态中的航向。
建模包括写入测量的磁场M包含航向信息,例如通过以下等式:
Figure BDA0002431191460000011
其中,R是旋转矩阵,使得可以从物体的参考系转到地球的惯性参考系成为可能。
Figure BDA0002431191460000012
是欧拉角,并且MEARTH是地球磁场。
该等式的相关性以测量方差为特征,也就是说,假设这个等式中的误差是具有零期望值的高斯随机变量。
该方差用于自动计算卡尔曼增益,所述卡尔曼增益对修正进行加权,以便考虑不同的噪声(与外部环境相关的动态噪声和磁力计的测量噪声)。
在又一种方法中,实施使用相同类型的建模和表征通过线性滤波进行的计算。
在这种方法中,表征等式相关性的是等式增益幅度的相对调整。一般来说,这个增益是手动调整的。本领域技术人员将知道如何根据建模的相关性对其进行加权。
因此,在这两种方法中,表征模型相关性的参数(在卡尔曼滤波器的情况下为高斯噪声的方差,在线性滤波的情况下为增益)通常是恒定的参数,与测量的磁场无关。
最近,已经提出了对考虑了所测量磁场的建模的相关性的表征,例如在出版物中:
-W.T.Faulkner,R.Alwood,W.T.David和J.Bohlin,“Gps-denied pedestriantracking in indoor environments using an imu and magnetic compass(在室内环境中使用惯性测量单元和磁罗盘的拒绝GPS的行人跟踪),”在Proceedings of the 2010International Technical Meeting of The Institute of Navigation(2010年导航学会国际技术会议会议录)中,(San Diego(圣地亚哥),CA),第198–204页,
-M.H.Afzal,V.Renaudin和G.Lachapelle,“Magnetic field based headingestimation for pedestrian navigation environments(基于磁场的行人导航环境航向估计),”在2011International Conference on Indoor Positioning and IndoorNavigation(2011年室内定位和室内导航国际会议)中,(Guimaraes,Portugal)(葡萄牙,圭马拉斯)M.H.Afzal,V.Renaudin,and G.Lachapelle,2011年9月,
然而,由此进行的表征并没有完全令人满意。
在磁场受到相当大的干扰的环境中尤其如此,特别是在“室内”或城市地区,其中起源于建筑物中的金属元素或电流的磁干扰会叠加在地球磁场上。
发明内容
本发明的总体目标是提出一种能够更好地表征所使用的建模的相关性的解决方案,特别是在高度干扰的环境中。
值得注意的是,本发明提出了一种用于确定航向的方法,其中磁场由至少一个磁力计测量,并且磁航向由所述测量值确定,并且其中计算装置对于磁航向的每次测量都实施对由此确定的航向的表征值进行计算的处理,其特征在于,所述表征值是测量的磁场梯度的函数,所述值用于确定航向。
特别地,所述表征值是梯度的范数(特别是范数L2)的单调函数,优选地递增的单调函数(将采取直接正系数的线性函数的示例)。
应当理解,如果例如在航向受到越来越多干扰时,任意确定所述表征值具有越来越低的值,则所述表征值可以是梯度范数的递减函数。
计算装置对于给定的采样时间确定:
-对航向预测的估计,所述估计是在先前的采样时间确定的航向的函数,
-根据从磁力计测量值确定的磁航向对所估计的预测航向进行的修正。
并且其中修正根据所述表征值确定,并且有利地仍然是表征值的单调函数(尤其是递增函数)。
因此,所述修正有利地是所测量磁场的梯度的函数,并且特别是梯度范数的单调函数(通过单调函数的组合)。以这种方式,考虑了环境被干扰越多,梯度越高,因此航向测量值被篡改就越多的事实。
在优选实施例中,计算装置根据以下来估计修正的幅度:
-两个采样时间之间磁干扰变化的模型,以及
-干扰幅度的先验估计。
以这种方式,处理考虑了两个采样时间之间的磁航向干扰的变化,并且考虑了在连续采样时间的情况下由于环境而引起的干扰的潜在空间相关性。
由此产生的计算比现有技术更可靠。
值得注意的是,对于给定的采样时间k+1,计算装置通过计算以下公式来估计相关联的航向干扰:
Figure BDA0002431191460000041
其中
Figure BDA0002431191460000042
并且vψ (d)是方差的高斯随机变量
a[k]2(1-α[k]2)
-k是先前的采样时间,
-a[k]是表示磁干扰的先验幅度的参数,
-
Figure BDA0002431191460000043
和σu是由计算装置计算的作为对干扰变化的期望值和方差的估计的两个参数。
有利地,计算装置将参数a[k]估计为磁场梯度的范数的线性函数。
参数
Figure BDA0002431191460000044
本身例如由计算装置计算为:直接从对于时间k+1和k的磁力计输出确定的磁航向之间的差值,从该差值中减去预测航向(ω[k]dt)的变化量。
参数σu可以被估计为两个连续采样时间之间的位移或位移速度的函数。
该处理实现例如卡尔曼滤波,该卡尔曼滤波的状态至少具有实际航向和磁航向干扰(ψ,ψ(d))来作为参数。
航向预测可以根据惯性单元的一个或多个传感器的测量值来确定。
本发明还涉及一种用于通过磁传感器确定航向的设备,该设备包括磁力计和用于根据由所述磁力计测量到的磁场来计算航向的计算装置,其特征在于,所述计算装置针对不同的连续采样时间实施上述处理。
还提出了一种包括至少一个这样的航向测量设备的磁惯性导航系统。
这种系统有利地用于城市环境或建筑物内部。
本发明还涉及:
-一种计算机程序产品,包括用于当所述程序在计算机上执行时执行前述类型的方法的代码指令;
-一种计算机设备可读的存储装置,其上的计算机程序产品包括用于执行这种方法的代码指令。
附图说明
从下面的描述中,本发明的其他特征和优点将变得更加清楚,这些描述纯粹是说明性的而非限制性的,并且应该参照附图来阅读,在附图中:
-图1是用于实施根据本发明的方法的设备图;
-图2更详细地表示了根据本发明的方法的实施的示例;
-图3示出了符合本发明实施例的方法的主要步骤。
具体实施方式
一般性说明–测量设备
参考图1,所提出的测量设备例如用于估计在环境磁场(通常是地球磁场,其可能会被附近的金属物体轻微改变)中移动的物体1的运动,标注为
Figure BDA0002431191460000051
众所周知,磁场是三维矢量场,也就是说,将三维矢量与其中物体可移动的每个三维点相关联。
该物体1可以是需要位置信息的任何可移动物体,例如轮式车辆、无人机等,而且也可以是行人。
物体1在壳体2(支架)中包括多个磁测量传感器20,即轴向磁力计20。轴向磁力计是指能够测量所述磁场的分量(例如,所述磁场矢量
Figure BDA0002431191460000052
在所述磁力计20的水平上沿其轴的投影)的元件。
更准确地说,磁力计20与壳体2是一体的。它们的运动在地面参考系中与壳体2和物体1基本相同。
以优选的方式,物体1的参考系设置有正交笛卡尔参考点,其中磁力计20在该参考点上具有预定位置。
在图2中,壳体2通过附接装置23被固定在物体1(例如行人的肢体)上。这些附接装置23例如由手镯,例如夹紧肢体并实现整体连接的自夹紧带(self-gripping strap)组成。
显然,本发明不限于估计行人的运动,而是在这种用途中特别有利,因为它能够大大减小体积,这对于人体以人体工程学的方式而便携的壳体是必需的。
壳体2可以包括用于直接实时实现本方法的处理操作的计算装置21(通常是处理器),或者取而代之的是,将测量值经由通信装置25传输到诸如移动终端(智能手机)3之类的外部设备,或者甚至是远程服务器4,或者取而代之的是,可以将测量值记录在本地数据存储存储器装置22(例如,闪存型存储器)中,用于例如服务器4上的后验处理。
通信装置25可以实现短程无线通信,例如蓝牙或Wi-Fi(特别是在具有移动终端3的实施例中),或者甚至是用于连接到移动网络(通常为UMTS/LTE)以进行长距离通信的装置。应当注意,通信装置25可以是例如有线连接件(通常是USB),用于将数据从本地数据存储装置22传递到移动终端3或服务器4的数据存储装置。
如果是承载“智能”的移动终端3(对应的是服务器4),则其包括计算装置31(对应的是41),诸如用于实现将要描述的本方法的处理操作的处理器。当所使用的计算装置是壳体2的计算装置21时,它还可以包括用于传送估计位置的通信装置25。例如,佩戴者的位置可以被发送到移动终端3以在导航软件界面中显示该位置。
壳体2、智能手机3和远程服务器4的数据计算装置21、31、41可以无区别地并且根据应用执行该方法的全部或部分步骤。
为此,它们每个都包括存储装置,其中存储了用于执行该方法的全部或部分代码指令序列。
预测和修正
计算装置实现(图3)滤波100,滤波100一方面通过预测来计算对航向值的估计(步骤101),另一方面根据误差估计来实现修正(步骤102)。
值得注意的是,已知例如由陀螺仪给出的角速度ω,步骤101在时间k+1计算航向ψk+1等于
ψk+1=ψk+ω.Δt
其中ψk是前一时间k的航向,并且其中Δt是分隔这两个采样时间的持续时间。
修正102考虑了磁力计20执行的测量。
此后,磁航向针对给定时间k的测量值(从磁力计的测量值导出)记为zψ[k],其中
zψ[k]=ψk
通常,磁航向由下式给出
Figure BDA0002431191460000071
其中,My和Mx是地球参考系中磁场的两个水平分量,这两个分量是根据物体1的惯性单元为物体1确定的姿态来计算的。
计算装置将与时间k+1相对应的修正航向
Figure BDA0002431191460000072
计算为等于此时所估计的航向ψk+1和修正的总和,所述修正有利地是以下的函数:
-相对于先前时间K计算或调整的增益Kk
-预测航向ψk+1和测量值zψ之间误差Err(ψk+1,zψ)的估计值
因此,计算装置使用从磁力计的测量值中导出的磁航向,以通过计算
Figure BDA0002431191460000073
来修正状态,特别是航向。
通常,误差Err(ψk+1,zψ)可能是在时间k+1(修正之前)估计的预测航向和从磁力计测量值中导出的磁航向zψ之间的简单差值。
然而,其他误差函数也是可能的,特别是在非线性滤波的情况下。
特别地,修正Kk.Err(ψk+1,zψ),并且更特别地,增益Kk可以有利地取决于由磁力计20测量的磁场梯度。如前所述,增益Kk优选地是梯度范数的单调函数,尤其是递增的单调函数。
这样,当磁场变化很大并因此很容易在航向测量中引起重大误差时,修正校正尤为重要。
由此获得的修正航向值由所述计算装置21、31、41存储和/或由所述装置用于其余的处理以及用于磁惯性导航信息(线速度、角速度、位置、航向等)的计算。
此外,它可以由计算装置传输到接口装置,例如在电话上以显示在其屏幕上。
与环境相关的干扰变化模型
下面详述一个可能的修正计算的示例,其中通过估计与环境相关的磁场干扰来校正磁航向zψ的测量。
从磁力计的测量值中导出的磁航向zψ实际上可以认为被分解如下:
Figure BDA0002431191460000081
其中:
-ψ对应于旨在确定的实际航向,
-
Figure BDA0002431191460000082
对应于高斯测量误差,以及
(d)(“d”代表干扰)对应于与环境有关的磁航向干扰(通常是与城市环境或建筑物中的金属基础设施和电缆有关的干扰)。
与环境相关联的磁航向干扰在空间上高度相关。与环境干扰相关联的磁场实际上是一个连续的矢量场,并且当两个给定的点A和B在空间上相邻时,这两个给定点A和B处的磁场都更接近。
在一个实施例中,与环境相关联的航向干扰ψ(d)由计算装置借助于以下公式来估计:考虑了两个连续采样时间之间的磁场干扰可能预期的变化相关性(时间、空间相关性或更复杂)的公式。
构建此估计使得实现用于ψd变化的Markovian模型(可用于递归滤波器),以及
·使得有可能考虑磁航向干扰变化(时间、空间相关性或更复杂)的模型,
·使得构建航向可观测的滤波模型成为可能(否则,构建航向估计器的任何希望都将落空)
在所提出的估计中,修正的幅度被估计为两个采样时间之间磁干扰变化的模型和干扰幅度的先验估计的函数。
发明人已经发现(并且在数学上得到验证),在空间相关性的情况下(在城市地区或室内的干扰的情况下),由时间步长分开的两个采样时间之间的干扰ψ(d)变化的合适的估计有利地如下:
Figure BDA0002431191460000091
其中
Figure BDA0002431191460000092
并且vψ(d)是方差的随机高斯变量
a[k]2(1-α[k]2)
-a[k]是代表磁干扰的先验幅度的参数,
-
Figure BDA0002431191460000093
和σu是随机变量的期望值和方差,该随机变量是两个采样时间之间磁航向干扰变化的模型(以随机高斯变量的形式),
-vψ(d)是考虑到磁力计噪声的估计。
该估计值可从zψ中减去以进行误差计算,这可以在ψk+1和(zψ(d)[k+1])之间直接进行计算,或者限制于在其状态下既有ψ(d)又有ψ的卡尔曼型滤波器中。
参数a[k]的确定
参数a[k]表示磁航向干扰的先验幅度。它表征了航向,并且例如被计算为磁场梯度范数的线性函数,其是单调函数的一个示例。
例如
a[k]=a0+a1N[k]
其中,
Figure BDA0002431191460000101
即范数L2。然而,也可以使用其他范数,诸如范数L1、Frobenius范数或任何其他范数。
a0和a1是在实施滤波处理之前预先固定的两个参数,
Figure BDA0002431191460000102
是时间k处的磁场。从a1为正的那一刻起,函数a[k]就是范数的递增函数。
可替换地,可以使用其他单调函数,例如从梯度分量中适当选择的参数的二次函数。
以此方式,如上所述,修正项Kk.Err(ψk+1,zψ)是增益Kk的单调函数,其本身是由磁力计20测量的磁场梯度的范数的单调函数,这使得能够进行高效的校正。
其他方法也可以确定参数a[k]。
值得注意的是,a[k]也可以通过与地球磁场模型进行比较来确定,例如通过实施诸如以下出版物中提出的技术来确定:
“Assessment of Indoor Magnetic Field Anomalies using MultipleMagnetometers Assessment of Indoor(使用室内的多台磁力计评估来评估室内磁场异常)”-M.H.Afzal,V.Renaudin,G.Lachapelle-Proceedings of the 23rd InternationalTechnical Meeting of The Satellite Division of the Institute of Navigation(ION GNSS 2010)(第23届导航研究所卫星分部国际技术会议会议录(ION GNSS 2010)),-M.H.Afzal,V.Renaudin,G.Lachapelle-2010年9月21日至24日,,2010年9月21日至24日。
其中参数a[k]是磁场梯度的函数的解决方案可以实现更好的校正,并且不需要地球磁场的模型。
参数
Figure BDA0002431191460000111
的确定
在第一近似中,术语vψ(d)可以被认为是可忽略的,并且在计算中被认为是零。
参数
Figure BDA0002431191460000112
代表干扰变化的最可能值。
例如,使用在短期内假定正确的陀螺仪来计算该参数。然后,通过磁航向测量值的变化与陀螺仪航向的变化之间的差来给出干扰最可能的变化:
Figure BDA0002431191460000113
其中:E[.]是数学期望,ω[k]是陀螺仪的转速,并且dt是采样间距。
因此,参数
Figure BDA0002431191460000114
被计算装置计算为直接从对于时间k+1和k的磁力计输出确定的磁航向之间的差值,从该差值中减去预测旋转量ω[k]dt。
参数σu的确定
参数σu表示在步骤k处的干扰与步骤k+1处的干扰之间的相关性的图像。
可以有利地选择以考虑到磁干扰的空间相关性的这种方式根据位移速度(或两个连续采样时间之间的位移,这是相似的)对其进行索引。
可以举个示例
σu=c.||v||,其中c是调整系数。
其他实施例
在另一个实施例中,组合不同的磁航向值,以便获得取决于航向的变量值,该变量值尽可能少地受到磁干扰相关联的误差的损害。
旨在确定的变量例如是轨迹的一部分的平均航向,其他航向相关变量显然能够被设想。
为此,计算装置针对每个采样间距k记录所测量的磁航向的值zψ[k]和称为a[k]的磁航向的表征值,该表征值被计算为梯度范数的单调函数。
例如,
Figure BDA0002431191460000121
因此,可以以优化的形式提出在两个采样时间k=k0和k=k1之间定义的轨迹部分的最可能航向
Figure BDA0002431191460000122
Figure BDA0002431191460000123
其中||.||是范数(L1,L2或其他)。
由计算装置实现的优化算法可以是任何合适的类型:梯度下降法、Levenberg-Marquardt算法或Metropolis-Hastings算法。
此外,在范数L2的情况下,然后以最小二乘法的形式提出问题,并且计算装置直接确定
Figure BDA0002431191460000131
容易理解的是,该值比未通过表征值加权的平均值更具优势,因为磁干扰最强的时间对加权平均值的影响要比未加权平均值的影响小得多。

Claims (15)

1.一种用于确定航向的方法,其中通过至少一个磁力计测量磁场,并且根据所述测量值确定磁航向,并且其中计算装置针对磁航向的每次测量,实施对由此确定的航向的表征值进行计算的处理,其特征在于,所述表征值是所测量的磁场梯度的范数的单调函数,所述值被用于确定所述航向。
2.根据权利要求1所述的方法,其中,对于给定的采样时间,所述计算装置确定:
-对航向预测进行估计,所述估计是在前一个采样时间确定的航向的函数,
-根据从磁力计测量值确定的磁航向,对所估计的预测航向进行修正,
并且其中所述修正根据所述表征值确定。
3.根据权利要求1和2中的一项所述的方法,其中,所述计算装置根据多个测量的磁航向和表征值来确定取决于所述航向的变量值。
4.根据权利要求1至3中的一项所述的方法,其中,在所述修正期间,所述计算装置根据以下来估计所述修正的幅度:
-两个采样时间之间的磁干扰的变化的模型,以及
-对干扰幅度的先验估计。
5.根据权利要求2所述的方法,其中,对于给定的采样时间k+1,所述计算装置通过计算以下公式来估计关联的航向干扰:
Figure FDA0002431191450000011
其中,
Figure FDA0002431191450000012
其中,vψ(d)是方差a[k]2(1-α[k]2)的高斯随机变量,
k是前一个采样时间,
a[k]是代表所述磁干扰的先验幅度的参数,
Figure FDA0002431191450000021
和σu是由所述计算装置计算的作为对干扰变化的期望值和方差的估计的两个参数。
6.根据权利要求5的方法,其中,所述计算装置将所述参数a[k]估计为所述磁场梯度的范数的线性函数。
7.根据权利要求5和6中的一项所述的方法,其中,由所述计算装置将所述参数
Figure FDA0002431191450000022
计算为:直接从对于时间k+1和k的磁力计输出确定的所述磁航向之间的差值,从所述差值中减去预测的航向变化量(ω[k]dt),其中ω[k]是回转测量旋转速度,并且dt是采样间距。
8.根据权利要求5到7中的一项所述的方法,其中,所述参数σu被估计为两个连续采样时间之间的位移或位移速度的函数。
9.根据权利要求5至8中的一项所述的方法,其中,所述处理实现卡尔曼滤波,所述卡尔曼滤波的状态至少具有真实航向和磁航向干扰(ψ,ψ(d))来作为参数。
10.根据前述权利要求中的一项所述的方法,其中,所述航向预测根据惯性单元的一个或多个传感器的测量值来确定。
11.一种用于通过磁传感器确定航向的设备,包括磁力计和用于根据由所述磁力计测量到的磁场来计算航向的计算装置,其特征在于,所述计算装置对于不同的连续采样时间实施根据前述权利要求中的一项所述的方法的处理。
12.一种磁惯性导航系统,包括至少一个根据权利要求11所述的磁场测量设备。
13.一种根据权利要求12所述的系统在城市环境或建筑物内的应用中的用途。
14.一种计算机程序产品,包括当在计算机上执行所述程序时用于执行根据权利要求1至10中的一项所述的磁场测量方法的代码指令。
15.一种计算机设备可读的存储装置,其上的计算机程序产品包括用于执行根据权利要求1至10中的一项所述的磁场测量方法的代码指令。
CN201880063738.9A 2017-07-28 2018-07-27 一种确定航向的方法和通过磁传感器确定航向的设备 Active CN111183335B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR1757220 2017-07-28
FR1757220A FR3069634B1 (fr) 2017-07-28 2017-07-28 Procede et dispositif de caracterisation d'un cap determine a partir de la mesure du champ magnetique
PCT/FR2018/051940 WO2019020961A1 (fr) 2017-07-28 2018-07-27 Procédé et dispositif de caractérisation d'un cap déterminé à partir de la mesure du champ magnétique

Publications (2)

Publication Number Publication Date
CN111183335A true CN111183335A (zh) 2020-05-19
CN111183335B CN111183335B (zh) 2022-02-11

Family

ID=60382327

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201880063738.9A Active CN111183335B (zh) 2017-07-28 2018-07-27 一种确定航向的方法和通过磁传感器确定航向的设备

Country Status (7)

Country Link
US (1) US11415420B2 (zh)
EP (1) EP3658853B1 (zh)
JP (1) JP7404227B2 (zh)
KR (1) KR102547926B1 (zh)
CN (1) CN111183335B (zh)
FR (1) FR3069634B1 (zh)
WO (1) WO2019020961A1 (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114046785A (zh) * 2021-11-10 2022-02-15 广东微电科技有限公司 磁探测信号线性噪声滤波方法及系统、计算机可读存储介质、磁导航传感器、agv小车
FR3137764B1 (fr) 2022-07-08 2024-07-12 Sysnav Procédé de recalage d’une attitude fournie par un système de navigation à l’estime au moyen d’un système de positionnement relatif

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5319561A (en) * 1991-01-21 1994-06-07 Sumitomo Electric Industries, Ltd. Apparatus for estimating current heading using magnetic and angular velocity sensors
EP2264485A1 (en) * 2009-06-04 2010-12-22 Tracedge Method for analysing moving object continuous trajectory based on sampled GPS position
CN102252689A (zh) * 2010-05-19 2011-11-23 北京国浩传感器技术研究院(普通合伙) 一种基于磁传感器的电子罗盘校准方法
CN202177381U (zh) * 2011-05-06 2012-03-28 微迈森惯性技术开发(北京)有限公司 运动姿态数据获取装置与人体运动姿态追踪系统
JP5061264B1 (ja) * 2012-03-23 2012-10-31 国立大学法人 千葉大学 小型姿勢センサ
CN103791905A (zh) * 2012-10-30 2014-05-14 雅马哈株式会社 姿态估计方法和装置
CN103994766A (zh) * 2014-05-09 2014-08-20 北京航空航天大学 一种抗gps失效固定翼无人机定向方法
CN106017459A (zh) * 2016-05-19 2016-10-12 极翼机器人(上海)有限公司 一种抗磁干扰的方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140297212A1 (en) * 2013-04-02 2014-10-02 Kionix, Inc. Systems and Methods for Compensating for a Misalignment Angle Between an Accelerometer and a Magnetometer
JP6372751B2 (ja) 2014-09-22 2018-08-15 カシオ計算機株式会社 電子機器及びオフセット値取得方法、オフセット値取得プログラム
JP6384662B2 (ja) 2014-09-22 2018-09-05 カシオ計算機株式会社 電子機器及びセンサ較正方法、センサ較正プログラム
EP3325990B1 (en) * 2015-07-21 2022-12-21 Israel Aerospace Industries Ltd. Gradiometer system and method
US9958274B2 (en) 2015-12-29 2018-05-01 Fujitsu Limited Magnetic field based location determination

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5319561A (en) * 1991-01-21 1994-06-07 Sumitomo Electric Industries, Ltd. Apparatus for estimating current heading using magnetic and angular velocity sensors
EP2264485A1 (en) * 2009-06-04 2010-12-22 Tracedge Method for analysing moving object continuous trajectory based on sampled GPS position
CN102252689A (zh) * 2010-05-19 2011-11-23 北京国浩传感器技术研究院(普通合伙) 一种基于磁传感器的电子罗盘校准方法
CN202177381U (zh) * 2011-05-06 2012-03-28 微迈森惯性技术开发(北京)有限公司 运动姿态数据获取装置与人体运动姿态追踪系统
JP5061264B1 (ja) * 2012-03-23 2012-10-31 国立大学法人 千葉大学 小型姿勢センサ
CN103791905A (zh) * 2012-10-30 2014-05-14 雅马哈株式会社 姿态估计方法和装置
CN103994766A (zh) * 2014-05-09 2014-08-20 北京航空航天大学 一种抗gps失效固定翼无人机定向方法
CN106017459A (zh) * 2016-05-19 2016-10-12 极翼机器人(上海)有限公司 一种抗磁干扰的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
AFZAL,MH ET.AL.: "Assessment of Indoor Magnetic Field Anomalies using Multiple Magnetometers", 《 INSTITUTE OF NAVIGATION SATELLITE DIVISION PROCEEDINGS OF THE INTERNATIONAL TECHNICAL MEETING》 *
BERGAMINI, ELENA 等: "Estimating Orientation Using Magnetic and Inertial Sensors and Different Sensor Fusion Approaches: Accuracy Assessment in Manual and Locomotion Tasks", 《SENSORS》 *
牛德森 等: "高动态精度抗磁干扰胶囊内镜定位方法", 《激光与光电子学进展》 *

Also Published As

Publication number Publication date
EP3658853B1 (fr) 2024-06-05
FR3069634A1 (fr) 2019-02-01
JP2020528560A (ja) 2020-09-24
FR3069634B1 (fr) 2021-06-11
KR20200037333A (ko) 2020-04-08
WO2019020961A1 (fr) 2019-01-31
US11415420B2 (en) 2022-08-16
EP3658853C0 (fr) 2024-06-05
KR102547926B1 (ko) 2023-06-26
CN111183335B (zh) 2022-02-11
US20200208980A1 (en) 2020-07-02
JP7404227B2 (ja) 2023-12-25
EP3658853A1 (fr) 2020-06-03

Similar Documents

Publication Publication Date Title
CN110887480B (zh) 基于mems传感器的飞行姿态估计方法及系统
JP7025215B2 (ja) 測位システム及び測位方法
CN106370178B (zh) 移动终端设备的姿态测量方法及装置
CN105716610A (zh) 一种地磁场模型辅助的载体姿态和航向计算方法和系统
CN113465628A (zh) 惯性测量单元数据补偿方法及系统
JP2021527214A (ja) 物体の磁力計を較正する方法
CN111183335B (zh) 一种确定航向的方法和通过磁传感器确定航向的设备
US20170074689A1 (en) Sensor Fusion Method for Determining Orientation of an Object
CN111148966B (zh) 一种用磁传感器确定航向的方法和设备
US11519728B2 (en) Method for estimating the movement of an object moving in a magnetic field
EP3227634B1 (en) Method and system for estimating relative angle between headings
JP2013122384A (ja) カルマンフィルタ、及び、状態推定装置
Qian et al. Optical flow based step length estimation for indoor pedestrian navigation on a smartphone
Qian et al. Optical flow-based gait modeling algorithm for pedestrian navigation using smartphone sensors
WO2021245640A1 (en) Magnetic indoor positioning with magnetometer calibration errors equalization
CN111133274B (zh) 用于估计在环境和磁场中运动的对象的运动的方法
JPWO2019020962A5 (zh)
CN110579212B (zh) 室内定位方法及装置
Kuusniemi et al. Evaluation of Bayesian approaches for multi-sensor multi-network seamless positioning
Yang et al. Magnetic Distortion Resistant Orientation Estimation
CN115690909A (zh) 变电站施工人员自主定位方法及装置
CN114322991A (zh) 设备的姿态检测方法、装置、电子设备及计算机存储介质

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