第五章GPS卫星定位基本原理.docx

1、第五章GPS卫星定位基本原理5.1概述测距交会确定点： 无线电导航定位系统 卫星激光测距 定位系统无线电导航定位：三已知点三维定位，两个已知点平面 定位 .卫星大地测量中的卫星激光测距定位。 利用地面上三个 已知点上的卫星激光测距仪同时测定某一时刻至卫星的空 间距离，从而来确定卫星的空间位置。卫星定位的基本原理： 依据测距的原理：伪距法定位，载波相位测量定位，以及差分 GPS 定位。 根据待定点的状态分为：静态定位（绝对定位）和动态定位（至少有一台接收机处于运动状态）和相对定位。 利用测距码或载波相位均可进行静态定位， 实际为减少误差，可利用载波相位观测值的各种线性组合（即差分）作 为观测值，

2、获得两点之间高精度的 GPS 基线向量（即坐标 差）。5.2伪距测量 伪距测量：由卫星发射的测距码信号到达 GPS 接收机的传 播时间乘以光速所得出的量测距离。由于卫星钟、接收机钟 的误差以及无线电信号经过电离层和对流层中的延迟， 实际 测出距离 与卫星到接收机的几何距离 有一定差值，因此一般称量测出的距离为伪距。 C/A码伪距，P码伪距。伪距法定位测量定位精度不高(P码定位误差约为10m, C/A码 定位误差为20-30m)，但因其具有定位速度快， 是GPS定位系统中进行导航定位的基本方法。 作为载波相位测量中解决整波数不确定(模糊度)的辅助资料。5.2.1伪距测量伪距测量的基本原理：为什么

3、采用码相关技术来确定伪距？GPS卫星发射的测距码是按照一定规律排列的， 在一个周期内，每个码对应某一特定的时间。应该说识别出每个码的形 状特征，即用每个码的某一标志即可推算出时延值 进行伪距测量。但实际上每个码在产生过程中都带有随机误差，并 且信号经过长距离传送后也会产生变形。 所以根据码的某一标志来推算时延值 就会产生很大的误差。因此采用码相关 技术，在自相关系数R( ) max的情况下来确定信号的传播时 间。由于测距码和信号在产生的过程中不可避免地带有误 差，而且测距码在传播过程中还有变形，因而自相关系数往 往不可能达到“ 1”，只能在自相关系数为最大的情况下确定 伪距，此时基本对齐。1R

4、() Ta(t )a(t t )dt伪距测量原理GPS定位的茶础就是测距即通过测托信号从卫星传播到接收机所用的时间 获得卫星和接收机之I司的更离= cxa；）.卫厦fl才刻产生梢躅相任r心f后到达ij1 LLI !-；ra图2利用伪前机码的自相耒特件如图弟令本地码片以一左速度移动如图牡当 本地复现码与城卫星传到的码相关值达到最大时，本地码穆动了的时间即为所求 f a调整本地码延迟 ，可使相关输出达到最大值R(t) Rmax(t)t t t5-3可得机钟差tk与卫星钟差 tj，即ttk t j ，若再考虑到信号传播经电离层的延迟和大气 对流层的延迟，则（ 5-5）式改写为1 2 c tk c t

5、5.2.2 伪距定位的观测方程2(X S X)2(YSY)2(Zs Z)2(XsjX)2 (YsjY)2(Zsj2 1/2Z)2 1/2 c tkjjj12c t j5.3载波相位测量 利用测距码进行伪距测量是全球定位系统基本测距方法， 然 后由于测距码的码元长度较大， 对于一些高精度应用来说其 测距精度还显得多低无法满足需要。 如果观测精度均取至测 距码波长的百分之一，则伪距测量对 P 码而言量测精度为30cm,对C/A而言3m左右。而如果把载波作为量测信号 ，由 于载波的波长短， L1 19cm, L2 24cm ，所以 就可达到很高精度。 目前的大地型接收机的载波相位测量精 度一般为1-

6、2mm载波信号是一种周期性的正弦信号，而相 位测量又只能测定其不足一个波长的部分， 因而存在着整周 期不确定性的问题，使解算过程比较复杂。重建载波定义。GPS卫星采用L频带的两种不同频率的电磁波作为高频信 号，L1，L2。L1上调制C/A码、P码以及导航电文；L2载波 上仅调制P码与导航电文。GPS卫星发射信号的频率，受到 卫星上原子钟的基准频率的控制。 P 码采用基准频率， C/A 码仅取基准频率的 1/10 。而 L1 载波为基准频率的 154 倍， L2 载波取基准频率的 120 倍。调制信号，已调波。GPS卫星的测距码和数据码信号，又是怎样调制到载波上的 呢?GPS信号调制采用调相技术

7、实现的。 GPS的测距码信号和数据码信号，都是以二进制数为码元的时间序列，它具有信号 波形和信号序列两种表述形式。信号波形称为码状态，以符 号 u（t ）表示。信号序列通常以符号 u 表示，信号序列中 每一个元素取值为 0 或者 1，称为码值，并且约定，当码值 为 0 时，对应的码状态为 +1，当码值取 1 时，对应码状态为 -1. 阐述信号波形与信号序列的对应关系。实现码信号与载波信号的调制，只需取码状态与载波相乘就可以了。载波是一种电磁波，由于 GPS卫星上原子钟的振荡器产生，其数学表达式为一正弦波。 因此，当码状态 +1 与载波相乘时，显然不会改变载波的相位；当码状态取 -1 时与载波相

8、乘时， 载波相位改变 180 度。这样，当码值由 0 变为 1，或由 1 变 为 0 时，都会使调制后的载波相位改变 180 度，称为相位跃 迁。的角频率( j=1,2 )； i 为载波 Li 的初相( i=1,2 )。利用图来说明一正弦波加载测距码信号或数据码信号后，在码值由 0变为 1 或由 1 变为 0 的交替处，调制后的载波出现相位跃迁 码相关解调技术 平方解调技术5.3.1 载波相位测量原理载波相位测量的观测量是 GPS接收机所接收的卫星载波信号与接收机本振参考信号的相位差。在初始 to时刻，测得小于一周的相位差 0 ，其整周数为 N0j ，此时包含整周数的相位观测值为：k (t0)

9、 0 N0k (t0) k (t0) N0接收机连续跟踪卫星信号，不断测定小于一周的相位差 i ， 并利 用 整 周 计数 器 记 录 从 t0 到 ti 时 间 内的 整 周 数 变化 量 Int(),只要卫星Sj从to到t之间卫星的信号没有中断，贝y初始 时刻整周模糊度 N0j 就为一个常数，任意一时刻 ti 卫星 Sj 到 k接收机的相位差为：kj (ti ) k(ti) kj (ti ) N0j Int( ) (5-11 )532载波相位测量的观测方程 载波相位观测量是接收机(天线)和卫星位置的函数，只有 得到了它们之间的函数关系，才能从观测量中求解接收机 (或卫星)的位置。设在GPS

10、标准时刻Ta (卫星钟面时刻ta )卫星Sj发射的 载波信号相位为 (ta)，经过传播延迟 后，在GPS标准时刻 Tb (接收机钟面时刻tb)到达接收机。(tb) j (ta)考虑到卫星钟差和接收机钟差有(5-12 )Ta ta ta,Tb 如则有(Tb tb) (Ta ta)对于卫星钟和接收机钟，其振荡器频率一般稳定良好，所以其信号的相位与频率的关系可表示为：（t t） （t） f t （ 5-13）f为信号的频率，t为微小的时间间隔， 以2为单位。设fj为j卫星发射的载波频率，fi为接收机本振产生的固定 参考频率，且=fj=f，同时考虑到Tb Ta ，则有：（Tb） j（Ta） f （ 5

11、-14）顾及（5-13 ）和（5-14 ）两式，（5-12 ）式可以改写为：(5-15 )(Tb) f tb j(Ta) f taf tb f ta传播延迟 中考虑到电离层和对流层的影响 1和2则:1（ 1 2） （ 5-16）c（5-16 ）代入到（5-15 ）中：1 2) f ta f tb(5-17 )考虑至U （ 5-11 ）式，即顾及载波相位整周数Nk N0 Int()后，有k f ta f tb 1 2 Nkc c c（5-18 ）（5-18 ）即为接收机k对卫星j的载波相位测量的观测方程。533整周未知数No的确定载波相位相对定位具有很好的测量精度，是目前 GPS定位的主要方法，

12、但是这种高精度是以正确求定整周未知数 N0和彻底消除周跳为前提。因为无论整周未知数确定得不正确， 还是周跳没有消除干净，一个整周数值的错误，就可以产生 0.2m的误差，所以，整周未知数的确定、 整周跳变的探测与 消除，在利用GPS载波相位进行精密定位中，具有非常重要 的意义。1.伪距法伪距法是在进行载波相位测量的同时又进行了伪距测量， 将伪距观测值减去载波相位测量的实际观测值 （化为以距离为单位）后即可得到 No。但由于伪距测量的精度较低，所以有较多的 N o取平均值才能获得正确的整波段数。2.整周未知数当做平差的待定参数一经典方法把整周未知数作为基线向量平差计算中的待定参数， 平差过程中与其

13、他参数一起求解确定。静态相对定位中常用方法。（1）整数解（固定解）整周未知数从理论上讲应该是一个整数，但是，由于各种误差的影响，平差求得的整周未知数往往不是一个整数，而是 一个实数。对于短基线，当进行 1h 以上的静态相对定位，由于测站间 星历误差、大气折射等误差具有较强相关性，相对定位可以 使这些误差大大削弱；同时也由于较长的观测时间，观测卫 星的几何分布 会产生较大变化，因此能以较高的精度来求 定整周未知数。此时，平差求出的整周未知数一般为较接近 于邻近整数的实数，且如果整周未知数估值的中误差甚小， 则可直接取邻近的整数为整周未知数由整周未知数的整数解获得的基线解也称为固定解。 对于短基线

14、，由于这种方法顾及了整周未知数的整数特性， 因此能够改善相对定位的精度。（2）实数解（浮动解） 对于长基线，误差的相关性降低，因此卫星星历、大气折射 等误差的影响难以消除，求解的整周未知数精度较低。事实 上，整周未知数的实数解中往往包含了一些系统误差。 此时， 在将其取为某一整数， 实际上对于相对定位精度只会有损无益。所以通常对于 20Km以上的长基线通常不在考虑整周未 知数的整数性质，直接将实数作为整周未知数的解。由实数 整周未知数获得的基线解，也称为浮动解。3.三差法（多普勒法） 由于连续跟踪的所有载波相位观测量观测值中均含有相同 的整周未知数 N0 ，所以将相邻两个观测历元的载波相位相

15、减，就将该未知参数消去，从而直接解出坐标参数。这就是 多普勒法。 但两个观测历元之间的载波相位观测值之差受到 此期间接收机及卫星钟的随机误差的影响，所以精度不好， 往往用来解算未知数参数的初值。 三差法可以消除掉许多误 差，所以使用较广泛。4.快速确定整周未知数法 基本思想：利用初始平差的解向量（接收机点的坐标及整周 未知数的实数解）及其精度信息（单位权中误差和方差协方 差阵），以数理统计理论的参数估计和统计假设检验为基础， 确定在某一置信区间整周未知数可能的整数解组合， 然后依 次将整周未知数的每一组合作为已知值，重复进行平差计 算。其中使估值的验后方差或方差和为最小的一组整周未知 数，即为

16、整周未知数的最佳估值。从统计检验的角度，取整周未知数估值加上 3 倍的中误差 （ Nr 3 Nr ）为整周未知数的整数取值范围， 该范围内包含 的所有整数均作为整周未知数的候选值。当所有的整周未知数都取了整数后，可作为已知值代 入观测方程，再进行最小二乘平差求基线解。如果整周未知 数整数侯选值不止一个， 则应将所有卫星的候选值构成不同 组合，逐一代入进行平差计算，最后取能使基线解方差最小 的那一组整数作为整周未知数。5.4整周跳变的修复 计数器累计整周数 周跳（定义） 整周跳变的探测与修复： 探测出何时发生了周跳并求出丢失 的周跳。何时有意义，何时无意义。5.4.1屏幕扫描法 作业人员在计算机

17、屏幕前依次对每个站、每个时段、每个卫 星的相位观测值变化率的图像进行逐段检查， 观测其变化率 是否连续。如果出现不规则的突然变化时，就说明在相应的 相位观测中出现了整周跳变现象。 然后用手工编辑的方法逐 点、逐段修复。5.4.2高次差或多项式拟合法 此方法是根据周跳现象的发生将会破坏载波相位测量的观 测值 Int（ ） 随时间而有规律变化的特殊性来探测。 GPS 卫星的径向速度最大可达 0.9km/s（ 一般指物体运动速度在观察者视线方向的速度 分量，即速矢量在视线方向的投影。 ） 。因而整周计数每秒 钟可变化千周。因此 15S 输出一个观测值，相邻观测值间的 差值可达数万周，那么对于十几周的

18、跳变不易发现。但是如 果在相邻的观测值间求差，这些变化值就小很多，四次差、 五次差接近 0. 由于时钟振荡器的随机误差给相邻 L1 载波相 位造成影响为 2.4 周，所以采用求差法一般难以探测出只有几周的小周跳。 利用几个相位值拟合一个 n 阶多项式， 据此多项式来估计下 一个观测值并与实际观测值比较，从而发现周跳并修复。5.4.3在卫星间求差法 卫星间求差探测小周跳。 双差相位观测值高阶差来发现周跳。5.4.4用双频观测值修复周跳对 流层改正模型只与， 测站绝对温度，测站 气压，测站水气压，测站高程，对流层外边 缘高度。双频载波相位观测值求差法探测周跳 假设由测站 Ti 观测卫星 Sj ，相

19、位观测值为i ，测相伪距观测方程写为：i, (t)ij (t) ij (t) c ti(t) t j (t) Nij (t0) ij,f (t)几何距离，接收机钟差，卫星钟差，初始整 周模糊度，电离层，对流层。 对于观测历元 t 和 t+1 的相位观测值求差，如果不发生周跳，由于两个观测历元比较接 近，大气的变化不大，电离层和对流层的影 响可以消除，因此有：j1 i 丄2 2022 j令21i j(t 1) 协/(t 1)/(t) gt 1)ti(t) (j(t 1) j(t) : H(t 1)t(t) ( j(t 1 j(t) 对于双频观测值可得：ji,Li由于在相邻历元、两个频率的相位观测

20、值间求差，可以消除 电离层、对流层、卫星钟差、接收机钟差等的影响，剩余的 主要是相位观测误差，因此，检验 量可以看成服从正态分2L2v 0.023周，因此，可以根据双频载波相位观测值解算则可以判断历元t到t+1没有发生周跳，否则或L2发生周跳。理论上检验量 能探测出大于0.07周的周跳。如果Li，L2同时发生周跳，且周跳之比接近频率比时， 检验量V 0.07,双频相位求差法将无法检测到周跳。5.4.5根据平差后的残差发现和修复整周跳变经过上述处理的观测值中还可能存在一些未被发现的小周跳。修复后的观测值中也可能引入 1-2周的偏差。用这些观测值来进行平差计算，求得各观测值的残差。由于载波相位 观

21、测量的精度很高，因而这些残差的数值一般均很小。有周 跳的观测值上会出现很大的残差，据此可以发现和修复周 跳。5.5 GPS绝对定位与相对定位GPS色对定位也叫单点定位GPS相对定位5.5.1 静态绝对定位可 D/(t) c tij lj(t) T/(t)(伪距定位观测方程)上式分别为站星几何距离，卫星接收机相对钟差，电离层等 效距离误差，对流层等效距离误差其中Dij(t)为非线性项，表示测站与卫星间几何距离显然有 ：Dij (t) (Xj(t) Xi)2 (Yj(t) Yi)2 (Zj(t)乙)21/2这里X j(t), Yj(t),Zj(t)为t时刻卫星Sj的三维地心坐标，Xi ,Y,Zi则

22、是测站Ti的三维地心坐标.如果设Zi0(3-45 )电文提供的卫星瞬时坐标为固定值，那么，对Di（t）以式子中:Zi1 nJ;(ZJ(t)Zi0)卅(3-46)站星几何距离的线性表达式为：(DiJ(t)o (XJ(t) X0)2 (Y J(t) Y)2 (ZJ(t)乙)2)( 3-47)可得线性化的伪距观测方程为:Di j (t) (Dij (t)0 kij(t) Xi lij (t) Yi mij(t) Zic tij Iij(t) Ti j (t)( 3-48 ) 伪距法绝对定位解：在方程(3-48 )中，含有3个测站未知数 Xi，Y，Zj以及一个钟差未知数t,因此接收机至少要跟踪 4颗G

23、PS卫星，才能组成 4 个伪距观测方程， 由此推算出测站的三维地心坐标。 电离层和对流层延迟等效距离误差可通过适当的数学模型 解算。 (3-48) 式子可令 :Rij (t) Dij(t) Iij (t) Ti j (t) (3-49)并记 Di c tij (3-50)称为卫星钟与接收机钟相对钟差等效距离误差 . 于是 , 伪距 观测方程可以改写成为 :Ri j (t) (Dij(t)0 kij (t) Xi lij(t) Yi mij(t) Zi Di(3-51) 式中 j=1,2,3,4. 采用矩阵形式 , 则有上式可以简化为 :Ai (t) GiLi(t)4 4 4 14 1 (3-5

24、2)由此 , 伪距法绝对定位解可以表示为 :Gi41Ai(t)441 Li4(1t) (3-53)当卫星的颗数多于4 颗, nj4 , 可用最小二乘法组成误差方程 :vi(t) Ai (t) Gi Li(t) n j 1 n j 4 4 1 n j 1 (3-54)相应的最小二乘解为 :T 1 TGi (Ai (t)T Ai (t) 1 Ai (t)T Li (t)4 1 4 nj nj 4 4 nj n j 1 (3-55)对于 nt 个历元 , 可以形成 nt 组误差方程组 (3-55)Ti (心)施4 1 4 4 4 1 (3-57)未知数协因数矩阵则为且参数向量各分量的中误差 ：(mT

25、i )k oV(QOkk (3-59)式中,为伪距测量中误差,(Q)kk为主对角线上第k个元 素.当观测时间较长时，接收机钟差随时间变化不容忽视 .这时可用下述方法处理：一：是将钟差表示成多项式，平差时同时求出系数，另一种方法是不同的观测历元 ，分别取独立的钟差参数.卫星几何分布精度因子：影响GPS绝对定位的精度有两个：第一个：单位权中误差0, 它由码相关伪距测量的精度 ，卫星星历精度以及大气折射影 响等许多因素；另一个因素是未知数的协因数矩阵 QTi,它由卫星的空间几何分布确定.几何精度因子 DOP以此作为衡量卫星空间几何分布对定位 精度的标准.未知数的协因数矩阵为：qiiqi2qi3qi4

26、q21q22q23q24QTiq3iq32q33q34 (3-61)q41q42q43q44式子中各个元素反映出在特定的卫星空间几何分布下 ，不同参数的定位精度及其相关性信息 .因此利用这些元素的不同组合,即可定义出若干从不同侧面描述卫星空间几何分布对 定位精度影响的精度因子钟差精度因子TDOPTDOP= 、q44 (3-62)相应的中误差为mT oTDOP（3-63）三维位置精度因子 PDOPPDOP= Jqii q22 q33 (3-64)相应三维位置中误差为：mp 0 PDOP结合TDOP和 PDOP可定义反映卫星空间几何分布对接收机钟差和位置综合影响的精度因子 -GDOP:GDOP=相

27、应的时空精度中误差为mG 0 GDOP垂直分量精度因子 VDOPVDOP= q33 (3-68),相应表示卫星空间几何分布对接收机位置垂直分量的影响的垂直分量中误差为：mv 0 VDOP (3-69)VDOP的另一种定义称为高程精度因子 ：VDOP= Jr q/|r| (3-70)式子中，r x, y,z为测站概略位置向量；q qn.q22.q33为三 维精度因子向量.由(3-70)定义的VDOP计算mv,即得高程定 位中误差.水平分量精度因子 HDOPHDOP= qn q?2 (3-71)相应水平分量的中误差为：mH HDOP (3-72)HDOP也有另外一种定义称为水平位置精度因子HDOP

28、= PDOP2VDOP2 (3-73)由上式定义的HDOF计算mH ,即得水平位置中误差.从上式子中可看出,GPS绝对定位的误差和精度因子 (DOP)的 大小成正比，因此在精度确定的情况下，应尽量采用精度 因子小的一组卫星进行观测 .换句话说,当接收机跟踪的卫星多于4颗时,便可选择其中 GDOP最小的一组卫星进行观 测，这项工作称为选星，通常接收机可以自动完成。GDOP和体积的关系(70页)。5.5.2 静态相对定位相对定位是用两台接收机分别安置在基线的两端， 同步观测相同的GPS卫星，以确定基线端点的相对位置或基线向量。同样，多台接收机安置在若干条基线的端点，通过同步观测 GPS卫星可以确定

29、多条基线向量。在一个端点坐标已知的情 况下，可以用基线向量推求另一待定点的坐标。相对定位分为静态相对定位和动态相对定位之分。1.观测值的线性组合两个或多个观测站在同步观测相同卫星下，卫星的轨道误差、卫星钟差、接收机钟差以及电离层和对流层的折射误差 等对观测量的影响具有一定的相关性， 利用这些测量的不同 组合(求差)进行相对定位，可有效地消除或减弱相关误差 的影响，从而提高定位精度。卫星间求差，接收机间求差，不同历元间求差。 一次差：接收机间求一次差。SD1k2(ti)k2(ti ) 1k(ti ) (5-39 )SD1j2 (ti )2(ti ) 1 (ti ) (5-40 )二次差：在一次差

30、分观测值继续求差，所得的结果仍可以被当做虚拟观测值，叫做载波相位观测值的二次差或双差。常 用的求二次差是在接收机间求一次差后再在卫星间求二次 差，叫做星站二次差分DD1k2j SD1j2(ti ) SD1k2(ti) 2j (ti) 1j(ti) k2(ti ) 1k(ti) (5-41) 同样在 ti 1时刻，对 k、 j 卫星的测站间单差观测值求差也可 求得双差观测值。对二次差继续求差称为求三次差。三次差是在接收机、卫星 间和历元之间求三次差。TD1k2j (ti ,ti 1) DD1k2j (ti 1) DD1k2j (ti ) 单差观测值可以消除与卫星有关的载波相位及其钟差项， 双 差观测值中可以消除与接收机有关的载波相位及其钟差项，三差观测值可以消除卫星和接收机有关的初始整周模糊度项N。因而差分观测值模型是 GPS测量应用中广泛采用的平差模型，特别是双差观测值即星站二次差分模型，是大多数GPS基线向量处理软件包中必选的模型。观测方程线性化及平差模型载波相位观测方程的线性化：Dij(t) (Dij(t) Kj(