《GPS相对定位基本原理》由会员分享，可在线阅读，更多相关《GPS相对定位基本原理（31页珍藏版）》请在金锄头文库上搜索。
1、G P S 相 对 定 位 原 理1. 相对定位原理概述不论是测码伪距绝对定位还是测相伪距绝对定位，由于卫星星历误差、接收机钟与卫 星钟同步差、大气折射误差等各种误差的影响，导致其定位精度较低。虽然这些误差已作 了一定的处理，但是实践证明绝对定位的精度仍不能满足精密定位测量的需要。为了进一 步消除或减弱各种误差的影响，提高定位精度，一般采用相对定位法。相对定位，是用两台GPS接收机，分别安置在基线的两端，同步观测相同的卫星，通 过两测站同步采集GPS数据，经过数据处理以确定基线两端点的相对位置或基线向量（图 1-1）。这种方法可以推广到多台GPS接收机安置在若干条基线的端点，通过同步观测相同
2、的GPS卫星，以确定多条基线向量。相对定位中，需要多个测站中至少一个测站的坐标值 作为基准，利用观测出的基线向量，去求解出其它各站点的坐标值。图1-1 GPS相对定位在相对定位中，两个或多个观测站同步观测同组卫星的情况下，卫星的轨道误差、卫 星钟差、接收机钟差以及大气层延迟误差，对观测量的影响具有一定的相关性。利用这些 观测量的不同组合，按照测站、卫星、历元三种要素来求差，可以大大削弱有关误差的影 响，从而提高相对定位精度。根据定位过程中接收机所处的状态不同，相对定位可分为静态相对定位和动态相对定 位（或称差分GPS定位）。2. 静态相对定位原理设置在基线两端点的接收机相对于周围的参照物固定不
3、动，通过连续观测获得充分的 多余观测数据，解算基线向量，称为静态相对定位。静态相对定位，一般均采用测相伪距观测值作为基本观测量。测相伪距静态相对定位 是当前GPS定位中精度最高的一种方法。在测相伪距观测的数据处理中，为了可靠的确定 载波相位的整周未知数，静态相对定位一般需要较长的观测时间（1.0h3.0h）,称为经典 静态相对定位。可见，经典静态相对定位方法的测量效率较低，如何缩短观测时间，以提高作业效率 便成为广大GPS用户普遍关注的问题。理论与实践证明，在测相伪距观测中，首要问题是 如何快速而精确的确定整周未知数。在整周未知数确定的情况下，随着观测时间的延长， 相对定位的精度不会显着提高。
4、因此提高定位效率的关键是快速而可靠的确定整周未知数为此，美国的Remondi B.W提出了快速静态定位方法。其基本思路是先利用起始基线 确定初始整周模糊度（初始化），再利用一台GPS接收机在基准站T静止不动的对一组卫星0 进行连续的观测，而另一台接收机在基准站附近的多个站点 T 上流动，每到一个站点则停i 下来进行静态观测，以便确定流动站与基准站之间的相对位置，这种“走走停停”的方法 称为准动态相对定位。其观测效率比经典静态相对定位方法要高，但是流动站的GPS接收 机必须保持对观测卫星的连续跟踪，一旦发生失锁，便需要重新进行初始化工作。这里将 讨论静态相对定位的基本原理。2.1观测值的线性组合
5、假设安置在基线端点的GPS接收机T (i = 1,2 ),相对于卫星Sj和Sk，于历元t (i = 1,2 )进ii行同步观测(如图2-1)，则可获得以下独立的载波相位观测量：e j (t)，e)(t)，o k (t)，e k (t)，1 1 1 2 1 1 1 2ej(t )， ej(t )，2 1 2 2ek(t )，21ek(t22S)图2-1 GPS相对定位的观测量Sk (t )2T2Sj(t )1k2k2e j2)j (t)2)5 k 0- /1 2在静态相对定位中，利用这些观测量的不同组合求差进行相对定位，可以有效地消除这些观测量中包含的相关误差，提高相对定位精度。目前的求差方式有
6、三种：单差、双差、三差，定义如下： 单差(Single-Difference):不同观测站同步观测同一颗卫星所得观测量之差2-1)Ae j =e j (t )-0 j (t)21 双差(Double-Difference)：不同观测站同步观测同组卫星所得的观测量单差之差VA k (t )= k (t )- j (t)k (t)-k (t)21- j (t)- j (t)212-2)三差(Triple-Difference)：不同历元同步观测同组卫星所得的观测量双差之差8VAk (t)= VAk (t )-VAk (t )2-3)2.2观测方程2.2.1单差观测方程1- Ak (t )-A j
7、(t ) 11 j (t )- j (t )卩2 2 1 221=Ak (t )-A j (t )22=Tk (t )-k (t ) ”22 1 2 -T k (t )- k (t)图2-2单差示意图测相伪距观测方程为：九 j (t)=p j (t)+ c rs t (t )-8 tj (t)-九 Nj (t )+Aj (t )+Aj (t)(2-4)i i i i 0i,I pi ,T参见图2-2，将(2-4)式的测相伪距观测方程应用于测站T、T，并代入(2-1)式,12可得:九Mj (t)=p j (t)p j (t) + crt (t)-61 (t)九Nj (t) Nj (t)2 1 2
8、 1 2 1+ Aj (t)+ Aj (t) Taj (t) + Aj (t)_|2, I2,T1,I1,T令 At (t) = 61 (t)-61 (t), ANj = Nj (t)-Nj (t)2 1 2 1AAj(t)=Aj (t)Aj (t), AAj (t)=Aj (t)Aj (t)I2, I1,IT2,T1,T则单差观测方程可写为:九Aq j C)= p j C) p j C)+ cAtC)九AN j + AA j C)+ AAj C)21IT2-5)2-6)由（2-6）式可见：卫星的钟差影响可以消除。同时由于两测站相距较近（100km）, 同一卫星到两个测站的传播路径上的电离层、
9、对流层延迟误差的相近，取单差可进一步明 显的减弱大气延迟的影响。Sj(t)qj (t )1申j C)2Sk(t)q k C)2q k C)12.2.2双差观测方程图 2-3 双差示意图参见图2-3,两台GPS接收机安置在测站T、T，对卫星Sj的单差为AQ j （t），对卫星Sk12的单差为AQk (t)，则由(2-6)式，双差观测方程可表示为:M k (t )xvaq k (t )=(p k (t) p k (t )(p j (t )一 p j (t)-XVANj(27)2 1 2 1在上式中可见，接收机的钟差影响完全消除，大气折射残差取二次差可以略去不计。这是双差模型的突出优点。2.2.3三
10、差观测方程参见图21,分别以t和t两个观测历元，对上述的双差观测方程求三次差，可得三差1 2观测方程为壮 VAQ j(t)=(pk (t )-pk2 21-(pk (t )-pkL 2 11(t )-(p j (t )-p j (t M22 21 2 (t )-(p j (t )-p j (t )112 11128)从三差观测方程中可见，三差模型进一步消除了整周模糊度的影响。2.2.4准动态相对定位观测方程准动态相对定位方法是将一台GPS接收机固定在基准站不动，而另一台接收机在其周 围的观测站流动，在每个流动站静止观测几分钟，以确定流动站与基准站之间的相对位置。 准动态相对定位的数据处理是以载
11、波相位观测量为依据的，其中的整周未知数在初始化的 过程中已经预先解算出来。因此，准动态相对定位可以在非常短的时间内获得与经典静态 相对定位精度相当的定位结果。根据(2-4)式的测相伪距观测方程，若整周模糊度Nj (t )已经确定，将其移到等式左 i0端，则测相伪距观测方程可以写为Rj (t )=p j (t)+ c L t (t )-8 tj (t ) + Aj (t )+Aj (t)(2-9)i i ii ,Ii ,T式中：Rj (t)=九心(t)+XNj (t )。i i i 0若忽略大气折射残差影响，则上式求取站间单差观测方程可得ARj (t)= p j2(t) - pj (t) + c
12、At (t)2-10)若采用双差模型进行准动态相对定位，则由(2-9)式，再对卫星间取双差可得：VARk (t)= pk (t)- pk (t)- p j (t)+ p j (t)(2-11)2 12 12.3静态相对定位观测方程的线性化及平差模型为了求解测站之间的基线向量，首先就应该将观测方程线性化，然后列出相应的误差 方程式，应用最小二乘法平差原理求解观测站之间的基线向量。下面我们根据间接平差原 理来讨论载波相位观测量的不同线性组合的平差模型。假设，在协议地球坐标系中，观测站T的待定坐标近似值向量为i其改正数向量为q观测站T至卫星S j的测相伪距方程是非线性的，必须将其线性化。i2.3.1
13、单差模型取两个观测站T和T，其中T为基准站，其坐标已知。线性化的载波相位单差观测方1 2 1程：A j (t)=mj (t)nj (t)5 y + f At (t)ANj25 z2-12)AAj (t)+AAj (t)1+ Ip2o2式中，大气折射延迟误差的残差很小，忽略。于是相应的误差方程可写成如下形式：AV)j2l1 _九一一)j2j25 x 25 y25 z2-f At (t )+ANj +Alj(t)2-13)(tj1P-)j 20P上述情况是两观测站同时观测同一颗卫星Sj的情况，可以将其推广到两观测站于历元 t时刻同时观测数颗卫星的情况，设同步观测的卫星数为nj颗，则相应的方程组为：
14、或者写为v(t)= a(t)5 X + b(t)AN + c(t)At(t)+Al(t)(2-14)2若进一步考虑到观测的历元次数为n，则相应的误差方程为: tv(t)v (t1)2a(t )_ a (t1)25 X +b (t )1 b (t1)2AN +c (t )i00c (t )20 一0At(t )_At (t1)2+Al(t )_Al (t1)2上式可写v (t )a (t )2b (t)0 0 *c (t )At()Al (t )一n -t一n -t一n -t1- n 一t一n 一t一n 一tV = A5 X + BAN + CAt + L(2-15)2或者5 X+L2-16)、 2C J ANAt按最小二乘法求解：_5 X A_2ANBAtCP a-1A_C
