c语言计算gps卫星位置.pdfVIP

C 语言计算GPS 卫星位置 1 概述 在用 GPS 信号进行导航定位以及制订观测计划时，都必须已知GPS 卫星在空间的瞬间 位置。卫星位置的计算是根据卫星电文所提供的轨道参数按一定的公式计算的。本节专门 讲解观测瞬间GPS 卫星在地固坐标系中坐标的计算方法。 2 卫星位置的计算 1. 计算卫星运行的平均角速度n 根据开普勒第三定律，卫星运行的平均角速度n0 可以用下式计算： 3    n  G /a  411 0 M ( a)3 14 3 2 式中μ 为 WGS-84 坐标系中的地球引力常数，且μ =3.986005×10 m /s 。平均角速度 n 加上卫星电文给出的摄动改正数Δ n，便得到卫星运行的平均角速度n 0 n=n +Δ n (4-12) 0 2. 计算归化时间t 首先对观测时刻t′作卫星钟差改正 t=t′-Δ t ta a (tt )a (tt )2 0 1 oc 2 oc 然后对观测时刻t 归化 GPS 时系 t =t-t (4-13) oc t t 式中 称作相对于参考时刻 的归化时间（读者注意：toc≠toe）。 oe 3. 观测时刻卫星平近点角M 的计算 M =M +n (4-14) 0 t 式中M 是卫星电文给出的参考时刻toe 的平近点角。 0 4. 计算偏近点角E E =M +esinE (E,M 以弧度计) （4-15） 上述方程可用迭代法进行解算，即先令 E =M ,代入上式，求出 E 再代入上式计算，因 为GPS 卫星轨道的偏心率e 很小，因此收敛快，只需迭代计算两次便可求得偏近点角E 。 5. 真近点角V 的计算 由于:   cosV (cosE e) (1ecosE ) 416 k k k  2    sinV arctg 1e sinE (cosE e)k  k k  417 因此:  2    V  arctg (1e ) sin E (cosE e)k  k k  418 6．升交距角Φ 的计算 k   V   (4  19) k k ω为卫星电文给出的近地点角距。 7. 摄动改正项δ u,δ r,δ i 的计算      C