c语言计算gps卫星位置,C语言学习知识-详细计算GPS卫星位置.doc

\\

C语言计算GPS卫星位置

1 概述

在用GPS信号进行导航定位以及制订观测计划时，都必须已知GPS卫星在空间的瞬间位置。卫星位置的计算是根据卫星电文所提供的轨道参数按一定的公式计算的。本节专门讲解观测瞬间GPS卫星在地固坐标系中坐标的计算方法。

2 卫星位置的计算

1. 计算卫星运行的平均角速度n

根据开普勒第三定律，卫星运行的平均角速度n0可以用下式计算：

式中μ为WGS-84坐标系中的地球引力常数，且μ=3.986005×1014m3/s2

n=n0+Δn (4-12)

2. 计算归化时间tk

首先对观测时刻t′作卫星钟差改正

t=t′-Δt

然后对观测时刻t归化到GPS时系

tk=t-toc (4-13)

式中tk称作相对于参考时刻toe的归化时间(读者注意：toc≠toe)。

3. 观测时刻卫星平近点角Mk的计算

Mk=M0+ntk (4-14)

式中M0是卫星电文给出的参考时刻toe的平近点角。

4. 计算偏近点角Ek

Ek=Mk+esinEk(Ek,Mk以弧度计) (4-15)

上述方程可用迭代法进行解算，即先令Ek=Mk,代入上式，求出Ek再代入上式计算，因为GPS卫星轨道的偏心率e很小，因此收敛快，只需迭代计算两次便可求得偏近点角Ek。

5. 真近点角Vk的计算

由于:

因此:

6．升交距角Φk的计算

ω为卫星电文给出的近地点角距。

7. 摄动改正项δu,δr,δi的计算

<