- 博客(25)
- 收藏
- 关注
RSS订阅原创 STK12 RPO模块学习（3）
Maintain Natural Motion Circumnavigation RPO序列在目标星和追踪星经历不同的力的情况下保持NMC。通常这种差异是由于阻力和太阳光压造成的。这些是主要不同力当执行接近任务的时候，因为重力和相对三体摄动力非常小当相对距离在100km以内的时候。目标卫星和追逐卫星之间的阻力和SRP差异很大程度上是由于不同的尺寸、面积、质量和形状。要有效地使用此序列，必须在正确的位置使用Update_Spacecraft_parameters序列，特别注意这些参数。
2024-05-14 09:36:38
646
原创 STK12 RPO模块学习（2)
这个序列运行卫星直到它达到了下面三个条件之一。1）截至时间。2)圈数到达了限制。3）其他条件，比如近地点。默认情况下，Astrogator使用“Earth HPOP Default v10”预报器。你能够修改呈其他修改器。下面图片展示了一个Coast序列从RBar序列的初始位置500m的地方运行一圈的结果包含四个要素Coast_Definition：从下行菜单中选择停止的方式。包括Duration，Revs，Condition。
2024-05-13 22:35:42
620
原创 STK12 RPO模块学习 （1）
在STK12中，在Astrogator的模块上开发了新的模块（Rendezvous and proximity operations)。轨道交会接近通常来说是一个很复杂的过程。RPO实现需要对轨道动力学有一个清晰的理解，并且对于Astrogator模块具备很强的背景和经验，同时对于数值积分求解器有一个很深的理解。同时掌握基本的编程能力去和Astrogator的脚本交互。AGI公司为了让整个RPO任务更容易上手，推出了该模块。
2024-05-13 21:16:10
928
原创 STK与matlab交互 Astrogator模块（14）
高轨卫星的轨道保持。与任何其它轨道状态一样，地球同步轨道也会受到各种扰动力的影响，这些摄动力会影响GEO卫星在位置方面的稳定性。摄动的主要来源是地球的非地球位势、太阳辐射压力和第三体效应（主要是月球和太阳）。总体影响是经度漂移（其方向取决于经度本身）和倾斜度随时间的变化。GEO卫星运营商将每个航天器放在一个定义明确的盒子里，以免干扰周围空间中的其它卫星。这项活动被称为“站岗”。在这里，我们将了解如何使用自动序列设计和执行台位保持。我们将创建一个东西向测站保持序列（基于经度控制）和一个南北向测站维持序列。
2024-05-11 23:42:30
584
1
原创 STK与matlab交互 Astrogator模块（13）
可以看到前5分钟的卫星按照给定的初始轨道条件进行预报，并且这六颗卫星的数据值均记录在SatJ2000这个6*6的数组里，蓝方1星在8点05时的在J2000坐标系下位置速度为第一行对应的数据，蓝方2星为第二行，依此类推，最后一行为红星的位置速度。第三步，根据人工智能等算法计算6颗卫星这个时刻沿着速度方向施加的脉冲大小，这里由于人工智能算法没有搭建出来，将23次施加脉冲策略提前给出，导入数据DeltaV，其中第一行对应的数据为蓝方1，2，3，4，5星，红星在8点05时刻施加的脉冲大小。下面开始创建仿真场景。
2024-04-27 14:32:31
438
原创 STK与matlab交互 Astrogator模块 （12）
本文结合控制的背景，采用小推力模型的条件下，实现每个时刻给卫星施加一个有限推力的脉冲，即发动机的采用恒推力发动机，同时考虑卫星的质量损耗，控制的变量为推力的方向以及发动机的推力。根据（11）中的数据作为文章的开头，首先获取初始化卫星，分析整个区间的时间序列。注意最后的LVLH报告STK是没有自带的，需要自己定义，整个上面代码的目的是获得目标在LVLH坐标系下的时间序列，为下一步时间段的分析打好基础。
2024-04-13 21:27:47
203
1
原创 STK与matlab交互 Astrogator模块 （11）
前面由于定轨的大作业，关于Astrogator模块的学习有所滞后，在本节将重新聚焦Astrogator模块，在本节中，首先解决的问题是已知两个卫星的轨道六根数，求解其中某一颗卫星LVLH坐标下另一颗卫星的位置速度。这里首先需要自己定义一个报表的内容，例如A，B星的轨道六根数分别为。该报告生成的是卫星B在卫星A的LVLH坐标系下初始的位置和速度。第二步，在STK操作，生成一个B在A星的LVLH坐标系下的坐标，其内容如下。首先利用MATLAB与STK互联生成两颗卫星。
2024-04-13 10:49:29
347
1
原创 使用STK求出坐标转换矩阵
使用STK可以获得一颗卫星的地心惯性系J2000，地心固联坐标系，测站坐标系的x,y,z轴的指向分量，例如通过报告得到了测站坐标系（南东天）三个分量在J2000坐标系下的状态分量。如下图所示，那么从J2000坐标系到测站坐标系的坐标转换矩阵为那么同样的通过测站坐标系在地心固联系下的状态分量，那么从地心固联系到测站坐标系的坐标转换矩阵为那么根据坐标转换的传递性，可以求出J2000坐标系到地固系的坐标转换矩阵。
2024-04-11 16:26:17
516
原创 J2摄动条件下的轨道转移方程
在进行精密定轨的过程中，通常涉及到求解卫星自身的状态转移方程，根据考虑不同的摄动力会对应着不同的状态转移方程，在本文中仅考虑J2项摄动影响，讨论卫星的状态转移方程。根据牛顿第二定律以及J2项摄动情况下的势函数，可以得到下式式中，为地球的势函数，为地球引力常数，为该位置离地球质心的距离，为地球赤道半径，为该位置对应的纬度根据状态转移方程，设卫星的状态空间为，其中为系统观测误差向量，则系统运动方程满足：得到的状态转移方程的微分形式为由于动力学模型的非线性，系统状态转移方程。
2024-04-09 22:52:58
459
原创 最小二乘法中的状态转移方程 （4）
根据（3）的介绍，得到了加速度的变化量，如果采用数值法，可以模拟出轨道的预报，在本节中，将讨论将轨道预报模型线性化，求出轨道的状态转移方程。根据前面的公式，整个状态转移方程可以写为。，而瞬时真天球坐标系下的状态转移方程为，根据（2）中已经求得两个坐标系之间的转换矩阵。通过上式可以得到瞬时真地球坐标系下的状态转移方程。根据上式求得加速度关于位置的关系式为。
2024-04-07 22:24:39
188
原创 最小二乘法的状态转移方程 （3）
在关于未受扰动的开普勒运动的导言一章中，假设地球的总质量集中在坐标系的中心和引力因此，可以用来计算卫星在处所感受到的加速度。为了讨论一个更现实的模型，方便地使用涉及相应重力势梯度的等价表示。通过总结单个质量元素的贡献，这个势的表达式可以很容易地推广到任意的质量分布，取一体积微元，其表达式为，则得到每一体积微元的势函数表达式可以写为这里表示地球内部某个点的密度，是卫星到这个地方的距离球形谐波的展开为了求出上述方程的积分，将距离的导数可以展开为一系列的勒让德多项式。对于。
2024-04-07 20:19:12
366
原创 最小二乘法的状态转移方程 （2）
通常给出的地球指向参数文本格式为：1962 01 01：公历日期37665 ：儒略日-0.0127 ：x极移，单位：角秒0.213 ：y极移，单位：角秒0.0326338：世界时UT1与协调世界时UTC的差0.0017230：LOD 地球日长变化 单位：s0.064261 ：x向的章动值 单位：角秒0.006067 ：y向的章动值 单位：角秒0.000000 ：x向极移的变化量 单位：角秒。
2024-04-07 11:42:40
688
原创 最小二乘法的状态转移方程 （1）
根据精密定轨的介绍，要想使用最小二乘法来对轨道实行微分修正，必须求出卫星的状态转移矩阵，不同的轨道预报模型对应着不同的状态转移矩阵，本文讨论精确轨道预报模型的状态转移方程。首先考虑的力模型是地球自身的引力模型，考虑非球形摄动，地球非球形部分的引力位为。取到71阶常数，同时考虑WGS84地球模型，加载地球各项常数，同时考虑三体摄动，以及太阳光压等常数。
2024-04-06 17:09:43
217
1
原创 轨道定轨中的最小二乘法思想（1）
系统动力学模型（状态模型或状态方程）：其中，为系统状态变量;为系统模型噪声，通常假设为零均值高斯白噪声。系统观测模型（观测方程）其中，为系统观测量，为系统观测噪声，也通常假设为零均值高斯白噪声。假设已获得系统在时刻观测量，就可以估计时刻的系统状态。假设初始估计值为,根据系统微分方程，得到的状态变量，然后根据观测量与状态变量的关系，得到该估计的系统状态方程，得到估计观测量的值残差：实际观测量与估计观测量之间偏差，即。
2024-04-02 21:32:26
813
1
原创 STK与matlab交互 雷达探测和轨道确定 （1）
在本节中，将利用STK软件构建太空目标监视雷达的仿真模型，并且通过给定的空间目标，能够分析优选探测雷达搜索屏设置方案，能够模拟生成雷达跟踪测量数据，之后根据初轨确定和轨道改进的原理，针对给定的探测数据，编写MATLAB程序，实现目标轨道的确定。
2024-03-24 15:44:57
399
1
原创 STK与matlab交互 Astrogator模块（10）
本节重点介绍STK中新的求解器，Sparse Nonlinear Optimizaer(SNOPT)，该求解器能完成一个既定的目标，被花费的函数和一系列约束条件的函数。这个文件迭代修改决策变量或者控制变量来寻找一个可行的、最优的求解方案。SNOPT是一个用于解决大规模优化问题（线性和非线性程序）的软件包。它采用了一种稀疏序列二次规划（SQP）算法，该算法具有有限内存的拉格朗日Hessian拟牛顿近似。SNOPT对于函数和梯度评估昂贵的非线性问题尤其有效。函数应该是平滑的，但不必是凸的。
2024-03-21 11:51:57
439
1
原创 STK与matlab交互 Astrogator模块（9）
前面八节讨论的轨道机动模型均为脉冲形式，既使卫星的速度在瞬时发生改变，在本节中，将调整卫星的机动方式改为有限推力脉冲机动的方式，即Finite形式的脉冲。为了便于理解，本文中结合STK自带帮助文档的教程，简单介绍一下有限推力模型情况下的机动。
2024-03-19 16:45:32
728
1
原创 STK与matlab交互 Astrogator模块（8)
入轨点的示意图，根据CW方程可以知道，一个绕飞椭圆的时间为被绕飞星的一个周期，我们取一个周期内，以60s为步长，离散得到入轨点，如图所示，每个蓝色的点在VVLH坐标系的位置和速度都可以由初始给出的位置和速度递推出来，以椭圆最上方的蓝点为第一个顶点，逆时针开始计数，直到再次回到椭圆最上方的位置，总共有1436个入轨点，按照这个顺序，得到了不同入轨点与脉冲大小之间的关系图，本文选取的绕飞椭圆大小半长轴为10km。根据观察，耗费脉冲的最大值和最小值相差0.3m/s，在实际工程应用中也需要相应的重视。
2024-03-19 09:43:01
498
1
原创 STK与matlab交互 Astrogator模块（7）
根据前六节的介绍，均采用STK自带的Lambert工具箱求解双脉冲转移问题。从前文可以发现，采用Lambert制导求出的脉冲过大，形成的轨迹与实际工程能力严重不符。从这一节开始，将采用STK中的Differential Corrector对双脉冲转移问题进行求解。
2024-03-18 21:21:39
506
原创 STK与matlab交互 Astrogator模块（6）
在（5）的结尾中，给出了一个关于入轨点优化的情况，即只要在绕飞椭圆上的某一点的位置和速度不受外力的情况下，将会一直保持绕飞状态，一个自然绕飞周期通常对应的是目标星的周期，因此在本文中，将在matlab中使用CW方程递推，VVLH坐标系下的相对位置和速度，然后利用matlab与STK交互，将该数据放入STK中实现从相对坐标系到绝对坐标系的转化，即该任务时刻，5 Feb 2024 04:00:00.000 UTCG能够实现绕飞的到达位置有很多种。现在将该位置的选择和转移的时间作为同时考虑的两个变量进行优化。
2024-03-17 20:01:28
522
1
原创 STK与matlab交互 Astrogator模块（5）
通过（5）的操作，最终实现了绕飞的整个过程，但在本文中开始绕飞的入轨点是事先给定的，理论上来说，只要在VVLH坐标系的绕飞椭圆里的任一状态都能够保持绕飞状态，因此，最终的状态点也能据此进行优化，就意味着初状态和末状态都是变化的，初状态与实施机动的时间有关，而末状态与其进入绕飞的轨迹点有关。下一节将会围绕这个问题进行进一步的讨论。创作不易，感兴趣的朋友欢迎咨询闲鱼账号：爱stk的鱼酱，将提供1对1答疑服务。
2024-03-17 15:44:09
682
2
原创 STK与matlab交互 Astrogator模块（4）
由于整个操作过程涉及到与STK交互，使用常规的遍历找最优和遗传算法产生大量的数据通常是不合适的，本文主要介绍了fminsearch函数，使用无导数法计算无约束的多变量函数的最小值。fminsearch有几个输入参数，结构体通常包含的参数：PlotFcns：绘制算法执行过程种的各个进度测量值。
2024-03-16 22:14:03
407
原创 STK与matlab交互 Astrogator模块（3）
结合（2）介绍的STK中自带的LAMBERT求解器，本文主要介绍matlab与STK交互，自动化完成上述任务需求。
2024-03-15 10:34:59
954
原创 STK与matlab交互 Astrogator模块（2）
（3）注意Solution Option的选择，其有三种选项，一种是固定时间约束(Fixed time)，即给定整个拦截的时间，另外两种分别是最小能量椭圆(Minimum Energyt)和基本椭圆(Minimum eccentric)方式,在这里参照Battin书籍里面的《航天动力学的数学方法》（博主的闲鱼有电子版出售）补充说明一下基本椭圆（Fundamental Ellipse)和最小能量椭圆。为轨道的半长轴，即最小能量椭圆对应着半长轴尽可能小的椭圆轨道，如果P1和P2都在实轴为2a的轨道上，那么。
2024-03-14 21:10:24
1067
原创 STK与matlab交互 Astrogator模块（1）
在本文中，将详细介绍matlab2022b与STK11.6版本的交互，本文中尝试着将STK对指定卫星的轨道数据导入matlab中，以数组的形式进行保存，以便于在下一步利用matlab的轨道数据进行处理。（1）打开STK（2）启动matlab，新建脚本，输入命令，这样end%% 根据你的需要设定场景的名称（3）接下来在matlab脚本中写入下列命令，设置场景的分析时间% 场景开始时间% 场景结束时间。
2024-03-13 21:38:44
606
4
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人