GPS卫星坐标位置计算

GPS卫星坐标位置计算，以北京房山 2019/1/2 上午6:01:00 GPS数据为例 勘误：观测时刻计算时少加了两整天，即应改为：t = 2*24*60*60+20*60*60#观测时刻t

import numpy as np
from numpy import *
from math import *

#bjfs0010.19n 32.0 kB	2019/1/2 上午6:01:00
#北京房山 2019/1/2 上午6:01:00 GPS数据
#GPS广播星历
 #1 19  1  1 20  0  0.0 -1.383270137012D-04 -6.139089236967D-12  0.000000000000D+00
 #   6.700000000000D+01 -3.762500000000D+01  3.919091817334D-09  2.254689281430D+00
 #  -1.918524503708D-06  8.239355636761D-03  1.081638038158D-05  5.153665266037D+03
 #   2.448000000000D+05 -9.685754776001D-08 -3.484545068489D-01 -1.117587089539D-07
 #   9.742196594735D-01  1.775312500000D+02  6.903077333794D-01 -7.651032981949D-09
 #   4.271606500985D-10  1.000000000000D+00  2.034000000000D+03  0.000000000000D+00
 #   2.000000000000D+00  0.000000000000D+00  5.587935447693D-09  6.700000000000D+01

#1.计算卫星运行的平均角速度n
muon = double（3.986005*pow（1014））#WGS