CN107228683B - 一种多星敏感器间慢变误差实时在轨修正方法 - Google Patents
一种多星敏感器间慢变误差实时在轨修正方法 Download PDFInfo
- Publication number
- CN107228683B CN107228683B CN201710501158.0A CN201710501158A CN107228683B CN 107228683 B CN107228683 B CN 107228683B CN 201710501158 A CN201710501158 A CN 201710501158A CN 107228683 B CN107228683 B CN 107228683B
- Authority
- CN
- China
- Prior art keywords
- head
- quaternion
- heads
- star sensor
- star
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Manufacturing & Machinery (AREA)
- Automation & Control Theory (AREA)
- Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
一种多星敏感器间慢变误差实时在轨修正方法,修正星敏感器各头部之间的慢变误差,包含如下步骤:S1、计算星敏感器各头部两两之间的坐标变换四元数滤波实测值QRotOHiOHjCAL;S2、由地面控制中心指定星敏感器上某一头部为基准头部OHref;S3、计算基准头部OHref到星敏感器各头部OHi的坐标变换四元数地面精测值QOHiOHref;S4、计算星敏感器各头部OHi到基准头部OHref的坐标变换四元数实测值QRotOHiOHref;S5、根据S3和S4得到的QOHiOHref和QRotOHiOHref计算修正头部间慢变误差后的各头部姿态四元数qOHiRc。本发明能在轨实时估计并修正,通过简单的计算方式提高多星敏感器联合定姿的姿态确定精度,并且简单有效的修正非基准头部与基准头部之间的慢变误差,并获得修正头部间慢变误差后的头部姿态四元数。
Description
技术领域
本发明涉及卫星星敏感器姿态修正技术领域,具体涉及一种多星敏感器间慢变误差实时在轨修正方法。
背景技术
随着航天科技与应用领域的不断发展与开拓,卫星的研制与实现技术呈现出两种趋势:一是以高精度遥感卫星为代表,任务要求决定了卫星星体的结构组成复杂、系统控制精度指标极高;另一个方向是朝着小型、快捷、低成本、结构与配置简化方向发展,实现单一任务,或者通过发射较多数量实现多任务或某项复杂工作。针对前一种趋势的复杂卫星,控制系统通常配置星敏感器作为姿态敏感器部件。
星敏感器是目前卫星上使用的精度最高的姿态测量敏感器,它通过测量恒星的观测矢量在卫星坐标系中的方位以及恒星亮度,再利用星历表得到这些恒星在惯性坐标系中的方位,经姿态确定算法计算即可提供卫星在惯性系中的姿态信息,精度可达到角秒级。但由于星敏感器光轴方向的测量噪声较大,单星敏感器姿态确定精度较低,不能满足高精度姿态确定要求。因此,目前主流的高精度控制系统都采用多星敏感器联合定姿技术来提高姿态确定的精度,利用一个星敏感器垂直光轴方向的姿态信息来补偿另一个星敏感器光轴方向姿态测量精度的不足。
卫星在轨运行过程中受太阳照射角度影响呈现温度周期性慢变,因此,星敏感器安装结构会受其影响而产生形变,导致星敏感器测量输出相对本体坐标系基准发生动态偏离。此外,星敏感器观测恒星过程中,其光轴所指向的天区随卫星轨道运动发生周期性变化,导致光学系统误差、标定误差和星表误差随恒星进出视场而变,也会造成周期性慢变误差,这些变化对星敏感器基准头部及非基准头部的姿态四元数均产生影响,进而对多星敏感器联合定姿的姿态确定精度造成影响。因此,为了获取高精度的可信的卫星姿态,需要对多星敏感器间慢变误差的实时估计及修正方法进行研究。
发明内容
本发明的目的在于提供一种多星敏感器间慢变误差实时在轨修正方法,应用于卫星上多个星敏感器之间修正慢变误差。
一种多星敏感器间慢变误差实时在轨修正方法,修正星敏感器各头部之间的慢变误差,包含如下步骤:
S1、计算星敏感器各头部两两之间的坐标变换四元数滤波实测值QRotOHiOHjCAL;
S2、由地面控制中心指定星敏感器上某一头部为基准头部OHref;
S3、计算基准头部OHref到星敏感器各头部OHi的坐标变换四元数地面精测值QOHiOHref;
S4、计算星敏感器各头部OHi到基准头部OHref的坐标变换四元数实测值QRotOHiOHref;
S5、根据S3和S4得到的QOHiOHref和QRotOHiOHref计算修正头部间慢变误差后的各头部姿态四元数qOHiRc。
上述的一种多星敏感器间慢变误差实时在轨修正方法,其中,所述步骤S1包含:
S11、计算任意星敏感器头部OHi到另一任意头部OHj的坐标变换四元数Q’RotOHiOHjCAL,其中,i=1,2,3……N-1,j=i+1,i+2,i+3……N,N为头部的个数;
S12、对S1得到的Q’RotOHiOHjCAL进行滤波,得到任意头部OHi到另一任意头部OHj的坐标变换四元数滤波实测值QRotOHiOHjCAL,其初值为地面安装精测值;
S13、重复步骤S11和S12直至遍历所有头部。
上述的一种多星敏感器间慢变误差实时在轨修正方法,其中,所述步骤S12中的滤波算法如下:
其中T为控制系统的控制周期,L为滤波时间常数,k为滤波算法计算时刻。
上述的一种多星敏感器间慢变误差实时在轨修正方法,其中,所述步骤S3中,
其中,QbOHref表示基准头部OHref姿态测量坐标系到卫星本体坐标系的转换四元数。
上述的一种多星敏感器间慢变误差实时在轨修正方法,其中,所述步骤S4中,
当i<ref时,该头部与基准头部的坐标变换四元数实测值为:
QRotOHiOHref=QRotOHiOHjCAL
当i>ref时,该头部与基准头部的坐标变换四元数实测值为:
本发明的优点和有益效果是:
(1)本发明可由星上软件运行,能在轨实时估计并修正,通过简单的计算方式,即可提高多星敏感器联合定姿的姿态确定精度。
(2)本发明简单有效的修正非基准头部与基准头部之间的慢变误差,并获得修正头部间慢变误差后的头部姿态四元数。
附图说明
图1是本发明的流程示意图。
图2是本发明中步骤S1的细化流程示意图。
具体实施方式
以下结合附图,通过详细说明一个较佳的具体实施例,对本发明做进一步阐述。
如图1所示,是一种多星敏感器间慢变误差实时在轨修正方法,包含如下步骤:
S1、计算星敏感器各头部两两之间的坐标变换四元数滤波实测值QRotOHiOHjCAL;
S2、由地面控制中心指定星敏感器上某一头部OHj为基准头部OHref;
S3、计算基准头部OHref到星敏感器各头部OHi的坐标变换四元数地面精测值QOHiOHref:
其中,QbOHref表示基准头部OHref姿态测量坐标系到卫星本体坐标系的转换四元数。若地面控制中心指定以OH1为基准头部,则QbOHref=QbOH1;若地面控制中心指定以OH2为基准头部,则QbOHref=QbOH2;若地面控制中心指定以OHj为基准头部,则QbOHref=QbOHj,各头部与本地系的转换四元数均由地面进行精测得到。
S4、计算星敏感器各头部OHi到基准头部OHref的坐标变换四元数实测值QRotOHiOHref:
当i<ref时,该头部与基准头部的坐标变换四元数实测值为:
QRotOHiOHref=QRotOHiOHjCAL
当i>ref时,该头部与基准头部的坐标变换四元数实测值为:
S5、根据S3和S4得到的QOHiOHref和QRotOHiOHref计算修正头部间慢变误差后的各头部姿态四元数qOHiRc:
其中,QRotOHiOHref为星敏感器头部OHi测量坐标系转换到基准头部OHref测量坐标系的四元数实际值;QOHiOHref为基准头部OHref测量坐标系转换到星敏感器头部OHi测量坐标系的安装四元数精测值。
如图2所示,所述步骤S1具体包含:
S11、计算任意星敏感器头部OHi到另一任意头部OHj的坐标变换四元数Q’RotOHiOHjCAL:
其中,i=1,2,3……N-1,j=i+1,i+2,i+3……N,N为头部的个数,QOHiRc为任意头部OHi输出的姿态四元数,QOHjRc为另一任意头部OHj输出的姿态四元数。;
S12、对S1得到的Q’RotOHiOHjCAL进行滤波,得到任意头部OHi到另一任意头部OHj的坐标变换四元数滤波实测值QRotOHiOHjCAL,其初值为地面安装精测值;
S13、重复步骤S11和S12直至遍历所有头部。
其中,所述步骤S12中的滤波算法如下:
式中,T为控制系统的控制周期,L为滤波时间常数,k为滤波算法计算时刻。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (2)
1.一种多星敏感器间慢变误差实时在轨修正方法,修正星敏感器各头部之间的慢变误差,其特征在于,包含如下步骤:
S1、计算星敏感器各头部两两之间的坐标变换四元数滤波实测值QRotOHiOHjCAL;
S2、由地面控制中心指定星敏感器上的其中一头部为基准头部OHref;
S3、计算基准头部OHref到星敏感器各头部OHi的坐标变换四元数地面精测值QOHiOHref;
其中,QbOHref表示基准头部OHref姿态测量坐标系到卫星本体坐标系的转换四元数;
S4、计算星敏感器各头部OHi到基准头部OHref的坐标变换四元数实测值QRotOHiOHref;
当i<ref时,该头部与基准头部的坐标变换四元数实测值为:
QRotOHiOHref=QRotOHiOHjCAL
当i>ref时,该头部与基准头部的坐标变换四元数实测值为:
S5、根据S3和S4得到的QOHiOHref和QRotOHiOHref计算修正头部间慢变误差后的各头部姿态四元数qOHiRc;
其中,S1进一步包含:
S11、计算任意星敏感器头部OHi到另一任意头部OHj的坐标变换四元数Q′RotOHiOHjCAL,其中i=1,2,3……N-1,j=i+1,i+2,i+3……N,N为头部的个数;
S12、对S1得到的Q′RotOHiOHjCAL进行滤波,得到任意头部OHi到另一任意头部OHj的坐标变换四元数滤波实测值QRotOHiOHjCAL;
S13、重复步骤S11和S12直至遍历所有头部。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710501158.0A CN107228683B (zh) | 2017-06-27 | 2017-06-27 | 一种多星敏感器间慢变误差实时在轨修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710501158.0A CN107228683B (zh) | 2017-06-27 | 2017-06-27 | 一种多星敏感器间慢变误差实时在轨修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107228683A CN107228683A (zh) | 2017-10-03 |
CN107228683B true CN107228683B (zh) | 2020-04-10 |
Family
ID=59936491
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710501158.0A Active CN107228683B (zh) | 2017-06-27 | 2017-06-27 | 一种多星敏感器间慢变误差实时在轨修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107228683B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109649692B (zh) * | 2018-12-28 | 2021-10-01 | 上海航天控制技术研究所 | 一种基于推力器偏差估计的轨控策略在轨修正方法和系统 |
CN109858151B (zh) * | 2019-01-31 | 2023-04-18 | 北京空间飞行器总体设计部 | 一种适用于惯性空间观测卫星的热环境保障设计方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102735260A (zh) * | 2012-06-18 | 2012-10-17 | 航天东方红卫星有限公司 | 一种星敏感器在轨测量误差的确定方法 |
WO2014086340A1 (de) * | 2012-12-04 | 2014-06-12 | Jena Optronik Gmbh | Verfahren zur automatischen korrektur von ausrichtungsfehlern in sternsensorsystemen |
CN104792340A (zh) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法 |
CN105277195A (zh) * | 2015-11-04 | 2016-01-27 | 上海新跃仪表厂 | 一种星敏感器单机间相对安装误差在轨辨识方法 |
CN105318871A (zh) * | 2015-11-09 | 2016-02-10 | 中国人民解放军63680部队 | 双星敏感器载体安装矩阵动态标定方法 |
CN105698764A (zh) * | 2016-01-30 | 2016-06-22 | 武汉大学 | 一种光学遥感卫星影像时变系统误差建模补偿方法及系统 |
CN106338296A (zh) * | 2016-11-04 | 2017-01-18 | 上海航天控制技术研究所 | 一种双星敏感器在轨实时交互的修正方法 |
CN106441283A (zh) * | 2016-10-12 | 2017-02-22 | 上海航天控制技术研究所 | 一种自主标较的星矢量级多视场星敏感器数据融合方法 |
-
2017
- 2017-06-27 CN CN201710501158.0A patent/CN107228683B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102735260A (zh) * | 2012-06-18 | 2012-10-17 | 航天东方红卫星有限公司 | 一种星敏感器在轨测量误差的确定方法 |
WO2014086340A1 (de) * | 2012-12-04 | 2014-06-12 | Jena Optronik Gmbh | Verfahren zur automatischen korrektur von ausrichtungsfehlern in sternsensorsystemen |
CN104792340A (zh) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法 |
CN105277195A (zh) * | 2015-11-04 | 2016-01-27 | 上海新跃仪表厂 | 一种星敏感器单机间相对安装误差在轨辨识方法 |
CN105318871A (zh) * | 2015-11-09 | 2016-02-10 | 中国人民解放军63680部队 | 双星敏感器载体安装矩阵动态标定方法 |
CN105698764A (zh) * | 2016-01-30 | 2016-06-22 | 武汉大学 | 一种光学遥感卫星影像时变系统误差建模补偿方法及系统 |
CN106441283A (zh) * | 2016-10-12 | 2017-02-22 | 上海航天控制技术研究所 | 一种自主标较的星矢量级多视场星敏感器数据融合方法 |
CN106338296A (zh) * | 2016-11-04 | 2017-01-18 | 上海航天控制技术研究所 | 一种双星敏感器在轨实时交互的修正方法 |
Non-Patent Citations (2)
Title |
---|
四元数理论及其在坐标转换中的应用;张捍卫;《大地测量与地球动力学》;20151031;第35卷(第5期);第807-810页 * |
星敏感器慢变误差校准方法研究;徐樱;《上海航天》;20160921;第33卷(第4期);第63-69页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107228683A (zh) | 2017-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108827310B (zh) | 一种船用星敏感器辅助陀螺仪在线标定方法 | |
CN103994763B (zh) | 一种火星车的sins/cns深组合导航系统及其实现方法 | |
CN102506893B (zh) | 一种基于地标信息的星敏感器低频误差补偿方法 | |
CN110926468B (zh) | 基于传递对准的动中通天线多平台航姿确定方法 | |
CN111323050B (zh) | 一种捷联惯导和多普勒组合系统标定方法 | |
CN105806369B (zh) | 一种星敏感器在轨光行差修正方法 | |
CN107747953B (zh) | 一种多敏感器数据与轨道信息时间同步方法 | |
CN112504298B (zh) | 一种gnss辅助的dvl误差标定方法 | |
CN110285815B (zh) | 一种可在轨全程应用的微纳卫星多源信息姿态确定方法 | |
US20180346154A1 (en) | Method and Apparatus for Spacecraft Gyroscope Scale Factor Calibration | |
JP2008216062A (ja) | 移動体姿勢計測装置 | |
CN110228605B (zh) | 一种基于太阳敏感器的卫星安全对日控制方法 | |
CN110285834B (zh) | 基于一点位置信息的双惯导系统快速自主重调方法 | |
CN111189474A (zh) | 基于mems的marg传感器的自主校准方法 | |
CN109489661B (zh) | 一种卫星初始入轨时陀螺组合常值漂移估计方法 | |
CN107228683B (zh) | 一种多星敏感器间慢变误差实时在轨修正方法 | |
CN108072387B (zh) | 一种低精度敏感器在轨偏差的地面校正方法及系统 | |
CN110388942B (zh) | 一种基于角度和速度增量的车载姿态精对准系统 | |
CN115754861A (zh) | 一种卫星磁强计在轨标定方法 | |
CN116448145A (zh) | 一种基于偏振矢量空间差分的航姿确定方法 | |
CN113447043B (zh) | 一种基于gnss的卫星天文导航系统误差自主标定方法及系统 | |
CN204255368U (zh) | 一种适用于火星车的sins/cns深组合导航系统 | |
CN109655080B (zh) | 一种数字式太阳敏感器在轨标定方法 | |
CN110160530B (zh) | 一种基于四元数的航天器姿态滤波方法 | |
CN109506645B (zh) | 一种星敏感器安装矩阵地面精测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |