摘要：为了更充分利用低轨重力卫星的高精度观测数据，根据卫星轨道的扰动理论，导出了应用卫星轨道与星间距离观测值联合反演地球重力场模型的算法．该算法的实质是将牛顿运动方程在卫星轨道处进行展开，转化为第二类Volterra积分方程，并采用基于移动窗口的9次多项式内插公式进行数值求解．给出了该算法的观测方程，用QR分解法消去局部参数矩阵，最后采用预条件共轭梯度法求解法方程．利用GRACE卫星2008-01-01～200808-01时间段内的轨道及星间距离观测数据，解算了120阶次的地球重力场模型SWJTU-GRACE01S，该模型在120阶处的阶方差为1．58×10，大地水准面差距累计误差为22．29cm，与美国GPS水准网比较的标准差为0．793m，结果表明：SWJTU—GRACE01S模型精度介于EIGEN—GRACE01S与EIGEN—GRACE02S模型之间，从而验证了该算法的有效性．

关键词：卫星轨道扰动 地球重力场模型 牛顿运动方程 星间距离 多项式内插 

单位：西南交通大学测量系 成都610031 同济大学测量与国土信息工程系 上海200092 中国科学院测量与地球物理研究所 武汉430077