CN108038290A - 一种多重约束下的多维度小行星探测目标筛选方法 - Google Patents

一种多重约束下的多维度小行星探测目标筛选方法 Download PDF

Info

Publication number
CN108038290A
CN108038290A CN201711264888.XA CN201711264888A CN108038290A CN 108038290 A CN108038290 A CN 108038290A CN 201711264888 A CN201711264888 A CN 201711264888A CN 108038290 A CN108038290 A CN 108038290A
Authority
CN
China
Prior art keywords
asteroid
constraint
detector
detection
target
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.)
Pending
Application number
CN201711264888.XA
Other languages
English (en)
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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201711264888.XA priority Critical patent/CN108038290A/zh
Publication of CN108038290A publication Critical patent/CN108038290A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开的一种多重约束下的多维度小行星探测目标筛选方法，属于航空航天技术领域。本发明根据探测任务的多重约束强弱情况对约束进行分类和排序，确定筛选顺序；根据探测器在轨时间寿命，测控站所支持的观测通信距离，确定探测时间范围约束和探测器与地球的测控通信约束，并对小行星目标进行初选；根据探测器成像约束及小行星物理特性对目标进行进一步筛选；基于日‑地‑月系统转移轨道两点边值问题求解，根据速度增量约束进一步筛选目标小行星，确定筛选出的目标小行星数量，并完成相应的小行星探测任务。本发明在考虑多重约束条件下实现小行星探测目标筛选，并通过对约束强弱等级和筛选难易程度进行判断，提高筛选速度和筛选效率。

Description

一种多重约束下的多维度小行星探测目标筛选方法
技术领域
本发明涉及一种多重约束下的多维度小行星探测目标筛选方法，尤其适用于从日地平衡点轨道出发、满足多重约束下的多维度小行星探测目标筛选方法，属于航空航天技术领域。
背景技术
自20世纪后期，美国国家航空航天局(NASA)Galileo号探测器首次开展小行星探测活动以来，小行星探测成为深空探测的热点。而小行星探测的目标选择是开展小行星探测的难点，如何更加快捷方便并且准确有效地找到符合探测任务要求的小行星目标从一定程度上决定了任务的可实现性及任务的价值。
在已发展的关于航天器探测小行星的目标选择方法中在先技术[1](参见：焦维新，钟俊.近地小行星探测目标选择[J].地球物理学报，2016，59(11)：3955-3959)给出基于小行星交会所需能量大小、小行星半径大小、小行星自旋周期、小行星可能撞击地球的概率等方面的目标小行星选择方法。这种方法需要对所有近地小行星进行速度增量计算，目标选择耗时多，没有考虑实际任务中测控和任务时间等约束情况，且仅适用于从地球出发的小行星探测任务。此外，这种方法需要进行大量从近地轨道到小行星的转移速度增量计算，效率低，速度慢。
发明内容
本发明公开的一种多重约束下的多维度小行星探测目标筛选方法要解决的技术问题为，在考虑多重约束条件下实现小行星探测目标筛选，并通过对约束强弱等级和筛选难易程度进行判断，提高筛选速度和筛选效率。所述的多重约束包括探测时间范围约束、探测器与地球的测控通信约束、探测器成像约束、探测器速度增量约束。
本发明的目的是通过下述技术方案实现的：
本发明公开的一种多重约束下的多维度小行星探测目标筛选方法，根据探测任务的多重约束强弱情况对约束进行分类和排序，确定筛选顺序。根据探测器在轨时间寿命，测控站所支持的观测通信距离，确定探测时间范围约束和探测器与地球的测控通信约束。根据确定的探测时间范围约束和探测器与地球的测控通信约束，对小行星目标进行初选。根据探测器成像约束及小行星物理特性对目标进行进一步筛选。基于日-地-月系统转移轨道两点边值问题求解，根据速度增量约束进一步筛选目标小行星，并确定筛选出的目标小行星数量。根据筛选出飞越探测的目标小行星及目标小行星数量，完成相应的小行星探测任务。
本发明公开的一种多重约束下的多维度小行星探测目标筛选方法，包括如下步骤：
步骤一：根据探测任务的多重约束强弱情况对约束进行分类和排序，确定筛选顺序。
步骤一优选如下方法实现：
根据探测器实际情况和任务要求，给出探测任务存在的约束类型和约束量级，并对多重约束强弱程度进行分类和排序。所述的多重约束包括探测时间范围约束、探测器与地球的测控通信约束、探测器成像约束、探测器速度增量约束。其中探测时间范围约束为强约束，决定任务执行的时间；其次是探测器与地球的测控通信约束，仅与探测器自身能力有关；探测器观测条件约束和探测器速度增量约束均与目标小行星状态有关。从计算复杂度分析，探测器速度增量约束最复杂，因此按照探测时间范围约束、探测器与地球的测控通信约束、探测器成像约束顺序进行筛选，最后以探测器速度增量为约束对目标进行筛选。
步骤二：确定探测时间范围约束和探测器与地球的测控距离约束。
给定探测器在任务轨道运行所携带的燃料重量，评估探测器寿命，从而选定小行星飞越探测时间范围，保证探测的可靠性。
根据探测器测控站的天线长度、所采用的码速率确定目标探测器飞越小行星时与地球的最大距离约束。
步骤三：根据步骤二确定的探测时间范围约束和探测器与地球的测控距离约束，对小行星目标进行初选。
近地小行星分为四种类型，分别为Amor型、Aten型、Apollo型和Atira型。采用小行星定轨星历，考虑的各大天体摄动，对轨道进行递推，搜索满足探测时间范围约束中与地球最近距离小于测控距离约束的目标小行星，轨道递推的动力学模型为
其中：r为探测器实时位置，μi为各摄动天体的引力常数，χi＝r-ri，ri为轨道积分过程中摄动天体实时的位置矢量，由小行星定轨星历提供，积分终止时间为时间范围约束的终点。轨道递推后，筛选得到在满足时间约束的区间内小行星距离地球距离满足测控通信约束的目标小行星。
步骤三所述的小行星定轨星历优选DE421行星星历，递推算法优选龙格库塔78法。
步骤四：根据探测器成像约束对目标进行进一步筛选。
根据探测器所携带的相机配置及成像质量要求，确定探测小行星的最小直径约束。
小行星的直径计算由公式(2)给出
D＝1329×10-H/5pv -1/2 (2)
其中，D为小行星直径，H为绝对星等，pv为反照率。所述的绝对星等H、反照率pv通过数据库确定。
根据小行星最小直径约束，从步骤三筛选出的目标星中进一步筛选出少量符合要求的小行星。
步骤五：根据速度增量约束进一步筛选目标小行星，并确定筛选出的目标小行星数量。
给定探测器所携带的剩余推进剂重量，计算估计能够用于转移的总速度增量Δv大小。根据速度增量约束对步骤四筛选出的剩余小行星的飞越探测机会以及转移所需总速度增量Δv进行计算分析，最后筛选出飞越探测的目标小行星，并确定筛选出的目标小行星数量。
步骤五所述的计算估计能够用于转移的总速度增量Δv大小具体实现方法如下：基于日-地-月以及探测器四体星历动力学模型，建立探测器飞越探测目标小行星时与小行星的最近距离、切向脉冲速度增量以及飞行时间之间的关系。通过探测时间范围约束和探测器与地球测控距离约束计算用于转移的总速度增量Δv大小。
还包括步骤六：根据步骤五筛选出的飞越探测的目标小行星及目标小行星数量，完成相应的小行星探测任务。
有益效果：
1、在已有的探测小行星目标选择方法中，都需要进行大量速度增量计算，造成筛选速度慢、效率低。本发明公开的一种多重约束下的多维度小行星探测目标筛选方法，根据探测任务的多重约束强弱情况对约束进行分类和排序，确定筛选顺序，按照确定的筛选顺序对小行星探测目标进行筛选，提高筛选速度和筛选效率。
2、本发明公开的一种多重约束下的多维度小行星探测目标筛选方法，根据探测器实际情况和任务要求，给出探测任务存在的约束类型和约束量级，并对多重约束强弱程度进行分类和排序；所述的多重约束包括探测时间范围约束、探测器与地球的测控通信约束、探测器成像约束、探测器速度增量约束，因此，适用于具有探测时间范围约束、探测器与地球的测控通信约束、探测器成像约束、探测器速度增量约束的小行星探测目标筛选问题。
附图说明
图1为本发明公开的一种多重约束下的多维度小行星探测目标筛选方法流程图。
图2为本发明公开的一种多重约束下的多维度小行星探测目标筛选方法转移轨道图。
具体实施方式
为了更好地说明本发明的目的和优点，下面结合具体实施示例对本发明做出详细解释。
实施例1：
如图1所示，本实施例公开的一种多重约束下的多维度小行星探测目标筛选方法，具体实现步骤如下：
步骤一：根据探测任务的多重约束强弱情况对约束进行分类和排序，确定筛选顺序。
根据探测器实际情况和任务要求，给出探测任务存在的约束类型和约束量级，并对多重约束强弱程度进行分类和排序。在本实施例中，多重约束包括探测时间范围约束、探测器与地球的测控通信约束、探测器成像约束、探测器速度增量约束。其中探测时间范围约束为强约束，决定任务执行的时间；其次是探测器与地球的测控通信约束，仅与探测器自身能力有关；探测器观测条件约束和探测器速度增量约束均与目标小行星状态有关。从计算复杂度分析，探测器速度增量约束最复杂，因此按照探测时间范围约束、探测器与地球的测控通信约束、探测器成像约束顺序进行筛选，最后以探测器速度增量为约束对目标进行筛选。
步骤二：确定探测时间范围约束和探测器与地球的测控距离约束。
对日-地拉格朗日L2点Lissajous轨道运行的探测器寿命进行评估，假设探测器初始时刻为2011年10月，探测器寿命小于2年，选定小行星飞越探测时间范围2011年10月至2013年6月。同时通过监测探测器在日-地拉格朗日L2点附近的Lissajous轨道(距离地球约170万千米)飞行的状态，根据探测器与地面测控站数据传输所采用的码速率以及测控站天线长度，确定能够保证探测器与测控站通信可靠性的小行星距离地球的最远距离rmax＜1500万km，因此目标小行星与地球的距离约束为r＜rmax。其中，r为探测器飞越探测目标小行星时与地球的距离。
步骤三：根据步骤二确定的探测时间范围约束和探测器与地球的测控距离约束，对小行星目标进行初选。
对目前发现的一万多颗近地小行星进行轨道递推，搜索满足在探测时间范围约束[t1,t2]时间段内与地球最近距离小于rmax的小行星，其中，rmax为步骤二中确定的小行星探测距离约束。小行星的星历采用美国喷气动力实验室(JPL)公布的观测定轨星历。轨道递推的动力学模型为
其中r为探测器实时位置，μi为各摄动天体的引力常数，χi＝r-ri，ri为轨道积分过程中摄动天体实时的位置矢量，由DE421星历得到，积分终止时间为时间范围约束的终点。递推算法选择龙格库塔78法，递推时间为DE421星历的历元时刻至2013年6月。
轨道递推后，可以得到满足探测时间约束和测控通信约束的目标小行星。这一步筛选出的小行星数量为100颗左右。
步骤四：根据探测器成像约束对目标进行进一步筛选。
根据探测器任务要求，综合考虑星载相机的性能指标以及成像质量要求，确定飞越探测的目标小行星的最小直径为Dmin＝1km，确定小行星直径约束满足约束的目标共3颗为Toutatis，1998QE2，2005NZ6。
步骤五：根据速度增量约束进一步筛选目标小行星，并确定筛选出的目标小行星数量。
给定探测器所携带的剩余推进剂重量，计算估计可用于转移的总速度增量大小Δv＜120m/s。根据该约束对优选出的3颗小行星的飞越探测机会以及转移所需总速度增量进行计算分析，最后筛选和确定飞越探测的目标星。
步骤五所述的计算估计能够用于转移的总速度增量Δv大小具体实现方法如下：基于日-地-月以及探测器四体星历动力学模型，求解两点边值问题，计算得出探测器飞越探测小行星的飞越机会以及转移所需的总速度增量。在t0 -时刻探测器在Lissajous轨道上的运行状态为x0＝[r0,v0]T，在t0 +时刻沿轨道切向施加脉冲速度增量Δv，使得探测器从Lissajous轨道出发进行轨道转移。基于日-地-月-卫四体星历动力学模型，建立探测器飞越探测目标小行星时与小行星的最近距离xf、切向脉冲速度增量大小Δv以及飞行时间tf之间的关系。通过探测时间范围约束和探测器与地球测控距离约束计算用于转移的总速度增量Δv大小，搜索方法为遍历算法或优化算法。优选出的3颗小行星的飞越探测机会以及转移所需总速度增量见表1。从日地L2点Lissajous轨道到目标小行星转移轨道图见图2。
表1优选出的目标小行星转移总速度增量表
备选目标中Toutatis飞越探测所需速度增量为Δv＝87.4m/s，且飞越时间在2年之内，满足所有约束，可作为小行星探测任务的备选目标。
步骤六：根据步骤五筛选出的飞越探测的目标小行星及目标小行星数量，完成此次小行星探测任务。
以上所述的具体描述，对发明的目的、技术方案和优点进行了进一步详细说明。需要理解的是，以上所述仅为本发明的具体实施示例，用于解释本发明，并不用于限定本发明的保护范围。凡在本发明的精神和原则之内所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种多重约束下的多维度小行星探测目标筛选方法，其特征在于：包括如下步骤，
步骤一：根据探测任务的多重约束强弱情况对约束进行分类和排序，确定筛选顺序；
步骤二：确定探测时间范围约束和探测器与地球的测控距离约束；
步骤三：根据步骤二确定的探测时间范围约束和探测器与地球的测控距离约束，对小行星目标进行初选；
步骤四：根据探测器成像约束对目标进行进一步筛选；
步骤五：根据速度增量约束进一步筛选目标小行星，并确定筛选出的目标小行星数量。
2.如权利要求1所述的一种多重约束下的多维度小行星探测目标筛选方法，其特征在于：步骤一具体实现方法为，
根据探测器实际情况和任务要求，给出探测任务存在的约束类型和约束量级，并对多重约束强弱程度进行分类和排序；所述的多重约束包括探测时间范围约束、探测器与地球的测控通信约束、探测器成像约束、探测器速度增量约束；其中探测时间范围约束为强约束，决定任务执行的时间；其次是探测器与地球的测控通信约束，仅与探测器自身能力有关；探测器观测条件约束和探测器速度增量约束均与目标小行星状态有关；从计算复杂度分析，探测器速度增量约束最复杂，因此按照探测时间范围约束、探测器与地球的测控通信约束、探测器成像约束顺序进行筛选，最后以探测器速度增量为约束对目标进行筛选。
3.如权利要求2所述的一种多重约束下的多维度小行星探测目标筛选方法，其特征在于：步骤二具体实现方法为，
给定探测器在任务轨道运行所携带的燃料重量，评估探测器寿命，从而选定小行星飞越探测时间范围，保证探测的可靠性；
根据探测器测控站的天线长度、所采用的码速率确定目标探测器飞越小行星时与地球的最大距离约束。
4.如权利要求3所述的一种多重约束下的多维度小行星探测目标筛选方法，其特征在于：步骤三具体实现方法为，
近地小行星分为四种类型，分别为Amor型、Aten型、Apollo型和Atira型；采用小行星定轨星历，考虑的各大天体摄动，对轨道进行递推，搜索满足探测时间范围约束中与地球最近距离小于测控距离约束的目标小行星，轨道递推的动力学模型为
<mrow> <mover> <mi>r</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>=</mo> <mo>-</mo> <mfrac> <mi>r</mi> <msup> <mi>r</mi> <mn>3</mn> </msup> </mfrac> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>11</mn> </munderover> <msub> <mi>&amp;mu;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mfrac> <msub> <mi>&amp;chi;</mi> <mi>i</mi> </msub> <msubsup> <mi>&amp;chi;</mi> <mi>i</mi> <mn>3</mn> </msubsup> </mfrac> <mo>+</mo> <mfrac> <msub> <mi>r</mi> <mi>i</mi> </msub> <msubsup> <mi>r</mi> <mi>i</mi> <mn>3</mn> </msubsup> </mfrac> <mo>)</mo> </mrow> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mn>11</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
其中：r为探测器实时位置，μi为各摄动天体的引力常数，χi＝r-ri，ri为轨道积分过程中摄动天体实时的位置矢量，由小行星定轨星历提供，积分终止时间为时间范围约束的终点；轨道递推后，筛选得到在满足时间约束的区间内小行星距离地球距离满足测控通信约束的目标小行星。
5.如权利要求4所述的一种多重约束下的多维度小行星探测目标筛选方法，其特征在于：步骤四具体实现方法为，
根据探测器所携带的相机配置及成像质量要求，确定探测小行星的最小直径约束；
小行星的直径计算由公式(2)给出
D＝1329×10-H/5pv -1/2 (2)
其中，D为小行星直径，H为绝对星等，pv为反照率；所述的绝对星等H、反照率pv通过数据库确定；
根据小行星最小直径约束，从步骤三筛选出的目标星中进一步筛选出少量符合要求的小行星。
6.如权利要求5所述的一种多重约束下的多维度小行星探测目标筛选方法，其特征在于：步骤五具体实现方法为，
给定探测器所携带的剩余推进剂重量，计算估计能够用于转移的总速度增量Δv大小；根据速度增量约束对步骤四筛选出的剩余小行星的飞越探测机会以及转移所需总速度增量Δv进行计算分析，最后筛选出飞越探测的目标小行星，并确定筛选出的目标小行星数量。
7.如权利要求6所述的一种多重约束下的多维度小行星探测目标筛选方法，其特征在于：步骤五所述的计算估计能够用于转移的总速度增量Δv大小具体实现方法如下，
基于日-地-月以及探测器四体星历动力学模型，建立探测器飞越探测目标小行星时与小行星的最近距离、切向脉冲速度增量以及飞行时间之间的关系；通过探测时间范围约束和探测器与地球测控距离约束计算用于转移的总速度增量Δv大小。
8.如权利要求7所述的一种多重约束下的多维度小行星探测目标筛选方法，其特征在于：步骤三所述的小行星定轨星历选DE421行星星历，递推算法选龙格库塔78法。
9.如权利要求1、2、3、4、5、6、7或8所述的一种多重约束下的多维度小行星探测目标筛选方法，其特征在于：还包括步骤六，根据步骤五筛选出的飞越探测的目标小行星及目标小行星数量，完成相应的小行星探测任务。
CN201711264888.XA 2017-12-05 2017-12-05 一种多重约束下的多维度小行星探测目标筛选方法 Pending CN108038290A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711264888.XA CN108038290A (zh) 2017-12-05 2017-12-05 一种多重约束下的多维度小行星探测目标筛选方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711264888.XA CN108038290A (zh) 2017-12-05 2017-12-05 一种多重约束下的多维度小行星探测目标筛选方法

Publications (1)

Publication Number Publication Date
CN108038290A true CN108038290A (zh) 2018-05-15

Family

ID=62095300

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711264888.XA Pending CN108038290A (zh) 2017-12-05 2017-12-05 一种多重约束下的多维度小行星探测目标筛选方法

Country Status (1)

Country Link
CN (1) CN108038290A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112861373A (zh) * 2021-03-04 2021-05-28 中国科学院国家天文台 近地小行星撞击轨道的生成方法和装置
CN113536536A (zh) * 2021-06-08 2021-10-22 北京航天飞行控制中心 一种多天体飞越探测轨道设计方法
CN114997043A (zh) * 2022-05-18 2022-09-02 中国科学院国家空间科学中心 一种撞击近地小行星的方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103158891A (zh) * 2013-03-04 2013-06-19 北京理工大学 一种从动平衡点轨道出发飞越小天体的目标选择方法
CN105083596A (zh) * 2015-09-21 2015-11-25 上海卫星工程研究所 近地小行星探测目标星选择方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103158891A (zh) * 2013-03-04 2013-06-19 北京理工大学 一种从动平衡点轨道出发飞越小天体的目标选择方法
CN105083596A (zh) * 2015-09-21 2015-11-25 上海卫星工程研究所 近地小行星探测目标星选择方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CUI PINGYUAN 等: "Target selection and transfer trajectories design for exploring asteroid mission", 《SCIENCE CHINA TECHNOLOGICAL SCIENCES》 *
乔栋 等: "嫦娥二号卫星飞越探测小行星的目标选择", 《中国科学：技术科学》 *
熊旭 等: "小行星采样返回探测任务发射机会搜索", 《中国宇航学会深空探测技术专业委员会第十届学术年会论文集》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112861373A (zh) * 2021-03-04 2021-05-28 中国科学院国家天文台 近地小行星撞击轨道的生成方法和装置
CN112861373B (zh) * 2021-03-04 2023-09-01 中国科学院国家天文台 近地小行星撞击轨道的生成方法和装置
CN113536536A (zh) * 2021-06-08 2021-10-22 北京航天飞行控制中心 一种多天体飞越探测轨道设计方法
CN113536536B (zh) * 2021-06-08 2023-10-24 北京航天飞行控制中心 一种多天体飞越探测轨道设计方法
CN114997043A (zh) * 2022-05-18 2022-09-02 中国科学院国家空间科学中心 一种撞击近地小行星的方法

Similar Documents

Publication Publication Date Title
Sibeck et al. Imaging plasma density structures in the soft X-rays generated by solar wind charge exchange with neutrals
CN108038290A (zh) 一种多重约束下的多维度小行星探测目标筛选方法
CN111680455B (zh) 基于搭载形式的撞击探测轨道设计方法和系统
CN104792299B (zh) 一种基于测角数据的小行星轨道识别方法
Pu et al. Optimal small satellite orbit design based on robust multi-objective optimization method
CN109117543B (zh) 载人航天器对近地小行星探测并返回的轨道设计方法
CN115314101A (zh) 一种基于并行计算的低轨通信卫星星座快速建模方法
CN108628345B (zh) 一种电磁航天器编队悬停协同控制方法及系统
CN112800532B (zh) 一种网格化发射弹道规划方法
Barbee et al. Methodology and results of the near-earth object (NEO) human space flight (HSF) accessible targets study (NHATS)
CN104864875A (zh) 一种基于非线性h∞滤波的航天器自主定位方法
Vassev et al. Autonomy requirements engineering: a case study on the BepiColombo mission
Fraser Adaptive extended Kalman filtering strategies for autonomous relative navigation of formation flying spacecraft
Tolson et al. Atmospheric modeling using accelerometer data during Mars Atmosphere and Volatile Evolution (MAVEN) flight operations
Bellome Trajectory design of multi-target missions via graph transcription and dynamic programming.
Luyi et al. Fast solution to the free return orbit's reachable domain of the manned lunar mission by deep neural network
CN116151052B (zh) 某型卫星大幅变轨丢失重捕方法及其装置、电子设备
Chu et al. Mars base technology program overview
CN113536536B (zh) 一种多天体飞越探测轨道设计方法
Phillips et al. Confirmation of the April Rho Cygnids (ARC, IAU# $348)
CN109211225A (zh) 获取大椭圆轨道空间物体剩余轨道寿命的方法、系统及设备
Smith CubeSat De-Orbit Point Targeting using Drag Modulations
Hughes Modeling the Effects of $\mathrm {J} _ {2} $ Perturbation on Orbital Debris and Satellite Motion
Ruan et al. Preliminary Safety Analysis of Megaconstellations in Low Earth Orbit: Assessing Short-Term and Long-Term Collision Risks
Zhu et al. Autonomous navigation algorithm for spacecraft based on sun-earth-moon astronomical information

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20180515