发明名称 双星敏感器联合确定载体姿态方法
摘要 本发明涉及一种双星敏感器联合确定船体姿态的方法,属于天文导航领域。方法包括以下步骤:一是星敏感器基于恒星参考矢量的蒙气差修正方法,即通过修正恒星参考矢量实现蒙气差的修正;二是基于联合双星敏感器恒星观测矢量的船体姿态确定方法,即将双星敏感器两个视场中每颗可识别恒星的观测矢量联合为观测矢量矩阵,最终与对应的恒星参考矢量矩阵完成船体姿态解算。本发明在简化星敏感器蒙气差修正算法的基础上,同时提高双星敏感器获取船体姿态的能力。主要内容涉及基于恒星参考矢量的蒙气差修正方法和基于联合双星敏感器恒星观测矢量的船体姿态确定方法两个方面。
申请公布号 CN104406583B 申请公布日期 2017.02.22
申请号 CN201410728741.1 申请日期 2014.12.05
申请人 中国人民解放军63680部队 发明人 潘良;朱伟康;刘冰;张同双;周海渊;赵李健;茅永兴;钟德安;焦宏伟
分类号 G01C21/02(2006.01)I 主分类号 G01C21/02(2006.01)I
代理机构 江阴市同盛专利事务所(普通合伙) 32210 代理人 唐纫兰;沈国安
主权项 一种双星敏感器联合确定载体姿态方法,其特征在于所述方法包括以下步骤:步骤一、基于恒星参考矢量的蒙气差修正1、将恒星S<sub>i</sub>在J2000.0坐标系的单位参考矢量V<sub>i</sub>,转换到船体地平坐标系,转换后的单位矢量记为V<sub>DPi</sub>,V<sub>i</sub>转换到V<sub>DPi</sub>的转换关系为:<img file="514843dest_path_image001.GIF" wi="510" he="79" />式(1)中,La、Lo分别为观测点的地理纬度和经度,x<sub>p</sub>、y<sub>p</sub>分别为极移参数,GAST为格林尼治真恒星时,ε<sub>A</sub>为平黄赤交角,ΔΨ、Δε分别为黄经章动和交角章动,Z<sub>A</sub>、θ<sub>A</sub>和ξ<sub>A</sub>为赤道岁差参数;2、计算恒星的地平方位角A<sub>DPi</sub>和地平仰角E<sub>DPi</sub>根据式(1)计算得到V<sub>DPi</sub>表示为[x<sub>DPi</sub> y<sub>DPi</sub> z<sub>DPi</sub>]<sup>T</sup>,V<sub>DPi</sub>的地平方位角A<sub>DPi</sub>、E<sub>DPi</sub>与x<sub>DPi</sub>、y<sub>DPi</sub>、z<sub>DPi</sub>有以下关系:<img file="743568dest_path_image002.GIF" wi="202" he="79" />3、根据地平仰角对V<sub>DPi</sub>进行蒙气差修正,得到蒙气差修正后的参考单位矢量V'<sub>DPi</sub><img file="523305dest_path_image003.GIF" wi="270" he="87" />其中:<img file="918514dest_path_image004.GIF" wi="122" he="35" />,ρ即为蒙气差修正量,<img file="795204dest_path_image005.GIF" wi="144" he="30" />,其中A、B的计算公式分别为:<img file="617666dest_path_image006.GIF" wi="128" he="91" />式中,t温度为(℃)、P为气压(hpa),ρ<sub>0</sub>计算公式为:<img file="697749dest_path_image007.GIF" wi="470" he="51" />当E<sub>DPi</sub>&gt;45°时,a为1;当14°≤E<sub>DPi</sub>≤45°时,a计算公式如下为:<img file="896649dest_path_image008.GIF" wi="458" he="58" />;步骤二、基于联合双星敏感器观测矢量的船体姿态确定1、采用如下公式将星敏感器观测矢量转换到船体甲板系<img file="565528dest_path_image009.GIF" wi="109" he="67" />(4)其中,R<sub>S1</sub><sup>b</sup>、R<sub>S2</sub><sup>b</sup>分别为星敏感器1、星敏感器2的安装矩阵,W<sub>1</sub>、W<sub>2</sub>分别为星敏感器1、星敏感器2的观测矢量矩阵;2、将星敏感器参考矢量转换到船体地平系将星敏感器1、星敏感器2识别恒星的参考矢量矩阵V<sub>1</sub>、V<sub>2</sub>,转换到船体地平系并分别经公式(1)、(2)、(3)的蒙气差修正重组后得到参考矢量分别为V<sub>1</sub>'、V<sub>2</sub>';3、组成联合观测和参考矢量矩阵解算船体姿态角将星敏感器1、星敏感器2相对船体甲板系的观测矢量矩阵W<sub>b1</sub>、W<sub>b2</sub>,组成联合观测矢量矩阵W<sub>b</sub>,将V<sub>1</sub>'、V<sub>2</sub>'在船体地平系下的参考矢量矩阵,组成联合参考矢量矩阵V<sub>DP</sub>',根据式(5)得到船体甲板系到船体地平系的姿态转换矩阵,最后根据式(6)得到船体姿态角,<img file="621208dest_path_image011.GIF" wi="482" he="190" />其中,<img file="375538dest_path_image012.GIF" wi="499" he="47" />为姿态矩阵R<sup>b</sup><sub>DP</sub>中的第i行、第j列元素,K为航向角、θ为纵摇角、Ψ为横摇角,航向K需要进行象限判断。
地址 214400 江苏省无锡市江阴市103信箱850号