CSDN上还有一个,但是结果算出来有问题,公式修改过,和精密星历对比误差单位在米级。