CN111259337B - 一种基于统计的重残骸实时落点预报方法 - Google Patents
一种基于统计的重残骸实时落点预报方法 Download PDFInfo
- Publication number
- CN111259337B CN111259337B CN202010042658.4A CN202010042658A CN111259337B CN 111259337 B CN111259337 B CN 111259337B CN 202010042658 A CN202010042658 A CN 202010042658A CN 111259337 B CN111259337 B CN 111259337B
- Authority
- CN
- China
- Prior art keywords
- component
- debris
- acceleration
- velocity
- calculation
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于统计的重残骸实时落点预报方法,重残骸落区规划与回收救援是运载火箭执行飞行方案的重要组成部分,目前运载火箭重残骸的实时计算受多种不确定测量因素影响,其预报落点与实际落点还存在较大的误差,导致误差产生的原因主要包括运载器特征点弹道信息获取误差、重残骸飞行信息计算误差、大气阻滞计算误差以及运动模型计算误差。本发明计算方法引进了大数据综合分析理念,立足航天重残骸测控现状,充分利用历次发射宝贵实测数据资源,对影响落点计算的风阻、速度修正量进行反演计算,形成了基于统计原理的实时落点计算模型,通过历次数据验证表明,新模型落点预测精度较高。
Description
技术领域
本发明涉及残骸分离后的落地预测技术,具体是一种基于统计原理的实时积分外推落点预报方法。
背景技术
一般卫星发射不同于载人航天回收任务,各子级残骸分离后是没有外测测量设备继续跟踪的,因此,在进行落点计算时,残骸初始运动信息是无法通过直接测量得出的,而是通过分离时刻运载器弹道信息估算得出的。其原因是运载器分离时刻箭体运动矢量过于复杂(除了在射向方向加速运动外,箭体自身相对运动方向还有方位、俯仰和滚动等矢量),同时分离时的相对作用力大小、方向和作用时间都是一个模糊的控制量。既然精确计算比较困难,任务实施中就只能利用经验值的模式进行修正。
在落点计算中,航天器被动段的运动可视为质点运动,我们仅从物体运动学原理方面理解,航天器被动段飞行主要受地球引力、哥氏力和空气阻力的影响。被动段飞行从受力不同上可分为自由飞行段和再入段,自由飞行段航天器只在地球引力和自转引起的哥氏力的作用下运动,再入段与自由飞行段的差别仅仅增加了空气阻力的影响。运动方程如下所示:
上式运动方程中,空气阻力项:XD
其中SM为残骸迎风面的面积,v为残骸综合速度,g0为地面重力加速度,ρ/ρ0为跟高度有关的空气密度,可通过查落区大气密度表得到,CD为空气阻力系数,其值与残骸速度有关,可通过查马赫数对应的气阻系数表得到。从公式可以看出,在任务实际执行中,残骸迎风面的面积的精确值不容易获取,随高度变化的空气密度获取也不够精准,这些都是导致计算大气阻力不精确的因素。
落点计算是由运动方程构成的一个微分方程,在进行数值计算中需设计一定的步长进行迭代计算,选择适当的步长是非常重要的,特别是实时任务中,步长过小会导致计算量过大而影响实时性需求,过大则造成计算误差过大的情况发生。
发明内容
为解决上述问题,本发明的目的在于提供一种基于历史数据统计反演的实时落点预报方法。具体涉及残骸分离信息计算、大气阻滞补偿计算、无动力陨落加速度计算、数值迭代计算,以及基于历史真实数据的大气阻滞综合系数、修正参数ΔV和ΔK的数值反演。
本发明具体包括以下,(1)在分离点弹道计算时采用单向高阶最小二乘非线性拟合来消除分离时刻弹道为拐点的影响,同时为进一步提高分离时刻弹道的精度,对残骸初始运动速度信息进行了修正;(2)在常规大气阻力模型的基础之上,通过分析大气密度随高度h的分布规律,用最小二乘拟合了大气密度随高度变化的指数,结合火箭残骸的截面、重量、马赫数的影响,建立了含大气阻滞综合系数、马赫数影响系数以及大气密度指数的大气阻滞模型,同时为进一步提高大气阻滞综合系数精度,对大气阻滞综合系数进行补偿修正;(3)综合考虑地球引力、哥氏力和空气阻力等因素的影响,计算无动力陨落加速度;(4)进行数值迭代计算落点,实现落点的预报,采用5ms步长的“龙格-库塔”模型进行迭代计算,迭代计算终止标识一般设置为落点相对地球椭球体表面的高程H小于500米;(5)结合分离作用力影响和历史经验,设计了修正参数的取值范围,ΔV取值范围为[2,15]和ΔK取值范围为[0.5,4.5];(6)基于历史数据统计反演大气阻滞系数、大气阻滞修正系数ΔK、修正速度ΔV。
鉴于此,本发明采用的技术方案是,一种基于统计的重残骸实时落点预报方法,包括以下步骤:
基于历史数据进行反演,获取最佳修正速度ΔV和最佳修正系数ΔK。
从遥测时间中获取最佳残骸分离时间信息。
从运载器众多飞行弹道中挑选最优飞行弹道。
计算残骸分离时刻弹道信息,包括弹道速度信息和弹道位置信息。
根据最佳修正速度ΔV修正残骸分离模型中残骸分离时刻弹道速度信息。
根据最佳修正系数ΔK修正大气阻滞补偿模型。
计算大气阻力。
根据无动力陨落加速度模型计算无动力加速度。
按迭代算法和时间迭代步长进行数值迭代。
判断迭代是否满足终止条件,若满足则输出残骸预示落点,若不满足则执行计算大气阻力步骤。
采用上述技术的原因是落点计算时无法准确获取残骸初始运动信息和大气阻力。
本发明的有益技术效果是:
1、充分考虑了运载器分离时刻运动过程复杂性和弹道特性,采用单向高阶最小二乘非线性拟合来消除分离时刻弹道为拐点的影响。
2、对残骸初始运动速度信息进行修正,进一步提高残骸初始运动速度信息精度。
3、综合考虑多方面因素,建立大气阻滞模型,基于历史数据统计反演大气阻滞系数。
4、对大气阻滞系数进行修正,进一步提高大气阻滞系数精度。
5、结合分离作用力影响和历史经验,设计了修正参数(修正系数和修正速度)的取值范围。
6、基于历史数据统计反演大气阻滞修正系数和修正速度。
7、结合计算精度和实时计算需求,设计了合适的迭代方法。
8、本发明编程易实现,提高了实时计算落点的效率和精度。
附图说明
图1为本发明的处理流程图。
具体实施方式
为使本发明的目的、技术方案及优点得到更好的理解,现以具体实施例并结合附图,具体说明如下:
基于统计原理的实时落点预报方法是以常规运动模型为基础建立,根据综合误差模糊计算理论,在基础计算模型中增加两个模糊修正量,一个用来补偿计算与实际运动的速度误差;一个用来调节模型计算与实际风阻的关系,同时根据重残骸运动特性,设定值域减少迭代计算次数,包括:
(1)建立残骸分离信息计算模型
进行残骸落点计算第一步是需要知道残骸分离时刻的弹道信息,然后才能根据加速度情况进行数值计算。这里之所以采用高阶拟合计算,主要原因是因为在残骸分离前后,其主体火箭加速度变化较大,如果利用线性拟合计算将导致计算结果误差较大。模型利用分离前、后弹道信息分别拟合,最后利用两组结果求均值得出。
设yj,j=1,2,……n为(n为多项式回归拟合段落点数)弹道参数(包括位置信息和速度信息)时间序列t的第j个采样点测量值;bj为yj序列的多项式回归系数;p为拟合多项式的阶数(取为2,3,4……;),解算模型为:
B=(DTD)-1DTY (1)
其中:
bj,j=1,2,……p为待求多项式回归系数,p为待求多项式阶数;ti,i=1,2,……n为弹道序列对应的时间,yi,i=1,2,……n为为弹道信息中时间序列t的第j个采样点测量值。
y=b0+b1x+…+bpxp,p=2,3… (3)
公式(3)即为最小二乘非线性拟合方程,将分离点时刻带入上式即可求出拟合弹道信息,本方法仿真计算中拟合多项式阶数取4。x表示数据拟合时间,y表示x时刻的拟合值。
假设经模型公式(3)解算的残骸分离的速度、位置信息分别为:V0=(vx0 vy0 vz0),P0=(x0 y0 z0)。vx0 vy0 vz0分别表示速度在X方向的分量、速度在Y方向的分量、速度在Z方向的分量,x0 y0 z0分别表示位置在X方向的分量、位置在Y方向的分量、位置在Z方向的分量。加入补偿修正量(ΔV)后,其分离信息(根据分离特性,只做速度补偿,位置补偿不做修正)为:
式中
xd=x0;yd=y0;zd=z0。
(2)建立大气阻滞补偿计算模型
大气阻滞加速度矢量模型推导过程较为复杂,严格意义上需要考虑残骸陨落过程空域的大气密度、风向等指标,实际任务中,往往采用地区通用大气密度数据进行计算,本发明中模型是在常规大气阻力模型的基础之上,通过分析大气密度随高度h的分布规律,用最小二乘拟合了大气密度随高度变化的指数,结合火箭残骸的截面、重量、马赫数的影响,建立了含大气阻滞综合系数、马赫数影响系数以及大气密度指数的大气阻滞模型:
gbx表示阻力加速度在X方向的分量,gby表示阻力加速度在Y方向的分量,gbz表示阻力加速度在Z方向的分量,vx表示速度在X方向的分量,vy表示速度在Y方向的分量,vz表示速度在Z方向的分量,V表示综合速度。
指数id是一个跟高度(公里)有关的一个高阶拟合值,拟合方程如下:
a1……a8,b1……b6表示指数系数,如表1所示。hk表示高度。
参数kv是与残骸飞行的马赫数有关的一个拟合值。
cv表示马赫数;k为大气阻滞综合系数。
表1指数系数表
补偿系数ΔK带入公式(5),从而形成带补偿系数的大气阻滞计算模型。
(3)建立无动力陨落加速度模型
残骸陨落过程中,自身已无动力***工作,下坠过程中除受大气阻力影响外,就是受引力影响,这里的引力一般指地球引力与因地球自转而产生的哥氏力,此力是残骸下坠过程中影响最大,计算模型如下:
其中gxd表示重力加速度在X方向的分量,gyd表示重力加速度在Y方向的分量,gzd表示重力加速度在Z方向的分量。
上式中里斯常数J=0.00108263,u=3.986004415e14为地球引力常数,we=0.000072921151467为地球自转角速率,地球半径Re=6378245.0。
残骸无动力陨落过程中,其运动加速度由两部分构成,一个是引力加速度,一个是大气阻滞加速度,因此残骸无动力陨落加速度模型可由公式(6)和公式(7)得出:
(4)进行落点数值迭代计算方法
残骸陨落轨迹和落点的计算,根据物体运动学原理,已知物体运动初始速度和位置,在加速度可计算的条件下,即可计算出下一时刻的速度和位置,在进行积分迭代计算中,时间步长不宜取的过大,一般不超过10毫秒,实际计算中所取步长越小,其迭代计算精度越高,但同时计算所用资源也随着步长的减少而成倍增加,实际应用中需结合计算需要和计算机运算速度来综合确定(仿真计算t取5ms)。迭代计算公式见公式(9)。在结合分离作用力影响和历史经验,在进行计算中ΔV取值范围为[2,15],步长取0.24,ΔK取值范围为[0.5,4.5],步长取1/50.0。
参见图1,基于历史数据统计反演修正系数流程:
1、历史数据收集整理,包括残骸实际落点、残骸分离时间、运载器飞行弹道等。
2、结合分离作用力影响和历史经验,设计ΔV取值范围[2,15],迭代步长0.24,ΔK取值范围[0.5,4.5],迭代步长1/50.0。
3、设计迭代算法、时间迭代步长、迭代终止条件等。
4、计算残骸分离时刻弹道信息。
5、对残骸分离时刻弹道速度信息按修正步长进行修正,并判断修正量是否超出取值范围,若超出取值范围执行步骤11,若未超出取值范围执行步骤6。
6、对大气阻滞系数按修正步长进行修正,并判断修正量是否超出取值范围,若超出取值范围执行步骤5,若未超出取值范围执行步骤7。
7、计算大气阻力。
8、计算无动力加速度。
9、按迭代算法和时间迭代步长进行数值迭代。
10、判断迭代是否满足终止条件,若满足则执行步骤6,若不满足执行步骤7。
11、将计算结果与历史实际落点进行对比分析。
12、选取最优修正系数。
基于统计原理的实时落点预报:
1、设计迭代算法、时间迭代步长、迭代终止条件等。
2、装填最优修正系数。
3、从遥测时间中获取最佳残骸分离时间信息。
4、从运载器众多飞行弹道中挑选最优飞行弹道。
5、计算残骸分离时刻弹道信息。
6、修正残骸分离时刻弹道速度信息。
7、修正大气阻滞系数。
8、计算大气阻力。
9、计算无动力加速度。
10、按迭代算法和时间迭代步长进行数值迭代。
11、判断迭代是否满足终止条件,若满足则执行步骤12,若不满足执行步骤7。
12、输出残骸预示落点。
Claims (7)
1.一种基于统计的重残骸实时落点预报方法,其特征在于,包括以下步骤:
基于历史数据进行反演,获取最佳修正速度ΔV和最佳修正系数ΔK;
从遥测时间中获取最佳残骸分离时间信息;
从运载器众多飞行弹道中挑选最优飞行弹道;
计算残骸分离时刻弹道信息,包括弹道速度信息和弹道位置信息;
根据最佳修正速度ΔV修正残骸分离模型中残骸分离时刻弹道速度信息;
根据最佳修正系数ΔK修正大气阻滞补偿模型,修正后的大气阻滞补偿模型为
式中,gbx表示阻力加速度在X方向的分量,gby表示阻力加速度在Y方向的分量,gbz表示阻力加速度在Z方向的分量,vx表示速度在X方向的分量,vy表示速度在Y方向的分量,vz表示速度在Z方向的分量,V表示综合速度;
指数id的拟合方程如下:
a1……a8,b1……b6表示指数系数,hk表示高度;
参数kv是与残骸飞行的马赫数有关的一个拟合值;
cv为马赫数;k为大气阻滞综合系数;
计算大气阻力;
根据无动力陨落加速度模型计算无动力加速度;
按迭代算法和时间迭代步长进行数值迭代;
判断迭代是否满足终止条件,若满足则输出残骸预示落点,若不满足则执行计算大气阻力步骤;
所述基于历史数据进行反演,包括以下步骤:
1)历史数据收集整理,包括残骸实际落点、残骸分离时间、运载器飞行弹道;
2)结合分离作用力影响和历史经验,设计修正速度和修正系数;
3)设计迭代算法、时间迭代步长、迭代终止条件;
4)计算残骸分离时刻弹道信息;
5)对残骸分离时刻弹道速度信息按修正速度进行修正,并判断修正量是否超出取值范围,若超出取值范围执行步骤11),若未超出取值范围执行步骤6);
6)对大气阻滞系数按修正系数进行修正,并判断修正量是否超出取值范围,若超出取值范围执行步骤5),若未超出取值范围执行步骤7);
7)计算大气阻力;
8)计算无动力加速度;
9)按迭代算法和时间迭代步长进行数值迭代;
10)判断迭代是否满足终止条件,若满足则执行步骤6),若不满足执行步骤7);
11)将计算结果与历史实际落点进行对比分析;
12)选取最佳修正速度ΔV和最佳修正系数ΔK。
2.根据权利要求1所述一种基于统计的重残骸实时落点预报方法,其特征在于:所述残骸分离时刻弹道信息利用分离前、后弹道信息分别拟合,最后利用两组结果求均值得出。
3.根据权利要求2所述一种基于统计的重残骸实时落点预报方法,其特征在于:所述弹道信息拟合采用单向高阶最小二乘非线性拟合。
7.根据权利要求1所述一种基于统计的重残骸实时落点预报方法,其特征在于:结合分离作用力影响和历史经验,设计的修正速度ΔV和修正系数ΔK,ΔV取值范围为[2,15]和ΔK取值范围为[0.5,4.5]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010042658.4A CN111259337B (zh) | 2020-01-15 | 2020-01-15 | 一种基于统计的重残骸实时落点预报方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010042658.4A CN111259337B (zh) | 2020-01-15 | 2020-01-15 | 一种基于统计的重残骸实时落点预报方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111259337A CN111259337A (zh) | 2020-06-09 |
CN111259337B true CN111259337B (zh) | 2023-03-31 |
Family
ID=70950566
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010042658.4A Active CN111259337B (zh) | 2020-01-15 | 2020-01-15 | 一种基于统计的重残骸实时落点预报方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111259337B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112113471B (zh) * | 2020-08-18 | 2022-08-16 | 中国人民解放军92941部队 | 基于最优模糊***的惯导测量弹道修正方法 |
CN113449460B (zh) * | 2021-05-20 | 2023-06-20 | 北京理工大学 | 破片阻力系数的测量方法、装置及计算机可读存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105203110A (zh) * | 2015-08-28 | 2015-12-30 | 中国科学院空间应用工程与技术中心 | 一种基于大气阻力模型补偿的低轨卫星轨道预报方法 |
CN107679655A (zh) * | 2017-09-15 | 2018-02-09 | 中国人民解放军63816部队 | 一种航天发射火箭落点预测*** |
CN109323698A (zh) * | 2018-12-03 | 2019-02-12 | 西安四方星途测控技术有限公司 | 一种空间目标陨落多模型跟踪引导技术 |
CN109798902A (zh) * | 2019-03-11 | 2019-05-24 | 北京星际荣耀空间科技有限公司 | 一种适用于运载火箭入轨修正的迭代制导方法 |
US10302325B1 (en) * | 2016-04-25 | 2019-05-28 | American Air Filter Company, Inc. | Air filter operation and use modification according to identify a filter resistance of an air handling unit |
CN110489879A (zh) * | 2019-08-22 | 2019-11-22 | 中国人民解放军32035部队 | 一种适用于空间环境扰动情况下的空间目标陨落预报方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10372575B1 (en) * | 2018-01-31 | 2019-08-06 | Dell Products L.P. | Systems and methods for detecting and removing accumulated debris from a cooling air path within an information handling system chassis enclosure |
-
2020
- 2020-01-15 CN CN202010042658.4A patent/CN111259337B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105203110A (zh) * | 2015-08-28 | 2015-12-30 | 中国科学院空间应用工程与技术中心 | 一种基于大气阻力模型补偿的低轨卫星轨道预报方法 |
US10302325B1 (en) * | 2016-04-25 | 2019-05-28 | American Air Filter Company, Inc. | Air filter operation and use modification according to identify a filter resistance of an air handling unit |
CN107679655A (zh) * | 2017-09-15 | 2018-02-09 | 中国人民解放军63816部队 | 一种航天发射火箭落点预测*** |
CN109323698A (zh) * | 2018-12-03 | 2019-02-12 | 西安四方星途测控技术有限公司 | 一种空间目标陨落多模型跟踪引导技术 |
CN109798902A (zh) * | 2019-03-11 | 2019-05-24 | 北京星际荣耀空间科技有限公司 | 一种适用于运载火箭入轨修正的迭代制导方法 |
CN110489879A (zh) * | 2019-08-22 | 2019-11-22 | 中国人民解放军32035部队 | 一种适用于空间环境扰动情况下的空间目标陨落预报方法 |
Non-Patent Citations (3)
Title |
---|
Research on neural network based debris flow prediction;Zhang Wenjiang;《2014 International Conference on Information Science, Electronics and Electrical Engineering》;20141106;全文 * |
于古胜 ; 李连登 ; 翟丽丽.航天器实时落点计算误差修正方法.《舰船电子工程》.2010, * |
落点预估法实现舰空导弹飞行试验安全判定;张远;《现代防御技术》;20170415;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111259337A (zh) | 2020-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107679655B (zh) | 一种航天发射火箭落点预测*** | |
CN108646555B (zh) | 基于航天器回收***空投试验数据的参数辨识方法 | |
CN111259337B (zh) | 一种基于统计的重残骸实时落点预报方法 | |
CN113479347B (zh) | 一种火箭垂直回收着陆段轨迹控制方法 | |
CN107544067A (zh) | 一种基于高斯混合近似的高超声速再入飞行器跟踪方法 | |
CN108549785B (zh) | 一种基于三维飞行剖面的高超声速飞行器精准弹道快速预测方法 | |
CN109323698B (zh) | 一种空间目标陨落多模型跟踪引导方法 | |
CN110015446B (zh) | 一种半解析的火星进入制导方法 | |
WO2021102669A1 (zh) | 超低轨道卫星轨道自主维持方法 | |
CN109446582B (zh) | 一种考虑地球自转的高精度降阶平稳滑翔动力学建模方法 | |
CN107368085B (zh) | 基于模型预测的风场中平流层飞艇高度控制方法 | |
CN103708045A (zh) | 一种探月飞船跳跃式再入的在线参数辨识方法 | |
Wenz et al. | Estimation of wind velocities and aerodynamic coefficients for UAVs using standard autopilot sensors and a moving horizon estimator | |
CN107101636A (zh) | 一种使用卡尔曼滤波器辨识多旋翼动力学模型参数的方法 | |
CN109901387A (zh) | 一种航空器自动近地防撞***自适应飞行轨迹预测方法 | |
CN112683261A (zh) | 一种基于速度预测的无人机鲁棒性导航方法 | |
CN113468828B (zh) | 一种飞机空中飞行颠簸强度指数计算方法 | |
CN114462293B (zh) | 一种高超声速目标中长期轨迹预测方法 | |
CN107796401B (zh) | 跳跃式再入飞行器线性伪谱参数修正横向制导方法 | |
CN109765141A (zh) | 一种基于swarm-c卫星提取大气密度的方法 | |
Halaas et al. | Formation flight for aerodynamic benefit simulation development and validation | |
CN116341377A (zh) | 一种基于lstm神经网络的下投式探测组件轨迹预测方法 | |
CN113359861B (zh) | 一种无人飞艇编队飞行控制方法及*** | |
CN111651860B (zh) | 一种可重复使用运载器再入段的预测校正鲁棒制导方法 | |
CN111382514A (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 |