CN106628257B - 地球摄动引力场中近地航天器相对运动轨道的保持方法 - Google Patents
地球摄动引力场中近地航天器相对运动轨道的保持方法 Download PDFInfo
- Publication number
- CN106628257B CN106628257B CN201610860433.3A CN201610860433A CN106628257B CN 106628257 B CN106628257 B CN 106628257B CN 201610860433 A CN201610860433 A CN 201610860433A CN 106628257 B CN106628257 B CN 106628257B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- msub
- mtr
- mover
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 230000000737 periodic effect Effects 0.000 claims abstract description 43
- 238000011217 control strategy Methods 0.000 claims abstract description 17
- 239000013598 vector Substances 0.000 claims description 30
- 239000011159 matrix material Substances 0.000 claims description 21
- 239000006185 dispersion Substances 0.000 claims description 18
- 230000005484 gravity Effects 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000008707 rearrangement Effects 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 230000001174 ascending effect Effects 0.000 claims 1
- 230000015572 biosynthetic process Effects 0.000 abstract description 2
- 230000000694 effects Effects 0.000 abstract description 2
- 238000011160 research Methods 0.000 abstract description 2
- 230000006835 compression Effects 0.000 abstract 1
- 238000007906 compression Methods 0.000 abstract 1
- 238000011438 discrete method Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 239000000446 fuel Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000012456 homogeneous solution Substances 0.000 description 1
- 238000011423 initialization method Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
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
- B64G1/242—Orbits and trajectories
-
- 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)
- External Artificial Organs (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种地球摄动引力场中近地航天器相对运动轨道的保持方法,通过对地球扁率摄动影响下近地轨道卫星的相对运动进行建模,并且提出一种基于时间离散法的周期解估计方案,能够简单有效地对周期性相对运动的初始条件进行较为精确的估计,以及定量地研究近地飞行器的相对运动在受地心引力摄动的作用下动力学特性的变化。为飞行器的编队飞行任务及其他任务中飞行器的绕飞阶段的轨道设计提供精确可靠的数学模型和数据参考。并且将一种闭合投影轨道引入相对运动轨道的保持方案,通过与DLQR控制方法结合进而设计了一种高效率、低燃耗的控制策略。
Description
技术领域
本发明属于航空航天技术领域,涉及一种近地航天器相对运动轨道的保持方法,特别涉及地球摄动引力场中近地航天器相对运动轨道的保持方法。
背景技术
近地轨道航天器的相对运动问题在过去的几十年中受到了大量相关研究人员的关注。许多二体问题的动力学和运动学模型被相继提出,最为著名的相对运动数学模型是C-W方程。C-W方程提供了一种求解近圆轨道航天器相对运动周期轨道的方法。然而,C-W方程中采用的一些假设条件(假设圆轨道、忽略地球扁率等)会导致模型的误差在某些情况下严重偏离实际。由于周期性相对轨道在相对轨道保持,以及航天器编队飞行方面具有重要意义,迫切需要更加精确的相对运动模型以及求解相对运动周期轨道的方法。迄今,已经有过许多对C-W方程进行一般化的尝试。通过研究得到了相对运动线性方程组的齐次解,得到了T-H方程的周期轨道初始条件。该结果能用于对参考轨道为椭圆轨道的大型航天器集群飞行进行初始化。然而,当同时考虑重力场的非线性、地球J2摄动、和大偏心率参考轨道时,初始条件就无法再满足要求了。当考虑非线性项的存在时,非线性动力学模型都没有解析解,因为模型方程不具备线性系统的可加性和齐次性。总的来说,当前还未能针对近地轨道航天器的相对运动提出一种精确、可靠、适应性强的周期轨道初始化方法,也未针对航天器相对运动轨道的保持提供一种普遍有效、节省燃耗的控制策略。
发明内容
本发明的目的在于提供一种任意形状翻滚卫星的非接触式惯量系数辨识方法,以克服上述现有技术存在的缺陷,本发明提出了基于时间离散法的相对运动轨道初值快速解算方法和轨道保持控制方法,计算效率高,并且节省控制燃料,因此具有明显的优势。
为达到上述目的,本发明采用如下技术方案:
地球摄动引力场中近地航天器相对运动轨道的保持方法,包括如下步骤:
1)建立地球引力摄动下主从星的精确相对运动模型;
2)使用时间离散法对主从星相对运动模型进行求解;
3)获取闭合投影轨道;
4)设计闭环控制策略。
进一步地,步骤1)中建立地球引力摄动下主从星的精确相对运动模型具体为:
同时采用地心惯性坐标系和星体坐标系建立从航天器相对于主航天器的椭圆轨道相对运动模型:
其中,xr,yr,zr为从星在主星星体坐标系下的位置坐标,ωx,ωz,αx,αz分别为主星星体坐标系绕自身x轴、z轴的旋转角速度及相应的旋转角加速度,η由如下关系式给出:
r,i,和θ均为描述主星轨道的周期性时变参数并由如下主星轨道方程给出:
其中,s表示三角函数sin,c表示三角函数cos,r为主星相对于地心的距离,vr为r相对于时间的导数,h为主星的角动量,θ为主星与升交点之间的夹角,i为轨道平面倾角,Ω为升交点赤经,kJ2为J2相关项;
Fx,Fy,和Fz是从星的控制力在主星星体坐标系下的分量,ηr和是相对位置坐标xr,yr,zr的函数并由如下关系给出:
其中,
rrZ=(r+xr)sisθ+yrsicθ+zrci,
在上式中,rr为从星相对于地心的距离,rrZ为rr在地心惯性系Z轴上的投影,J2为地心引力场的摄动项,μ为引力常数,Re为地球半径;
然后,通过降阶增维的方式,将主从星相对运动方程中的二阶导数消去,从而将其转化为如下的一阶状态空间形式:
其中,X=(xr yr zr vx vy vz)T为状态空间向量,G为映射函数。
进一步地,步骤2)中使用时间离散法对主从星相对运动模型进行求解具体为:
将状态空间向量X假设为傅里叶级数的周期解形式,以函数f(t)为例,将其表示为:
其中a0,an,bn为谐波系数,ω为周期运动的频率,n为谐波的阶数,t表示时间;
根据傅里叶级数本身的特性,这种形式的周期解存在如下所示的变换关系:
其中f(tj)和分别为时间tj处函数f(t)及其一阶时间导数的值,矩阵EAE-1为相应的变换矩阵,ωf为假设的周期运动频率,假设相对运动的周期为时间T,在该周期内选取K个时间点,即得到如下关于离散点的方程组:
其中,t1至tK为在时间周期内选取的配点,根据傅里叶级数形式的周期解所固有的变换关系,将上式的左端项转化为:
其中,ωxr,ωyr,和ωzr分别为x,y,z轴的运动分量的频率,ExAEx -1,EyAEy -1和EzAEz -1分别为x,y,z轴的运动分量的变换矩阵,Xr,Yr,Zr和Vx,Vy,Vz分别为对应位置和速度分量在时间配点处的值;
由此,得到该模型对应的时间离散法代数方程组,简写如下
其中,Q=(Xr,Yr,Zr,Vx,Vy,Vz)T,是对(G(X(t1)),G(X(t2)),...G(X(tK)))T重新排列后得到的矢量,为对E重新排列后得到的矩阵,通过对上述方程组进行求解即得到周期解的离散点信息,利用该离散点信息即能够得到周期解的估计值。
进一步地,步骤3)中获取闭合投影轨道具体为:
假设闭合投影轨道能够展开为傅里叶级数的形式,通过在慢漂移轨道的某一段上选择均匀分布的M个离散点,得到余量函数其中,矢量为闭合投影轨道的谐波估计系数,为相应的谐波矩阵,Qc为慢漂移轨道在配点处的位置和速度信息,通过最小二乘法,寻找使RTR取最小值的闭合投影轨道的频率ωc *和离散信息即求解下式:
其中,R为余量函数,ωc *为最优的闭合投影轨道的频率,为该最优闭合投影轨道在配点处的位置及速度信息,将余量函数代入上式,得到:
该非线性代数方程的解ωc *和即用于构造闭合投影轨道,
上式中,矩阵和T如下所示:
其中,A的具体形式为:
t为配点处的时间t1至tK构成的对角阵。
进一步地,步骤4)中设计闭环控制策略具体为:
对相对运动模型进行线性化处理之后得到下式
其中Xd=X-Xp是真实轨迹X相对投影轨道的漂移量,Xp是投影轨道的状态矢量,矩阵为控制矩阵,如下所示;
令控制矢量u为固定时间间隔的脉冲推力,即:
其中tk为脉冲推力开始作用的瞬时时间,为脉冲间隔Ts的整数倍:tk=kTs,d为脉冲推力作用的持续时间,uk为离散控制器给出的控制信号,由此将相对运动模型离散化为:
其中,dτ为对时间τ的微分;
使用DLQR控制策略,使如下的二次型性能指标最小化
其中,Xdk为在时间域内离散化的状态矢量;
通过Ricatti方程得到最优控制律
uk=-K(X(tk)-Xp(tk))
其中,X(tk)为真实状态矢量在tk处的值,Xp(tk)为闭合投影轨道的状态矢量在tk处的值,K为增益矩阵。
与现有技术相比,本发明具有以下有益的技术效果:
本发明方法基于时间离散法,通过在时间域内离散获得航天器相对运动的代数方程,避免了一般的非线性系统分析方法中大量的符号运算,同时简化了系统方程,使得求解效率得到较大提高。本发明所提出的时间离散法求解方案扩展性较强,能够应用到考虑其他摄动因素(如太阳光压和大气阻力)的更加精确的相对运动模型的求解中。另外,提出了一种基于闭合投影轨道的闭环控制策略,由于闭合投影轨道与真实的相对运动周期轨道近似程度较高,这使得该控制策略能够有效地降低轨道保持的燃料消耗,同时闭合投影轨道的求解过程简单,不会过多占用星载计算机的计算资源。
附图说明
图1为本发明方法的流程图;
图2为航天器相对运动模型示意图;
图3为时间离散法求解方案得到的相对运动轨道;
图4为为无控条件下航天器的慢漂移轨迹和闭合投影轨道;
图5为DLQR控制策略下的相对运动轨迹。
具体实施方式
下面对本发明作进一步详细描述:
采用时间离散法对主从星相对运动模型进行求解,将相对运动状态矢量的三个分量假设为傅里叶级数的周期解形式。然后通过周期内的时间离散,得到离散后的系统的代数方程组,通过选取合适的不动点迭代计算方法及迭代初值,获得周期解的估计解,即为相对周期运动的精确初值。使用时间离散法对系统模型进行求解后,能够得到精度较高的估计解,但是由于求解过程中存在计算误差,由该初始条件得到的相对运动轨迹并不完全闭合,而是一条漂移速率较慢的近似周期轨道。我们采用最小二乘原理将慢漂移轨道投影到一个周期轨道上,得到闭合投影轨道。基于此闭合投影轨道,通过与DLQR控制方法结合进而设计了一种高效率、低燃耗的控制策略。
具体步骤包括:
1)基于拉格朗日方程,建立地球引力场在J2摄动下航天器的精确相对运动模型:
同时采用地心惯性坐标系(ECI)和星体坐标系(LVLH)建立从航天器相对于主航天器的椭圆轨道相对运动模型
其中,xr,yr,zr为从星在主星星体坐标系下的位置坐标,ωx、ωz,αx,αz分别为主星星体坐标系绕自身x轴、z轴的旋转角速度及相应的旋转角加速度,为η由如下关系式给出
r,i,和θ均为描述主星轨道的周期性时变参数并由如下主星轨道方程给出:
其中,s表示三角函数sin,c表示三角函数cos,r为主星相对于地心的距离,vr为r相对于时间的导数,h为主星的角动量,θ为主星与升交点之间的夹角,i为轨道平面倾角,Ω为升交点赤经,kJ2为J2相关项;
Fx,Fy,和Fz是从星的控制力在主星星体坐标系下的分量,ηr和是相对位置坐标xr,yr,zr的函数并由如下关系给出:
其中,
rrZ=(r+xr)sisθ+yrsicθ+zrci,
在上式中,rr为从星相对于地心的距离,rrZ为rr在地心惯性系Z轴上的投影,J2为地心引力场的摄动项,μ为引力常数,Re为地球半径;
然后,通过降阶增维的方式,将主从星相对运动方程中的二阶导数消去,从而将其转化为如下的一阶状态空间形式:
其中,X=(xr yr zr vx vy vz)T为状态空间向量,G为映射函数。
2)建立时间离散法求解方案对航天器相对运动模型进行求解:
将状态空间向量X假设为傅里叶级数的周期解形式,以函数f(t)为例,将其表示为:
其中a0,an,bn为谐波系数,ω为周期运动的频率,n为谐波的阶数,t表示时间;
根据傅里叶级数本身的特性,这种形式的周期解存在如下所示的变换关系
其中f(tj)和分别为时间tj处函数f(t)及其一阶时间导数的值,矩阵EAE-1为相应的变换矩阵,ωf为假设的周期运动频率。假设相对运动的周期为时间T,在该周期内选取K个时间点,可以得到如下关于离散点的方程组
其中,t1至tK为在时间周期内选取的配点,根据傅里叶级数形式的周期解所固有的变换关系,上式的左端项可以转化为
其中,ωxr,ωyr,和ωzr分别为x,y,z轴的运动分量的频率。ExAEx -1,EyAEy -1和EzAEz -1分别为x,y,z轴的运动分量的变换矩阵,Xr,Yr,Zr和Vx,Vy,Vz分别为对应位置和速度分量在时间配点处的值;
由此,能够得到该模型对应的时间离散法代数方程组,简写如下
其中,Q=(Xr,Yr,Zr,Vx,Vy,Vz)T。是对(G(X(t1)),G(X(t2)),...G(X(tK)))T重新排列后得到的矢量,为对E重新排列后得到的矩阵。上式为一组非线性代数方程,通过选取合适的不动点迭代计算方法及迭代初值,就能进行求解,所得到的解Q为周期解的离散点信息,利用该信息能够获得周期解的估计解,也能够对相对周期运动进行初始化。
3)获取闭合投影轨道
使用时间离散法对系统模型进行求解后,能够得到精度较高的估计解,同时也能获得周期轨道的初始化条件,但是由于求解过程中存在计算误差,由该初始条件得到的相对运动轨迹并不完全闭合,而是一条漂移速率较慢的近似周期轨道。为了获得相对运动控制所需的参考轨迹,我们将其投影到一个周期性的闭合轨道上,并以该轨道为基础,设计闭环控制策略。
仍然假设闭合投影轨道能够展开为傅里叶级数的形式,通过在慢漂移轨道的某一段上选择均匀分布的M个离散点,得到余量函数其中,矢量为闭合投影轨道的谐波估计系数,为相应的谐波矩阵,Qc为慢漂移轨道在配点处的位置和速度信息,通过最小二乘法,寻找使RTR取最小值的闭合投影轨道的频率ωc *和离散信息也就是求解下式
其中,R为余量函数,ωc *为最优的闭合投影轨道的频率,为该最优闭合投影轨道在配点处的位置及速度信息,将余量函数代入上式,可以得到
该非线性代数方程的解ωc *和即用于构造闭合投影轨道。
其中,矩阵和T如下所示:
其中,A的具体形式为:
t为配点处的时间t1至tK构成的对角阵。
4)设计闭环控制策略:
对相对运动模型进行线性化处理之后得到下式
其中Xd=X-Xp是真实轨迹X相对投影轨道的漂移量,Xp是投影轨道的状态矢量,矩阵为控制矩阵,如下所示;
令控制矢量u为固定时间间隔的脉冲推力,即:
其中tk为脉冲推力开始作用的瞬时时间,为脉冲间隔Ts的整数倍:tk=kTs,d为脉冲推力作用的持续时间,uk为离散控制器给出的控制信号,由此将相对运动模型离散化为:
其中,dτ为对时间τ的微分。
使用DLQR控制策略,使如下的二次型性能指标最小化
其中,Xdk为在时间域内离散化的状态矢量;
通过Ricatti方程得到最优控制律
uk=-K(X(tk)-Xp(tk))。
其中,X(tk)为真实状态矢量在tk处的值,Xp(tk)为闭合投影轨道的状态矢量在tk处的值,K为增益矩阵。
为了更好地说明本发明的目的和优点,下面结合附图和实例对本发明内容做进一步说明:
本发明针对近地轨道近距离交会任务,提出了一种近地轨道航天器相对运动轨道的保持方法,该方法的操作流程如图1所示,其具体实施方式包括如下步骤:
第一步:建立航天器的相对运动模型。该模型中采用了两组笛卡尔坐标系,如图2所示。其中地心惯性坐标系(ECI)由一组单位向量X,Y,Z表示,体坐标系(LVLH)固联在参考航天器S0上。建模过程中考虑地球引力场的J2摄动项,建立了一阶标准化动力学方程其中X=(xr yr zr vx vy vz)T,也就是相对运动的状态矢量。
第二步:相对运动模型的时间离散法求解方案。取周期轨道的近似解为傅里叶级数的谐波形式,其中傅里叶级数的系数为未知量。将近似解带入动力学方程然后在的一个周期上进行离散,得到代数方程组。使用牛顿迭代法求解代数方程组,得到近似解。这样,得到的估计解X(t)与真实的周期解足够接近,是一个慢漂移的相对运动轨道,如图3。
第三步:获取相对运动的精确闭合投影轨道。
上一步中,使用时间离散法对系统模型进行求解后,能够得到精度较高的估计解,同时也能获得周期轨道的初始化条件,但是由于求解过程中存在计算误差,由该初始条件得到的相对运动轨迹并不完全闭合,而是一条漂移速率较慢的近似周期轨道。为了获得相对运动控制所需的参考轨迹,我们将其投影到一个周期性的闭合轨道上,并以该轨道为基础,设计闭环控制策略。仍然假设闭合投影轨道能够展开为傅里叶级数的形式,通过在慢漂移轨道的某一段上选择均匀分布有限个离散点,带入控制方程得到余量函数通过最小二乘法,寻找使RTR取最小值的闭合投影轨道的高精度闭合解,如图4。
第四步:设计闭环控制策略。基于获取的相对运动的精确闭合投影轨道,并采用离散化线性二次型最优控制器进行周期轨道控制,在消耗燃料很低的情况下,得到了精确的相对运动周期轨道,如图5所示。
Claims (4)
1.地球摄动引力场中近地航天器相对运动轨道的保持方法,其特征在于,包括如下步骤:
1)建立地球引力摄动下主从星的精确相对运动模型;具体为:
同时采用地心惯性坐标系和星体坐标系建立从航天器相对于主航天器的椭圆轨道相对运动模型:
其中,xr,yr,zr为从星在主星星体坐标系下的位置坐标,ωx,ωz,αx,αz分别为主星星体坐标系绕自身x轴、z轴的旋转角速度及相应的旋转角加速度,η由如下关系式给出:
r,i,和θ均为描述主星轨道的周期性时变参数并由如下主星轨道方程给出:
<mrow>
<mover>
<mi>r</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>v</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mover>
<mi>v</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>r</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mi>&mu;</mi>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<mfrac>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
<msup>
<mi>r</mi>
<mn>3</mn>
</msup>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>k</mi>
<mrow>
<mi>J</mi>
<mn>2</mn>
</mrow>
</msub>
<msup>
<mi>r</mi>
<mn>4</mn>
</msup>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>3</mn>
<msubsup>
<mi>s</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>s</mi>
<mi>&theta;</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<mover>
<mi>h</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>J</mi>
<mn>2</mn>
</mrow>
</msub>
<msubsup>
<mi>s</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>s</mi>
<mi>&theta;</mi>
<mn>2</mn>
</msubsup>
</mrow>
<msup>
<mi>r</mi>
<mn>3</mn>
</msup>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<mover>
<mi>&theta;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mi>h</mi>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mrow>
<mi>J</mi>
<mn>2</mn>
</mrow>
</msub>
<msubsup>
<mi>c</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>s</mi>
<mi>&theta;</mi>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<msup>
<mi>hr</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<mover>
<mi>i</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>J</mi>
<mn>2</mn>
</mrow>
</msub>
<msubsup>
<mi>s</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>s</mi>
<mi>&theta;</mi>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>hr</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<mover>
<mi>&Omega;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mrow>
<mi>J</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>s</mi>
<mi>&theta;</mi>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<msup>
<mi>hr</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,s表示三角函数sin,c表示三角函数cos,r为主星相对于地心的距离,vr为r相对于时间的导数,h为主星的角动量,θ为主星与升交点之间的夹角,i为轨道平面倾角,Ω为升交点赤经,kJ2为J2相关项;
Fx,Fy,和Fz是从星的控制力在主星星体坐标系下的分量,ηr和是相对位置坐标xr,yr,zr的函数并由如下关系给出:
其中,
<mrow>
<msub>
<mi>r</mi>
<mi>r</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>+</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>y</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>z</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>,</mo>
</mrow>
rrZ=(r+xr)sisθ+yrsicθ+zrci,
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>J</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msub>
<mi>J</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>&mu;R</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>,</mo>
</mrow>
在上式中,rr为从星相对于地心的距离,rrZ为rr在地心惯性系Z轴上的投影,J2为地心引力场的摄动项,μ为引力常数,Re为地球半径;
然后,通过降阶增维的方式,将主从星相对运动方程中的二阶导数消去,从而将其转化为如下的一阶状态空间形式:
<mrow>
<mover>
<mi>X</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,X=(xr yr zr vx vy vz)T为状态空间向量,G为映射函数;
2)使用时间离散法对主从星相对运动模型进行求解;
3)获取闭合投影轨道;
4)设计闭环控制策略。
2.根据权利要求1所述的地球摄动引力场中近地航天器相对运动轨道的保持方法,其特征在于,步骤2)中使用时间离散法对主从星相对运动模型进行求解具体为:
将状态空间向量X假设为傅里叶级数的周期解形式,以函数f(t)为例,将其表示为:
其中a0,an,bn为谐波系数,ω为周期运动的频率,n为谐波的阶数,t表示时间;
根据傅里叶级数本身的特性,这种形式的周期解存在如下所示的变换关系:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>f</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>f</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>f</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>K</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msub>
<mi>&omega;</mi>
<mi>f</mi>
</msub>
<msup>
<mi>EAE</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>K</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中f(tj)和分别为时间tj处函数f(t)及其一阶时间导数的值,矩阵EAE-1为相应的变换矩阵,ωf为假设的周期运动频率,假设相对运动的周期为时间T,在该周期内选取K个时间点,即得到如下关于离散点的方程组:
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mover>
<mi>X</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mover>
<mi>X</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mover>
<mi>X</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>K</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>G</mi>
<mo>(</mo>
<mrow>
<mi>X</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>,</mo>
<mi>G</mi>
<mo>(</mo>
<mrow>
<mi>X</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>,</mo>
<mn>...</mn>
<mi>G</mi>
<mo>(</mo>
<mrow>
<mi>X</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>K</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
其中,t1至tK为在时间周期内选取的配点,根据傅里叶级数形式的周期解所固有的变换关系,将上式的左端项转化为:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>X</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>r</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>Y</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>r</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>Z</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>r</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>V</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>x</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>V</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>y</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>V</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>z</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>x</mi>
<mi>r</mi>
</mrow>
</msub>
<msub>
<mi>E</mi>
<mi>x</mi>
</msub>
<msup>
<msub>
<mi>AE</mi>
<mi>x</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>y</mi>
<mi>r</mi>
</mrow>
</msub>
<msub>
<mi>E</mi>
<mi>y</mi>
</msub>
<msup>
<msub>
<mi>AE</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>z</mi>
<mi>r</mi>
</mrow>
</msub>
<msub>
<mi>E</mi>
<mi>z</mi>
</msub>
<msup>
<msub>
<mi>AE</mi>
<mi>z</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>x</mi>
<mi>r</mi>
</mrow>
</msub>
<msub>
<mi>E</mi>
<mi>x</mi>
</msub>
<msup>
<msub>
<mi>AE</mi>
<mi>x</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>y</mi>
<mi>r</mi>
</mrow>
</msub>
<msub>
<mi>E</mi>
<mi>y</mi>
</msub>
<msup>
<msub>
<mi>AE</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>z</mi>
<mi>r</mi>
</mrow>
</msub>
<msub>
<mi>E</mi>
<mi>z</mi>
</msub>
<msup>
<msub>
<mi>AE</mi>
<mi>z</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>X</mi>
<mi>r</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Y</mi>
<mi>r</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mi>r</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>x</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>y</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>z</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,ωxr,ωyr,和ωzr分别为x,y,z轴的运动分量的频率,ExAEx -1,EyAEy -1和EzAEz -1分别为x,y,z轴的运动分量的变换矩阵,Xr,Yr,Zr和Vx,Vy,Vz分别为对应位置和速度分量在时间配点处的值;
由此,得到该模型对应的时间离散法代数方程组,简写如下
<mrow>
<mover>
<mi>G</mi>
<mo>~</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mi>E</mi>
<mo>~</mo>
</mover>
<mi>Q</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
其中,Q=(Xr,Yr,Zr,Vx,Vy,Vz)T,是对(G(X(t1)),G(X(t2)),...G(X(tK)))T重新排列后得到的矢量,为对E重新排列后得到的矩阵,通过对上述方程组进行求解即得到周期解的离散点信息,利用该离散点信息即能够得到周期解的估计值。
3.根据权利要求2所述的地球摄动引力场中近地航天器相对运动轨道的保持方法,其特征在于,步骤3)中获取闭合投影轨道具体为:
假设闭合投影轨道能够展开为傅里叶级数的形式,通过在慢漂移轨道的某一段上选择均匀分布的M个离散点,得到余量函数其中,矢量为闭合投影轨道的谐波估计系数,为相应的谐波矩阵,Qc为慢漂移轨道在配点处的位置和速度信息,通过最小二乘法,寻找使RTR取最小值的闭合投影轨道的频率ωc *和离散信息即求解下式:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>R</mi>
<mi>T</mi>
</msup>
<mi>R</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<msub>
<mi>&omega;</mi>
<mi>c</mi>
</msub>
<mo>*</mo>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>R</mi>
<mi>T</mi>
</msup>
<mi>R</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mover>
<mi>Q</mi>
<mo>^</mo>
</mover>
<mo>*</mo>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
</mrow>
其中,R为余量函数,ωc *为最优的闭合投影轨道的频率,为该最优闭合投影轨道在配点处的位置及速度信息,将余量函数代入上式,得到:
该非线性代数方程的解ωc *和即用于构造闭合投影轨道,
上式中,矩阵和T如下所示:
其中,A的具体形式为:
t为配点处的时间t1至tK构成的对角阵。
4.根据权利要求1所述的地球摄动引力场中近地航天器相对运动轨道的保持方法,其特征在于,步骤4)中设计闭环控制策略具体为:
对相对运动模型进行线性化处理之后得到下式
<mrow>
<msub>
<mover>
<mi>X</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mo>=</mo>
<mover>
<mi>A</mi>
<mo>&OverBar;</mo>
</mover>
<msub>
<mi>X</mi>
<mi>d</mi>
</msub>
<mo>+</mo>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mi>u</mi>
</mrow>
其中Xd=X-Xp是真实轨迹X相对投影轨道的漂移量,Xp是投影轨道的状态矢量,矩阵为控制矩阵,如下所示;
<mrow>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
令控制矢量u为固定时间间隔的脉冲推力,即:
<mrow>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<msub>
<mi>u</mi>
<mi>k</mi>
</msub>
<mi>d</mi>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo><</mo>
<mi>t</mi>
<mo><</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<mi>d</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<mi>d</mi>
<mo><</mo>
<mi>t</mi>
<mo><</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中tk为脉冲推力开始作用的瞬时时间,为脉冲间隔Ts的整数倍:tk=kTs,d为脉冲推力作用的持续时间,uk为离散控制器给出的控制信号,由此将相对运动模型离散化为:
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<msub>
<mi>X</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<mover>
<mi>B</mi>
<mo>~</mo>
</mover>
<msub>
<mi>u</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
</mrow>
其中,dτ为对时间τ的微分;
使用DLQR控制策略,使如下的二次型性能指标最小化
<mrow>
<mi>J</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<mo>&lsqb;</mo>
<msup>
<msub>
<mi>X</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
</msub>
<mi>T</mi>
</msup>
<msub>
<mi>QX</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<msup>
<msub>
<mi>u</mi>
<mi>k</mi>
</msub>
<mi>T</mi>
</msup>
<msub>
<mi>Ru</mi>
<mi>k</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
其中,Xdk为在时间域内离散化的状态矢量;
通过Ricatti方程得到最优控制律
uk=-K(X(tk)-Xp(tk))
其中,X(tk)为真实状态矢量在tk处的值,Xp(tk)为闭合投影轨道的状态矢量在tk处的值,K为增益矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610860433.3A CN106628257B (zh) | 2016-09-28 | 2016-09-28 | 地球摄动引力场中近地航天器相对运动轨道的保持方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610860433.3A CN106628257B (zh) | 2016-09-28 | 2016-09-28 | 地球摄动引力场中近地航天器相对运动轨道的保持方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106628257A CN106628257A (zh) | 2017-05-10 |
CN106628257B true CN106628257B (zh) | 2018-01-19 |
Family
ID=58853568
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610860433.3A Active CN106628257B (zh) | 2016-09-28 | 2016-09-28 | 地球摄动引力场中近地航天器相对运动轨道的保持方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106628257B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490966B (zh) * | 2018-01-31 | 2021-02-05 | 中国人民解放军国防科技大学 | 基于微分代数的静止轨道摄动相对轨迹高阶制导方法 |
CN109145391A (zh) * | 2018-07-26 | 2019-01-04 | 西北工业大学 | 一种面向任务的航天器连续推力轨道设计方法 |
CN110321598B (zh) * | 2019-06-10 | 2023-04-18 | 南京航空航天大学 | 一种j2项摄动条件下航天器相对运动解析求解方法 |
CN110209053A (zh) * | 2019-06-10 | 2019-09-06 | 西北工业大学深圳研究院 | 集群卫星系统的控制器设计方法及控制方法 |
CN111591469B (zh) * | 2020-03-03 | 2021-12-07 | 航天科工空间工程发展有限公司 | 一种低轨星座系统相位保持方法、系统、设备和存储介质 |
CN111731513B (zh) * | 2020-06-15 | 2022-03-04 | 航天东方红卫星有限公司 | 一种基于单脉冲轨控的高精度引力场中回归轨道维持方法 |
CN112046794B (zh) * | 2020-07-16 | 2022-02-25 | 中国人民解放军军事科学院国防科技创新研究院 | 基于混合高斯模型的固定时间约束的航天器集群控制方法 |
CN112415896B (zh) * | 2020-11-24 | 2022-05-06 | 西北工业大学 | 一种基于自适应超扭滑模控制的平动点轨道保持方法 |
CN113220003B (zh) * | 2021-03-31 | 2022-08-26 | 西北工业大学 | 一种非合作柔性组合体航天器姿态稳定混合非脆弱控制方法 |
CN113282097B (zh) * | 2021-06-04 | 2022-07-29 | 中国人民解放军战略支援部队航天工程大学 | 一种geo博弈航天器相对位置非球形摄动误差的控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102880184A (zh) * | 2012-10-24 | 2013-01-16 | 北京控制工程研究所 | 一种静止轨道卫星自主轨道控制方法 |
CN102878995A (zh) * | 2012-10-24 | 2013-01-16 | 北京控制工程研究所 | 一种静止轨道卫星自主导航方法 |
CN105718659A (zh) * | 2016-01-21 | 2016-06-29 | 西北工业大学 | 一种高面质比航天器轨道动力学分析方法 |
-
2016
- 2016-09-28 CN CN201610860433.3A patent/CN106628257B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102880184A (zh) * | 2012-10-24 | 2013-01-16 | 北京控制工程研究所 | 一种静止轨道卫星自主轨道控制方法 |
CN102878995A (zh) * | 2012-10-24 | 2013-01-16 | 北京控制工程研究所 | 一种静止轨道卫星自主导航方法 |
CN105718659A (zh) * | 2016-01-21 | 2016-06-29 | 西北工业大学 | 一种高面质比航天器轨道动力学分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106628257A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106628257B (zh) | 地球摄动引力场中近地航天器相对运动轨道的保持方法 | |
CN103728980B (zh) | 航天器相对轨道的控制方法 | |
CN101320524B (zh) | 多处理器实时仿真平台 | |
CN102819266B (zh) | 一种拟周期j2不变相对轨道编队飞行控制方法 | |
CN102073755A (zh) | 近空间高超声速飞行器运动控制仿真方法 | |
CN104898642A (zh) | 一种用于航天器姿态控制算法的集成测试仿真系统 | |
CN104309822B (zh) | 一种基于参数优化的航天器单脉冲水滴形绕飞轨迹悬停控制方法 | |
CN103970964A (zh) | 一种挠性卫星模态参数在轨辨识方法 | |
CN103439975A (zh) | 一种分布式指数时变滑模姿态协同跟踪控制方法 | |
Hua et al. | Effect of elastic deformation on flight dynamics of projectiles with large slenderness ratio | |
CN105930305B (zh) | 一种三脉冲交会接近制导方法 | |
CN105631099A (zh) | 一种小天体探测器着陆动力学模拟系统 | |
CN106156414B (zh) | 一种卫星轨道仿真方法及装置 | |
Liu et al. | Collision-free trajectory design for long-distance hopping transfer on asteroid surface using convex optimization | |
Lee et al. | Lagrangian mechanics and Lie group variational integrators for spacecraft with imbalanced reaction wheels | |
Manchester et al. | Recursive inertia estimation with semidefinite programming | |
CN105116905A (zh) | 一种飞行器姿态控制方法 | |
CN106709161A (zh) | 一种获取航天器贮箱内液体燃料大幅晃动作用力的近似方法 | |
CN103514362A (zh) | 基于模型误差补偿的两行根数生成方法 | |
Zhao et al. | An adaptive optimization algorithm based on clustering analysis for return multi-flight-phase of VTVL reusable launch vehicle | |
CN105912020A (zh) | 一种弱稳定边界区域内航天器编队飞行方法 | |
CN110321598B (zh) | 一种j2项摄动条件下航天器相对运动解析求解方法 | |
Traub et al. | Assessment of the dependencies of realistic differential drag controlled in-plane reconfiguration maneuvers on relevant parameters | |
CN107300861B (zh) | 一种航天器动力学分布式计算方法 | |
CN113741197B (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 |