卫星定位测量 卫星的位置计算 卫星的位置计算.docVIP

卫星的位置计算 利用GPS卫星进行导航和定位，就是根据已知的卫星轨道参数计算出卫星瞬时位置，通过观测和数据处理，确定接收机的位置和载体的运动速度。所以，获取准确的卫星轨道参数，计算出卫星在观测瞬间的位置，是GPS导航定位的基础。因为GPS系统坐标系统采用WGS-84坐标系统。为了计算卫星在WGS-84大地坐标系中的位置，首先需要计算卫星在其轨道平面内的位置。此时定义：原点与地心相重合，轴指向升交点，轴在轨道平面内垂直于轴，我们称其为轨道平面直角坐标系，它是一种过渡性的坐标系。再进行坐标系的转换，将卫星在其轨道的坐标转换到地面直角坐标系下。 1 ）卫星在其轨道平面内的坐标 1.1 计算卫星运行的平均角速度n 卫星运行的平均角速度 式中为卫星轨道的长半轴，为地球引力常数（含大气层），它们都有确定值。利用导航电文中给出的摄动改正数，计算观测时刻卫星的平均角速度 （1） 1.2计算归化时间 导航电文提供的轨道参数是相应于参考时刻时的数值，为了求得观测时刻的参数，需求出此时刻相对于参考时刻的时间差： （2 ） 在GPS时间系统中，时间是从一周的开始（星期日子夜）连续以秒计算的，所以计算归化时间时应顾及到一个星期（604 800s）的开始或结束。亦即当＞302 400s时，应减去604 800s。当＜－302 400s时，应加上604 800s。 1.3计算真近点角 首先计算观测瞬间的卫星平近点角。因为导航电文中已经给出参考时刻的平近点角，所以 （3） 其次计算偏近点角。根据导航电文中给出的偏心率和上面计算的平近点角，根据公式根据开普勒方程 可以写出： （4） 上述开普勒方程仍按前述的迭代法计算，即先令，代入上式。通常只需两次迭代计算即可求得偏近点角。 然后再计算出真近点角。由式可得： （5） 1.4计算升交距角 根据计算出的其近点角和导航电文提供的近地点角距，按下式计算升交距角 （6） 1.5计算经过摄动改正的升交距角、卫星的地心距离及轨道倾角 在卫星导航电文中，给出了6个卫星轨道摄动的修正参数。根据这些参数，计算出轨道参数s的修正量： （7） 经过摄动改正的升交距角、卫星的地心距离、轨道倾角则为： （8） 1.6计算卫星的轨道平面直角坐标 （9） 2 ）卫星在地心空间直角坐标系中的位置计算 2.1 计算观测时刻的升交点经度 为了计算卫星在地心空间直角坐标系中的位置，必须知道卫星升交点在观测时刻的精度。 图1广播星历参数 假若观测时刻的升交点赤经为（升交点与春分点所对应的地心夹角），观测时刻时的春分点格林尼治恒星时为GAST，由图1可知，此时的升交点经度应为 （10） 如果参考时刻的升交点赤经为，其变率为，则观测时刻时的升交点赤经为 （11） 其次，式（10 2-26）中的GAST将随地球自转而不断增加，其增值速率即为地球自转的速率。设一个星期开始（星期日子夜）时刻的格林尼治恒星时为GAST（），则 （12） 式中：——710-5rad/s； 导航电文给我们提供的不是时的升交点赤经，而是始于格林尼治起始子午线到升交点的准经度，它们之间的关系是： （13） 将式（11）、（12）、（13）一并代入式（9），则得升交点的经度为： （14） 2.2 计算卫星在地心空间直角坐标系中的坐标 由式（9)可求出卫星在其轨道平面直角坐标系中的坐标，在该坐标系统中卫星的空间位置可表示为： （15） 此时轴过地心指向轨道平面的垂直方向。 这样，轨道直角坐标系（，，）与 地心直角坐标系（X，Y，Z）具有相同的原点（参阅图2和图3），其差别仅是坐标轴的定向不同。为了使两坐标系取得一致，可将坐标系（，，）作如下旋转：（1）绕轴方向顺转角度，使轨道平面与赤道平面重合，轴与Z轴相重合；（2）绕轴顺转角度，使、轴与X、Y轴同时重合。 图2 卫星轨道参数