发明名称 一种飞行器星敏感器与捷联惯组安装误差的测量方法
摘要 本发明公开了一种飞行器星敏感器与捷联惯组安装误差的测量方法,方法包括以下步骤:S1、建立中间转换坐标系O<sub>m</sub>X<sub>m</sub>Y<sub>m</sub>Z<sub>m</sub>,其原点与惯组坐标系原点重合,O<sub>m</sub>X<sub>m</sub>Y<sub>m</sub>平面位于当地水平面内,O<sub>m</sub>X<sub>m</sub>轴方向沿飞行器纵向方向,O<sub>m</sub>Y<sub>m</sub>轴方向与捷联惯组棱镜准直;S2、计算捷联惯组在中间转换坐标系下的姿态矩阵<img file="DDA0000850411250000011.GIF" wi="55" he="67" />;S3、计算星敏感器测量坐标系在中间转换坐标系中的姿态矩阵<img file="DDA0000850411250000012.GIF" wi="55" he="64" />;S4、解算星敏感器测量坐标系在捷联惯组坐标系下的转换矩阵,确定星敏感器测量坐标系上各轴与捷联惯组坐标系各轴的投影关系,获得星敏感器与捷联惯组的安装误差。本发明方便快捷,准确度高,可在室内进行快速的星敏感器与惯组安装误差测量。
申请公布号 CN105424060A 申请公布日期 2016.03.23
申请号 CN201510794566.0 申请日期 2015.11.18
申请人 湖北航天技术研究院总体设计所 发明人 王亚军;周海;李春权;王勇;陆壮志
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 武汉东喻专利代理事务所(普通合伙) 42224 代理人 宋业斌
主权项 一种飞行器星敏感器与捷联惯组安装误差的测量方法,其特征在于,所述方法包括以下步骤:S1、建立中间转换坐标系O<sub>m</sub>X<sub>m</sub>Y<sub>m</sub>Z<sub>m</sub>,所述中间转换坐标系O<sub>m</sub>X<sub>m</sub>Y<sub>m</sub>Z<sub>m</sub>原点O<sub>m</sub>与惯组坐标系O<sub>t</sub>X<sub>t</sub>Y<sub>t</sub>Z<sub>t</sub>原点O<sub>t</sub>重合,O<sub>m</sub>X<sub>m</sub>Y<sub>m</sub>平面位于当地水平面内,O<sub>m</sub>X<sub>m</sub>轴方向沿飞行器纵向方向,O<sub>m</sub>Y<sub>m</sub>轴方向通过第一经纬仪与捷联惯组棱镜准直,根据右手法则确定O<sub>m</sub>Z<sub>m</sub>轴方向;S2、计算捷联惯组在中间转换坐标系下O<sub>1</sub>的姿态矩阵<img file="FDA0000850411220000011.GIF" wi="87" he="85" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msubsup><mi>C</mi><mi>b</mi><mi>t</mi></msubsup><mo>=</mo></mrow></mtd></mtr><mtr><mtd><mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>c</mi><mi>o</mi><mi>s</mi><mi>&psi;</mi><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi><mo>+</mo><mi>sin</mi><mi>&psi;</mi><mi>sin</mi><mi>&theta;</mi><mi>sin</mi><mi>&phi;</mi></mrow></mtd><mtd><mrow><mo>-</mo><mi>c</mi><mi>o</mi><mi>s</mi><mi>&phi;</mi><mi>sin</mi><mi>&psi;</mi></mrow></mtd><mtd><mrow><mi>sin</mi><mi>&theta;</mi><mi>cos</mi><mi>&psi;</mi><mo>-</mo><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi><mi>sin</mi><mi>&psi;</mi><mi>sin</mi><mi>&phi;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi><mi>sin</mi><mi>&psi;</mi><mo>-</mo><mi>sin</mi><mi>&phi;</mi><mi>sin</mi><mi>&theta;</mi><mi>cos</mi><mi>&psi;</mi></mrow></mtd><mtd><mrow><mi>cos</mi><mi>&psi;</mi><mi>cos</mi><mi>&phi;</mi></mrow></mtd><mtd><mrow><mi>sin</mi><mi>&theta;</mi><mi>sin</mi><mi>&psi;</mi><mo>+</mo><mi>c</mi><mi>o</mi><mi>s</mi><mi>&psi;</mi><mi>sin</mi><mi>&phi;</mi><mi>cos</mi><mi>&theta;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mi>sin</mi><mi>&theta;</mi><mi>cos</mi><mi>&phi;</mi></mrow></mtd><mtd><mrow><mo>-</mo><mi>sin</mi><mi>&phi;</mi></mrow></mtd><mtd><mrow><mi>c</mi><mi>o</mi><mi>s</mi><mi>&phi;</mi><mi>cos</mi><mi>&theta;</mi></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000850411220000012.GIF" wi="1814" he="318" /></maths>其中,θ,Φ,<img file="FDA0000850411220000013.GIF" wi="55" he="69" />分别为惯组坐标系下的俯仰角、滚动角和偏航角;偏航角<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&psi;</mi><mo>=</mo><mn>0</mn><mo>,</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>&theta;</mi><mo>=</mo><mi>arcsin</mi><mrow><mo>(</mo><mfrac><mrow><mo>-</mo><msubsup><mi>f</mi><mi>x</mi><mi>b</mi></msubsup></mrow><mi>g</mi></mfrac><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>&phi;</mi><mo>=</mo><mi>arcsin</mi><mrow><mo>(</mo><mfrac><mrow><mo>-</mo><msubsup><mi>f</mi><mi>y</mi><mi>b</mi></msubsup></mrow><mrow><mi>g</mi><mi>cos</mi><mi>&theta;</mi></mrow></mfrac><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000850411220000014.GIF" wi="791" he="271" /></maths>其中<img file="FDA0000850411220000015.GIF" wi="202" he="156" />为飞行器坐标系中的加速度矢量在x、y轴的分量,g为重力加速度;S3、采集第二经纬仪模拟的星矢量O<sub>m</sub>X<sub>m</sub>Y<sub>m</sub>Z<sub>m</sub>坐标系中的坐标值(X<sub>i</sub>,Y<sub>i</sub>,Z<sub>i</sub>),所述第二经纬仪的观测镜头在星敏感器的视场中央,并与第一经纬仪准直;采集CCD上对应像素点在星敏感器测量坐标系中的坐标值(u<sub>i</sub>,v<sub>i</sub>),根据<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>u</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>v</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>=</mo><mfrac><mn>1</mn><mrow><mi>S</mi><mi>x</mi></mrow></mfrac><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>p</mi><mn>11</mn></msub></mtd><mtd><msub><mi>p</mi><mn>12</mn></msub></mtd><mtd><msub><mi>p</mi><mn>13</mn></msub></mtd><mtd><msub><mi>p</mi><mn>14</mn></msub></mtd></mtr><mtr><mtd><msub><mi>p</mi><mn>21</mn></msub></mtd><mtd><msub><mi>p</mi><mn>22</mn></msub></mtd><mtd><msub><mi>p</mi><mn>23</mn></msub></mtd><mtd><msub><mi>p</mi><mn>24</mn></msub></mtd></mtr><mtr><mtd><msub><mi>p</mi><mn>31</mn></msub></mtd><mtd><msub><mi>p</mi><mn>32</mn></msub></mtd><mtd><msub><mi>p</mi><mn>33</mn></msub></mtd><mtd><msub><mi>p</mi><mn>34</mn></msub></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>X</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Y</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Z</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000850411220000016.GIF" wi="1028" he="270" /></maths>求解系数矩阵<img file="FDA0000850411220000021.GIF" wi="734" he="199" />利用系数矩阵p计算星敏感器测量坐标系在中间转换坐标系O<sub>m</sub>X<sub>m</sub>Y<sub>m</sub>Z<sub>m</sub>中的姿态矩阵<img file="FDA0000850411220000022.GIF" wi="95" he="79" /><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>C</mi><mi>s</mi><mi>t</mi></msubsup><mo>=</mo><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>r</mi><mn>1</mn></msub></mtd><mtd><msub><mi>r</mi><mn>2</mn></msub></mtd><mtd><msub><mi>r</mi><mn>3</mn></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000850411220000026.GIF" wi="529" he="70" /></maths>其中<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>r</mi><mn>1</mn><mo>=</mo><mo>|</mo><mi>p</mi><mn>3</mn><msup><mo>|</mo><mn>2</mn></msup><mo>&CenterDot;</mo><mi>p</mi><mn>1</mn><mo>-</mo><mi>p</mi><mn>1</mn><mo>&CenterDot;</mo><mi>p</mi><mn>3</mn><mo>&CenterDot;</mo><mi>p</mi><mn>3</mn><mo>/</mo><mo>|</mo><mi>p</mi><mn>1</mn><mo>&times;</mo><mi>p</mi><mn>3</mn><mo>|</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>r</mi><mn>2</mn><mo>=</mo><mo>|</mo><mi>p</mi><mn>3</mn><msup><mo>|</mo><mn>2</mn></msup><mo>&CenterDot;</mo><mi>p</mi><mn>2</mn><mo>-</mo><mi>p</mi><mn>2</mn><mo>&CenterDot;</mo><mi>p</mi><mn>3</mn><mo>&CenterDot;</mo><mi>p</mi><mn>3</mn><mo>/</mo><mo>|</mo><mi>p</mi><mn>2</mn><mo>&times;</mo><mi>p</mi><mn>3</mn><mo>|</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>r</mi><mn>3</mn><mo>=</mo><mi>p</mi><mn>3</mn></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000850411220000023.GIF" wi="1134" he="229" /></maths>p1、p2、p3为系数矩阵p中的元素组成的矩阵,p1=[p<sub>11</sub> p<sub>12</sub> p<sub>13</sub>],p2=[p<sub>21</sub> p<sub>22</sub> p<sub>23</sub>],p3=[p<sub>31</sub> p<sub>32</sub> p<sub>33</sub>],Sx为CCD横坐标和纵坐标之间的比例系数;S4、解算星敏感器测量坐标系在捷联惯组坐标下的转换矩阵<img file="FDA0000850411220000024.GIF" wi="374" he="86" />根据所述转换矩阵<img file="FDA0000850411220000025.GIF" wi="74" he="78" />确定星敏感器的测量坐标系上各轴与捷联惯组坐标系各轴的投影关系,获得星敏感器与捷联惯组的安装误差。
地址 430040 湖北省武汉市金山大道9号