矢量化双体传播器（开普勒宇宙变量...）

简介：
Keplerian 通用变量二体传播方法采用 Newton-Raphson 迭代算法，从而构建了一个高效的传播程序，它避免了将坐标系从笛卡尔坐标系转换成 Keplerian 坐标系，然后再将其转换回笛卡尔坐标系的繁琐步骤。 这种方法能够适应各种轨道形状，包括圆形、偏心、抛物线以及双曲线轨道，并且能够精确地提供小于 0.5 公里的位置误差，该误差值是通过与解析二体传播例程的 STK 实现进行比较所获得的。 具体的使用示例如下： >> [r,v] = keplerUniversal(r0,v0,t,mu); 其中： 输入： r0 是在 ECI 参考坐标系中定义的初始位置向量，单位为千米 (km)，其维度为 [3xN]； v0 是在 ECI 参考坐标系中定义的初始速度矢量，单位为千米/秒 (km/s)，其维度也为 [3xN]； t 是以秒为单位的时间向量，其维度为 [1xN]； mu 是中心天体的重力常数。若未指定该参数，则默认值为地球的重力常数（双）。 输出： r 是在 ECI 参考坐标系中计算出的传播位置向量，单位为千米 (km)，其维度同样为 [3xN]； v 是在 ECI 参考坐标系中计算出的传播速度矢量。