一种基于牛顿迭代的Lambert变轨问题半长轴迭代空间转移轨道计算方法技术

本发明专利技术公开了一种基于牛顿迭代思想的Lambert变轨问题改进求解方法。步骤1：计算转移角θ、辅助变量c和辅助变量s；步骤2：根据步骤1的转移角和辅助变量，确定轨道形状；步骤3：当Δt＞Δt

【技术实现步骤摘要】
一种基于牛顿迭代的Lambert变轨问题半长轴迭代空间转移轨道计算方法
本专利技术属于航天器轨道设计的
；具体涉及一种基于牛顿迭代的Lambert变轨问题半长轴迭代空间转移轨道计算方法。
技术介绍
Lambert变轨问题，也被称为两点边界值问题，是指航天器在指定转移时间飞过两个预定位置矢量的轨道确定。该方法求解空间转移轨道不受轨道高度，轨道倾角等条件的影响，可求解任意两个位置矢量间的圆锥曲线轨道，具有广泛适用性。由Lagrange转移时间方程可以知道，转移时间和转移轨道的轨道要素之间关系复杂，无法通过数值方法直接确定满足要求的转移轨道曲线。目前，Lambert变轨问题求解的基本思路为选择迭代变量对转移时间求解进行转化，确定其初值，最后通过对迭代结果进行转换，得到转移轨道起始点速度和交汇点的速度，进而确定转移轨道的轨道要素。Gauss最早提出了利用迭代变量求解Lambert变轨问题的方法，Battin-Vaughan进一步完善Gauss算法，通过对辅助变量x的迭代计算得到转移轨道半长轴、半通径、离心率等轨道根数，并得到了转移轨道上起始时刻和交会时刻的速度矢量,该方法只对于360°的等径转移存在奇异，但是计算繁琐，耗时很长。Klumpp进一步改进Battin算法，消除了奇异性。但该方法仍需求解超几何函数，计算复杂，时效性差。
技术实现思路
本专利技术的目的是为了提高航天器交会中的Lambert变轨问题的求解速度，提出了基于牛顿迭代思想的一种Lambert变轨问题的半长轴迭代求解方法。r>本专利技术通过以下技术方案实现：一种基于牛顿迭代思想的Lambert变轨问题的半长轴迭代求解方法，所述半长轴迭代空间转移轨道计算方法包括以下步骤：步骤1：已知初始时刻t1任务飞行器在起始轨道上的位置矢量为r1，交会时刻t2目标飞行器在目标轨道的位置矢量为r2，计算转移角θ、辅助变量c和辅助变量s；步骤2：根据步骤1计算的转移角θ、辅助变量c和辅助变量s，确定飞行器的空间转移轨道形状；步骤3：当转移时间Δt大于抛物线轨道转移时间Δtp时，则转移轨道为椭圆轨道，对椭圆轨道求解；引入最小能量椭圆半长轴当a＜am时，两点间不存在椭圆转移轨道，即对于椭圆轨道有am＜a＜+∞,0≤e＜1；步骤4：当转移时间Δt小于抛物线轨道转移时间Δtp时，则转移轨道为双曲线轨道，对双曲线轨道求解；对于双曲线轨道有-∞＜a＜0,e＞1；步骤5：当转移时间Δt等于抛物线轨道转移时间Δtp时，则转移轨道为抛物线轨道，对抛物线轨道求解；对于抛物线轨道有半长轴az＝∞,e＝0；步骤6：将步骤3-5中的将Lambert变轨问题解析方程表示成圆锥曲线半长轴的函数，计算圆锥曲线半通径；步骤7：根据步骤3至步骤6知为初始时刻t1任务飞行器在起始轨道上的速度矢量，为交会时刻t2目标飞行器在目标轨道的速度矢量，计算任务飞行器和目标飞行器变轨点速度和交汇点速度，从而完成空间转移轨道计算过程。进一步的，所述步骤1中计算转移角θ、辅助变量c和辅助变量s，具体为，转移角辅助变量辅助变量式中，r1为初始时刻t1任务飞行器在起始轨道上的位置矢量，r2为交会时刻t2目标飞行器在目标轨道的位置矢量，r1＝||r1||,r2＝||r2||；设Δt为转移轨道的转移时间，Δt＝t2-t1。进一步的，所述步骤2中根据步骤1计算的转移角θ、辅助变量c和辅助变量s，确定轨道形状；其中μ为地球引力常数，Δtp为抛物线轨道转移时间。进一步的，所述步骤3包括以下步骤：步骤3.1：利用固定大步长迭代法得出搜索闭区间：令搜索区间上限a0＝am,搜索区间下限a1＝10am，由和得到辅助变量α,β，当θ≤π时，取0≤α≤2π,0≤β≤π；当θ＞π时，取0≤α≤2π,-π≤β≤0；步骤3.2：取并将Lambert变轨问题解析方程表示成圆锥曲线半长轴的函数，即分别计算f(a0),f(az),f(a1)；步骤3.3：如果f(a0)f(a1)＜0，则进行步骤3.4；如果f(a0)f(a1)＞0，则a0＝a1,a1＝2a1，重复步骤3.2；步骤3.4：引入牛顿迭代法，记迭代变量判断是否满足f(ai+1)→0，若不满足则循环执行步骤3.4，直到满足条件为止；步骤3.5：当f(ai+1)→0时，此时取a＝ai+1即为利用牛顿迭代法求得的椭圆转移轨道的半长轴。进一步的，所述步骤4中包括以下步骤：步骤4.1：选定固定大步长迭代法的搜索闭区间：令搜索区间上限a0＝-106，搜索区间下限a1＝-108，由和得到辅助变量α,β，当θ≤π时，取0≤α,0≤β；当θ＞π时，取0≤α,β≤0；步骤4.2：取并将Lambert变轨问题解析方程表示成圆锥曲线半长轴的函数，即分别计算f(a0),f(az),f(a1)；步骤4.3：如果f(a0)f(a1)＜0，则进行步骤4.4；如果f(a0)f(a1)＞0，则a0＝a1,a1＝2a1，重复步骤4.2；步骤4.4：引入牛顿迭代法，记迭代变量判断是否满足f(ai+1)→0，若不满足则循环执行步骤4.4，直到满足条件为止；步骤4.5：当f(ai+1)→0时，此时取a＝ai+1即为利用牛顿迭代法求得的双曲线转移轨道的半长轴。进一步的，所述步骤6中将步骤3-5中的将Lambert变轨问题解析方程表示成圆锥曲线半长轴的函数，计算圆锥曲线半通径，椭圆轨道抛物线轨道其中A＝r12-2r1r2cosΔf+r22，B＝-2r1r2(r1+r2)(1-cosΔf)，C＝r12r22(1-cosΔf)2，Δf为飞行转移角，其中双曲线轨道进一步的，所述步骤7中计算变轨点速度和交汇点速度具体为，变轨点速度矢量交会点速度矢量本专利技术的有益效果是：本专利技术解决了已知初始点位置、交会点位置和转移时间的条件下，求解空间转移轨道的问题，提高Lambert变轨问题的求解效率，具有广阔的应用前景。附图说明附图1本专利技术的总体框图。具体实施方式下面将结合本专利技术实施例中的附图对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。实施例1一种基于牛顿迭代思想的Lambert变轨问题的半长轴迭代求解方法，所述半长轴迭代求解方法包括以下步骤：步骤1：已知初始时刻t1任务飞行器在起始轨道上的位置矢量为r1，交会时刻t2目标飞行器在目标轨道的位置矢量为r2，计算转移角θ、辅助变量c和辅助变量s；步骤2：根据步骤1计算的转移角θ、辅助变量c和辅助变量s，确定飞行器的空间转移轨道形状；步骤3：当转本文档来自技高网...

【技术保护点】
1.一种基于牛顿迭代的Lambert变轨问题半长轴迭代空间转移轨道计算方法，其特征在于，所述半长轴迭代空间转移轨道计算方法包括以下步骤：/n步骤1：已知初始时刻t

【技术特征摘要】
1.一种基于牛顿迭代的Lambert变轨问题半长轴迭代空间转移轨道计算方法，其特征在于，所述半长轴迭代空间转移轨道计算方法包括以下步骤：
步骤1：已知初始时刻t1任务飞行器在起始轨道上的位置矢量为r1，交会时刻t2目标飞行器在目标轨道的位置矢量为r2，计算转移角θ、辅助变量c和辅助变量s；
步骤2：根据步骤1计算的转移角θ、辅助变量c和辅助变量s，确定飞行器的空间转移轨道形状；
步骤3：当转移时间Δt大于抛物线轨道转移时间Δtp时，则转移轨道为椭圆轨道，对椭圆轨道求解；
引入最小能量椭圆半长轴当a＜am时，两点间不存在椭圆转移轨道，即对于椭圆轨道有am＜a＜+∞,0≤e＜1；
步骤4：当转移时间Δt小于抛物线轨道转移时间Δtp时，则转移轨道为双曲线轨道，对双曲线轨道求解；
对于双曲线轨道有-∞＜a＜0,e＞1；
步骤5：当转移时间Δt等于抛物线轨道转移时间Δtp时，则转移轨道为抛物线轨道，对抛物线轨道求解；
对于抛物线轨道有半长轴az＝∞,e＝0；
步骤6：将步骤3至步骤5中的将Lambert变轨问题解析方程表示成圆锥曲线半长轴的函数，计算圆锥曲线半通径；
步骤7：根据步骤3至步骤6知为初始时刻t1任务飞行器在起始轨道上的速度矢量，为交会时刻t2目标飞行器在目标轨道的速度矢量，计算任务飞行器和目标飞行器变轨点速度和交汇点速度，从而完成空间转移轨道计算过程。


2.根据权利要求1所述半长轴迭代空间转移轨道计算方法，其特征在于，所述步骤1中计算转移角θ、辅助变量c和辅助变量s，具体为，
转移角
辅助变量
辅助变量
式中，r1为初始时刻t1任务飞行器在起始轨道上的位置矢量，r2为交会时刻t2目标飞行器在目标轨道的位置矢量，r1＝||r1||,r2＝||r2||，设Δt为转移轨道的转移时间，Δt＝t2-t1。


3.根据权利要求1所述半长轴迭代空间转移轨道计算方法，其特征在于，所述步骤2中根据步骤1计算的转移角θ、辅助变量c和辅助变量s，确定轨道形状；



其中μ为地球引力常数，Δtp为抛物线轨道转移时间。


4.根据权利要求1所述半长轴迭代空间转移轨道计算方法，其特征在于，所述步骤3包括以下步骤：
步骤3.1：利用固定大步长迭代法得出搜索闭区间：令搜索区间上限a0...

【专利技术属性】
技术研发人员：王松艳，晁涛，蒋瑞晔，杨明，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江;23