CN109655080B - 一种数字式太阳敏感器在轨标定方法 - Google Patents
一种数字式太阳敏感器在轨标定方法 Download PDFInfo
- Publication number
- CN109655080B CN109655080B CN201811527132.4A CN201811527132A CN109655080B CN 109655080 B CN109655080 B CN 109655080B CN 201811527132 A CN201811527132 A CN 201811527132A CN 109655080 B CN109655080 B CN 109655080B
- Authority
- CN
- China
- Prior art keywords
- sun sensor
- orbit
- digital sun
- coordinate
- digital
- 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
Abstract
本发明涉及一种数字式太阳敏感器在轨标定方法,包含:S1、根据星上时间和太阳星历数据,计算J2000惯性系下的太阳矢量;S2、根据星敏感器四元数,计算卫星本体系下的太阳矢量投影;S3、计算理论测量系下的太阳矢量投影;S4、结合数字式太阳敏感器的参数,计算数字式太阳敏感器的理论质心坐标;S5、修正数字式太阳敏感器的在轨质心坐标;S6、计算数字式太阳敏感器相对理论测量系的安装偏差角;S7、采用伪逆修正数字式太阳敏感器的参数;S8、优化数字式太阳敏感器的质心原点坐标。本发明利用实际在轨测量数据标定数字式太阳敏感器的参数,有效提高在轨数字式太阳敏感器的测量精度。
Description
技术领域
本发明涉及一种在轨标定方法,具体是指一种数字式太阳敏感器的在轨标定方法。
背景技术
数字式太阳敏感器作为空间飞行器的太阳方位测量姿态敏感器,其测量精度直接决定了在轨姿态确定的精度,对姿态确定精度要求较高的空间飞行器则尤为重要。
由于数字式太阳敏感器的探测器的玻璃盖片厚度H、折射率n、成像焦距f、坐标原点、安装矩阵等均会影响太阳角度解算精度。因此,为达到高精度测量需求,需要对这些参数进行精确标定。
现有技术中,标定技术主要通过地面太阳模拟器将探测器每一个像素点遍历测量,并通过地面测量数据进行标定。但是,太阳模拟器由于光源稳定性的原因,与真实太阳之间仍然存在一定的差别,并且转台精度无法精确评估,因此现有方法不能真实的对数字式太阳敏感器进行标定。
此外,在实际应用中存在产品安装偏差、平台变形、环境等多种因素的影响,导致在轨应用中会影响数字式太阳敏感器的测量精度。因此,利用在轨数据对影响太阳角度解算的参数进行标定,是提升数字式太阳敏感器测量精度的一项至关重要的工作。
基于上述,本发明提出一种数字式太阳敏感器在轨标定方法,以解决现有技术中存在的缺点和限制。
发明内容
本发明的目的是提供一种数字式太阳敏感器在轨标定方法,利用实际在轨测量数据标定数字式太阳敏感器的参数,有效提高在轨数字式太阳敏感器的测量精度。
为实现上述目的,本发明提供一种数字式太阳敏感器在轨标定方法,包含以下步骤:
S1、根据星上时间和太阳星历数据,计算J2000惯性系下的太阳矢量;
S2、根据星敏感器四元数,计算卫星本体系下的太阳矢量投影;
S3、计算理论测量系下的太阳矢量投影;
S4、结合数字式太阳敏感器的参数,计算数字式太阳敏感器的理论质心坐标;
S5、修正数字式太阳敏感器的在轨质心坐标;
S6、计算数字式太阳敏感器相对理论测量系的安装偏差角;
S7、采用伪逆修正数字式太阳敏感器的参数;
S8、优化数字式太阳敏感器的质心原点坐标。
所述的S1之前,还包含采集遥测数据的步骤,所述的遥测数据为:星上时间T(i);对应T(i)时刻的数字式太阳敏感器的在轨质心坐标XS(i)、YS(i);星敏感器四元数Qj(i)。
所述的S2中,具体为:根据J2000惯性系下的太阳矢量Si以及星敏感器四元数Qj(i),结合星敏感器的安装矩阵,计算得到卫星本体系下的太阳矢量投影Sb为:
Sb=Abm×Aj(i)×Si
其中,Abm为星敏感器由J2000惯性系转换至卫星本体系的转换矩阵,Aj(i)为与星敏感器四元数Qj(i)对应的矩阵;Si为J2000惯性系下的太阳矢量。
所述的S3中,具体为:根据卫星本体系下的太阳矢量投影Sb,结合数字式太阳敏感器的安装矩阵,计算得到理论测量系下的太阳矢量投影Ss为:
Ss=Atb×Sb
其中,Atb为数字式太阳敏感器由卫星本体系转换至理论测量系的转换矩阵。
所述的S5中,具体为:将数字式太阳敏感器的质心原点坐标的标定误差Δx0、Δy0作为优化量,修正遥测的数字式太阳敏感器的在轨质心坐标XS(i)、YS(i),得到修正后的在轨质心坐标X′S(i)、Y′S(i)为:
X′S(i)=Xs(i)+Δx0
Y′S(i)=YS(i)+Δy0
其中,Δx0、Δy0的初始值均为0。
所述的S6中,具体包含以下步骤:
其中,ΔX、ΔY为修正后的在轨质心坐标X′S(i)、Y′S(i)与理论质心坐标XL(i)、YL(i)之间的偏差量;
S63、分别将修正后的在轨质心坐标X′S(i)、Y′S(i),与理论质心坐标XL(i)、YL(i)转换成极坐标ρ(θS,RS)和ρ(θL,RL);
所述的S7中,具体包含以下步骤:
S72、计算理论质心坐标和修正后的在轨质心坐标的偏差ΔR,采用伪逆估算成像焦距f、玻璃盖片厚度H、折射率n的偏差值Δf、ΔH、Δn,具体为:
ΔR=RL-RS
其中,通过重复执行S4~S72的步骤迭代至少5次,计算得到Δf、ΔH、Δn;
S73、修正成像焦距f、玻璃盖片厚度H、折射率n,为:
f'=f+Δf
H'=H+ΔH
n'=n+Δn
得到修正后的数字式太阳敏感器的成像焦距f’、玻璃盖片厚度H’、折射率n’。
所述的S8中,以修正后的在轨质心坐标X′S(i)、Y′S(i)与理论质心坐标XL(i)、YL(i)之间的偏差量ΔX、ΔY的测量偏差最小量作为优化量,采用单纯形法优化质心原点坐标,得到优化后的Δx0’、Δy0’。
综上所述,本发明所提供的数字式太阳敏感器在轨标定方法,利用数字式太阳敏感器的实际在轨测量数据,对其成像焦距、玻璃盖片厚度、折射率、质心原点坐标、安装偏差角等参数进行标定,有效提高在轨数字式太阳敏感器的测量精度,即提高在轨姿态确定精度,克服了地面标定存在测量误差以及无法考虑安装偏差、平台变形、环境等因素引起的偏差等缺点,为数字式太阳敏感器在轨标定应用提供了有效的方法。
附图说明
图1为本发明中的数字式太阳敏感器在轨标定方法的流程图。
具体实施方式
以下结合图1,通过优选实施例对本发明的技术内容、构造特征、所达成目的及功效予以详细说明。
如图1所示,为本发明所提供的数字式太阳敏感器在轨标定方法,包含以下步骤:
S1、根据星上时间和太阳星历数据,计算J2000惯性系下的太阳矢量;
S2、根据星敏感器四元数,计算卫星本体系下的太阳矢量投影;
S3、计算理论测量系下的太阳矢量投影;
S4、结合数字式太阳敏感器的参数,计算数字式太阳敏感器的理论质心坐标;
S5、修正数字式太阳敏感器的在轨质心坐标;
S6、计算数字式太阳敏感器相对理论测量系的安装偏差角;
S7、采用伪逆修正数字式太阳敏感器的参数;
S8、优化数字式太阳敏感器的质心原点坐标。
所述的S1之前,还包含采集遥测数据的步骤,所述的遥测数据为:星上时间T(i);对应T(i)时刻的数字式太阳敏感器的在轨质心坐标XS(i)、YS(i);星敏感器四元数Qj(i)。
所述的S2中,具体为:根据J2000惯性系下的太阳矢量Si以及星敏感器四元数Qj(i),结合星敏感器的安装矩阵,计算得到卫星本体系下的太阳矢量投影Sb为:
Sb=Abm×Aj(i)×Si
其中,Abm为星敏感器由J2000惯性系转换至卫星本体系的转换矩阵,Aj(i)为与星敏感器四元数Qj(i)对应的矩阵;Si为J2000惯性系下的太阳矢量。
所述的S3中,具体为:根据卫星本体系下的太阳矢量投影Sb,结合数字式太阳敏感器的安装矩阵,计算得到理论测量系下的太阳矢量投影Ss为:
Ss=Atb×Sb
其中,Atb为数字式太阳敏感器由卫星本体系转换至理论测量系的转换矩阵。
其中,所述的成像焦距f、玻璃盖片厚度H、折射率n的初始值为地面标定值。
所述的S5中,具体为:将数字式太阳敏感器的质心原点坐标的标定误差Δx0、Δy0作为优化量,修正遥测的数字式太阳敏感器的在轨质心坐标XS(i)、YS(i),得到修正后的在轨质心坐标X′S(i)、Y′S(i)为:
X′S(i)=Xs(i)+Δx0
Y′S(i)=YS(i)+Δy0
其中,Δx0、Δy0的初始值均为0。
所述的S6中,具体包含以下步骤:
其中,ΔX、ΔY为修正后的在轨质心坐标X′S(i)、Y′S(i)与理论质心坐标XL(i)、YL(i)之间的偏差量;
S63、分别将修正后的在轨质心坐标X′S(i)、Y′S(i),与理论质心坐标XL(i)、YL(i)转换成极坐标ρ(θS,RS)和ρ(θL,RL);
所述的S7中,具体包含以下步骤:
S72、计算理论质心坐标和修正后的在轨质心坐标的偏差ΔR,采用伪逆估算成像焦距f、玻璃盖片厚度H、折射率n的偏差值Δf、ΔH、Δn,具体为:
ΔR=RL-RS
其中,通过重复执行S4~S72的步骤迭代至少5次,计算得到Δf、ΔH、Δn;
S73、修正成像焦距f、玻璃盖片厚度H、折射率n,为:
f'=f+Δf
H'=H+ΔH
n'=n+Δn
得到修正后的数字式太阳敏感器的成像焦距f’、玻璃盖片厚度H’、折射率n’。
所述的S8中,以修正后的在轨质心坐标X′S(i)、Y′S(i)与理论质心坐标XL(i)、YL(i)之间的偏差量ΔX、ΔY的测量偏差最小量作为优化量,采用单纯形法优化质心原点坐标,得到优化后的Δx0’、Δy0’。
综上所述,本发明所提供的数字式太阳敏感器在轨标定方法,利用数字式太阳敏感器的实际在轨测量数据,对其成像焦距、玻璃盖片厚度、折射率、质心原点坐标、安装偏差角等参数进行标定,有效提高在轨数字式太阳敏感器的测量精度,即提高在轨姿态确定精度,克服了地面标定存在测量误差以及无法考虑安装偏差、平台变形、环境等因素引起的偏差等缺点,为数字式太阳敏感器在轨标定应用提供了有效的方法。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (9)
1.一种数字式太阳敏感器在轨标定方法,其特征在于,包含以下步骤:
S1、根据星上时间和太阳星历数据,计算J2000惯性系下的太阳矢量;
S2、根据星敏感器四元数,计算卫星本体系下的太阳矢量投影;
S3、计算理论测量系下的太阳矢量投影;
S4、结合数字式太阳敏感器的参数,计算数字式太阳敏感器的理论质心坐标;
S5、修正数字式太阳敏感器的在轨质心坐标;
S6、计算数字式太阳敏感器相对理论测量系的安装偏差角;
S7、采用伪逆修正数字式太阳敏感器的参数;
S8、优化数字式太阳敏感器的质心原点坐标。
2.如权利要求1所述的数字式太阳敏感器在轨标定方法,其特征在于,所述的S1之前,还包含采集遥测数据的步骤,所述的遥测数据为:星上时间T(i);对应T(i)时刻的数字式太阳敏感器的在轨质心坐标XS(i)、YS(i);星敏感器四元数Qj(i)。
3.如权利要求2所述的数字式太阳敏感器在轨标定方法,其特征在于,所述的S2中,具体为:根据J2000惯性系下的太阳矢量Si以及星敏感器四元数Qj(i),结合星敏感器的安装矩阵,计算得到卫星本体系下的太阳矢量投影Sb为:
Sb=Abm×Aj(i)×Si
其中,Abm为星敏感器由J2000惯性系转换至卫星本体系的转换矩阵,Aj(i)为与星敏感器四元数Qj(i)对应的矩阵;Si为J2000惯性系下的太阳矢量。
4.如权利要求3所述的数字式太阳敏感器在轨标定方法,其特征在于,所述的S3中,具体为:根据卫星本体系下的太阳矢量投影Sb,结合数字式太阳敏感器的安装矩阵,计算得到理论测量系下的太阳矢量投影Ss为:
Ss=Atb×Sb
其中,Atb为数字式太阳敏感器由卫星本体系转换至理论测量系的转换矩阵。
6.如权利要求5所述的数字式太阳敏感器在轨标定方法,其特征在于,所述的S5中,具体为:将数字式太阳敏感器的质心原点坐标的标定误差Δx0、Δy0作为优化量,修正遥测的数字式太阳敏感器的在轨质心坐标XS(i)、YS(i),得到修正后的在轨质心坐标X′S(i)、Y′S(i)为:
X′S(i)=Xs(i)+Δx0
Y′S(i)=YS(i)+Δy0
其中,Δx0、Δy0的初始值均为0。
7.如权利要求6所述的数字式太阳敏感器在轨标定方法,其特征在于,所述的S6中,具体包含以下步骤:
其中,ΔX、ΔY为修正后的在轨质心坐标X′S(i)、Y′S(i)与理论质心坐标XL(i)、YL(i)之间的偏差量;
S63、分别将修正后的在轨质心坐标X′S(i)、Y′S(i),与理论质心坐标XL(i)、YL(i)转换成极坐标ρ(θS,RS)和ρ(θL,RL);
8.如权利要求7所述的数字式太阳敏感器在轨标定方法,其特征在于,所述的S7中,具体包含以下步骤:
S72、计算理论质心坐标和修正后的在轨质心坐标的偏差ΔR,采用伪逆估算成像焦距f、玻璃盖片厚度H、折射率n的偏差值Δf、ΔH、Δn,具体为:
ΔR=RL-RS
其中,通过重复执行S4~S72的步骤迭代至少5次,计算得到Δf、ΔH、Δn;
S73、修正成像焦距f、玻璃盖片厚度H、折射率n,为:
f'=f+Δf
H'=H+ΔH
n'=n+Δn
得到修正后的数字式太阳敏感器的成像焦距f’、玻璃盖片厚度H’、折射率n’。
9.如权利要求8所述的数字式太阳敏感器在轨标定方法,其特征在于,所述的S8中,以修正后的在轨质心坐标X′S(i)、Y′S(i)与理论质心坐标XL(i)、YL(i)之间的偏差量ΔX、ΔY的测量偏差最小量作为优化量,采用单纯形法优化质心原点坐标,得到优化后的Δx0’、Δy0’。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811527132.4A CN109655080B (zh) | 2018-12-13 | 2018-12-13 | 一种数字式太阳敏感器在轨标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811527132.4A CN109655080B (zh) | 2018-12-13 | 2018-12-13 | 一种数字式太阳敏感器在轨标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109655080A CN109655080A (zh) | 2019-04-19 |
CN109655080B true CN109655080B (zh) | 2020-05-29 |
Family
ID=66114398
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811527132.4A Active CN109655080B (zh) | 2018-12-13 | 2018-12-13 | 一种数字式太阳敏感器在轨标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109655080B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114777730B (zh) * | 2022-06-16 | 2022-09-02 | 航天宏图信息技术股份有限公司 | 一种地基太阳高度角的计算方法和装置 |
CN117516551A (zh) * | 2024-01-05 | 2024-02-06 | 中国西安卫星测控中心 | 一种卫星本体系太阳矢量的综合获得方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6152403A (en) * | 1998-11-11 | 2000-11-28 | Hughes Electronics Corporation | Gyroscopic calibration methods for spacecraft |
CN100414253C (zh) * | 2007-07-06 | 2008-08-27 | 北京航空航天大学 | 一种数字式太阳敏感器标定方法及装置 |
CN102175260B (zh) * | 2010-12-31 | 2012-11-14 | 北京控制工程研究所 | 一种自主导航系统误差校正方法 |
CN105928527B (zh) * | 2016-04-25 | 2019-04-30 | 航天东方红卫星有限公司 | 基于可见光地球敏感器和太阳敏感器的卫星姿态确定方法 |
CN106096148B (zh) * | 2016-06-14 | 2019-05-24 | 中国空间技术研究院 | 一种简单姿态控制下的大倾角轨道卫星太阳帆板指向方法 |
CN108072387B (zh) * | 2016-11-14 | 2021-06-22 | 上海微小卫星工程中心 | 一种低精度敏感器在轨偏差的地面校正方法及系统 |
CN106767906B (zh) * | 2016-11-30 | 2020-05-19 | 上海航天控制技术研究所 | 一种标定数字太阳敏感器主点、焦距的方法 |
CN108657470B (zh) * | 2018-05-14 | 2019-06-25 | 上海微小卫星工程中心 | 航天器帆板一维驱动机构停滞位置在轨辨识及修正方法 |
-
2018
- 2018-12-13 CN CN201811527132.4A patent/CN109655080B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109655080A (zh) | 2019-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103674063B (zh) | 一种光学遥感相机在轨几何定标方法 | |
CN100414253C (zh) | 一种数字式太阳敏感器标定方法及装置 | |
CN104729537B (zh) | 一种星敏感器低频误差在轨实时补偿方法 | |
CN107564057B (zh) | 顾及大气折光校正的高轨面阵光学卫星在轨几何标定方法 | |
CN105806369B (zh) | 一种星敏感器在轨光行差修正方法 | |
CN104807477B (zh) | 一种基于靶标控制点的卫星ccd阵列影像几何检校方法 | |
CN104154928A (zh) | 一种适用于惯性平台内置式星敏感器的安装误差标定方法 | |
CN107607127B (zh) | 一种基于外场的星敏感器内部参数标定及精度快速验证系统 | |
CN109655080B (zh) | 一种数字式太阳敏感器在轨标定方法 | |
CN107967700B (zh) | 大视场宽工作距双目相机的在轨几何校正及精度验证方法 | |
CN111798523A (zh) | 星相机在轨定标定姿及遥感影像几何定位方法、系统 | |
CN107024228A (zh) | 一种星敏感器非高频误差在轨修正方法 | |
CN106679695B (zh) | 一种基于星敏感器的焦平面热变形的测试方法 | |
CN109855652B (zh) | 星载激光测高仪指向角误差为非常数时的在轨标定方法 | |
CN109990801B (zh) | 基于铅垂线的水平仪装配误差标定方法 | |
CN112097794B (zh) | 遥感卫星载荷平台标定方法及系统 | |
CN107860400B (zh) | 遥感卫星图像定位整星级综合优化设计与分析方法 | |
CN113447043A (zh) | 一种基于gnss的卫星天文导航系统误差自主标定方法及系统 | |
CN107228683B (zh) | 一种多星敏感器间慢变误差实时在轨修正方法 | |
CN113776560B (zh) | 一种基于高精度三轴转台的陀螺标度因数测试方法 | |
CN109506645B (zh) | 一种星敏感器安装矩阵地面精测方法 | |
CN110674888A (zh) | 一种基于数据融合的头部姿态识别方法 | |
CN110702096B (zh) | 一种太阳敏感器的高阶误差补偿方法 | |
CN109489656B (zh) | 基于旋量的星敏定姿方法 | |
CN109581428B (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 |