CN111854764A - 基于星间测量信息的航天器姿态确定方法及系统 - Google Patents

基于星间测量信息的航天器姿态确定方法及系统 Download PDF

Info

Publication number
CN111854764A
CN111854764A CN202010696750.2A CN202010696750A CN111854764A CN 111854764 A CN111854764 A CN 111854764A CN 202010696750 A CN202010696750 A CN 202010696750A CN 111854764 A CN111854764 A CN 111854764A
Authority
CN
China
Prior art keywords
satellite
inter
vector
information
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.)
Granted
Application number
CN202010696750.2A
Other languages
English (en)
Other versions
CN111854764B (zh
Inventor
陈婷
容建刚
张科科
刘善伍
王俊
孙国文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Engineering Center for Microsatellites
Innovation Academy for Microsatellites of CAS
Original Assignee
Shanghai Engineering Center for Microsatellites
Innovation Academy for Microsatellites of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Engineering Center for Microsatellites, Innovation Academy for Microsatellites of CAS filed Critical Shanghai Engineering Center for Microsatellites
Priority to CN202010696750.2A priority Critical patent/CN111854764B/zh
Publication of CN111854764A publication Critical patent/CN111854764A/zh
Application granted granted Critical
Publication of CN111854764B publication Critical patent/CN111854764B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/24Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • Automation & Control Theory (AREA)
  • General Physics & Mathematics (AREA)
  • Navigation (AREA)

Abstract

本发明提供了一种基于星间测量信息的航天器姿态确定方法及系统,包括:测量本星和目标星的轨道信息,获取本星与目标星之间的相对位置信息;将与星间测量光电设备的视线轴对准相关的参数转换为表征航天器姿态特征的星间矢量;直接获取本星星上姿态敏感器的测量信息;通过双矢量定姿,将所述相对位置信息、所述星间矢量和所述测量信息进行融合,形成航天器姿态信息,根据所述航天器姿态信息进行航天器三轴姿态确定。

Description

基于星间测量信息的航天器姿态确定方法及系统
技术领域
本发明涉及航天器定姿技术领域,特别涉及一种基于星间测量信息的航天器姿态确定方法及系统。
背景技术
在现有文献中,有一部分文献对通过星间通讯对合作目标进行类GPS的多航天器协同导航进行了研究,也有一部分文献针对非合作目标,运用陀螺和星间测量信息,结合滤波方法,对非合作目标的相对导航和相对姿态确定进行了研究。可以看出,现有文献都着重对如何运用星间测量信息和星上其它敏感器信息来确定本星和目标星之间的相对姿态进行了研究,而忽略了星间测量信息其实就和星上其它敏感器信息一样,能够在一定程度上反映本星的姿态特征的基本事实,从而忽视了对充分将星间测量信息和星上其它敏感器信息联合起来确定本星的三轴姿态的研究,更没有充分利用本星和目标星之间的星间测量信息确定本星三轴姿态的工程实例。
发明内容
本发明的目的在于提供一种基于星间测量信息的航天器姿态确定方法及系统,以解决现有的姿态确定精度方法没有充分利用星上有效资源的问题。
为解决上述技术问题,本发明提供一种基于星间测量信息的航天器姿态确定方法,该方法包括:
测量本星和目标星的轨道信息,获取本星与目标星之间的相对位置信息;
将与星间测量光电设备的视线轴对准相关的参数转换为表征航天器姿态特征的星间矢量;
直接获取本星星上姿态敏感器的测量信息;
通过双矢量定姿,将所述相对位置信息、所述星间矢量和所述测量信息进行融合,形成航天器姿态信息,根据所述航天器姿态信息进行航天器三轴姿态确定。
可选的,在所述的基于星间测量信息的航天器姿态确定方法中,所述基于与星间测量光电设备的视线轴对准相关的参数转换为表征航天器姿态特征的星间矢量包括:
航天器本体上的星间测量光电伺服机构测量得到能够使自身视线轴对准目标星的视线俯仰角α和偏航角β,结合航天器自身视线轴在尚未执行对准操作时的在星间激光载荷测量坐标系中的方向矢量Pinit,以及星间测量光电伺服机构的驱动顺序为先j轴后i轴,j轴带动i轴转动,得到视线轴对准目标星后的视线轴在星间激光载荷测量坐标系中的方向矢量Pm
Pm=Ri(anglei)Rj(anglej)Pinit
其中,i=y,z,j=y,z,angley=α,anglez=β,且
Figure BDA0002591464610000021
Figure BDA0002591464610000022
可选的,在所述的基于星间测量信息的航天器姿态确定方法中,所述基于与星间测量光电设备的视线轴对准相处的参数转换为表征航天器姿态特征的星间矢量还包括:
设星间激光载荷测量坐标系的X轴、Y轴、Z轴在卫星本体坐标系下的方向矢量分别为Mx、My、Mz,所述星间激光载荷测量坐标系到卫星本体坐标系的转换矩阵为激光载荷的安装矩阵:
Mbm=[Mx My Mz]T
获得卫星本体坐标系下的星间矢量为:
Pb=MbmPm
可选的,在所述的基于星间测量信息的航天器姿态确定方法中,测量本星和目标星的轨道信息,获取本星与目标星之间的相对位置信息包括:
计算轨道系下本星和目标星之间的归一化相对位置矢量,包括:
由星上定轨模块提供本星惯性系下本星的轨道位置RAi和速度VAi,以及本星惯性系下目标星的轨道位置RBi,计算本星惯性系下本星指向目标星的归一化位置矢量:
PABi=[RBi-RAi];
对所述本星惯性系下本星指向目标星的归一化位置矢量PABi进行归一化,得到uPABi
对本星惯性系下本星的轨道位置RAi和速度VAi进行归一化,得到相应的单位矢量uRAi=[rx ry rz]T和uVAi=[vx vy vz]T,并通过计算得到本星惯性系到轨道系的转换矩阵MAoi
MAoi[0][0]=-(uRAi·uVAi)rx+(uRAi·uRAi)vx
MAoi[0][1]=-(uRAi·uVAi)ry+(uRAi·uRAi)vy
MAoi[0][2]=-(uRAi·uVAi)rz+(uRAi·uRAi)vz
MAoi[1][0]=vyrz-vzry
MAoi[1][1]=vzrx-vxrz
MAoi[1][2]=vxry-vyrx
MAoi[2][0]=-rx
MAoi[2][1]=-ry
MAoi[2][2]=-rz
由uPABi和MAoi计算轨道系下本星指向目标星的归一化位置矢量:uPABo=MAoiuPABi
可选的,在所述的基于星间测量信息的航天器姿态确定方法中,直接获取本星星上姿态敏感器的测量信息包括:
由星上其它敏感器信息计算本体系下其它矢量信息,以及计算轨道系下其它环境矢量信息;
由星上其它敏感器信息计算本体系下其它矢量信息包括:
所述星上其它敏感器包括磁强计、太阳敏感器及地平仪,所述星上其它敏感器输出相应敏感器单机坐标系中的第一单位矢量uHm,所述第一单位矢量包括地磁矢量、太阳矢量及地球矢量;
由所述星上其它敏感器的安装方式得到相应的敏感器单机安装矩阵MbH,根据所述敏感器单机安装矩阵MbH得到所述本体系下其它矢量信息:
uHb=MbHuHm
计算轨道系下其它环境矢量信息包括:
由星上环境矢量确定模块提供轨道系下的第二单位矢量uHo,所述第二单位矢量包括归一化的地磁矢量、归一化的太阳矢量及归一化的地心矢量。
可选的,在所述的基于星间测量信息的航天器姿态确定方法中,将所述相对位置信息、所述星间矢量和所述测量信息进行融合包括:
融合星间矢量和星上其它敏感器信息确定航天器三轴姿态,包括:
Vo1=uPABo
Figure BDA0002591464610000041
Vo3=Vo1×Vo2
Figure BDA0002591464610000042
Vb1=Pb
Figure BDA0002591464610000043
Vb3=Vb1×Vb2
Figure BDA0002591464610000044
轨道系到本星的本体坐标系的转换矩阵为:
Figure BDA0002591464610000045
本星的三轴姿态角为滚动角
Figure BDA0002591464610000046
俯仰角θ及偏航角ψ,按照以下步骤解算得出:
若MAbo[1,2]=1,ψ与θ角取前一时刻角度值近似当前角度,则
Figure BDA0002591464610000047
θ=θk,ψ=ψk
若MAbo[1,2]=-1,ψ与θ取前一时刻角度值近似当前角度,则
Figure BDA0002591464610000048
θ=θk,ψ=ψk
若MAbo[1,2]≠±1,则
Figure BDA0002591464610000051
θ=arctan 2(-MAbo[0,2],MAbo[2,2]),θ∈(-π,π],
ψ=arctan 2(-MAbo[1,0],MAbo[1,1]),ψ∈(-π,π];
当所述星上其它敏感器为地平仪时,直接输出滚动角
Figure BDA0002591464610000052
和俯仰角θ,偏航角通过以下计算得出:
ψ=arctan 2[(k1·uPABoy-k2·uPABox),(k1·uPABox+k2·uPABoy)]
其中,k1=cosθ·Pbx+sinθ·Pbz
Figure BDA0002591464610000053
本发明还提供一种基于星间测量信息的航天器姿态确定系统,所述基于星间测量信息的航天器姿态确定系统包括:
第一载荷,被配置为测量本星和目标星的轨道信息,获取本星与目标星之间的相对位置信息;
星间矢量模块,被配置为将与星间测量光电设备的视线轴对准相关的参数转换为表征航天器姿态特征的星间矢量;
第二载荷,被配置为直接提供本星星上姿态敏感器的测量信息;
融合模块,被配置为通过双矢量定姿,将所述相对位置信息、所述星间矢量和所述测量信息进行融合,形成航天器姿态信息,根据所述航天器姿态信息进行航天器三轴姿态确定。
在本发明提供的基于星间测量信息的航天器姿态确定方法及系统中,通过给出了将星上载荷输出的星间测量信息转换为可以表征航天器姿态特征的星间矢量的方法,并给出了将星间矢量与星上其它敏感器输出信息进行融合以确定航天器三轴姿态的方法,使得星间测量信息能够充分应用到航天器姿态在轨实时确定的工程实际中。本发明可应用于带有星间激光通信载荷的卫星的三轴姿态确定,为整星提供一种备份的较高精度的三轴姿态确定方法。
航天器姿态确定为航天器姿态控制提供必需的输入,是进行航天器姿态控制的基础和前提,多一种的姿态确定方法,就能为航天器的姿态有效控制提供多一层的保障,也就为航天器姿控系统增加一些可靠性,从而,充分利用星上资源,例如星上载荷资源的输出信息进行航天器的三轴姿态确定就显得尤其重要。
目前对星上载荷资源的输出信息的利用着重于对不同航天器之间的相对姿态的确定上,而忽略了在本星相对于参考坐标系(通常为质心轨道坐标系或惯性坐标系)的三轴姿态确定方面的应用。为此,本发明提供一种可充分应用星上载荷输出的星间测量信息实现航天器的较高精度的、而且可以用于工程实际的三轴姿态确定方法。
本发明充分应用星间测量光电伺服机构输出的视线俯仰角和偏航角,基于星间视线轴对准原理,获得了可以表征星体姿态特征的星间矢量;再融合星间矢量和星上其它敏感器的输出信息,给出了能够唯一确定航天器三轴姿态的具体方法。
本发明无需额外增加其它敏感器,只是充分利用星上载荷的输出信息获取能够表征航天器姿态的星间矢量,通过与星上已有敏感器的输出信息进行有效融合,即可实现航天器三轴姿态的较高精度确定,很好地适应了现代商业航天重量轻价格低的要求,而且整套算法原理简单清晰,运算量小,易于在航天器机载计算机中实现,为航天器在轨姿态的实时高精度确定增加了一层备份,为航天器的姿态有效控制增加了一道保障。
附图说明
图1是本发明一实施例的基于星间测量信息的航天器姿态确定方法示意图。
具体实施方式
以下结合附图和具体实施例对本发明提出的基于星间测量信息的航天器姿态确定方法及系统作进一步详细说明。根据下面说明和权利要求书,本发明的优点和特征将更清楚。需说明的是,附图均采用非常简化的形式且均使用非精准的比例,仅用以方便、明晰地辅助说明本发明实施例的目的。
发明人发现,在现有文献中,有一部分文献对通过星间通信对合作目标进行类GPS的多航天器协同导航进行了研究,也有一部分文献针对非合作目标,运用陀螺和星间测量信息,结合滤波方法,对非合作目标的相对导航和相对姿态确定进行了研究。可以看出,现有文献都着重对如何运用星间测量信息和星上其它敏感器信息来确定本星和目标星之间的相对姿态进行了研究,而忽略了星间测量信息其实就和星上其它敏感器信息一样,能够在一定程度上反映本星的姿态特征的基本事实,从而忽视了对充分将星间测量信息和星上其它敏感器信息联合起来确定本星的三轴姿态的研究,更没有充分利用本星和目标星之间的星间测量信息确定本星三轴姿态的工程实例。
因此,如何将星间通信载荷获得的星间测量信息转换为姿态确定需要的星间矢量信息,以及如何利用星间矢量信息和星上其它敏感器信息进行航天器的三轴姿态确定,这块领域的研究非常少,并未查阅或检索到相关专利库/文献库有类似资料。
现有技术着重进行了应用星间测量信息进行不同卫星之间相对姿态确定的研究,而忽略了对单星三轴姿态确定的研究,从而未能提出如何应用星间测量信息确定航天器三轴姿态的方法,更缺乏如何将此方法应用到工程实际在轨实时应用中的考虑。
目前技术缺乏如何将星上的星间通信载荷的输出信息转换为可表征航天器姿态特征的星间矢量的方法,难以进一步充分利用星上载荷的输出信息为航天器提供多一种的三轴姿态确定方法。
本发明的核心思想在于提供一种基于星间测量信息的航天器姿态确定方法及系统,以解决现有的姿态确定精度方法没有充分利用星上有效资源的问题。
为实现上述思想,本发明提供了一种基于星间测量信息的航天器姿态确定方法及系统,所述基于星间测量信息的航天器姿态确定系统包括:第一载荷,被配置为测量本星和目标星的轨道信息,获取本星与目标星之间的相对位置信息;星间矢量模块,被配置为将与星间测量光电设备的视线轴对准相关的参数转换为表征航天器姿态特征的星间矢量;第二载荷,被配置为直接提供本星星上姿态敏感器的测量信息;融合模块,被配置为通过双矢量定姿,将所述相对位置信息、所述星间矢量和所述测量信息进行融合,形成航天器姿态信息,根据所述航天器姿态信息进行航天器三轴姿态确定。
<实施例一>
本发明的实施例提供一种基于星间测量信息的航天器姿态确定方法,如图1所示,所述基于星间测量信息的航天器姿态确定方法包括:测量本星和目标星的轨道信息,获取本星与目标星之间的相对位置信息;将与星间测量光电设备的视线轴对准相关的参数转换为表征航天器姿态特征的星间矢量;直接获取本星星上姿态敏感器的测量信息;通过双矢量定姿,将所述相对位置信息、所述星间矢量和所述测量信息进行融合,形成航天器姿态信息,根据所述航天器姿态信息进行航天器三轴姿态确定。
具体的,在所述的基于星间测量信息的航天器姿态确定方法中,所述基于与星间测量光电设备的视线轴对准相关的参数转换为表征航天器姿态特征的星间矢量包括:航天器本体上的星间测量光电伺服机构测量得到能够使自身视线轴对准目标星的视线俯仰角α和偏航角β,结合航天器自身视线轴在尚未执行对准操作时的在星间激光载荷测量坐标系中的方向矢量Pinit,以及星间测量光电伺服机构的驱动顺序为先j轴后i轴,j轴带动i轴转动,得到视线轴对准目标星后的视线轴在星间激光载荷测量坐标系中的方向矢量Pm
Pm=Ri(anglei)Rj(anglej)Pinit
其中,i=y,z,j=y,z,angley=α,anglez=β,且
Figure BDA0002591464610000081
Figure BDA0002591464610000091
进一步的,在所述的基于星间测量信息的航天器姿态确定方法中,所述基于与星间测量光电设备的视线轴对准相处的参数转换为表征航天器姿态特征的星间矢量还包括:设星间激光载荷测量坐标系的X轴、Y轴、Z轴在卫星本体坐标系下的方向矢量分别为Mx、My、Mz,所述星间激光载荷测量坐标系到卫星本体坐标系的转换矩阵为激光载荷的安装矩阵:
Mbm=[Mx My Mz]T
获得卫星本体坐标系下的星间矢量为:
Pb=MbmPm
另外,在所述的基于星间测量信息的航天器姿态确定方法中,测量本星和目标星的轨道信息,获取本星与目标星之间的相对位置信息包括:计算轨道系下本星和目标星之间的归一化相对位置矢量,包括:由星上定轨模块提供本星惯性系下本星的轨道位置RAi和速度VAi,以及本星惯性系下目标星的轨道位置RBi,计算本星惯性系下本星指向目标星的归一化位置矢量:
PABi=[RBi-RAi];
对所述本星惯性系下本星指向目标星的归一化位置矢量PABi进行归一化,得到uPABi
对本星惯性系下本星的轨道位置RAi和速度VAi进行归一化,得到相应的单位矢量uRAi=[rx ry rz]T和uVAi=[vx vy vz]T,并通过计算得到本星惯性系到轨道系的转换矩阵MAoi
MAoi[0][0]=-(uRAi·uVAi)rx+(uRAi·uRAi)vx
MAoi[0][1]=-(uRAi·uVAi)ry+(uRAi·uRAi)vy
MAoi[0][2]=-(uRAi·uVAi)rz+(uRAi·uRAi)vz
MAoi[1][0]=vyrz-vzry
MAoi[1][1]=vzrx-vxrz
MAoi[1][2]=vxry-vyrx
MAoi[2][0]=-rx
MAoi[2][1]=-ry
MAoi[2][2]=-rz
由uPABi和MAoi计算轨道系下本星指向目标星的归一化位置矢量:uPABo=MAoiuPABi
具体的,在所述的基于星间测量信息的航天器姿态确定方法中,直接获取本星星上姿态敏感器的测量信息包括:由星上其它敏感器信息计算本体系下其它矢量信息,以及计算轨道系下其它环境矢量信息;由星上其它敏感器信息计算本体系下其它矢量信息包括:所述星上其它敏感器包括磁强计、太阳敏感器及地平仪,所述星上其它敏感器输出相应敏感器单机坐标系中的第一单位矢量uHm,所述第一单位矢量包括地磁矢量、太阳矢量及地球矢量;由所述星上其它敏感器的安装方式得到相应的敏感器单机安装矩阵MbH,根据所述敏感器单机安装矩阵MbH得到所述本体系下其它矢量信息:
uHb=MbHuHm
计算轨道系下其它环境矢量信息包括:由星上环境矢量确定模块提供轨道系下的第二单位矢量uHo,所述第二单位矢量包括归一化的地磁矢量、归一化的太阳矢量及归一化的地心矢量。
进一步的,在所述的基于星间测量信息的航天器姿态确定方法中,将所述相对位置信息、所述星间矢量和所述测量信息进行融合包括:融合星间矢量和星上其它敏感器信息确定航天器三轴姿态,包括:
Vo1=uPABo
Figure BDA0002591464610000101
Vo3=Vo1×Vo2
Figure BDA0002591464610000102
Vb1=Pb
Figure BDA0002591464610000103
Vb3=Vb1×Vb2
Figure BDA0002591464610000104
轨道系到本星的本体坐标系的转换矩阵为:
Figure BDA0002591464610000111
本星的三轴姿态角为滚动角
Figure BDA0002591464610000112
俯仰角θ及偏航角ψ,按照以下步骤解算得出:
若MAbo[1,2]=1,ψ与θ角取前一时刻角度值近似当前角度,则
Figure BDA0002591464610000113
θ=θk,ψ=ψk
若MAbo[1,2]=-1,ψ与θ取前一时刻角度值近似当前角度,则
Figure BDA0002591464610000114
θ=θk,ψ=ψk
若MAbo[1,2]≠±1,则
Figure BDA0002591464610000115
θ=arctan2(-MAbo[0,2],MAbo[2,2]),θ∈(-π,π],
ψ=arctan2(-MAbo[1,0],MAbo[1,1]),ψ∈(-π,π];
当所述星上其它敏感器为地平仪时,直接输出滚动角
Figure BDA0002591464610000116
和俯仰角θ,偏航角通过以下计算得出:
ψ=arctan 2[(k1·uPABoy-k2·uPABox),(k1·uPABox+k2·uPABoy)]
其中,k1=cosθ·Pbx+sinθ·Pbz
Figure BDA0002591464610000117
本实施例还提供一种基于星间测量信息的航天器姿态确定系统,所述基于星间测量信息的航天器姿态确定系统包括:第一载荷,被配置为测量本星和目标星的轨道信息,获取本星与目标星之间的相对位置信息;星间矢量模块,被配置为将与星间测量光电设备的视线轴对准相关的参数转换为表征航天器姿态特征的星间矢量;第二载荷,被配置为直接提供本星星上姿态敏感器的测量信息;融合模块,被配置为通过双矢量定姿,将所述相对位置信息、所述星间矢量和所述测量信息进行融合,形成航天器姿态信息,根据所述航天器姿态信息进行航天器三轴姿态确定。
在本发明提供的基于星间测量信息的航天器姿态确定方法及系统中,通过给出了将星上载荷输出的星间测量信息转换为可以表征航天器姿态特征的星间矢量的方法,并给出了将星间矢量与星上其它敏感器输出信息进行融合以确定航天器三轴姿态的方法,使得星间测量信息能够充分应用到航天器姿态在轨实时确定的工程实际中。本发明可应用于带有星间激光通信载荷的卫星的三轴姿态确定,为整星提供一种备份的较高精度的三轴姿态确定方法。
航天器姿态确定为航天器姿态控制提供必需的输入,是进行航天器姿态控制的基础和前提,多一种的姿态确定方法,就能为航天器的姿态有效控制提供多一层的保障,也就为航天器姿控系统增加一些可靠性,从而,充分利用星上资源,例如星上载荷资源的输出信息进行航天器的三轴姿态确定就显得尤其重要。
目前对星上载荷资源的输出信息的利用着重于对不同航天器之间的相对姿态的确定上,而忽略了在本星相对于参考坐标系(通常为质心轨道坐标系或惯性坐标系)的三轴姿态确定方面的应用。为此,本发明提供一种可充分应用星上载荷输出的星间测量信息实现航天器的较高精度的、而且可以用于工程实际的三轴姿态确定方法。
本发明充分应用星间测量光电伺服机构输出的视线俯仰角和偏航角,基于星间视线轴对准原理,获得了可以表征星体姿态特征的星间矢量;再融合星间矢量和星上其它敏感器的输出信息,给出了能够唯一确定航天器三轴姿态的具体方法。
本发明无需额外增加其它敏感器,只是充分利用星上载荷的输出信息获取能够表征航天器姿态的星间矢量,通过与星上已有敏感器的输出信息进行有效融合,即可实现航天器三轴姿态的较高精度确定,很好地适应了现代商业航天重量轻价格低的要求,而且整套算法原理简单清晰,运算量小,易于在航天器机载计算机中实现,为航天器在轨姿态的实时高精度确定增加了一层备份,为航天器的姿态有效控制增加了一道保障。
综上,上述实施例对基于星间测量信息的航天器姿态确定方法及系统的不同构型进行了详细说明,当然,本发明包括但不局限于上述实施中所列举的构型,任何在上述实施例提供的构型基础上进行变换的内容,均属于本发明所保护的范围。本领域技术人员可以根据上述实施例的内容举一反三。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
上述描述仅是对本发明较佳实施例的描述,并非对本发明范围的任何限定,本发明领域的普通技术人员根据上述揭示内容做的任何变更、修饰,均属于权利要求书的保护范围。

Claims (7)

1.一种基于星间测量信息的航天器姿态确定方法,其特征在于,该方法包括:
测量本星和目标星的轨道信息,获取本星与目标星之间的相对位置信息;
将与星间测量光电设备的视线轴对准相关的参数转换为表征航天器姿态特征的星间矢量;
直接获取本星星上姿态敏感器的测量信息;
通过双矢量定姿,将所述相对位置信息、所述星间矢量和所述测量信息进行融合,形成航天器姿态信息,根据所述航天器姿态信息进行航天器三轴姿态确定。
2.如权利要求1所述的基于星间测量信息的航天器姿态确定方法,其特征在于,所述基于与星间测量光电设备的视线轴对准相关的参数转换为表征航天器姿态特征的星间矢量包括:
航天器本体上的星间测量光电伺服机构测量得到能够使自身视线轴对准目标星的视线俯仰角α和偏航角β,结合航天器自身视线轴在尚未执行对准操作时的在星间激光载荷测量坐标系中的方向矢量Pinit,以及星间测量光电伺服机构的驱动顺序为先j轴后i轴,j轴带动i轴转动,得到视线轴对准目标星后的视线轴在星间激光载荷测量坐标系中的方向矢量Pm
Pm=Ri(anglei)Rj(anglej)Pinit
其中,i=y,z,j=y,z,angley=α,anglez=β,且
Figure FDA0002591464600000011
Figure FDA0002591464600000012
3.如权利要求2所述的基于星间测量信息的航天器姿态确定方法,其特征在于,所述基于与星间测量光电设备的视线轴对准相处的参数转换为表征航天器姿态特征的星间矢量还包括:
设星间激光载荷测量坐标系的X轴、Y轴、Z轴在卫星本体坐标系下的方向矢量分别为Mx、My、Mz,所述星间激光载荷测量坐标系到卫星本体坐标系的转换矩阵为激光载荷的安装矩阵:
Mbm=[Mx My Mz]T
获得卫星本体坐标系下的星间矢量为:
Pb=MbmPm
4.如权利要求3所述的基于星间测量信息的航天器姿态确定方法,其特征在于,测量本星和目标星的轨道信息,获取本星与目标星之间的相对位置信息包括:
计算轨道系下本星和目标星之间的归一化相对位置矢量,包括:
由星上定轨模块提供本星惯性系下本星的轨道位置RAi和速度VAi,以及本星惯性系下目标星的轨道位置RBi,计算本星惯性系下本星指向目标星的归一化位置矢量:
PABi=[RBi-RAi];
对所述本星惯性系下本星指向目标星的归一化位置矢量PABi进行归一化,得到uPABi
对本星惯性系下本星的轨道位置RAi和速度VAi进行归一化,得到相应的单位矢量uRAi=[rx ry rz]T和uVAi=[vx vy vz]T,并通过计算得到本星惯性系到轨道系的转换矩阵MAoi
MAoi[0][0]=-(uRAi·uVAi)rx+(uRAi·uRAi)vx
MAoi[0][1]=-(uRAi·uVAi)ry+(uRAi·uRAi)vy
MAoi[0][2]=-(uRAi·uVAi)rz+(uRAi·uRAi)vz
MAoi[1][0]=vyrz-vzry
MAoi[1][1]=vzrx-vxrz
MAoi[1][2]=vxry-vyrx
MAoi[2][0]=-rx
MAoi[2][1]=-ry
MAoi[2][2]=-rz
由uPABi和MAoi计算轨道系下本星指向目标星的归一化位置矢量:uPABo=MAoiuPABi
5.如权利要求4所述的基于星间测量信息的航天器姿态确定方法,其特征在于,直接获取本星星上姿态敏感器的测量信息包括:
由星上其它敏感器信息计算本体系下其它矢量信息,以及计算轨道系下其它环境矢量信息;
由星上其它敏感器信息计算本体系下其它矢量信息包括:
所述星上其它敏感器包括磁强计、太阳敏感器及地平仪,所述星上其它敏感器输出相应敏感器单机坐标系中的第一单位矢量uHm,所述第一单位矢量包括地磁矢量、太阳矢量及地球矢量;
由所述星上其它敏感器的安装方式得到相应的敏感器单机安装矩阵MbH,根据所述敏感器单机安装矩阵MbH得到所述本体系下其它矢量信息:
uHb=MbHuHm
计算轨道系下其它环境矢量信息包括:
由星上环境矢量确定模块提供轨道系下的第二单位矢量uHo,所述第二单位矢量包括归一化的地磁矢量、归一化的太阳矢量及归一化的地心矢量。
6.如权利要求5所述的基于星间测量信息的航天器姿态确定方法,其特征在于,将所述相对位置信息、所述星间矢量和所述测量信息进行融合包括:
融合星间矢量和星上其它敏感器信息确定航天器三轴姿态,包括:
Vo1=uPABo
Figure FDA0002591464600000031
Vo3=Vo1×Vo2
Figure FDA0002591464600000032
Vb1=Pb
Figure FDA0002591464600000033
Vb3=Vb1×Vb2
Figure FDA0002591464600000034
轨道系到本星的本体坐标系的转换矩阵为:
Figure FDA0002591464600000041
本星的三轴姿态角为滚动角
Figure FDA0002591464600000042
俯仰角θ及偏航角ψ,按照以下步骤解算得出:
若MAbo[1,2]=1,ψ与θ角取前一时刻角度值近似当前角度,则
Figure FDA0002591464600000043
θ=θk,ψ=ψk
若MAbo[1,2]=-1,ψ与θ取前一时刻角度值近似当前角度,则
Figure FDA0002591464600000044
θ=θk,ψ=ψk
若MAbo[1,2]≠±1,则
Figure FDA0002591464600000045
θ=arctan2(-MAbo[0,2],MAbo[2,2]),θ∈(-π,π],
ψ=arctan2(-MAbo[1,0],MAbo[1,1]),ψ∈(-π,π];
当所述星上其它敏感器为地平仪时,直接输出滚动角
Figure FDA0002591464600000046
和俯仰角θ,偏航角通过以下计算得出:
ψ=arctan2[(k1·uPABoy-k2·uPABox),(k1·uPABox+k2·uPABoy)]
其中,k1=cosθ·Pbx+sinθ·Pbz
Figure FDA0002591464600000047
7.一种基于星间测量信息的航天器姿态确定系统,其特征在于,所述基于星间测量信息的航天器姿态确定系统包括:
第一载荷,被配置为测量本星和目标星的轨道信息,获取本星与目标星之间的相对位置信息;
星间矢量模块,被配置为将与星间测量光电设备的视线轴对准相关的参数转换为表征航天器姿态特征的星间矢量;
第二载荷,被配置为直接提供本星星上姿态敏感器的测量信息;
融合模块,被配置为通过双矢量定姿,将所述相对位置信息、所述星间矢量和所述测量信息进行融合,形成航天器姿态信息,根据所述航天器姿态信息进行航天器三轴姿态确定。
CN202010696750.2A 2020-07-20 2020-07-20 基于星间测量信息的航天器姿态确定方法及系统 Active CN111854764B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010696750.2A CN111854764B (zh) 2020-07-20 2020-07-20 基于星间测量信息的航天器姿态确定方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010696750.2A CN111854764B (zh) 2020-07-20 2020-07-20 基于星间测量信息的航天器姿态确定方法及系统

Publications (2)

Publication Number Publication Date
CN111854764A true CN111854764A (zh) 2020-10-30
CN111854764B CN111854764B (zh) 2022-05-27

Family

ID=73000917

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010696750.2A Active CN111854764B (zh) 2020-07-20 2020-07-20 基于星间测量信息的航天器姿态确定方法及系统

Country Status (1)

Country Link
CN (1) CN111854764B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117416531A (zh) * 2023-11-17 2024-01-19 中国科学院微小卫星创新研究院 基于光学载荷的姿态控制方法、系统及计算机可读介质

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5959576A (en) * 1997-08-12 1999-09-28 Honeywell Inc. Satellite attitude determination using GPS and intersatellite line of sight communications
CN101381004A (zh) * 2008-08-20 2009-03-11 南京航空航天大学 基于大气阻力的微小卫星编队飞行控制方法及控制装置
CN102494686A (zh) * 2011-10-17 2012-06-13 北京国科环宇空间技术有限公司 一种卫星姿态轨道确定系统及方法
CN102679985A (zh) * 2012-05-11 2012-09-19 北京航空航天大学 一种应用星间跟踪的航天器星座分散化自主导航方法
CN104034334A (zh) * 2014-06-05 2014-09-10 哈尔滨工程大学 一种小视场星敏感器的单星及双星定姿方法
CN104296752A (zh) * 2014-09-24 2015-01-21 上海卫星工程研究所 航天器天文测角测速组合自主导航系统及其导航方法
US20160047878A1 (en) * 2011-10-25 2016-02-18 The Boeing Company Combined orbit and attitude determination system and methods
CN108279010A (zh) * 2017-12-18 2018-07-13 北京时代民芯科技有限公司 一种基于多传感器的微小卫星姿态确定方法
CN108427427A (zh) * 2018-03-16 2018-08-21 北京控制工程研究所 一种航天器对地表定向目标姿态角计算方法
CN110044361A (zh) * 2019-03-07 2019-07-23 中国人民解放军63921部队 基于目标投影位置的光学载荷星上自主调度方法
CN111307139A (zh) * 2019-12-09 2020-06-19 北京航空航天大学 一种基于偏振/天文信息融合的航向与姿态确定方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5959576A (en) * 1997-08-12 1999-09-28 Honeywell Inc. Satellite attitude determination using GPS and intersatellite line of sight communications
CN101381004A (zh) * 2008-08-20 2009-03-11 南京航空航天大学 基于大气阻力的微小卫星编队飞行控制方法及控制装置
CN102494686A (zh) * 2011-10-17 2012-06-13 北京国科环宇空间技术有限公司 一种卫星姿态轨道确定系统及方法
US20160047878A1 (en) * 2011-10-25 2016-02-18 The Boeing Company Combined orbit and attitude determination system and methods
CN102679985A (zh) * 2012-05-11 2012-09-19 北京航空航天大学 一种应用星间跟踪的航天器星座分散化自主导航方法
CN104034334A (zh) * 2014-06-05 2014-09-10 哈尔滨工程大学 一种小视场星敏感器的单星及双星定姿方法
CN104296752A (zh) * 2014-09-24 2015-01-21 上海卫星工程研究所 航天器天文测角测速组合自主导航系统及其导航方法
CN108279010A (zh) * 2017-12-18 2018-07-13 北京时代民芯科技有限公司 一种基于多传感器的微小卫星姿态确定方法
CN108427427A (zh) * 2018-03-16 2018-08-21 北京控制工程研究所 一种航天器对地表定向目标姿态角计算方法
CN110044361A (zh) * 2019-03-07 2019-07-23 中国人民解放军63921部队 基于目标投影位置的光学载荷星上自主调度方法
CN111307139A (zh) * 2019-12-09 2020-06-19 北京航空航天大学 一种基于偏振/天文信息融合的航向与姿态确定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
XUELIANG BAI; PHILIPP HAGEL; XIAOFENG WU; SIZE XIAO: "Improved model predictive control for virtual satellite attitude control", 《PROCEEDINGS OF THE 31ST CHINESE CONTROL CONFERENCE》 *
杏建军,都晓宁,王威,高玉东: "星间相对测量在三星编队中的应用", 《空间科学学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117416531A (zh) * 2023-11-17 2024-01-19 中国科学院微小卫星创新研究院 基于光学载荷的姿态控制方法、系统及计算机可读介质
CN117416531B (zh) * 2023-11-17 2024-05-31 中国科学院微小卫星创新研究院 基于光学载荷的姿态控制方法、系统及计算机可读介质

Also Published As

Publication number Publication date
CN111854764B (zh) 2022-05-27

Similar Documents

Publication Publication Date Title
CN101788296B (zh) 一种sins/cns深组合导航系统及其实现方法
CN112082574B (zh) 星敏感器的校正方法及系统
CN106382927A (zh) 一种基于卫星识别的星敏感器自主导航方法
CN103017760B (zh) 一种大椭圆轨道火星探测器自主对火定向方法
CN104880192B (zh) 一种基于偏振罗盘的载体航向角计算方法
Rad et al. Optimal attitude and position determination by integration of INS, star tracker, and horizon sensor
CN110285815A (zh) 一种可在轨全程应用的微纳卫星多源信息姿态确定方法
CN107449420B (zh) 一种星箭一体化惯性姿态敏感器
CN110672130A (zh) 一种大失准角下惯性/偏振光组合导航系统ekf对准方法
CN110296719A (zh) 一种在轨标定方法
CN112129288A (zh) 一种基于偏振光/地磁航向约束的位姿估计方法及系统
CN111854764B (zh) 基于星间测量信息的航天器姿态确定方法及系统
CN109945854B (zh) 一种基于四元数差值融合星敏感器数据的方法
Somov et al. Guidance, navigation and control of a surveying satellite when an area imagery for disaster management.
Lu et al. A high-accuracy two-position alignment inertial navigation system for lunar rovers aided by a star sensor with a calibration and positioning function
CN115855063A (zh) 基于绝对姿态递推修正的交会对接敏感器数据预处理方法
CN114802818A (zh) 晨昏轨道卫星及其对日姿态计算方法、导引方法
CN115655284A (zh) 一种高动态下相对导航的时间对准与补偿方法
Mahdi et al. Attitude determination and control system design of KufaSat
CN115112117A (zh) 一种基于轨道罗盘的卫星姿态确定方法
Li et al. Attitude and heading measurement based on adaptive complementary Kalman filter for PS/MIMU integrated system
Hong et al. Application of EKF for missile attitude estimation based on “SINS/CNS” integrated guidance system
Fujita et al. Lessons Learned from On-orbit Gyroscope Malfunction and Recovery Operation of Microsatellite RISESAT
Zhang et al. A new method of single celestial-body sun positioning based on theory of mechanisms
CN113525720B (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