CN103112604A - 一种卫星轨道控制方法 - Google Patents
一种卫星轨道控制方法 Download PDFInfo
- Publication number
- CN103112604A CN103112604A CN2013100363994A CN201310036399A CN103112604A CN 103112604 A CN103112604 A CN 103112604A CN 2013100363994 A CN2013100363994 A CN 2013100363994A CN 201310036399 A CN201310036399 A CN 201310036399A CN 103112604 A CN103112604 A CN 103112604A
- Authority
- CN
- China
- Prior art keywords
- momentum
- attitude
- track
- satellite
- attitude control
- 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
Images
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种卫星轨道控制方法,能够在偏航轴和滚动轴的其中一个轴的姿态控制发动机工作不正常情况下实现卫星的轨道控制,首先对姿态控制发动机工作正常的两个轴进行动量轮卸载条件判断,当满足动量轮卸载条件时,对姿态控制发动机工作正常的两个轴进行动量轮卸载,对姿态控制发动机工作不正常的轴采用动量轮进行姿态控制;当不满足动量轮卸载条件时,采用动量轮进行三轴姿态控制直至卫星三轴姿态角和角速度均小于设定的角度阈值和角速度阈值,然后计算轨道控制时间,轨道控制发动机按照所计算的轨道控制时间进行轨道控制。本发明的方法可实现欠驱动情况下的快速轨道机动。
Description
技术领域
本发明涉及一种卫星轨道控制方法,可直接应用于具有欠驱动推力器和动量轮的卫星轨道控制领域。
背景技术
以往卫星的轨道控制大都是基于三轴姿态控制发动机喷气控制实现姿态稳定,或者在轨道控制发动机干扰力矩较小且轨道控制时间较短时,利用动量轮可吸收轨道控制干扰力矩时,利用轮控维持姿态稳定。但是,在对地定向的卫星姿态控制发动机故障(即出现欠驱动的情况)、轨道控制发动机干扰力矩较大且要求轨道调整量较大时,无法仅利用姿态控制发动机或者动量轮实现姿态稳定。
目前没有使用欠驱动推力器和动量轮实现大干扰力矩情况下的快速轨道控制方法。
发明内容
本发明所要解决的技术问题是:针对轨道控制发动机工作时会产生较大干扰力矩,姿态控制发动机为欠驱动的情况且要求轨道调整量较大时,提出了一种基于姿态控制发动机和动量轮结合的卫星轨道控制方法,可实现欠驱动情况下的快速轨道机动。
本发明包括如下技术方案:一种卫星轨道控制方法,能够在偏航轴和滚动轴的其中一个轴的姿态控制发动机工作不正常情况下实现卫星的轨道控制,包括如下步骤:
(1)对姿态控制发动机工作正常的两个轴进行动量轮卸载条件判断,当满足动量轮卸载条件时,转入步骤(5);当不满足动量轮卸载条件时,转入步骤(2);
(2)采用动量轮进行三轴姿态控制;
(3)判断卫星三轴姿态角和角速度是否均小于设定的角度阈值和角速度阈值;如果卫星三轴姿态角和角速度均小于设定的角度阈值和角速度阈值则转入步骤(4);否则转入步骤(2);
(4)计算轨道控制时间,轨道控制发动机按照所计算的轨道控制时间进行轨道控制,在轨道控制结束后返回步骤(1);
(5)对姿态控制发动机工作正常的两个轴进行动量轮卸载,对姿态控制发动机工作不正常的轴采用动量轮进行姿态控制;然后返回步骤(1)。
所述步骤(1)中,若当前角动量与标称角动量的差值大于设定的角动量偏差阈值,则满足动量轮卸载条件;若当前角动量与标称角动量的差值不大于设定的角动量偏差阈值,则不满足动量轮卸载条件。
所述轨道控制时间的根据如下公式计算:
t△=(Hz△-Hz-Iz*wz)/Tgz
其中:t△为轨道控制时间;Hz△为姿态控制发动机工作不正常的轴的动量轮的角动量阈值;Hz为姿态控制发动机工作不正常的轴的实际角动量;Iz为姿态控制发动机工作不正常的轴的转动惯量;wz为姿态控制发动机工作不正常的轴的角速度;Tgz为轨道控制发动机在姿态控制发动机工作不正常的轴上产生的干扰力矩。
本发明与现有技术相比具有如下优点:
本发明针对在偏航轴和滚动轴姿态控制发动机的其中一个轴工作不正常情况下、轨道控制发动机工作时产生较大干扰力矩的卫星,利用轨道罗盘原理,将偏航轴(或滚动轴)的干扰角动量转移到滚动轴(或偏航轴)上,并采用姿态控制发动机的推力器对工作正常两个轴的动量轮卸载,转移并消除干扰角动量。
本发明利用姿态控制发动机和动量轮结合进行轨道控制,在偏航轴有较大干扰力矩、姿态控制发动机为欠驱动且要求轨道调整量较大的情况下实现了在较短时间内的连续轨道控制。
附图说明
图1为本发明的卫星轨道控制方法的流程图。
图2为卫星控制系统的组成示意图。
图3为三轴的姿态角曲线示意图;图3a为滚动轴姿态角曲线示意图;图3b为俯仰轴姿态角曲线示意图,图3c为偏航轴姿态角曲线示意图;其中横坐标为时间(s),纵坐标为姿态角(度)。
图4为三轴的姿态角速度曲线示意图;图4a为滚动轴姿态角速度曲线示意图,图4b为俯仰轴姿态角速度曲线示意图,图4c为偏航轴姿态角速度曲线示意图;其中横坐标为时间(s),纵坐标为姿态角速度(度/秒)。
具体实施方式
下面就结合附图对本发明做进一步介绍。
所述欠驱动的情况为偏航或滚动姿态控制发动机工作不正常的情况,下面以偏航轴姿态控制发动机工作不正常、且轨道控制过程中产生偏航干扰力矩的情况为例进行说明,对于滚动轴的姿态控制发动机工作不正常的情况可类同。
如图1所示,本发明的卫星轨道控制方法步骤如下:
(1)对滚动、俯仰轴进行动量轮卸载条件判断;当满足动量轮卸载条件时,转入步骤(5);当不满足动量轮卸载条件时,转入步骤(2);
若当前角动量与标称角动量的差值大于设定的角动量偏差阈值,则认为满足动量轮卸载条件;若当前角动量与标称角动量的差值不大于设定的角动量偏差阈值,则认为不满足动量轮卸载条件。所述阈值根据动量轮可工作的范围确定。
(2)采用动量轮进行俯仰、偏航和滚动三轴姿态控制;
一般根据卫星姿态角偏差和角速度偏差,采用一定的控制律(如PID控制律)计算出三轴姿态控制所需的控制力矩,将控制力矩分配到三轴动量轮上,由动量轮产生控制力矩实现对三轴姿态的控制。
(3)判断卫星三轴姿态角和角速度是否均小于设定的角度阈值和角速度阈值;如果卫星三轴姿态角和角速度均小于设定的角度阈值和角速度阈值则转入步骤(4);否则转入步骤(2)。
角度阈值根据轨道控制时的姿态要求(保证一定的轨道控制方向和效率)确定,一般取为1-3度的范围;
角速度阈值根据轨道控制时所需要的姿态稳定情况和系统可容忍的角动量偏差确定,例如角速度阈值可以为0.001rad/s。
(4)计算轨道控制时间,轨道控制发动机按照所计算的时间进行轨道控制,在轨道控制结束后返回步骤(1)。
轨道控制时间计算公式如下:
t△=(Hz△-Hz-Iz*wz)/Tgz
其中:t△为本次轨道控制时长;Hz△为偏航轴动量轮正常工作的角动量阈值;Hz为当前偏航轴动量轮角动量;Iz为卫星偏航轴转动惯量;wz为卫星偏航轴角速度;Tgz为轨道控制发动机轨道控制时产生的偏航干扰力矩。
以偏航姿态控制发动机故障为例,设轨道控制发动机轨道控制时产生的偏航干扰力矩Tgz为0.5Nm,若偏航轴动量轮当前角动量Hz为4Nms,动量轮正常工作的角动量阈值Hz△为10Nms,卫星偏航轴转动惯量为1000kgm2,卫星偏航角速度为0.001rad/s,则偏航轴可吸收的干扰角动量Hz△-Hz-Iz*wz为5Nms,计算出本次轨道控制发动机可工作时长为10s。
(5)对滚动、俯仰轴进行动量轮卸载,偏航轴采用动量轮进行偏航轴的姿态控制;然后返回步骤(1)。
对滚动轴的动量轮进行卸载的过程如下:根据该滚动轴的实际角动量与标称角动量的偏差采用如PI等控制方式计算动量轮的控制力矩,并送给动量轮,以期动量轮角动量向标称角动量接近;然后,将该滚动轴的姿态角偏差及角速度偏差作为喷气相平面控制的两个输入,计算出星体的喷气控制力矩,将上述喷气控制力矩送至滚动轴的姿态控制发动机,以抑制由上述动量轮角动量变化引起的星体姿态变化。对俯仰轴的动量轮进行卸载的过程同上。
由于偏航轴无姿态控制发动机,因此偏航轴的动量轮进行正常控制,保证偏航姿态和角速度稳定。
实施例
卫星控制系统组成如图2所示:执行机构由姿态控制发动机(16个1N推力器)、轨道控制发动机(2个20N推力器)、动量轮(5个25Nms的动量轮)组成,本发明中所述的卫星轨道控制方法在姿轨控计算机上运行。
某卫星入轨后,出现偏航轴姿态控制发动机不可用,导致卫星处于姿态控制发动机欠驱动的状态。同时该卫星的轨道控制发动机轨道控制时受羽流影响,产生卫星偏航轴较大的干扰力矩,约为0.5Nm。由于载荷正常工作的需要,该卫星必需抬高轨道55km至目标轨道,因此需要卫星在偏航欠驱动、轨道控制发动机干扰力矩较大的情况下进行连续轨道控制。
按照前面所述的方法进行连续自主轨道控制,选取动量轮卸载条件判断中的角动量与标称角动量的偏差阈值为5Nms,满足轨道控制的角速度阈值为0.005度/秒,角度阈值为2度,动量轮进行三轴姿态控制选取PID控制律。经仿真,在一个轨道周期内可连续进行16次轨道控制,每次轨道控制时长为10s,一个轨道周期内进行滚动方向的动量轮卸载次数为3次,动量轮卸载过程中卫星姿态角控制在2度以内,角速度在0.1度/秒以内,姿态稳定,轨道控制阶段姿态角及角速度曲线参见图3和图4。
本发明未详细说明部分属本领域技术人员公知常识。
Claims (3)
1.一种卫星轨道控制方法,能够在偏航轴和滚动轴的其中一个轴的姿态控制发动机工作不正常情况下实现卫星的轨道控制,包括如下步骤:
(1)对姿态控制发动机工作正常的两个轴进行动量轮卸载条件判断,当满足动量轮卸载条件时,转入步骤(5);当不满足动量轮卸载条件时,转入步骤(2);
(2)采用动量轮进行三轴姿态控制;
(3)判断卫星三轴姿态角和角速度是否均小于设定的角度阈值和角速度阈值;如果卫星三轴姿态角和角速度均小于设定的角度阈值和角速度阈值则转入步骤(4);否则转入步骤(2);
(4)计算轨道控制时间,轨道控制发动机按照所计算的轨道控制时间进行轨道控制,在轨道控制结束后返回步骤(1);
(5)对姿态控制发动机工作正常的两个轴进行动量轮卸载,对姿态控制发动机工作不正常的轴采用动量轮进行姿态控制,然后返回步骤(1)。
2.如权利要求1所述的卫星轨道控制方法,其特征在于:
所述步骤(1)中,若当前角动量与标称角动量的差值大于设定的角动量偏差阈值,则满足动量轮卸载条件;若当前角动量与标称角动量的差值不大于设定的角动量偏差阈值,则不满足动量轮卸载条件。
3.如权利要求1所述的卫星轨道控制方法,其特征在于:
所述步骤(4)中的轨道控制时间根据如下公式计算:
t△=(Hz△-Hz-Iz*wz)/Tgz
其中:t△为轨道控制时间;Hz△为姿态控制发动机工作不正常的轴的动量轮的角动量阈值;Hz为姿态控制发动机工作不正常的轴的实际角动量;Iz为姿态控制发动机工作不正常的轴的转动惯量;wz为姿态控制发动机工作不正常的轴的角速度;Tgz为轨道控制发动机在姿态控制发动机工作不正常的轴上产生的干扰力矩。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013100363994A CN103112604B (zh) | 2013-01-30 | 2013-01-30 | 一种卫星轨道控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013100363994A CN103112604B (zh) | 2013-01-30 | 2013-01-30 | 一种卫星轨道控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103112604A true CN103112604A (zh) | 2013-05-22 |
CN103112604B CN103112604B (zh) | 2013-11-20 |
Family
ID=48411006
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2013100363994A Active CN103112604B (zh) | 2013-01-30 | 2013-01-30 | 一种卫星轨道控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103112604B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103332301A (zh) * | 2013-06-08 | 2013-10-02 | 北京航空航天大学 | 利用充液变惯量飞轮控制航天器姿态的方法及其执行机构 |
CN103424116A (zh) * | 2013-07-23 | 2013-12-04 | 中国西安卫星测控中心 | 一种适应轨道机动的地球同步卫星精密定轨方法 |
CN103472726A (zh) * | 2013-09-22 | 2013-12-25 | 上海无线电设备研究所 | 一种提高空间活动部件寿命的非线性跟踪控制方法 |
CN103676954A (zh) * | 2013-11-18 | 2014-03-26 | 航天东方红卫星有限公司 | 一种卫星入轨成功评判方法 |
CN103818564A (zh) * | 2014-03-04 | 2014-05-28 | 中国人民解放军国防科学技术大学 | 一种采用小推力的航天器轨道维持与对地定向姿态保持一体化控制方法 |
CN103941740A (zh) * | 2014-04-15 | 2014-07-23 | 北京控制工程研究所 | 一种考虑地球椭率的多轴机动成像卫星偏航姿态控制方法 |
CN104267735A (zh) * | 2014-09-02 | 2015-01-07 | 上海新跃仪表厂 | 一种倾斜轨道卫星偏航机动时机判断方法 |
CN107255570A (zh) * | 2017-05-18 | 2017-10-17 | 上海卫星工程研究所 | 适合于大部件翻身式展开试验的姿态调整方法 |
CN107776915A (zh) * | 2017-09-15 | 2018-03-09 | 深圳市盛路物联通讯技术有限公司 | 设备系统控制方法及相关产品 |
CN108657468A (zh) * | 2018-04-20 | 2018-10-16 | 北京控制工程研究所 | 一种具有最大角动量包络的动量轮驱动力矩分配方法 |
CN108761507A (zh) * | 2018-05-21 | 2018-11-06 | 中国人民解放军战略支援部队信息工程大学 | 基于短弧定轨和预报的导航卫星轨道快速恢复方法 |
CN108983799A (zh) * | 2018-07-24 | 2018-12-11 | 北京电子工程总体研究所 | 一种满足快速绕飞要求的航天器绕飞观测方法 |
CN109625329A (zh) * | 2018-12-10 | 2019-04-16 | 上海航天控制技术研究所 | 一种基于离散喷气的飞轮角动量自主卸载方法 |
CN113602532A (zh) * | 2021-08-25 | 2021-11-05 | 重庆零壹空间科技集团有限公司 | 一种固体运载火箭入轨修正方法 |
CN115384811A (zh) * | 2022-09-06 | 2022-11-25 | 北京控制工程研究所 | V型轮控且单轮掉电卫星的三轴姿态快速稳定控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101066706A (zh) * | 2007-07-03 | 2007-11-07 | 北京控制工程研究所 | 减小姿控喷气对轨道影响的估计与补偿方法 |
EP2181923A1 (fr) * | 2008-10-31 | 2010-05-05 | Thales | Procédé et système de désaturation des roues d'inertie d'un engin spatial |
CN102880183A (zh) * | 2012-10-15 | 2013-01-16 | 北京控制工程研究所 | 一种用于偏航机动的控制力矩陀螺角动量管理方法 |
-
2013
- 2013-01-30 CN CN2013100363994A patent/CN103112604B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101066706A (zh) * | 2007-07-03 | 2007-11-07 | 北京控制工程研究所 | 减小姿控喷气对轨道影响的估计与补偿方法 |
EP2181923A1 (fr) * | 2008-10-31 | 2010-05-05 | Thales | Procédé et système de désaturation des roues d'inertie d'un engin spatial |
CN102880183A (zh) * | 2012-10-15 | 2013-01-16 | 北京控制工程研究所 | 一种用于偏航机动的控制力矩陀螺角动量管理方法 |
Non-Patent Citations (1)
Title |
---|
王振枫: "《某型微小卫星姿态控制系统设计及仿真研究》", 《南京航空航天大学》 * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103332301A (zh) * | 2013-06-08 | 2013-10-02 | 北京航空航天大学 | 利用充液变惯量飞轮控制航天器姿态的方法及其执行机构 |
CN103332301B (zh) * | 2013-06-08 | 2015-09-23 | 北京航空航天大学 | 利用充液变惯量飞轮控制航天器姿态的方法及其执行机构 |
CN103424116A (zh) * | 2013-07-23 | 2013-12-04 | 中国西安卫星测控中心 | 一种适应轨道机动的地球同步卫星精密定轨方法 |
CN103424116B (zh) * | 2013-07-23 | 2015-09-23 | 中国西安卫星测控中心 | 一种适应轨道机动的地球同步卫星精密定轨方法 |
CN103472726A (zh) * | 2013-09-22 | 2013-12-25 | 上海无线电设备研究所 | 一种提高空间活动部件寿命的非线性跟踪控制方法 |
CN103676954A (zh) * | 2013-11-18 | 2014-03-26 | 航天东方红卫星有限公司 | 一种卫星入轨成功评判方法 |
CN103676954B (zh) * | 2013-11-18 | 2016-07-06 | 航天东方红卫星有限公司 | 一种卫星入轨成功评判方法 |
CN103818564A (zh) * | 2014-03-04 | 2014-05-28 | 中国人民解放军国防科学技术大学 | 一种采用小推力的航天器轨道维持与对地定向姿态保持一体化控制方法 |
CN103818564B (zh) * | 2014-03-04 | 2015-11-25 | 中国人民解放军国防科学技术大学 | 一种采用小推力的航天器轨道维持与对地定向姿态保持一体化控制方法 |
CN103941740B (zh) * | 2014-04-15 | 2016-07-06 | 北京控制工程研究所 | 一种考虑地球椭率的多轴机动成像卫星偏航姿态控制方法 |
CN103941740A (zh) * | 2014-04-15 | 2014-07-23 | 北京控制工程研究所 | 一种考虑地球椭率的多轴机动成像卫星偏航姿态控制方法 |
CN104267735B (zh) * | 2014-09-02 | 2017-06-06 | 上海新跃仪表厂 | 一种倾斜轨道卫星偏航机动时机判断方法 |
CN104267735A (zh) * | 2014-09-02 | 2015-01-07 | 上海新跃仪表厂 | 一种倾斜轨道卫星偏航机动时机判断方法 |
CN107255570A (zh) * | 2017-05-18 | 2017-10-17 | 上海卫星工程研究所 | 适合于大部件翻身式展开试验的姿态调整方法 |
CN107776915A (zh) * | 2017-09-15 | 2018-03-09 | 深圳市盛路物联通讯技术有限公司 | 设备系统控制方法及相关产品 |
CN107776915B (zh) * | 2017-09-15 | 2020-06-16 | 深圳市盛路物联通讯技术有限公司 | 应用于太空设备的控制方法及装置 |
CN108657468B (zh) * | 2018-04-20 | 2020-08-14 | 北京控制工程研究所 | 一种具有最大角动量包络的动量轮驱动力矩分配方法 |
CN108657468A (zh) * | 2018-04-20 | 2018-10-16 | 北京控制工程研究所 | 一种具有最大角动量包络的动量轮驱动力矩分配方法 |
CN108761507A (zh) * | 2018-05-21 | 2018-11-06 | 中国人民解放军战略支援部队信息工程大学 | 基于短弧定轨和预报的导航卫星轨道快速恢复方法 |
CN108983799A (zh) * | 2018-07-24 | 2018-12-11 | 北京电子工程总体研究所 | 一种满足快速绕飞要求的航天器绕飞观测方法 |
CN109625329A (zh) * | 2018-12-10 | 2019-04-16 | 上海航天控制技术研究所 | 一种基于离散喷气的飞轮角动量自主卸载方法 |
CN109625329B (zh) * | 2018-12-10 | 2022-06-24 | 上海航天控制技术研究所 | 一种基于离散喷气的飞轮角动量自主卸载方法 |
CN113602532A (zh) * | 2021-08-25 | 2021-11-05 | 重庆零壹空间科技集团有限公司 | 一种固体运载火箭入轨修正方法 |
CN115384811A (zh) * | 2022-09-06 | 2022-11-25 | 北京控制工程研究所 | V型轮控且单轮掉电卫星的三轴姿态快速稳定控制方法 |
CN115384811B (zh) * | 2022-09-06 | 2024-05-14 | 北京控制工程研究所 | V型轮控且单轮掉电卫星的三轴姿态快速稳定控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103112604B (zh) | 2013-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103112604B (zh) | 一种卫星轨道控制方法 | |
CN103072702B (zh) | 卫星轨道和姿态控制方法 | |
CN105539884B (zh) | 一种卫星偏航控制导引方法 | |
CN106660641A (zh) | 控制卫星在地球轨道中的轨道的方法、卫星和控制这种卫星的轨道的系统 | |
CN103112602B (zh) | 一种基于推力器和动量轮联合的三轴角速度阻尼控制方法 | |
CN104880946A (zh) | 一种基于鲁棒预见控制的舰载机自动着舰控制方法 | |
US20180194495A1 (en) | Energy efficient satellite maneuvering | |
CN104181941A (zh) | 一种适应倾斜轨道卫星的帆板双向控制方法 | |
CN103955224B (zh) | 一种用于相对运动视线跟踪的姿态控制方法 | |
JPH06135396A (ja) | 宇宙船の姿勢制御のための補償転換 | |
CN109625329B (zh) | 一种基于离散喷气的飞轮角动量自主卸载方法 | |
CN103963966A (zh) | 用于滑行的飞行器的侧向控制的方法和装置 | |
CN102880183A (zh) | 一种用于偏航机动的控制力矩陀螺角动量管理方法 | |
CN104058104A (zh) | 无加速度计情况下一种基于关调制的高精度轨控方法 | |
CN104567545A (zh) | Rlv大气层内主动段的制导方法 | |
CN103076809B (zh) | 一种利用干扰积累角动量自平衡的轨道控制方法 | |
CN103818564B (zh) | 一种采用小推力的航天器轨道维持与对地定向姿态保持一体化控制方法 | |
US20140138491A1 (en) | Spacecraft Momentum Unload and Station-keeping Techniques | |
CN103253382B (zh) | 一种高精度发动机联合变轨方法 | |
CN102880059B (zh) | 一种采用正弦偏航导引律的偏航机动控制方法 | |
US20150375854A1 (en) | Differential steering control of electric taxi landing gear | |
CN108839824B (zh) | 一种基于合作博弈的混合执行机构动量优化管理方法 | |
CN109032158B (zh) | 一种直线规划的头对日—尾对日交替连续偏航姿态控制方法 | |
EP2390189B1 (en) | Methods and systems for reducing angular velocity using a gyroscope array | |
CN113569391A (zh) | 地月转移轨道的参数确定方法、装置、设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |