发明名称 用于定位水下运动磁目标的磁传感器阵列校准方法
摘要 本发明公开了一种用于定位水下运动磁目标的磁传感器阵列校准方法。所述校准方法,包括N个已校准三维量的磁传感器,按照一定的阵型布放于海底;各个磁传感器的检测数据由海底电缆传输至水下或水上基站进行处理;在每个磁传感器上集成有压力传感器和罗经,磁传感器和罗经集成后,同时完成罗经与磁传感器坐标系角度差的校准。通过校准每个磁传感器的空间坐标以及三维姿态,给出布放的磁传感器的阵列布局,从而达到较宽海域的磁目标监测,对运动磁目标进行精度更高的定位。
申请公布号 CN103792595A 申请公布日期 2014.05.14
申请号 CN201410069499.1 申请日期 2014.02.27
申请人 中国电子科技集团公司第四十一研究所 发明人 杨保国;韩顺利
分类号 G01V13/00(2006.01)I 主分类号 G01V13/00(2006.01)I
代理机构 济南舜源专利事务所有限公司 37205 代理人 王连君
主权项 1.用于定位水下运动磁目标的磁传感器阵列校准方法,其特征在于,包括N个已校准三维量的磁传感器,按照一定的阵型布放于海底;各个磁传感器的检测数据由海底电缆传输至水下或水上基站进行处理,N为整数,且N≥2;在每个磁传感器上集成有压力传感器和罗经,磁传感器和罗经集成后,完成罗经与磁传感器坐标系角度差的校准;所述校准方法包括如下校准步骤:a磁传感器的姿态校准设磁目标的坐标为<img file="FDA0000470627490000011.GIF" wi="326" he="90" />第i个磁传感器的大地坐标为<img file="FDA0000470627490000012.GIF" wi="314" he="90" />1≤i≤N,磁传感器三维姿态角为<img file="FDA0000470627490000013.GIF" wi="318" he="88" />对磁目标的定向数据方位角为<img file="FDA0000470627490000014.GIF" wi="334" he="89" />定向坐标系为(X′<sub>a</sub>,Y′<sub>a</sub>,Z′<sub>a</sub>);磁目标在以磁传感器为原点的大地坐标系的方位角为<img file="FDA0000470627490000015.GIF" wi="328" he="89" />坐标系为(X<sub>a</sub>,Y<sub>a</sub>,Z<sub>a</sub>);上述两个坐标系的旋转矩阵表示为:<maths num="0001"><![CDATA[<math><mrow><mi>R</mi><mrow><mo>(</mo><mi>&Psi;</mi><mo>)</mo></mrow><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mi></mi><mi>&kappa;</mi><mi>cos</mi><mi>A</mi></mtd><mtd><mi>cos</mi><mi></mi><mi>&kappa;</mi><mi>sin</mi><mi>A</mi></mtd><mtd><mi>sin</mi><mi>&kappa;</mi></mtd></mtr><mtr><mtd><mo>-</mo><mi>cos</mi><msup><mi>&phi;</mi><mo>`</mo></msup><mi>sin</mi><mi>A</mi><mo>-</mo><mi>sin</mi><msup><mi>&phi;</mi><mo>`</mo></msup><mi>sin</mi><mi></mi><mi>&kappa;</mi><mi>cos</mi><mi>A</mi></mtd><mtd><mi>cos</mi><msup><mi>&phi;</mi><mo>`</mo></msup><mi>cos</mi><mi>A</mi><mo>-</mo><mi>sin</mi><msup><mi>&phi;</mi><mo>`</mo></msup><mi>sin</mi><mi></mi><mi>&kappa;</mi><mi>sin</mi><mi>A</mi></mtd><mtd><mi>sin</mi><msup><mi>&phi;</mi><mo>`</mo></msup><mi>cos</mi><mi>&kappa;</mi></mtd></mtr><mtr><mtd><mi>sin</mi><msup><mi>&phi;</mi><mo>`</mo></msup><mi>sin</mi><mi>A</mi><mo>-</mo><mi>cos</mi><msup><mi>&phi;</mi><mo>`</mo></msup><mi>sin</mi><mi></mi><mi>&kappa;</mi><mi>cos</mi><mi>A</mi></mtd><mtd><mo>-</mo><mi>sin</mi><msup><mi>&phi;</mi><mo>`</mo></msup><mi>cos</mi><mi>A</mi><mo>-</mo><mi>cos</mi><msup><mi>&phi;</mi><mo>`</mo></msup><mi>sin</mi><mi>&kappa;A</mi></mtd><mtd><mi>cos</mi><msup><mi>&phi;</mi><mo>`</mo></msup><mi>cos</mi><mi>&kappa;</mi></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths>其中,φ`=acsin(sinφ/cosκ),磁目标的坐标表示为(X<sub>a</sub>,Y<sub>a</sub>,Z<sub>a</sub>)<sup>T</sup>=R(Ψ)(X′<sub>a</sub>,Y′<sub>a</sub>,Z′<sub>a</sub>)<sup>T</sup>;设罗经与磁传感器坐标系角度差的校准数据为Ψ′,布放后罗经的角度数据为Ψ″,得到:R(Ψ)=R(Ψ″)R(Ψ′),因此,磁目标的定向数据<img file="FDA0000470627490000017.GIF" wi="62" he="73" />和大地坐标数据的<img file="FDA0000470627490000018.GIF" wi="54" he="73" />校准公式为:(cosα<sub>i</sub>,cosβ<sub>i</sub>,cosγ<sub>i</sub>)<sup>T</sup>=R(Ψ)(cosα′<sub>i</sub>,cosβ′<sub>i</sub>,cosγ′<sub>i</sub>)<sup>T</sup>;b磁传感器阵列阵型的校准压力传感器提供每个磁传感器的深度数据或相对深度数据,得到磁传感器的zi,进一步对磁传感器的大地平面坐标(x<sub>i</sub>,y<sub>i</sub>)进行估计;采用水面目标舰船轨迹航行的方法进行阵列校准,舰船装有卫星定位系统GPS实时测量舰船的大地坐标,记GPS坐标数据为<img file="FDA0000470627490000019.GIF" wi="384" he="90" />同时装有同步芯片,使用通信系统发送坐标数据至基站;通过采用舰船盲航行获得各个阵元的初始值,即当第i个磁传感器定向为正上方时,记录此时的GPS平面坐标为第i个磁传感器的平面坐标,作为校准的初始值;获得初始值后,舰船按照一定的轨迹航行,此时包含测量数据和待估量的观测方程为:<maths num="0002"><![CDATA[<math><mrow><mfrac><msub><mrow><mi>cos</mi><mi>&beta;</mi></mrow><mn>1</mn></msub><msub><mrow><mi>cos</mi><mi>&alpha;</mi></mrow><mn>1</mn></msub></mfrac><mo>=</mo><mfrac><mrow><msub><mi>y</mi><mi>G</mi></msub><mo>-</mo><msub><mi>y</mi><mn>1</mn></msub></mrow><mrow><msub><mi>x</mi><mi>G</mi></msub><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub></mrow></mfrac><mo>+</mo><mi>V</mi></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mfrac><msub><mrow><mi>cos</mi><mi>&beta;</mi></mrow><mn>2</mn></msub><msub><mrow><mi>cos</mi><mi>&alpha;</mi></mrow><mn>2</mn></msub></mfrac><mo>=</mo><mfrac><mrow><msub><mi>y</mi><mi>G</mi></msub><mo>-</mo><msub><mi>y</mi><mn>2</mn></msub></mrow><mrow><msub><mi>x</mi><mi>G</mi></msub><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub></mrow></mfrac><mo>+</mo><mi>V</mi><mo>,</mo></mrow></math>]]></maths>...<maths num="0004"><![CDATA[<math><mrow><mfrac><msub><mrow><mi>cos</mi><mi>&beta;</mi></mrow><mi>N</mi></msub><msub><mrow><mi>cos</mi><mi>&alpha;</mi></mrow><mi>N</mi></msub></mfrac><mo>=</mo><mfrac><mrow><msub><mi>y</mi><mi>G</mi></msub><mo>-</mo><msub><mi>y</mi><mi>N</mi></msub></mrow><mrow><msub><mi>x</mi><mi>G</mi></msub><mo>-</mo><msub><mi>x</mi><mi>N</mi></msub></mrow></mfrac><mo>+</mo><mi>V</mi></mrow></math>]]></maths>以及<maths num="0005"><![CDATA[<math><mrow><msub><mrow><mi>cot</mi><mi>&gamma;</mi></mrow><mn>1</mn></msub><mo>=</mo><mfrac><msub><mi>z</mi><mi>i</mi></msub><msqrt><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>G</mi></msub><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>G</mi></msub><mo>-</mo><msub><mi>y</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac><mo>+</mo><mi>V</mi></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mrow><mi>cot</mi><mi>&gamma;</mi></mrow><mn>2</mn></msub><mo>=</mo><mfrac><msub><mi>z</mi><mi>i</mi></msub><msqrt><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>G</mi></msub><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>G</mi></msub><mo>-</mo><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac><mo>+</mo><mi>V</mi><mo>,</mo></mrow></math>]]></maths>...<maths num="0007"><![CDATA[<math><mrow><msub><mrow><mi>cot</mi><mi>&gamma;</mi></mrow><mi>N</mi></msub><mo>=</mo><mfrac><msub><mi>z</mi><mi>i</mi></msub><msqrt><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>G</mi></msub><mo>-</mo><msub><mi>x</mi><mi>N</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>G</mi></msub><mo>-</mo><msub><mi>y</mi><mi>N</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac><mo>+</mo><mi>V</mi></mrow></math>]]></maths>其中,V为系统的观测残差,使用最小二乘准则估计各个磁传感器的位置,由于方程带有一定的非线性,因此需采用非线性最小二乘方法,利用初始值进行估计,得第i个磁传感器的平面坐标(x<sub>i</sub>,y<sub>i</sub>)。
地址 266555 山东省青岛市经济技术开发区香江路98号