最新3 GPS卫星的坐标计算汇总Word格式.docx

1、远地点fs（2）开普勒第二定律 卫星的地心向径在单位时间内所扫过的面积相等。表明卫星在椭圆轨道上的运行速度是不断变化的，在近地点处速度最大，在远地点处速度最小。（3）开普勒第三定律 卫星运行周期的平方与轨道椭圆长半径的立方之比为一常量，等于GM的倒数。 假设卫星运动的平均角速度为n，则n=2/Ts，可得 当开普勒椭圆的长半径确定后，卫星运行的平均角速度也随之确定，且保持不变。3.1.3卫星运动的轨道参数由开普勒定律可知，卫星运动的轨道，是通过地心平面上的一个椭圆，且椭圆的一个焦点与地心相重合。而确定椭圆的形状和大小至少需要两个参数，即椭圆的长半径as及其偏心率es（或椭圆的短半径bs）。另外，

2、为确定任意时刻卫星在轨道上的位置，需要一个参数，一般取真近点角fs，即在轨道平面上，卫星与近地点之间的地心角距，该参数为时间的函数，它确定了卫星在轨道上的瞬时位置。参数as，es，fs唯一地确定了卫星轨道的形状、大小以及卫星在轨道上的瞬时位置。如图3-1所示，称之为轨道椭圆形状参数。图3-1 开普勒轨道参数为了确定该椭圆在上述坐标系中的方向，尚需三个参数。一般采用开普勒轨道参数（图31），或称开普勒轨道根数。升交点的赤经，即在地球赤道平面上，升交点与春分点之间的地心夹角（升交点，即当卫星由南向北运行时轨道与地球赤道面的一个交点）。i轨道面的倾角，即卫星轨道平面与地球赤道面之间的夹角。上两个参数

3、，唯一地确定了卫星轨道平面与地球体之间的相对定向，称之为轨道平面定向参数。s近地点角距，即在轨道平面上，升交点与近地点之间的地心夹角，这一参数表达了开普勒椭圆在轨道面上的定向，称之为轨道椭圆定向参数。在此，参数as、es、i、s和fs所构成的坐标系统，通常称为轨道坐标系统。其中，参数as、es、i、s的大小，则是由卫星的发射条件决定。在该系统中，当6个轨道参数一经确定后，卫星在任一瞬间相对地球体的空间位置及其速度，便可唯一地确定。真近点角fs的计算 在6个轨道参数中，参数as、es、i、s的大小，则是由卫星的发射条件决定。所以计算卫星瞬时位置的关键在于计算参数fs。图33 偏近点角与真近点角如

4、图，为了计算真近点角，引入两个辅助参数Es偏近点角 Ms平近点角。 Ms是一个假设量，当卫星运动的平均角速度为n，则Ms = n （ t - t0 ），t0为卫星过近地点的时刻，t为观测卫星时刻。平近点角与偏近点角间存在如下关系（开普勒方程）： Es = Ms + es sinEs。 真近点角fs和偏近点角Es之间的关系：又可得：或者：3.2 GPS卫星的导航电文（数据码）GPS卫星的导航电文（简称卫星电文又叫数据码（D码）：是用户用来定位和导航的数据基础。它主要包括：卫星星历、时钟改正、电离层时延码转换到捕获P码的信息。它的基本单位是长1500bit的一个主帧（如图4-1所示），传输速率是5

5、0bit/s,30秒钟传送完毕一个主帧。一个主帧包括5个子帧，第1、2、3子帧每30秒钟重复一次，内容每小时更新一次。第4、5子帧的全部信息则需要750秒钟才能够传送完。即第4、5子帧是12.5分钟播完一次，然后再重复之，其内容仅在卫星注入新的导航数据后才得以更新。3.2.1遥测码（TLW，即TelemetryWord）遥测码位于各子帧的开头，它用来表明卫星注入数据的状态，以次指示用户是否选用该颗卫星。3.2.2转换码（HOW,即HandOverWord）转换码位于每个子帧的第二个字码，其作用是提供帮助用户从所捕获的C/A码转换到捕获P码的Z计数，它表示从每星期天零时到星期六24小时，P码子码

6、X1的周期（1.5秒）重复数。3.2.3第一数据块第一数据块位于第1子帧的第310字码，它的主要内容包括：标识码，时延差改正；星期序号；卫星健康状况；数据龄期；卫星时钟改正系数等。1.时延差改正时延差改正就是载波L1、L2的电离层时延差。当使用单频接收机时，为了减小电离层的影响，提高定位精度，要用改正观测结果，双频接收机可通过L1，L2两频率的组合来消除电离层的影响，不需要此项改正。2.数据龄期AODC卫星时钟的数据龄期AODC是时钟改正数的外推时间间隔，它指明卫星时钟改正数的置信度。AODC=Toc-tl （4-1）式中：toc为第一数据块的参考时刻；tl是计算时钟改正参数所用数据的最后观测

7、时间。3.星期序号WNWN表示从1980年1月6日子夜零点（UTC）起算的星期数，即GPS星期数。4.卫星时钟改正GPS时间系统是以地面主控站的主原子钟为基准。由于主控站主钟的不稳定性，使得GPS时间和UTC时间之间存在着差值。地面监控系统通过临测确定出这种差值，并用导航电文播发给广大用户。每一颗GPS卫星的时钟相对GPS时系存在着差值，需加以改正，这是卫星时钟改正：ts=a0+a1（t-toc）+a2（t-toc）2（4-2）式中，a0,a1,a2含义见3.4节。3.2.4第二数据块包含第2和第3子帧，其内容表示GPS卫星的星历，这些数据为用户提供了有关计算卫星运动位置的信息。描述卫星的运行

8、及其轨道的参数包括下列三类。（如图4-2所示）1.开普勒六参数这6个参数为：，e,i0,0,M0（其含义同3.4）。2.轨道摄动九参数这9个参数为：n，,，Cuc,Cus,Crc,Crs,Cic,Cis（其含义同3.4）3.时间二参数（1）从星期日子夜零点开始度量的星历参考时刻toe，变化与0604800秒；（2）星历表的数据龄期AODE，有：AODEtoe-t1 （）式中t1为作预报星历测量的最后观测时间，因此AODE就是预报星历的外推时间长度。3.2.5第三数据块第三数据块包括第4和第5两个子帧，其内容包括了所有GPS卫星的历书数据。当接收机捕获到某颗GPS卫星后，根据第三数据块提供的其他

9、卫星的概略星历、时钟改正、卫星工作状态等数据，用户可以选择工作正常和位置适当的卫星，并且较快地捕获到所选择的卫星。1.第4子帧（1）第2，3，4，5，7，8，9，10页面提供第2532颗卫星的历书；（2）第17页面提供专用电文，第18页面给出电离层改正模型参数和UTC数据；（3）第25页面提供所以卫星的型号、防电子对抗特征符和第2532颗卫星的健康状况；（4）第1，6，11，12，16，19，20，21，22，23，24页作备用，第13，14，15页为空闲页。2.第5子帧（1）第124页面给出第124颗卫星的历书；（2）第25页面给出第124颗卫星的健康状况和星期编号。在第三数据块中，第4和第

10、5子帧的每个页面的第3字码，其开始的8个比特是识别字符，且分成两种形式：（a）第1和第2比特为电文识别（DATAID）；（b）第38比特为卫星识别（SVID）。GPS卫星广播星历预报参数及其定义参数参数定义Toe（s）星历表参考历元（秒）IODE（AODE）星历表数据量（N）M0（rad）按参考历元toe计算的平近点角（弧度）n（rad/s）由精密星历计算得到的卫星平均角速度与按给定参数计算所得的平均角速度之差（弧度）e轨道偏心率（m）轨道长半径的平方根（0.5m）0（rad）按参考历元toe计算的升交点赤径（弧度）i0（rad）按参考历元知计算的轨道倾角（弧度）（rad）近地点角距（弧度）（

11、rad/s）升交点赤径变化率（弧度/秒）I（rad/s）轨道倾角变化率（弧度/秒）Cuc（rad）升交距角（+fS）的余弦调和项改正的振幅（弧度）Cus（rad）升交距角（+fS）的正弦调和项改正的振幅（弧度）Crc（m）轨道半径的余弦调和项改正的振幅（米）Crs（m）轨道半径的正弦调和项改正的振幅（米）Cic（rad）轨道倾角的余弦调和项改正的振幅（弧度）Cis（rad）轨道倾角的正弦调和项改正的振幅（弧度）GPD周数（周）Tgd电离层延迟改正（秒），IODC星钟的数据量（N）a0卫星钟差时间偏差（秒）a1卫星钟速（秒/秒）频率偏差系数a2卫星钟速变率（秒/秒2）漂移系数卫星精度（N）卫星健

12、康GPS n文件的文件头格式说明文件头标签说明格式RINEX VERSION/TYPE- RINEX格式的版本号（如2.10）- 文件类型（在本文件中为“N”）：F9.2,11XA1,19XPGM/RUN BY/DATE- 创建本数据文件所采用程序的名称- 创建本数据文件单位的名称- 创建本数据文件的日期A20COMMENT- 注释行A60ION ALPHA- 历书中的电离层参数A0A3（第四子帧的第18页）2X,4D12.4ION BETA- 历书中的电离层参数B0B33X,4D12.4DELTA-UTC:A0,A1,T,W 用于计算UTC时间的历书参数（第四子帧的第18页）- A0，A1：

13、多项式系数- T：UTC数据的参考时刻- W：UTC参考周数，为连续计数，不是1024的余数3X,2D19.12I9LEAP SECONDS- 由于跳秒而造成的时间差I6END OF HEADER- 文件头的最后一个记录60X GPS n文件的数据记录节说明观测值记录PRN号/历元/卫星钟- 卫星的PRN号- 历元：TOC（卫星钟的参考时刻） 年（2个数字，如果需要可补0） 月，日，时，分 秒- 卫星钟的偏差（s）- 卫星钟的漂移（s/s）- 卫星钟的漂移速度（s/） I21X,I2.24（1X,I2） F5.13D19.12广播轨道1- IODE（Issue of Data, Ephemer

14、is/数据、星历发布时间）- （rad/s）（rad）3X,4D19.12广播轨道2- e轨道偏心率（radians）- sqrt（A）（广播轨道3- TOE星历的参考时刻（GPS周内的秒数）（rad）（OMEGA）广播轨道4（rad）（OMEGA DOT）广播轨道5- i（rad/s）（IDOT）上的码- GPS周数（与TOE一同表示时间）。为连续计数，不是1024的余数P码数据标记广播轨道6- 卫星精度（m）- 卫星健康状态（第一子帧第三字第1722位）- TGD（sec）- IODC钟的数据龄期广播轨道7- 电文发送时刻（单位为GPS周的秒，通过交接字（HOW）中的Z计数得出）- 拟合区

15、间（h），如未知则为零- 备用3.3 GPS卫星的坐标计算3.3.1用广播星历计算卫星位置第一步：计算卫星运动的平均角速度n首先根据广播星历中给出的参数计算出参考时刻Toe的平均角速度no：式中，GM为万有引力常数G与地球总质量M之乘积，其值为GM3.98600471014m3s2。然后根据广播星历中给定的摄动参数n计算观测时刻卫星的平均角速度n：第二步：计算观测瞬间卫星的平近点角MS由于卫星的运行周期为12小时左右，采用卫星过近地点时刻to来计算平近点角M时，外推间隔最大有可能达6小时。而广播星历每2小时更新一次，将参考时刻设在中央时刻时，外推间隔1小时。所以用Toe来取代卫星过近地点时刻t

16、o后，外推间隔将大大减小，用较简单的模型也能获得精度较高的结果。式中，Mo为参考时刻Toe时的平近点角，由广播星历给出。第三步：计算偏近点角依据（28）式，用弧度表示的开普勒方程为：第四步：计算真近点角根据开普勒轨道方程，可得近点角fs与偏近点角ES之间的关系式中，e为卫星轨道的偏心率，由广播星历给出。由此可得真近点角计算常用公式第五步：计算升交距角为近地点角距，由广播星历给出。第六步：计算摄动改正项广播星历中给出了Cuc、Cus、Crc、Crs、Cic、Cis6个摄动参数，据此可求出由于地球引力场位函数的二阶带谐系数项而引起的升交距角u的摄动改正项、卫星矢径r的摄动改正项r，和卫星轨道倾角i

17、的摄动改正项i。计算公式如下：第七步：对升交距角、卫星矢径r、轨道倾角i进行摄动改正aS为卫星轨道的长半径，a（）2，由广播星历给出。：io为TOe时刻的轨道倾角，由广播星历中的开普勒六参数给出。完为i的变化率，由广播星历中的摄动九参数给出。第八步：计算卫星在轨道面坐标系中的位置在轨道平面直角坐标系中（坐标原点位于地心，X轴指向升交点）卫星的平面直角坐标为：第九步：计算观测瞬间升交点的经度L若参考时刻Toe时升交点的赤经为Toe，升交点对时间的变化率为，那么观测瞬间t的升交点赤经应为：可从广播星历的摄动参数中给出。设本周开始时刻（星期日0时）格林尼治恒星时为GASTweek，则观测瞬间的格林尼

18、治恒星时为：e为地球自转角速度，其值为e7.292115l0-5rads；t为本周内的时间（s）。这样就可求得观测瞬间升交点的经度值为：令则有：注意：广播星历中给出的L并不是参考时刻Toe的升交点赤经Toe，而是该值与本周起始时刻的格林尼治恒星时GASTweek之差。说明：第十步：计算卫星在瞬时地球坐标系中的位置已知升交点的大地经度L以及轨道平面的倾角i后，就可通过两次旋转求得卫星在地固坐标系中的位置：第十一步：计算卫星在协议地球坐标系中的位置观测瞬间卫星在协议地球坐标系中的位置在GPS定位导航中，需要多次连续计算卫星位置和速度，如用上述方法计算需占用较多的内存空间和计算时间。为此，常将卫星星

19、历用一个时间多项式来表示，在内存中仅保存该多项式的系数，供计算时调用。在各种多项式中切贝雪夫多项式的拟合效果最佳，即使在该时间段的两端近似性也很好。用n阶切贝雪夫多项式来逼近时间段t0，t0+t中的卫星星历时，先将变量变换为变量于是卫星坐标可表示为：n为多项式的阶数；Cxi为切贝雪夫多项式的系数。根据已知的卫星坐标，用最小二乘法拟合出多项式系数Cxi后，就可用式（225）计算出该时段中任一时刻的卫星位置。切贝雪夫多项式Ti，的递推公式如下：卫星的运动速度和加速度也可用类似方法计算。3.3.2、用精密星历计算卫星位置精密星历是按一定的时间间隔（通常为15min）来给出卫星在空间的三维坐标、三维运

20、动速度及卫星钟改正数等信息。著名的IGS（国际GPS地球动力学服务局）综合精密星历需12周后才能获得。由IGS提出的格式被广泛采用，其中ASCII格式的SPI和二进制格式的E（F）格式不仅给出了卫星的三维位置信息（km），也给出了卫星的三维运动速度信息（kms）。而SP2（ASCII格式）和ECF2（二进制格式）则仅给出了卫星的三维位置信息。速度信息需通过位置信息用数值微分的方法来求出。采用这种格式时存储量可减少一半左右。在SP3（ASCII格式）和ECF3（二进制格式）中增加了卫星钟的改正数信息。观测瞬间的卫星位置及运动速度可采用内插法求得。其中拉格朗日（Lagrange）多项式内插法被广泛采用，因为这种内插法速度快且易于编程。拉格朗日插值公式十分简单：已知函数yf（x）的n+1个节点x0，xl，x2，xn及其对应的函数值y0，yl，y2，yn对插值区间内任一点x，可用下面的拉格朗日插值多项式来计算函数值：研究表明：对GPS卫星而言，如果要精确至108，用30min的历元间隔和9阶内插已足够保证精度。用40min的历元间隔以及17阶多项式就能获得毫米级（10l0）的精度。例：精密星历sp3数据格式