05、卫星信号

1、第四章 电磁波的传播与GPS卫 星信号 GPS定位的基本观测量是观测站（用户 接收天线）至GPS卫星（信号发射天线） 的距离（或称信号传播路径），它是通 过测定卫星信号在该路径上的传播时间 （时间延迟）或测定卫星载波信号相位 在该路径上的变化周数（相位延迟）来 导出的。 tv 1.电磁波的传播速度与大气折射 假设电磁波在真空中的传播速度为cvac，则有 cvac= vacf= vac/T= /kvac。在卫星大地测量中， 国际上当前采用的真空光速为 c=2.99782458108(m/s)。 对GPS而言，卫星发射信号传播到接收机天线的 时间约0.1秒，当光速值的最后一位含有一个 单位的误差，

2、将会引起0.1m的距离误差。表明 准确确定电磁波传播速度的重要意义。实际的 电磁波传播是在大气介质中，在到达地面接收 机前要穿过性质、状态各异且不稳定的若干大 气层，这些因素可能改变电磁波传播的方向、 速度和强度，这种现象称为 大气折射对GPS观测结果的影响，往往超过了 GPS精密定位所容许的精度范围。 电磁波在大气中的传播速度可以用折射率n来表 示，n=c/v。折射率与大气的组成和结构密切 相关，其实际值接近于1，故常用折射数N0来 表示，N0=(n-1)106。 根据大气物理学，如果电磁波在某种介质中的 ，则该介质成为。 介质的弥散现象是由于传播介质的内电场和入 射波的外电场之间的电磁转换

3、效应而产生的。 当介质的原子频率与入射波的频率接近一致时， 将发生共振，由此而影响电磁波的传播速度。 通常称。如果把具有不同频 率的多种波叠加，所形成的复合波称为群波， 则在具有速度弥散现象的介质中，单一频率正 弦波的传播与群波的传播是不同的。 假设单一正弦波的相位传播速度为相速vp，群波 的传播速度为群速vg，则有 式中为通过大气层的电磁波波长。 若取通过大气层的电磁波频率为f，则相应的折射 率为 在GPS定位中， p pg v vv f n fnn p pg 1.大气层的结构与性质 地球表面被一层很厚的大气所包围，大气的总质 量约为3.9 1018（kg），约为地球总质量的百 万分之一。由

4、于地球引力的作用，大气质量在 垂直方向上分布极不均匀，主要集中在大气底 部，其中75%的质量分布在10km以下， 90% 的以上质量分布在30km以下。同时大气在垂直 方向上的物理性质差异也很大，根据温度、成 分和荷电等物理性质的不同，大气可分为性质 各异的若干大气层。按不同标准有不同的分层 方法，根据对电磁波传播的不同影响，一般分 为对流层和电离层。 是指从地面向上约40km范围内的大气底 层，占整个大气质量的99%。对流层与地面接 触，从地面得到辐射热能，温度随高度的上升 而降低，平均每升高1km降低6.50C，而在水平 方向（南北方向）上，温差每100km 一般不超 过10C。对流层虽仅

5、有少量带电离子，但却具 有很强的对流作用，云、雾、雨、雪、风等主 要天气现象均出现其中。该层大气中除了含有 各种气体元素外，还含水滴、冰晶和尘埃等杂 质，对电磁波的传播有很大影响。 分布于地球大气层的顶部，约在地面向 上70km以上范围。由于原子氧吸收了太阳紫 外线的能量，该大气层的温度随高度上升而迅 速升高，同时由于太阳和其它天体的各种射线 作用，使大部分大气分子发生电离，具有密度 较高的带电粒子。电离层中电子的密度决定于 太阳辐射强度和大气密度，因而导致电离层的 电子密度不仅随高度而异，而且与太阳黑子的 活动密切相关。电磁波在电离层中的传播速度 与频率相关，电离层属于 l折射数随高度的变化

6、。 某一瞬间全球电子密度与测站间的关系 2.对流层的影响与改正 在对流层中，折射率略大于1，随着高度的增加 逐渐减小，当接近对流层顶部时，其值接近于 1。对流层的折射影响，在天顶方向（高度角 900）可产生2.3m的电磁波传播路径误差，当 高度角为100时，传播路径误差可达20m。在 精密定位中，对流层的影响必须顾及。 通常将对流层的大气折射分为两 部分，Nd和Nw分别表示干、湿分量的折射数， 则N0= Nd+Nw。 Nd和Nw与大气的压力、温度和 湿度有如下近似关系 (A) 式中P为大气压力（mbar），Tk为绝对温度 （Tk=0C+273.2）,e0为水汽分压（mbar）。沿 天顶方向，对

7、流层大气对电磁波传播路径的影 响，可表示为 2 0 5 1073. 3 6 .77 k w k d T e N T P N wd SSS 干分量引起的电磁波传播路径距离差主要与地面 的大气压力和温度有关；湿分量引起的电磁波 传播路径距离差主要与传播路径上的大气状况 密切相关。 由地球表面向上沿天顶方向的电磁波传播路径为 考虑干、湿分量的折射数，则有 S0为电磁波在真空中的传播路径，Hd为当Nd趋近 于0时的高程值（约40km）， Hw为当Nw趋近 于0时的高程值（约10km）. HH dHNSndHS 06 0 10 wd H w H d dHNdHNSS 66 0 1010 于是沿天顶方向电

8、磁波传播路径的距离差 为 dHNS dHNS SSSSS w d H ww H dd wd 6 6 0 10 10 在卫星大地测量中，不可能沿电磁波传播路线直 接测定对流层的折射数，一般可以根据地面的 气象数据来描述折射数与高程的关系。 根据理论分析，折射数的干分量与高程H的关系 为 Nd0为按前（A）式计算的地面大气折射数的干 分量，对于参数 Hd，H. Hopfield通过分析全球 高空气象探测资料，推荐了如下经验公式。 4 0 d d dd H HH NN )16.273(72.14840136 kd TH 由于大气湿度随地理纬度、季节和大气状况而变 化，尚难以建立折射数湿分量的理论模型

9、，一 般采用与干分量相似的表示方法 式中Nw0为按（A）式计算的地面大气折射数的湿 分量，高程的平均值取为Hw=11000m 4 0 w w ww H HH NN 积分可得沿天顶方向对流层对电磁波传播路径影 响的近似关系： 数字分析表明，在大气的正常状态下，沿天顶方 向，折射数干分量对电磁波传播路径的影响约 为2.3m，约占天顶方向距离总误差的90%，湿 分量的影响远较干分量影响小。 实际观测时，观测站接收的卫星信号往往不是来 自天顶，此时在考虑对流层影响时必须顾及电 磁波传播方向的高度角。 w k w d k d H T e S H T P S 2 0 5 5 4810 10552. 1 1

10、0552. 1 假设GPS卫星相对观测站的高度角为hs，可得 实践表明，上式中含有较大的模型误差，当hs大 于100时，改正量的估算误差可达0.5m。许多 学者先后推荐了改正模型，比较精确的一种形 式如下，其中HT为观测站的高程，单位m。 sww sdd S S sinh/ sinh/ )11000(1046512. 7 )16.273(72.1484013610552. 1 )25. 6sin(/ )25. 6sin(/ 2 0 2 5 212 212 T k w Tk k d sww sdd H T e S HT T P S hS hS 目前采用的各种对流层模型，即使应用实时测量的 气象资

11、料，电磁波的传播路径，经过对流层折射 改正后的残差，仍保持在对流层影响的5%左右。 减弱对流层折射改正项残差影响主要措施： 尽可能充分地掌握观测站周围地区的实时气象资料。 利用水汽辐射计，准确地测定电磁波传播路径上的 水汽积累量，以便精确的计算大气湿分量的改正 项。但设备庞大价格昂贵，一般难以普遍采用。 当基线较短时（20km），在稳定的大气条件下，利 用相对定位的差分法来减弱大气折射的影响。 完善对流层大气折射的改正模型。 3.电离层的影响与改正 在电离层中，由于太阳和其它天体的强烈辐射，大部分 气体分子被电离，产生了密度很高的自由电子，在离 子化的大气中，单一频率正弦波相折射率的弥散公式：

12、 式中et为电荷量，me为电荷质量，Ne为电子密度，0为真 空介质常数。当取常数值et=1.602110-19，me=9.11 10-31，0=8.859 10-12，并略去二次微小项，可得： 21 0 22 2 4 1 e te p mf eN n 2 28.401 f N n e p 根据群折射率与相折射率的关系，可得 可见，在电离层中，相折射率和群折射率是不同的， GPS定位中，对于码相位测量和载波相位测量的修正 量，应采用群折射率和相折射率分别计算。当电磁波 沿天顶方向通过电离层时，由于折射率的变化而引起 的传播路径距离差和相位延迟，一般可写为： 2 28.401 f N n f n

13、fnn e g p pg s s dsn c f dsn ) 1( ) 1( 由相折射率和群折射率引起的路径传播误差(m)和时间 延迟(ns)分别为 2 7 2 2 7 2 103436. 1 28.40 103436. 1 28.40 f N t f N f N t f N g g p p 在电离层中产生的各种延迟量，对确定的电磁波频率， 只有电子密度是唯一的独立变量。 实际资料分析表明，电离层的电子密度，白天约为夜间 的5倍，一年中冬季与夏季相差4倍，太阳活动高峰期 约为低峰期的4倍。电离层电子密度的大致变化范围在 109-31012电子数/m3。沿天顶方向电子密度总量，日 间为51017

14、电子数/m2，夜间为51016电子数/m2。此 外，电子密度在不同高度、不同时间都有明显差别。 当电磁波的传播方向偏离天顶方向时，电子总量会明显 增加，在倾角为hs方向上，电子总数Nh有如下近似： 电离层对不频率电磁波沿天顶方向传播路径的影响如下 sNNhsinh/ 单频400MHz1600MH z 2000MH z 8000MH z 平均50m3m2m0.12m 90%小 于 250m15m10m0.6m 最大500m30m20m1.2m 由于影响电离层电子密度的因素复杂（时间、高 度、太阳辐射及黑子活动、季节和地区等）， 难以可靠地确定观测时刻沿电磁波传播路线的 电子总量。对GPS单频接收

15、用户，一般均利用 电离层模型来近似计算改正量，但目前有效性 不会优于75%。即当电离层的延迟为50m，经 过模型改正后，仍含有约12.5m的残差。 为减弱电离层的影响，比较有效的措施为： （1）利用两种不同的频率进行观测 两种频率电磁波同步观测时电离层对传播路径的影响分 别为 可得消除电离层折射影响的距离： 2 2 1 2 2 2 1 12 2 1 28.40 28.40 f f f N f N ff f f 2 2 2 1 00 0 122 11 f f fff ff 2 1 2 2 2 2 0 1 ff f f 2 2 2 1 2 2 121 f ff fff 2 1 2 2 2 2 1

16、ff f f 同理可得不同频率电磁波的相位延迟关系以及经过电离 层折射改正后的相位值： 目前，为进行高精度卫星定位，普遍采用双频观测技术， 以便有效减弱电离层折射影响；不同的双频组合，对 电离层影响的改善程度也不同，改正后的残差如下： 2 1 28.40 28.40 2 1 cf N cf N f f 2 1 2 2 2 2 2 1 0 2111 ff f f f ffff 双频150/400MHz400/2000MH z 1227/1572MH z 2000/8000MH z 平均0.6m0.9cm0.3cm0.04cm 90%小于10m6.6cm1.7cm0.21cm 最大36m22cm4

17、.5cm0.43cm (2)两观测站同步观测量求差 用两台接收机在基线的两端进行同步观测，取其观测量 之差。因为当两观测站相距不太远时，卫星至两观测 站电磁波传播路径上的大气状况相似，大气状况的系 统影响可通过同步观测量的差分而减弱。 该方法对小于20km的短基线效果尤为明显，经过电离层 折射改正后，基线长度的相对残差约为10-6。故在短基 线相对定位中，即使使用单频接收机也能达到相当高 精度。但随着基线长度的增加，精度将明显降低。 4.3GPS卫星的测距码信号卫星的测距码信号 1.关于GPS卫星信号 GPS卫星所发射的信号包括载波信号、P码（或Y码）、 C/A码和数据码（或D码）等多种信号分

18、量，其中P码 和C/A码统称为。 GPS卫星信号的产生与构成主要考虑了如下因素； （1）适应多用户系统要求。 （2）满足实时定位要求。 （3）满足高精度定位需要。 （4）满足军事保密要求。 2.码与码的产生 （1）码的概念 在现代数字通信中，广泛使用二进制数（0和1）及其 组合，来表示各种信息。表达不同信息的二进制数及 其组合，称为。一位二进制数叫一个码元码元或一比特比特。 比特为码和信息量的度量单位。 如果将各种信息例如声音、图象和文字等通过量化， 并按某种预定规则，表示成二进制数的组合形式，则 这一过程称为编码编码。 在二进制数字化信息的传输中，每秒传输的比特数称 为数码率，表示数字化信息

19、的传输速度，单位为bit/s。 （2）随机噪声码 既然码是用以表达各种信息的二进制数的组合，是一组 二进制的数码序列，则这一序列就可以表达成以0和1 为幅度的时间函数。假设一组码序列u(t)，对某一时刻 来说，码元是0或1完全是随机的，但出现的概率均为 1/2。这种码元幅度的取值完全无规律的码序列，称为 序列（或序列）。它是一种非周期 性序列，无法复制，但其自相关性好。而相关性的好 坏，对提高利用GPS卫星码信号测距精度，极其重要。 为了说明随机码的自相关性，现将随机序列u(t)平移k 个码元，得到一个新的随机序列u(t)，如果两随机序列 u(t)和u(t)所对应的码元中，相同的码元数（同为0

20、或1） 为Au，相异的码元数为Bu，则随机序列u(t)的自相关 系数R(t)定义为： R(t) = (Au - Bu) / (Au + Bu) 当平移的码元数k=0，说明两个结构相同的随机码序列， 相应的码元相互对齐， Bu=0，自相关系数R(t)=1。当 k0时，由于码序列的随机性，当序列中码元数充分大 时，则AuBu，即自相关系数R(t) 0。于是，根据码 序列自相关系数的取值，可以判断两个随机码序列的 相应码元是否对齐。 假设GPS卫星发射的是一个随机码序列u(t)， 而GPS接收机若能同时复制出结构与之相同的 随机码序列u(t)，则由于卫星信号时间传播延 迟的影响，被接收的u(t)与u

21、(t)之间产生了平移， 即相应的码元错开，因而R(t) 0。如果通过一 个时间延迟器来调整u(t)，使之与u(t)的码元相 互完全对齐，即有R(t) =1。则可以从接收机的 时间延迟器中测出卫星信号到达用户接收机的 准确传播时间，从而准确测定站星距离。 （3）伪随机噪声码及其产生 尽管随机码具有良好的自相关性，但却是一种非周期 序列，不服从任何编码规则，实际中无法复制和利用。 GPS采用了一种伪随机噪声码（Pseudo Random NoicePRN）简称。它的特点是： 具有随机码的良好自相关性，又具有某种确定的编码 规则，是周期性的，容易复制。 伪随机码是由一个“多极反馈移位寄存器”的装置产

22、 生的。移位寄存器由一组连接在一起的存储单元组成， 每个存储单元只有0或1两种状态。移位寄存器的控制 脉冲有两个：钟脉冲和置1脉冲。移位寄存器是在钟脉 冲的驱动和置1脉冲的作用下而工作的。 假设移位寄存器是由4个存储单元组成的四级反馈移位 寄存器，当钟脉冲加到该移位寄存器后，每个存储单 元的内容，都顺序地由上一单元转移到下一单元，与 此同时，将其中某几个单元，如单元3和单元4的内容 进行模2相加，反馈给第一个单元。 1234 钟脉冲 输出 置1脉冲 模2相加 移位寄存器在工作开始时，置1脉冲的作用，使各级存储单元的内容 全部处于1状态，此后在钟脉冲驱动下，将经历15种状态。 状态编号各级状态

23、4 3 2 1 模2加反馈 3+4 末极输出的 二进制数 1 2 3 4 1 1 1 1 1 1 1 0 1 1 0 0 1 0 0 0 0 0 0 1 1 1 1 1 5 6 7 8 0 0 0 1 0 0 1 0 0 1 0 0 1 0 0 1 0 0 1 1 0 0 0 1 9 10 11 12 0 0 1 1 0 1 1 0 1 1 0 1 1 0 1 0 0 1 0 1 0 0 1 1 13 14 15 0 1 0 1 1 0 1 1 0 1 1 1 1 1 1 0 1 0 移位寄存器在经历了上表所列的15种状态后，再重复 全1状态，完成一个最大周期。同时，从第四级存储单 元也输出一

24、个最大周期为15tu的二进制数序列， tu为两 个钟脉冲的时间间隔。这种周期最大的二进制数序列， 通常称为。 上例中是一个码长包含有15个码元的周期序列，其中 任意4个连续的二进制数所构成的码都不相同，任何一 个码在周期序列中都有相应确定的位置和时刻。 一般情况下，对于一个r级反馈移位寄存器来说，将产 生更复杂的周期性m序列，此时移位寄存器可能经历 的状态（码长）为Nu=2r-1，最多包含Nu个码元，最大 周期为Tu=Nutu。 在r级反馈移位寄存器所产生的周期性m序列中，有时 可以截取其中的一部分，组成一个新的周期性序列加 以利用，这种新的周期较短的序列，称为截短序列或截短序列或 截短码截短

25、码。实际中有时还需要将多个周期较短的m序列， 按预定的规则，构成一个周期较长的序列，称为 由于移位寄存器不容许出现全“0”状态，因此2r-1个码元中， “1”的个数总比“0”的个数多一个。这样当两个周期相同 的m序列其相应码完全相齐时，自相关系数R(t)=1,而在其他 情况有： R(t)-1/Nu=-1/(2r-1) 当r足够大时，就有R(t)=0,所以伪随机噪声码与随机噪声码一 样，具有良好的自相关性，又是一种结构稳定，可以复制 的周期性序列。用户接收机在接受到卫星发来的伪随机码 信号后，当即复制相同结构的伪随机码，然后通过和接受 到的码信号比较（相关），精确测定信号的传播时间延迟， 进一步

26、计算某一时刻测站和卫星的距离。 3.GPS的测距码 GPS卫星所采用的两种测距码，即C/A码和P码（或Y 码），均属于伪随机码。 （1）C/A码：是由两个10级反馈移位寄存器组合而产生。 码长Nu=210-1=1023比特，码元宽为tu=1/f1=0.97752s, （ f1为基准频率f0的10分之1，1.023 MHz）,相应的距 离为293.1m。周期为Tu= Nutu=1ms，数码率为 1.023Mbit/s。 C/A码的码长短，共1023个码元，若以每秒50码元的速 度搜索，只需20.5s，易于捕获，称。 码元宽度大，假设两序列的码元对齐误差为为码元宽度 的100分之1，则相应的测距误

27、差为2.9m。由于精度低， 又称。 （2）P码 P码产生的原理与C/A码相似，但更复杂。发生电路采 用的是两组各由12级反馈移位寄存器构成。码长 Nu2.351014比特，码元宽为tu=1/f0=0.097752s，相 应的距离为29.3m。周期为Tu= Nutu 267d，数码率为 10.23Mbit/s。 P码的周期长，267天重复一次，实际应用时P码的周 期被分成38部分，（每一部分为7天，码长约6.19 1012比特），其中1部分闲置，5部分给地面监控站使 用，32部分分配给不同卫星，每颗卫星使用P码的不 同部分，都具有相同的码长和周期，但结构不同。P码 的捕获一般是先捕获C/A码，再

28、根据导航电文信息，捕 获P码。由于P码的码元宽度为C/A码的1/10，若取码 元对齐精度仍为码元宽度的1/100，则相应的距离误差 为0.29m，故P码称为。 1.导航电文及其格式 导航电文是包含有关卫星的星历、卫星工作状态、时 间系统、卫星钟运行状态、轨道摄动改正、大气折射 改正和由C/A码捕获P码等导航信息的数据码（或D 码）。 导航电文也是二进制码，依规定格式组成，按帧向外 播送。每帧电文含有1500比特，播送速度50bit/s，每 帧播送时间30s。 每帧导航电文含5个子帧，每个子帧分别含有10个字， 每个字30比特，故每个子帧共300比特，播发时间6s。 为记载多达25颗卫星，子帧4

29、、5各含有25页。子帧1、 2、3和子帧4、5的每一页构成一个主帧。主帧中1、2、 3的内容每小时更新一次，4、5的内容仅当给卫星注入 新的导航电文后才得以更新。 导航电文的格式： 123 45 12345678910 30s 6s 0.02s 0.6s 25页 10个字 30比特 一帧导航电文的内容 TLM HOW数据块1时钟修正参数 TLM HOW数据块2星历表 TLM HOW数据块2星历表继续 TLM HOW数据块3卫星历书等 TLM HOW数据块3卫星历书等 1 2 3 4 5 一个子帧6s长，10个字，每字30比特 1帧 30s 1500 比特 （1）遥测字（TLMTelemetry

30、 WORD） 位于个子帧的开头，作为捕获导航电文的前导。 （2）交接字（HOWHand Over Word） 紧接各子帧的遥测字，主要向用户提供用于捕获P码的 Z记数。所谓Z记数是从每个星期六/星期日子夜零时起 算的时间记数，表示下一子帧开始瞬间的GPS时。 （3）数据块1：含有卫星钟改正参数及数据龄期、星期 的周数编号、电离层改正参数、和卫星工作状态等信 息。 卫星钟改正参数a0、a1、a2分别表示该卫星的钟差 、钟 速和钟速变化率。任意时刻t的钟改正数为t=a0+a1(t- t0e)+a2(t-t0e)2。 参考历元t0e为数据块1的基准时间，从GPS时星期六/星 期日子夜零时起算，变化于

31、0-604800s之间。 数据龄期AODA表示基准时间t0e与最近一次更新星历的 时间之差，主要用于评价钟改正数的可信程度。 现时星期编号WN：表示从1980年1月6日协调时零点起 算的GPS时星期数。 （4）数据块2：包含在2、3两个子帧里，主要向用户提 供有关计算卫星运行位置的信息。该数据一般称为卫 星星历。 （5）数据块3：包含在4、5两个子帧中，主要向用户提 供GPS卫星的概略星历及卫星的工作状态信息，称为 卫星的历书。 4.5GPS卫星信号的构成 1卫星的载波信号与调制 GPS卫星信号包含三种信号分量：载波、测距码和数据 码。信号分量的产生都是在同一个基本频率 f0=10.23MHz

32、的控制下产生，GPS卫星信号示意图如 下 基本频率 10.23MHz L1载波 1575.42MHz L2载波 1227.60MHz C/A码 1.023MHz P码 10.23MHz P码 10.23MHz 数据码 50BPS 数据码 50BPS 154 120 10 204600 从上图中可见，GPS卫星取L波段的两种不同电磁波频率 为载波，L1载波频率为1575.42MHz，波长为19.03cm； L2载波频率为1227.60MHz，波长为24.42cm。在L1载 波上，调制有C/A码、P码（或Y码）和数据码； L2载 波上，只调制有P码（或Y码）和数据码。 在无线电通信中，为有效地传播信息，一般将频率较低 的信号加载到频率较高的载波上，此时频率较低的信 号称为。 GPS卫星的测距码和数据码是采用调相技术调制到载波 上，且调制码的幅值只取0或1。如果码值取0，则对应 的码状态取+1；而码值取1时，对应码状态为-1，载波 和相应的码状态相乘后，即实现了载波的调制。 2.卫星信号的解调 为进行载波相位测量，当用户接收到卫星发播的信号 后，可通过以下两种解调技术来恢复载波相位。 （1）复制码与卫星信号相乘：由于调制码的码值是用1 的码状态来表示的，当把接收的卫星码信号与用户接 收机产生的复制码（结构与卫星测距码信号完全相同 的测距码），在两码同步的条件下相乘，即可去