动态绝对定位原理
GPS绝对定位又叫单点定位,即以GPS卫星和用户接收机之间的距离观测值为基础,并根据卫星星历确定的卫星瞬时坐标,直接确定用户接收机天线在WGS-84坐标系中相对于坐标原点(地球质心)的绝对位置。
根据用户接收机天线所处的状态不同,绝对定位又可分为动态绝对定位和静态绝对定位。因为受到卫星轨道误差、钟差以及信号传播误差等因素的影响,静态绝对定位的精度约为米级,而动态绝对定位的精度约为10~40m。因此静态绝对定位主要用于大地测量,而动态绝对定位只能用于一般性的导航定位中。
将GPS用户接收机安装在载体上,并处于动态情况下,确定载体的瞬时绝对位置的定位方法,称为动态绝对定位。一般,动态绝对定位只能获得很少或者没有多余观测量的实时解,因而定位精度不是很高,主要被广泛应用于飞机、船舶、陆地车辆等运动载体的导航。另外在航空物探和卫星遥感领域也有着广阔的应用前景。
定位解算中,根据采用的距离观测原理的不同,可以分为测码伪距动态绝对定位和测相伪距动态绝对定位。
一、测码伪距动态绝对定位 测码伪距观测方程
????0???l为了推导方便,取:
??x???c?t?????? (1)
?m?n???yk12?????z??R???????1???2 (2)
则测码伪距观测方程可写为
R???0???l??x???c?t (3)
?m?n???yk?????z??j
式中的电离层和对流层延迟改正数可从卫星发射的导航电文中获得,而卫星S在地球协议坐标系中的坐标也可从卫星星历中得到。
显然,式中在某个历元t只有测站Ti在协议地球坐标系中的坐标向量?x,y,z?和接收机钟的钟差?tk这4个未知参数,正是我们需要求解的。为此,至少需要建立4个类似的方程。所以,用户至少需要同步观测4颗卫星以便获得4个以上测码伪距观测方程。
在动态绝对定位的情况下,由于测站是运动的,所以获得的伪距观测量很少。但为了获得实时定位结果,必须至少同步观测4颗卫星。
假设GPS接收机在测站Ti于某一历元t同步观测j颗卫星(j?1,2,3,?n),则由(3)式可得:
1?R?1???0??l1?2??2??2R????0??l ????????????n??n??n???0????R?????lm1m2?mnn1n2?nn?1???x????1???y?? (4 ????z?????1?c?t??k?其误差方程为:
?v1??l1?2??2vlV???????????n??n??v????lm1m2?mnn1n2?nn1?1???x??R?1??0?????22??1???yR??0??AX?L (5) ?????????z???????nn??1?c?t??R??0????k??当同时跟踪卫星数刚好为4颗(即n?4)时,无多余观测量,此时上式中
AX?L?0 (6) 直接解此方程组得唯一定位解:
X??A?1L (7)
很明显,当同时观测卫星数多于4颗时,则观测量的个数超过待求参数的个数,此时要利用最小二乘法平差求解误差方程(5),得:
X??ATA解的精度为:
???1ATL (8)
mX??0qii (9) 式中 mT为解的中误差;?0为伪距测量中误差;qii为权系数阵QZ主对角线的相应元素,
QZ?AiTAi???1。
上述测码伪距绝对定位模型(7)、(8),已被广泛应用于实时动态单点定位。因为通过卫星星历中获得的卫星瞬时坐标是WGS-84坐标,因此解求得到的接收机位置坐标也是WGS-84坐标系中的坐标。
实际应用中,有时给定的近似坐标偏差较大,而且线性化过程中略去二次及二次以上项对平差结果也有影响,在解算过程中往往一次平差不能达到理想的解算结果,因此常常采用迭代法。
顺便要指出,这里在解算载体位置时,不是直接求出它的三维坐标,而是求各个坐标分量的修正分量,也就是给定用户的三维坐标初始值,而求解三维坐标的改正数。在解算运动载体的实时点位时,前一个点的点位坐标可作为后续点位的初始坐标值。
二、测相伪距动态绝对定位
载波相位方程(或测相伪距方程)的线性化形式:
??x????N?t?????0???l?m?n???y0?? (10) ????z??c?tk???1???2令
R?????t????1???2 (11)
代入(10)式,则测相伪距观测方程可写为
R???0???l??x????N?t??c?t (12)
?m?n???y0k?????z??由于测相伪距法中引入了另外的未知参数——整周未知数,因此,若和测码伪距法一样,
观测4颗卫星无法解算出测站的三维坐标。
假设GPS接收机在测站Ti于某一历元t同步观测n颗以上卫星(j?1,2,3,4,?n),则由(5-28)式可得误差方程组为
?v1??l1?2??2v??l?V????????nj??n??v????lm1m2?mnn1???1??x????n2?????1?y????c?tk?????????z?????nn???1??1?0???N1?t0???R?1??0?1?1???N2?t???R?2??2?0?0??????????????????n???njn?????NtR??01???0?0???????AX?B?T?CN?L (13)
可见,误差方程中的未知参数有:三个测站点坐标,一个接收机钟差,n个整周未知数。这样误差方程中总未知参数为4+n个,而观测方程的总数只有n个,如此则不可能实时求解。
如果在载体运动之前,GPS接收机在t0时刻锁定卫星S后,先保持载体静止,求出整周模糊度Njj(j?1,2,3,4,?n)。据前述分析,只要在初始历元t0之后的后续时间里?t0?,
没有发生卫星失锁现象,它们仍然是只与初始历元t0有关的常数,在载体运动过程中当成常数来处理。
则(13)式可写为
1?1???xi??R?1??0??N1?t0?????R?2??2??N2?t???1???y00?????? (14) ????z????????nnnj????1???c?tk???R??0??Nt0???v1?2vV?????n?v???l1??2???l?????n???l?m1m2?m1n1n2?nn?AX?L
这样,就与(5)式在形式上完全一致。此时,同步观测4颗以上卫星,就可得到(8)是完全一样的实时解,只是解方程过程中采用的是测相伪距观测值,因此定位解的精度较之测码伪距法要高。
值得注意的是,采用测相伪距动态绝对定位时,载体上的GPS接收机在运动之前必须初始化,而且运动过程中不能发生信号失锁,否则就无法实现实时定位。然而载体在运动过程中,要始终保持对所观测卫星的连续跟踪,目前在技术上尚有一定困难,一旦发生周跳,则须在动态条件下重新初始化。因此,在实时动态绝对定位中,寻找快速确定动态整周模糊度的方法是非常关键的问题。