形状参数(Shape Parameter)
本文作者:天疆说
定义
在轨迹设计的形状基方法(shape-based method)中,形状参数是定义轨迹几何曲线的自由参数。不对应于特定的物理量,而是通过假设轨迹具有某种解析形状(如指数正弦曲线、多项式曲线等)来参数化航天器的运动路径;推力剖面、飞行时间和各段燃料消耗由形状参数通过逆动力学反推得到(Petropoulos 2002;Izzo 2006)。
最知名的形状参数体系是 Petropoulos (2002) 提出的指数正弦曲线(exponential sinusoid),由四个参数定义:
其中 为极径, 为极角。四个参数的角色(Izzo 2006):
:控制轨道总体尺寸(距离尺度),决定极大/极小极径。
:控制每个周期内极径振荡的幅度—— 越大,轨迹在径向的螺旋波动越显著。
:控制极角频率,决定每个公转周期内轨迹绕行圈数。在 Petropoulos 的原公式中, 就是形状参数(shape parameter)的狭义含义。
:控制相位,即轨迹在 时的初始相位角。
动力学可行性
假设航天器始终沿轨迹切线方向施加推力(tangential thrust),则推力剖面和极角时间历程均由上述四个参数唯一确定(Petropoulos 2002;Izzo 2006):
关键:轨迹物理可行(推力始终为正、轨迹不遭回折)当且仅当 (Izzo 2006)。
多圈指数正弦 Lambert 问题
在出发行星速度已知的条件下,四个参数退化为一个自由参数 ——即形状参数(shape parameter)。给定 (出发位置)、(目标位置)、(转移极角)和 (绕行圈数),所有可能连接两点的指数正弦曲线构成一个由出发点航迹角 参数化的单参数族 (Izzo 2006)。
该族中满足动力学可性()的轨迹,其飞行时间 是 的单调函数——这对多圈指数正弦 Lambert 问题的适定性和单值性至关重要。Lambert 问题通过数值求积和二分法则(Regula Falsi)求解,给定 和飞行时间即可找到唯一对应的指数正弦轨迹(Izzo 2006)。在行星际全局优化中, 是核心决策变量之一:例如日地到火星的 N=0 最优解对应 ,N=1 对应 ,N=2 对应 (Izzo 2006 表 1)。
与小推力转移设计的关联
形状基方法的吸引力在于:它将通常需要间接法(庞特里亚金极大值原理 + 多点边值问题)求解的复杂最优控制问题,简化为对少数形状参数的搜索。缺点在形状参数数量有限(指数正弦线最多 4 个)时不可避免——所表达的轨迹仅是全部可行轨迹的一个子集("shape-restricted optimal" 而非 "globally optimal")。然而对初步设计阶段的方案筛选、全局优化算法的初值生成、以及大量轨迹方案的快速扫描而言,这是目前最实用的折中之一(Izzo 2006;Vasile et al. 2007)。
对于涉及多项式变化的振幅和相位假设:平动点附近的小推力转移中,可假设振幅和相位按多项式随时间变化(线性或更高阶),以此拟合螺旋转移特性,但这仅是对特定轨道族的一种启发式建模——其通用性不及指数正弦曲线的形状参数(Petropoulos 2002)。
相关概念
参考文献
Petropoulos, 2002, Ph.D. dissertation, Purdue University(指数正弦曲线形状基方法的原始提出)
Izzo, 2006, Lambert's problem for exponential sinusoids, JGCD(多圈指数正弦 Lambert 问题的系统化求解,形状参数 的角色与最优值的确定)
Vasile, Schütze et al., 2007(用更多自由参数的形状基方法扩展指数正弦方法,权衡曲线保真度与搜索维度)
