1)在test_EKF_GEOS3中，我们定义了全局变量:const(天文和数学常数)，CS(非归一化重力场系数)，AuxParam(表示轨道传播器的模型参数和时间的结构)，eopdata(地球方向参数)，PC (JPL星历计算的行星系数)。
2)从下面的链接下载地球方向参数，删除标题和eop19620101.txt之类的文本，然后重命名它(如果它过时了)。
3)从eop19620101.txt中读取地球方位参数。
4)观测数据来自GEOS3.txt，减去了“Fundamentals of Astrodynamics and Applications, 4th Edition, 2013”一书中表4-4中的偏差。测量噪声设置如表4-4所示。在此之后，使用position计算站点的ecef位置。M，并选择两组观测值(方位角、仰角和距离)。然后用LTCMatrix计算LTC矩阵(从格林威治子午线系统到局部切线坐标的转换)。m，由Elements求得该历元(第一次测量时间)卫星状态向量的初始猜测值。m和State.m。
5)扩展卡尔曼滤波器的初始参数和轨道传播器参数(AuxParam。Mjd_UTC = Mjd_UTC;AuxParam。N = 20;AuxParam。M = 20;AuxParam。太阳= 1;AuxParam。月亮= 1;AuxParam。行星= 1;)设置。
6)将epoch的状态向量传播到第二次测量的时间并更新。在此之后，更新的状态向量被传播到下一次测量的时间并更新。这个过程一直持续到最后一次测量。最后，在最后一次测量时更新的状态向量被传播到第一次测量时。将结果从TOD转换为ECI坐标系，并与真实状态向量进行比较。
请注意数值积分是在TOD坐标系下完成的。所以，Accel中只考虑地球自转。m和VarEqn.m。
引用:
O. Montenbruck, E. Gill，“卫星轨道:模型、方法和应用”，施普林格Verlag，海德堡，2005年。
D. Vallado，“天体动力学基础与应用”，2013年第4版。
引用
Meysam Mahooti(2023)。初始定轨(扩展卡尔曼滤波)(//www.ru-cchi.com/matlabcentral/fileexchange/55627-initial-orbit-determination-extended-kalman-filter)， MATLAB中央文件交换。取得済み．
初始定轨(扩展卡尔曼滤波)
| バ，ジョン | 公開済み | リリ，スノ，ト | |
|---|---|---|---|
| 2.2.1 | 修订日期为2022-10-25。 |
||
| 2.2.0 | test_EKF_GEOS3.mwas modified. |
||
| 2.1.2 | 添加了文档，test_EKF_GEOS3。M被修改。 |
||
| 2.1.1 | 修订日期为2022-10-22。 |
||
| 2.1.0.0 | 香波-戈登积分器被Radau IIA积分器取代。 |
||
| 2.0.0.0之间 | 添加DE430全矩阵。 |
||
| 1.1.1.1 | 2020年3月16日修订。 |
||
| 1.0.1.0 | 添加文档。 |
||
| 1.0.0.0 | 完成镜像的添加。 |
您也可以从以下列表中选择一个网站: