发明名称 一种室内纯磁场导航的方法
摘要 本发明公开了一种室内纯磁场导航的方法。使用与导航区域一边等长的通电螺线管平行放置在室内的最长边,通上电流产生空间磁场;设定室内导航区域中目标的空间坐标点,计算得到空间目标点两个垂直分量的磁场大小;设定导航物体所在导航区域坐标点为初始坐标点;通过测量经过点的磁场大小信息及基于卡尔曼滤波的导航方法来预测将要行进的下一点的坐标位置,逐步靠近目标点;当在行进至某一点测得的两个垂直分量磁场大小信息与之前设定的目标点计算得到的磁场大小一致,则导航结束,导航物体到达设定的目标点。本发明公开的室内磁场导航方法不需要基于事先存储的磁场信息,不需要与任何其他设备通信,不需要依赖其他导航手段。
申请公布号 CN103344235B 申请公布日期 2015.08.26
申请号 CN201310244459.1 申请日期 2013.06.18
申请人 浙江大学 发明人 赵钟伟;冉立新
分类号 G01C21/00(2006.01)I 主分类号 G01C21/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 林怀禹
主权项 一种室内纯磁场导航的方法,其特征在于:使用与导航区域一边等长的通电螺线管平行放置在室内的最长边,通上电流产生空间磁场;设定室内导航区域中目标的空间坐标点,计算得到空间目标点两个垂直分量的磁场大小;设定导航物体所在导航区域坐标点为初始坐标点;通过测量经过点的磁场大小信息及基于卡尔曼滤波的纯磁场导航方法来预测将要行进的下一点的坐标信息,逐步靠近目标点;当在某一目标点测得的两个垂直分量磁场大小信息与之前设定的目标点计算得到的两个垂直分量的磁场大小一致,则导航结束,导航物体到达设定的目标点;所述的基于卡尔曼滤波的纯磁场导航方法为:(1)根据目标点坐标(x<sub>e</sub>,y<sub>e</sub>)和通电螺线管长度l、半径a、每米的匝数n、通电电流大小I,计算出目标点两个垂直分量的磁场大小信息,记为Bx<sub>e</sub>和By<sub>e</sub>;(2)设定导航物体所在导航区域坐标点为初始坐标点,记为(x<sub>0</sub>,y<sub>0</sub>),用双轴磁传感器测得初始位置的磁场大小信息,即位Bx<sub>0</sub>和By<sub>0</sub>;记初始误差均方值P<sub>0</sub>=0,初始磁场测量值Z<sub>0</sub>=[Bx<sub>0</sub>,By<sub>0</sub>];(3)从初始坐标点开始,一边测量经过路径上的磁场信息一边靠近目标点,整个过程分为k个步进,通过在第k步经过点(x<sub>k</sub>,y<sub>k</sub>)磁场大小及梯度信息来预测行进下一步(第k+1步)即将经过点的空间坐标点位置,记第k步测得的磁场大小信息为Bx<sub>k</sub>和By<sub>k</sub>,记Zm(k)=[Bx<sub>k</sub>,By<sub>k</sub>]<sup>T</sup>,记第k步所在点测得的磁场梯度信息为X<sub>dotk</sub>和Y<sub>dotk</sub>,记第k步X<sub>dotk</sub>和Y<sub>dotk</sub>与y轴的夹角分别为θ<sub>1</sub>和θ<sub>2</sub>,记下一步行进的测量前预测方向与y轴夹角为<img file="FDA0000698459580000011.GIF" wi="121" he="78" />根据方程<img file="FDA0000698459580000012.GIF" wi="926" he="174" />求出;然后求出下一步测量前的坐标预测值<img file="FDA0000698459580000013.GIF" wi="953" he="108" />其中m为每一步测量的步长距离;求出测量前下一步磁场预测值<img file="FDA0000698459580000014.GIF" wi="1681" he="96" />(4)行进至测量前的坐标预测点,根据下一步测量前的坐标预测值<img file="FDA0000698459580000015.GIF" wi="126" he="96" />及该点磁场测量值Z<sub>m</sub>(k+1),计算得到k+1步位置测量后预测状态值<img file="FDA0000698459580000021.GIF" wi="153" he="97" />根据公式<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>M</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msubsup><mover><mi>M</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>-</mo></msubsup><mo>+</mo><msub><mi>K</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>Z</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>Z</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000698459580000022.GIF" wi="1028" he="103" /></maths>其中:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>K</mi><mi>k</mi></msub><mo>=</mo><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><msup><mi>N</mi><mi>T</mi></msup><msup><mrow><mo>(</mo><msubsup><mi>NP</mi><mi>k</mi><mo>-</mo></msubsup><msup><mi>N</mi><mi>T</mi></msup><mo>+</mo><mi>R</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>,</mo></mrow>]]></math><img file="FDA0000698459580000023.GIF" wi="655" he="91" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><mo>=</mo><msub><mi>AP</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msup><mi>A</mi><mi>T</mi></msup><mo>+</mo><mi>Q</mi><mo>,</mo></mrow>]]></math><img file="FDA0000698459580000024.GIF" wi="457" he="88" /></maths>初始时,P<sub>0</sub>=0;A=1;N=Zm(k)对M<sub>k</sub>的导数;Q为处理噪声均方值,设为[100,100;100,100]R为测量噪声均方值,设为[100;100]*[100,100];得到下一步(k+1步)测量后坐标预测点后,行进至该点,重复之前的步骤,以此类推,直到行进至点时,测量到的两个垂直分量的磁场与目标点的两个垂直分量的磁场相同,即认为到达了目标点,导航方法结束。
地址 310027 浙江省杭州市西湖区浙大路38号