CN103853047B - 一种基于状态量反馈的小推力跟踪制导方法 - Google Patents
一种基于状态量反馈的小推力跟踪制导方法 Download PDFInfo
- Publication number
- CN103853047B CN103853047B CN201410073905.1A CN201410073905A CN103853047B CN 103853047 B CN103853047 B CN 103853047B CN 201410073905 A CN201410073905 A CN 201410073905A CN 103853047 B CN103853047 B CN 103853047B
- Authority
- CN
- China
- Prior art keywords
- orbit
- state
- time
- msubsup
- optimal
- 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.)
- Expired - Fee Related
Links
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种基于状态量反馈的小推力跟踪制导方法,该方法包括有构建航天器从初始轨道转移到目标轨道的时间最优转移轨道步骤;在时间最优转移轨道基础上设计时间最优状态量参考轨道的步骤;构造单圈内最优控制问题,设计出单圈控制律的步骤。本发明方法是将多圈的连续推力轨道机动优化问题转化为一系列单圈的状态量参考轨道跟踪控制问题来解决。通过构造与时间无关的状态量参考轨道,设计出单圈控制律,从而实现状态量参考轨道的跟踪。
Description
技术领域
本发明属于航天飞行器轨道控制技术领域,尤其是涉及一种基于状态量的小推力跟踪制导方法。
背景技术
小推力离子发动机相较于传统化学推进发动机来,以其高比冲、长寿命、高可靠性的特点,在星际航行、轨道提升和轨道转移、地球静止轨道卫星的站位保持等领域有着广泛的应用前景。
与小推力行星际轨道转移不同,地球静止轨道转移通常包含上百圈的中途转移轨道,对最优转移轨道设计和制导都提出了很大的挑战。最优转移轨道设计通常是以时间最优或燃料最优为优化目标,通过直接法或间接法求解得到开环的最优控制律。如果直接将开环最优控制律作用于航天器,由于航天器控制系统本身及外界环境的不确定性干扰,将导致航天器偏离原来设计的最优转移轨道,即便在设计最优控制律时考虑精细的动力学模型,但实际过程中必然会遇到无法预知的扰动,依然无法保证最后能到达预定的目标轨道。对于长时间多圈的地球静止轨道射入问题来讲,该问题尤甚,设计一套反馈的制导方法尤为必要。
传统的跟踪制导方法都是基于时间相关的,即在设计参考轨道阶段充分考虑推力发动机的冗余,在制导跟踪阶段跟踪特定时间下参考轨道参数的变化量;传统跟踪控制思路简单,未全面考虑控制能力约束,且无法保证最优性。
发明内容
为了解决小推力长时间多圈的地球静止轨道射入问题,本发明提出了一种基于状态量反馈的小推力跟踪制导方法。该方法是将多圈的连续推力轨道机动优化问题转化为一系列单圈的参考轨道跟踪控制问题来解决。通过构造与时间无关的状态量参考轨道,设计出单圈控制律,从而实现状态量参考轨道的跟踪。
本发明是一种基于状态量反馈的小推力跟踪制导方法,该方法包括有构建航天器从初始轨道转移到目标轨道的时间最优转移轨道步骤;在时间最优转移轨道基础上设计时间最优状态量参考轨道的步骤;构造单圈内最优控制问题,设计出单圈控制律的步骤。
本发明基于状态量反馈的小推力跟踪制导方法的优点在于:
①本发明方法消除了传统跟踪制导方法与时间的相关性,构造了与轨道状态量相关的跟踪参考轨道。
②本发明方法将多圈的连续推力轨道转移问题转换成一系列单圈的状态量参考轨道跟踪问题,从而实现了闭环的轨道跟踪控制。
③本发明通过跟踪状态量参考轨道,能够充分考虑推进器的变轨能力,保证了跟踪方案的最优性。
④本发明在求解单圈跟踪控制问题时,能够通过多种方法找到最优控制律,具有跟踪状态量参考轨道多元性和单圈最优控制设计多样性的特点。
附图说明
图1是航天器飞行动力学原理中对轨道要素定义的示意图。
图1A是卫星轨道姿态动力学与控制中对轨道要素定义的示意图。
图2A是本发明中轨道伪半长轴与状态量参考轨道示意图。
图2B是本发明中轨道倾角与状态量参考轨道示意图。
图2C是本发明中轨道偏心率与状态量参考轨道示意图。
图3是应用本发明基于状态量反馈的小推力跟踪制导进行轨道转移的流程图。
具体实施方式
下面将结合附图和实施例对本发明做进一步的详细说明。
在1995年12月第1版《航天器飞行动力学原理》,肖业伦编著,第44页的图2-13轨道要素定义中,一般卫星轨道包括有六个要素,即轨道偏心率e、轨道半长轴a、近地点幅角ω、纬度幅角u、轨道倾角i和轨道升交点赤经Ω。
在本发明中,在时间最优的参考轨道O下的轨道转移任务时间记为T最优,在实际转移轨道Q下的轨道转移任务时间记为T实际,轨道转移的任意一时刻记为ts。
实际转移轨道记为Q,则构建所述Q所涉及到的轨道要素分别记为轨道半长轴aQ、偏心率eQ、轨道倾角iQ、升交点赤经ΩQ、近地点幅角ωQ和真近点角θQ;任意一个时刻ts下的实际转移轨道为Q(ts)={aQ,eQ,iQ,ΩQ,ωQ,θQ},实际转移轨道的状态量的横轴为T实际。
参见图3所示,应用本发明提出的基于状态量反馈的小推力跟踪制导方法进行轨道转移的执行步骤为:
(A)依据时间最优的参考轨道TO,获取状态量参考轨道TR;
(B)依据状态量参考轨道TR设计单圈的最优控制律uk;k表示圈数;
(C)航天飞行器按照单圈的最优控制律u进行轨道转移;
(D)通过当前时刻的实际轨道状态信息QS和状态量参考轨道TR计算下一圈的最优控制律uk+1;
(E)航天飞行器按照uk+1进行轨道转移;
(F)重复(D)和(E)步骤,直至轨道状态设计量为零,完成飞行任务。
本发明是一种基于状态量反馈的小推力跟踪制导方法,该方法包括有下列步骤:
步骤1:构建航天器从初始轨道转移到目标轨道的时间最优转移轨道;
参见图1、图1A所示,航天器的轨道要素通常包括有轨道半长轴a(单位为千米)、偏心率e(单位为无量纲)、轨道倾角i(单位为度)、升交点赤经Ω(单位为度)、近地点幅角ω(单位为度)和真近点角θ(单位为度)。轨道半长轴a和偏心率e的定义参见1998年8月第1版《卫星轨道姿态动力学与控制》,章仁为编著,第4页图1.1-4相关内容。
在本发明中,目标轨道则为地球静止轨道E,所述目标轨道轨道要素包括有轨道半长轴aE、偏心率eE、轨道倾角iE、升交点赤经ΩE、近地点幅角ωE和真近点角θE;由于不考虑所述目标轨道的定点问题,则ΩE、ωE和θE的变量取值自由,无约束;而aE取值为42165km,eE取值为0,iE取值为0度。
初始轨道I的轨道要素包括有轨道半长轴aI、偏心率eI、轨道倾角iI、升交点赤经ΩI、近地点幅角ωI和真近点角θI。
在小推力控制律u作用下,航天器从初始轨道I转移到目标轨道(即地球静止轨道E)上,并且转移时间最短(单位为秒)的转移轨道记为时间最优转移轨道TO。所述时间最优转移轨道TO所涉及到的轨道要素分别是轨道半长轴aTO、偏心率eTO、轨道倾角iTO、升交点赤经ΩTO、近地点幅角ωTO和真近点角θTO;任意一个时刻ts下的时间最优转移轨道的轨道要素为TO(ts)={aTO,eTO,iTO,ΩTO,ωTO,θTO}。
在任务轨道转移时间T任务范围内的第一个时刻记为t1,第二个时刻记为t2,……,最后一个时刻记为ts,s表示时刻的标识号,为了方便说明ts也称为任意一个时刻。
本发明中初始轨道I到目标轨道地球静止轨道E的时间最优转移轨道的构造可以通过直接法或间接法进行设计。直接法和间接法的相关理论参见《航天器轨迹优化理论,方法及应用》,唐国金著。
步骤2:在时间最优转移轨道基础上设计时间最优状态量参考轨道;
在时间最优转移轨道TO上,以任意一个时刻ts下的轨道状态设计量作为横轴,以跟踪量AA为纵轴,构造与轨道状态量相关的时间最优状态量最优轨道TOS,如图2A、图2B、图2C所示。
所述轨道状态设计量
所述跟踪量AA包括有eTO(ts)和iTO(ts),采用集合形式表达为AA={bTO(ts),eTO(ts),iTO(ts)}。
在任务轨道转移时间T任务结束时刻,各跟踪量和状态量都将同时到达原点(图2A、图2B、图2C中的“0”点)。为了在保留时间最优状态量最优轨道TOS的长期变化规律的同时,同时消除时间最优状态量最优轨道TOS的小幅振荡,可以采用三次多项式曲线拟合生成时间最优状态量参考轨道。
时间最优状态量参考轨道记为TR,其所涉及到的轨道状态量分别是轨道伪半长轴bTR、偏心率eTR、轨道倾角iTR;状态量的横轴为应用集合表达形式为TR={bTR,eTR,iTR}。由于时间最优状态量参考轨道TR与时间最优转移轨道TO在状态时的横轴是相同的,则有
在本发明中,三次多项式曲线记为f(x)=p1x3+p2x2+p3x+p4;其中,x表示横轴轨道状态设计量,p1为三次多项式的第一个系数,p2为三次多项式的第二个系数,p3为三次多项式的第三个系数,p4为三次多项式的第四个系数。为了保证三次多项式过原点,p4为0。
步骤3:构造单圈内最优控制问题,设计单圈控制律
本发明中航天器实际转移轨道记为Q,则构建所述Q所涉及到的轨道要素分别是轨道半长轴aQ、偏心率eQ、轨道倾角iQ、升交点赤经ΩQ、近地点幅角ωQ和真近点角θQ;任意一个时刻ts下的实际转移轨道为Q(ts)={aQ,eQ,iQ,ΩQ,ωQ,θQ},实际转移轨道的状态量的横轴为T实际。
航天器实际轨道状态量记为QS,其所涉及到的轨道状态量分别是伪半长轴bQS、偏心率eQS、轨道倾角iQS;状态量的横轴为与半长轴和偏心率相关的状态设计量应用集合表达形式为QS={bQS,eQS,iQS}。所述轨道状态设计量伪半长轴偏心率eQS(ts)=eQ(ts);轨道倾角iQS(ts)=iQ(ts)。
初始轨道I到目标轨道(即地球静止轨道E)的小推力轨道转移过程通常包含上百圈的中途转移轨道,对每一圈单独设计最优控制律,该最优控制率需保证航天器的实际转移轨道QS的状态量伪半长轴bQS、偏心率eQS、轨道倾角iQS跟踪时间最优状态量参考轨道TR的同时,使得与半长轴和偏心率相关的状态设计量最快的逼近原点。
在本发明中,构造单圈内最优控制问题,从初始轨道到目标轨道将形成多圈的中途转移轨道,为了将多圈的连续推力轨道机动优化问题转化为一系列单圈的参考轨道跟踪控制问题来解决,则航天器轨道中的单圈轨道要素记为X,X={a,e,i,Ω,ω,θ},对于轨道转移过程中任意一圈轨道标识号记为k,第k圈的轨道要素通用表达记为Xk,Xk={ak,ek,ik,Ωk,ωk,θk}。第k圈的轨道状态量伪半长轴记为bk,与半长轴和偏心率相关的状态设计量记为所述伪半长轴
第k圈的起始时刻记为时刻航天器的实际轨道要素记为X0 k,所述X0 k涉及到的轨道要素分别是轨道半长轴a0 k、偏心率e0 k、轨道倾角i0 k、升交点赤经Ω0 k、近地点幅角ω0 k和真近点角θ0 k。写为集合形式有X0 k={a0 k,e0 k,i0 k,Ω0 k,ω0 k,θ0 k}。
第k圈起始时刻航天器的密切轨道周期为其中μ是地球的引力常数;第k圈终点时刻记为即通过设计第k圈的最优控制律uk,使得在时刻的性能指标Jk最小,且满足轨道动力学关系、质量导数关系、初始时刻约束关系以及跟踪时间最优状态量参考轨道的约束关系。
所述的性能指标
所述轨道动力学关系
所述质量导数关系
所述初始时刻约束关系
所述跟踪时间最优状态量参考轨道的约束关系
其中:第k圈系统的性能指标记为Jk,航天器的质量记为m,发动机推力大小记为F,发动机喷气速度记为c。分别为时间最优状态量参考轨道跟踪量AA的各元素在状态量的横轴为等于时的值。
在本发明中,将所述第k圈的最优控制律uk应用到航天器第k圈的实际轨道转移过程中,航天器实际轨道标识记为QS,实际轨道QS中的状态量伪半长轴bQS、偏心率eQS、轨道倾角iQS以及与半长轴和偏心率相关的状态设计量同时到达原点(图2A、图2B、图2C中的“0”点)时,在该原点对应的轨道要素即为目标轨道(即地球静止轨道E)的轨道要素,航天器的跟踪制导过程结束。
如果实际轨道QS中的状态量伪半长轴bQS、偏心率eQS、轨道倾角iQS以及与半长轴和偏心率相关的状态设计量未能同时到达原点,则计算第k+1圈的最优控制律uk+1,然后将所述第k+1圈的最优控制律uk+1应用到航天器第k+1圈的实际轨道转移过程中,再次判断实际轨道QS中的状态量伪半长轴bQS、偏心率eQS、轨道倾角iQS以及与半长轴和偏心率相关的状态设计量是否同时到达原点(图2A、图2B、图2C中的“0”点),若到达,则跟踪制导过程结束,若未同时到达,再次重复获取下一圈的最优控制律。
步骤4:制导策略实施及仿真验证
将第k圈获得的最优控制uk作用于航天器,在第k圈终点处实际转移轨道的状态量与时间最优状态量参考轨道状态之间的偏差,将在第k+1圈中得到修正。由此可知,通过一系列开环的最优控制过程,最终形成了一个闭环的制导过程。
以一个静止转移轨道GTO到GEO轨道转移的实例,对上述制导过程进行仿真验证。仿真程序采用MATLAB语言编写,运行环境为Windows XP操作系统,MATLAB7.0。
初始GTO轨道的经典轨道要素为
航天器质量为2000kg;电推进系统的最大推力为0.55N,发动机喷气速度为(2400×9.8)s。
若考虑轨道动力学方程为二体问题模型,利用最优控制理论采用间接法得到开环时间最优转移需要105.5天,消耗燃料212.98kg。
在本实施例中,间接法引用《航天器轨迹优化理论、方法及应用》,唐国金等著,2012年1月第一次印刷,科学出版社。第37、38页,3.1.1节内容。
利用时间最优转移轨道采用三次多项式插值,构造状态反馈的参考轨道。
在参考轨道基础上,设计单圈跟踪最优控制。同样利用最优控制理论求解单圈最优控制,单圈内轨道动力学方程考虑为二体问题模型。得到单圈控制律后,作用于航天器,此时轨道动力学模型考虑有摄动模型的精细动力学模型,同时考虑地影的影响,在地影区内推力大小为零。整个轨道转移制导过程需要110.99天,燃料消耗为224.11kg。相对于时间最优的开环解来讲,任务时间和燃料消耗都多消耗了5.2%左右。但本发明中的制导方法具有较强的鲁棒性,能够保证航天器最后进入目标轨道,最终轨道半长轴偏差小于1km,偏心率和轨道倾角偏差都在万分之一的量级;而开环时间最优解只是一个具有参考意义的理论解,无法直接满足实际工程需求。
本发明中提出的一种基于状态量反馈的小推力轨道转移制导方法,该消除传统制导方案与时间的相关性,构造了与轨道状态量相关的参考轨道;提出了一种状态量反馈形式的闭环小推力轨道机动跟踪控制方案,通过求解一系列开环单圈跟踪的最优控制问题,实现闭环轨道制导控制。该闭环反馈控制策略能够有效地消除各种摄动、地影等因素引起的轨道转移误差,保证航天器进入目标轨道,具有一定的工程操作性。
本发明中应用字母的物理意义见下表:
Claims (2)
1.一种基于状态量反馈的小推力跟踪制导方法,该方法包括有构建航天器从初始轨道转移到目标轨道的时间最优转移轨道步骤;在时间最优转移轨道基础上设计时间最优状态量参考轨道的步骤;构造单圈内最优控制问题,设计出单圈控制律的步骤;其特征在于:
构建航天器从初始轨道转移到目标轨道的时间最优转移轨道,即在小推力控制律u作用下,航天器从初始轨道I转移到目标轨道E上,并且转移时间最短的时间最优转移轨道TO;
所述时间最优转移轨道TO所涉及到的轨道要素分别是轨道半长轴aTO、偏心率eTO、轨道倾角iTO、升交点赤经ΩTO、近地点幅角ωTO和真近点角θTO;
任意一个时刻ts下的时间最优转移轨道的轨道要素为TO(ts)={aTO,eTO,iTO,ΩTO,ωTO,θTO};
在时间最优转移轨道基础上设计时间最优状态量参考轨道的步骤是指,在时间最优转移轨道TO上,以任意一个时刻ts下的设计量作为横轴,以跟踪量AA为纵轴,构造与轨道状态量相关的时间最优状态量最优轨道TOS;
所述设计量
所述跟踪量AA包括有eTO(ts)和iTO(ts),采用集合形式表达为AA={bTO(ts),eTO(ts),iTO(ts)};aE为轨道半长轴;
对时间最优状态量最优轨道TOS采用三次多项式曲线拟合生成时间最优状态量参考轨道TR,所述的TR涉及到的轨道状态量分别是轨道伪半长轴bTR、偏心率eTR、轨道倾角iTR;状态量的横轴为应用集合表达形式为TR={bTR,eTR,iTR};
构造单圈内最优控制问题,设计出单圈控制律的步骤是指,航天器轨道中的单圈轨道要素记为X,X={a,e,i,Ω,ω,θ},对于轨道转移过程中任意一圈轨道标识号记为k,第k圈的轨道要素通用表达记为Xk,Xk={ak,ek,ik,Ωk,ωk,θk};第k圈的轨道状态量伪半长轴记为bk,与半长轴和偏心率相关的状态设计量记为所述伪半长轴
第k圈的起始时刻记为时刻航天器的实际轨道要素记为X0 k,所述X0 k涉及到的轨道要素分别是轨道半长轴a0 k、偏心率e0 k、轨道倾角i0 k、升交点赤经Ω0 k、近地点幅角ω0 k和真近点角θ0 k;写为集合形式有X0 k={a0 k,e0 k,i0 k,Ω0 k,ω0 k,θ0 k};
第k圈起始时刻航天器的密切轨道周期为其中μ是地球的引力常数;第k圈终点时刻记为即通过设计第k圈的最优控制律uk,使得在时刻的性能指标Jk最小,且满足轨道动力学关系、质量导数关系、初始时刻约束关系以及跟踪时间最优状态量参考轨道的约束关系;
所述的系统性能指标
所述轨道动力学关系
所述质量导数关系
所述初始时刻约束关系
所述跟踪时间最优状态量参考轨道的约束关系
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msup>
<mi>b</mi>
<mi>k</mi>
</msup>
<mo>(</mo>
<msubsup>
<mi>t</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
<mo>)</mo>
<mo>-</mo>
<msup>
<mi>b</mi>
<mrow>
<mi>T</mi>
<mi>R</mi>
</mrow>
</msup>
<mo>&lsqb;</mo>
<msubsup>
<mi>x</mi>
<mi>G</mi>
<mi>k</mi>
</msubsup>
<mo>(</mo>
<msubsup>
<mi>t</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
<mo>)</mo>
<mo>&rsqb;</mo>
<mo>=</mo>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>e</mi>
<mi>k</mi>
</msup>
<mo>(</mo>
<msubsup>
<mi>t</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
<mo>)</mo>
<mo>-</mo>
<msup>
<mi>e</mi>
<mrow>
<mi>T</mi>
<mi>R</mi>
</mrow>
</msup>
<mo>&lsqb;</mo>
<msubsup>
<mi>x</mi>
<mi>G</mi>
<mi>k</mi>
</msubsup>
<mo>(</mo>
<msubsup>
<mi>t</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
<mo>)</mo>
<mo>&rsqb;</mo>
<mo>=</mo>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>i</mi>
<mi>k</mi>
</msup>
<mo>(</mo>
<msubsup>
<mi>t</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
<mo>)</mo>
<mo>-</mo>
<msup>
<mi>i</mi>
<mrow>
<mi>T</mi>
<mi>R</mi>
</mrow>
</msup>
<mo>&lsqb;</mo>
<msubsup>
<mi>x</mi>
<mi>G</mi>
<mi>k</mi>
</msubsup>
<mo>(</mo>
<msubsup>
<mi>t</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
<mo>)</mo>
<mo>&rsqb;</mo>
<mo>=</mo>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
其中:第k圈的系统性能指标记为Jk,航天器的质量记为m,发动机推力大小记为F,发动机喷气速度记为c;分别为时间最优状态量参考轨道跟踪量AA的各元素在状态量的横轴为等于时的值。
2.根据权利要求1所述的基于状态量反馈的小推力跟踪制导方法,其特征在于:应用所述的基于状态量反馈的小推力跟踪制导方法进行轨道转移的执行步骤为:
(A)依据时间最优的参考轨道TO,获取状态量参考轨道TR;
(B)依据状态量参考轨道TR设计单圈的最优控制律uk;k表示圈数;
(C)航天飞行器按照单圈的最优控制律u进行轨道转移;
(D)通过当前时刻的实际轨道状态信息QS和状态量参考轨道TR计算下一圈的最优控制律uk+1;
(E)航天飞行器按照uk+1进行轨道转移;
(F)重复(D)和(E)步骤,直至设计量为零,完成飞行任务。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410073905.1A CN103853047B (zh) | 2014-02-28 | 2014-02-28 | 一种基于状态量反馈的小推力跟踪制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410073905.1A CN103853047B (zh) | 2014-02-28 | 2014-02-28 | 一种基于状态量反馈的小推力跟踪制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103853047A CN103853047A (zh) | 2014-06-11 |
CN103853047B true CN103853047B (zh) | 2018-01-09 |
Family
ID=50860846
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410073905.1A Expired - Fee Related CN103853047B (zh) | 2014-02-28 | 2014-02-28 | 一种基于状态量反馈的小推力跟踪制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103853047B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105138011B (zh) * | 2015-08-31 | 2017-10-03 | 哈尔滨工业大学 | 一种航天器在轨服务观测空间目标局部范围的时间和燃料脉冲最优遍历方法 |
CN108303874B (zh) * | 2017-12-28 | 2020-06-23 | 北京航空航天大学 | 一种针对绳系空间拖船系统摆振的小推力切换控制方法 |
CN109164829B (zh) * | 2018-10-23 | 2021-08-27 | 哈尔滨工业大学(深圳) | 一种基于力反馈装置和vr感知的飞行机械臂系统及控制方法 |
CN113955153B (zh) * | 2021-12-21 | 2022-04-08 | 南京航空航天大学 | 一种燃料最优的连续小推力轨道转移方法 |
CN114384806B (zh) * | 2022-01-12 | 2023-05-23 | 北京理工大学 | 多摄动地影约束下电推进航天器多圈变轨的分段优化方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6341249B1 (en) * | 1999-02-11 | 2002-01-22 | Guang Qian Xing | Autonomous unified on-board orbit and attitude control system for satellites |
CN103226631A (zh) * | 2013-03-29 | 2013-07-31 | 南京航空航天大学 | 一种小推力转移轨道快速设计与优化方法 |
-
2014
- 2014-02-28 CN CN201410073905.1A patent/CN103853047B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6341249B1 (en) * | 1999-02-11 | 2002-01-22 | Guang Qian Xing | Autonomous unified on-board orbit and attitude control system for satellites |
CN103226631A (zh) * | 2013-03-29 | 2013-07-31 | 南京航空航天大学 | 一种小推力转移轨道快速设计与优化方法 |
Non-Patent Citations (5)
Title |
---|
Low-Thrust Guidance Scheme for the Spacecraft using B-plane Targeting and Lyapunov Feedback control;Dong-Hyun Cho等;《International Conference on Control, Automation and Systems》;20121021(第12期);第1969-1973页 * |
基于最优状态反馈的小推力转移轨道制导策略;尚海滨等;《吉林大学学报(工学版)》;20070731;第37卷(第4期);第949-954页 * |
小推力轨道保持方法;吕秋杰等;《上海航天》;20101231(第4期);第23-28,42页 * |
无需初值猜测的间接法小推力轨道设计;林厚源等;《天文学报》;20120531;第53卷(第3期);第197-204页 * |
有限推力下时间最优轨道转移;梁新刚等;《航天控制》;20070228;第25卷(第1期);第46-51页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103853047A (zh) | 2014-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hu et al. | Robust fault-tolerant tracking control for spacecraft proximity operations using time-varying sliding mode | |
CN103853047B (zh) | 一种基于状态量反馈的小推力跟踪制导方法 | |
CN104527994B (zh) | 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法 | |
Hu et al. | Disturbance observer based finite-time attitude control for rigid spacecraft under input saturation | |
Lu et al. | Optimal aerocapture guidance | |
Morgan et al. | Model predictive control of swarms of spacecraft using sequential convex programming | |
Naasz | Classical element feedback control for spacecraft orbital maneuvers | |
CN106814746A (zh) | 一种航天器姿轨一体化反步跟踪控制方法 | |
CN104309822B (zh) | 一种基于参数优化的航天器单脉冲水滴形绕飞轨迹悬停控制方法 | |
CN102819266B (zh) | 一种拟周期j2不变相对轨道编队飞行控制方法 | |
CN105930305B (zh) | 一种三脉冲交会接近制导方法 | |
CN106114910A (zh) | 一种航天器飞行轨道滚动时域控制方法 | |
Wu et al. | Null motion strategy for spacecraft large angle agile maneuvering using hybrid actuators | |
Zhao et al. | Maneuver control of the hypersonic gliding vehicle with a scissored pair of control moment gyros | |
Morgan et al. | Decentralized model predictive control of swarms of spacecraft using sequential convex programming | |
Wang et al. | Control analysis of an underactuated spacecraft under disturbance | |
Brown et al. | Energetics of control moment gyroscopes as joint actuators | |
Mu et al. | Trajectory planning for terminal area energy management phase of reusable launch vehicles | |
Bucchioni et al. | Optimal time-fixed impulsive non-Keplerian orbit to orbit transfer algorithm based on primer vector theory | |
CN104793613A (zh) | 一种航天器在日地系统不稳定平动点轨道间交会控制方法 | |
Lei et al. | A steering method with multiobjective optimizing for nonredundant single-gimbal control moment gyro systems | |
CN108891625A (zh) | 固体微推进器阵列与磁力矩器联合控制方法 | |
Wang et al. | Controllability of an underactuated spacecraft with one thruster under disturbance | |
Xu et al. | Coordinated control method of space-tethered robot system for tracking optimal trajectory | |
Sun et al. | Piecewise attitude tracking control of a gravity gradient microsatellite for coplanar orbital transfer |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180109 Termination date: 20190228 |