CN108490973B - 航天器编队相对轨道确定方法及装置 - Google Patents

航天器编队相对轨道确定方法及装置 Download PDF

Info

Publication number
CN108490973B
CN108490973B CN201810355799.4A CN201810355799A CN108490973B CN 108490973 B CN108490973 B CN 108490973B CN 201810355799 A CN201810355799 A CN 201810355799A CN 108490973 B CN108490973 B CN 108490973B
Authority
CN
China
Prior art keywords
spacecraft
relative
orbit
motion parameter
formation
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
Application number
CN201810355799.4A
Other languages
English (en)
Other versions
CN108490973A (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201810355799.4A priority Critical patent/CN108490973B/zh
Publication of CN108490973A publication Critical patent/CN108490973A/zh
Application granted granted Critical
Publication of CN108490973B publication Critical patent/CN108490973B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10—Simultaneous control of position or course in three dimensions
    • G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • G05D1/104—Simultaneous control of position or course in three dimensions specially adapted for aircraft involving a plurality of aircrafts, e.g. formation flying

Landscapes

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

Abstract

本发明实施例公开一种航天器编队相对轨道确定方法及装置。所述方法包括：测量围绕目标对象飞行的航天器编队中的第一航天器相对于参考对象的第一运动参数；根据所述第一运动参数，确定出所述第一航天器绕所述目标对象飞行的第一轨道根数；测量所述航天器编队中第二航天器相对于所述第一航天器的第二运动参数；根据所述第一轨道根数及所述第二运动参数，确定出所述第二航天器绕所述目标对象飞行的第二轨道根数；结合所述第一轨道根数及所述第二轨道根数，确定出所述第二航天器相对于所述第一航天器运动的平均相对轨道根数。

Description

航天器编队相对轨道确定方法及装置
技术领域
本发明涉及航天技术领域，尤其涉及一种航天器编队相对轨道确定方法及装置。
背景技术
由多个航天器按照一定顺序编队飞行的组合，可以称之为航天器编队。航天器编队中的各个航天器保持预定的相对轨道，是维持航天器编队的飞行队列的前提，故航天器编队相对轨道确定直接影响航天器编队飞行控制的成败。例如，月球轨道航天器编对是绕月飞行的航天器编队。月球轨道航天器编队相对轨道确定采用地面站跟踪测量。对于月球大椭圆轨道航天器编队，由于轨道周期很长，且在近月点航天器速度很快以及测量资源的限制等问题，总是存在着航天器编队的飞行控制精度低的问题。
发明内容
有鉴于此，本发明实施例期望提供一种航天器相对轨道确定方法及装置，至少部分解决上述飞行控制精度低的问题。
为达到上述目的，本发明的技术方案是这样实现的：
第一方面，本发明实施例提供一种航天器编队相对轨道确定方法，包括：
测量围绕目标对象飞行的航天器编队中的第一航天器相对于参考对象的第一运动参数；
根据所述第一运动参数，确定出所述第一航天器绕所述目标对象飞行的第一轨道根数；
测量所述航天器编队中第二航天器相对于所述第一航天器的第二运动参数；
根据所述第一轨道根数及所述第二运动参数，确定出所述第二航天器绕所述目标对象飞行的第二轨道根数；
结合所述第一轨道根数及所述第二轨道根数，确定出所述第二航天器相对于所述第一航天器运动的平均相对轨道根数。
可选地，所述测量围绕目标对象飞行的航天器编队中的第一航天器相对于参考对象的第一运动参数，包括：
测量围绕月球飞行的航天器编队中的第一航天器相对于地球的第一运动参数。
可选地，所述第一运动参数包括：所述第一航天器运行预设时长所经过的位置和对应位置处的速度；
所述根据所述第一运动参数，确定出所述第一航天器绕所述目标对象飞行的第一轨道根数，包括：
根据所述第一航天器运行预设时长所经过的位置和对应位置处的速度，计算出所述第一航天器瞬时轨道根数。
可选地，所述第二运动参数包括：所述第二航天器运行一周相对于所述第一航天器的相对位置处的相对速度；
所述根据所述第一轨道根数及所述第二运动参数，确定出所述第二航天器绕所述目标对象飞行的第二轨道根数，包括：
根据所述第二运动参数，确定出所述第二航天器运行一周相对于所述参考对象的位置和对应位置处的速度；
根据所述第二航天器运行一周相对于所述参考对象的位置和对应位置处的速度，计算出所述第二航天器瞬时轨道根数。
可选地，所述测量所述航天器编队中第二航天器相对于所述第一航天器的第二运动参数，包括：
利用无线电测量所述第二航天器相对于所述第一航天器的相对运动参数，其中，所述相对运动参数包括以下至少之一：相对位置及相对速度；相对角度及相对速度；相对位置、相对角度及相对速度；
和/或，
利用图像采集测量所述第二航天器相对于所述第一航天器的相对运动参数，其中，所述相对运动参数包括以下至少之一：相对位置及相对速度；相对角度及相对速度；相对位置、相对角度及相对速度。
第二方面，本发明实施例提供一种航天器编队相对轨道确定装置，包括：
第一测量模块，用于测量围绕目标对象飞行的航天器编队中的第一航天器相对于参考对象的第一运动参数；
第一确定模块，用于根据所述第一运动参数，确定出所述第一航天器绕所述目标对象飞行的第一轨道根数；
第二测量模块，用于测量所述航天器编队中第二航天器相对于所述第一航天器的第二运动参数；
第二确定模块，用于根据所述第一轨道根数及所述第二运动参数，确定出所述第二航天器绕所述目标对象飞行的第二轨道根数；
第三确定模块，用于结合所述第一轨道根数及所述第二轨道根数，确定出所述第二航天器相对于所述第一航天器运动的平均相对轨道根数。
可选地，所述第一测量模块，用于测量围绕月球飞行的航天器编队中的第一航天器相对于地球的第一运动参数。
可选地，所述第一运动参数包括：所述第一航天器运行预设时长所经过的位置和对应位置处的速度；
所述第一确定模块，用于根据所述第一航天器运行预设时长所经过的位置和对应位置处的速度，计算出所述第一航天器瞬时轨道根数。
可选地，所述第二运动参数包括：所述第二航天器运行一周相对于所述第一航天器的相对位置处的相对速度；
所述第二确定模块，用于根据所述第二运动参数，确定出所述第二航天器运行一周相对于所述参考对象的位置和对应位置处的速度；根据所述第二航天器运行一周相对于所述参考对象的位置和对应位置处的速度，计算出所述第二航天器瞬时轨道根数。
可选地，所述第二测量模块，具体用于：
利用无线电测量所述第二航天器相对于所述第一航天器的相对运动参数，其中，所述相对运动参数包括以下至少之一：相对位置及相对速度；相对角度及相对速度；相对位置、相对角度及相对速度；
和/或，
利用图像采集测量所述第二航天器相对于所述第一航天器的相对运动参数，其中，所述相对运动参数包括以下至少之一：相对位置及相对速度；相对角度及相对速度；相对位置、相对角度及相对速度。
本发明实施例提供的航天器相对轨道确定方法及装置，仅对航天器编队中部分航天器进行相对于参考对象的远距离测量，而其他航天器则通过与进行了远距离测量的航天器之间的相对运动参数的测量，从而可以换算出其他航天器相对于参考对象的测量；如此减少了远距离测量，减少远距离测量导致的精确度低的问题，同时还减少了因为远距离测量导致的功耗大及计算繁琐等问题，从而具有测量精确度高、平均相对轨道根数预测精确的特点，同时还具有实现简单的特点。
附图说明
图1为本发明实施例提供的第一种航天器编队相对轨道确定方法的流程示意图；
图2为本发明实施例提供的第一种航天器编队相对轨道确定装置的结构示意图；
图3为本发明实施例提供的第二种航天器编队相对轨道确定方法的流程示意图；
图4为本发明实施例提供的一种用于测量第二运动参数的坐标系的示意图；
图5为本发明实施例提供的一种相对位置的误差效果示意图；
图6为本发明实施例提供的一种相对速度的误差效果示意图。
具体实施方式
以下结合说明书附图及具体实施例对本发明的技术方案做进一步的详细阐述。
如图1所示，本实施例提供一种航天器编队相对轨道确定方法，包括：
步骤S110：测量围绕目标对象飞行的航天器编队中的第一航天器相对于参考对象的第一运动参数；
步骤S120：根据所述第一运动参数，确定出所述第一航天器绕所述目标对象飞行的第一轨道根数；
步骤S130：测量所述航天器编队中第二航天器相对于所述第一航天器的第二运动参数；
步骤S140：根据所述第一轨道根数及所述第二运动参数，确定出所述第二航天器绕所述目标对象飞行的第二轨道根数；
步骤S150：结合所述第一轨道根数及所述第二轨道根数，确定出所述第二航天器相对于所述第一航天器运动的平均相对轨道根数。
本实施例中所述航天器编队可包括2个或2个以上的航天器，在本实施例中第一航天器和第二航天器可以泛指任意一个航天器。可选地为，所述第一航天器可为航天器编队中的核心航天器，第二航天器可为分布在核心航天器旁边的其他航天器。所述核心航天器通常为主航天器，其他航天器可为主航天器提供支撑和/或辅助服务的航天器。
在本实施例中，航天器编队绕目标对象飞行。例如，航天器编队为月球的卫星，则目标对象为月球；若航天器编队为地球的卫星，则目标对象为地球，若航天器为火星的卫星，则目标对象为火星。在本实施例中参考对象可为进行第一航天器的运动参数的基准点所在星球。所述参考对象可为月球的航天器编队参考的地球，也可以火星的航天器编队的地球。
在本实施例中可以通过地磁计等对地测量或者地面对卫星的远程监控，得到第一航天器相对于参考对象的第一运行参数，例如，第一航天器绕月球运动一周相对于地球的各个位置的速度。根据第一运行参数就可以按照任意一种方法计算出第一航天器的轨道根数，即所述第一轨道根数。轨道根数(或称轨道要素或轨道参数)是用来描述天体在其轨道运行状态的一组参数。通常情况下指的是用经典万有引力定律描述天体按圆锥曲线运动时多个参数，例如，以下提供6个参考的参数：
轨道半长轴a，椭圆轨道的半长轴，并不是长轴与短轴的算术平均数。若圆轨道，则为轨道高度加上地球半径。
轨道偏心率e，为椭圆扁平程度的一种量度，定义是椭圆两焦点间的距离与长轴长度的比值。
轨道倾角i，行星轨道面对黄道面的倾角；在升交点处从黄道面逆时针方向量到行星轨道面的角度。
升交点赤经Ω，行星轨道升交点的黄道经度。自X轴(春分点)方向在赤道平面内沿逆时针方向度量到升交点的地心夹角。从春分点向东度量，0至360°。
近地点幅角ω，从升交点沿行星运动轨道逆时针量到近地点的角度。自轨道升交点在轨道平面内沿卫星运动方向度量到近地点的角度。从升交点顺轨道运行方向度量，0-360°。
指定历元的平近点角，行星对应于初始时刻(t0)时该的平近点角。
位于同一个航天编队中的两个航天器的距离通常是较近的，然后测量第二航天器相对于第一航天器的第二运动参数，如此则不用采用远距离的对参考对象的远距离测量，如此可以提升测量的精度，从而获得高精度的第二运动参数，且相对于远距离测量还节省测量的功耗。
在本实施例中可以根据第二运动参数，同时计算出第二航天器的轨道根数，即所述第二轨道根数。
在结合第一轨道根数和第二轨道根数，计算出两个航天器之间的平均相对轨道根数。
如此，在进行平均相对轨道根数确定时，仅需确定出航天编队中一个或多个(即部分)航天器相对于参考对象的运动参数，其他航天器利用与已经相对于参考对象对应的航天器的相对运动测量，得到其运动参数；减少远距离测量导致的精度低及功耗高的问题，故本申请提供的平均相对轨道根数的确定具有测量精度高及功耗低的特点。
可选地，所述步骤S130可包括：测量围绕月球飞行的航天器编队中的第一航天器相对于地球的第一运动参数。
在本实施例中，月球为目标对象；地球为参考对象，即航天编队为绕月飞行的月球航天器编队。
可选地，所述第一运动参数包括：所述第一航天器运行预设时长所经过的位置和对应位置处的速度；所述步骤S120可包括：
根据所述第一航天器运行预设时长所经过的位置和对应位置处的速度，计算出所述第一航天器瞬时轨道根数。
例如，第一航天器绕月球一段时间所经过的位置及该位置的速度，从而可以计算出第一航天瞬时轨道根数。
可选地，所述第二运动参数包括：所述第二航天器运行一周相对于所述第一航天器的相对位置处的相对速度；所述步骤S140可包括：根据所述第二运动参数，确定出所述第二航天器运行一周相对于所述参考对象的位置和对应位置处的速度；根据所述第二航天器运行一周相对于所述参考对象的位置和对应位置处的速度，计算出所述第二航天器瞬时轨道根数。
在本实施例中，首先将测量出第二航天器相对于第一航天器的相对运动参数，然后结合第一航天器的第一运动参数，计算出第二航天器相对于参考对象的运动参数，例如，计算出月球航天编队中的第二航天器相对于地球的运动参数；然后基于第二航天器相对于地球的运动参数，确定出第二航天器的瞬时轨道根数。
所述步骤S150可包括：根据第一航天器的瞬时轨道根数及第二航天器的瞬时轨道根数，可以确定出第一航天器及第二航天器之间的平均相对轨道根数。平均相对轨道根数为相对轨道根数的均值。
可选地，所述步骤S130可包括：
利用无线电测量所述第二航天器相对于所述第一航天器的相对运动参数，其中，所述相对运动参数包括以下至少之一：相对位置及相对速度；相对角度及相对速度；相对位置、相对角度及相对速度；
和/或，
利用图像采集测量所述第二航天器相对于所述第一航天器的相对运动参数，其中，所述相对运动参数包括以下至少之一：相对位置及相对速度；相对角度及相对速度；相对位置、相对角度及相对速度。
例如，基于无线电测量可包括：利用各种交变电磁场进行的距离测量、角度测量或速度测量等。
利用图像采集测量包括采集图像，通过被采集的航天器在图像中尺寸和位置的变化，可以确定出两个航天器之间的相对距离、相对角度及相对速度等各种相对运动参数。例如，利用位于第一航天器上的各种测量仪器，例如，无线电测量仪器和/或照相机来测量相对运动参数。再例如，也可以利用位于第二航天器上的各种测量仪器，例如，无线电测量仪器和/或照相机来测量相对运功参数。
如图2所示，本实施例提供一种航天器编队相对轨道确定装置，包括：
第一测量模块110，用于测量围绕目标对象飞行的航天器编队中的第一航天器相对于参考对象的第一运动参数；
第一确定模块120，用于根据所述第一运动参数，确定出所述第一航天器绕所述目标对象飞行的第一轨道根数；
第二测量模块130，用于测量所述航天器编队中第二航天器相对于所述第一航天器的第二运动参数；
第二确定模块140，用于根据所述第一轨道根数及所述第二运动参数，确定出所述第二航天器绕所述目标对象飞行的第二轨道根数；
第三确定模块150，用于结合所述第一轨道根数及所述第二轨道根数，确定出所述第二航天器相对于所述第一航天器运动的平均相对轨道根数。
第一测量模块110、第一确定模块120、第二测量模块130、第二确定模块140及第三确定模块150均可为程序模块，被处理器或处理电路执行后，可以实现前述一个或多个实施例提供的功能，从而可以利用一个航天器对参考对象的远距离测量，并结合航天器之间的相对测量，计算出航天器之间的平均相对轨道根数，具有结果精确及功耗低等特点。
可选地，所述第一测量模块110，用于测量围绕月球飞行的航天器编队中的第一航天器相对于地球的第一运动参数。
可选地，所述第一运动参数包括：所述第一航天器运行预设时长所经过的位置和对应位置处的速度；
所述第一确定模块120，用于根据所述第一航天器运行预设时长所经过的位置和对应位置处的速度，计算出所述第一航天器瞬时轨道根数。
可选地，所述第二运动参数包括：所述第二航天器运行一周相对于所述第一航天器的相对位置处的相对速度；
所述第二确定模块140，用于根据所述第二运动参数，确定出所述第二航天器运行一周相对于所述参考对象的位置和对应位置处的速度；根据所述第二航天器运行一周相对于所述参考对象的位置和对应位置处的速度，计算出所述第二航天器瞬时轨道根数。
可选地，所述第二测量模块130，具体用于：利用无线电测量所述第二航天器相对于所述第一航天器的相对运动参数，其中，所述相对运动参数包括以下至少之一：相对位置及相对速度；相对角度及相对速度；相对位置、相对角度及相对速度；和/或，利用图像采集测量所述第二航天器相对于所述第一航天器的相对运动参数，其中，所述相对运动参数包括以下至少之一：相对位置及相对速度；相对角度及相对速度；相对位置、相对角度及相对速度。
以下结合上述任意实施例提供几个具体示例：
示例1：
如图3所示，本示例提供一种航天器编队的相对轨道确定的方法，可包括：
第一步：由地面站测得A星(相当于前述的第一航天器)初始时刻的绝对位置/速度，通过轨道预报得到一圈内A星的绝对位置/速度，并由此计算出一圈内A星的瞬时轨道根数。
第二步：由航天器自身携带的测量仪器短时间测得两个航天器的相对距离和角度的测量值，再通过相对测量获得其估计值。由A星初始时刻的绝对位置/速度和航天器相对距离和角度的估计值，得到初始时刻B星(相当于前述的第二航天器)相对位置/速度。再通过轨道预报得到一圈内B星的绝对位置/速度，并由此计算出一圈内B星的瞬时轨道根数。此时的绝对位置和速度都是对地位置和对地速度，即以地球为参考点测量的位置及速度。
第三步：由A星和B星的瞬时轨道根数计算出一圈的瞬时相对轨道根数；最终，由瞬时相对轨道根数计算出平均相对轨道根数。该方式计算平均相对轨道根数简单易行，可以准确、可靠地确定月球大椭圆轨道航天器编队相对轨道参数。
示例2：
本示例提供一种基于星上短时间相对测量的月球大椭圆轨道航天器编队相对轨道确定方法，包括：
联合地面站短时间绝对轨道测量和星上短时间相对测量，高精度确定月球大椭圆轨道航天器编队相对轨道；所述短时间可为一个航天器绕月球飞行一周的时间或数周的时间。
利用地面站较短时间测量所测得主星的初始轨道根数，通过轨道递推得到一个周期内主星的轨道根数；
采用无线电测距和测角相机测角的信息测量方式，仅需较短时间内相对测量信息，通过相对测量和相关计算得到双星的相对平均轨道根数。
联合地面站精度较低的较短时间绝对轨道定轨数据和较短时间星上相对距离、相对角度测量数据，高精度的确定月球大椭圆轨道航天器编队相对平均轨道根数。
示例3：
相对测量：
采用无线电测距加测角相机测角的信息测量方式，测距精度为2m(1σ)，测角精度为0.02°(1σ)。
如图4所示，定义测量坐标系OcXcYcZc：坐标原点Oc为测角相机光心，可以认为与卫星质心重合，Xc轴指向光轴方向，Yc轴与Zc轴对应测角相机成像的水平方向和竖直方向。Yc与Zc轴和Xc轴构成右手直角坐标系。
相对测量采用扩展卡尔曼滤波，其计算过程如下：
时间更新方程:
Figure GDA0002907323540000111
其中
Figure GDA0002907323540000112
f(X,t)是如下的非线性相对动力学方程：其中，x为在Xc轴的位置；y为在Yc轴上的位置；z为在Zc轴上的位置；
Figure GDA0002907323540000113
为在Xc轴的速度；
Figure GDA0002907323540000114
为在Yc轴上的速度；
Figure GDA0002907323540000115
为Zc轴的速度；其中，
Figure GDA0002907323540000116
中T表示转置。
Figure GDA0002907323540000117
Figure GDA0002907323540000118
Figure GDA0002907323540000119
其中，
Figure GDA00029073235400001110
为轨道角速度；
Figure GDA00029073235400001111
为轨道角加速度；
Figure GDA00029073235400001112
为在Xc轴的加速度分量；
Figure GDA00029073235400001113
为在Yc轴的加速度分量；
Figure GDA00029073235400001114
为在Zc轴的加速度分量；μ为引力常数；r为航天器到参考对象的距离，例如，绕月飞行的航天器到地球中心的距离。
根据上一时刻的相对状态估计值，通过四阶龙哥库塔法求解非线性相对动力学方程(1)，得到当前时间的状态预测值。
Figure GDA0002907323540000121
式中状态转移矩阵如下
Figure GDA0002907323540000122
其中τ为滤波步长，本发明实施例选为1秒。k＝1,2,..为滤波常数；f为真近点角。
Figure GDA0002907323540000123
为
Figure GDA0002907323540000124
的转置矩阵；Qk-1过程激励噪声的协方差；
测量更新方程：
Figure GDA0002907323540000125
Figure GDA0002907323540000126
Pk＝(I-KkHk)Pk,k-1 (7)
其中，Rk为测量噪声协方差；Zk为相对坐标系下的相对位置的测量值，Zk通过如下的方程求解；
Figure GDA0002907323540000127
其中，Cc2b为测量坐标系到卫星本体坐标系的转换矩阵，Cb2LVLH为卫星本体坐标系到相对坐标系的转换矩阵。此处的本体坐标系可为航天器自身的坐标系，通常以航天器的质心为原点的坐标系；此处的相对坐标系可为前述的测量坐标系OcXcYcZc。
测量矩阵Hk为：
Figure GDA0002907323540000131
卫星运行于远地点附近时，相对速度较小，导航收敛快。而卫星运行于近地点附近，相对速度较大，导航收敛慢。因此相对测量选择在远地点附近进行。
如图5和图6所示，运用现有的Matlab软件进行仿真，得到相对位置和速度估计误差均在要求之内。其中，相对测量时间可为10分钟等预设时间长度。
由上述相对测量算法可估计出相对位置、相对速度，考虑到地面测定轨***可测量参考卫星的绝对位置速度。通过计算便可得到从卫星的绝对位置和速度。把两颗卫星的绝对位置速度分别转换成瞬时轨道根数，通过相减便可得瞬时相对轨道根数的估计。再通过如下相关计算可得平均相对轨道根数。
平均相对轨道根数计算：
Figure GDA0002907323540000132
的估计：
对于
Figure GDA0002907323540000133
的估计采用平均法。首先，将两颗卫星月惯系下绝对位置和速度转换成轨道根数的六个参数；然后，一个轨道周期数据取平均；最后，两个卫星的平均轨道根数作差得到平均相对轨道根数的轨道偏心率
Figure GDA0002907323540000134
平均相对轨道根数的轨道倾角
Figure GDA0002907323540000135
平均相对轨道根数的近地点幅角
Figure GDA0002907323540000136
平均相对轨道根数的升交点赤经
Figure GDA0002907323540000137
以平均相对偏心率
Figure GDA0002907323540000138
为例，
Figure GDA0002907323540000139
估计如下：
Figure GDA00029073235400001310
式中eA和eB分别代表A星和B星的轨道偏心率。j表示一个轨道周期内数据序号。N表示一个轨道周期内的数据总个数。假定δt为数据的间隔时长，则Nδt等于一个轨道周期。
平均相对半长轴
Figure GDA0002907323540000141
估计：
根据相对运动学理论知道，由于相对半长轴Δa偏差存在，随着时间增加相对平近点角ΔM也不断增大。相对平近点角ΔM和相对半长轴Δa的关系如下：
Figure GDA0002907323540000142
式中n为A星轨道角速度，a为A星的半长轴。Δt为时间长度。
因此，可通过如下公式间接估计出Δa
Figure GDA0002907323540000143
式中，T为1个轨道周期，
Figure GDA0002907323540000144
为A星的平均半长轴。
为了进一步提高精度，对通过(12)计算的Δaest(t)进行100个数据点取平均，得到平均相对半长轴，如下式所示：
Figure GDA0002907323540000145
平均相对半长轴
Figure GDA0002907323540000146
估计：
估计出
Figure GDA0002907323540000147
计算得到
Figure GDA0002907323540000148
是一个轨道周期内的ΔM的平均值。由于Δa的存在，ΔM是随着Δa呈线性变化的。因此需要对ΔM进行补偿。补偿方法如下
Figure GDA0002907323540000149
(14)仿真中A星初始轨道根数和双星初始平均相对轨道根数偏差可如表1和表2所示：
表1：A星初始轨道根数
参数 a e i Ω ω M
参数值 6451Km 0.68 20deg 91.7deg 130deg 180deg
表2：初始相对轨道根数偏差
Figure GDA0002907323540000151
地面站所提供的A星位置确定精度6Km(3σ)，速度确定精度为3m/s(3σ)。星上相对测量测量时间为10分钟。利用STK中月球HPOP轨道动力学模型得到两颗卫星的轨道根数，通过计算得到相对平均轨道根数作为参考值。采用本发明所提方法的平均相对轨道根数估计如表3所示。通过表3可以看出，采用本发明实施例所发明的方法估计相对平均轨道根数精度很高。
表3：平均相对轨道根数估计精度
Figure GDA0002907323540000152
在本申请所提供的几个实施例中，应该理解到，所揭露的设备和方法，可以通过其它的方式实现。以上所描述的设备实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，如：多个单元或组件可以结合，或可以集成到另一个***，或一些特征可以忽略，或不执行。另外，所显示或讨论的各组成部分相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口，设备或单元的间接耦合或通信连接，可以是电性的、机械的或其它形式的。
上述作为分离部件说明的单元可以是、或也可以不是物理上分开的，作为单元显示的部件可以是、或也可以不是物理单元，即可以位于一个地方，也可以分布到多个网络单元上；可以根据实际的需要选择其中的部分或全部单元来实现本实施例方案的目的。
另外，在本发明各实施例中的各功能单元可以全部集成在一个处理模块中，也可以是各单元分别单独作为一个单元，也可以两个或两个以上单元集成在一个单元中；上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。
本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储于一计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介质包括：移动存储设备、只读存储器、随机存取存储器、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。

Claims (10)

1.一种航天器编队相对轨道确定方法，其特征在于，包括：
测量围绕目标对象飞行的航天器编队中的第一航天器相对于参考对象的第一运动参数；
根据所述第一运动参数，确定出所述第一航天器绕所述目标对象飞行的第一轨道根数；
测量所述航天器编队中第二航天器相对于所述第一航天器的第二运动参数；
根据所述第一轨道根数及所述第二运动参数，确定出所述第二航天器绕所述目标对象飞行的第二轨道根数；
结合所述第一轨道根数及所述第二轨道根数，确定出所述第二航天器相对于所述第一航天器运动的平均相对轨道根数；
所述平均相对轨道根数，包括平均相对偏心率
Figure FDA0002907323530000011
以及平均相对半长轴
Figure FDA0002907323530000012
其中，
Figure FDA0002907323530000013
式中eA和eB分别代表A星和B星的轨道偏心率，j表示一个轨道周期内数据序号，N表示一个轨道周期内的数据总个数，假定δt为数据的间隔时长，则Nδt等于一个轨道周期；
通过相对平近点角ΔM和相对半长轴Δa的关系，获得所述相对半长轴Δa，包括：
Figure FDA0002907323530000014
式中n为A星轨道角速度，a为A星的半长轴，Δt为时间长度，
Figure FDA0002907323530000015
式中，T为1个轨道周期，
Figure FDA0002907323530000016
为A星的平均半长轴，
通过Δaest(t)进行100个数据点取平均，得到平均相对半长轴Δa，
Figure FDA0002907323530000021
其中，Δaest(t)为相对半长轴Δa的估计值。
2.根据权利要求1所述的方法，其特征在于，
所述测量围绕目标对象飞行的航天器编队中的第一航天器相对于参考对象的第一运动参数，包括：
测量围绕月球飞行的航天器编队中的第一航天器相对于地球的第一运动参数。
3.根据权利要求1或2所述的方法，其特征在于，
所述第一运动参数包括：所述第一航天器运行预设时长所经过的位置和对应位置处的速度；
所述根据所述第一运动参数，确定出所述第一航天器绕所述目标对象飞行的第一轨道根数，包括：
根据所述第一航天器运行预设时长所经过的位置和对应位置处的速度，计算出所述第一航天器瞬时轨道根数。
4.根据权利要求3所述的方法，其特征在于，
所述第二运动参数包括：所述第二航天器运行一周相对于所述第一航天器的相对位置处的相对速度；
所述根据所述第一轨道根数及所述第二运动参数，确定出所述第二航天器绕所述目标对象飞行的第二轨道根数，包括：
根据所述第一轨道根数及所述第二运动参数，确定出所述第二航天器运行一周相对于所述参考对象的位置和对应位置处的速度；
根据所述第二航天器运行一周相对于所述参考对象的位置和对应位置处的速度，计算出所述第二航天器瞬时轨道根数。
5.根据权利要求1或2所述的方法，其特征在于，
所述测量所述航天器编队中第二航天器相对于所述第一航天器的第二运动参数，包括：
利用无线电测量所述第二航天器相对于所述第一航天器的相对运动参数，其中，所述相对运动参数包括以下至少之一：相对位置及相对速度；
和/或，
利用图像采集测量所述第二航天器相对于所述第一航天器的相对运动参数，其中，所述相对运动参数包括以下至少之一：相对位置及相对速度。
6.一种航天器编队相对轨道确定装置，其特征在于，包括：
第一测量模块，用于测量围绕目标对象飞行的航天器编队中的第一航天器相对于参考对象的第一运动参数；
第一确定模块，用于根据所述第一运动参数，确定出所述第一航天器绕所述目标对象飞行的第一轨道根数；
第二测量模块，用于测量所述航天器编队中第二航天器相对于所述第一航天器的第二运动参数；
第二确定模块，用于根据所述第一轨道根数及所述第二运动参数，确定出所述第二航天器绕所述目标对象飞行的第二轨道根数；
第三确定模块，用于结合所述第一轨道根数及所述第二轨道根数，确定出所述第二航天器相对于所述第一航天器运动的平均相对轨道根数；
所述平均相对轨道根数，包括平均相对偏心率
Figure FDA0002907323530000031
以及平均相对半长轴
Figure FDA0002907323530000032
其中，
Figure FDA0002907323530000033
式中eA和eB分别代表A星和B星的轨道偏心率，j表示一个轨道周期内数据序号，N表示一个轨道周期内的数据总个数，假定δt为数据的间隔时长，则Nδt等于一个轨道周期；
通过相对平近点角ΔM和相对半长轴Δa的关系，获得所述相对半长轴Δa，包括：
Figure FDA0002907323530000034
式中n为A星轨道角速度，a为A星的半长轴，Δt为时间长度，
Figure FDA0002907323530000035
式中，T为1个轨道周期，
Figure FDA0002907323530000041
为A星的平均半长轴，
通过Δaest(t)进行100个数据点取平均，得到平均相对半长轴Δa，
Figure FDA0002907323530000042
其中，Δaest(t)为相对半长轴Δa的估计值。
7.根据权利要求6所述的装置，其特征在于，
所述第一测量模块，用于测量围绕月球飞行的航天器编队中的第一航天器相对于地球的第一运动参数。
8.根据权利要求6或7所述的装置，其特征在于，
所述第一运动参数包括：所述第一航天器运行预设时长所经过的位置和对应位置处的速度；
所述第一确定模块，用于根据所述第一航天器运行预设时长所经过的位置和对应位置处的速度，计算出所述第一航天器瞬时轨道根数。
9.根据权利要求8所述的装置，其特征在于，
所述第二运动参数包括：所述第二航天器运行一周相对于所述第一航天器的相对位置处的相对速度；
所述第二确定模块，用于根据所述第二运动参数，确定出所述第二航天器运行一周相对于所述参考对象的位置和对应位置处的速度；根据所述第二航天器运行一周相对于所述参考对象的位置和对应位置处的速度，计算出所述第二航天器瞬时轨道根数。
10.根据权利要求6或7所述的装置，其特征在于，
所述第二测量模块，具体用于：
利用无线电测量所述第二航天器相对于所述第一航天器的相对运动参数，其中，所述相对运动参数包括以下至少之一：相对位置及相对速度；
和/或，
利用图像采集测量所述第二航天器相对于所述第一航天器的相对运动参数，其中，所述相对运动参数包括以下至少之一：相对位置及相对速度。
CN201810355799.4A 2018-04-19 2018-04-19 航天器编队相对轨道确定方法及装置 Active CN108490973B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810355799.4A CN108490973B (zh) 2018-04-19 2018-04-19 航天器编队相对轨道确定方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810355799.4A CN108490973B (zh) 2018-04-19 2018-04-19 航天器编队相对轨道确定方法及装置

Publications (2)

Publication Number Publication Date
CN108490973A CN108490973A (zh) 2018-09-04
CN108490973B true CN108490973B (zh) 2021-04-13

Family

ID=63313619

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810355799.4A Active CN108490973B (zh) 2018-04-19 2018-04-19 航天器编队相对轨道确定方法及装置

Country Status (1)

Country Link
CN (1) CN108490973B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109543292B (zh) * 2018-11-21 2022-04-29 哈尔滨工业大学 基于轨道预测的星间链路信号捕获过程的仿真方法
CN109189102B (zh) * 2018-11-23 2021-10-01 上海航天控制技术研究所 一种星上高精度计算双星半长轴偏差的方法
CN110471432B (zh) * 2019-07-04 2020-09-08 中国科学院电子学研究所 一种卫星编队构型的方法、装置及存储介质
CN113204250B (zh) * 2021-04-29 2022-03-08 西安电子科技大学 强动态条件下的卫星编队相对位置鲁棒高精度估计方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004210032A (ja) * 2002-12-27 2004-07-29 Mitsubishi Electric Corp 編隊飛行衛星
CN102322862A (zh) * 2011-06-29 2012-01-18 航天东方红卫星有限公司 一种编队飞行卫星绝对和相对轨道确定方法
CN103954281A (zh) * 2014-04-28 2014-07-30 北京控制工程研究所 一种基于共面星间测向和异面星间测距的星座导航方法
CN104142686A (zh) * 2014-07-16 2014-11-12 北京控制工程研究所 一种卫星自主编队飞行控制方法
JP2015000703A (ja) * 2013-06-18 2015-01-05 三菱電機株式会社 宇宙機の軌道制御方法およびその装置
CN105068546A (zh) * 2015-07-31 2015-11-18 哈尔滨工业大学 一种卫星编队相对轨道自适应神经网络构形包含控制方法
CN105138010A (zh) * 2015-08-31 2015-12-09 哈尔滨工业大学 一种编队卫星分布式有限时间跟踪控制方法
CN105138005A (zh) * 2015-06-16 2015-12-09 北京航空航天大学 一种基于星间距离的相对轨道要素确定方法
CN105737834A (zh) * 2014-12-09 2016-07-06 上海新跃仪表厂 一种基于轨道平根数的相对导航鲁棒滤波方法
CN106767768A (zh) * 2016-12-06 2017-05-31 东南大学 一种双星编队的自主导航方法
CN107450578A (zh) * 2017-07-26 2017-12-08 中国人民解放军空军工程大学 一种绕地球同步轨道卫星分布式共轨飞行的卫星群轨道设计方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016111317A1 (ja) * 2015-01-09 2016-07-14 三菱電機株式会社 軌道制御装置および衛星
US9919813B2 (en) * 2015-04-15 2018-03-20 The United States Of America, As Represented By The Secretary Of The Navy Control system and method for a plane change for satellite operations

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004210032A (ja) * 2002-12-27 2004-07-29 Mitsubishi Electric Corp 編隊飛行衛星
CN102322862A (zh) * 2011-06-29 2012-01-18 航天东方红卫星有限公司 一种编队飞行卫星绝对和相对轨道确定方法
JP2015000703A (ja) * 2013-06-18 2015-01-05 三菱電機株式会社 宇宙機の軌道制御方法およびその装置
CN103954281A (zh) * 2014-04-28 2014-07-30 北京控制工程研究所 一种基于共面星间测向和异面星间测距的星座导航方法
CN104142686A (zh) * 2014-07-16 2014-11-12 北京控制工程研究所 一种卫星自主编队飞行控制方法
CN105737834A (zh) * 2014-12-09 2016-07-06 上海新跃仪表厂 一种基于轨道平根数的相对导航鲁棒滤波方法
CN105138005A (zh) * 2015-06-16 2015-12-09 北京航空航天大学 一种基于星间距离的相对轨道要素确定方法
CN105068546A (zh) * 2015-07-31 2015-11-18 哈尔滨工业大学 一种卫星编队相对轨道自适应神经网络构形包含控制方法
CN105138010A (zh) * 2015-08-31 2015-12-09 哈尔滨工业大学 一种编队卫星分布式有限时间跟踪控制方法
CN106767768A (zh) * 2016-12-06 2017-05-31 东南大学 一种双星编队的自主导航方法
CN107450578A (zh) * 2017-07-26 2017-12-08 中国人民解放军空军工程大学 一种绕地球同步轨道卫星分布式共轨飞行的卫星群轨道设计方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Autonomous Orbit Determination and Observability Analysis for Formation Satellites;Yangwei Ou等;《2016 35th Chinese Control Conference》;20160829;全文 *
Relative States Autonomous Determination of Satellites Formation Flying;Xi-bin Cao,等;《2006 International Conference on Machine Learning and Cybernetics》;20091231;全文 *
卫星编队相对轨道的分布式估计方法;雪丹等;《吉林大学学报（工学报）》;20061231;第36卷(第6期);全文 *
基于模型预测的卫星编队队形机动控制;吴宝林等;《吉林大学学报(工学版)》;20071231;第219-222段及图1 *
基于轨道根数差的卫星编队自主定轨研究;雪丹等;《***仿真学报》;20061231;第18卷(第10期);全文 *

Also Published As

Publication number Publication date
CN108490973A (zh) 2018-09-04

Similar Documents

Publication Publication Date Title
CN108490973B (zh) 航天器编队相对轨道确定方法及装置
Vetter Fifty years of orbit determination
CN109931955B (zh) 基于状态相关李群滤波的捷联惯性导航***初始对准方法
CN104713559B (zh) 一种高精度sins模拟器的设计方法
Abdelrahman et al. Sigma-point Kalman filtering for spacecraft attitude and rate estimation using magnetometer measurements
RU2318188C1 (ru) Способ автономной навигации и ориентации космических аппаратов
CN104048664A (zh) 一种导航卫星星座自主定轨的方法
CN112129300A (zh) 位置间动力学约束的低轨卫星星载gnss精密定轨方法及***
CN113091731A (zh) 一种基于恒星视线相对论效应的航天器自主导航方法
CN116105730A (zh) 基于合作目标卫星甚短弧观测的仅测角光学组合导航方法
Somodi et al. Application of numerical integration techniques for orbit determination of state-of-the-art LEO satellites
Liu et al. Autonomous orbit determination and timekeeping in lunar distant retrograde orbits by observing X‐ray pulsars
CN106643726B (zh) 一种统一惯性导航解算方法
Wesam et al. Kalman filter implementation for small satellites using constraint GPS data
Kramlikh et al. Estimation method for nanosatellite orbital parameters in case of abnormal operation of navigation equipment
Leonard et al. Liaison-supplemented navigation for geosynchronous and lunar l1 orbiters
İpek Satellite orbit estimation using kalman filters
Olson Sequential estimation methods for small body optical navigation
Um et al. GPS attitude determination for the SOAR experiment
Searcy Observability-enhanced dual-filter design for attitude estimation with minimum observations
Tomasicchio et al. Multi-Sensor Fusion for Improved Navigation in Lunar Landing Missions
Zuiani et al. Orbit determination and control for the European Student Moon Orbiter
CN116552817B (zh) 优化卫星轨道控制的方法、设备和存储介质
CN115096317B (zh) 一种地月空间dro航天器编队相对导航方法和***
Amert et al. Hardware Demonstration and Improvements of the Stellar Positioning System

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