CN108128484A - 一种基于线性二次型调节器的双星系统轨道保持方法 - Google Patents
一种基于线性二次型调节器的双星系统轨道保持方法 Download PDFInfo
- Publication number
- CN108128484A CN108128484A CN201711366153.8A CN201711366153A CN108128484A CN 108128484 A CN108128484 A CN 108128484A CN 201711366153 A CN201711366153 A CN 201711366153A CN 108128484 A CN108128484 A CN 108128484A
- Authority
- CN
- China
- Prior art keywords
- binary
- track
- control
- star system
- star
- 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
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Chemical & Material Sciences (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Combustion & Propulsion (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及一种双星系统的轨道保持方法,涉及一种采用线性二次型调节器跟踪控制实现双星系统轨道保持的方法,属于航空航天技术领域。本发明实现方法如下:利用简化双体动力学模型设计标称轨道,然后在精确的双体动力学方程下以标称轨道的初值得真实轨道,计算真实轨道与标称轨道的偏差,基于线性二次型调节器设计最优控制律,获得最优加速度,施加连续控制,使真实轨道收敛在标称轨道附近,实现双星系统下的轨道跟踪和保持。本发明适用于存在初始误差和模型误差情况下的双星系统轨道保持,具有收敛性好,保持精度高等特点。
Description
技术领域
本发明涉及一种双星系统的轨道保持方法,特别涉及一种采用线性二次型调节器跟踪控制实现双星系统轨道保持的方法,属于航空航天技术领域。
背景技术
双小行星系统是特殊的一类小行星,双星系统质心在围绕太阳进行公转的同时,两颗小行星还围绕质心进行圆形或椭圆运动。相比单小行星探测,双小行星系统探测具有更多的科学回报,同时双星系统的形成和演化仍然存在争论,因此对双星系统开展探测是未来小行星探测的热点之一。与单小行星相同,双小行星系统也存在形状不规则,自旋状态不确定等特点,导致双星系统内的动力学环境复杂,而双星之间的相对运动使探测器的运动非线性混沌线性明显,基于简化动力学设计的轨道通常无法稳定运行,给探测轨道带来了困难,需要通过轨道控制实现轨道维持。
目前双星系统的轨道保持方法中,在先技术【1】(参见Loic Chappaz,KathleenC.Howell.Trajectory design for bounded motion near uncertain binary systemscomprised of small irregular bodies exploiting sliding control modes[J].ActaAstronautica,2015,115,226-240)采用高阶滑模制导的方法对有界轨道进行控制,将标称轨道的若干点作为目标点,并将实际轨道分为无控段和受控段,在受控段利用滑模制导使轨道逼近目标的从而实现轨道保持。但该方法由于仅约束末端状态,实际轨道在除目标的附近的区域外可能和标称轨道差异较大,跟踪效果较差。
在先技术【2】(参见Pamela Woo,Arun K.Misra等,Bounded trajectories of aspacecraft near an equilibrium point of a binary asteroid system,2015,110,313-323)通过构建合适的Lyapunov函数,实现探测器在双星系统平衡点附近轨道的反馈控制,但该方法仅能应用与平衡点附近轨道,不适合双星系统内的其他轨道。
发明内容
本发明公开的一种基于线性二次型调节器的双星系统轨道保持方法要解决的技术问题是,提出一种适用于双星系统周期轨道的轨道保持方法,基于线性二次型调节器的控制律施加连续控制,实现真实轨道收敛在标称轨道附近,能够适用于存在初始误差和模型误差情况下的双星系统轨道保持,具有收敛性好,保持精度高等特点。
所述真实轨道指考虑初始误差和精确引力场模型后的轨道。
本发明的目的是通过下述技术方案实现的。
本发明一种基于线性二次型调节器的双星系统轨道保持方法,包括以下步骤:
步骤一:根据确定的目标双星系统参数,通过忽略双体的自旋角速度差别,将小行星简化为三轴椭球体,在质心旋转坐标系下建立建立如公式(1)所示的简化双体动力学模型。所述的双星系统参数包括双星的形状,质量比μ。
其中r=[x,y,z]表示探测器到质心的位置矢量,ω为系统的旋转角速度。U1,2表示双星的引力势能
其中α,β,γ为椭球体的三个主轴长度,初状态下星体的最长主轴与X轴重合。
r1,r2分别表示探测器到双星系统两颗天体的位置矢量,μ为质量比。
步骤二:基于步骤一简化的双体动力学模型进行标称轨迹设计。
在简化的双体动力学模型下设计标称轨道,选择周期轨道作为标称轨道,利用双体动力学模型的对称性,周期轨道应满足 面对称轨道满足初状态且半个轨道周期后的状态为轴对称的轨道满足初状态半个轨道周期后的状态满足根据已有的搜索方法对系统内的区域进行搜索得到初值,利用微分修正方法得到精确的周期轨道作为标称轨道Γ0。
步骤二所述的已有的搜索方法包括利用速度庞加莱映射或网格搜索。
步骤三:考虑双星系统的自旋速度Ω和精确引力场模型,建立精确的双体动力学方程。
令双星系统的两个天体自旋速度分别为Ω1,Ω2,初始状态天体固连系与质心旋转系的夹角分别为θ1,θ2,则在质心旋转系下的两个天体的自旋速度分别为ω1′=Ω1-ω0,ω2′=Ω2-ω0,自旋角随时间的变化为φ1=θ1+ω1′t,φ2=θ2+ω2′t,则动力学方程变为
其中u=[acx,acy,acz]为控制加速度,Rz为绕Z轴的旋转矩阵。
步骤四:在精确的双体动力学方程下以标称轨道Γ0的初值作为初始值考虑误差进行动力学积分,得真实轨道Γ。计算真实轨道与标称轨道的偏差,根据线性二次型调节器设计最优控制律。
在精确的双体动力学方程下以标称轨道Γ0的初值作为初始值考虑误差进行动力学积分时间T0,得真实轨道Γ。记真实轨道Γ的状态为标称轨道的状态为记真实轨道和标称轨道的误差为e=X1-X0。
根据双体动力学方程(1)得
其中
控制问题表述为通过控制使误差e消除为0,考虑在[t0,tf]时间内,使跟踪误差和消耗的燃料最优,将问题转化为最优控制问题,取二次型性能指标函数为
式中Q为6×6正定常矩阵,R为3×3正定常矩阵,分别为轨道误差e和控制加速度u的加权系数。Q主要决定控制精度和动态效果,R主要决定控制的能量消耗,对于控制成本较大的情况,选择较大的R和较小的Q;控制成本较低的情况选择较大的Q和较小的R以提高控制精度。由于探测器携带的燃料有限,需考虑控制成本选择Q和R。优选Q=diag(1,1,1,1,1,1),R=diag(1010,1010,1010)。
反馈增益矩阵为
K=R-1BTP (10)
式中P通过求解黎卡提代数方程得到
PA+ΑTP-PBR-1BTP+Q=0 (11)
从而得到最优控制律为
u=-Ke=-K(X1-X0) (12)
根据最优控制律得到最优控制加速度u代入方程(6),继续积分。
步骤五选择积分时长T0,重复步骤四,计算真实轨道Γ与标称轨道Γ0的误差并根据线性二次型调节器最优控制律得到新的最优控制加速度u,再次积分直至任务结束,保证真实轨道始终在标称轨道附近,实现双星系统下的轨道跟踪和保持。
本发明公开的一种基于线性二次型调节器的双星系统轨道保持方法,利用简化双体动力学模型设计标称轨道,然后在精确的双体动力学方程下以标称轨道的初值得真实轨道Γ,计算真实轨道与标称轨道的偏差,基于线性二次型调节器设计最优控制律,获得最优加速度,施加连续控制,使真实轨道收敛在标称轨道附近,实现双星系统下的轨道跟踪和保持。本发明能够适用于存在初始误差和模型误差情况下的双星系统轨道设计,具有收敛性好,保持精度高等特点。
所述真实轨道指考虑初始误差和精确引力场模型后的轨道。
有益效果:
1、本发明公开的一种基于线性二次型调节器的双星系统轨道保持方法,以简化的双体动力学模型计算得到的标称轨道,通过计算真实轨道与标称轨道的误差施加控制加速度,实现在考虑初始误差和引力场模型的高精度动力学下实现双体系统轨道保持,计算效率高,适用于考虑初始误差和模型误差的精确动力学模型下的轨道保持。
2、本发明公开的一种基于线性二次型调节器的双星系统轨道保持方法,通过推导合适的线性二次型控制律调节器生成最优控制律,综合考虑轨道误差精度和能量消耗,具有保持精度高,收敛性好,燃料消耗低的优点。
附图说明
图1本发明的一种基于线性二次型调节器的双星系统轨道保持方法流程示意图;
图2本发明的一种基于线性二次型调节器的双星系统轨道保持方法步骤二简化动力学模型下的标称轨道。
图3本发明的一种基于线性二次型调节器的双星系统轨道保持方法步骤四真实动力学模型下的跟踪轨道。
图4本发明的一种基于线性二次型调节器的双星系统轨道保持方法步骤四真实轨道的位置跟踪误差。
图5本发明的一种基于线性二次型调节器的双星系统轨道保持方法步骤四真实轨道的速度跟踪误差。
图6本发明的一种基于线性二次型调节器的双星系统轨道保持方法步骤四真实轨道的控制加速度变化。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合附图和实例对发明内容做进一步说明。
实施例1:
本实施例的基于线性二次型调节器的双星系统轨道保持方法,包括以下步骤:
步骤一:根据确定的目标双星系统参数,通过忽略双体的自旋角速度差别,将小行星简化为三轴椭球体,在质心旋转坐标系下建立简化双体动力学模型。所述的双星系统参数包括双星的形状,质量比μ。
其中r=[x,y,z]表示探测器到质心的位置矢量,ω为系统的旋转角速度,。U1,2表示双星的引力势能
其中α,β,γ为椭球体的三个主轴长度,假设初状态下星体的最长主轴与X轴重合。
r1,r2分别表示探测器到双星系统两颗天体的位置矢量。
以双星系统Hermes为例,简化模型中主星的三轴长分别为0.6124km,0.6km,0.5880km;子星的三轴长分别为0.5530km,0.5364km,0.5300km,系统的质量比μ=0.4083,轨道周期13.8936h,轨道角速度ω=1.2562×10-4rad/s。
步骤二:基于步骤一简化的双体动力学模型进行标称轨迹设计。
在简化的双体动力学模型下设计标称轨道,选择周期轨道作为标称轨道,利用双体动力学模型的对称性,周期轨道应满足 面对称轨道满足初状态且半个轨道周期后的状态为轴对称的轨道满足初状态半个轨道周期后的状态满足根据已有的搜索方法对系统内的区域进行搜索得到初值,利用微分修正方法得到精确的周期轨道作为标称轨道Γ0。
图2给出一组利用网格搜索得到的标称轨道所示,轨道周期为11.0067h。
步骤三:考虑双星系统的自旋速度Ω和精确引力场模型,建立精确的双体动力学方程。
取小行星的自旋周期分别为6.9468h和12.0545h,初始的自旋角均为θ1=0°,θ2=0°,则在质心旋转系下的自旋速度分别为ω′1=ω1-ω0=1.2562rad/s,ω′2=ω2-ω0=1.9167×10-5rad/s,自旋角随时间的变化为φ1=θ1+ω′1t,φ2=θ2+ω′2t,则动力学方程变为
其中u=[acx,acy,acz]为控制加速度,Rz为绕Z轴的旋转矩阵。
步骤四:在精确的双体动力学方程下以标称轨道Γ0的初值作为初始值考虑误差进行动力学积分,得真实轨道Γ。计算真实轨道与标称轨道的偏差,根据线性二次型调节器设计最优控制律。
在真实动力学下以标称轨道Γ0的初值作为初始值,考虑初始位置误差为10m,速度误差为0.001m/s进行动力学积分,积分时间T=400s,得真实轨道Γ。
令真实轨道的状态为标称轨道的状态为记真实轨道和标称轨道的误差为e=X1-X0
根据动力学方程(1)可得
其中
u=[acx,acy,acz]表示控制加速度。
控制问题可以表述为通过控制使误差e消除为0,考虑在[t0,tf]时间内,使跟踪误差和消耗的燃料最优,因此可以考虑为最优控制问题,取二次型性能指标函数为
式中Q为6×6正定常矩阵,R为3×3正定常矩阵,分别为轨道误差e和控制加速度u的加权系数。Q主要决定控制精度和动态效果,R主要决定控制的能量消耗,对于控制成本较大的情况,选择较大的R和较小的Q;控制成本较低的情况选择较大的Q和较小的R以提高控制精度。由于探测器携带的燃料有限,需考虑控制成本选择Q和R。优选Q=diag(1,1,1,1,1,1),R=diag(1010,1010,1010)。
反馈增益矩阵为
K=R-1BTP (10)
式中P通过求解黎卡提代数方程得到
PA+ATP-PBR-1BTP+Q=0 (11)
从而得到最优控率为
u=-Ke=-K(X1-X0) (12)
根据最优控制律得到最优控制加速度u代入方程(6),继续积分。
步骤五选择积分时长T=400s,重复步骤四,计算真实轨道Γ与标称轨道Γ0的误差并根据线性二次型调节器最优控制律得到新的最优控制加速度u,再次积分直至任务结束,保证真实轨道始终在标称轨道附近,实现双星系统下的轨道跟踪和保持。
图3给出了真实动力学模型下的跟踪轨道。图4和图5分别给出轨道的位置和速度跟踪误差,图6给出了真实轨道的控制加速度变化。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于线性二次型调节器的双星系统轨道保持方法,其特征在于:包括以下步骤,
步骤一:根据确定的目标双星系统参数,通过忽略双体的自旋角速度差别,将小行星简化为三轴椭球体,在质心旋转坐标系下建立如公式(1)所示的简化双体动力学模型;所述的双星系统参数包括双星的形状,质量比μ;
其中r=[x,y,z]表示探测器到质心的位置矢量,ω为系统的旋转角速度;U1,2表示双星的引力势能
其中α,β,γ为椭球体的三个主轴长度,初状态下星体的最长主轴与X轴重合;
r1,r2分别表示探测器到双星系统两颗天体的位置矢量,μ为质量比;
步骤二:基于步骤一简化的双体动力学模型进行标称轨迹设计;
在简化的双体动力学模型下设计标称轨道,选择周期轨道作为标称轨道,利用双体动力学模型的对称性,周期轨道应满足 面对称轨道满足初状态且半个轨道周期后的状态为轴对称的轨道满足初状态半个轨道周期后的状态满足根据已有的搜索方法对系统内的区域进行搜索得到初值,利用微分修正方法得到精确的周期轨道作为标称轨道Γ0;
步骤三:考虑双星系统的自旋速度Ω和精确引力场模型,建立精确的双体动力学方程;
令双星系统的两个天体自旋速度分别为Ω1,Ω2,初始状态天体固连系与质心旋转系的夹角分别为θ1,θ2,则在质心旋转系下的两个天体的自旋速度分别为ω′1=Ω1-ω0,ω′2=Ω2-ω0,自旋角随时间的变化为φ1=θ1+ω′1t,φ2=θ2+ω′2t,则动力学方程变为
其中u=[acx,acy,acz]为控制加速度,Rz为绕Z轴的旋转矩阵;
步骤四:在精确的双体动力学方程下以标称轨道Γ0的初值作为初始值考虑误差进行动力学积分,得真实轨道Γ;计算真实轨道与标称轨道的偏差,根据线性二次型调节器设计最优控制律;
在精确的双体动力学方程下以标称轨道Γ0的初值作为初始值考虑误差进行动力学积分时间T0,得真实轨道Γ;记真实轨道Γ的状态为标称轨道的状态为记真实轨道和标称轨道的误差为e=X1-X0;
根据双体动力学方程(1)得
其中K=Urr+C,
控制问题表述为通过控制使误差e消除为0,考虑在[t0,tf]时间内,使跟踪误差和消耗的燃料最优,将问题转化为最优控制问题,取二次型性能指标函数为
式中Q为6×6正定常矩阵,R为3×3正定常矩阵,分别为轨道误差e和控制加速度u的加权系数;Q主要决定控制精度和动态效果,R主要决定控制的能量消耗,对于控制成本较大的情况,选择较大的R和较小的Q;控制成本较低的情况选择较大的Q和较小的R以提高控制精度;由于探测器携带的燃料有限,需考虑控制成本选择Q和R;
反馈增益矩阵为
K=R-1BTP (10)
式中P通过求解黎卡提代数方程得到
PA+ΑTP-PBR-1BTP+Q=0 (11)
从而得到最优控制律为
u=-Ke=-K(X1-X0) (12)
根据最优控制律得到最优控制加速度u代入方程(6),继续积分;
步骤五选择积分时长T0,重复步骤四,计算真实轨道Γ与标称轨道Γ0的误差并根据线性二次型调节器最优控制律得到新的最优控制加速度u,再次积分直至任务结束,保证真实轨道始终在标称轨道附近,实现双星系统下的轨道跟踪和保持。
2.如权利要求1所述的一种基于线性二次型调节器的双星系统轨道保持方法,其特征在于:步骤二所述的已有的搜索方法包括利用速度庞加莱映射或网格搜索。
3.如权利要求1或2所述的一种基于线性二次型调节器的双星系统轨道保持方法,其特征在于:步骤四中选Q=diag(1,1,1,1,1,1),R=diag(1010,1010,1010)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711366153.8A CN108128484B (zh) | 2017-12-18 | 2017-12-18 | 一种基于线性二次型调节器的双星系统轨道保持方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711366153.8A CN108128484B (zh) | 2017-12-18 | 2017-12-18 | 一种基于线性二次型调节器的双星系统轨道保持方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108128484A true CN108128484A (zh) | 2018-06-08 |
CN108128484B CN108128484B (zh) | 2020-08-28 |
Family
ID=62390604
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711366153.8A Active CN108128484B (zh) | 2017-12-18 | 2017-12-18 | 一种基于线性二次型调节器的双星系统轨道保持方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108128484B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110405764A (zh) * | 2019-07-25 | 2019-11-05 | 北京理工大学 | 基于仿生眼机器人的轨迹追踪方法及装置、仿生眼、仿生眼机器人 |
CN111731513A (zh) * | 2020-06-15 | 2020-10-02 | 航天东方红卫星有限公司 | 一种基于单脉冲轨控的高精度引力场中回归轨道维持方法 |
CN112859877A (zh) * | 2021-02-01 | 2021-05-28 | 上海同普电力技术有限公司 | 一种agv四舵轮同步控制算法 |
US11614719B2 (en) | 2019-07-25 | 2023-03-28 | Beijing Institute Of Technology | Wide-field-of-view anti-shake high-dynamic bionic eye |
CN116817932A (zh) * | 2022-12-13 | 2023-09-29 | 北京航空航天大学 | 一种轨道维持、轨道确定和重力场实时测绘的一体化方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103257352A (zh) * | 2013-05-10 | 2013-08-21 | 航天东方红卫星有限公司 | 一种双星gps闭环控制测试方法 |
CN106338296A (zh) * | 2016-11-04 | 2017-01-18 | 上海航天控制技术研究所 | 一种双星敏感器在轨实时交互的修正方法 |
CN106597856A (zh) * | 2016-12-29 | 2017-04-26 | 北京理工大学 | 基于网格搜索的双星系统空间轨道族搜索方法 |
CN106650278A (zh) * | 2016-12-29 | 2017-05-10 | 北京理工大学 | 非同步双星系统的稳定轨道计算方法 |
CN106682274A (zh) * | 2016-12-05 | 2017-05-17 | 北京理工大学 | 考虑振幅约束的一种Halo轨道在轨保持方法 |
-
2017
- 2017-12-18 CN CN201711366153.8A patent/CN108128484B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103257352A (zh) * | 2013-05-10 | 2013-08-21 | 航天东方红卫星有限公司 | 一种双星gps闭环控制测试方法 |
CN106338296A (zh) * | 2016-11-04 | 2017-01-18 | 上海航天控制技术研究所 | 一种双星敏感器在轨实时交互的修正方法 |
CN106682274A (zh) * | 2016-12-05 | 2017-05-17 | 北京理工大学 | 考虑振幅约束的一种Halo轨道在轨保持方法 |
CN106597856A (zh) * | 2016-12-29 | 2017-04-26 | 北京理工大学 | 基于网格搜索的双星系统空间轨道族搜索方法 |
CN106650278A (zh) * | 2016-12-29 | 2017-05-10 | 北京理工大学 | 非同步双星系统的稳定轨道计算方法 |
Non-Patent Citations (1)
Title |
---|
杏建军等: "基于改进线性二次型调节器的近地轨道编队卫星鲁棒控制", 《国防科技大学学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110405764A (zh) * | 2019-07-25 | 2019-11-05 | 北京理工大学 | 基于仿生眼机器人的轨迹追踪方法及装置、仿生眼、仿生眼机器人 |
US11614719B2 (en) | 2019-07-25 | 2023-03-28 | Beijing Institute Of Technology | Wide-field-of-view anti-shake high-dynamic bionic eye |
CN111731513A (zh) * | 2020-06-15 | 2020-10-02 | 航天东方红卫星有限公司 | 一种基于单脉冲轨控的高精度引力场中回归轨道维持方法 |
CN112859877A (zh) * | 2021-02-01 | 2021-05-28 | 上海同普电力技术有限公司 | 一种agv四舵轮同步控制算法 |
CN116817932A (zh) * | 2022-12-13 | 2023-09-29 | 北京航空航天大学 | 一种轨道维持、轨道确定和重力场实时测绘的一体化方法 |
CN116817932B (zh) * | 2022-12-13 | 2024-03-15 | 北京航空航天大学 | 一种轨道维持、轨道确定和重力场实时测绘的一体化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108128484B (zh) | 2020-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108128484A (zh) | 一种基于线性二次型调节器的双星系统轨道保持方法 | |
CN106682274B (zh) | 考虑振幅约束的一种Halo轨道在轨保持方法 | |
CN106672266B (zh) | 考虑时间约束的平衡点Halo轨道调相轨道转移方法 | |
CN109283934B (zh) | 基于旋转路径质量的航天器多约束姿态机动优化方法 | |
CN106679674B (zh) | 基于星历模型的地月L2点Halo轨道阴影分析方法 | |
CN111306989A (zh) | 一种基于平稳滑翔弹道解析解的高超声速再入制导方法 | |
CN109911249B (zh) | 低推重比飞行器的星际转移有限推力入轨迭代制导方法 | |
WO2021063073A1 (zh) | 一种指定发射仰角的自由弹道构造方法 | |
CN104898642A (zh) | 一种用于航天器姿态控制算法的集成测试仿真系统 | |
CN105574261A (zh) | 一种月球借力约束的地月平动点转移轨道设计方法 | |
CN109539903A (zh) | 一种固体运载火箭椭圆转移轨道迭代制导控制方法 | |
CN109033604B (zh) | 含旋转载荷的卫星动力学建模及轴承处受力的确定方法 | |
CN112572835B (zh) | 一种具有姿态切换的卫星在轨角动量管理及控制方法 | |
CN113343442B (zh) | 一种求解固定时间有限燃料多脉冲转移轨道的方法及系统 | |
CN107589756A (zh) | 一种奔月卫星编队初始化方法 | |
CN104249816A (zh) | 非合作目标绕飞悬停的姿轨协同控制方法 | |
CN105253329B (zh) | 一种基于弱稳定边界的两脉冲行星捕获轨道方法 | |
CN105930305B (zh) | 一种三脉冲交会接近制导方法 | |
CN113815903B (zh) | 一种用于遥感卫星的飞轮过零规避方法 | |
CN107367941A (zh) | 基于非线性增益的高超声速飞行器攻角观测方法 | |
CN113093776A (zh) | 一种航天器的离轨参数确定方法及装置 | |
CN116331525B (zh) | 一种卫星飞轮转速过零规避方法 | |
CN106092096A (zh) | 高精度轨道仿真中基于迭代逼近方法的卫星位置确定方法 | |
CN115113638A (zh) | 一种燃料最优主动漂移三维成像轨道控制方法 | |
CN111897355B (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 |