发明名称 一种基于CCD星敏感器的标定方法
摘要 本发明提供的是一种基于CCD星敏感器的标定方法。(1)采集CCD星敏感器的输出:CCD星敏感器的坐标系相对于i系之间的姿态信息(2)采集当地位置信息,得到地球坐标系e系相对于导航坐标系n系的转换矩阵(3)求解e系相对于i系之间的转换矩阵(4)通过(1)、(2)、(3)步骤中所给出的信息,解算得到姿态矩阵;(5)将步骤(4)中得到的姿态矩阵经过换算得到失准角,将其作为观测方程,代入卡尔曼滤波器进行滤波估计;(6)通过步骤(5)估计出陀螺的常值漂移和加速度计零偏。本发明的方法,短时间内可以达到稳定的标定结果。不需要进行任何机动措施,便可以估计出陀螺常值漂移和加速度计零偏。
申请公布号 CN101943585B 申请公布日期 2013.07.31
申请号 CN201010215400.6 申请日期 2010.07.02
申请人 哈尔滨工程大学 发明人 高伟;付建楠;张鑫;奔粤阳;徐博;周广涛;于强;张永刚;吴晓;胡丹
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 代理人
主权项 1.一种基于CCD星敏感器的标定方法,其特征是:(1)采集CCD星敏感器的输出:CCD星敏感器的坐标系相对于惯性坐标系即i系之间的姿态信息<img file="FSB00001071923300011.GIF" wi="77" he="55" />所述采集CCD星敏感器的输出的方法为:i系与船舶载体坐标系即b系之间的转换矩阵为<img file="FSB00001071923300012.GIF" wi="249" he="63" />其中:<img file="FSB00001071923300013.GIF" wi="53" he="63" />为CCD星敏感器坐标系即s系与b系之间的转换矩阵,在导航设备装船时通过光学瞄准精确获得;将天球坐标系O-UVW按照先绕W轴逆时针转ω角,得到O-U<sub>1</sub>V<sub>1</sub>W<sub>1</sub>坐标系,再绕U<sub>1</sub>逆时针转u角,使W<sub>1</sub>轴与Z<sub>s</sub>重合,得到O-U<sub>2</sub>V<sub>2</sub>W<sub>2</sub>坐标系,最后再绕W<sub>2</sub>轴逆时针旋转v角,得到O<sub>s</sub>-U<sub>s</sub>V<sub>s</sub>W<sub>s</sub>坐标系;<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>i</mi><mi>s</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mi></mi><mi>&omega;</mi><mi>cos</mi><mi>v</mi><mo>-</mo><mi>sin</mi><mi></mi><mi>&omega;</mi><mi>sin</mi><mi>v</mi><mi>cos</mi><mi>u</mi></mtd><mtd><mi>sin</mi><mi></mi><mi>&omega;</mi><mi>cos</mi><mi>v</mi><mo>+</mo><mi>cos</mi><mi></mi><mi>&omega;</mi><mi>sin</mi><mi>v</mi><mi>cos</mi><mi>u</mi></mtd><mtd><mi>sin</mi><mi>v</mi><mi>sin</mi><mi>u</mi></mtd></mtr><mtr><mtd><mo>-</mo><mi>cos</mi><mi></mi><mi>&omega;</mi><mi>sin</mi><mi>v</mi><mo>-</mo><mi>sin</mi><mi></mi><mi>&omega;</mi><mi>cos</mi><mi>v</mi><mi>cos</mi><mi>u</mi></mtd><mtd><mo>-</mo><mi>sin</mi><mi></mi><mi>&omega;</mi><mi>sin</mi><mi>v</mi><mo>+</mo><mi>cos</mi><mi></mi><mi>&omega;</mi><mi>cos</mi><mi>v</mi><mi>cos</mi><mi>u</mi></mtd><mtd><mi>cos</mi><mi>v</mi><mi>sin</mi><mi>u</mi></mtd></mtr><mtr><mtd><mi>sin</mi><mi></mi><mi>&omega;</mi><mi>sin</mi><mi>u</mi></mtd><mtd><mo>-</mo><mi>cos</mi><mi></mi><mi>&omega;</mi><mi>sin</mi><mi>u</mi></mtd><mtd><mi>cos</mi><mi>u</mi></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>(2)采集当地位置信息即经度和纬度,得到地球坐标系即e系相对于导航坐标系即n系的转换矩阵<img file="FSB00001071923300015.GIF" wi="78" he="56" /><img file="FSB00001071923300016.GIF" wi="845" he="258" />(3)求解e系相对于i系之间的转换矩阵<img file="FSB00001071923300017.GIF" wi="78" he="55" /><maths num="0002"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>i</mi><mi>e</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mrow><mo>(</mo><msub><mi>A</mi><mi>j</mi></msub><mo>+</mo><msub><mi>&omega;</mi><mi>ie</mi></msub><mo>&CenterDot;</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><mi>sin</mi><mrow><mo>(</mo><msub><mi>A</mi><mi>j</mi></msub><mo>+</mo><msub><mi>&omega;</mi><mi>ie</mi></msub><mo>&CenterDot;</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>-</mo><mi>sin</mi><mrow><mo>(</mo><msub><mi>A</mi><mi>j</mi></msub><mo>+</mo><msub><mi>&omega;</mi><mi>ie</mi></msub><mo>&CenterDot;</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><mi>cos</mi><mrow><mo>(</mo><msub><mi>A</mi><mi>j</mi></msub><mo>+</mo><msub><mi>&omega;</mi><mi>ie</mi></msub><mo>&CenterDot;</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>ω<sub>ie</sub>为地球自转角速度,t是世界标准时间系统提供的具体时间,A<sub>j</sub>是初始位置与春分点之间的夹角;(4)通过(1)、(2)、(3)步骤中所给出的信息,解算得到姿态矩阵,所述姿态矩阵<img file="FSB00001071923300019.GIF" wi="52" he="66" />为:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>n</mi><mi>b</mi></msubsup><mo>=</mo><msubsup><mi>C</mi><mi>i</mi><mi>b</mi></msubsup><mo>&CenterDot;</mo><mo>(</mo><msubsup><mi>C</mi><mi>e</mi><mi>n</mi></msubsup><mo>&CenterDot;</mo><msubsup><mi>C</mi><mi>i</mi><mi>e</mi></msubsup><msup><mo>)</mo><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>;</mo></mrow></math>]]></maths>(5)将步骤(4)中得到的姿态矩阵经过换算得到失准角,将其作为观测方程,代入卡尔曼滤波器进行滤波估计;具体方法为:<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>n</mi><msup><mi>n</mi><mo>&prime;</mo></msup></msubsup><mo>=</mo><msubsup><mi>C</mi><mi>b</mi><msup><mi>n</mi><mo>&prime;</mo></msup></msubsup><mo>&CenterDot;</mo><msubsup><mi>C</mi><mi>i</mi><mi>b</mi></msubsup><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><msubsup><mi>C</mi><mi>e</mi><mi>n</mi></msubsup><mo>&CenterDot;</mo><msubsup><mi>C</mi><mi>i</mi><mi>e</mi></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup></mrow></math>]]></maths>φ<sub>x</sub>=C(2,3)φ<sub>y</sub>=C(3,1)φ<sub>z</sub>=C(1,2)其中C表示<img file="FSB00001071923300023.GIF" wi="111" he="70" />量测量为:<maths num="0005"><![CDATA[<math><mrow><mi>Z</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&delta;V</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&phi;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中:δV(t)为惯导解算的速度与真实的速度之间的差值;(6)通过步骤(5)估计出陀螺的常值漂移和加速度计零偏。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室