太原建设厅官方网站,上传户型图生成效果图,专注wordpress主题,个人网站建设平台目录
1.计算卫星运动的平均角速度 n
2.计算观测瞬间卫星的近地点角
3.计算偏近点角
4.计算真近点角 f
5.计算升交角距
6.计算摄动改正项
7.进行摄动改正
8.计算卫星在轨道面坐标系中的位置
9.计算观测瞬间升交点的经度 L
10.计算卫星在瞬时地球坐标系中的位置
11.…目录
1.计算卫星运动的平均角速度 n
2.计算观测瞬间卫星的近地点角
3.计算偏近点角
4.计算真近点角 f
5.计算升交角距
6.计算摄动改正项
7.进行摄动改正
8.计算卫星在轨道面坐标系中的位置
9.计算观测瞬间升交点的经度 L
10.计算卫星在瞬时地球坐标系中的位置
11.计算卫星在协议地球坐标系中的位置 使用广播星历进行GPS卫星位置的计算的详细步骤如下
关于在广播星历中各参数的分布情况可以参考下图 1.计算卫星运动的平均角速度 n
首先根据广播星历中给出的参数 计算参考时刻(TOE) 的平均角速度 式中GM为万有引力常数 G 与地球总质量 M 之乘积其值为 然后根据广播星历中给定的摄动参数 计算观测时刻卫星的平均角速度 n 注观测时刻地面站观测的时间参考时刻卫星发出信号的时间
2.计算观测瞬间卫星的近地点角 式中 为参考时刻 时的平近点角由广播星历给出。
问为什么要用参考时刻 来替代卫星过近地点时刻 来计算呢?
答因为广播星历每 2h 更新一次将参考时刻设在中央时刻时外推间隔小于等于 1h。而卫星的运行周期为 12h 左右采用卫星过近地点时刻 来计算时外推间隔最大有可能达 6h 。用 来取代卫星过近地点时刻 后外推间隔将大大减小用较简单的模型也能获得精度较高的结果。
3.计算偏近点角
用弧度表示的开普勒方程为 用角度表示的开普勒方程为 解上述方程可用迭代法或微分改正法。
4.计算真近点角 f 式中e 为卫星轨道的偏心率由广播星历给出。
所以真近点角 f 可以表示为 5.计算升交角距 式中 为近地点角距由广播星历给出。
6.计算摄动改正项 、 、
广播星历中给出了下列6个摄动参数 据此可以求出引起的升交角距 u 的摄动改正项 、卫星矢径 r 的摄动改正项 和卫星轨道倾角 i 的摄动改正项 。
卫星矢径 r 接收机与卫星之间的绝对距离直线距离以米为单位。
计算公式如下: 7.对 、 、 进行摄动改正 式中a 为卫星轨道的长半径a 由广播星历给出 为 时刻的轨道倾角由广播星历中的开普勒六参数给出 为轨道倾角 i 的变化率由广播星历中的摄动九参数给出。
8.计算卫星在轨道面坐标系中的位置
在轨道平面直角坐标系中(坐标原点位于地心, 轴指向升交点)
卫星的平面直角坐标为:
卫星矢径 r 、升交角距 u 9.计算观测瞬间升交点的经度 L
设参考时刻 时升交点的赤经为 , 升交点对时间的变化率为
那么观测瞬间 t 的升交点赤经 应为: 可从广播星历的摄动参数中给出。
设本周开始时刻(星期日 0 时)格林尼治恒星时为 则观测瞬间的格林尼治恒星时为 其中 为地球自转角速度其值为
t 为本周内的时间(s)从本周开始时刻起算至本周观测时刻瞬间的秒数
显然上述算法中把地球自转看成是完全匀速的未顾及地球自转的不均匀性。
这样就可以求得观测瞬间升交点的经度值 L 令 则有 注意广播星历中给出的并不是参考时刻 的升交点赤经 而是该值与本周起始时刻的格林尼治恒星时 之差(即上述的 )。
10.计算卫星在瞬时地球坐标系中的位置
已知升交点的大地经度 L 以及轨道平面的倾角 i 后
就可通过两次旋转方便地求得卫星在地固坐标系中的位置 11.计算卫星在协议地球坐标系中的位置
观测瞬间卫星在协议地球坐标系中的位置为 在GPS定位中常常需要多次计算卫星的位置和速度如用上述方法计算需占用较多的内存和计算时间。为此常将卫星星历用一个时间多项式来表示在内存中仅保存该多项式的系数供计算时调用。
在各种多项式中切比雪夫多项式的拟合效果最佳即使在该时间段的两端近似性也很好。用 n 阶切比雪夫多项式来逼近时间段 中的卫星星历时先将变量 变换为变量 于是卫星坐标就可以表示为 式中n 为多项式的阶数 为切比雪夫多项式的系数。 根据已知的卫星坐标用最小二乘法拟合出多项式系数 就可用上式计算出该时段中任一时刻的卫星位置。
切比雪夫多项式 的递推公式如下
