GPS测量基本原理

1＞概述测量学中有测距交会确定点位的方法。与其相似，无线电导航定位系统、卫星激光测距定位系统，其定位原理也是利用测距交会的原理定位。就无线电导航定位来说，设想在地面上有三个无线电发射台，其坐标为已知，用户接收机在某一时刻采用无线电测距的方法分别测得了接收机至三个发射台的距离d1，d2,d3。只需以三个发射台为球心，以d1，d2,d3为半径作出三个定位球面，即可交会出用户接收机的空间位置。如果只有两个无线电发射台的话，则可根据用户接收机的概略位置交会出接收机的平面位置。这种无线电导航定位系统是迄今为止仍在使用的飞机船舶的的中导航定位方法。近代卫星大地测量中的卫星激光测距定位也是应用了测距交会定位的原理和方法。虽然用于测距的卫星（表面安装有激光反射镜）是在不停的运动中，但总可以利用固定于地面上三个已知点上的卫星激光测距仪同时测定某一时刻至卫星的距离d1，d2,d3,应用测距交会的原理便可确定该时刻卫星的空间位置。如此，可以确定三可以上卫星的空间位置。如果第四个地面点上（坐标未知）也有一台卫星测距仪同时参与了测定改点到三颗卫星的空间距离，则利用所测定的三个空间距离可交会出该地面点的空间位置。将无线电信号发射台从地面搬到卫星上，组成一颗卫星导航定位系统，应用无线电测距交会的原理，便可利用三个以上地面已知点（控制站）交会处卫星的位置，反之利用三颗以上的卫星的已知空间位置又可交会出地面未知点（用户接收机）的位置。这便是GPS卫星定位的基本原理。GPS卫星发射测距信号和导航电文，导航电文中含有卫星的位置信息。用户用GPS接收机在某一时刻同时接收三个以上的GPS卫星信号，测量出测站点（接收机天线中心）P至三颗以上GPS卫星的距离并解算出该时刻GPS卫星的空间位置坐标，据此利用距离交会法解算出测站P的位置坐标，如下图所示，设在时刻ti在在测站P用GPS接收机同时测出P点至三颗GPS卫星的距离P1，P2,P3，通过GPS电文解释出该时刻三颗GPS卫星的三维坐标分别为（Xi，Yi，Zi），j=l，2，3。用距离交会的方法求解出P点的三维坐标（X，Y，Z）的观测方程为在GPS定位中，GPS卫星是高速运动的卫星，其坐标随时间在快速变化着，需要实时地由GPS卫星信号测量测站至卫星之间的距离，实时地由卫星的导航电文解算出卫星的坐标值，并标定测站点的定位，依据测距的原理，其定位原理与方法主要由伪距法定位，载波相位测量定位以及差分GPS定位等，对于待定点来说，根据其运动状态可以将GPS定位分为静态定位和动态定位。静态定位指的是对于固定不动的待定点，将GPS接收机安置于其上，观测数分钟乃至更长的时间，以确定该点的三维坐标，又叫绝对定位，若以两台GPS接收机分别置于两个不同的固定不动的待定点上，则通过一定时间的观测，可以确定两个待定点之间的相对位置，又叫相对定位。而动态定位则至少有一台GPS接收机处于运动状态，测定的是各观测时刻（观测历元）运动中的接收机的点位（绝对点位或相对定位）利用接收到的卫星信号（测距码）或载波相位，均可进行静态定位。实际应用中，为了减弱卫星的轨道误差、卫星钟差、接收机钟差以及电离层和对流层的折射误差的影响常采用载波相位观测值的各种线性组合（即差分值）作为观测值，获得两点之间高精度的GPS基线向量（即坐标差）。2＞伪距测量伪距法定位是由GPS接收机在某一时刻测出得到四颗以上GPS卫星的伪距以及已知点的卫星位置，采用距离交会的方法求定接收机天线所在点的三维坐标。所测伪距就是由卫星发射的测距码信号到达接收机的传播时间乘以光速所得的量测距离。由于卫星钟、接收机钟的误差以及无线电信号经过电离层和对流层中的延迟，实际测出的距离P’与卫星到接收机的几何距离有一定的差距，因此一般称量测出的距离为伪距。用C/A码进行测量的伪距为C/A码伪距，用P码测量出来的伪距为P码伪距。伪距法定位虽然一次定位精度不高（P码定位误差约为10cm，C/A码定位误差为20~30m），但因其有定位速度快，且无多值性问题等优点，仍然是GPS定位系统进行导航的最基本的方法。同时，所测伪距又可以作为载波相位测量中解决整波数不确定性问题（模糊度）的辅助资料。因此，有必要了解伪距测量以及伪距法定位的基本原理。2.1伪距测量GPS卫星依据自己的时钟发出某一结构的测距码，该测距码经过t时间的传播后的到达接收机。接收机在自己的时钟控制下产生一组结构完全相同的测距码一一复制码，并通过时延器使其延迟时间T’将这两组测距码进行相关处理，若自相关系数R（T’）尹1，则继续调整延迟时间T’直至自相关系数R（t'）=1为止。使接收机所产生的复制码与接收到的GPS卫星测距码完全对齐，那么其延迟时间T’即为GPS卫星信号从卫星传播到接收机所用的时间t。GPS卫星信号的传播时一种无线电信号的传播，其速度等于光速c，卫星至接收机的距离即为T’与c的乘积。为什么采用码相关技术来确定伪距呢？GPS卫星发射出的测距码是按照某一规律排列的，在一个周期内每个码对应某一特定的时间。应该说识别出每一个码的形状特征，即用每个码的某一标志推算出时延值T进行伪距测量。但实际上每个码在产生过程中都带有随机误差，并且信号经过长距离传播后也会产生在自相关系数R(t’)=MAX的情况下来确定信号的传播时间t。这样就排除了随机误差的影响，实质上就是采用了多个码特征来确定T的方法。由于测距码和复制码在产生的过程中均不可避免地带有误差，而且测距码在传播过程中还会由于各种外界干扰而产生形变，因而自相关系数往往不可避免地带有误差，而其自相关系数不可能达到“1”，只能在自相关系数为最大的情况下来确定伪距，也就是本地码和接收码基本上对齐了。这样可以最大幅度地消除各种随机误差的影响，以达到提高精度的目的。测定自相关系数R(t’)的工作由接收机锁相环路中的相关器和积分器来完成。如下图由卫星钟控制的测距码a(t)在GPS时间t时刻自卫星天线发出，经传播延迟T到达GPS接收机，接收机所接收到的信号为a(t-T)。由接收机钟控制的本地码发生器产生一个与卫星发播相同的本地码a?(t+At)，At为接收机钟与卫星钟的钟差。经过码移位电路将本地码延迟t’，送至相关器与所接收到的卫星发播信号进行相关运算，经过积分器后，即可得到自相关系数R(t’)输出：伪距测量原理图调整本地码延迟T',可使相关输出达到最大值'R(t)=R(t)< maxt-T=t+△t—1(1)可得T'=T+△t+nT1一 .:p=p+cAt+n人⑵式中：p'为伪距测量值，「为卫星到接收机的几何距离，T为测距码的周期，入二cT为相应测距码的波长，n=0,1,2,3……是整数值，c为信号传播速度。式(2)中即为伪距测量的基本方程。式中n入称为测距模糊度。如果已知待测距离小于测距码的波长(如用P码测距)，则n=0,具有p'=P+cAt(3)称为无模糊度测距。由式（3）可知，伪距观测值P’的待测距离与钟差等效距离之和。钟差△《包含接收机钟差8tkK与卫星钟差如，即At=-5tK+5tj，若考虑到信号传播经电离层和大气对流层的延迟，则（3）式改写为：p—p'+8p+Sp+cSt-cStj8p..一...、一一.、一（4）式即为所测伪距与真正的几何距离之间的关系式。式中Sp1， 2分别为电离层和对流层的改正项。8tK的下标表示接收机号，5tj的上标j表示卫星号。2.2伪距定位观测方程从（4）式中可以看出，电离层和对流层改正可以按照一定的模型进行计算，卫星钟差atj可以自导航电文中取得。而几何距离P与卫星坐标（X，Y，Z）与接收机坐标（X，Y，Z）之间有如下关系：p2—（X-X）2+（Y-Y）2+（Z-Z）2s s s （5）式中，卫星坐标可根据导航电文求得，所以式中只包含接收机坐标三个未知数。如果将接收机钟差8tK也作为未知数，则共有四个未知数，接收机必须同时至少测定四颗卫星的距离才能解算出接收机的三维坐标值。为此，将（5）式代入（4）式，有：—p'j+Sp/+Spj一cStj式中，j=1,2,3……。式（6）即为伪距定位的观测方程组。3＞载波相位测量利用测距码进行伪距测量是全球定位系统的基本测距方法。然而由于测距码的码元长度较大，对于一些高精度应用来讲其测距精度还显得过低无法满足需要。如果观测精度均取至测距码波长的百分之一，则伪距测量对P码而言量测精度为30cm，对C/A码而言为3m左右。而如果把载波作为量测信号，由于载波的波长短，X=19cm，久=24cm，所以就可以达到很高的精度。目前的大地型L1 L2接收机的载波相位测量精度一般为1〜2mm，有的精度更高。但载波信号是一种周期性的正弦信号，而相位相关有只能测定不足一个波长的部分，因而存在着整周期数不确定的问题，使解算过程变得比较复杂。仅供个人学习参考在GPS信号中由于己用相位调整的方法在载波上调制了测距码和导航电文，因而接收到的载波的相位一不在连续，所以在进行载波相位测量以前，首先要进行解调工作，设法将调制在载波上的测距码和卫星电文去掉，重新获得载波，这一工作称为重建载波，重建载波一般可以采用两种方法，一种是码相关法，另一种是平方法。采用前者，用户可同时提取测距信号和卫星电文，但用户必须知道测距码的结构；采用后者，用户无须知道测距码的结构，但只能获得载波信号而无法获得测距码和导航电文。3.1载波相位测距原理载波相位测量的观测量是GPS接收机所接收的卫星载波信号，与接收机本振参考信号的相位差相位差，以w](七)表示k接收机在接收机钟面时刻'时所接收到的j卫星载波信号的相位值，vk(七)表示k接收机在面钟时刻\时所产生的本地参考信号的相位值，则k接收机在接收机钟面时刻'k时观测j卫星所取得的相位观测可写为wj(t)=w(t)—wj(t)kkkkkk(7)通常的相位或者相位差测量只是测出一周以内的相位值。实际观测中，如果对整周进行计数，则自某一初始取样时刻'0以后就可以取得的相位观测值。如下图，在初始'0时刻，测得小于一周的相位差为Av0，其整周数为N0，此时包含整周数的相位观测值应为：=wj+w(t)+Njkk0 0(8)接收机继续跟