发明名称 一种用于三轴磁传感器简易校正的两步校正方法
摘要 本发明公开了一种用于三轴磁传感器简易校正的两步校正方法,将方型工装载体做任意角度的旋转,得到一系列地磁场强度三轴分量的测量值。利用基于椭球拟合的三轴磁传感器误差补偿方法,对上述的地磁场强度三轴分量的测量值进行校正,求得误差校正系数矩阵和综合零偏误差。然后将方型工装载体置于水平大理石平台上,令方型工装载体坐标系的X、Y、Z轴分别朝上,让方型工装载体绕朝上的坐标轴缓慢旋转至少一周得到测量数据,利用最小二乘的原理求得正交坐标系转换矩阵。本发明对三轴磁传感器自身坐标系与载体坐标系不重合导致的安装误差进行了有效的校正,同时整个校正过程简捷、省时、精度高,不依赖于精密仪器提供准确的方向基准、水平基准等。
申请公布号 CN106569150A 申请公布日期 2017.04.19
申请号 CN201610945745.4 申请日期 2016.11.02
申请人 南京理工大学 发明人 吴志强;项建梁;姜磊;朱欣华;王宇;苏岩
分类号 G01R33/02(2006.01)I;G01R33/00(2006.01)I 主分类号 G01R33/02(2006.01)I
代理机构 南京理工大学专利中心 32203 代理人 朱沉雁
主权项 一种用于三轴磁传感器简易校正的两步校正方法,其特征在于,方法步骤如下:步骤1、基于椭球拟合对三轴磁传感器进行标定:步骤1‑1、将方型工装载体做任意角度的旋转,从而得到一系列做任意角度旋转时的地磁场强度三轴分量的测量值;步骤1‑2、利用基于椭球拟合的三轴磁传感器误差补偿方法,对上述的地磁场强度三轴分量的测量值进行校正,求得误差校正系数矩阵C<sub>1</sub>和综合零偏误差<img file="FDA0001141160490000018.GIF" wi="75" he="67" />步骤2、基于水平面内三轴旋转对三轴磁传感器的安装误差进行校正:步骤2‑1、将方型工装载体置于水平大理石平台上,并将方型工装载体坐标系的X轴朝上,让方型工装载体绕X轴缓慢旋转至少1周,求得向量[c<sub>11</sub> c<sub>12</sub> c<sub>13</sub>]<sup>Τ</sup>:步骤2‑1‑1、将方型工装载体置于水平大理石平台上,并将方型工装载体坐标系的X轴朝上,让方型工装载体绕X轴缓慢旋转至少1周,得到一系列绕X轴旋转时的地磁场强度三轴分量的测量值;步骤2‑1‑2、规定方型工装载体坐标系与三轴磁传感器自身坐标系并不完全重合,即存在安装误差,因此,存在一个正交的坐标系转换矩阵C<sub>2</sub>,<img file="FDA0001141160490000011.GIF" wi="558" he="221" />使得以下式成立:<img file="FDA0001141160490000012.GIF" wi="1038" he="231" />在式(1)中,<img file="FDA0001141160490000013.GIF" wi="50" he="55" />表示地磁场在方型工装载体坐标系中的真实值,<img file="FDA0001141160490000014.GIF" wi="284" he="70" />分别是<img file="FDA0001141160490000015.GIF" wi="53" he="54" />沿方型工装载体坐标系X、Y、Z轴方向的分量值,<img file="FDA0001141160490000016.GIF" wi="51" he="54" />表示经过步骤1校正后地磁场在磁传感器自身坐标系中的真实值,角度ψ、θ、<img file="FDA0001141160490000017.GIF" wi="35" he="46" />是三个欧拉误差角,c<sub>ij</sub>是正交坐标系转换矩阵C<sub>2</sub>中的第i行、第j列元素,i=1,2,3,j=1,2,3;步骤2‑1‑3、利用步骤1‑2中得到的误差校正系数矩阵C<sub>1</sub>和综合零偏误差<img file="FDA0001141160490000021.GIF" wi="58" he="71" />对步骤2‑1‑1中得到的绕X轴旋转时的地磁场强度三轴分量的测量值进行校正,得到<img file="FDA0001141160490000022.GIF" wi="91" he="55" />则有以下关系式成立:<maths num="0001"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>c</mi><mn>11</mn></msub></mtd><mtd><msub><mi>c</mi><mn>12</mn></msub></mtd><mtd><msub><mi>c</mi><mn>13</mn></msub></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>B</mi><mi>x</mi><mrow><mi>x</mi><mi>s</mi></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>B</mi><mi>y</mi><mrow><mi>x</mi><mi>s</mi></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>B</mi><mi>z</mi><mrow><mi>x</mi><mi>s</mi></mrow></msubsup></mtd></mtr></mtable></mfenced><mo>=</mo><msubsup><mi>B</mi><mi>x</mi><mi>b</mi></msubsup><mo>=</mo><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001141160490000023.GIF" wi="1158" he="236" /></maths>式(2)中<img file="FDA0001141160490000024.GIF" wi="71" he="55" />表示将绕X轴旋转时得到的测量值进过步骤1校正后,地磁场在磁传感器自身坐标系中的真实值,<img file="FDA0001141160490000025.GIF" wi="326" he="63" />分别表示<img file="FDA0001141160490000026.GIF" wi="67" he="55" />沿磁传感器自身坐标系X、Y、Z轴方向上的分量值;步骤2‑1‑4、根据式(2),将步骤2‑1‑3中得到的N<sub>x</sub>个<img file="FDA0001141160490000027.GIF" wi="62" he="53" />展开为方程组如下:<maths num="0002"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>B</mi><mrow><mi>x</mi><mn>1</mn></mrow><mrow><mi>x</mi><mi>s</mi></mrow></msubsup></mtd><mtd><msubsup><mi>B</mi><mrow><mi>y</mi><mn>1</mn></mrow><mrow><mi>x</mi><mi>s</mi></mrow></msubsup></mtd><mtd><msubsup><mi>B</mi><mrow><mi>z</mi><mn>1</mn></mrow><mrow><mi>x</mi><mi>s</mi></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>B</mi><mrow><mi>x</mi><mn>2</mn></mrow><mrow><mi>x</mi><mi>s</mi></mrow></msubsup></mtd><mtd><msubsup><mi>B</mi><mrow><mi>y</mi><mn>2</mn></mrow><mrow><mi>x</mi><mi>s</mi></mrow></msubsup></mtd><mtd><msubsup><mi>B</mi><mrow><mi>z</mi><mn>2</mn></mrow><mrow><mi>x</mi><mi>s</mi></mrow></msubsup></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr></mtable></mtd><mtd><mtable><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr></mtable></mtd><mtd><mtable><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr></mtable></mtd></mtr><mtr><mtd><msubsup><mi>B</mi><mrow><msub><mi>xN</mi><mi>x</mi></msub></mrow><mrow><mi>x</mi><mi>s</mi></mrow></msubsup></mtd><mtd><msubsup><mi>B</mi><mrow><msub><mi>yN</mi><mi>x</mi></msub></mrow><mrow><mi>x</mi><mi>s</mi></mrow></msubsup></mtd><mtd><msubsup><mi>B</mi><mrow><msub><mi>zN</mi><mi>x</mi></msub></mrow><mrow><mi>x</mi><mi>s</mi></mrow></msubsup></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>c</mi><mn>11</mn></msub></mtd></mtr><mtr><mtd><msub><mi>c</mi><mn>12</mn></msub></mtd></mtr><mtr><mtd><msub><mi>c</mi><mn>13</mn></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001141160490000028.GIF" wi="1166" he="309" /></maths>步骤2‑1‑5、由式(3)利用最小二乘拟合,求得向量[c<sub>11</sub> c<sub>12</sub> c<sub>13</sub>]<sup>Τ</sup>:<maths num="0003"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>c</mi><mn>11</mn></msub></mtd></mtr><mtr><mtd><msub><mi>c</mi><mn>12</mn></msub></mtd></mtr><mtr><mtd><msub><mi>c</mi><mn>13</mn></msub></mtd></mtr></mtable></mfenced><mo>=</mo><msup><mrow><mo>(</mo><msup><mi>B</mi><mrow><mi>x</mi><mi>T</mi></mrow></msup><msup><mi>B</mi><mi>x</mi></msup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>B</mi><mrow><mi>x</mi><mi>T</mi></mrow></msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001141160490000029.GIF" wi="1141" he="294" /></maths>其中矩阵<img file="FDA00011411604900000210.GIF" wi="537" he="309" />步骤2‑2、将方型工装载体置于水平大理石平台上,并将方型工装载体坐标系的Y轴朝上,让方型工装载体绕Y轴缓慢旋转至少1周,求得向量[c<sub>21</sub> c<sub>22</sub> c<sub>23</sub>]<sup>Τ</sup>;步骤2‑3、将方型工装载体置于水平大理石平台上,并将方型工装载体坐标系的Z轴朝上,让方型工装载体绕Z轴缓慢旋转至少1周,求得向量[c<sub>31</sub> c<sub>32</sub> c<sub>33</sub>]<sup>Τ</sup>;步骤3、磁传感器测量结果的综合校正模型<img file="FDA0001141160490000031.GIF" wi="54" he="55" />为:<maths num="0004"><math><![CDATA[<mrow><msup><mover><mi>B</mi><mo>&RightArrow;</mo></mover><mi>b</mi></msup><mo>=</mo><msub><mi>C</mi><mn>2</mn></msub><msub><mi>C</mi><mn>1</mn></msub><mrow><mo>(</mo><msup><mover><mover><mi>B</mi><mo>&RightArrow;</mo></mover><mo>^</mo></mover><mi>s</mi></msup><mo>-</mo><msub><mover><mi>B</mi><mo>&RightArrow;</mo></mover><msup><mi>z</mi><mo>&prime;</mo></msup></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001141160490000032.GIF" wi="1214" he="135" /></maths>式(11)中<img file="FDA0001141160490000033.GIF" wi="53" he="70" />是地磁场矢量在磁传感器自身坐标系下的实际测量结果。
地址 210094 江苏省南京市玄武区孝陵卫200号