STK与matlab交互 Astrogator模块（2）

一、LAMBERT问题

       给定转移时间以及两个位置矢量来确定轨道的问题称为兰伯特问题。这在今天是引导航天器与导弹的基本方法。求解这个问题是非常广泛的，各种各样的算法存在来求解这个问题。兰伯特问题是与二体轨道动力学内在联系的，需要一个中心的天体。本文结合STK自带的design tool里面的兰伯特求解器，来利用matlab和STK互联实现一个航天器拦截另一个航天器的问题。

二、STK的LAMBERT求解器

       首先，介绍一下STK中自带LAMBERT求解器的手动操作的过程。

（1）打开STK界面，【Utilities menu】——【Component Browser】——【Design Tools】

(2)点击【Lambert Solver】，然后点击上面的【Duplicate】，重新命名My Lambert Solver（只要是英文即可），再点击My Lambert Solver，打开界面如下所示

（3）注意Solution Option的选择，其有三种选项，一种是固定时间约束(Fixed time)，即给定整个拦截的时间，另外两种分别是最小能量椭圆(Minimum Energyt)和基本椭圆(Minimum eccentric)方式,在这里参照Battin书籍里面的《航天动力学的数学方法》（博主的闲鱼有电子版出售）补充说明一下基本椭圆（Fundamental Ellipse)和最小能量椭圆

(3.1)最小能量椭圆

      从P1到P2点存在着许多种可能的轨道，而最小能量椭圆，根据二体轨道的总能量是,其中为轨道的半长轴，即最小能量椭圆对应着半长轴尽可能小的椭圆轨道，如果P1和P2都在实轴为2a的轨道上，那么

显然只要的最小值是弦长c。因此虚焦点位于弦P1P2上并且对应的半长轴值为

(3.2)基本椭圆

       基本椭圆也被称为最小偏心率椭圆，具有对称性，其图像如下图所示

根据几何关系，基本椭圆的半长轴为

在这里，我们介绍完了上述两种椭圆轨道的基本性质，具体的公式推导及其应用，可以参考《航天动力学的数学方法》，不做过多叙述。

（4）设置求解参数：把历元时间设定为任务分析开始的时刻，因为LAMBER库只能在地球惯性系下进行（Earth Inertial )计算，由于地球章动和平动的影响，不同时刻的地球惯性系的XYZ轴指向在绝对惯性空间是不同的。本文设置为任务开始时刻26 Jan 04:00:00.000 UTCG。

设置任务想定是：蓝星经过三天的时间机动到指定的目标星实行交会，那么我们只要把蓝星的初始时刻的地球惯性系下的坐标作为初始坐标和目标星三天后到达的位置在地球惯性系下的坐标作为结束坐标，这样我们就能顺利实现拦截了。点击【Compute】，给自己的这个机动序列命名，这里命名为one1,点击【Construct Sequence】——【Add to MCS Segment】，然后打开blue的MCS中删除默认的任务序列，鼠标左击【Insert Before】——【one1】——【Run MCS】，这样整个机动过程就完成了，为了验证结果的正确性，左击blue，点击【Report&Graphic Manager】——【Installed Style】——【Inertial Position Velocity】，查看29 Jan 2024 04:00:00.000 UTCG时刻的位置是否到达给定的位置。

下一篇文章，将利用matlab代码的形式完成上述操作，制作不易，希望有所收获的小伙伴，点赞收藏。