CN103662092A - 一种衔接主减速和接近段的预测校正方法 - Google Patents

Abstract

一种衔接主减速和接近段的预测校正方法，步骤为：计算主减速段制导参数；基于主减速段制导参数和探测器当前状态进行快速调整段终端预测；基于接近段入口条件修正主减速段制导目标；判断是否满足主减速切换条件，如果满足，则切换到快速调整制导；否则重新计算步骤（1）～（3），直到满足主减速切换条件；利用切换时刻状态确定快速调整段制导参数；根据快速调整段制导参数确定当前制导指令；利用当前状态规划接近段入口制导指令；快速调整段当前制导指令与接近段入口制导指令夹角小于设定值或制导时间完成，切换到接近段制导。本发明保证了接近段入口对姿态、高度、速度和加速度的需求，满足了探测器从主减速段平缓过渡到接近段的要求。

Description

一种衔接主减速和接近段的预测校正方法

技术领域

本发明涉及一种衔接主减速和接近段的预测校正方法，主要应用于深空着陆探测器，属于航天器制导、导航与控制技术领域。可应用于月球以及火星、小行星等深空天体探测任务，具有广泛的应用价值和市场前景。

背景技术

深空天体软着陆动力下降过程一般可分为六个任务段：（1）主减速段：距天体表面高度从约15km到约3km，该段主要任务是软着陆制动，减小探测器的速度至预设值，高度下降至约3km。（2）快速调整段：距天体表面高度从约3km到约2.4km，该段主要任务是快速衔接主减速和接近段。（3）接近段：距天体表面高度从约2.4km到约100m，该段主要任务是粗避障。（4）悬停段：距天体表面高度约100m，该段主要任务是对着陆区域的精障碍检测。保持探测器处于悬停状态，利用三维成像敏感器对着陆区进行观测，选择出安全着陆点。（5）避障段：距天体表面高度从约100m到约30m，该段主要任务是精避障和下降。（6）缓速下降段：距天体表面高度从约30m到0m，该段主要任务是保证探测器平稳缓速下降到天体表面。

Luna和Surveyor系列月球探测器分别采用了标称轨迹和重力转弯制导；Apollo主减速段和接近段都采用四次多项式制导，通过在地面优化设计参数来满足推力变化过程的衔接，但是对于姿态的衔接是没有考虑；美国的重返月球计划提出，在主减速段采用PEG制导，接近段采用多项式制导，主减速过渡到接近段是通过快速姿态调整实现，没有提出具有衔接功能的自适应修正制导方法。

对于动力下降的主减速段末期，探测器姿态仍接近水平，主发动机仍工作在最大推力段，推力加速度也达到最大；而动力下降的接近段要求探测器姿态接近垂直，主发动机工作在低推力水平上，高度、速度和加速度满足一定关系；可见，主减速段末端状态和接近段初始状态很难直接衔接上。

发明内容

本发明技术解决问题：克服现有技术的不足，提供一种衔接主减速和接近段的预测校正方法，考虑到主减速段末端状态和接近段初始状态衔接的需求，设计了主减速段带有制导目标修正功能的制导律和基于切换状态确定制导参数的自适应制导律，保证了接近段入口对姿态、高度、速度和加速度的需求，满足了探测器从主减速段平缓过渡到接近段的要求。

本发明技术解决方案：一种衔接主减速和接近段的预测校正方法，实现步骤如下：

（1）计算主减速段制导参数

利用探测器当前状态和主发动机比冲和推力参数，计算主减速段制导参数：剩余主减速时间、参考时间、推力方向变化率和参考速度增量方向；

（2）基于主减速段制导参数和探测器当前状态进行快速调整段终端预测

利用步骤（1）中计算主减速段制导参数和探测器当前状态，预测主减速终端状态和制导指令；根据接近段入口对姿态和加速度的要求，确定快速调整段终端状态；由此确定快速调整段制导参数，进而利用快速调整段制导参数预测快速调整过程产生的位置和速度变化量，最终预测得到快速调整终端状态，即预测得到接近段入口状态；

（3）基于接近段入口条件修正主减速段制导目标

利用预测的接近段初始状态与接近段入口目标状态进行比较，得到接近段入口状态的差，所述接近段入口目标状态即接近段入口对高度和速度的要求；再利用状态转移矩阵递推到主减速末端，得到主减速段制导参数的修正量，进而达到反馈修正的目的；经过主减速末期的多次预测修正后，可以保证接近段入口目标状态的实现；

（4）判断是否满足主减速切换条件，如果满足，则切换到快速调整制导；否则重新计算步骤（1）～（3），直到满足主减速切换条件；

（5）利用切换时刻的探测器推力加速度大小、方向和接近段入口目标状态确定快速调整段制导参数；

（6）根据快速调整段制导参数确定快速调整段当前制导指令；

（7）利用探测器当前状态，规划接近段入口制导指令；

（8）快速调整段当前制导指令与接近段入口制导指令夹角小于设定值或制导时间完成，切换到接近段制导。

所述步骤（2）具体实现如下：

（1）定义主减速段末端加速度为

接近段初始加速度为

快速调整段初始加速度为

快速调整段时间为Δt_p，其受姿态机动最大角速度约束，最大姿态角速度可以取为5°/s，制导参数计算公式如下

$\begin{array}{c}{\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{Ff}}={\stackrel{\‾}{\mathrm{\λ}}}_v+(t_{\mathrm{go}}-K)\·\stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}};m_f=m_0-\frac{F\·t_{\mathrm{go}}}{I_{\mathrm{sp}}};{\stackrel{\‾}{a}}_{b,f}=\frac{F}{m_f}\·\frac{{\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{Ff}}}{\left|\right|{\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{Ff}}\left|\right|}\\ a_{a,i}=\left|\right|{\stackrel{\‾}{a}}_{a,i}\left|\right|;a_{b,f}=\left|\right|{\stackrel{\‾}{a}}_{b,f}\left|\right|\\ {\stackrel{\·}{a}}_p=(a_{a,i}-a_{b,f})/{\mathrm{\Δt}}_p\\ {\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{pv}}=({\stackrel{\‾}{a}}_{a,i}+{\stackrel{\‾}{a}}_{b,f})/2\\ {\stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}}}_p=({\stackrel{\‾}{a}}_{a,i}-{\stackrel{\‾}{a}}_{b,f})/\mathrm{\Δ}{t}_p\\ t_{\mathrm{\λ}}={\mathrm{\Δt}}_p/2\\ {\stackrel{\‾}{\mathrm{\λ}}}_p={\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{pv}}+(t-t_{\mathrm{\λ}}){\stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}}}_p\\ a_p=a_{b,f}+{\stackrel{\·}{a}}_p\·t\end{array}---\left(1\right)$

式中，

为主减速段末端推力方向，

为主减速段制导参数—参考速度增量方向，t_go为主减速段制导参数—主减速剩余时间，K为主减速段制导参数—参考时间，

为主减速段制导参数—推力方向变化率，m₀探测器当前质量，m_f为主减速段末端探测器质量，F和I_sp分别为主减速段主发动机推力和比冲，a_a,i为接近段初始加速度大小，a_b,f为主减速段末端加速度大小，为加速度变化率；

为快速调整段推力指令方向；

为快速调整段参考推力方向；t_λ为参考时间；

为推力方向的变化率；t为制导时间；a_p为快速调整段指令推力。

（2）由式（1）得到推力加速度方向和大小，于是可以预测快速调整段推力所产生的位置和速度变化量（

）；基于主减速末端状态，预测快速调整段月球引力加速度所产生的位置和速度变化量（

），组合这两项变化量，得到快速调整段共能产生的位置和速度变化量分别为

$\begin{array}{c}{\stackrel{\‾}{r}}_{\mathrm{pbias}}={\stackrel{\‾}{r}}_{\mathrm{thrust}}+{\stackrel{\‾}{r}}_{\mathrm{grav}}\\ {\stackrel{\‾}{v}}_{\mathrm{pbias}}={\stackrel{\‾}{v}}_{\mathrm{thrust}}+{\stackrel{\‾}{v}}_{\mathrm{grav}}\end{array}---\left(2\right)$

（3）定义主减速末端位置和速度分别为

和

则预测进入接近段的初始位置和速度，分别为

$\begin{array}{c}{\stackrel{\widehat{\‾}}{r}}_{\mathrm{ai}}={\stackrel{\‾}{r}}_{\mathrm{bf}}+{\mathrm{\Δt}}_p\·{\stackrel{\‾}{v}}_{\mathrm{bf}}+{\stackrel{\‾}{r}}_{\mathrm{pbias}}\\ {\stackrel{\widehat{\‾}}{v}}_{\mathrm{ai}}={\stackrel{\‾}{v}}_{\mathrm{bf}}+{\stackrel{\‾}{v}}_{\mathrm{pbias}}\end{array}---\left(3\right).$

所述步骤（3）具体实现步骤如下：

（1）定义接近段初始状态目标高度为h_ai，其由地面根据任务需要设定；利用（2）预测的接近段初始位置

计算预测的接近段初始高度

r_m为天体参考半径；并计算需要反馈的高度差

（2）定义接近段初始状态目标速度为

利用（2）预测的接近段初始速度

计算需要反馈的速度差

（3）利用反馈的高度差、速度差确定修正的主减速制导目标，修正的主减速制导位置目标为

式右侧的

为上一周期修正的结果，初值为地面预先设定主减速制导位置目标；速度目标修正为

C_h和C_v分别为高度和速度修正系数，由地面根据状态间关系确定，一般可取C_h=0.01，C_v=0.01；

（4）通过多个周期的多次修正，达到高度和速度差收敛到零的要求。

所述步骤（5）具体实现步骤如下：

（1）根据接近段初始状态要求以及主减速末端加速度，计算需要的速度增量方向矢量

制导参考时间t_λ

$\begin{array}{c}{\stackrel{\‾}{\mathrm{\λ}}}_v=({\stackrel{\‾}{a}}_{a,i}+{\stackrel{\‾}{a}}_{b,f})/2\\ t_{\mathrm{\λ}}={\mathrm{\Δt}}_p/2\end{array}---\left(4\right)$

其中，Δt_p为快速调整段时间，其受姿态机动最大角速度约束，由

和的夹角确定；

（2）确定快速调整指令加速度大小变化率

指令推力矢量变化率

$\begin{array}{c}{\stackrel{\·}{a}}_p=(a_{a,i}-a_{b,f})/{\mathrm{\Δt}}_p\\ \stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}}=({\stackrel{\‾}{a}}_{a,i}-{\stackrel{\‾}{a}}_{b,f})/{\mathrm{\Δt}}_p\end{array}---\left(5\right).$

本发明与现有技术相比的优点在于：

（1）本发明设计了主减速段带有制导目标修正功能的制导律和基于切换状态确定制导参数的自适应制导律，解决了主减速段末端状态和接近段初始状态衔接问题，保证了探测器姿态、高度、速度和加速度从主减速段过渡到接近段的平缓性，提高了深空天体探测器软着陆的安全性。

（2）本发明属于深空探测器的导航、制导与控制技术领域，非常适合深空天体探测器的高安全着陆制导的需要，可以应用于着陆或附着深空天体探测器的导航、制导与控制，也可以应用到返回地球任务的导航与制导，具有广泛的应用价值和市场前景。

附图说明

图1一种衔接主减速和接近段的预测校正方法流程图。

具体实施方式

如图1所示，本发明具体实现如下：

（1）计算主减速段制导参数

利用探测器当前状态和主发动机比冲和推力等参数，计算制导参数。主减速制导律采用PEG制导律（详见文章Space Shuttle Ascent Guidance，Navigation,and Control.The Journal of the Astronautical Science,Vol.XXVII,No.1,pp1-38,January-March,1979），具体计算过程简单描述如下。

1）确定参考速度增量方向其中，

和

分别为主减速制导速度目标和探测器当前速度；

2）定义制导参数

m₀为当前探测器质量，I_sp和F分别为发动机比冲和推力，计算剩余主减速时间为

3）计算积分参数

S=L(t_go-τ)+I_spt_go，J=Lt_go-S，

4）计算参考时间

预测推力产生的位移

其中，

为探测器当前位置，主减速制导位置目标，

为预测重力引起的位移，

为需要补偿的偏差；

5）计算推力方向变化率

为

除去航向位移部分后剩余部分；

6）计算推力产生的速度增量和位移：

和 ${\stackrel{\‾}{r}}_{\mathrm{thrust}}=S{\stackrel{\‾}{\mathrm{\λ}}}_v+(Q-\mathrm{KS})\stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}};$

7）预测主减速制导目标：速度 ${\stackrel{\‾}{v}}_{\mathrm{pd}}=\stackrel{\‾}{v}+{\stackrel{\‾}{v}}_{\mathrm{thrust}}+{\stackrel{\‾}{v}}_{\mathrm{grav}}$ 和 ${\stackrel{\‾}{r}}_{\mathrm{pd}}=\stackrel{\‾}{r}+\stackrel{\‾}{v}{t}_{\mathrm{go}}+{\stackrel{\‾}{r}}_{\mathrm{thrust}}+{\stackrel{\‾}{r}}_{\mathrm{grav}},$ 其中，

和

为天体引力产生的速度增量和位移；

8）计算制导速度偏差如果条件

利用确定的制导参数计算推力方向

其中，

为表达制导推力方向，t为制导时间；否则，修正需要推力速度增量

（ρ_g为阻尼因数，通常ρ_g=1），从1）开始新的计算，直至满足上述条件。

（2）基于制导参数和状态的快速调整终端预测

利用步骤（1）中计算制导参数和探测器当前状态，预测主减速终端状态和制导指令；根据接近段入口对姿态和加速度的要求，确定快速调整段终端状态；由此可以确定快速调整段制导参数，进而利用快速调整段制导参数预测快速调整过程产生的位置和速度变化量；最终预测得到快速调整终端状态（即预测得到接近段入口状态）。

在主减速段，采用预测/修正的制导策略。定义主减速段末端加速度为接近段初始加速度为

快速调整段初始加速度为

快速调整段时间为Δt_p，其受姿态机动最大角速度约束，这里最大姿态角速度取为5°/s。制导参数计算公式如下

$\begin{array}{c}{\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{Ff}}={\stackrel{\‾}{\mathrm{\λ}}}_v+(t_{\mathrm{go}}-K)\·\stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}}\\ m_f=m_0-\frac{F\·t_{\mathrm{go}}}{I_{\mathrm{sp}}}\\ {\stackrel{\‾}{a}}_{b,f}=\frac{F}{m_f}\·\frac{{\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{Ff}}}{\left|\right|{\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{Ff}}\left|\right|}\\ a_{a,i}=\left|\right|{\stackrel{\‾}{a}}_{a,i}\left|\right|\\ a_{b,f}=\left|\right|{\stackrel{\‾}{a}}_{b,f}\left|\right|\\ {\stackrel{\·}{a}}_p=(a_{a,i}-a_{b,f})/{\mathrm{\Δt}}_p\\ {\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{pv}}=({\stackrel{\‾}{a}}_{a,i}+{\stackrel{\‾}{a}}_{b,f})/2\\ {\stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}}}_p=({\stackrel{\‾}{a}}_{a,i}-{\stackrel{\‾}{a}}_{b,f})/{\mathrm{\Δt}}_p\\ t_{\mathrm{\λ}}={\mathrm{\Δt}}_p/2\\ {\stackrel{\‾}{\mathrm{\λ}}}_p={\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{pv}}+(t-t_{\mathrm{\λ}}){\stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}}}_p\\ a_p=a_{b,f}+{\stackrel{\·}{a}}_p\·t\end{array}---\left(1\right)$

式中，为主减速段末端推力方向，

为主减速段制导参数—参考速度增量方向，t_go为主减速段制导参数—主减速剩余时间，K为主减速段制导参数—参考时间，

为主减速段制导参数—推力方向变化率，m₀探测器当前质量，m_f为主减速段末端探测器质量，F和I_sp分别为主减速段主发动机推力和比冲，a_a,i为接近段初始加速度大小，a_b,f为主减速段末端加速度大小，为加速度变化率；

为快速调整段推力指令方向；

为快速调整段参考推力方向；t_λ为参考时间；

为推力方向的变化率；t为制导时间；a_p为快速调整段指令推力。

由式（1）可得到推力加速度方向和大小，于是，可以预测快速调整段推力所产生的位置和速度变化量（

）。基于主减速末端状态，可以预测快速调整段月球引力加速度所产生的位置和速度变化量（

）。组合这两项变化量，可以得到快速调整段共能产生的位置和速度变化量分别为

$\begin{array}{c}{\stackrel{\‾}{r}}_{\mathrm{pbias}}={\stackrel{\‾}{r}}_{\mathrm{thrust}}+{\stackrel{\‾}{r}}_{\mathrm{grav}}\\ {\stackrel{\‾}{v}}_{\mathrm{pbias}}={\stackrel{\‾}{v}}_{\mathrm{thrust}}+{\stackrel{\‾}{v}}_{\mathrm{grav}}\end{array}---\left(2\right)$

定义主减速末端位置和速度分别为

和

则可预测进入接近段的初始位置和速度，分别为

$\begin{array}{c}{\stackrel{\widehat{\‾}}{r}}_{\mathrm{ai}}={\stackrel{\‾}{r}}_{\mathrm{bf}}+{\mathrm{\Δt}}_p\·{\stackrel{\‾}{v}}_{\mathrm{bf}}+{\stackrel{\‾}{r}}_{\mathrm{pbias}}\\ {\stackrel{\widehat{\‾}}{v}}_{\mathrm{ai}}={\stackrel{\‾}{v}}_{\mathrm{bf}}+{\stackrel{\‾}{v}}_{\mathrm{pbias}}\end{array}---\left(3\right)$

（3）基于接近段入口条件修正主减速段制导目标

利用预测的接近段初始状态与接近段入口目标状态（接近段入口对高度和速度的要求）进行比较，得到接近段入口状态的差；再利用状态转移矩阵递推到主减速末端，得到主减速段制导参数的修正量，进而达到反馈修正的目的。对于主减速末期，需要经过多次预测/修正来保证接近段入口目标状态的实现。具体实现过程如下。

①定义接近段初始状态目标高度为h_ai，其由地面根据任务需要设定；利用（2）预测的接近段初始位置

计算预测的接近段初始高度

r_m为天体参考半径；并计算需要反馈的高度差

②定义接近段初始状态目标速度为

利用（2）预测的接近段初始速度

计算需要反馈的速度差

③利用反馈的高度差、速度差确定修正的主减速制导目标，修正的主减速制导位置目标为

式右侧的

上一周期修正的结果，初值为地面预先设定主减速制导位置目标；速度目标修正为

C_h和C_v分别为高度和速度修正系数，由地面根据状态间关系确定，一般可取C_h=0.01，C_v=0.01。

④通过多个周期的多次修正，达到高度和速度差收敛到零的要求。

（4）判断是否满足主减速切换条件（如制导时间小于设定值，例如0.128s），如果满足，则切换到快速调整制导；否则重新计算步骤（1）～（3），直到满足主减速切换条件。

（5）切换到快速调整制导模式后，利用切换状态确定快速调整段制导参数

1）根据接近段初始状态要求以及主减速末端加速度，计算需要的速度增量方向矢量

制导参考时间t_λ

$\begin{array}{c}{\stackrel{\‾}{\mathrm{\λ}}}_v=({\stackrel{\‾}{a}}_{a,i}+{\stackrel{\‾}{a}}_{b,f})/2\\ t_{\mathrm{\λ}}={\mathrm{\Δt}}_p/2\end{array}---\left(4\right)$

其中，Δt_p为快速调整段时间，其受姿态机动最大角速度约束，由

和

的夹角确定。

2）确定快速调整指令加速度大小变化率

指令推力矢量变化率

$\begin{array}{c}{\stackrel{\·}{a}}_p=(a_{a,i}-a_{b,f})/{\mathrm{\Δt}}_p\\ \stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}}=({\stackrel{\‾}{a}}_{a,i}-{\stackrel{\‾}{a}}_{b,f})/{\mathrm{\Δt}}_p\end{array}---\left(5\right)$

（6）根据快速调整段制导参数确定快速调整段当前制导指令

取快速调整段推力加速度a_pi初值为a_b,f，其中，a_b,f为主减速末端推力加速度大小，于是

$\begin{array}{c}\stackrel{\‾}{\mathrm{\λ}}={\stackrel{\‾}{\mathrm{\λ}}}_v+(t-t_{\mathrm{\λ}})\stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}}\\ a_{\mathrm{pi}}=a_{\mathrm{pi}}+{\stackrel{\·}{a}}_p\·{\mathrm{\Δt}}_{\mathrm{guid}}\\ {\stackrel{\‾}{a}}_{\mathrm{IC}}=a_{\mathrm{pi}}\·\frac{\stackrel{\‾}{\mathrm{\λ}}}{\left|\right|\stackrel{\‾}{\mathrm{\λ}}\left|\right|}\end{array}---\left(6\right)$

式中，t为快速调整段当前制导时间，Δt_guid为快速调整段制导周期。

（7）切换到快速调整制导模式后，利用当前状态规划出接近段入口制导指令

利用当前状态规划制导目标，调用四次多项式制导律或其它制导律确定接近段入口制导加速度指令，四次多项式制导律的具体实现见Apollolunar-descent guidance（JPL report R-695）。

（8）快速调整段当前制导指令与接近段入口制导指令夹角小于设定值或制导时间Δt_p完成，切换到接近段制导模式。指令夹角设定值一般根据设定的最大姿态角速度确定，这里取为3°。

本发明未详细阐述部分属于本领域技术人员的公知技术。

Claims (4)

1.一种衔接主减速和接近段的预测校正方法，其特征在于实现步骤如下：

（1）计算主减速段制导参数

利用探测器当前状态和主发动机比冲和推力参数，计算主减速段制导参数：剩余主减速时间、参考时间、推力方向变化率和参考速度增量方向；

（2）基于主减速段制导参数和探测器当前状态进行快速调整段终端预测

利用步骤（1）中计算主减速段制导参数和探测器当前状态，预测主减速终端状态和制导指令；根据接近段入口对姿态和加速度的要求，确定快速调整段终端状态；由此确定快速调整段制导参数，进而利用快速调整段制导参数预测快速调整过程产生的位置和速度变化量，最终预测得到快速调整终端状态，即预测得到接近段入口状态；

（3）基于接近段入口条件修正主减速段制导目标

利用预测的接近段初始状态与接近段入口目标状态进行比较，得到接近段入口状态的差，所述接近段入口目标状态即接近段入口对高度和速度的要求；再利用状态转移矩阵递推到主减速末端，得到主减速段制导参数的修正量，进而达到反馈修正的目的；经过主减速末期的多次预测修正后，可以保证接近段入口目标状态的实现；

（4）判断是否满足主减速切换条件，如果满足，则切换到快速调整制导；否则重新计算步骤（1）～（3），直到满足主减速切换条件；

（5）利用切换时刻的探测器推力加速度大小、方向和接近段入口目标状态确定快速调整段制导参数；

（6）根据快速调整段制导参数确定快速调整段当前制导指令；

（7）利用探测器当前状态，规划接近段入口制导指令；

（8）快速调整段当前制导指令与接近段入口制导指令夹角小于设定值或制导时间完成，切换到接近段制导。

2.根据权利要求1所述的一种衔接主减速和接近段的预测校正方法，其特征在于：所述步骤（2）具体实现如下：

（1）定义主减速段末端加速度为接近段初始加速度为

快速调整段初始加速度为

快速调整段时间为Δt_p，其受姿态机动最大角速度约束，制导参数计算公式如下

$\begin{array}{c}{\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{Ff}}={\stackrel{\‾}{\mathrm{\λ}}}_v+(t_{\mathrm{go}}-K)\·\stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}};m_f=m_0-\frac{F\·t_{\mathrm{go}}}{I_{\mathrm{sp}}};{\stackrel{\‾}{a}}_{b,f}=\frac{F}{m_f}\·\frac{{\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{Ff}}}{\left|\right|{\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{Ff}}\left|\right|}\\ a_{a,i}=\left|\right|{\stackrel{\‾}{a}}_{a,i}\left|\right|;a_{b,f}=\left|\right|{\stackrel{\‾}{a}}_{b,f}\left|\right|\\ {\stackrel{\·}{a}}_p=(a_{a,i}-a_{b,f})/{\mathrm{\Δt}}_p\\ {\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{pv}}=({\stackrel{\‾}{a}}_{a,i}+{\stackrel{\‾}{a}}_{b,f})/2\\ {\stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}}}_p=({\stackrel{\‾}{a}}_{a,i}-{\stackrel{\‾}{a}}_{b,f})/\mathrm{\Δ}{t}_p\\ t_{\mathrm{\λ}}={\mathrm{\Δt}}_p/2\\ {\stackrel{\‾}{\mathrm{\λ}}}_p={\stackrel{\‾}{\mathrm{\λ}}}_{\mathrm{pv}}+(t-t_{\mathrm{\λ}}){\stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}}}_p\\ a_p=a_{b,f}+{\stackrel{\·}{a}}_p\·t\end{array}---\left(1\right)$

式中，

为主减速段末端推力方向，为主减速段制导参数—参考速度增量方向，t_go为主减速段制导参数—主减速剩余时间，K为主减速段制导参数—参考时间，为主减速段制导参数—推力方向变化率，m₀探测器当前质量，m_f为主减速段末端探测器质量，F和I_sp分别为主减速段主发动机推力和比冲，a_a,i为接近段初始加速度大小，a_b,f为主减速段末端加速度大小，

为加速度变化率；

为快速调整段推力指令方向；

为快速调整段参考推力方向；t_λ为参考时间；

为推力方向的变化率；t为制导时间；a_p为快速调整段指令推力;

（2）由式（1）得到推力加速度方向和大小，于是可以预测快速调整段推力所产生的位置和速度变化量（

）；基于主减速末端状态，预测快速调整段月球引力加速度所产生的位置和速度变化量（

），组合这两项变化量，得到快速调整段共能产生的位置和速度变化量分别为

$\begin{array}{c}{\stackrel{\‾}{r}}_{\mathrm{pbias}}={\stackrel{\‾}{r}}_{\mathrm{thrust}}+{\stackrel{\‾}{r}}_{\mathrm{grav}}\\ {\stackrel{\‾}{v}}_{\mathrm{pbias}}={\stackrel{\‾}{v}}_{\mathrm{thrust}}+{\stackrel{\‾}{v}}_{\mathrm{grav}}\end{array}---\left(2\right)$

（3）定义主减速末端位置和速度分别为和

则预测进入接近段的初始位置和速度，分别为

$\begin{array}{c}{\stackrel{\widehat{\‾}}{r}}_{\mathrm{ai}}={\stackrel{\‾}{r}}_{\mathrm{bf}}+{\mathrm{\Δt}}_p\·{\stackrel{\‾}{v}}_{\mathrm{bf}}+{\stackrel{\‾}{r}}_{\mathrm{pbias}}\\ {\stackrel{\widehat{\‾}}{v}}_{\mathrm{ai}}={\stackrel{\‾}{v}}_{\mathrm{bf}}+{\stackrel{\‾}{v}}_{\mathrm{pbias}}\end{array}---\left(3\right).$

3.根据权利要求1所述的一种衔接主减速和接近段的预测校正方法，其特征在于：所述步骤（3）具体实现步骤如下：

（1）定义接近段初始状态目标高度为h_ai，其由地面根据任务需要设定；利用（2）预测的接近段初始位置

计算预测的接近段初始高度

为天体参考半径；并计算需要反馈的高度差

（2）定义接近段初始状态目标速度为

利用（2）预测的接近段初始速度

计算需要反馈的速度差

（3）利用反馈的高度差、速度差确定修正的主减速制导目标，修正的主减速制导位置目标为式右侧的为上一周期修正的结果，初值为地面预先设定主减速制导位置目标；速度目标修正为

C_h和C_v分别为高度和速度修正系数，由地面根据状态间关系确定；

（4）通过多个周期的多次修正，达到高度和速度差收敛到零的要求。

4.根据权利要求1所述的一种衔接主减速和接近段的预测校正方法，其特征在于：所述步骤（5）具体实现步骤如下：

（1）根据接近段初始状态要求以及主减速末端加速度，计算需要的速度增量方向矢量

制导参考时间t_λ

$\begin{array}{c}{\stackrel{\‾}{\mathrm{\λ}}}_v=({\stackrel{\‾}{a}}_{a,i}+{\stackrel{\‾}{a}}_{b,f})/2\\ t_{\mathrm{\λ}}={\mathrm{\Δt}}_p/2\end{array}---\left(4\right)$

其中，Δt_p为快速调整段时间，其受姿态机动最大角速度约束，由和

的夹角确定；

（2）确定快速调整指令加速度大小变化率指令推力矢量变化率

$\begin{array}{c}{\stackrel{\·}{a}}_p=(a_{a,i}-a_{b,f})/{\mathrm{\Δt}}_p\\ \stackrel{\·}{\stackrel{\‾}{\mathrm{\λ}}}=({\stackrel{\‾}{a}}_{a,i}-{\stackrel{\‾}{a}}_{b,f})/{\mathrm{\Δt}}_p\end{array}---\left(5\right).$

Images