1、第四章 GPS 定位的观测量和观测方程GPS 定位要解决两个问题：一是观测瞬间 GPS 卫星的空间位置。 解决方法：通过导航电文中的卫星星历来确定。 （前已述）二是测站点卫星之间的距离。 解决方法：通过测定信号传播时间计算。传播时间通过 GPS 的观测量计算的。4.1 GPS 的主要观测量主要观测量：码伪距： C1、 P1、P2载波相位：L1、L2多普勒 D1、D2即：1)L1 载波相位观测值（L1 ）2)L2 载波相位观测值（L2 ）3)调制在 L1 上的 C/A-code 伪距（C1） 4)调制在 L1 上的 P-code 伪距（P1 ） 5)调制在 L2 上的 P-code 伪距（P2
2、） 6) L1 载波 Dopple 观测值（D1 ）7) L2 载波 Dopple 观测值（D1 ）在 RINEX 2.10 中, 定义了下列观测值类型 :L1，L2:L1 和 L2 上的相位观测值;C1：采用 L1 上 C/A 码所测定的伪距;P1, P2:采用 L1、L2 上的 P 码所测定的伪距;D1 , D 2:L1 和 L2 上的多普勒频率 ;T1, T2:子午卫星的 150( T1) 和 400 MHz( T2) 信号上的多普勒积分;S1, S2: 接收机所给出的 L1、L2 相位观测值的原始信号强度或 SNR 值。在反欺骗( AS) 之下所采集的观测值将被转换为“L2 ”或“P2
3、”, 并将失锁指示符( 见表 9-2) 的第二位置 1。观测值的单位: 载波相位为周 , 伪距为 m, 多普勒为 Hz,子午卫星为周, SNR 等则与接收机有关。【例】载波GPS 使用两种载波：L1 载波：fL1=154f0=1575.42MHz,波长 1=19.032cm,L2 载波：fL2=120f0=1227.6MHz,波长 2=24.42cm 。载波的作用：1、 在无线电通信技术中，为了有效地传播信息，都是将频率较低的信号加载在频率较高的载波上，此过程称为调制。然后载波携带着有用信号传送出去，到达用户接收机。2、 载波也可以用作测距信号来使用。3、 选择这两个载频，目的在于测量出或消除
4、掉由于电离层而引起的延迟误差。数据码数据码即为导航电文，D 码。测距码GPS 卫星的测距码属于伪随机噪声码，简称 PRN，是一个具有一定周期的取值 0和 1 的离散符号串。它不仅具有高斯噪声所有的良好的自相关特征，而且具有某种确定的编码规则，使我们人工能复制出来。C/A 码用于粗测距和捕获 GPS 卫星信号的伪随机码。选择 32 个码以 PRN1PRN32 命名各种 GPS 卫星。C/A 码的码元宽度较大。假设两个序列的码元对齐误差为码宽的 1/101/100，则此进相应的测距误差为 29.32.93m 。P 码是卫星的精测码，码率为 10.23MHz，码长约为 619x1012 比特。码元宽
5、度为 CA 码的 110，这时若取码元的对齐精度仍为码元宽度的l10l100，则由此引起的相应距离误差约为 293-029m，仅为 CA 码的110。所以 P 码可用于较精密的导航和定位，故通常也称之为精码。根据美国国际部规定，P 码是专为军用的。目前只有极少数高档次测地型接收机才能接收 P 码，且价格昂贵。即使如此，美国国防部又宣布实施 AS 政策，即在 P 码上增加一个极度保密的 W 码，形成新的 Y 码，绝对禁止非特许用户应用。载波、测距码和数据码共同构成了 GPS 卫星信号，如图所示。4.2 测码伪距观测方程伪距定义：含有接收机卫星钟的误差及大气传播误差的站星观测距离，故称为伪距。观测
6、量为时间延迟： =ti(GPS)-tj(GPS) 真值ji符号约定：tj(GPS)为卫星 sj 发射信号时的理想 GPS 时刻，ti(GPS)为接收机 Ti 收到该卫星信号时的理想 GPS 时刻，tj 为卫星 sj 发射信号时的卫星钟时刻，ti 为接收机 Ti 收到该卫星信号时的接收机钟时刻。tij 为卫星信号到达观测站的传播时间。tj 为卫星钟相对理想 GPS 时的钟差，ti 为接收机钟相对理想 GPS 时的钟差。则有tj=tj(GPS)+tj，t i=ti(GPS)+ti信号从卫星传播到观测站的时间为tij=ti-tj=ti(GPS)-tj(GPS)+ti-tj= jijit假设卫星至观测
7、站的几何距离为 ij，在忽略大气影响的情况下可得相应的伪距：当卫星钟与接收机钟严格同步时，上式所确定的伪距即为站星几何距离。通常 GPS 卫星的钟差可从卫星发播的导航电文中获得，经钟差改正后，各卫星之间的时间同步差可保持在 20ns 以内。如果忽略卫星钟差影响，并考虑电离层、对流层折射影响，可得测码伪距观测方程的常用形式测码伪距方程的线性化：上式中假设： 为卫星 于历元 t 的坐标近似值向量为观测站的坐标 近似值向量为卫星坐标的改正数向量为观测站坐标的改正数向量测站到卫星的方向余弦：所以：jijijijijiji tctct )()()()( tTtIttctt jigjijijiji ijj
8、itt)( 2122 )()()( ijijij ztytx0Xjsi iTjjjj zytxt )(,),()(TiiizyxX,)()(1)(00tltxt jiijjiji )(tmytyt jiijjiji )(100tnztzt jiijjiji )()(tlxtjiijmyjiij)()(tnztjiij2102020 )(ijijijji tzytxtt )()()()(,),(-)(0 tTtItcXtnmtlt jigjiiijijijijiji ijjijijijiji误差方程为：4.3 测相伪距观测方程载波相位观测值：载波相位应被称为载波拍频相位，它是收到的受多普勒频移影
9、响的卫星信号载波相位与接收机本机振荡产生信号相位之差。载波信号量测精度优于波长的 1/100，载波波长（L1=19cm,L2=24cm）比 C/A 码波长(C/A=293m)短得多，所以 GPS 测量采用载波相位观测值可以获得比伪距（C/A 码或P 码）定位高得多的成果精度。观测量为相位延迟： ijt=iti-jtj=Nij（t 0）+N ij（t-t 0）+ j（t ）i=Nij（t 0）+解释：N ij（t 0）为卫星被锁定后就确定不变的量，称为初始历元整周未知数，是一个未知的常数。Nij（t-t 0）为从初始历元到后续观测历元之间相位变化的整周数，由接收机自动连续记数，可测定。j（t）是
10、后续观测历元 t 时刻不足一周的小数部分，可测定。iNij（tt 0）+ j（t） ，是载波相位的实际观测值，即户接收机相位观测输出值。i测量点 P 到任一卫星在 ti 时刻相应测量有下式 :i=i+Ni+N0(1)式中,i 为在 ti 时刻测量点到任一卫星间相位测量的相位数值;i 为在 ti 时刻相应测量的小数值 (1 周内);Ni 为信号从 t0 时锁定到 ti 时刻接收到达的载波信号进行拍频计数的累计整周数;N0 为在 t0 时刻前载波相位从卫星到测量点在空间传播的整周数。i是瞬时可以测量到的值,Ni 是在相位锁定后可通过累计得到,唯 N0 是个未知数,这是无法通过观测得到的值,称初始整
11、周未知数,也称整周模糊度或多值性。 N0FrN0Int()iit0 ti)(ji)(ti )(-)(-)()(-)(,),(tV0ji tTtItttcXtnmtl jigjijijiiijijiji ）（同时， ijt=iti-jtj= iti(GPS)+ ti(ti)-jtj(GPS)+tj(ti)=iti(GPS)-jtj(GPS)+f ti(ti)-tj(ti)=f ij +f ti(ti)-tj(ti) （(t+ t)= (t)+f t，f j=fi =f）式中 ij= ijti(GPS) , tj(GPS)/c由于 tj(GPS) =ti(GPS) - ij ，将上式按级数展开得：
12、上式中二次项影响很小可忽略，并考虑接收机的钟差，可得以观测历元 ti 为根据的表达式：上式的计算可采用迭代法，并略去二次项如果顾及大气折射影响，则卫星信号的传播时间最终表达为代入上式：代入第一式，得：考虑关系式 =c/f，可得：上式中上标项对伪距的影响为米级。在相对定位中，如果基线较短（小于 20km） ，则有关项可忽略，载波相位观测方程和测相伪距观测方程分别为:.)(21)(1)(1 2jiijjiijijji GPStcGPStcPStc )()()( iijjiijijji ttt )(1)(1)( iijijijji tcttc )()(1)(1)(1)()(ijipjiiijijij
13、ji tTtIctcttct )()()()()(1)( ijipjiiijijijji tTtItttt )()()()(1)()()( ttIcftftcftctft jipjijijijijiji )()()(10tNtTtIcf jijipji jijijijiji 0tNtTtIjijipjiii )()()()()()( 0tTtIcftttftcft jipjijijijiji Njijijijijiji载波相位观测方程和测相伪距观测方程线性化：误差方程为：整周跳变修复整周跳变：卫星信号被障碍物挡住而暂时中断，或受无线电信号干扰造成失锁，计数器无法连续计数，当信号重新被跟踪后，使
14、整周计数不正确，但不到一整周的相位观测值仍是正确的。这种现象称为周跳。整周跳变的探测与修复常用的方法有下列几种：1.屏幕扫描法此种方法是由作业人员在计算机屏幕前依次对每个站、每个时段、每个卫星的相位观测值变化率的图像进行逐段检查，观测其变化率是否连续。如果出现不规则的突然变化时，就说明在相应的相位观测中出现了整周跳变现象。然后用手工编辑的方法逐点、逐段修复。2.用高次差或多项式拟合法此种方法是根据有周跳现象的发生将会破坏载波相位测量的观测值 Int()+ 随时间而有规律变化的特性来探测的。整周计数每秒钟可变化数千周，那么对于几十周的跳变就不易发现。但如果在相邻的两个观测值间依次求差而求得观测值
15、的一次差的话，这些一次差的变化就要小得多。在一次差的基础上再求二次差，三次差、四次差、五次差时，其变化就小的更多了。此时就能发现有周跳现象的时段来。四次、五次差已趋近于零。由于振荡器的随机误差而给相邻的 L1 载波相位造成的影响为 24 周，所以用求差的方法一般难以探测出只有几周的小周跳。通常也采用曲线拟合的方法进行计算。根据几个相位测量观测值拟合一个 n 阶多项式，据此多项式来预估下一个观测值并与实测值比较，从而来发现周跳并修正整周计数。3.在卫星间求差法在 GPS 测量中，每一瞬间要对多颗卫星进行观测，因而在每颗卫星的载波相位测量观测值中，所受到的接收机振荡器的随机误差的影响是相同的。在卫
16、星间求差后即可消除此项误差的影响。4.根据平差后的残差发现和修复整周跳变经过上述处理的观测值中还可能存在一些未被发现的小周跳。修复后的观测值中也可能引入 12 周的偏差。用这些观测值来进行平差计算，求得的各观测值的残差。由于载波相位测量的精度很高，因而这些残差的数值一般均很小。有周跳的观测值上则会出现很大的残差，据此可以发现和修复周跳。5.用双频观测值修复周跳采用双频载波相位观测值的组合，并考虑电离层折射改正有：)()()(-)()()( 00 tTtItctNtntmtltt jigjiijiijijijijiji X)(-)(-)()()(-)(,),(tV 00ji ttIttttctt
17、l jigjijijijiiijijiji ）（ 1-1 0f jijijiiijijijijiji 122121 fcAfNff 由于组合后右边只剩下整周数之差和电离层折射的残差项，由于电离层折射的残差项很小，故可用 来探测周跳。此法又称电离层残差法。整周未知数 N0 的确定确定整周未知数 N0 是载波相位测量的一项重要工作。常用的方法有下列几种：1伪距法伪距法是在进行载波相位测量的同时又进行了伪距测量，将伪距观测值减去载波相位测量的实际观测值(化为以距离为单位 )后即可得到 N0。但由于伪距测量的精度较低，所以要有较多的 N0 取平均值后才能获得正确的整波段数。2.将整周未知数当作平差中的
18、待定参数经典方法把整周未知数当作平差计算中的待定参数来加以估计和确定有两种方法。(1)整数解整周未知数从理论上讲应该是一个整数，利用这一特性能提高解的精度。短基线定位时一般采用这种方法。具体步骤如下：首先根据卫星位置和修复了周跳后的相位观测值进行平差计算，求得基线向量和整周未知数。由于各种误差的影响，解得的整周未知数往往不是一个整数，称为实数解。然后将其固定为整数(通常采用四舍五入法 )，并重新进行平差计算。在计算中整周未知数采用整周值并视为已知数，以求得基线向量的最后值。(2)实数解当基线较长时，误差的相关性将降低，许多误差消除得不够完善。所以无论是基线向量还是整周未知数，均无法估计得很准确
19、。在这种情况下再将未知数固定为某一整数往往无实际意义，所以通常将实数解作为最后解。采用经典方法解算整周未知数时，为了能正确求得这些参数，往往需要一个小时甚至更长的观测时间，从而影响了作业效率，所以只有在高精度定位领域中才应用。3.多普勒法(三差法)由于连续跟踪的所有载波相位测量观测值中均含有相同的整周未知数 No，所以将相邻两个观测历元的载波相位相减，就将该未知参数消去，从而直接解出坐标参数。这就是多普勒法。但两个历元之间的载波相位观测值之差受到此期间接收机钟及卫星钟的随机误差的影响，所以精度不太好，往往用来解算未知参数的初始值。三差法可以消除许多误差，所以使用较广泛。4.快速确定整周未知数法1990 年 E.Frei 和 G.Beutler 提出了利用快速模糊度（即整周未知数）解算法进行快速定位的方法。采用这种方法进行短基线定位时，利用双频接收机只需观测一分钟便能成功地确定整周未知数。
