第三章-卫星运动及GPS卫星信号.ppt

第三章第三章 卫星运动及信号卫星运动及信号研究研究GPS卫星卫星在协议地球坐标系中的瞬时位置在协议地球坐标系中的瞬时位置掌握掌握GPS卫星信号卫星信号1.卫星轨道在卫星轨道在GPS定位中的意义定位中的意义v 卫星在空间运行的轨迹称为卫星在空间运行的轨迹称为轨道轨道，描述卫星轨道位置和，描述卫星轨道位置和状态的参数称为状态的参数称为轨道参数轨道参数v 由于利用由于利用GPS进行导航和测量时，卫星作为位置已进行导航和测量时，卫星作为位置已知的高空观测目标，在进行绝对定位时，卫星轨道误差将知的高空观测目标，在进行绝对定位时，卫星轨道误差将直接影响用户接收机位置的精度；而在相对定位时，尽管直接影响用户接收机位置的精度；而在相对定位时，尽管卫星轨道误差的影响将会减弱，但当基线较长或精度要求卫星轨道误差的影响将会减弱，但当基线较长或精度要求较高时，轨道误差影响不可忽略较高时，轨道误差影响不可忽略v 此外，为了制订此外，为了制订GPS测量的观测计划和便于捕获卫测量的观测计划和便于捕获卫星发射的信号，也需要知道卫星的轨道参数星发射的信号，也需要知道卫星的轨道参数2.影响卫星轨道的因素及其研究方法影响卫星轨道的因素及其研究方法 卫星在空间绕地球运行时，除了受地球重力场的卫星在空间绕地球运行时，除了受地球重力场的引力作用外，还受到太阳、月亮和其它天体的引引力作用外，还受到太阳、月亮和其它天体的引力影响，以及太阳光压、大气阻力和地球潮汐力力影响，以及太阳光压、大气阻力和地球潮汐力等因素影响。

卫星实际运行轨道十分复杂，难以等因素影响卫星实际运行轨道十分复杂，难以用简单而精确的数学模型加以描述用简单而精确的数学模型加以描述v 在各种作用力对卫星运行轨道的影响中，在各种作用力对卫星运行轨道的影响中，地地球引力场球引力场的影响为主，其它作用力的影响相对要的影响为主，其它作用力的影响相对要小的多若假设地球引力场的影响为小的多若假设地球引力场的影响为1 1，其它引力，其它引力场的影响均小于场的影响均小于1010-5-5人造卫星运动的轨道参数取决于人造卫星运动的轨道参数取决于所受的力的作用所受的力的作用•1.地球质心引力（中心引力）地球质心引力（中心引力）•2.摄动力摄动力•与之对应的无摄与受摄运动与之对应的无摄与受摄运动3.1 概述概述卫星的受力分为两类•中心力：中心力：是假设地球为均质球体的引力（质量集是假设地球为均质球体的引力（质量集中于球体的中心）决定着卫星运动的基本规律中于球体的中心）决定着卫星运动的基本规律和特征，由此决定的卫星轨道，可视为和特征，由此决定的卫星轨道，可视为理想轨道理想轨道，，是分析卫星实际轨道的基础是分析卫星实际轨道的基础 摄动力或非中心力：摄动力或非中心力：包括地球非球形对称的作包括地球非球形对称的作用力、日月引力、大气阻力、光辐射压力以及地用力、日月引力、大气阻力、光辐射压力以及地球潮汐力等。

摄动力使卫星的运动产生一些小的球潮汐力等摄动力使卫星的运动产生一些小的附加变化而偏离理想轨道，同时偏离量的大小也附加变化而偏离理想轨道，同时偏离量的大小也随时间而改变随时间而改变受摄运动，受摄轨道）（受摄运动，受摄轨道）研究流程研究流程二体问题二体问题： 忽略所有的摄动力，仅考虑地球质心引力忽略所有的摄动力，仅考虑地球质心引力研究卫星相对于地球的运动，在天体力学研究卫星相对于地球的运动，在天体力学中，称之为二体问题中，称之为二体问题 ①①它是卫星运动的第一近似描述；它是卫星运动的第一近似描述； ②②它是至今惟一能得到的严密分析解的运动；它是至今惟一能得到的严密分析解的运动； ③③它是全部作用力下的卫星运动更精确解的它是全部作用力下的卫星运动更精确解的基础二体问题二体问题1.根据万有引力定律及牛顿第二定律导出卫根据万有引力定律及牛顿第二定律导出卫根据万有引力定律及牛顿第二定律导出卫根据万有引力定律及牛顿第二定律导出卫星与地球的运动方程星与地球的运动方程星与地球的运动方程星与地球的运动方程§卫星卫星卫星卫星S S 相对于地球相对于地球相对于地球相对于地球质心质心质心质心O O 的加速度的加速度的加速度的加速度：：：：§卫星卫星卫星卫星S S 在万有引力在万有引力在万有引力在万有引力作用下的加速度：作用下的加速度：作用下的加速度：作用下的加速度：§地球地球地球地球O O 在万有引力作在万有引力作在万有引力作在万有引力作用下的加速度：用下的加速度：用下的加速度：用下的加速度：§卫星与地球间卫星与地球间卫星与地球间卫星与地球间引力的大小：引力的大小：引力的大小：引力的大小： G为引力常数，为引力常数，M为地球质量，为地球质量，ms为卫星质为卫星质量，量，r为卫星的地心向径。

根据上式来研究为卫星的地心向径根据上式来研究地球和卫星之间的相对运动问题，在天体地球和卫星之间的相对运动问题，在天体力学中称为力学中称为二体问题二体问题引力加速度决定了引力加速度决定了卫星绕地球运动的基本规律卫星在上述卫星绕地球运动的基本规律卫星在上述地球引力场中的无摄运动，也称地球引力场中的无摄运动，也称开普勒运开普勒运动动，其规律可通过开普勒定律来描述其规律可通过开普勒定律来描述椭圆轨道的发现•约翰尼斯约翰尼斯.开普勒开普勒 (1571 - 1630)•利用第谷多年积累的观测资料，仔细分利用第谷多年积累的观测资料，仔细分析研究，发现了行星沿椭圆轨道运行析研究，发现了行星沿椭圆轨道运行•提出行星运动三定律（即开普勒定律）提出行星运动三定律（即开普勒定律）u第谷第谷 (1546 - 1601): 著名的天文观测者著名的天文观测者u对天体进行了精确，细致的观测对天体进行了精确，细致的观测3.2.2 开普勒定律3-73-103.2.3 无摄卫星运动的轨道参数无摄卫星运动的轨道参数Y 升交点赤经Ω轨道倾角i近地点角距 真近点角Vf轨道长半轴a 近地点近地点赤道平面赤道平面轨道偏心率升交点升交点 Ω卫星卫星rX ωZiV Vf fb§通过开普勒轨道通过开普勒轨道6个参数可以确定出卫星在个参数可以确定出卫星在轨道平面上的瞬间位置和速度。

轨道平面上的瞬间位置和速度 § a a、、、、 e e、、、、  、、、、 i i、、、、    、、、、 V V 由卫星发射条由卫星发射条件决定，已知件决定，已知为为时间的函数时间的函数，，需计算出需计算出3.2.4 真近点角计算真近点角计算•由平近点角由平近点角MS计算偏近点角计算偏近点角Es ，再由偏近，再由偏近点角计算真近点角，而平近点角是时间的点角计算真近点角，而平近点角是时间的函数Es见下图见下图真近点角的计算真近点角的计算§以椭圆轨道以椭圆轨道O’ O’ 为圆心，以椭圆半径为圆心，以椭圆半径a a为半径作为半径作一辅助圆一辅助圆O’-AS’A’O’-AS’A’§偏近点角偏近点角( (E E ) )§平近点角平近点角( (M M ) ) M=n M=n ( (t- t-  ) )§真近点角真近点角( (V V ) )卫星平均卫星平均角速度角速度卫星过近地点卫星过近地点时刻，已知时刻，已知真近点角的计算真近点角的计算偏近点角偏近点角( (E E ) )平近点角平近点角( (M M ) )真近点角真近点角( (V V ) )开普勒方程开普勒方程: : M=M= E - esinE E - esinE平近点角与偏近点角•由开普勒方程：•ES=MS+eSsinES•偏近点角采用迭代法求解。

真近点角fs•ascosEs=rcosfs+ases•最后得到3-203.3 卫星的瞬时位置与瞬时速度•瞬时位置：瞬时位置：•1.卫星在轨道直角坐标系的位置卫星在轨道直角坐标系的位置3-22•2.卫星在天球坐标系的位置（对轨道直角坐卫星在天球坐标系的位置（对轨道直角坐标系进行旋转）标系进行旋转）3-23•3.卫星在地球坐标系中的位置卫星在地球坐标系中的位置3-251.卫星在轨道直角坐标系的位置卫星在轨道直角坐标系的位置 取直角坐标系的原点与地球质心相重合，取直角坐标系的原点与地球质心相重合， s轴指向近地点、轴指向近地点、 s轴垂直于轨道平面向上轴垂直于轨道平面向上 ,  s轴在轨道平面上垂直于轴在轨道平面上垂直于 s轴构成右手系，轴构成右手系，则卫星在任意时刻的坐标为则卫星在任意时刻的坐标为:2.卫星在天球坐标系的位置卫星在天球坐标系的位置在轨道平面直角坐标系中只确定了卫星在轨在轨道平面直角坐标系中只确定了卫星在轨道平面上的位置，而轨道平面与地球体的道平面上的位置，而轨道平面与地球体的相对定向尚需由轨道参数相对定向尚需由轨道参数 、、i和和 s确定天球坐标系（天球坐标系（x,y,z）与轨道坐标系）与轨道坐标系( s,  s,  s)具有相同的原点，差别在于坐标系的定具有相同的原点，差别在于坐标系的定向不同，为此需将轨道坐标系作如下旋转：向不同，为此需将轨道坐标系作如下旋转：用旋转矩阵表示如下：3.卫星在地球坐标系中的位置卫星在地球坐标系中的位置利用利用GPS定位时，应使观测卫星和观测站的位置处于统一的定位时，应使观测卫星和观测站的位置处于统一的坐标系统。

由于瞬时地球空间直角坐标系与瞬时天球空间坐标系统由于瞬时地球空间直角坐标系与瞬时天球空间直角坐标系的差别在于直角坐标系的差别在于x轴的指向不同，若取其间的夹角轴的指向不同，若取其间的夹角为春分点的格林尼治恒星时为春分点的格林尼治恒星时GAST，则在地球坐标系中卫，则在地球坐标系中卫星的瞬时坐标（星的瞬时坐标（X，，Y，，Z）与天球坐标系中的瞬时坐标）与天球坐标系中的瞬时坐标（（x，，y，，z）存在如下关系：）存在如下关系：•瞬时速度：瞬时速度：•1.轨道直角坐标系中三维坐标量对时间求导轨道直角坐标系中三维坐标量对时间求导得得3.31式式•其余两坐标系同理求解其余两坐标系同理求解3.4 受摄运动（了解）受摄运动（了解）v对于卫星精密定位来说，必须考虑对于卫星精密定位来说，必须考虑地球引地球引力场摄动力、日月摄动力、大气阻力、光力场摄动力、日月摄动力、大气阻力、光压摄动力、潮汐摄动力对卫星运动状态的压摄动力、潮汐摄动力对卫星运动状态的影响影响考虑了摄动力作用的卫星运动称为考虑了摄动力作用的卫星运动称为卫星的卫星的受摄运动受摄运动卫星的受摄运动轨道卫星受到的摄动力卫星受到的摄动力§地球引力场摄动力影响约为地球引力场摄动力影响约为地球引力场摄动力影响约为地球引力场摄动力影响约为1010-3-3量级，其他摄量级，其他摄量级，其他摄量级，其他摄动力影响大多小于或接近于动力影响大多小于或接近于动力影响大多小于或接近于动力影响大多小于或接近于1010-6-6量级。

量级•由于在由于在GPS星历电文中必须包含卫星位置和速度星历电文中必须包含卫星位置和速度的非常精确的信息的非常精确的信息,除非在非常接近于历元的时间除非在非常接近于历元的时间内，否则仅用开普勒密切轨道元素以计算内，否则仅用开普勒密切轨道元素以计算GPS卫卫星的位置是不够的解决这个问题的一个方法是星的位置是不够的解决这个问题的一个方法是十分十分频繁地更新频繁地更新GPS星历电文，但是不现实星历电文，但是不现实的•因此用因此用“校正校正”参数对参数对GPS星历电文中的开普勒星历电文中的开普勒密切轨道参数进行了扩充，以使用户能在卫星星密切轨道参数进行了扩充，以使用户能在卫星星历电文的两次更新之间的时间内能十分精确地估历电文的两次更新之间的时间内能十分精确地估计开普勒轨道元素对于一个特定的星历电文，计开普勒轨道元素对于一个特定的星历电文，在历元之后的任何时候，在历元之后的任何时候，GPS接受机用校正参数接受机用校正参数以估计在所希望时刻的轨道元素以估计在所希望时刻的轨道元素3.4 受摄运动（了解）受摄运动（了解）•无摄运动轨道为开普勒轨道无摄运动轨道为开普勒轨道•但是，地球受到很多摄动力的影响。

但是，地球受到很多摄动力的影响•图图3-8•表表3-1•注：受摄运动方程是位置、速度、时间的注：受摄运动方程是位置、速度、时间的函数，及其复杂，在此不作讨论函数，及其复杂，在此不作讨论•其单独影响见～，自己阅读其单独影响见～，自己阅读3.5 GPS星历星历•卫星星历是描述卫星运动轨道的信息描述卫星运动轨道的信息•卫星星历包含某一时刻的卫星轨道参数（也叫轨道根数）及其变化率•根据卫星星历，可以计算出任一时刻的卫星位置及其速度•精确的轨道信息是精密定位的基础3.5 GPS星历星历•GPS卫星星历可以分为：分为：•预报星历预报星历（广播星历）（广播星历）•后处理星历后处理星历（精密星历）（精密星历）3.5.1 预报星历（广播星历）预报星历（广播星历）§预报星历：预报星历：是是通过通过卫星发射的含有轨道信息卫星发射的含有轨道信息的的导航电文导航电文传输给用户的，用户接收机接收传输给用户的，用户接收机接收到这些信号，经过解码便可获得所需要的卫到这些信号，经过解码便可获得所需要的卫星星历，这种星历也叫星星历，这种星历也叫广播星历广播星历. .§( (用户在观测时可以用户在观测时可以实时实时地得到地得到) )。

3.5.1 预报星历（广播星历）预报星历（广播星历）•预报星历通常包括相对某一参考历元某一参考历元的开普勒轨道参数和必要的轨道摄动改正项参数：–开普勒轨道的6个参数；–反映摄动力影响的9个参数；–参考时刻；–星历数据龄期；•预报星历包含的参数是根据GPS监测站约一周的观测资料推算的，也称参考星历轨道参数：星历的基准时间：星历的基准时间：半长轴的平方根：半长轴的平方根：偏心率：偏心率：参考时刻的倾角：参考时刻的倾角：升交点赤经：升交点赤经：近地点角距：近地点角距：平近点角：平近点角：轨道倾角变化率：轨道倾角变化率：星历表数据龄期：星历表数据龄期：升交点赤经变化率：升交点赤经变化率：平均角速度的校正：平均角速度的校正：纬度幅角余弦的校正：纬度幅角余弦的校正：纬度幅角正弦的校正：纬度幅角正弦的校正：轨道半径余弦的校正：轨道半径余弦的校正：轨道半径正弦的校正：轨道半径正弦的校正：倾角余弦的校正：倾角余弦的校正：倾角正弦的校正：倾角正弦的校正3.5.1 预报星历预报星历•预报星历（参考星历）只代表卫星在参考历元的瞬时轨道参数，在摄动力的影响下，卫星的实际轨道将偏离参考轨道•偏离程度主要取决于：观测历元与参考历元的时间差。

3.5.1 预报星历（广播星历）预报星历（广播星历）•用轨道参数的摄动项对已知的卫星星历加以改正，可以外推出任意观测历元的卫星星历•只要保证外推时间间隔不太长，可以保证卫星预报星历的精度•地面的监控站时刻观测卫星的运行轨道；主控站每天更新卫星的参考星历，由注入站每天向卫星注入新的参考星历3.5.1 预报星历（广播星历）预报星历（广播星历）•GPS卫星发射的广播星历，每小时更新一次•那么外推星历的时间间隔将不会超过0.5小时•预报星历的精度一般为20～40m预报星历的编制和传送过程：监测站监测卫星的运行状态 主控站编制卫星星历注入站将卫星星历注入卫星卫星向用户发送卫星星历用户接收机接收卫星星历（监测站接收机）周期1小时循环闭合的过程 用广播星历计算卫星位置用广播星历计算卫星位置在在用用GPS信信号号进进行行导导航航定定位位以以及及制制定定观观测测计计划划时时，，都都必必须须已已知知GPS卫卫星星在在空空间间的的瞬瞬时时位位置置卫卫星星位位置置的的计计算算是是根根据据卫卫星星电电文文所所提提供供的的轨轨道道参参数数按按一一定的公式计算的定的公式计算的1 1、计算卫星运行的平均速度、计算卫星运行的平均速度n根据开普勒第三定律，卫星运行的平均角速度根据开普勒第三定律，卫星运行的平均角速度n0可以用下式计算可以用下式计算 式式 中中 ，， u为为 WGS-84WGS-84坐坐 标标 系系 中中 的的 地地 球球 引引 力力 常常 数数 , ,且且 u=3.986 3.986 005×10005×1014m m3/s/s2。

平平均均角角速速度度n0 0加加上上导导航航电电文文中中给给出出的的摄摄动动改改正正数数∆n，，便得到卫星运行的平均角速度便得到卫星运行的平均角速度n然后将观测时刻然后将观测时刻t归化到归化到GPSGPS时系时系式中，式中，tk称作相对于参考时刻称作相对于参考时刻t0e的归化时间的归化时间3 3、观测时刻卫星平近点角、观测时刻卫星平近点角Mk的计算的计算4 4、采用迭代算法计算偏近点角、采用迭代算法计算偏近点角Ek式中式中M0是导航电文给出的参考时刻是导航电文给出的参考时刻t0e的平近点角的平近点角2 2、计算归化时间、计算归化时间tk首先对观测时刻首先对观测时刻t作卫星钟差改正作卫星钟差改正5 5、计算真近点角、计算真近点角Vk 6 6、计算升交距角、计算升交距角Φk 其中，其中，w为近地点角距，由卫星星历给出为近地点角距，由卫星星历给出 7 7、计算摄动改正数、计算摄动改正数δu、、δr、、δi 式中，式中，δδu、、δδr、、δδi分别为升交距角分别为升交距角u的摄动量、卫星矢径的摄动量、卫星矢径r的摄的摄动量和轨道倾角动量和轨道倾角i i的摄动量的摄动量Cuc , Cus，，Crc , Crs，，Cic , Cis由导航电由导航电文给出。

文给出 8 8、计算经摄动改正后的升交距角、计算经摄动改正后的升交距角uk、卫星矢径、卫星矢径rk和轨道倾角和轨道倾角ik 9 9、计算卫星在轨道平面坐标系中的坐标、计算卫星在轨道平面坐标系中的坐标 1010、计算观测时刻升交点经度、计算观测时刻升交点经度Ωk we=7.29211567×107.29211567×10-5-5rad/s,rad/s,为地球自转的速率，其它有关参数从电为地球自转的速率，其它有关参数从电文中得到文中得到1111、计算卫星在地固坐标系中的空间直角坐标、计算卫星在地固坐标系中的空间直角坐标 1212、计算卫星在协议地球坐标系中的空间直角坐标、计算卫星在协议地球坐标系中的空间直角坐标 其中，其中，XP、、YP为地极瞬时坐标为地极瞬时坐标 在GPS定位中，轨道平面坐标系的x轴指向升交点，y轴垂直于x轴指向地极北方向，原点位于地球质心 3.5.2 后处理星历后处理星历•后处理星历，是根据地面跟踪站所获得的精密观测资料计算而得到的星历，是一种不包含外推误差不包含外推误差的实测星历•由精密观测资料计算而得，不含外推误差，由精密观测资料计算而得，不含外推误差，精度可达精度可达米米级。

级用精密星历计算卫星位置用精密星历计算卫星位置§精密星历精密星历是按一定的时间间隔（通常为是按一定的时间间隔（通常为1515分钟）来分钟）来给出卫星在空间的三维坐标、三维运动速度及卫星给出卫星在空间的三维坐标、三维运动速度及卫星钟改正数等信息钟改正数等信息§观测瞬间的卫星位置及运动速度可采用内插法求得观测瞬间的卫星位置及运动速度可采用内插法求得•后处理星历不是通过卫星广播的•需要通过无线电、网络等通信方式向用户传递•是有偿服务优点： 轨道参数非常准确，也称精密星历缺点： 不能做到实时后处理星历的编制和传送过程：建立卫星跟踪系统，随时监测卫星运动状态计算卫星星历用磁带或通过电视、电传、卫星通讯等方式，向用户提供以往观测时刻的星历 建立和维持一个独立的跟踪系统来精密测定建立和维持一个独立的跟踪系统来精密测定GPS卫星的轨道，技术复杂，投资大，因此，利用卫星的轨道，技术复杂，投资大，因此，利用GPS预报星历进行精密定位工作仍是目前一个重预报星历进行精密定位工作仍是目前一个重要的研究和开发领域要的研究和开发领域 中国的中国的GPS卫星跟踪站在卫星跟踪站在“八五八五”期间已经建成，期间已经建成， “九五九五”期间已向国内用户提供服务。

现在正在期间已向国内用户提供服务现在正在进行数据处理的全国进行数据处理的全国GPS2000网，其点数可达网，其点数可达2000余个，其中余个，其中GPS卫星永久性跟踪站增加为卫星永久性跟踪站增加为25个，将来还要继续增加到个，将来还要继续增加到250个左右，以满足全国个左右，以满足全国动动‘态或高精度定位用户的需要态或高精度定位用户的需要3.6 GPS卫星的伪随机噪声码•GPS卫星向广大用户发送的导航定位信号，是经过伪随机噪声码扩频扩频的调相信号•导航信号采用L波段（22cm）做载波L波段的优点：–减少拥挤，避免“撞车”目前L波段的频率占用率较低–适应扩频，传送宽带信号在占用率较低的L波段，易于传送宽带信号3.6 GPS卫星的伪随机噪声码–卫星高轨运行能获得较大的多普勒频移载波频率越高、多普勒频移越大，有利于测量用户的行驶速度–大气衰减小，有益于研制用户设备大气对某些频率的电磁波产生谐振吸收使用L波段避免了谐振吸收，有利于用低成本设备测量接收GPS信号3.6 GPS卫星的伪随机噪声码•码，表达信息的二进制数及其组合码，表达信息的二进制数及其组合3.6.1 3.6.1 有关码的基本知识有关码的基本知识n（（1 1）码的基本概念）码的基本概念• 在现代数字化通讯中，广泛使用在现代数字化通讯中，广泛使用二进制数二进制数（即（即“0”和和“1”）及其组合表示各种信息，）及其组合表示各种信息，这这些表达不同信息的二进制数及其组合称为码些表达不同信息的二进制数及其组合称为码。

• 一位二进制数叫做一个一位二进制数叫做一个码元码元或或一比特一比特；；• 比特比特（（bit））被取为码的度量单位；被取为码的度量单位；• 数码率数码率：在二进制数字化信息的传输中，每：在二进制数字化信息的传输中，每秒中传输的比特数，记为秒中传输的比特数，记为bit/s或或BPS（（2 2））随机噪声码随机噪声码§定义定义：假设一组码序列：假设一组码序列u(t)u(t)，，对某一时刻来说，对某一时刻来说，码元是码元是0 0或或1 1完全是随机的，但其出现的概率均完全是随机的，但其出现的概率均为为1/21/2这种码元幅度的取值完全无规律的码序这种码元幅度的取值完全无规律的码序列，称为列，称为随机码序列随机码序列1 1 10 0 0 01 10 01 1 1 1随机码序列随机码序列n随机码的随机码的自相关性自相关性 设将随机码序列设将随机码序列u(t)u(t)平移平移k k个码元个码元，而得到一，而得到一个新的随机码序列个新的随机码序列u u1 1(t)(t)如果两随机序列所对应如果两随机序列所对应的码元中，相同的码元数（的码元中，相同的码元数（0 0或或1 1）为）为SuSu，相异的码，相异的码元数为元数为DuDu，则随机序列，则随机序列u(t)u(t)的自相关系数的自相关系数R(t)R(t)定义定义为：为：•自相关函数是时间的函数，自相关函数是时间的函数，GPS通过调整自相关通过调整自相关函数值来计算函数值来计算 信号传播时间。

信号传播时间•随机码的特点：随机码的特点：非周期序列，无法复制非周期序列，无法复制两种结构相同的码序列完全对齐时R（t）=1两种结构相同的码序列平移k个码元时，R（t）≈03.6.2 伪伪随机噪声码随机噪声码n伪随机噪声码伪随机噪声码：看似一组杂乱无章的随机噪声码，：看似一组杂乱无章的随机噪声码，其实是其实是按照一定规律编排起来按照一定规律编排起来的、的、可以复制可以复制的的周周期性期性的二进制序列，且的二进制序列，且具有类似于随机噪声码的具有类似于随机噪声码的自相关性自相关性n GPS采用了一种伪随机噪声码（采用了一种伪随机噪声码（Pseudo Random Noise-PRN)，，简称伪随机码或伪码简称伪随机码或伪码伪随机码的产生•伪随机码是使用移位寄存器产生的，这种移位寄存器叫做抽头式反馈移位寄存器（tapped feedback shift register）•产生的伪随机码也称为m序列算法实例见表算法实例见表3.4四级反馈移位寄存器D1D2D3D4x0111110111100111000111000001000001001001111000011001011101011101001101111100一个一个r级级移位寄存器所产生的移位寄存器所产生的m序列，在一个周序列，在一个周期内码元最大个数期内码元最大个数Nu=2r-1 Nu称为码长称为码长。

M序列的最大周期为：序列的最大周期为：Tu=（（2r-1））tu= Nu tutu为为码元宽度码元宽度（钟脉冲的时间间隔）（钟脉冲的时间间隔）移位寄存器不允许出现全移位寄存器不允许出现全0状态，因此状态，因此1的个数的个数总比总比0的个数多的个数多1两个周期相同的两个周期相同的m序列其相应码元完全对齐时，序列其相应码元完全对齐时，R（（t））=1其他情况：其他情况： R（（t））=-1/Nu3.6.3 测距码信号两种测距码信号：两种测距码信号：C/A码，码，p码码均为均为伪随机噪声码伪随机噪声码1、粗码（C/A码）§定义定义：用于进行：用于进行粗略测距粗略测距和和捕获精码捕获精码的测距码的测距码§产生产生：：由两个由两个10级反馈移位寄存器相组合而产生出级反馈移位寄存器相组合而产生出1023种不同结构的种不同结构的C/A码，但其码长和数码率均相同码，但其码长和数码率均相同§不同的不同的GPS卫星采用结构相异的卫星采用结构相异的C/A码§特性特性：：码长码长 L=1023比特；比特；码元宽度码元宽度tu=1/f1=0.97752us(相应距离为相应距离为293.1m)；；§特点特点: : 测距精度低，测距误差可达测距精度低，测距误差可达29.3m~2.9m,也称粗码；也称粗码；码结构公开；目前，码结构公开；目前，C/A码只调制在码只调制在L1载波载波上，故无法精确地上，故无法精确地消除电离层延迟。

消除电离层延迟C/A码•C/A码码长较短，易于捕获通过捕获C/A码的信息，可以方便地捕获P码–GPS接收机为了捕获C/A码，通常使用逐个搜索法C/A码总共有1023个码元，如果每秒搜索50个码元，仅需20.5秒便可完成•C/A码的码元宽度较大定位精度较低2 2、、P P码码§定义定义：：用于精确测定从用于精确测定从GPS卫星至接收机距卫星至接收机距离的测距码离的测距码§产生产生：其发生电路采用两组各有两个：其发生电路采用两组各有两个12级反级反馈移位寄存器所构成，情况更为复杂馈移位寄存器所构成，情况更为复杂§特点特点：测距精度高（一般为：测距精度高（一般为0.3米米)，测距误差，测距误差约为约为2.93m~0.29m；；§ 码元宽度仅为码元宽度仅为C/A码的码的1/10，也称精码；，也称精码；• 码结构不公开，保密码结构不公开，保密P码•为了捕获和跟踪P码，采用两项措施：–将P码的时间周期“截短”把每周日子夜零时作为起点，周六午夜24时作为终点，周期为7天–分步捕获法先根据C/A码解译出卫星导航电文，根据电文的HOW信息捕获P码。

3.7 3.7 导航电文导航电文§定义定义：导航电文是由：导航电文是由GPSGPS卫星卫星向向用户用户播发的一播发的一组反映卫星在空间的组反映卫星在空间的位置位置、、卫星的状态卫星的状态、、卫卫星钟的修正参数星钟的修正参数、、电离层延迟修正参数电离层延迟修正参数等重等重要数据的二进制代码，也称要数据的二进制代码，也称数据码数据码（（D D码）§ 3.7.1 3.7.1 电文格式电文格式§传输速率传输速率：：50bit/s50bit/s，以，以““帧帧””为单位向外发送；为单位向外发送；§每帧的长度为每帧的长度为1500bit1500bit，播发完一帧需，播发完一帧需30s30s；；§一帧包括一帧包括5 5个子帧每个子帧包含个子帧每个子帧包含300bit300bit，播发时，播发时间为间为6s6s§每个子帧又分为每个子帧又分为1010个字个字§每个字由每个字由30bit30bit组成§其中第四、五两个子帧各有其中第四、五两个子帧各有2525个页面，每个页面，每2525帧组帧组成一个主帧，需要成一个主帧，需要750s750s（（12.512.5分钟）才能将分钟）才能将2525个个页面全部播发完。

页面全部播发完导航电文的格式导航电文的格式3.7.2 电文内容电文内容•1、遥测码（TLM – Telemetry Word）–每一子帧的第1个字–用作捕获导航电文的前导遥测码（telemetry word, TLW）•遥测码位于每个子帧的开头，用于表明卫星注入数据的状态–第1～8bit是同步码（10001001），为各子帧编码脉冲提供一个同步起点–第9～22bit是遥测电文，包括地面监控系统注入数据时的状态信息、诊断信息等–第22、24bit是连接码–第25～30bit是奇偶校验码85奇偶校验码•一个字码（30bit）包含6bit奇偶校验码每一个bit的奇偶校验码对应前面24bit的数据码•例如：1110 0101 1101 0001 1011 1010 10111086数据码文数据码文校验码文校验码文111010101011011000111011110100•2、转换码（交接字）（HOW – Hand Over Word）–每一子帧的第2个字–主要内容：Z计数转换码（hand over word, HOW）•转换码位于每个子帧的第二个子码•作用：提供用户从捕获的C/A码转换到捕获P码的Z计数。

Z计数位于转换码的第1～17bit，从每周六/周日零时起算的时间计数通过Z计数，可以知道观测瞬间在P码周期中所处的准确位置，以便迅速捕获P码•转换码的第25～30bit为奇偶校验码88–内容：内容：•WN – GPS周，从子夜零点起算的周期数周，从子夜零点起算的周期数•L2所调制所调制测距码标识符测距码标识符 – “10”表示表示C/A码，码，“01”表示表示P(Y)码码•传输参数传输参数N – URA•Tgd – 信号在卫星内部的时延（电离层时延差）信号在卫星内部的时延（电离层时延差）•星钟数据龄期星钟数据龄期AODC，卫星时钟改正数的置信度，卫星时钟改正数的置信度•星钟改正参数（星钟改正参数（GPS卫星时钟，相对卫星时钟，相对GPS时系的差值）时系的差值）•a0（钟偏），（钟偏），a1（钟速），（钟速），a2（钟漂）（钟漂）§3 3、第一数据块、第一数据块：：§第一子帧中的第第一子帧中的第3-103-10个字4、第二数据块、第二数据块：：由第二子帧和第三子帧构成，给出了该由第二子帧和第三子帧构成，给出了该GPS卫星的卫星的广播星历参数广播星历参数–第第2、、3子帧的第子帧的第3~10个字个字–内容内容•该发送信号卫星的星历该发送信号卫星的星历 －－ 广播星历广播星历–星历参数星历参数开普勒轨道系数，轨道摄动参数，时间参数开普勒轨道系数，轨道摄动参数，时间参数5、第三数据块、第三数据块：：包含在第四、第五子帧中。

向用户提供了其余包含在第四、第五子帧中向用户提供了其余GPS卫星的概略星历及其工作状态的信息卫星的概略星历及其工作状态的信息–第4、5子帧的第3~10个字内容：所有卫星历书（概略星历）第三数据块的内容每12.5分钟重复一次第4子帧•1、6、11、12、16、19-24，备用•2-10提供25-32颗卫星历书•13-15空闲•17专用电文•18电离层改正模型数据•25提供卫星型号及健康状况第5子帧•1-24提供1-24颗卫星历书•25给出1-24颗卫星的健康状况星期编号。