test_EKF_GEOS3 1)中,我们定义全局变量:常量(天文和数学常数),c(非规范重力场系数),AuxParam(结构表示模型参数和时间轨道传播算子),eopdata(地球定向参数)、电脑(喷气推进实验室星历表计算行星系数)。
2)从以下链接下载地球定向参数,像eop19620101删除标题和文本。txt,然后重命名它(如果它过时了)。
3)地球定向参数从eop19620101.txt读取。
4)观察从GEOS3读取。txt和偏见,从表4 - 4本书的“天体动力学原理及应用(第四版),2013年,“减去。然后从表4 - 4测量声音设置。在那之后,车站ecef位置计算使用位置。m,两组观察(方位角、仰角和范围)。随后,LTC矩阵(从格林威治子午线系统转换到本地切线坐标)由LTCMatrix计算。m,最初想卫星状态向量的时代(第一次测量的时间)发现的元素。m和State.m。
5)扩展卡尔曼滤波器的初始参数和轨道参数(AuxParam传播者。Mjd_UTC = Mjd_UTC;AuxParam。n = 20;AuxParam。m = 20;AuxParam。太阳= 1;AuxParam。月亮= 1;AuxParam。行星= 1;)。
6)时代的状态向量是传播到第二次测量的时间和更新。之后,状态向量更新传播到下一个测量的时间和更新。这一过程持续进行直到最后测量的时间。最后,更新状态向量的时候最后一个测量是传播到第一次测量的时间。结果被转换从TOD ECI坐标系统并与真实状态向量。
请注意,数值积分在托德坐标系统。所以,在Accel只考虑地球自转。m和VarEqn.m。
引用:
o . Montenbruck e·吉尔,“卫星轨道:模型、方法和应用,“施普林格-,海德堡,2005年。
d . Vallado“天体动力学原理和应用程序,第四版,2013年。
引用作为
Meysam Mahooti (2023)。初始轨道计算(扩展卡尔曼滤波)(//www.tatmou.com/matlabcentral/fileexchange/55627-initial-orbit-determination-extended-kalman-filter), MATLAB中央文件交换。检索。