CN105115512A - 一种火星大气进入段侧向预测校正制导方法 - Google Patents
一种火星大气进入段侧向预测校正制导方法 Download PDFInfo
- Publication number
- CN105115512A CN105115512A CN201510611464.0A CN201510611464A CN105115512A CN 105115512 A CN105115512 A CN 105115512A CN 201510611464 A CN201510611464 A CN 201510611464A CN 105115512 A CN105115512 A CN 105115512A
- Authority
- CN
- China
- Prior art keywords
- detector
- guidance
- angle
- heel
- energy
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Toys (AREA)
Abstract
本发明公开的一种火星大气进入段侧向预测校正制导方法,涉及侧向预测校正制导方法,属于深空探测技术领域。本发明在已有纵向预测校正制导方法的基础上,将预测校正方法引入探测器侧向运动的制导律设计中。首先,确定着陆探测任务所需的侧向运动约束条件,以确定探测器在进入段所需倾侧角反转次数;然后,求解侧向运动约束条件确定探测器倾侧角反转时刻的能量,在探测器能量超过倾侧角反转时刻能量时,进行倾侧角反转。在达到该约束条件后进行下一次倾侧角反转能量确定;最后,结合倾侧角反转时刻能量以及纵向制导律确定当前制导周期最终制导输出。本发明可在保证开伞点位置精度同时,根据任务需要对进入轨迹的侧向运动进行灵活规划,且节约燃料。
Description
技术领域
本发明涉及一种侧向预测校正制导方法,尤其涉及一种火星大气进入段侧向预测校正制导方法,属于深空探测技术领域。
背景技术
在未来火星着陆探测任务对实现火星表面的着陆精度提出了更高的要求。而火星大气进入段是整个火星进入、下降与着陆过程中历时最久,环境最恶劣的一个阶段,保证进入段任务的成功实施,是对于保证最终着陆的精度起着至关重要的作用。
到目前为止,共有7颗火星探测器成功着陆火星表面。最近一次着陆探测任务“火星科学实验室/好奇号”任务中,首次在进入段采用了制导律,并在预定的开伞区域内实现安全开伞。这表明随着人类火星探测的深入开展,在未来火星着陆探测任务中,尤其是诸如载人火星探测、采样返回等任务对着陆精度要求的提高,采用相应的进入制导律成为一种趋势。在火星大气进入段实施相应的制导方法,是保证进入段末端精度的重要手段。
目前所研究的火星大气进入段制导方法主要可分为标称轨迹法和预测校正制导法。其中,标称轨迹法通过离线或在线规划出一条能够满足着陆精度的标称轨迹,并由相应的轨迹跟踪律来跟踪该标称轨迹,以消除进入点偏差以及进入过程中各种干扰或不确定因素造成的偏差,从而保证开伞的位置精度。而预测校正制导法则是根据当前状态,在线预测落点偏差,根据落点偏差产生相应的制导指令,以修正各种因素造成的落点偏差,进而保证开伞精度。当前预测校正制导方法主要针对纵向航程设计相应的预测校正制导律,而侧向运动是通过相应的航向角误差走廊或者侧向航程误差走廊加以约束,当航向角或者侧向航程超出该走廊时,倾侧角产生反转,进而保证侧向运动的末端位置精度。而未来火星探测任务对进入轨迹的侧向运动部分提出了新的要求,例如倾侧角反转过于频繁会导致燃料的过多消耗,进而增加任务失败的几率;未来火星大气进入段无线电信标的布置,探测器飞越特定区域能够显著提高导航系统的可观测度,进而提高导航精度。对于这些需求,传统的通过规划误差走廊来约束侧向运动的方式缺少在规划侧向运动方面的灵活性。
为了满足未来火星探测任务对进入段侧向运动的需求,有必要针对火星大气进入段探测器的侧向运动,设计一种相应的预测校正制导律,以保证开伞点位置精度的同时,能够根据任务需要对进入轨迹的侧向运动进行灵活规划。
发明内容
本发明公开的一种火星大气进入段侧向预测校正制导方法所解决的技术问题是,实现在保证开伞点位置精度的同时,能够根据任务需要对进入轨迹的侧向运动进行灵活规划。
本发明的目的是通过下述技术方案实现的:
本发明公开的一种火星大气进入段侧向预测校正制导方法,包括如下步骤:
在已有纵向预测校正制导方法的基础上,将预测校正方法引入探测器侧向运动的制导律设计中。在每个制导周期内:首先,确定着陆探测任务所需的探测器的侧向运动的约束条件,以确定探测器在进入段所需倾侧角反转的次数;然后,利用数值方法求解侧向运动的约束条件,确定探测器倾侧角反转时刻的能量,在探测器的能量超过倾侧角反转时刻的能量时,进行倾侧角反转。在达到该约束条件后,开始进行下一次倾侧角反转能量的确定;最后,结合倾侧角反转时刻的能量以及纵向制导律确定当前制导周期最终的制导输出。进而实现在保证开伞点位置精度的同时,能够根据任务需要对进入轨迹的侧向运动进行灵活规划。
本发明公开的一种火星大气进入段侧向预测校正制导方法,包括如下步骤:
步骤1、确定纵向运动制导指令|σ|。
利用探测器动力学模型进行数值积分至满足开伞条件,得到开伞时刻的开伞剩余纵程与目标位置的开伞精度偏差sf。所述的开伞条件指探测器动压在区间[qmin,qmax]内和探测器马赫数在区间[Mamin,Mamax]内。具体实现方法为:
考虑火星自转影响的探测器对无量纲时间的三自由度无量纲进入动力学模型为
其中,s为剩余纵程,表征从探测器当前位置到目标开伞位置的火星表面大圆弧的距离,r为火星质心到探测器质心的距离,无量纲参数为火星半径R0,v探测器相对于火星的速度,无量纲参数为其中g0为火星表面重力加速度,γ为航迹角,σ为倾侧角,g为当地重力加速度,无量纲参数为g0。D和L分别阻力加速度和升力加速度
其无量纲参数均为g0,CD和CL分别为阻力系数和升力系数,S为探测器参考面积,m为探测器质量,q=ρv2/2为动压,β=m/SCD为探测器弹道系数,L/D为探测器升阻比。火星大气密度采用指数模型
其中ρ0为参考密度,h0为参考高度,hs为大气密度标高。
定义进入段飞行器的比能量
定义火星大气进入段的侧向航程
χ=R0sin-1(sinStogosin△ψ)(5)
其中,△ψ为航向角偏差,Stogo为剩余航程,由式(6)给出
Stogo=R0cos-1[sinφtsinφ+cosφtcosφcos(θt-θ)](6)
纵向运动的动力学由公式(7)进行描述
给出倾侧角剖面参数化形式
其中,σf为开伞时刻的倾侧角;e,e0和ef分别为当前时刻、初始时刻以及开伞时刻的能量。σ0>0为进入时刻的倾侧角,通过数值求解非线性方程(9)求取σ0。
再利用式(8)产生每一步的纵向制导指令|σ|,从而使探测器满足开伞点的位置精度。
步骤2、确定侧向运动的约束条件及相应倾侧角反转次数irev。
侧向运动的约束条件为方程组(10)
χ(etar,erev)=χ*(10)
其中,向量χ*为对应能量为向量etar时,侧向航程的取值;向量erev为倾侧角反转时刻的能量向量,所述的向量erev是侧向预测校正制导的待求参量。方程组(10)的分量形式为
倾侧角反转次数为irev=dim(erev)=dim(χ)时,方程组(10)有唯一解。侧向运动的相应约束条件由方程组(10)给定。
步骤3、求解约束条件求得倾侧角反转时刻的能量erev。
在求解方程组(10)的过程中,方程组(10)各个分量方程式(11)相互独立,能够分别独立求解。方程组(10)构成关于能量向量erev的非线性方程。通过采用数值方法对约束条件方程组(10)进行迭代求解,求得倾侧角反转时刻的能量erev。
步骤4、求解出侧向制导指令sign(σk(e))。
在探测器进入过程中,每当能量e依次超过制导律所求解出的能量向量erev的各个分量时,倾侧角σ便反转一次。数学表达式为式(12)
其中,sign(σ0(e))=-sign(χ0)。
步骤5、求解出第k次制导的制导指令输出
第k次制导的制导指令输出由第k次纵向制导指令|σ(k)|及第k次侧向制导指令sign(σ(k))得到
步骤6、重复步骤1-5实时更新制导指令,直至探测器的动压和马赫数满足相应的开伞条件,火星大气进入段制导过程结束。
有益效果:
本发明公开的一种火星大气进入段侧向预测校正制导方法,首次将预测校正的方法引入火星大气进入段侧向运动的制导律设计中,该制导方法根据任务需要,确定侧向运动所需满足的约束条件,进而确定倾侧角反转次数,实现在保证开伞点位置精度的同时,能够根据任务需要对进入轨迹的侧向运动进行灵活规划。
附图说明
图1为侧向预测校正制导指令生成流程图;
图2为倾侧角反转示意;
图3给出了基础算例和实施例1的轨迹;
图4给出了基础算例和实施例2的轨迹。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合两个实施例和相应附图对发明内容做进一步说明。
两个实施例中采用的大气进入的状态初值由表1给出。开伞点为(θf,φf)=(-40.1°,-43.8°)。
表1仿真初值的设定(归一化结果)
r0 | θ0 | φ0 | v0 | γ0 | ψ0 |
1.0368 | -1.016 | -0.723 | 1.5487 | -0.226 | 1.871 |
为方便与实施例中的结果进行对比,给出一个侧向运动制导方法的基础算例。
本基础算例所采用的侧向运动制导方法是,规划与相对速度成线性关系的侧向航程走廊式(14),当侧向航程超出该走廊时进行倾侧角反转。
χc=c1v+c0(14)
其中,c0=8.71×10-5,c1=5.21×10-3。
实施例1:
本实施例公开的一种火星大气进入段侧向预测校正制导方法,用以实现在火星大气进入段进行一次侧向机动,所述的侧向机动即倾侧角反转机动,从而达到节省燃料的目的。包括如下步骤:
利用探测器动力学模型进行数值积分至满足开伞条件(及开伞动压区间[qmin,qmax]和开伞马赫数区间[Mamin,Mamax]),计算开伞时刻的开伞剩余纵程与目标位置的开伞精度偏差sf,具体实现方法为:
步骤1、确定纵向运动制导指令|σ|。
考虑火星自转影响的探测器对无量纲时间的三自由度无量纲进入动力学模型为
其中,s为剩余纵程,表征从探测器当前位置到目标开伞位置的火星表面大圆弧的距离,r为火星质心到探测器质心的距离,无量纲参数为火星半径R0,v探测器相对于火星的速度,无量纲参数为其中g0为火星表面重力加速度,γ为航迹角,σ为倾侧角,g为当地重力加速度,其无量纲参数为g0。D和L分别阻力加速度和升力加速度
其无量纲参数均为g0,CD和CL分别为阻力系数和升力系数,S为探测器参考面积,m为探测器质量,q=ρv2/2为动压,β=m/SCD为探测器弹道系数,L/D为探测器升阻比。火星大气密度采用指数模型
其中ρ0=2×10-4kg/m3为参考密度,h0=40000m为参考高度,hs=7500km为大气密度标高。
定义进入段飞行器的比能量
定义火星大气进入段的侧向航程
χ=R0sin-1(sinStogosin△ψ)(19)
其中,△ψ为航向角偏差,Stogo为剩余航程,由式(20)给出
Stogo=R0cos-1[sinφtsinφ+cosφtcosφcos(θt-θ)](20)
纵向运动的动力学由公式(21)进行描述
给出倾侧角剖面参数化形式
其中,σf=70°为开伞时刻的倾侧角;e,e0和ef分别为当前时刻、初始时刻以及开伞时刻的能量,其中e0=-0.2348,ef=0.9875。σ0>0为进入时刻的倾侧角,通过数值求解非线性方程(23)求取σ0。
再利用式(22)产生每一步的纵向制导指令|σ|,从而使探测器满足开伞点的位置精度。
步骤2、确定侧向运动的约束条件及相应倾侧角反转次数irev。
本实施例所要实现的目标是在整个火星大气进入段中,探测器只进行一次侧向机动,则侧向运动的约束条件为方程组(24)
χ(etar,erev)=χ*(24)
可以写作
χ(ef,erev)=0(25)
即只在e=erev时进行一次倾侧角反转机动,满足开伞点处的位置精度要求χ(ef)=0(在etar=ef时,χ*=0)。由于只进行一次侧向机动,侧向运动的约束条件数为1,相应的倾侧角翻转次数irev=1。
步骤3、求解约束条件求得倾侧角反转时刻的能量erev。
约束方程(25)构成了关于倾侧角反转时刻能量erev的非线性方程。通过采用Newton-Raphson方法对约束条件方程(25)进行迭代求解,即
从而求得倾侧角反转时刻的能量erev。本实施例下,倾侧角反转时刻的能量为erev=0.482,对应倾侧角反转时刻为trev=115s。
步骤4、求解出侧向制导指令sign(σk(e))。
在探测器进入过程中,每当能量e依次超过制导律所求解出的能量向量erev的各个分量时,倾侧角σ便倾侧角反转一次。数学表达式为式(27)
其中,sign(σ0(e))=-sign(χ0)。
由于只进行一次侧向机动,因此在机动前,保持sign(σ(e))=-sign(χ0);在侧向机动后,保持sign(σ(e))=sign(χ0)。
步骤5、求解出第k次制导的制导指令输出
第k次制导的制导指令输出由第k次纵向制导指令|σ(k)|及第k次侧向制导指令sign(σ(k))得到
步骤6、重复步骤1-5实时更新制导指令,直至探测器的动压和马赫数满足相应的开伞条件,火星大气进入段制导过程结束。
图3给出了基础算例和本实施例的轨迹。可以看出,本实施例中只进行了一次倾侧角反转,从而达到节省能量的目的,并且保证了末端开伞时刻的位置精度。
实施例2:
本实施例公开的一种火星大气进入段侧向预测校正制导方法,针对未来火星大气进入段无线电导航方案时,飞行器飞越相应可观区域会显著提高导航系统的可观测度,进而提高导航精度。因此,期望飞行器飞越可观区域的轨迹尽可能长。为简单起见,此处假设飞行器飞越可观区域的形心(本实施例中选取可观区域的形心坐标为(θ1,φ1)=(-45.14°,-43.67°))时,其飞越可观区域的轨迹最长。
本实施例包括如下步骤:
利用探测器动力学模型进行数值积分至满足开伞条件(及开伞动压区间[qmin,qmax]和开伞马赫数区间[Mamin,Mamax]),计算开伞时刻的开伞剩余纵程与目标位置的开伞精度偏差sf,具体实现方法为:
步骤1、确定纵向运动制导指令|σ|。
考虑火星自转影响的探测器对无量纲时间的三自由度无量纲进入动力学模型为
其中,s为剩余纵程,表征从探测器当前位置到目标开伞位置的火星表面大圆弧的距离,r为火星质心到探测器质心的距离,无量纲参数为火星半径R0,v探测器相对于火星的速度,无量纲参数为其中g0为火星表面重力加速度,γ为航迹角,σ为倾侧角,g为当地重力加速度,其无量纲参数为g0。D和L分别阻力加速度和升力加速度
其无量纲参数均为g0,CD和CL分别为阻力系数和升力系数,S为探测器参考面积,m为探测器质量,q=ρv2/2为动压,β=m/SCD为探测器弹道系数,L/D为探测器升阻比。火星大气密度采用指数模型
其中ρ0=2×10-4kg/m3为参考密度,h0=40000m为参考高度,hs=7500km为大气密度标高。
定义进入段飞行器的比能量
定义火星大气进入段的侧向航程
χ=R0sin-1(sinStogosin△ψ)(33)
其中,△ψ为航向角偏差,Stogo为剩余航程,由式(34)给出
Stogo=R0cos-1[sinφtsinφ+cosφtcosφcos(θt-θ)](34)
纵向运动的动力学由公式(35)进行描述
给出倾侧角剖面参数化形式
其中,σf=70°为开伞时刻的倾侧角;e,e0和ef分别为当前时刻、初始时刻以及开伞时刻的能量,其中e0=-0.2348,ef=0.9875。σ0>0为进入时刻的倾侧角,通过数值求解非线性方程(37)求取σ0。
再利用式(36)产生每一步的纵向制导指令|σ|,从而使探测器满足开伞点的位置精度。
步骤2、确定侧向运动的约束条件及相应倾侧角反转次数irev。
本实施例所要实现的目标是在整个火星大气进入段中,探测器只进行一次侧向机动,则侧向运动的约束条件为方程组(38)
χ(etar,erev)=χ*(38)
可以写作
e1,分别为由纵向运动得到的可观区域形心处的能量及其侧向航程距离。本实施例中,可观区域形心对应的侧向航程为χ1=15km。能量向量erev的维数dim(erev)=2。本实施例中,ef=0.9875。此时侧向运动的约束条件数为2,相应的倾侧角翻转次数irev=2。
步骤3、求解约束条件求得倾侧角反转时刻的能量erev。
约束方程(25)构成了关于倾侧角反转时刻能量向量erev的非线性方程。通过采用Newton-Raphson方法对约束条件方程(39)进行迭代求解,即
求得倾侧角反转时刻的能量erev。本实施例下,倾侧角反转时刻的能量为erev=(0.189,0.891),对应倾侧角反转时刻为trev=(101s,165s)。
步骤4、求解出侧向制导指令sign(σk(e))。
在探测器进入过程中,每当能量e依次超过制导律所求解出的能量向量erev的各个分量时,倾侧角σ便反转一次。数学表达式为式(41)
其中,sign(σ0(e))=-sign(χ0)。
步骤5、求解出第k次制导的制导指令输出
第k次制导的制导指令输出由第k次纵向制导指令|σ(k)|及第k次侧向制导指令sign(σ(k))得到
步骤6、重复步骤1-5实时更新制导指令,直至探测器的动压和马赫数满足相应的开伞条件,火星大气进入段制导过程结束。
图4给出了基础算例和本实施例的轨迹。可以看出,经过两次倾侧角反转,实现了飞行器经过特定可观区域的形心,并且保证了末端开伞时刻的位置精度。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下,可对这些实施例进行多种变化、修改、替换和变形,本发明的范围由权利要求及其等同物限定。
Claims (3)
1.一种火星大气进入段侧向预测校正制导方法,其特征在于:在已有纵向预测校正制导方法的基础上,将预测校正方法引入探测器侧向运动的制导律设计中;在每个制导周期内,首先,确定着陆探测任务所需的探测器的侧向运动的约束条件,以确定探测器在进入段所需倾侧角反转的次数;然后,利用数值方法求解侧向运动的约束条件,确定探测器倾侧角反转时刻的能量,在探测器的能量超过倾侧角反转时刻的能量时,进行倾侧角反转;在达到约束条件后,开始进行下一次倾侧角反转能量的确定;最后,结合倾侧角反转时刻的能量以及纵向制导律确定当前制导周期最终的制导输出;进而实现在保证开伞点位置精度的同时,能够根据任务需要对进入轨迹的侧向运动进行灵活规划。
2.如权利要求1所述的一种火星大气进入段侧向预测校正制导方法,其特征在于:具体实现方法包括如下步骤,
步骤1、确定纵向运动制导指令|σ|;
利用探测器动力学模型进行数值积分至满足开伞条件,得到开伞时刻的开伞剩余纵程与目标位置的开伞精度偏差sf;所述的开伞条件指探测器动压在区间[qmin,qmax]内和探测器马赫数在区间[Mamin,Mamax]内;具体实现方法为,
考虑火星自转影响的探测器对无量纲时间的三自由度无量纲进入动力学模型为
其中,s为剩余纵程,表征从探测器当前位置到目标开伞位置的火星表面大圆弧的距离,r为火星质心到探测器质心的距离,无量纲参数为火星半径R0,v探测器相对于火星的速度,无量纲参数为其中g0为火星表面重力加速度,γ为航迹角,σ为倾侧角,g为当地重力加速度,无量纲参数为g0;D和L分别阻力加速度和升力加速度
其无量纲参数均为g0,CD和CL分别为阻力系数和升力系数,S为探测器参考面积,m为探测器质量,q=ρv2/2为动压,β=m/SCD为探测器弹道系数,L/D为探测器升阻比;火星大气密度采用指数模型
其中ρ0为参考密度,h0为参考高度,hs为大气密度标高;
定义进入段飞行器的比能量
定义火星大气进入段的侧向航程
χ=R0sin-1(sinStogosin△ψ)(5)
其中,△ψ为航向角偏差,Stogo为剩余航程,由式(6)给出
Stogo=R0cos-1[sinφtsinφ+cosφtcosφcos(θt-θ)](6)
纵向运动的动力学由公式(7)进行描述
给出倾侧角剖面参数化形式
其中,σf为开伞时刻的倾侧角;e,e0和ef分别为当前时刻、初始时刻以及开伞时刻的能量;σ0>0为进入时刻的倾侧角,通过数值求解非线性方程(9)求取σ0;
再利用式(8)产生每一步的纵向制导指令|σ|,从而使探测器满足开伞点的位置精度;
步骤2、确定侧向运动的约束条件及相应倾侧角反转次数irev;
侧向运动的约束条件为方程组(10)
χ(etar,erev)=χ*(10)
其中,向量χ*为对应能量为向量etar时,侧向航程的取值;向量erev为倾侧角反转时刻的能量向量,所述的向量erev是侧向预测校正制导的待求参量;方程组(10)的分量形式为
倾侧角反转次数为irev=dim(erev)=dim(χ)时,方程组(10)有唯一解;侧向运动的相应约束条件由方程组(10)给定;
步骤3、求解约束条件求得倾侧角反转时刻的能量erev;
在求解方程组(10)的过程中,方程组(10)各个分量方程式(11)相互独立,能够分别独立求解;方程组(10)构成关于能量向量erev的非线性方程;通过采用数值方法对约束条件方程组(10)进行迭代求解,求得倾侧角反转时刻的能量erev;
步骤4、求解出侧向制导指令sign(σk(e));
在探测器进入过程中,每当能量e依次超过制导律所求解出的能量向量erev的各个分量时,倾侧角σ便反转一次;数学表达式为式(12)
其中,sign(σ0(e))=-sign(χ0);
步骤5、求解出第k次制导的制导指令输出
第k次制导的制导指令输出由第k次纵向制导指令|σ(k)|及第k次侧向制导指令sign(σ(k))得到
步骤6、重复步骤1-5实时更新制导指令,直至探测器的动压和马赫数满足相应的开伞条件,火星大气进入段制导过程结束。
3.如权利要求2所述的一种火星大气进入段侧向预测校正制导方法,其特征在于:步骤3所述的通过采用数值方法对约束条件方程组(10)进行迭代求解,所述的数值方法采用Newton-Raphson方法对约束条件方程组(10)进行迭代求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510611464.0A CN105115512B (zh) | 2015-09-23 | 2015-09-23 | 一种火星大气进入段侧向预测校正制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510611464.0A CN105115512B (zh) | 2015-09-23 | 2015-09-23 | 一种火星大气进入段侧向预测校正制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105115512A true CN105115512A (zh) | 2015-12-02 |
CN105115512B CN105115512B (zh) | 2017-10-10 |
Family
ID=54663557
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510611464.0A Expired - Fee Related CN105115512B (zh) | 2015-09-23 | 2015-09-23 | 一种火星大气进入段侧向预测校正制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105115512B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105511493A (zh) * | 2015-12-28 | 2016-04-20 | 北京理工大学 | 一种基于火星大气辅助的低轨星座部署方法 |
CN105718660A (zh) * | 2016-01-21 | 2016-06-29 | 中国工程物理研究院总体工程研究所 | 临近空间大范围机动弹道三维包络计算方法 |
CN105867402A (zh) * | 2016-05-10 | 2016-08-17 | 北京航空航天大学 | 一种火星着陆器大气进入段抗干扰复合在线制导方法 |
CN107323691A (zh) * | 2017-07-04 | 2017-11-07 | 北京理工大学 | 一种多约束火星大气进入预测制导方法 |
CN107804487A (zh) * | 2017-10-16 | 2018-03-16 | 中国西安卫星测控中心 | 一种基于自适应偏差控制的跳跃式再入返回落点预报方法 |
CN107942673A (zh) * | 2017-12-11 | 2018-04-20 | 北京理工大学 | 一种开伞点高度跟踪的火星大气进入段解析制导方法 |
CN107966156A (zh) * | 2017-11-24 | 2018-04-27 | 北京宇航系统工程研究所 | 一种适用于运载火箭垂直回收段的制导律设计方法 |
CN111319794A (zh) * | 2020-02-25 | 2020-06-23 | 上海航天控制技术研究所 | 一种适用于火星探测制动捕获期间推进自主故障处理方法 |
CN113734474A (zh) * | 2021-08-09 | 2021-12-03 | 北京空间飞行器总体设计部 | 一种多约束下大底分离触发条件确定方法和系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100250031A1 (en) * | 2009-03-27 | 2010-09-30 | The Charles Stark Draper Laboratory, Inc. | Propulsive guidance for atmospheric skip entry trajectories |
CN103863579A (zh) * | 2014-03-31 | 2014-06-18 | 北京控制工程研究所 | 一种深空探测返回过程的预测校正制导方法 |
CN104035335A (zh) * | 2014-05-27 | 2014-09-10 | 北京航空航天大学 | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导律 |
CN104635488A (zh) * | 2014-12-12 | 2015-05-20 | 北京理工大学 | 一种火星大气进入段鲁棒预测制导方法 |
CN104648695A (zh) * | 2014-12-18 | 2015-05-27 | 北京控制工程研究所 | 一种基于倾侧角可用性的再入走廊最优规划方法 |
-
2015
- 2015-09-23 CN CN201510611464.0A patent/CN105115512B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100250031A1 (en) * | 2009-03-27 | 2010-09-30 | The Charles Stark Draper Laboratory, Inc. | Propulsive guidance for atmospheric skip entry trajectories |
CN103863579A (zh) * | 2014-03-31 | 2014-06-18 | 北京控制工程研究所 | 一种深空探测返回过程的预测校正制导方法 |
CN104035335A (zh) * | 2014-05-27 | 2014-09-10 | 北京航空航天大学 | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导律 |
CN104635488A (zh) * | 2014-12-12 | 2015-05-20 | 北京理工大学 | 一种火星大气进入段鲁棒预测制导方法 |
CN104648695A (zh) * | 2014-12-18 | 2015-05-27 | 北京控制工程研究所 | 一种基于倾侧角可用性的再入走廊最优规划方法 |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105511493A (zh) * | 2015-12-28 | 2016-04-20 | 北京理工大学 | 一种基于火星大气辅助的低轨星座部署方法 |
CN105511493B (zh) * | 2015-12-28 | 2018-04-03 | 北京理工大学 | 一种基于火星大气辅助的低轨星座部署方法 |
CN105718660B (zh) * | 2016-01-21 | 2019-03-01 | 中国工程物理研究院总体工程研究所 | 临近空间大范围机动弹道三维包络计算方法 |
CN105718660A (zh) * | 2016-01-21 | 2016-06-29 | 中国工程物理研究院总体工程研究所 | 临近空间大范围机动弹道三维包络计算方法 |
CN105867402A (zh) * | 2016-05-10 | 2016-08-17 | 北京航空航天大学 | 一种火星着陆器大气进入段抗干扰复合在线制导方法 |
CN105867402B (zh) * | 2016-05-10 | 2019-05-07 | 北京航空航天大学 | 一种火星着陆器大气进入段抗干扰复合在线制导方法 |
CN107323691A (zh) * | 2017-07-04 | 2017-11-07 | 北京理工大学 | 一种多约束火星大气进入预测制导方法 |
CN107323691B (zh) * | 2017-07-04 | 2020-01-07 | 北京理工大学 | 一种多约束火星大气进入预测制导方法 |
CN107804487A (zh) * | 2017-10-16 | 2018-03-16 | 中国西安卫星测控中心 | 一种基于自适应偏差控制的跳跃式再入返回落点预报方法 |
CN107966156A (zh) * | 2017-11-24 | 2018-04-27 | 北京宇航系统工程研究所 | 一种适用于运载火箭垂直回收段的制导律设计方法 |
CN107966156B (zh) * | 2017-11-24 | 2020-09-18 | 北京宇航系统工程研究所 | 一种适用于运载火箭垂直回收段的制导律设计方法 |
CN107942673A (zh) * | 2017-12-11 | 2018-04-20 | 北京理工大学 | 一种开伞点高度跟踪的火星大气进入段解析制导方法 |
CN111319794A (zh) * | 2020-02-25 | 2020-06-23 | 上海航天控制技术研究所 | 一种适用于火星探测制动捕获期间推进自主故障处理方法 |
CN111319794B (zh) * | 2020-02-25 | 2021-10-01 | 上海航天控制技术研究所 | 一种适用于火星探测制动捕获期间推进自主故障处理方法 |
CN113734474A (zh) * | 2021-08-09 | 2021-12-03 | 北京空间飞行器总体设计部 | 一种多约束下大底分离触发条件确定方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105115512B (zh) | 2017-10-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105115512A (zh) | 一种火星大气进入段侧向预测校正制导方法 | |
CN104615144B (zh) | 基于目标规划的高超声速飞行器再入轨迹在线优化方法 | |
US9671789B1 (en) | Method and system for anti-disturbance composite on-line guidance for atmosphere-entering phase of a Mars lander | |
CN110262495A (zh) | 可实现移动机器人自主导航与精确定位的控制系统及方法 | |
CN105005313B (zh) | 一种基于路径点规划的火星大气进入段预测制导方法 | |
CN106054604B (zh) | 基于模型预测控制理论的再入飞行器鲁棒最优制导方法 | |
CN109250153B (zh) | 火星大气进入段轨迹最优跟踪制导方法 | |
CN103926835B (zh) | 一种基于干扰观测器的着陆器动力下降段优化控制方法 | |
CN104035335A (zh) | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导律 | |
CN105573337B (zh) | 一种满足再入角和航程约束的离轨制动闭路制导方法 | |
CN106371312B (zh) | 基于模糊控制器的升力式再入预测-校正制导方法 | |
CN109839116B (zh) | 一种速度增量最小的交会接近方法 | |
CN104019818A (zh) | 一种基于预测轨迹的行星导航轨道器布局优化方法 | |
CN105807781A (zh) | 一种基于比例导引的空间近距飞越末制导方法 | |
CN110530374A (zh) | 一种无人机天钩回收阶段航迹偏差显示方法 | |
CN103955224A (zh) | 一种用于相对运动视线跟踪的姿态控制方法 | |
CN107792404A (zh) | 一种交会对接平移靠拢段的安全带确定方法 | |
Liang et al. | Decoupling trajectory tracking for gliding reentry vehicles | |
CN105678061B (zh) | 一种满足被动安全要求的空间站安全区域设计方法 | |
CN114138003B (zh) | 一种小型固定翼无人机双机协同区域侦察方法 | |
JP2015168315A (ja) | 誘導装置、誘導装置搭載宇宙機 | |
CN104898678A (zh) | 一种运载火箭飞行控制的冗余诊断方法 | |
Khuralay et al. | Computer simulation of intelligent control systems for high-precision cruise missiles | |
CN110562492B (zh) | 探测器火星大气进入轨迹快速生成方法 | |
Ostroumov | Navaids facility for aircraft positioning |
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: 20171010 Termination date: 20210923 |