开普勒状态转移矩阵：给定初始位置和速度，计算轨道的开普勒状态转移矩阵-matlab开发

Phi = keplerSTM(x0,dt,mu) 将返回状态转移矩阵 Phi 对于由初始位置和速度 (x0) 描述的一组对象具有引力参数 mu，沿开普勒轨道传播时间 dt。 时间 dt 的位置和速度由 Phi*x0 给出。 x0 必须是长度为 6n 的行向量，其中 n 是行星的数量。 对于每个行星，x0 必须包括 3 个位置和 3 个速度值顺序为 [r1,r2,r3,v1,v2,v3]。 这些位置和速度是在笛卡尔坐标系中以中心物体测量（即明星）在原点。 对于多个行星，只需堆叠这些向量中的几个相互叠加。 mu 必须包含 n 值等于 G(m+ms) 其中 G 是重力常数，m 是轨道物体的质量，ms 是中心物体的质量。 dt 是传播时间的标量值。 注：所有单位应一致。 如果位置在 AU 和速度以 AU/天为单位，那么 dt 必须以天为单位，而 mu 必须以AU^3/天^2。 这是使用 S