CN107402023B - 地球静止轨道卫星星敏慢变误差预报校正的方法 - Google Patents
地球静止轨道卫星星敏慢变误差预报校正的方法 Download PDFInfo
- Publication number
- CN107402023B CN107402023B CN201710681151.1A CN201710681151A CN107402023B CN 107402023 B CN107402023 B CN 107402023B CN 201710681151 A CN201710681151 A CN 201710681151A CN 107402023 B CN107402023 B CN 107402023B
- Authority
- CN
- China
- Prior art keywords
- error
- star sensor
- measurement
- attitude
- satellite
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Navigation (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种地球静止轨道卫星星敏慢变误差预报校正的方法,其步骤主要为:一、根据载荷标定的姿态确定误差,计算星敏感器的姿态测量误差;二、根据星敏感器的安装矩阵,计算星敏感器测量系下的测量误差;三、将测量误差数据通过分段多阶傅立叶函数来拟合星敏感器的慢变误差;四、根据地方时与恒星时角的转换关系,计算预报日期与标定日期的恒星时差;五、利用拟合的星敏感器慢变误差,预报星敏感器的姿态测量误差;六、根据星敏感器的测量模型,计算修正后的定姿角。本发明不仅可为图像导航与配准提供高精度的姿态基准,还可有效地提升光学遥感卫星的图像配准精度。
Description
技术领域
本发明涉及星敏感器慢变误差预报校正领域,特别涉及一种地球静止轨道卫星星敏慢变误差预报校正的方法。
背景技术
随着对地观测卫星的不断发展,人们对卫星的姿态控制精度和稳定性提出了更高的要求,特别是高精度光学遥感卫星,还要求卫星平台具有高姿态确定精度,为图像导航与配准提供高精度的姿态基准。由于空间环境温度及太阳光热效应的变化,对星敏感器的光学系统及成像光斑产生慢变影响,从而影响星敏感器的测量精度。星上常用的方法是采用高精度陀螺组合和高精度星敏感器联合卡尔曼滤波进行姿态确定实现高精度的姿态,但无法校正星敏感器在轨受空间环境温度及太阳光热效应影响产生的慢变误差。因此,为了提高图像导航与配准的姿态基准,有必要利用地球静止轨道卫星太阳光热效应重复性的优势,研究地球静止轨道卫星星敏感器慢变误差预报校正的方法。
发明内容
本发明的目的是提供一种地球静止轨道卫星星敏慢变误差预报校正的方法,它利用地球静止轨道卫星数据重复的优势,将在轨下传的星敏测量数据,根据载荷标定的姿态偏差,采用误差数据拟合方法得到星敏感器的慢变误差模型,将慢变误差上注,用于预报校正后期的星敏在轨慢变误差,为图像导航与配准提供高精度的姿态基准。
本发明的一种地球静止轨道卫星星敏慢变误差预报校正的方法,其步骤为:
S1:根据载荷标定的姿态确定误差,计算星敏感器的姿态测量误差;
S2:根据星敏感器的安装矩阵,计算星敏感器测量系下的测量误差;
S3:将所述星敏感器测量系下的测量误差数据通过分段多阶傅立叶函数拟合出星敏感器的慢变误差;
S4:根据地方时与恒星时角的转换关系,计算预报日期与标定日期的恒星时差;
S5:利用拟合的星敏感器的慢变误差,预报星敏感器的姿态测量误差;
S6:根据星敏感器的测量模型,计算修正后的定姿角。
优选地,所述步骤S1中,星敏感器的姿态测量误差矩阵Cbbr为:
其中,Cbo是星敏感器的姿态矩阵;Cbro是载荷标定的姿态矩阵。
优选地,星敏感器的姿态矩阵Cbo和载荷标定的姿态矩阵Cbro的计算公式分别为:
优选地,所述步骤S2中,星敏感器测量系下的误差矩阵Cssr为:
星敏感器测量系下的测量误差分别为:
优选地,所述星敏感器的安装矩阵Cbs的计算公式为:
其中,qbs=[qbs0qbs1qbs2qbs3]为标定的星敏测量安装四元数。
优选地,所述步骤S3中,采用分段多阶傅立叶函数将所述星敏感器测量系下的测量误差数据进行拟合,得到与时间对应的误差曲线。
优选地,采用分段四阶傅立叶函数将所述星敏感器测量系下的测量误差数据进行拟合,得到的误差曲线的表达式为:
其中,
(aθ0,aθ1,aθ2,aθ3,aθ4,bθ1,bθ2,bθ3,bθ4)是傅立叶函数拟合俯仰误差得到的系数;
ωθ是傅立叶函数拟合俯仰误差得到的基频;
θf(t)是俯仰误差拟合曲线函数;
(aψ0,aψ1,aψ2,aψ3,aψ4,bψ1,bψ2,bψ3,bψ4)是傅立叶函数拟合偏航误差得到的系数;
ωψ是傅立叶函数拟合偏航误差得到的基频;
ψf(t)是偏航误差拟合曲线函数;
t是距起始时刻的时间。
优选地,所述步骤S4中,预报日期t1与标定日期t0的恒星时差ΔG为:
其中,T1是t1的世纪数;T0是t0的世纪数。
星敏感器测量误差元数qerf的计算公式为:
其中,
θf为预报的星敏感器俯仰角误差,
ψf为预报的星敏感器偏航角误差。
优选地,所述步骤S6中,根据预报的星敏感器测量误差,修正星敏感器输出四元数;则星敏感器的测量模型为:
其中,
qbro=[qbro0 qbro1 qbro2 qbro3]T是修正后的星敏感器姿态四元数;
qio是轨道系到惯性系的转换四元数;
qsi是星敏感器的测量四元数;
qerf=[qerf0 qerf1 qerf2 qerf3]T是预报的星敏感测量误差四元数;
θc=atan2(2(qbro1qbro3+qbro2qbro0),-(qbro1)2-(qbro2)2+(qbro3)2+(qbro0)2);
ψc=atan2(2(qbro1qbro2+qbo3qbo0),-(qbro1)2+(qbro2)2-(qbro3)2+(qbro0)2)
其中,
θc为星敏感器修正的俯仰角,
ψc为星敏感器修正的偏航角。
与现有技术相比,本发明的有益效果为:本发明利用地球静止同步轨道卫星太阳光热效应重复性,采用载荷标定的姿态误差确定在轨下传数据的星敏感器慢变误差,根据重复规律预报在轨星敏的慢变误差,为图像导航与配准提供高精度的姿态基准,有效地提升光学遥感卫星的图像配准精度。
附图说明
图1为本发明的预报校正过程图。
具体实施方式
本发明提供了一种地球静止轨道卫星星敏慢变误差预报校正的方法,为使本发明更明显易懂,下面结合附图与具体实施方式对本发明做进一步说明。
如图1所示,本发明的一种地球静止轨道卫星星敏慢变误差预报校正的方法,其步骤具体如下:
步骤S1:根据载荷标定的姿态确定误差,计算星敏感器的姿态测量误差;
则星敏感器的姿态测量误差矩阵Cbbr为:
其中:
θr为载荷的俯仰角,
ψr为星敏感器的偏航角,
θ为星敏感器的俯仰角,
ψ为星敏感器的偏航角,
Cbo是星敏感器的姿态矩阵;
Cbro是载荷标定的姿态矩阵。
步骤S2:根据星敏感器安装矩阵计算星敏感器测量系下的测量误差;
设qbs=[qbs0qbs1qbs2qbs3]为标定的星敏测量安装四元数,则星敏感器的安装矩阵Cbs(即星敏感器测量系到本体系的转换矩阵)为:
星敏感器测量系下的误差矩阵Cssr为:
星敏感器测量系下的测量误差如下:
其中,
θse是星敏测量系下的俯仰角误差;
ψse是星敏测量系下的偏航角误差。
步骤S3:将星敏感器测量系下的测量误差数据采用分段多阶(四阶或四阶以上)傅立叶函数拟合出星敏感器的慢变误差;
将星敏感器测量系下的测量误差数据采用分段四阶傅立叶函数拟合,得到与时间对应的误差曲线表达式,如下:
其中,
(aθ0,aθ1,aθ2,aθ3,aθ4,bθ1,bθ2,bθ3,bθ4)是傅立叶函数拟合俯仰误差得到的系数;
ωθ是傅立叶函数拟合俯仰误差得到的基频;
θf(t)是俯仰误差拟合曲线函数;
(aψ0,aψ1,aψ2,aψ3,aψ4,bψ1,bψ2,bψ3,bψ4)是傅立叶函数拟合偏航误差得到的系数;
ωψ是傅立叶函数拟合偏航误差得到的基频;
ψf(t)是偏航误差拟合曲线函数;
t是距起始时刻的时间(单位:s)。
步骤S4:根据地方时与恒星时角的转换关系,计算预报日期与标定日期的恒星时差(地方时相同);
设标定的日期为t0,预报的日期为t1,可得:
其中,
ΔG是恒星时差,单位(s);
T1是t1的世纪数(天的整数);
T0是t0的世纪数(天的整数)。
步骤S5:利用拟合的星敏感器慢变误差预报星敏感器的姿态测量误差;
星敏感器测量误差元数qerf如下:
其中:
θf为预报的星敏感器俯仰角误差,
ψf为预报的星敏感器偏航角误差。
步骤S6:根据星敏感器的测量模型,计算修正后的定姿角;
根据预报的星敏感器测量误差修正星敏感器输出四元数,因设qbs为标定的星敏测量安装四元数,则星敏感器的测量模型为:
其中,
qbro=[qbro0 qbro1 qbro2 qbro3]T是修正后的星敏感器姿态四元数;
qio是轨道系到惯性系的转换四元数;
qsi是星敏感器的测量四元数;
qerf=[qerf0 qerf1 qerf2 qerf3]T是预报的星敏感测量误差四元数。
θc=atan2(2(qbro1qbro3+qbro2qbro0),-(qbro1)2-(qbro2)2+(qbro3)2+(qbro0)2);
ψc=atan2(2(qbro1qbro2+qbo3qbo0),-(qbro1)2+(qbro2)2-(qbro3)2+(qbro0)2)
其中:
θc为星敏感器修正的俯仰角,
ψc为星敏感器修正的偏航角。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (9)
1.一种地球静止轨道卫星星敏慢变误差预报校正的方法,其特征在于,其步骤为:
S1:根据载荷标定的姿态确定误差,计算星敏感器的姿态测量误差;
S2:根据星敏感器的安装矩阵,计算星敏感器测量系下的测量误差;
S3:将所述星敏感器测量系下的测量误差数据通过分段多阶傅立叶函数拟合出星敏感器的慢变误差;
S4:根据地方时与恒星时角的转换关系,计算预报日期与标定日期的恒星时差;
S5:利用拟合的星敏感器的慢变误差,预报星敏感器的姿态测量误差;
S6:根据星敏感器的测量模型,计算修正后的定姿角;
所述步骤S6中,根据预报的星敏感器测量误差,修正星敏感器输出四元数;则星敏感器的测量模型为:
其中,
qbro=[qbro0 qbro1 qbro2 qbro3]T是修正后的星敏感器姿态四元数;
qio是轨道系到惯性系的转换四元数;
qsi是星敏感器的测量四元数;
qerf=[qerf0 qerf1 qerf2 qerf3]T是预报的星敏感测量误差四元数;
θc=atan2(2(qbro1qbro3+qbro2qbro0),-(qbro1)2-(qbro2)2+(qbro3)2+(qbro0)2);
ψc=atan2(2(qbro1qbro2+qbo3qbo0),-(qbro1)2+(qbro2)2-(qbro3)2+(qbro0)2)
其中,
θc为星敏感器修正的俯仰角,
ψc为星敏感器修正的偏航角。
6.如权利要求1所述的一种地球静止轨道卫星星敏慢变误差预报校正的方法,其特征在于,
所述步骤S3中,采用分段多阶傅立叶函数将所述星敏感器测量系下的测量误差数据进行拟合,得到与时间对应的误差曲线。
7.如权利要求5所述的一种地球静止轨道卫星星敏慢变误差预报校正的方法,其特征在于,
采用分段四阶傅立叶函数将所述星敏感器测量系下的测量误差数据进行拟合,得到的误差曲线的表达式为:
其中,
(aθ0,aθ1,aθ2,aθ3,aθ4,bθ1,bθ2,bθ3,bθ4)是傅立叶函数拟合俯仰误差得到的系数;
ωθ是傅立叶函数拟合俯仰误差得到的基频;
θf(t)是俯仰误差拟合曲线函数;
(aψ0,aψ1,aψ2,aψ3,aψ4,bψ1,bψ2,bψ3,bψ4)是傅立叶函数拟合偏航误差得到的系数;
ωψ是傅立叶函数拟合偏航误差得到的基频;
ψf(t)是偏航误差拟合曲线函数;
t是距起始时刻的时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710681151.1A CN107402023B (zh) | 2017-08-10 | 2017-08-10 | 地球静止轨道卫星星敏慢变误差预报校正的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710681151.1A CN107402023B (zh) | 2017-08-10 | 2017-08-10 | 地球静止轨道卫星星敏慢变误差预报校正的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107402023A CN107402023A (zh) | 2017-11-28 |
CN107402023B true CN107402023B (zh) | 2020-06-23 |
Family
ID=60396544
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710681151.1A Active CN107402023B (zh) | 2017-08-10 | 2017-08-10 | 地球静止轨道卫星星敏慢变误差预报校正的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107402023B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108759821A (zh) * | 2018-06-14 | 2018-11-06 | 上海卫星工程研究所 | 一种geo卫星多个星敏感器构型 |
CN109443380B (zh) * | 2018-10-10 | 2022-02-18 | 中国科学院上海技术物理研究所 | 一种地球静止轨道光学相机热变形误差修正方法 |
CN111121787B (zh) * | 2019-12-06 | 2022-01-11 | 上海航天控制技术研究所 | 一种基于遥感图像的自主初轨确定方法 |
CN114088112A (zh) * | 2021-10-27 | 2022-02-25 | 中国空间技术研究院 | 一种卫星姿态确定精度评估方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1948085A (zh) * | 2005-10-12 | 2007-04-18 | 北京航空航天大学 | 一种基于星场的星敏感器校准方法 |
US7260456B2 (en) * | 2004-01-05 | 2007-08-21 | The Boeing Company | Pixel-frequency slews and filters for star data measurements |
CN104729537A (zh) * | 2015-03-19 | 2015-06-24 | 北京控制工程研究所 | 一种星敏感器低频误差在轨实时补偿方法 |
CN106940196A (zh) * | 2017-03-30 | 2017-07-11 | 上海航天控制技术研究所 | 一种星敏感器安装热变形修正方法 |
CN107024228A (zh) * | 2017-04-12 | 2017-08-08 | 上海航天控制技术研究所 | 一种星敏感器非高频误差在轨修正方法 |
-
2017
- 2017-08-10 CN CN201710681151.1A patent/CN107402023B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7260456B2 (en) * | 2004-01-05 | 2007-08-21 | The Boeing Company | Pixel-frequency slews and filters for star data measurements |
CN1948085A (zh) * | 2005-10-12 | 2007-04-18 | 北京航空航天大学 | 一种基于星场的星敏感器校准方法 |
CN104729537A (zh) * | 2015-03-19 | 2015-06-24 | 北京控制工程研究所 | 一种星敏感器低频误差在轨实时补偿方法 |
CN106940196A (zh) * | 2017-03-30 | 2017-07-11 | 上海航天控制技术研究所 | 一种星敏感器安装热变形修正方法 |
CN107024228A (zh) * | 2017-04-12 | 2017-08-08 | 上海航天控制技术研究所 | 一种星敏感器非高频误差在轨修正方法 |
Non-Patent Citations (1)
Title |
---|
星敏感器慢变误差校准方法研究;徐樱 等;《上海航天》;20160831;第33卷(第4期);第63-69页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107402023A (zh) | 2017-11-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107402023B (zh) | 地球静止轨道卫星星敏慢变误差预报校正的方法 | |
CN103323026B (zh) | 星敏感器和有效载荷的姿态基准偏差估计与修正方法 | |
CN108413887B (zh) | 光纤光栅辅助分布式pos的机翼形变测量方法、装置和平台 | |
Mazarico et al. | The gravity field, orientation, and ephemeris of Mercury from MESSENGER observations after three years in orbit | |
Tapley et al. | The gravity recovery and climate experiment: Mission overview and early results | |
US8433515B2 (en) | Method for measuring precision of star sensor and system using the same | |
CN104729537B (zh) | 一种星敏感器低频误差在轨实时补偿方法 | |
Olsen et al. | Magnetic observations from CryoSat-2: calibration and processing of satellite platform magnetometer data | |
CN104792340A (zh) | 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法 | |
CN103852086B (zh) | 一种基于卡尔曼滤波的光纤捷联惯导系统现场标定方法 | |
CN112082574B (zh) | 星敏感器的校正方法及系统 | |
CN102564455A (zh) | 星敏感器安装误差四位置标定与补偿方法 | |
Iwata | Precision attitude and position determination for the Advanced Land Observing Satellite (ALOS) | |
Lu et al. | Dual-filter transfer alignment for airborne distributed POS based on PVAM | |
CN109708663A (zh) | 基于空天飞机sins辅助的星敏感器在线标定方法 | |
CN108072387B (zh) | 一种低精度敏感器在轨偏差的地面校正方法及系统 | |
Schindelegger et al. | Atmospheric effects on Earth rotation | |
Iwata et al. | Precision attitude determination for the advanced land observing satellite (ALOS): design, verification, and on-orbit calibration | |
Li et al. | Airborne position and orientation system for aerial remote sensing | |
CN109655080A (zh) | 一种数字式太阳敏感器在轨标定方法 | |
Lai et al. | In-flight quality evaluation of attitude measurements from STECE APS-01 star tracker | |
Mears et al. | Estimated errors in retrievals of ocean parameters from SSMIS | |
CN103941593B (zh) | 低轨卫星姿态仿真方法 | |
Mo et al. | A study of lunar contamination and on‐orbit performance of the NOAA 18 Advanced Microwave Sounding Unit–A | |
Hwang et al. | Validation of geostationary satellite orbit determination using single-station antenna tracking data |
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 |