[发明专利]一种火星大气进入段最优预测制导方法

摘要

本发明公开的一种火星大气进入段最优预测制导方法，涉及一种最优制导方法，属于深空探测技术领域。本发明将最优控制的方法引入火星大气进入段的制导律设计中，通过确定末端高度最优条件下的剖面形式，最优倾侧角开关曲线的最大值|σ|_max和最小值|σ|_min取值，以及开关曲线的开关时机，能够在满足末端开伞位置精度的同时，最大限度地提升末端高度，从而为火星进入、下降和着陆任务的后续操作提供更大的时间裕度，保证火星着陆探测任务的安全实施。此外，本发明通过搜索使得性能指标取得极小值的方式来确定开关时机采用所述的性能指标能够避免在线数值求解非线性方程时可能出现的无解的情形，从而提高算法的稳定性和鲁棒性，能够保证末端位置精度。

主权项

1.一种火星大气进入段最优预测制导方法，其特征在于：包括如下步骤，步骤1、确定纵向运动在末端高度h_f最优条件下的剖面形式；利用探测器动力学模型进行数值积分至满足火星大气进入段末端条件，所述的火星大气进入段末端条件为开伞条件，得到末端航程s_f，进而得到开伞时刻的开伞剩余纵程与目标位置的开伞精度偏差；所述的开伞条件指探测器动压在区间[q_min,q_max]内和探测器马赫数在区间[Ma_min,Ma_max]内；利用探测器动力学模型进行数值积分至满足火星大气进入段末端条件具体实现方法为：考虑火星自转影响的探测器对无量纲时间的三自由度无量纲进入动力学模型为，其中，s为剩余纵程，表征从探测器当前位置到标称末端位置的火星表面大圆弧的距离，r为火星质心到探测器质心的距离，无量纲参数为火星半径R₀，v为探测器相对于火星的速度，速度v的无量纲参数为其中g₀为火星表面重力加速度，γ为航迹角，σ为倾侧角，纵向动力学只确定倾侧角σ的大小|σ|，而倾侧角σ的符号有相应的侧向逻辑确定，g为当地重力加速度，其无量纲参数为g₀；D和L分别为阻力加速度和升力加速度，其无量纲参数均为g₀，C_D和C_L分别为阻力系数和升力系数，S为探测器参考面积，m为探测器质量，q＝ρv²/2为动压，β＝m/SC_D为探测器弹道系数，L/D为探测器升阻比；火星大气密度采用如公式(3)所示的指数模型，其中ρ₀为参考密度，h为飞行器飞行高度，h₀为参考高度，h_s为大气密度标高；进入段飞行器的比能量由公式(4)给出，根据极小值原理，性能指标，在末端高度最优条件下的形式为，J＝‑h_f＝‑(r_f‑R₀)   (6)其中，Φ[x(t_f),t_f]＝‑r_f,L[x(t_f),u(t_f),t_f]＝0，t_f为末端时刻，下标“f”表示末端时刻对应的各物理量，h_f为末端时刻对应的高度，x(t_f)为末端时刻对应的状态向量，u(t_f)为末端时刻的控制变量，即倾侧角的大小|σ|；则根据动力学方程(1)，哈密顿函数H为，其中，为λ＝[λ_r,λ_v,λγ,λ_s]^T协状态变量，满足即根据极小值原理有即边界条件为横截条件为式(12)中，υ₁，υ₂为待定常数；式(10)表明，火星大气进入段过程中，末端高度最优条件下的倾侧角|σ|剖面为开关曲线，即bang‑bang控制律；步骤2、确定倾侧角开关曲线所需要满足的末端条件；当时，由公式(8)可知，的左极限为则由于在末端时刻有根据横截条件式(13)可知，υ₁＞0；由于均为v(t_f),g(t_f),cosγ(t_f)正数，则的极限由横截条件式(12)可知，λ_γ(t_f)＝0；结合式(16)可知，在t＝t_f的一个左邻域内，有λ_γ＜0   (17)则根据式(10)可知，在t＝t_f的一个左邻域内，有|σ|^*＝|σ|_min   (18)由式(18)可知，倾侧角开关曲线在经过最后一次开关后至进入段末端，倾侧角|σ|需满足取极小值|σ|_min的末端条件；步骤3、确定倾侧角开关曲线的最大值|σ|_max和最小值|σ|_min取值范围；为满足着陆位置精度要求，倾侧角开关曲线的最大值|σ|_max和最小值|σ|_min需位于常值倾侧角|σ|₀的两侧，即满足公式(19)0≤|σ|_min＜|σ|₀＜|σ|_max≤π   (19)步骤4、确定倾侧角开关曲线最大值|σ|_max和最小值|σ|_min的数值；在bang‑bang控制律下，倾侧角|σ|应在其上下边界上切换；为满足式(19)，取倾侧角开关曲线最大值|σ|_max和最小值|σ|_min分别为步骤5、确定倾侧角开关曲线的开关时机，进而保证末端位置精度；所述的开关时机指火星大气进入段轨迹的物理量ζ满足式(22)条件所对应的时间t_s所述的物理量ζ具体包括进入速度v、比能量e；对应的开关时机v_s、e_s分别称为开关速度、开关能量；通过确定开关时机能够确定倾侧角开关曲线执行开关操作的时刻，保证末端位置精度。