发明名称 动磁场下的滤波定位方法及装置
摘要 本发明属于磁场定位技术领域,公开了一种动磁场下的滤波定位方法及装置,本发明所提供的动磁场下的滤波定位方法及装置通过对采集的旋转磁场数据和姿态数据进行分段处理,将每段旋转磁场数据输入到旋转磁场模型进行解算获取每段旋转磁场数据的测量值,对上一段旋转磁场数据的滤波估计实施状态转移,得到上一段旋转磁场数据滤波的一步预测,将一步预测值和该段旋转磁场数据的测量值输入Kalman滤波器进行滤波,获取该段磁场数据的滤波估计,最终计算出目标物体的当前位置,大幅降低了定位结果的方差,具有更高的定位精度,实现在目标物体连续运动过程中的实时测量定位,提高了动磁场定位的抗干扰能力。
申请公布号 CN103630129B 申请公布日期 2016.09.21
申请号 CN201210297928.1 申请日期 2012.08.20
申请人 清华大学 发明人 张超;庞珂珂
分类号 G01C21/00(2006.01)I;G01C21/20(2006.01)I 主分类号 G01C21/00(2006.01)I
代理机构 北京路浩知识产权代理有限公司 11002 代理人 韩国胜
主权项 一种动磁场下的滤波定位方法,其特征在于,包括:100、在目标物体的运动过程中不断采集更新目标物体的旋转磁场数据和姿态数据并存储,同时根据预设规则生成一组时间节点,其中,旋转磁场数据包括x、y、z三个方向的磁场B<sub>x</sub>、B<sub>y</sub>、B<sub>z</sub>及采样时间,姿态数据包括俯仰角<img file="FDA00010129509100000111.GIF" wi="93" he="62" />方位角θ<sub>s</sub>及采样时间;101、k时刻,根据所述一组时间节点分别将存储的旋转磁场数据和姿态数据分为M段,分别记为s<sub>m</sub>和h<sub>m</sub>,其中,1≤m≤M,m、M为自然数;102、将第m段旋转磁场数据s<sub>m</sub>输入旋转磁场模型进行解算,得到第m个测量值<img file="FDA0001012950910000011.GIF" wi="604" he="127" />其中k为k时刻,<img file="FDA0001012950910000012.GIF" wi="74" he="71" />为k时刻第m个测点的相对距离,<img file="FDA0001012950910000013.GIF" wi="76" he="71" />为k时刻第m个测点的俯仰角,<img file="FDA0001012950910000014.GIF" wi="70" he="62" />为k时刻第m个测点的方位角;103、当2≤m≤M时,根据第m段姿态数据h<sub>m</sub>和第m个测量值<img file="FDA0001012950910000015.GIF" wi="94" he="79" />对第m‑1次Kalman滤波估计<img file="FDA0001012950910000016.GIF" wi="667" he="127" />进行状态转移,得到第m‑1次滤波估计的一步预测<img file="FDA0001012950910000017.GIF" wi="918" he="127" />即<img file="FDA0001012950910000018.GIF" wi="1254" he="278" /><img file="FDA0001012950910000019.GIF" wi="774" he="254" /><img file="FDA00010129509100000110.GIF" wi="1590" he="254" />或<img file="FDA0001012950910000021.GIF" wi="1377" he="255" />其中,<img file="FDA0001012950910000022.GIF" wi="1030" he="175" />;104、将第m次测量值和第m‑1次滤波估计的一步预测输入Kalman滤波方程进行第m次Kalman滤波,得到第m次滤波估计<img file="FDA0001012950910000023.GIF" wi="126" he="127" />其中,当m=1时,则将k‑1时刻第M次滤波估计作为k时刻第0次滤波估计的一步预测;105、重复102、103和104步骤,并最终得到k时刻的第M次滤波估计<img file="FDA0001012950910000024.GIF" wi="606" he="131" />则k时刻目标物体的绝对空间坐标为:<img file="FDA0001012950910000025.GIF" wi="542" he="95" /><img file="FDA0001012950910000026.GIF" wi="524" he="95" /><maths num="0001"><math><![CDATA[<mrow><msup><mi>z</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>=</mo><msubsup><mi>&gamma;</mi><mi>M</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><msubsup><mi>sin&theta;</mi><mi>M</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><mo>.</mo></mrow>]]></math><img file="FDA0001012950910000027.GIF" wi="370" he="92" /></maths>
地址 100084 北京市海淀区清华园北京100084-82信箱