星上下传的辅助数据由于各种原因经常发生误码，导致辅助数据存在错误值。为了剔除这些错值，需要稳健的剔除粗差的算法。现在针对星上直传的GPS测量数据，描述一种稳健的剔除错误点的算法。

算法输入原始的GPS文件，输出剔除错误点后的GPS文件。

（为何不直接选用卫星位置，而要选用星地距离来鉴别呢？这是因为，虽然卫星位置在地心直角坐标系中的三个坐标分量随时间变化也是很平滑的，但很多时候，即使坐标存在错误其随时间变化也还是那么平滑，常常不能有效地反应卫星测量数据正确与否，因此不被选做鉴别使用。）

算法原理：首先计算各时刻下(x)卫星与地心之间的距离(y)；再对测量时刻x和星地距离y进行多项式拟合，按照粗差点大于三倍中误差的原理剔除误差最大的点；再重复：对剩下的测量时刻(x)星地距离(y)进行拟合以及剔除误差最大的点，直到所有星地距离误差dy小于3倍中误差或达到指定的阈值为止。

需要工具软件请联系jsy96@qq.com