1、18GPS 卫星的测距码信号和 GPS 卫星的导航电文学习指南在这一章节中，主要讲述了 GPS 全球卫星定位系统的组成，简单地介绍了 GPS卫星信号、GPS 卫星星历和卫星运动理论基础。重点介绍了 GPS 卫星信号特点及其应用。 对本章的学习要重点突出 GPS 系统的组成、 GPS 卫星信号的应用，理解和掌握 GPS 卫星位置计算的各项参数物理意义和几何特点。本单元教学重点和难点1、GPS 编码的方法；2、导航电文的格式和内容。教学目标1、了解 GPS 卫星信号的作用；2、熟悉 GPS 编码的方法；3、熟悉导航电文的格式；4、熟悉导航电文的内容。一、 GPS 卫星信号1 GPS 卫星信号构成及
2、产生GPS 卫星发射的信号由载波、测距码和导航电文三部分组成。如图 25 所示：1.1 载波 L1、 、L 1 由卫星上的原于钟所产生的基准频率 f01.023MHz 倍频 154 倍和 120 倍产生。1.2 测距码1.2.1C/A 码C/A 码又称为粗捕获码，它被调制在 L1 载波上，是 1.023MHz 的伪随机噪声码（PRN 码） ，由卫星上的原子钟所产生的基准频率 f0 降频 10 倍产生，即：f C/Af 0101.023MHz 。由于每颗卫星的 C/A 码都不一样，因此，我们经常用它们的 PRN 号来区分它们。C/A 码是普通用户用以测定测站到卫星间的距离的一种主要的信号。1.2
3、.1 P 码P 码又被称为精码，它被调制在 L1 和 L2 载波上，是 10.23MHz 的伪随机噪声码，直接使用由卫星上的原于钟所产生的基准频率，即：f pf 01.023MHz，其周期为七天。在实施 AS 时，P码与 W 码进行模二相加生成保密的 Y 码，此时，一般用户无法利用 P 码来进行导航定位。1.2.1 L2C 码L2C 码称为城市码，它被调制在 L2 载波上，L2C 信号包括 2 个 PRN 码:即 CM 码和 CL 码。2005 年 9 月 23 日第一颗具有广播 L2C 信号功能的 GPS 卫星 SLC-17A 从 CapeCanaveral, Florida（佛罗里达）发射
4、升空。 L2 码同样可以提供高质量(低相噪，高灵敏度)的数据来进行导航定位。1.3 导航电文19导航信息被调制在 L1 载波上，其信号频率为 50Hz，包含有 GPS 卫星的轨道参数、卫星钟改正数和其它一些系统参数。用户一般需要利用此导航信息来计算某一时刻 GPS 卫星在地球轨道上的位置，导航信息也被称为广播星历。图 25 GPS 卫星信号构成及产生2 GPS 的测距码信号测距码是用于测定从卫星至接收机间距离的二进制码，如图 25 所示。GPS 卫星中所用的测距码从性质上讲属于伪随机噪声码。它们看似一组取值(0 或 1)完全无规律的随机噪声码序列，其实是具有确定编码规则编排起来的、可以复制的周
5、期性的二进制序列，且具有类似于随机噪声码的自相关性特性，结构相同的随机码序列 通过平移码元数，相应的码元相互对齐，(t)u易于测量。测距码是由若干个多级反馈移位寄存器所产生的 m 序列经平移、截短、求模二和等一系列复杂处理后形成的。根据性质和用途的不同，在 GPS 卫星发射的测距码信号中包含了CA 和 P（Y ）码两种伪随机噪声码信号，各卫星所用的测距码互不相同。下面将分别介绍其特点及作用。图 2-6 CA 码、P 码的特点2.1 CA 码(CoarseAcquisition Code)用于进行粗略测距和捕获精码的测距码称为粗码，也称捕获码。CA 码的测距精度一段为(2-3) m。C A 码是
6、一种结构公开的明码，供全世界所有的用户免费使用。CA 码的特征是：码长 Nu2 10-1=1023bit；码元宽度 tu0.97752us，相应长度 293.1m;周期 Tu Nutu=1ms；数码率 BPS1.023Mbit s ，如图 2-7 所示。GPS 星座中的不同卫星使用结构各异的 CA 码。这样既便于复制又易于区分。CA 码具有的特性：由于 CA 码的码长较短（周期1ms） ，在 GPS 导航和定位中，为了捕获 CA 码以测定卫星信号传播的时间延迟，通常对 CA 码金行逐个搜索，而 CA 码总共只有 1023 个码元，若以每秒 50 码元的速度搜索，仅需约 20.5s 便可完成，易
7、于捕获。而通过捕获 CA 码所得到的20卫星提供的导航电文信息，又可以方便地捕获 P 码，所以，通常称 CA 为捕获码。CA 码的码元宽度 t 0=1/f=0.97752s(码元持续时间 )，空间矢距：L=293.1m( 码元持续时间 与 c 乘积 )较大。若两个序列的码元相关误差为码元宽度的 l101100，则此时所对应的测距误差可达 29.32.9m。由于其精度较低，所以称 CA 码为粗捕获码。2.2 P 码(Precision Code)用于精确测定从 GPS 卫星至接收机距离的测距码称为精码。该测距码又同时调制在 L 1 和 L2两个载波上，可较完善地消除电离层延迟，故用它来测距可获得
8、较精确的结果。P 码是一种结构保密的军用码。目前，美国政府不提供一般 GPS 民用用户使用。P 码的特征是：码长 Nu=2.351014bit；码元宽度 tu=0.097752us，相应长度 29.3m ，数码率BPS10.23Mbit s ；如图 2-7 所示, 周期 TuNu tu=267 天，一个周期中约含 6.2 万亿个码元。实际上 P 码的一个整周期被分为 38 部分，每一部分周期为 7 天，码长约 6.19x1012bit。其中有 5部分由地面监控站使用，其他 32 部分分配给不同的卫星，1 个部分闲置。这样，每颗卫星所使用P 码便具有不同的结构，易于区分，但码长和周期相同。P 码
9、具有的特性：因为 P 码的码长较长（6.1910 12bit） ，在 GPS 导航和定位中，如果采用搜索 CA 码的办法来捕获 P 码，即逐个码元依次进行搜索，当搜索的速度仍为每秒 50 码元时，约需 14155天，那将是无法实现的，不易捕获。因此，一般都是先捕获 CA 码，然后根据导航电文中给出的有关信息，便可捕获 P 码。 P 码的码元宽度 t 0=1/f=0.097752s，每个码元所持续的时间为 CA 码的 110。空间矢距：L=29.3m( 码元持续时间 与 c 乘积 )较大。若两个序列的码元相关误差仍为码元宽度的l101100，则此时所引起的测距误差仅有 2.930.293m。仅为
10、 CA 码的 110。所以 P 码可用于较精密的导航和定位，称为精码。2.3 L2C 码目前，CA 码只调制在 Ll 载波上，故无法精确地消除电离层延迟。随着全球定位系统的现代化，在卫星上增设调制了 CA 码的第二民用频率码 L2C 码后，该问题将可得到解决。采用窄相关间隔(Narrow Correlator Spacing) 技术后测距精度可达分米级，与精码的测距精度大体相当。3 CPS 卫星的载波信号可运载调制信号的高频振荡波称为载波。GPS 卫星发射两种频率的载波信号，由于它们均位于微波的 L 波段，故分别称为 L 1 载波和 L2 载波。即频率为 1575.42MHz 的 L1 载波和
11、频率为1227.60HMz 的 L2 载波，它们的波长分别为 19.03cm 和 24.42cm，如图 2-6 所示。在无线电通信技术中，为了有效地高质量传播信息，都是将频率较低的信号加载在频率较高的载波上，此过程称为调制。GPS 卫星的 L1 和 L2 载波上携带着测距信号和导航电文传送出去，到达用户接收机。在一般的通信中，当调制波到达用户接收机解调出有用信息后，载波的作用便告完成。但在全球定位系统中，载波除了能更好地传送测距码和导航电文这些有用信息外，在载波相位测量中它又被当做一种测距信号来使用。其测距精度比伪距测量的精度高 23 个数量级。因此，载波相位测量在高精度定位中得到了广泛的应用
12、。采用两个不同频率载波的主要目的是为了较完善地消除电离层延迟。采用高频率载波的目的是为了更精确地测定多普勒频移和载波相位(对应的距离值 )，从而提高测速和定位的精度，减少信号的电离层延迟，因为电离层延迟与信号频率 f 的平方成反比。4 GPS 导航电文GPS 卫星导航电文是用户利用 GPS 定位和导航所必须的基础数据。它主要提供了卫星在空间的位置、卫星的工作状态、卫星钟的修正参数、电离层延迟修正参数等重要信息。这些信息是以二进制码的形式按规定格式编码，并按帧发给用户接收机，因此又称之为数据码(D 码) 。导航电文的传输速率为 50bits，以“帧”为单位向外发送。每帧的长度为 1500bit，
13、播发完21一个主帧需 30s。一个主帧包括 5 个子帧。每个子帧均包含 300bit，播发时间为 6s。每个子帧又可分为 10 个字，每个字都由 30bit 组成。其中第四、五两个子帧各有 25 个页面，需要 750s 才能将 25 个页面全部播发完。第一、 二、三子帧每 30s 重复一次，其内容每隔 2h 更新一次。第四、五子帧每 30s 翻转一页，12.5min 完整地播发一次，然后再重复。其内容仅在卫星注入新的导航数据后才得以更新。卫星电文的基本构成如图 27 所示。图 27 导航电文的组成在每帧导航电文中，各子帧电文的主要内容如图 27 所示，下面介绍电文各部分的基本含义。4.1 遥测
14、码(telemetry word，TLW)遥测码位于各子帧的开头，它用来表明卫星注入数据的状态。遥测码的 l8bit 是同步码(10001001)，为各子帧编码脉冲提供一个同步起点，接收机将从该起点开始顺序解译电文。第922bit 为遥测电文，包括地面监控系统注入数据时的状态信息、诊断信息及其他信息。第 23bit和第 24bit 是连接码，第 2530bit 为奇偶检查码，它用于发现和纠正错误。4.2 转换码(hand over word，HOW)转换码位于每个子帧的第二个字码。其作用是提供用户从捕获的 CA 码转换到捕获 P 码的Z 计数。 Z 计数值位于转换码的第 l17bit，是从每周
15、六周日零时起算的时间计数。因此，当知道了 Z 计数，即知道了观测瞬间在 P 码周期中所处的准确位置，这样便可迅速捕获 P 码。4.3 第一数据块第一数据块是位于第 1 子帧的第 310 字码，它的主要内容包括： (1)时延差改正 Tgd时延差改正 Tgd 就是载波 L1、L 2 的电离层时延差。当使用单频接收机时，为了减小电离层效应影响，提高定位精度，要用 Tgd 改正观测结果；双频接收机可通过 L1、L 2 两项频率的组合来消除电离层效应的影响，不需要此项改正。(2)数据龄期 AODC卫星时钟的数据龄期 AODC 是时钟改正数的外推时间间隔，它指明卫星时钟改正数的置信度。AODCt 0cti
16、 (21)式中，t 0c 为数据块 I 的参考时刻；t i 是计算时钟改正参数所用数据的最后观测时间。(3)星期序号 WN10.sec.611 1542512.min30sec2311 98765432 10secTLWHO第一数据块修正参数第二数据块星历表 第三数据块卫星历书(1500bit) (37500bit)22WN 表示从 1980 年 1 月 6 日子夜零点(UTC)起算的星期数，即 GPS 星期数。(4)卫星时钟改正GPS 时间系统是以地面主控站的主原子钟为基准。由于主控站主原子钟的不稳定性，使得GPS 时间和 UTC 时间之间存在差值。地面监控通过监测确定出这种差值，并用导航电
17、文播发给广大用户。GPS 卫星的时钟相对 GPS 时间系统存在着差值，需加以改正，这便是卫星时钟改正。(22)202010s )()(t cctt式中， 0 为卫星钟差(s)； 1 为卫星钟速(ss) ； 2 为卫星钟速变率(ss 2)4.4 第二数据块导航电文的第 2 和第 3 子帧组成第二数据块，其内容为 GPS 卫星星历，即描述卫星运行及其轨道参数的信息，提供有关计算卫星运行位置的数据，它是 GPS 卫星向导航、定位用户播发的主要电文，描述卫星的运行及其轨道参数包括以下三类：如图 28 所示，图 28GPS 卫星轨道参数(1)开普勒轨道六参数： 为卫星轨道椭圆长半轴的平方根； e 为卫星
18、轨道椭圆偏心率；i 0 为参考时刻 t0 的轨道面倾角； 0 为参考时刻 t0 的升交点赤经； 为近地点角距；M 0 为参考时刻 t0 的平近点角。(2)轨道摄动九参数： n 为平均角速度改正数 ; 为升交点赤经变化率;i 为卫星轨道平面倾角变化23率；C us、C uc 为升交角距的正余弦调和改正项振幅，C is、C ic 为轨道正面倾角的正余弦调和改正项振幅；C rs、C rc 为轨道向径正余弦调和改正项振幅。(3)时间二参数：从星期日子夜零点开始度量的星历参考时刻 t0e 及星历表的数据龄期 AODE。4.5 第三数据块第三数据块包括第 4 和第 5 两个子帧，其内容包括了所有 GPS
19、卫星的历书数据。当接收机捕获到某颗 GPS 卫星信号后，根据第三数据块提供的其他卫星的概略星历、时钟改正、卫星工作状态等数据，用户可以选择工作正常、位置适当的卫星，并较快地捕获到所选择的卫星。二、GPS 卫星星历1 卫星运动理论基础 GPS 卫星的星历是描述卫星运行及其轨道的参数，它的主要作用是利用 GPS 卫星系统进行导航定位时，计算卫星在空间的瞬时位置。而研究 GPS 卫星在协议地球坐标系中的瞬时位置，就是 GPS 卫星的轨道运动理论。本节首先对卫星运动理论作简要介绍。 卫星在空间运行的轨迹称为轨道，描述卫星位置及状态的参数，称为卫星轨道参数，而轨道参数取决于卫星所受到的各种力的作用。众所
20、周知，人造地球卫星在空间运行时，除了受地球重力场的引力作用外，还受到太阳、月亮及其他天体引力的影响，同时还受到大气的阻力、太阳光压力及地球潮汐的作用力等因素的影响。为了研究卫星运动的基本规律，一般将卫星受到的作用力分为两类：一类是地球质心引力，即将地球看作密度均匀并由无限多的同心球层所构成的圆球，它对球外一点的引力等效于质量集中于球心的质点所产生的引力，称为中心引力。另一类是摄动 力，也称非中心引力，它包括地球非球形对称的作用力、日月引力、大气阻力、光辐射压力及地球潮汐作用力。1.1 卫星无摄运动摄动力与中心引力相比，仅为 10-3 量级。所以，人造地球卫星在空间运行时，主要受地球重力场的引力
21、作用。它决定着卫星运动的基本规律和特征，由它所决定的卫星轨道是研究卫星实际轨道的基础。卫星在预定的轨道上运行，如果忽略摄动力的影响，地球可视为质量全部集中于质心的质点，卫星也可以看作是质量集中的质点。根据万有引力定律，地球受卫星的引力可表示为(23)rmGM2eF式中：M 为地球质量；m 为卫星质量；G6672x10 -8cm3(gs 2)为万有引力常数； 为r卫星在( 历元) 平天球坐标系中的位置向量； 为向量 的模，即卫星至地球的距离。rr根据牛顿第三定律，卫星受地球的引力 ，其大小与 相同而方向相反，即sFe(24) rmM2s根据牛顿第二定律，可得卫星及地球的运动方程为rGdt2224
22、(25)rmGMdtr22由此，可得卫星在无摄运动理想条件下，卫星相对地球的运动方程为(26)rtr22卫星在上述地球引力场中的无摄运动称为开普勒运动, 其规律可用开普勒定律来描述。由开普勒定律可知，卫星运动的轨道，是通过地心平面上的一个椭圆，且椭圆的一个焦点与地心相重合。而确定椭圆的形状和大小至少需要两个参数，即椭圆的长半径 as 及其偏心率es(或椭圆的短半径 bs)。另外，为确定任意时刻卫星在轨道上的位置，需要一个参数，一般取真近点角 fs，即在轨道平面上，卫星与近地点之间的地心角距，该参数为时间的函数，它确定了卫星在轨道上的瞬时位置。参数 as，e s，f s 唯一地确定了卫星轨道的形
23、状、大小以及卫星在轨道上的瞬时位置。如图 2-9 所示，称之为轨道椭圆形状参数。 但是，这时卫星轨道平面与地球体的相对位置和方向还无法确定。要确定卫星轨道与地球体之间的相互关系，亦可表达为确定开普勒椭圆在天球坐标系中的位置和方向。因为根据开普勒第一定律，轨道椭圆的一个焦点与地球的质心相重合，所以为了确定该椭圆在上述坐标系中的方向，尚需三个参数。 这三个参数的选择并不是唯一的。其中一组应用广泛的参数，称为开普勒轨道参数(图 29)，或称开普勒轨道根数。现将这组参数的惯用符号及其定义，综合介绍如下：升交点的赤经，即在地球赤道平面上，升交点与春分点之间的地心夹角(升交点，即当卫星由南向北运行时轨道与
24、地球赤道面的一个交点)。i轨道面的倾角，即卫星轨道平面与地球赤道面之间的夹角。上两个参数，唯一地确定了卫星轨道平面与地球体之间的相对定向，称之为轨道平面定向参数。s近地点角距，即在轨道平面上，升交点与近地点之间的地心夹角，这一参数表达了开普勒椭圆在轨道面上的定向，称之为轨道椭圆定向参数。在此，参数 as、e s、i、 s 和 fs 所构成的坐标系统，通常称为轨道坐标系统。其中，参数as、 es、 、i、 s 的大小，则是由卫星的发射条件决定。在该系统中，当 6 个轨道参数一经确定后，卫星在任一瞬间相对地球体的空间位置及其速度，便可唯一地确定。建立以地球质心为坐标原点，x 轴指向近地点，y 轴重
25、合于轨道的短轴，z 轴为轨道平面的图 210 开普勒轨道参数ases25法线方向，构成右手坐标系。在此坐标系内列出卫星运动的微分方程并求解，可以得出著名的开普勒轨道方程：(27)SS0sinEe-)t-n(式中， n 为观测时刻卫星的平均角速度。设卫星沿椭圆轨道运动的周期为 T，则 n=2/T。t0 为第六个积分常数，即 t0 是卫星过近地点时刻，它给出了辅助参数 E 与时间 t 的函数关系。ES 为偏近点角。如图 210 所示，假设过卫星质心 ms，作平行于椭圆短半袖的直线，则 为该直线与近地点至椭圆中心连线的交点， 为该直线与以椭圆中心为原点，并以 as 为m 半径的大圆的交点。E S 就
26、是在椭圆平面上，近地点 P 至 点的圆弧所对应的圆心角。由开普勒轨道方程知，当 t t0 时，E S0。顾及轨道方程式，可得 r=aS(1eS)。说明此时卫星正位于近地点处。从而证明了 t0 是卫星过近地点的时刻。令 MSn(tt 0)，则 M 随时间 t 以平均角速度 n 变化，故称 MS 为平近点角。又令 M0nt 0为过近地点的平近点角，则(28)00S-t)-(t开普勒轨道方程可写为：(29)SSsinEe-至此，我们得到了以轨道参数表示的六个积分常数参数 as、e s、i、 s 和 t0。若已知六个轨道参数，就可以唯一地确定卫星的运动状态。也就是说，已知六个轨道参数可以确定任意时刻的
27、卫星位置及其运动速度。1.2 . 卫星受摄运动在摄动力作用下，卫星的运动称为受摄运动，相应的卫星轨道称为受摄轨道。摄动力作用使卫星的运动产生一些小的附加变化而偏离理想轨道，同时，这种偏离量的大小也随时间而变化着。对于卫星精密定位来说，在只考虑地球质心引力情况下计算卫星的运动状态(即研究二体问题)是不能满足精度要求的。必须考虑地球引力场摄动力、日月摄动力、大气阻力、光压摄动力、潮汐摄动力对卫星运动状态的影响。讨论卫星相对于地球无摄运动的二体问题时，六个轨道参数均为常数。其中卫星过近地点的时刻 t0 也可用平近点角 M0 代替。在考虑了摄动力的作用后，卫星的受摄运动的轨道参数不再保持为常数，而是随
28、时间变化的轨道参数。卫星在地球质心引力和各种摄动力总的影响下的轨道参数称为瞬时轨道参数。卫星运动的真实轨道称为卫星的摄动轨道或瞬时轨道。瞬时轨道不是椭圆，轨道平面在空间的方向也不是固定不变的。 图 210 偏近点角与真近点角26在人造地球卫星所受的摄动力中，地球引力场摄动力最大，约为 10-3 量级，其他摄动力大多小于或近于是 106 量级。这些摄动力引起卫星位置的变化，引起轨道参数的变化。例如，考虑地球引力场摄动力中地球引力场位函数的二阶带谐系数项的影响，使轨道参数不断减小，即轨道平面不断西退，这种现象称为轨道面的进动。进动速度主要取决于轨道倾角 i 和轨道长半径as。对于两万公里高度，倾角
29、约为 55 度的 GPS 卫星来说，其进动速度约为 0.039 度天。轨道参数 的变化使得近地点在轨道面内不断旋转，或者说轨道椭圆以其不变的形状在轨道面内旋转。通过解算卫星受摄运动的微分方程，可以得到卫星轨道参数的变化规律。1.2 .1 用直角坐标表示的受摄运动方程在直角坐标系中，卫星的受摄运动方程形式简单。设作用于卫星上的摄动力位函数为 R，则受摄运动方程的分量形式可写为：(210) zRyxx/z)/r-(zyy/)/r-(x33式中 ， ， 分别为卫星在地球质心引力作用下产x)/r(-3)/r(-3生的加速度沿三个坐标轴的分量。这种形式的微分方程在求解的过程中不涉及卫星的轨道参数，可以用
30、数值方法求解。1.2 .2 用轨道参数表示的受摄运动方程用直角坐标表示的受摄运动方程难以得到关于卫星的运动轨道及其变化规律。而以轨道参数表示的受摄运动方程则既可以用于数值解法也可用于分析解法。其中如果摄动力的性质为非保守力时，例如太阳辐射压力、大气阻力因不存在位函数，具有代表性的卫星受摄运动方程是牛顿受摄运动方程。(211) )sin)1()2(cos1dt dtcositsin)co(dt1sti )cos(indte )1(si12ta2S222ss TfprSrefenaMiflfWeafrf TfESfefessssssssss式中，S 为沿卫星矢径 方向的分量，T 为在轨道平面上垂直
31、于矢径方向并指向卫星运动的分量，W 为沿轨道平面法线并按 S，T，W 组成右手坐标系取向的分量。此时，可将摄动力所产生的加速度分解为互相垂直的三个分量 S，T，W。不论摄动力的性质如何，都可以使用牛顿受摄运动方程解卫星的受摄运动。通过研究牛顿受摄运动方程可知，由于卫星在运动中受到各种摄动力作用的影响，其轨道参数随时间而变化。若已知某一初始时刻的轨道参数，通过分析解算含有轨道参数的受摄运动方程，可以求得轨道参数的变率，从而求得任一时刻的轨道参数。这样，利用二体问题的运动方程就可以求得任一时刻的27卫星位置和速度。2 GPS 卫星星历描述某一时刻卫星运动轨道的参数及其变率称卫星星历。根据卫星星历就
32、可以计算出任一时刻的卫星位置及其速度。GPS 卫星定位中，需要知道 GPS 卫星的位置。通过卫星的导航电文将已知的某一初始历元的轨道参数及其变率发给用户(接收机) ，即可计算出任一时刻的卫星位置。另外，通过在已知的地面站对 GPS 卫星进行观测，求得卫星在某一时刻的位置，可以反求出卫星的轨道参数，从而对卫星的轨道进行改进，实现精密定轨，用于 GPS 精密定位。因此，精确的轨道信息是精密定位的基础。GPS 卫星星历分为预报星历(广播星历)和后处理星历( 精密星历)。2.1 预报星历是通过卫星发射的含有轨道信息的导航电文传递给用户的，用户接收机接收到这些信号，经过解码便可获得所需要的卫星星历，所以
33、这种星历也叫广播历。卫星的预报星历，通常包括相对某一参考历元的开普勒轨道参数和必要的轨道摄动正项参数。参考历元的卫星开普勒轨道参数，也叫参考星历，它是根据 GPS 监测站约一周的观测资料推算的。参考星历只代表卫星在参考历元的瞬时轨道参数，但是在摄动力的影响下，卫星的实际轨道随后将偏离其参考轨道，偏离的程度主要决定于观测历元与所选参考历元间的时间差。一般来说，如果我们用轨道参数的摄动项对已知的卫星参考星历加以改正，就可以外推出任意观测历元的卫星星历。为了保持卫星预报星历的必要精度，一般采用限制预报星历外推时间间隔的方法。为此，GPS 跟踪站每天都利用其观测资料，更新用以确定卫星参考星历的数据，以
34、计算每天卫星轨道参数的更新值，并且，每天按时将其注入相应的卫星加以储存，以更新卫星的参考轨道之用。据此，GPS 卫星发射的广播星历，每小时更新一次，以供用户使用。这样，如果将上述计算参考星历的参考历元 toe 选在两次更新星历的中央时刻，则外推的时间间隔最大将不会超过 0.5 小时。从而可以在采用同样摄动力模型的情况下。有效地保持外推轨道参数的精度。预报星历的精度，一般约为 2040m。预报星历的内容包括：参考历元瞬间的开普勒 6 个参数，反映摄动力影响的 9 个参数，以及1 个参考时刻和星历数据龄期，共计 17 个星历参数。这些参数通过 GPS 卫星发射的含有轨道信息的导航电文传递给用户。
35、表 2-1 GPS 卫星广播星历预报参数及其定义参数 参数定义Toe（s） 星历表参考历元 (秒) IODE(AODE) 星历表数据量(N) M0（rad ） 按参考历元 toe 计算的平近点角(弧度) n（rad/s） 由精密星历计算得到的卫星平均角速度与按给定参数计算所得的平均角速度之差(弧度) e 轨道偏心率 （m）a轨道长半径的平方根(0.5m)0（rad） 按参考历元 toe 计算的升交点赤径(弧度) i0（rad ） 按参考历元知计算的轨道倾角 (弧度) （rad） 近地点角距 (弧度) （rad/s） 升交点赤径变化率( 弧度/ 秒) I（rad/s） 轨道倾角变化率( 弧度/
36、秒) Cuc（rad ） 升交距角 (+fS)的余弦调和项改正的振幅(弧度) Cus（rad ） 升交距角 (+fS)的正弦调和项改正的振幅(弧度) 28Crc（m） 轨道半径的余弦调和项改正的振幅 (米) Crs（m） 轨道半径的正弦调和项改正的振幅 (米) Cic（rad ） 轨道倾角的余弦调和项改正的振幅 (弧度) Cis（rad ） 轨道倾角的正弦调和项改正的振幅 (弧度) GPD 周数( 周)Tgd 电离层延迟改正 (秒)，IODC 星钟的数据量(N)a0 卫星钟差 时间偏差( 秒)a1 卫星钟速(秒/ 秒) 频率偏差系数a2 卫星钟速变率( 秒/ 秒 2)漂移系数卫星精度 (N)卫
37、星健康 (N)其中n 中包括了轨道根数 的常期摄动。n 中主要是二阶带谐项引起的的长期漂移，也包括了日、月引力摄动和太阳光压摄动。在中主要是二阶带谐项引起升交点赤径的长期漂移，也包括了极移的影响。星历表参考历元 toe 是从星期日子夜零点开始计算的参考时刻，星历表数据龄期 IODE 为从toe 时刻至为作预报星历测量的最后观测时刻之间的时间，故 IODE 是预报星历的外推时间间隔。为了对上述 3 类参数有一个具体的数据概念，表 2-2 列出一组用 Novatel RT2 双频接收机测量，并由其随机软件提供的 GPS 卫星广播星历：( 时间：03-Sep-00 11:16)。表 2-3 摘录了1
38、990 年 3 月 16 日用 Trimble 4000SST 双频接收机观测 3 颗卫星的轨道参数。表 2-2 Novatel RT2 双频接收机接收的 GPS 卫星广播星历2 NAVIGATION DATA RINEX VERSION / TYPESSUtilities v1.41 Your Company 03-Sep-00 11:16 PGM / RUN BY / DATE END OF HEADER 6 0 8 25 2 0 0.0 +.274740159512D-07 +.113686837722D-12 .000000000000D+00+.950000000000D+02 +.
39、888437500000D+02 +.454554648312D-08 -.170024978378D+01+.444985926151D-05 +.697577721439D-02 +.102780759335D-04 +.515367507172D+04+.439200000000D+06 +.838190317154D-07 -.717044767358D+00 -.707805156708D-07+.947643365081D+00 +.170531250000D+03 -.237102754179D+01 -.796068873750D-08-.203937066226D-09 .0
40、00000000000D+00 +.107600000000D+04 .000000000000D+00.000000000000D+00 .000000000000D+00 -.512227416039D-08 +.351000000000D+03.000000000000D+004 0 8 25 2 0 0.0 +.359198544174D-03 +.355839802069D-10 .000000000000D+00+.183000000000D+03 +.680000000000D+02 +.427839249795D-08 -.257055307047D+01+.357255339
41、622D-05 +.518933811691D-02 +.323541462421D-05 +.515354366875D+04+.439200000000D+06 -.465661287308D-07 +.346417413904D+00 +.186264514923D-08+.976317339170D+00 +.328218750000D+03 -.538595653048D+00 -.819355558019D-08+.137148569931D-09 .000000000000D+00 +.107600000000D+04 .000000000000D+00.000000000000
42、D+00 .000000000000D+00 -.605359673500D-08 +.439000000000D+03表 2-3 GPS 卫星广播星历轨道参数值CPS 卫星参数PRN02 PRN14 PRN19Toe（s） 5.184000E 十 005 5.184000E 十 005 5.184000E 十 005IODE（s） 2.120000E 十 002 7.500000E 十 001 1.980000E 十 002M0（rad ） 6.054760E001 2.098430E 一 001 7.666470E 一 00ln（rad/s） 1.439048E009 1.375952E0
43、09 1.514195E 一 009e 8.247306E 一 003 6.024992E003 3.710157E 一 00329（m）a5.153691E 十 003 5.153689E 十 003 5.153629E 十 0030（rad） 8.267319E001 1.815057E 一 001 8.471960E 一 001i0（rad ） 3.041002E 一 00l 3.060368E 一 001 3.038775E00l（rad） 9.525346E001 8.729396E00l 2.284949E001（rad/s） 一 2.569777E009 2.513957E 一
44、009 2.654474E 一 009I（rad/s） 2.273737E 一 012 2.023625E011 3.853984E011Cuc（rad ） 9.166209E007 5.02l 849E007 6.23l 362E 一 007Cus（rad ） 2.530490E006 3.095522E005 1.117021E 一 006Crc（rad ） 2.213125E 十 002 1.932500E 十 002 3.065000E 十 002Crs（rad ） 一 5.659375E 十 001 2.828125E 十 001 3.640625E 十 001Cic（m） 1.89
45、7275E 一 008 4.446738E 一 008 3.498100E 一 008Cis（m） 一 6.225433E 一 008 6.521882E 一 009 1.363666E008GPS 卫星向全球用户播发的星历，是用两种波码进行传送的。一种是用叫做 CA 码所传送的 GPS 卫星星历(简称 C A 码星历)，其星历精度为数十米。另一种用 P 码所传送的 GPS 卫星星历( 简称 P 码星历)精度提高到 5m 左右，只有工作于 P 码的接收机才能从 P 码中解译出精密的 P 码星历。精密的 P 码星历主要用于军事目的导航定位。C A 码星历交付民用。目前绝大多数的商品接收机，都是工
46、作于 CA 码的，只能使用降低了精度的 CA 码星历。CA 码星历精度的人为降低，给用户的 GPS 定位引入相应误差。这是非特许用户进行高精度的 GPS 测量时必须解决的一个问题。利用精密的后处理星历能够解决这一问题。 2.2 后处理星历由于 GPS 卫星的预报星历是根据跟踪站前一段时间的观测资料，外推的参考轨道参数，并加入轨道的摄动改正后得到的外推星历。因此，广播星历包含外推误差其精度必然受到限制，不能满足某些从事精密定位工作的用户要求。后处理星历，是根据地面跟踪站所获得的精密观测资料计算而得到的星历，它是一种不包含外推误差的实测星历，可为用户提供观测时刻的卫星精密星历，其精度可达米级，以后
47、其精度有望进一步提高到分米级。这种星历不是通过 GPS 卫星的导航电文向用户传递，一些国家某些部门，根据各自建立的卫星跟踪站所获得的对 GPS 卫星的精密观测资料。应用与确定广播星历相似的方法而计算卫星星历。然而，这种星历用户无法实时通过卫星信号而获得，而是利用磁带或通过电视、电传、卫星通讯等方式在事后有偿地向用户提供所需要的服务。3 GPS 卫星位置的计算3.1 用广播星历计算卫星位置第一步：计算卫星运动的平均角速度 n首先根据广播星历中给出的参数 计算出参考时刻 Toe 的平均角速度 no：a(212)30)(GM式中，GM 为万有引力常数 G 与地球总质量 M 之乘积，其值为 GM3.9
48、86004710 14m3s 2。然后根据广播星历中给定的摄动参数n 计算观测时刻卫星的平均角速度 n：(213)n0第二步：计算观测瞬间卫星的平近点角 M由于卫星的运行周期为 12 小时左右，采用卫星过近地点时刻 to 来计算平近点角 M 时，外推30间隔最大有可能达 6 小时。而广播星历每 2 小时更新一次，将参考时刻设在中央时刻时，外推间隔1 小时。所以用 Toe 来取代卫星过近地点时刻 to 后，外推间隔将大大减小，用较简单的模型也能获得精度较高的结果。(214)Toe)-n(tM0S式中，Mo 为参考时刻 Toe 时的平近点角，由广播星历给出。第三步： 计算偏近点角依据(28)式，用弧度表示的开普勒方程为：(215)SSsinEeE解上述方程可用迭代法或微分改正法。第四步： 计算真近点角根据开普勒轨道方程，可得近点
