STK与MATLAB联动第二期（2）

1. 轨道外推（J2、J4、HPOP）

由于地球并非理想的球形，且卫星运动过程受到各种引力作用，所以，卫星的运动与二体运动有所偏离，通常称为有摄运动。因此，在实际航天仿真中，常用的外推算法包括J2摄动外推、J4摄动外推、高精度轨道外推等。

（1）J2摄动外推(J2Pertubation)

J2摄动（一阶）外推算法考虑了由于地球非球形引起的轨道根数的长周期变化，这种算法不考虑大气阻力以及日月引力。

（2）J4摄动外推(J4Pertubation)

J4摄动（二阶）外推算法考虑了由于地球非球形引起的轨道根数的长周期变化，这种算法不考虑大气牵引力以及日月引力。J4外推模型包括了J2的一阶和二阶影响以及J4的一阶影响，不包括引起长周期变化的J3系数。

（3）高精度轨道外推(HPOP)

高精度轨道外推能够处理员、椭圆、抛物线、双曲线轨道，其距离涵盖地球表面到月球轨道以及更远的距离。HPOP与二体运动、J2、J4轨道外推一样利用相同的轨道根数设置航天器在某一历元时刻的状态，其具体设置以及matlab代码见上一期。

1. TLE轨道根数

1. TLE介绍

TLE轨道根数是一组用于确定空间目标位置和速度的轨道数据，它以文本格式给出，由两行组成，有效字符为数字0-9、大写字母A-Z、正负号、空格和句点。格式如下：

格式1：AAAAAU YYLLLPPP BBBBB.BBBBBBBB .CCCCCCCC DDDDD-D EEEEE-E F GGGGZ

• AAAA为空间目标编号；
• U为安全保密密级；
• YYLLLPPP为国际标志符，其中YY为发射年度，LLL为该年度发射序列，PPP为部件序列号；
• BBBBB. BBBBBBBB为历元时刻，其中前两位为年份，其次三位为天数，小数是不足一天的时间；
• .CCCCCCCC为轨道平动一阶导数；
• DDDDD-D为轨道平动二阶导数；
• EEEEE-E为BSTAR阻力系数，F为轨道计算模型；
• GGGG为星历数；
• Z为校验和

格式2：AAAAA HHH.HHHH III. IIII JJJJJJJ KKK. KKKK MMM. MMMMNN. NNNNNNNNRRRRRZ

• AAAA为空间目标编号；
• HHH.HHHH为轨道倾角；
• III. IIII为升交点赤经；
• JJJJJJJ为轨道偏心率，小数点在首位；
• KKK. KKKK为近地点幅角；
• MMM. MMMM为平近地角；
• NN. NNNNNNNN为平均运动角速度，单位为圈/天；
• RRRRR为历元时刻的圈数；
• Z为校验和

以中国空间站2023年8月31日公开轨道信息为例，具体可见：

中国空间站轨道参数-20230831_中国载人航天官方网站

1 48274U 21035A 23243.00000000 .00014080 00000-0 23907-3 0 9997

2 48274 41.4743 61.0218 0010871 268.9629 129.1488 15.61490428133535

STK中对应的参数设置

MATLAB代码示例

可通过直接在matlab指令中输入TLE数据或调用TLE文件两种方式进行卫星轨道的设定。

1. 指令输入

SetState */Satellite/Satellite1 TLE "1 10637U 78012A   04113.40484266  .00000000  00000-0 10000-3 0  9453" "2 10637 040.1939 040.6913 1575431 119.1582 257.4005 01.000538770048625"

1. TLE文件调用

TLE文件可从CelesTrak: Current GP Element Sets 网站中获取，本文选择的是https://celestrak.org/NORAD/elements/gp.php?GROUP=iridium&FORMAT=tle

ImportTLEFile * "C:\Users\raind\Desktop\gp.txt"