CN110489791A - 一种脉冲速度增量改变轨道根数的高阶解析计算方法 - Google Patents

Abstract

本发明提供一种脉冲速度增量改变轨道根数的高阶解析计算方法，包括以下步骤：1，输入脉冲速度增量改变轨道根数的计算方法所需要达到的精度阶数N；2，采用微分形式的高斯变分方程通过链式法则依次得到轨道根数的1～N阶导数；3，采用泰勒展开法得到轨道根数变化关于脉冲速度增量的N阶解析函数；4，利用所建立的N阶解析函数计算给定脉冲速度增量下的轨道根数改变量，完成高阶解析计算。本发明的计算方法是完全解析的，因而计算效率高；且能够达到任意阶次，因而计算精度高。这在航天器轨道控制方面具有广泛的使用价值。

Description

一种脉冲速度增量改变轨道根数的高阶解析计算方法

技术领域

本发明涉及航空航天技术领域，尤其涉及一种脉冲速度增量改变轨道根数的高阶解析计算方法。

背景技术

航天器的轨道控制通常是通过推力器瞬间产生的等效脉冲速度增量实现的。在任意给定的脉冲速度增量下，航天器的轨道根数变化量可通过高斯变分方程来评估。现在被广泛使用的积分型高斯变分方程(IGVEs)是在微分型高斯变分方程(DGVEs)基础上通过一阶近似和脉冲速度增量等效替换得到的，因此也称作一阶IGVEs。值得指出的是，尽管DGVEs是精确的，但它只能适用于连续力控制模式，而不能适用于脉冲速度增量控制模式，因而在工程上使用较少。但能够适用于脉冲速度增量控制模式的一阶IGVEs的精度不高，只能用于评估短时间或较小脉冲速度增量控制下的轨道改变，这大大限制了工程实践的应用。为此，有学者提出了构造二阶IGVEs方程的方法。该方法的实质是对轨道积分与位置速度矢量之间的解析函数进行二阶泰勒级数展开。相比一阶 IGVEs，二阶IGVEs能够以更高精度评估脉冲速度增量产生的轨道根数变化。但用于求出该二阶IGVEs的方法具有明显的缺点，即轨道积分与位置速度矢量之间的解析函数是非线性的、隐形式的，非常不便于进行高阶解的进一步求解。这也是上述方法仅得到了二阶IGVEs而没有得到三阶及更高阶IGVEs的原因。考虑到航天器轨道精确控制对更高阶IGVEs的需求，有必要采用一种简易便捷且通用的方法，从而构造出解析的任意阶IGVEs。

发明内容

为解决现有技术中存在的问题，本发明的目的在于提供一种脉冲速度增量改变轨道根数的高阶解析计算方法，本发明能够快速精确地确定给定脉冲速度增量下航天器轨道根数的变化。

为实现上述目的，本发明采用以下技术手段：

一种脉冲速度增量改变轨道根数的高阶解析计算方法，包括以下步骤：

S1，获取用于计算脉冲速度增量改变轨道根数的解析解的阶数N；其中N 为任意正整数；

S2，根据微分形式的高斯变分方程通过链式法则依次得到轨道根数的1～N 阶导数；

S3，根据泰勒展开法得到轨道根数变化关于脉冲速度增量的N阶解析函数；

S4，根据建立的解析函数计算给定脉冲速度增量下的轨道根数改变量，完成高阶解析计算。

作为本发明的进一步改进，步骤S2具体为：

S2.1，给出轨道根数的微分形式的高斯变分方程；

S2.2，采用微积分的链式法则，对上述微分形式的高斯变分方程进一步求导。

作为本发明的进一步改进，步骤S2具体为：

S2.1，给出微分形式的高斯变分方程如下：

其中，等式左边的t为时间，a、e、i、Ω、ω、f为轨道根数，p＝a(1-e²)、为与轨道根数相关的两个轨道参数，μ为地球引力系数，F_r、F_t、F_n为连续控制力；具体地：

a、轨道半长轴；e、轨道偏心率；i、轨道倾角；Ω、升交点赤经；ω、纬度幅角；f、真近点角；p、轨道通径；h、轨道角动量；F_r、径向控制力；F_t、切向控制力；F_n、法向控制力；

S2.2，采用微积分的链式法则，对上述微分形式的高斯变分方程进一步求导：对任意六个轨道根数中的一个α_k，若已得到其1～m-1阶的导数，则按照下述公式求解其第m阶的导数：

其中，的表达式采用微分型高斯变分方程；

当m＝N时，终止上述求导过程，否则继续。

作为本发明的进一步改进，步骤S3具体为：

S3.1，采用泰勒展开法得到轨道根数变化关于脉冲速度增量的N阶形式表达式；

S3.2，将表达式中的连续力和时间脉冲按照下述方式结合形成脉冲速度增量；

S3.3，将N阶形式表达式的连续力表达式变换为脉冲速度增量表达式，完成了N阶IGVEs的构建。

作为本发明的进一步改进，步骤S3具体为：

S3.1，采用泰勒展开法得到轨道根数变化关于脉冲速度增量的N阶形式表达式：

S3.2，将表达式中的连续力和时间脉冲按照下述方式结合形成脉冲速度增量：

Δv_r＝F_r·Δt,Δv_t＝F_t·Δt,Δv_n＝F_n·Δt (4)

S3.3，考虑公式(4)中的等价关系，将公式(3)中的连续力表达式变换为脉冲速度增量表达式，从而得到：

其中，η_kpq是经整理后得到的系数；完成了N阶IGVEs的构建。

与现有技术相比，本发明的有益效果是：

本发明先输入脉冲速度增量改变轨道根数的计算方法所需要达到的精度阶数，再采用微分形式的高斯变分方程通过链式法则依次得到1～N阶轨道根数导数；然后采用泰勒展开法得到轨道根数变化关于脉冲速度增量的N阶解析函数；最后利用所建立的解析函数计算给定脉冲速度增量下的轨道根数改变量，完成高阶解析计算。该方法能够逐次计算得到任意高阶的积分型高斯变分方程，因而可根据工程需要提高问题求解的精度。

优选地，所得结果具有解析表达式，因而计算代价小、计算效率高，适合在轨实时应用。

附图说明

图1是本发明的高阶解析计算方法的计算流程示意图；

图2是1阶、2阶、3阶IGVEs方法在10000次蒙特卡洛仿真中轨道半长轴改变量的计算误差统计结果。

图3是1阶、2阶、3阶IGVEs方法在10000次蒙特卡洛仿真中轨道偏心率改变量的计算误差统计结果。

图4是1阶、2阶、3阶IGVEs方法在10000次蒙特卡洛仿真中轨道倾角改变量的计算误差统计结果。

图5是1阶、2阶、3阶IGVEs方法在10000次蒙特卡洛仿真中升交点赤经改变量的计算误差统计结果。

图6是1阶、2阶、3阶IGVEs方法在10000次蒙特卡洛仿真中纬度幅角改变量的计算误差统计结果。

图7是1阶、2阶、3阶IGVEs方法在10000次蒙特卡洛仿真中真近点角改变量的计算误差统计结果。

具体实施方式

下面结合附图对本发明的结构和工作原理作进一步详细说明。

如图1所示，本发明提出的一种脉冲速度增量改变轨道根数的高阶解析计算方法，包括以下步骤：

S1，输入脉冲速度增量改变轨道根数的计算方法所需要达到的精度阶数N；其中N可以为1、2、3等任意正整数；

S2，采用微分形式的高斯变分方程通过链式法则依次得到轨道根数的1～N 阶导数；

S3，采用泰勒展开法得到轨道根数变化关于脉冲速度增量的N阶解析函数；

S4，利用所建立的解析函数计算给定脉冲速度增量下的轨道根数改变量，完成高阶解析计算。

优选的，S2具体包括以下步骤：

S2.1，给出微分形式的高斯变分方程如下：

其中，等式左边的t为时间，a、e、i、Ω、ω、f为轨道根数，p＝a(1-e²)、为与轨道根数相关的两个轨道参数，μ为地球引力系数，F_r、F_t、F_n为连续控制力。具体地：

a——轨道半长轴；

e——轨道偏心率；

i——轨道倾角；

Ω——升交点赤经；

ω——纬度幅角；

f——真近点角；

p——轨道通径；

h——轨道角动量；

F_r——径向控制力；

F_t——切向控制力；

F_n——法向控制力；

S2.2，采用微积分的链式法则，对上述微分形式的高斯变分方程进一步求导。即，对任意六个轨道根数中的一个α_k，若已得到其1～m-1阶的导数，则按照下述公式求解其第m阶的导数：

其中，的表达式采用微分型高斯变分方程。

当m＝N时，终止上述求导过程，否则继续。

优选的，S3具体包括以下步骤：

S3.1，采用泰勒展开法得到轨道根数变化关于脉冲速度增量的N阶形式表达式：

S3.2，将表达式中的连续力和时间脉冲按照下述方式结合形成脉冲速度增量：

Δv_r＝F_r·Δt,Δv_t＝F_t·Δt,Δv_n＝F_n·Δt (4)

S3.3，考虑公式(4)中的等价关系，将公式(3)中的连续力表达式变换为脉冲速度增量表达式，从而得到：

其中，η_kpq是经整理后得到的系数。这样能就完成了N阶IGVEs的构建。

下面列举一个具体实施例和附图，说明本发明的具体计算过程。

实施例1

如图1所示，本发明的计算方法包括：

S1，输入脉冲速度增量改变轨道根数的计算方法所需要达到的精度阶数 N＝3；

S2，采用微分形式的高斯变分方程通过链式法则依次得到轨道根数的1、2、 3阶导数；具体地，

S2.1给出微分形式的高斯变分方程作为轨道根数的1阶导数：

S2.2对上述方程进行求导，采用链式法则，得到轨道根数的2阶导数：

S2.3对上述方程进行求导，利用链式法则，得到轨道根数的3阶导数：

S3，采用泰勒展开法得到轨道根数变化关于脉冲速度增量的3阶解析函数；具体地，

S3.1，采用泰勒展开法得到轨道根数变化关于脉冲速度增量的3阶形式表达式：

S3.2，将表达式中的连续力和时间脉冲按照下述方式结合形成脉冲速度增量：

Δv_r＝F_r·Δt,Δv_t＝F_t·Δt,Δv_n＝F_n·Δt (5)

S3.3，考虑公式(5)中的等价关系，将公式(4)中的连续力表达式变换为脉冲速度增量表达式，从而得到：

其中上标(1)、(2)、(3)分别表示1阶、2阶、3阶项，具体表达式为：

S4，利用所建立的解析函数计算给定脉冲速度增量下的轨道根数改变量，完成高阶解析计算。

结果验证，采用10000次蒙特卡洛仿真说明本发明方法的有效性。

假设参考轨道的六个轨道根数按照如下方式随机选取：

a∈[a_min,a_max]，其中a_min＝R_E+300km,a_max＝R_E+20000km，R_E＝6378.14km；

e∈[e_min,e_max],其中e_min＝0.01,e_max＝1-a_min/a；

i∈[0,π],Ω∈[0,2π],ω∈[0,2π],f∈[0,2π].

假设脉冲速度增量按照如下方式随机选取：

Δv_r∈[-100,100]m/s、Δv_t∈[-100,100]m/s、Δv_n∈[-100,100]m/s。

则通过采用精确的非线性的位置速度到轨道根数的函数关系计算得到轨道根数的真实变化量，记作Δα_true。按照本发明提供的1阶、2阶、3阶解，得到的相应轨道根数变化量记作Δα_est。计算两者之间的偏差，并按照偏差量级进行统计，得到结果如图2～7所示。由图可知，对六个轨道根数中的任意一个，三阶解的精度总体上高于二阶解，而二阶解的精度总体上高于一阶解。

以上实施例仅用于说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细说明，所属领域的普通技术人员依然可以对本发明的具体实施方案进行修改或者等同替换，而这些并未脱离本发明精神和范围的任何修改或者等同替换，其均在本发明的权利要求保护范围之内。

Claims (5)

1.一种脉冲速度增量改变轨道根数的高阶解析计算方法，其特征在于，包括以下步骤：

S1，获取用于计算脉冲速度增量改变轨道根数的解析解的阶数N；其中N为任意正整数；

S2，根据微分形式的高斯变分方程通过链式法则依次得到轨道根数的1～N阶导数；

S3，根据泰勒展开法得到轨道根数变化关于脉冲速度增量的N阶解析函数；

S4，根据建立的解析函数计算给定脉冲速度增量下的轨道根数改变量，完成高阶解析计算。

2.根据权利要求1所述的脉冲速度增量改变轨道根数的高阶解析计算方法，其特征在于，步骤S2具体为：

S2.1，给出轨道根数的微分形式的高斯变分方程；

S2.2，采用微积分的链式法则，对上述微分形式的高斯变分方程进一步求导。

3.根据权利要求1或2所述的脉冲速度增量改变轨道根数的高阶解析计算方法，其特征在于，步骤S2具体为：

S2.1，给出微分形式的高斯变分方程如下：

其中，等式左边的t为时间，a、e、i、Ω、ω、f为轨道根数，p＝a(1-e²)、为与轨道根数相关的两个轨道参数，μ为地球引力系数，F_r、F_t、F_n为连续控制力；具体地：

a、轨道半长轴；e、轨道偏心率；i、轨道倾角；Ω、升交点赤经；ω、纬度幅角；f、真近点角；p、轨道通径；h、轨道角动量；F_r、径向控制力；F_t、切向控制力；F_n、法向控制力；

S2.2，采用微积分的链式法则，对上述微分形式的高斯变分方程进一步求导：对任意六个轨道根数中的一个α_k，若已得到其1～m-1阶的导数，则按照下述公式求解其第m阶的导数：

其中，的表达式采用微分型高斯变分方程；

当m＝N时，终止上述求导过程，否则继续。

4.根据权利要求1或2所述的脉冲速度增量改变轨道根数的高阶解析计算方法，其特征在于，步骤S3具体为：

S3.1，采用泰勒展开法得到轨道根数变化关于脉冲速度增量的N阶形式表达式；

S3.2，将表达式中的连续力和时间脉冲按照下述方式结合形成脉冲速度增量；

S3.3，将N阶形式表达式的连续力表达式变换为脉冲速度增量表达式，完成了N阶IGVEs的构建。

5.根据权利要求3所述的脉冲速度增量改变轨道根数的高阶解析计算方法，其特征在于，步骤S3具体为：

S3.1，采用泰勒展开法得到轨道根数变化关于脉冲速度增量的N阶形式表达式：

S3.2，将表达式中的连续力和时间脉冲按照下述方式结合形成脉冲速度增量：

Δv_r＝F_r·Δt,Δv_t＝F_t·Δt,Δv_n＝F_n·Δt (4)

S3.3，考虑公式(4)中的等价关系，将公式(3)中的连续力表达式变换为脉冲速度增量表达式，从而得到：

其中，η_kpq是经整理后得到的系数；完成了N阶IGVEs的构建。