CN112613223A - 基于径向基神经网络的月球引力辅助逃逸轨道递推方法 - Google Patents
基于径向基神经网络的月球引力辅助逃逸轨道递推方法 Download PDFInfo
- Publication number
- CN112613223A CN112613223A CN202011326683.1A CN202011326683A CN112613223A CN 112613223 A CN112613223 A CN 112613223A CN 202011326683 A CN202011326683 A CN 202011326683A CN 112613223 A CN112613223 A CN 112613223A
- Authority
- CN
- China
- Prior art keywords
- detector
- neural network
- gravity
- radial basis
- basis function
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000013528 artificial neural network Methods 0.000 title claims abstract description 47
- 230000005484 gravity Effects 0.000 claims abstract description 56
- 230000006870 function Effects 0.000 claims abstract description 47
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 238000012549 training Methods 0.000 claims abstract description 17
- 230000010354 integration Effects 0.000 claims abstract description 10
- 208000037280 Trisomy Diseases 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000003062 neural network model Methods 0.000 claims description 4
- 230000003631 expected effect Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 2
- 230000000630 rising effect Effects 0.000 claims description 2
- 210000002569 neuron Anatomy 0.000 abstract description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 5
- 238000013461 design Methods 0.000 description 5
- 239000000523 sample Substances 0.000 description 5
- 230000001174 ascending effect Effects 0.000 description 2
- 238000012938 design process Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- 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
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于径向基神经网络的月球引力辅助逃逸轨道递推方法，包括步骤如下：设置探测器初始状态参数；计算探测器在近心点处的状态及引力辅助前的轨道根数；通过圆形限制性三体问题动力学模型进行轨道递推，获得探测器在引力辅助之前的轨道根数及引力辅助之后的特征能量值；学习解析的圆锥曲线拼接法与三体问题结果之间的误差关系；将误差值与圆锥曲线拼接法计算结果进行代数相加运算，从而替代在轨道设计时三体问题所需的数值积分。本发明采用径向基神经网络对均匀采样的部分初始相空间进行网络训练，通过对训练时间、训练样本数及神经元个数的折衷权衡选择合适的神经网络结构；计算三体问题引力辅助时，计算效率高且精度高。
Description
技术领域
本发明属于航空航天技术领域，具体指代一种基于径向基神经网络的月球引力辅助逃逸轨道递推方法。
背景技术
利用月球引力辅助逃逸地月系统是低能耗深空探测的重要途径之一。对于地月逃逸问题而言，由于地月引力共同的作用，在轨道初始设计阶段若采用经典的圆锥曲线拼接方法设计月球引力辅助逃逸轨道则与高精度模型相比误差较大。因此，采用限制性三体问题模型来设计月球引力辅助十分必要。然而，受复杂、强非线性动力学环境的影响，三体问题没有解析解，需要大量的时间来进行数值积分计算，严重影响了复杂地月系统逃逸任务设计分析的效率。如何快速计算三体问题下的轨道也是当前科技人员关注的热点问题之一。
已发展的基于人工神经网络的三体问题快速计算方法中，在先文献De Smet S.Onthe Design of Solar Gravity Driven Planetocentric Transfers Using ArtificialNeural Networks[D].University of Colorado at Boulder,2018中公开了一种基于人工神经网络系统探索太阳引力驱动火星系统轨道转移问题的方法，可以高效地解决轨道设计过程中面临的大量数值积分问题。但是该方法面向的是从抛物线轨道捕获至火星环绕轨道的任务，其任务设计过程与地月系统逃逸问题不同，无法应用于月球引力辅助逃逸轨道递推计算中。
在先文献Breen P G,Foley C N,Boekholt T,et al.Newton vs the machine:solving the chaotic three-body problem using deep neural networks[J].arXivpreprint arXiv:1910.07291,2019中通过利用神经网络能够在固定计算成本内提供准确的三体问题的解并且速度提升1亿倍，但是其研究问题被限制为三个质量相等、初始速度为零的粒子在一个平面上的引力问题。三体问题并未被推广至一般情况，因此，无法在引力辅助问题中得到使用。
发明内容
针对于上述现有技术的不足，本发明的目的在于提供一种基于径向基神经网络的月球引力辅助逃逸轨道递推方法,以解决现有技术中难以用很短的时间实现限制性三体问题模型下月球引力辅助逃逸轨道快递精确递推计算的问题。
为达到上述目的，本发明采用的技术方案如下：
本发明的一种基于径向基神经网络的月球引力辅助逃逸轨道递推方法，包括步骤如下：
1)设置探测器初始状态参数；
2)利用解析的圆锥曲线拼接法计算探测器在近心点处的状态及引力辅助前的轨道根数；
3)在相同的近心点处，通过圆形限制性三体问题动力学模型进行轨道递推，从而获得探测器在引力辅助之前的轨道根数及引力辅助之后的特征能量C3值；
4)基于径向基神经网络学习解析的圆锥曲线拼接法与三体问题结果之间的误差关系；
5)基于上述步骤4)得到的误差值，将其与圆锥曲线拼接法计算结果进行代数相加运算，从而替代在轨道设计时三体问题所需的数值积分，提高计算效率。
进一步地，所述步骤1)具体包括：设定探测器的初始特征能量C3值；设定探测器初始阶段状态约束条件参数：轨道倾角i、近星距Rp与远星距Ra；两个主天体质量M1、M2，质量参数比μ＝M2/(M1+M2)，次天体影响球半径rSOI及两个主天体之间的距离REM。
进一步地，所述步骤2)具体包括：依次建立以两个主天体质心为原点的会合坐标系，再分别建立以主体为中心的惯性坐标系和以次体为中心的引力辅助坐标系，其中会合坐标系的坐标轴原点为两个主天体的质心，x轴由M1指向M2，z轴与地月系统的角动量方向一致，y轴满足右手定则；当t＝0时的会合坐标系为该地月系统的惯性坐标系；在引力辅助坐标系中，rp和vp分别描述在该坐标系下，探测器在近心点处的位置和速度；已知探测器的初始状态近星距Rp、远星距Ra和轨道倾角i时，得到此时轨道的半长轴，在主体引力作用下，利用二体能量积分公式计算得出探测器经过影响球边界时的速度大小，从而根据探测器的角动量计算出速度的侧向分量大小，进而得到惯性坐标系下探测器进入影响球时的速度矢量
且进入影响球时的位置矢量为
根据圆锥曲线模型的假设，引力辅助发生时得到行星的速度矢量是一个常矢量，表示为
从而得到引力辅助坐标系下探测器进入影响球时的剩余速度矢量；在已知引力辅助之后探测器具备预想能量C3值的前提下，根据式(1)计算出探测器飞出影响球时的速度大小：
从而得到引力辅助坐标系下，探测器在近心点处时的位置和速度矢量：
由此，根据已求得的进入影响球的位置和速度矢量计算出引力辅助前的轨道根数。
进一步地，所述步骤3)具体包括：
建立三体问题动力学模型，其运动方程在会合坐标系下表示为：
式中，r1和r2分别为探测器与天体M1和M2之间的距离，质量参数为μ，该数学模型的积分初始量为近心点处的位置和速度；
基于三体问题动力学模型，通过数值积分方法在时间上分别向前向后进行轨道递推获得引力辅助前的轨道根数及引力辅助之后的特征能量C3数据分布。
进一步地，所述步骤4)具体包括：
解析的圆锥曲线拼接法与三体问题之间的误差计算如下：
式中，a表示半长轴，e是偏心率，i是轨道倾角，Ω是升交点赤经，w是近地点幅角，C3表征三体问题下，特征能量的计算公式；
根据上述已经获得的误差，利用径向基(RBF)神经网络对其进行学习；确定神经网络的结构，径向基神经网络的网络层数为3层，分别为输入层、隐含层、输出层，输入层到隐含层的权值由输入进来的信号样本来确定，隐含层到输出层的权值通过求解线性方程组来获得；
在采用径向基神经网络预测时，使用MATLAB神经网络工具箱中的newrbe函数，该函数的调用格式为：net＝newrbe(P,T,spread)，P和T分别表示神经网络的输入向量和输出向量，spread表征径向基函数的分布密度，对spread进行不同数值的尝试，以确定一个最优值；spread的默认值为1，通过对均匀选取的训练样本进行模型的训练，然后通过调用sim函数测试网络性能；如果实际输出结果与真实值的误差相差大，则通过改变网络结构、数据样本和spread参数值进行调整，最终获得理想的网络模型；根据预期效果，在训练过程中需要对spread不断地修改与测试，最终确定spread值。
进一步地，所述步骤5)具体包括：根据上述已经获得的径向基神经网络的模型，在给定探测器的初始轨道倾角和远星距的情况下，带入径向基神经网络模型预测出两个方法的误差值，然后将预测得出的误差值与圆锥曲线拼接方法得到的结果进行代数相加从而获得三体问题下轨道高精度近似解。
本发明的有益效果：
本发明适用于三体问题引力辅助的快速计算；采用径向基神经网络对均匀采样的部分初始相空间进行网络训练，通过对训练时间、训练样本数及神经元个数的折衷权衡选择合适的神经网络结构；计算三体问题引力辅助时，计算效率高且精度高。
附图说明
图1为本发明计算方法流程示意图。
图2a表示会合坐标系示意图。
图2b表示惯性坐标系示意图。
图2c表示引力辅助坐标系示意图。
具体实施方式
为了便于本领域技术人员的理解，下面结合实施例与附图对本发明作进一步的说明，实施方式提及的内容并非对本发明的限定。
参照图1所示，本发明的一种基于径向基神经网络的月球引力辅助逃逸轨道递推方法，包括步骤如下：
1)设置探测器初始状态参数；
2)利用解析的圆锥曲线拼接法计算探测器在近心点处的状态及引力辅助前的轨道根数；
3)在相同的近心点处，通过圆形限制性三体问题动力学模型进行轨道递推，从而获得探测器在引力辅助之前的轨道根数及引力辅助之后的特征能量C3值；
4)基于径向基神经网络学习解析的圆锥曲线拼接法与三体问题结果之间的误差关系；
5)基于上述径向基神经网络模型预测得到的误差值，将其与圆锥曲线拼接法计算结果进行代数相加运算，从而替代在轨道设计时三体问题所需的数值积分，提高计算效率。
以地月引力辅助逃逸轨道递推为例：
1、设定探测器的初始特征能量C3＝0；探测器初始阶段状态约束条件参数轨道倾角i∈[0°,180°]、近星距Rp＝10000km与远星距Ra∈[400000,600000]km；M1代表地球的质量，M2代表月球的质量，质量参数比计算公式为μ＝M2/(M1+M2)，次天体影响球半径rSOI＝(μ/(1-μ))2/5及地月之间的距离为REM。
2、依次建立以地月质心为原点的会合坐标系，再分别建立以地球为中心的惯性坐标系和以月球为中心的引力辅助坐标系，如图2a-图2c所示，其中会合坐标系的坐标轴原点为地月的质心，x轴由M1指向M2，z轴与地月系统的角动量方向一致，y轴满足右手定则；当t＝0时的会合坐标系就是该地月系统的惯性坐标系；在引力辅助坐标系中，rp和vp分别描述了在该坐标系下，探测器在近心点处的位置和速度。
已知探测器的初始状态近地距Rp、远地距Ra和轨道倾角i时，得到此时轨道的半长轴，在地球引力作用下，利用二体能量积分公式计算得出探测器经过影响球边界时的速度大小，从而根据探测器的角动量计算出速度的侧向分量大小，进而得到惯性坐标系下探测器进入影响球时的速度矢量
且进入影响球时的位置矢量为
根据圆锥曲线模型的假设，引力辅助发生时得到行星的速度矢量是一个常矢量，表示为
其中V2的值为1.023km/s，从而得到引力辅助坐标系下探测器进入影响球时的剩余速度矢量；在已知引力辅助之后探测器具备预想能量C3值的前提下，根据式(1)计算出探测器飞出影响球时的速度大小：
从而可以得到引力辅助坐标系下，探测器在近心点处时的位置和速度矢量：
由此，根据已求得的进入影响球的位置和速度矢量计算出引力辅助前的轨道根数。
3、行星引力辅助本质上是一个三体问题，建立三体问题动力学模型，，其运动方程在会合坐标系下可以表示为：
其中，r1和r2分别为探测器与天体M1和M2之间的距离，质量参数为μ，该模型的积分初始量为近心点处的位置和速度。
基于三体问题动力学模型，通过数值积分方法在时间上分别向前向后进行轨道递推可以获得引力辅助前的轨道根数及引力辅助之后的特征能量C3数据分布。
4、解析的圆锥曲线拼接法与三体问题之间的误差计算如下：
其中，a是半长轴，e是偏心率，i是轨道倾角，Ω是升交点赤经，w表示近地点幅角，最后C3表征三体问题下，特征能量的计算公式。
根据上述已经获得的两个方法之间的误差，利用径向基(RBF)神经网络对其进行学习，首先确定神经网络的结构，RBF神经网络的网络层数为3层，分别为输入层、隐含层、输出层，其中输入层到隐含层的权值由输入进来的信号样本来确定，隐含层到输出层的权值可通过求解线性方程组来获得。
针对地月逃逸为例，将探测器初始状态的轨道倾角以步长为3°进行选取，远地距每隔2000km进行采样，所以训练样本数为61×101个。故将径向基神经网络的结构定为“2-6060-6”，轨道倾角和远地距Ra作为两个输入信号，根据6060个样本数量确定了6060个感知单元和6个输出信号，输入层到隐含层的权值个数为2×6060，隐含层到输出层的权值个数为6060×6，输出层为两个动力学模型之间轨道根数的误差及特征能量值C3。利用径向基神经网络进行模型的训练，只需要选择不到20％的三体问题所需的积分量，就可以达到很好的近似效果。
在采用径向基神经网络预测的时候，主要用到了MATLAB神经网络工具箱中的newrbe函数，该函数的调用格式为：net＝newrbe(P,T,spread)，P和T分别表示神经网络的输入向量和输出向量，spread表征径向基函数的分布密度，该值的大小直接决定了网络的性能，通常需要对spread进行不同数值的尝试，以确定一个最优值。spread的默认值为1，通过对均匀选取的训练样本进行模型的训练，然后通过调用sim函数测试网络性能。如果实际输出结果与真实值的误差相差较大，则可以通过改变网络结构、数据样本和spread参数值进行调整，最终获得理想的网络模型。根据预期效果，在训练过程中需要对spread不断地修改与测试，最终确定spread＝0.05。
由此，根据已给的数据样本可以训练出一个满足预期精度的径向基神经网络模型。
5、根据上述已经获得的径向基神经网络的模型，在给定探测器的初始轨道倾角和远地距的情况下，带入网络模型预测出两个方法的误差值，然后将预测得出的误差值与圆锥曲线拼接方法得到的结果进行代数相加从而获得三体问题下轨道高精度近似解。
根据近似结果发现半长轴和轨道倾角的预测值与真实值的相对百分比在-10-3～10-3量级之间，近乎可以忽略其误差，偏心率的相对误差百分比在-0.02％～0.03％区间范围内，升交点赤经呈现的相对误差较大，其百分比为-0.25％～1.50％，近地点幅角的相对误差百分比在-0.025％～0.10％之间。引力辅助之后，探测器所具备的特征能量C3的大小，其相对误差百分比在-10-3～10-3量级之内。利用训练好的网络模型计算6060个测试样本，径向基神经网络只需要0.49s，然而采用动力学模型积分所需要的平均时间为106.57s。
本发明具体应用途径很多，以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以作出若干改进，这些改进也应视为本发明的保护范围。
Claims (6)
1.一种基于径向基神经网络的月球引力辅助逃逸轨道递推方法，其特征在于，包括步骤如下：
1)设置探测器初始状态参数；
2)利用解析的圆锥曲线拼接法计算探测器在近心点处的状态及引力辅助前的轨道根数；
3)在相同的近心点处，通过圆形限制性三体问题动力学模型进行轨道递推，从而获得探测器在引力辅助之前的轨道根数及引力辅助之后的特征能量C3值；
4)基于径向基神经网络学习解析的圆锥曲线拼接法与三体问题结果之间的误差关系；
5)基于上述步骤4)得到的误差值，将其与圆锥曲线拼接法计算结果进行代数相加运算，从而替代在轨道设计时三体问题所需的数值积分。
2.根据权利要求1所述的基于径向基神经网络的月球引力辅助逃逸轨道递推方法，其特征在于，所述步骤1)具体包括：设定探测器的初始特征能量C3值；设定探测器初始阶段状态约束条件参数：轨道倾角i、近星距Rp与远星距Ra；两个主天体质量M1、M2，质量参数比μ＝M2/(M1+M2)，次天体影响球半径rSOI及两个主天体之间的距离REM。
3.根据权利要求1所述的基于径向基神经网络的月球引力辅助逃逸轨道递推方法，其特征在于，所述步骤2)具体包括：依次建立以两个主天体质心为原点的会合坐标系，再分别建立以主体为中心的惯性坐标系和以次体为中心的引力辅助坐标系，会合坐标系的坐标轴原点为两个主天体的质心，x轴由M1指向M2，z轴与地月系统的角动量方向一致，y轴满足右手定则；当t＝0时的会合坐标系为该地月系统的惯性坐标系；在引力辅助坐标系中，rp和vp分别描述在该坐标系下，探测器在近心点处的位置和速度；已知探测器的初始状态近星距Rp、远星距Ra和轨道倾角i时，得到此时轨道的半长轴，在主体引力作用下，利用二体能量积分公式计算得出探测器经过影响球边界时的速度大小，从而根据探测器的角动量计算出速度的侧向分量大小，进而得到惯性坐标系下探测器进入影响球时的速度矢量
且进入影响球时的位置矢量为
根据圆锥曲线模型的假设，引力辅助发生时得到行星的速度矢量是一个常矢量，表示为
从而得到引力辅助坐标系下探测器进入影响球时的剩余速度矢量；在已知引力辅助之后探测器具备预想能量C3值的前提下，根据式(1)计算出探测器飞出影响球时的速度大小：
从而得到引力辅助坐标系下，探测器在近心点处时的位置和速度矢量：
根据已求得的进入影响球的位置和速度矢量计算出引力辅助前的轨道根数。
4.根据权利要求1所述的基于径向基神经网络的月球引力辅助逃逸轨道递推方法，其特征在于，所述步骤3)具体包括：
建立三体问题动力学模型，其运动方程在会合坐标系下表示为：
式中，r1和r2分别为探测器与天体M1和M2之间的距离，质量参数为μ，该数学模型的积分初始量为近心点处的位置和速度；
基于三体问题动力学模型，通过数值积分方法在时间上分别向前向后进行轨道递推获得引力辅助前的轨道根数及引力辅助之后的特征能量C3数据分布。
5.根据权利要求1所述的基于径向基神经网络的月球引力辅助逃逸轨道递推方法，其特征在于，所述步骤4)具体包括：
解析的圆锥曲线拼接法与三体问题之间的误差计算如下：
式中，a表示半长轴，e是偏心率，i是轨道倾角，Ω是升交点赤经，w是近地点幅角，C3表征三体问题下，特征能量的计算公式；
根据上述已经获得的误差，利用径向基神经网络对其进行学习；确定神经网络的结构，径向基神经网络的网络层数为3层，分别为输入层、隐含层、输出层，输入层到隐含层的权值由输入进来的信号样本来确定，隐含层到输出层的权值通过求解线性方程组来获得；
在采用径向基神经网络预测时，使用MATLAB神经网络工具箱中的newrbe函数，该函数的调用格式为：net＝newrbe(P,T,spread)，P和T分别表示神经网络的输入向量和输出向量，spread表征径向基函数的分布密度，对spread进行不同数值的尝试，以确定一个最优值；spread的默认值为1，通过对均匀选取的训练样本进行模型的训练，然后通过调用sim函数测试网络性能；如果实际输出结果与真实值的误差相差大，则通过改变网络结构、数据样本和spread参数值进行调整，最终获得理想的网络模型；根据预期效果，在训练过程中需要对spread不断地修改与测试，最终确定spread值。
6.根据权利要求1所述的基于径向基神经网络的月球引力辅助逃逸轨道递推方法，其特征在于，所述步骤5)具体包括：根据已经获得的径向基神经网络的模型，在给定探测器的初始轨道倾角和远星距的情况下，带入径向基神经网络模型预测出两个方法的误差值，然后将预测得出的误差值与圆锥曲线拼接方法得到的结果进行代数相加从而获得三体问题下轨道高精度近似解。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011326683.1A CN112613223B (zh) | 2020-11-24 | 2020-11-24 | 基于径向基神经网络的月球引力辅助逃逸轨道递推方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011326683.1A CN112613223B (zh) | 2020-11-24 | 2020-11-24 | 基于径向基神经网络的月球引力辅助逃逸轨道递推方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112613223A true CN112613223A (zh) | 2021-04-06 |
| CN112613223B CN112613223B (zh) | 2024-04-16 |
Family
ID=75225544
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011326683.1A Active CN112613223B (zh) | 2020-11-24 | 2020-11-24 | 基于径向基神经网络的月球引力辅助逃逸轨道递推方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112613223B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113343561A (zh) * | 2021-05-25 | 2021-09-03 | 中国科学院国家空间科学中心 | 一种求解最优月球飞掠转移轨道的方法及系统 |
| CN113591403A (zh) * | 2021-09-28 | 2021-11-02 | 南京航空航天大学 | 一种航天器借力飞行轨道的设计方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102923324A (zh) * | 2012-11-29 | 2013-02-13 | 北京理工大学 | 基于不变流形与引力辅助的低能量行星逃逸轨道设计方法 |
| CN107451656A (zh) * | 2017-08-08 | 2017-12-08 | 上海航天控制技术研究所 | 一种火星探测器器上自主轨道计算方法 |
| CN107798380A (zh) * | 2017-09-12 | 2018-03-13 | 上海航天控制技术研究所 | 探测器器上星历计算方法和计算系统 |
| CN110096726A (zh) * | 2019-02-21 | 2019-08-06 | 上海卫星工程研究所 | 基于月球借力的geo卫星应急转移轨道快速优化设计方法 |
-
2020
- 2020-11-24 CN CN202011326683.1A patent/CN112613223B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102923324A (zh) * | 2012-11-29 | 2013-02-13 | 北京理工大学 | 基于不变流形与引力辅助的低能量行星逃逸轨道设计方法 |
| CN107451656A (zh) * | 2017-08-08 | 2017-12-08 | 上海航天控制技术研究所 | 一种火星探测器器上自主轨道计算方法 |
| CN107798380A (zh) * | 2017-09-12 | 2018-03-13 | 上海航天控制技术研究所 | 探测器器上星历计算方法和计算系统 |
| CN110096726A (zh) * | 2019-02-21 | 2019-08-06 | 上海卫星工程研究所 | 基于月球借力的geo卫星应急转移轨道快速优化设计方法 |
Non-Patent Citations (3)
| Title |
|---|
| 杨彬 等: "高精度火星大气制动轨迹智能高效优化方法", 中国科学:技术科学, no. 09, 13 September 2020 (2020-09-13) * |
| 杨洪伟 等: "精确动力学模型中的行星引力辅助轨道设计", 中国空间科学技术, no. 02, 25 April 2013 (2013-04-25) * |
| 贾建华 等: "三维引力辅助解析分析方法研究", 力学学报, no. 03, 18 May 2013 (2013-05-18) * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113343561A (zh) * | 2021-05-25 | 2021-09-03 | 中国科学院国家空间科学中心 | 一种求解最优月球飞掠转移轨道的方法及系统 |
| CN113343561B (zh) * | 2021-05-25 | 2021-12-28 | 中国科学院国家空间科学中心 | 一种求解最优月球飞掠转移轨道的方法及系统 |
| CN113591403A (zh) * | 2021-09-28 | 2021-11-02 | 南京航空航天大学 | 一种航天器借力飞行轨道的设计方法 |
| CN113591403B (zh) * | 2021-09-28 | 2022-01-04 | 南京航空航天大学 | 一种航天器借力飞行轨道的设计方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112613223B (zh) | 2024-04-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112613223B (zh) | 基于径向基神经网络的月球引力辅助逃逸轨道递推方法 | |
| CN112257343B (zh) | 一种高精度地面轨迹重复轨道优化方法及系统 | |
| CN110378012B (zh) | 一种考虑高阶重力场的严格回归轨道设计方法、系统及介质 | |
| CN107944084B (zh) | 一种深空探测器行星借力飞行轨道评估方法 | |
| CN108196449A (zh) | 最优着陆轨迹设计的协变量初值确定方法 | |
| CN104833466B (zh) | 一种航天器地面测试与在轨微振动力学环境映射方法 | |
| CN105631095A (zh) | 一种等间隔发射的多约束地月转移轨道簇搜索方法 | |
| CN102923324A (zh) | 基于不变流形与引力辅助的低能量行星逃逸轨道设计方法 | |
| de Almeida Prado | A comparison of the “patched-conics approach” and the restricted problem for swing-bys | |
| Khoury | Orbital rendezvous and spacecraft loitering in the earth-moon system | |
| CN115314101A (zh) | 一种基于并行计算的低轨通信卫星星座快速建模方法 | |
| McCarthy et al. | Cislunar transfer design exploiting periodic and quasi-periodic orbital structures in the four-body problem | |
| Sood et al. | L4, L5 solar sail transfers and trajectory design: Solar observations and potential Earth Trojan exploration | |
| Cho et al. | Optimal perilune altitude of lunar landing trajectory | |
| CN108100306B (zh) | 一种快速筛选日地平衡点小行星探测目标的方法 | |
| Rozek et al. | Multi-objective optimisation of NRHO-LLO orbit transfer via surrogate-assisted evolutionary algorithms | |
| CN107844462A (zh) | 一种行星际连续推力转移轨道评估方法 | |
| Yan et al. | ANN-based method for fast optimization of Jovian-moon gravity-assisted trajectories in CR3BP | |
| CN108082538A (zh) | 一种考虑始末约束的多体系统低能量捕获轨道方法 | |
| CN111428912A (zh) | 一种基于支持向量机的火星探测器轨道预测方法及系统 | |
| Song et al. | Fast estimation of gravitational field of irregular asteroids based on deep neural network and its application | |
| CN112580255B (zh) | 一种基于飞越映射参数学习的木卫引力辅助轨道递推方法 | |
| Song et al. | A lunar cargo mission design strategy using variable low thrust | |
| CN112009727B (zh) | 平动点轨道最优小推力转移分段设计方法 | |
| Biswas et al. | Evolved general relativity predicts “sought-after departure” of Space Time Asymmetry Research mission |
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 |