摘  要： 为了满足VANET(Vehicular Ad-hoc Network)中对车辆位置信息的高精度要求，本文研究了一种基于全球定位系统(Global Positioning System，GPS)伪距双差和有高精度里程仪的惯性导航(Inertial Navigation Systems，INS)的协作相对定位方法。通过数据融合技术，将GPS伪距双差、GPS信号的多普勒频移以及被高精度里程仪修正后的INS加速度等数据进行融合处理，获得具有良好精度的相对定位结果。结果表明，使用该方法的定位性能优于无里程仪和INS的定位性能。

关键词： GPS；惯性导航；里程仪；相对定位

0 引言

为满足现代车辆安全行驶以及车辆定位等应用对定位精度的高要求，迫切需要更高水平更高精度的车辆定位技术的研究和使用。在现有全球卫星导航系统[1-2](Global Navigation Satellite Systems，GNSS)与定位技术的基础上，针对原有定位技术存在定位精度不高，以及在GPS伪距双差相对定位时，GPS信号存在短时间中断的问题，本文研究了一种基于全球定位系统(GPS)伪距双差和有高精度里程仪的惯性导航(INS)的协作相对定位技术。该技术在原有GPS伪距双差相对定位技术的基础上，添加有高精度的里程仪惯性导航系统INS的测量器件。INS测量器件用来补充GPS短时间中断时的车辆位置信息，高精度的里程仪用来修正INS的积累误差，然后将GPS伪距双差、GPS信号的多普勒频移与被高精度里程仪修正后的INS加速度等数据信息进行融合分析，从而可以获得更高精度的车辆相对定位。

1 车辆协作相对定位方法

本文所研究的是有高精度里程仪的INS/GPS伪距双差的协作相对定位。先假设研究的车辆上已装有GPS接收机、INS设备以及高精度的里程仪，同时目标车辆所在的区域GPS信号覆盖良好，接收机可以观测到5颗以上的可见卫星，并且车辆之间能够进行车与车(Vehicle to Vehicle，V2V)通信，共享数据信息[3]。

1.1 现有的GPS伪距双差的相对定位方法

在t时刻，节点k处的GPS接收机与卫星i之间伪距观测值[2]为：

其中，Rki为节点k处的接收机与卫星i之间的几何距离，c为光速，δk为节点k处接收机的时钟误差，δi为卫星i的时钟误差，εi，ion为电离层误差，εi，trop为对流层误差，εi，sate为卫星i的轨迹误差，

ki为接收机k处的热噪声与来自卫星i的信号多径误差。

由于卫星的时间差别和高空层等误差对处于同一区域的车辆来说近似相同，因此可以通过对同一观测卫星的不同的GPS接收机进行差分消除这些公共误差。如使用不同车辆k、l对同一卫星i进行观测，其伪距单差为：

其中，SDikl为伪距单差，

为车辆k和车辆l到卫星i的几何距离之差。由于卫星与车辆之间的距离远远大于车辆之间的距离，因此可以将车辆到卫星的方向向量看作是平行的。

可以近似表示为：

其中，

为车辆k指向卫星i的单位方向向量，rkl=rl-rk，表示为两车辆之间的相对位置向量，ri、rk、rl分别是卫星i、车辆k、l的位置矢量。

对车辆间的伪距单差进一步求差分，可得伪距双差

，双差的作用是消除接收机时间误差的影响。

伪距双差可以表示为：

可以由车辆自身的伪距以及V2V通信从邻近的车辆中获取数据计算得到。

是不能通过伪距双差进行消除的非公共误差，假设

均值为0、方差相等，在后面的性能分析中，将其作为观测噪声处理。采用伪距双差

方法，在有4颗以上的可见卫星时，可以通过最小二乘法获得车辆间相对位置rkl(导航系下的相对位置信息)的结果[4]。

1.2 改进的GPS伪距双差相对定位方法

1.2.1 车辆间基于载波频移的相对速度的测量

在t时刻，车辆k的GPS接收机与卫星i之间GPS的载波相位伪距[5]可以参考上节中的GPS伪距双差表示为[2]：

相位伪距对时间求导，可得接收机k从卫星i处接收到的信号的载波频率偏移[2]，并能从中提取多普勒频移：

其中，f是GPS信号的载波频率，GPS信号使用L1=    1 575.42 MHz，λ是GPS L1信号的波长，

是观测噪声。

车辆k与车辆l对于卫星i的载波频移单差为：

其中，vkl=vk-vl是指车辆k、l之间的相对速度。由于卫星移动(卫星移动速度约3 km/s)引起的式(4)和式(10)中单位向量的变化较小，并且这种变化与车辆到卫星距离(20 000 km)相比，可以忽略不计，因此，ei和ej可以直接采用任何一辆车辆的接收机所在的最终位置，然后结合卫星星历进行计算得到。vkl可以通过多普勒频移解算得到[6]。

1.2.2 车辆间基于INS的相对加速度的测量

车辆k、l之间的相对加速度矢量akl=ak-al，其中ak、al分别是车辆k、l的加速度矢量。两车之间的相对加速度akl、相对速度vkl、相对位移rkl以及GPS卫星的定位信息都定义在导航坐标系下。

INS在本设计中主要为车辆提供加速度以及方位角信息。欧拉角用来构建姿态矩阵Cbn，然后将INS提供的载体坐标系下的车辆加速度信息，根据ak=Cbnakb，将其从载体坐标系转换到导航坐标系[2]中。

再通过接收到的GPS信号的多普勒频移来估测欧拉角，假设?驻fki是车辆k接收到的卫星i的多普勒频移，根据式(9)可得：

其中，vi是卫星i的速度矢量，当有4颗以上可见卫星时，车辆在导航系下的速度vk可以根据上述多普勒频移关系式计算得到。

假设目标车辆在道路行驶时不发生路边滑动或者偏离地面，也即在载体坐标系下，沿Y轴、Z轴的速度分量为0，车辆仅沿X轴方向行驶，车辆k在载体坐标系下的速度大小为vkb=|vkb|，可以由INS中的IMU提供的加速度通过一次对时间的积分得到，载体加速度与导航加速度满足如下关系式[7]：

其中，g是重力加速度，

z、

y是INS测量器件提供的车辆在载体系下沿各轴旋转的角速度，通过公式vk=Cbn[vkb 0 0]T，可以求得偏航角

，进而得到姿态矩阵Cbn，最终实现加速度由载体系到导航系的变换。考虑惯性导航技术存在误差积累的问题，再使用参考文献[7]中的高精度里程仪与惯性导航组合的方式，修正INS的误差，提高其信息精确度。

1.2.3 融合数据信息的分析方法

卡尔曼滤波器是协作定位技术进行数据融合分析的关键，在本文中，使用卡尔曼滤波器将本地数据与接收到的数据进行融合分析时的状态模型[8]为：

其中，

为状态向量，F是滤波器的传输函数模型，G是噪声处理模型，

是沿着各轴的高斯相对加速度噪声，假设其均值为零、抽样间隔为τ。对于车辆k、l，状态向量、滤波器传输模型以及噪声处理模型分别为：

其中，0n是n×n的零矩阵，In是n维单位矩阵。假设对于车辆k、l，可见卫星为m颗，则观测模型为：

H是观测模型，

是观测噪声，考虑伪距双差和GPS信号的多普勒频移以及INS提供的加速度，可得如下的观测模型：

其中，

是INS观测的相对加速度向量转变成的导航系下的数值，εkl为噪声。通常将具有最高仰角的卫星定为1，称之为参考卫星，由式(5)、(10)、(16)可得：

其中，O′为(m-1)×3的零矩阵，假设各观测量之间相互独立，则观测噪声的协方差可以表示为：

假设

a2是INS加速度噪声的协方差，三个坐标轴之间相互独立，观测到的加速度的协方差为

a2I3，因Cbn相互正交[9]，则导航系下的加速度协方差也为

a2I3。相对加速度的协方差为

。

ρ2和

2分别是GPS伪距和多普勒频移的观测误差，根据参考文献[10]，可以得到

，其中：

A=[1(m-1)×1  -Im-1  -1(m-1)×1  Im-1](20)

其中1表示一个全为1的(m-1)×1矩阵。已知H、F、Q、G时，就可以使用卡尔曼滤波器实现相对定位。

2 仿真结果与分析

本文采用测量距离均方根drms来分析INS/GPS伪距双差和GPS伪距双差相对定位的性能。对车辆i相对定位最佳为：

。当车辆(即是目标车辆)为1时，

。在此假设参与协作的车辆数为5，GPS伪距测量观测误差的标准差为

ρ=3 m，可见卫星数目为4~7颗，仿真结果如图1、图2所示。

采用GPS/INS或GPS伪距双差两种方法进行相对定位时，相对定位距离均方根drms都随着卫星数目的增加而减小，也即可见卫星数目越多，定位性能越优。如图2为在GPS伪距双差基础上改进的GPS/INS伪距双差，对比图1中改进前的相对定位均方根误差drms，在高精度里程仪与INS器件的作用下，GPS伪距双差相对定位的drms有所减少，也即是采用本文方法优于GPS伪距双差的定位性能。

3 结论

本文主要介绍了有高精度里程仪的INS/GPS伪距双差的协作相对定位技术。在现有伪距双差定位技术的基础上，添加INS测量器件和高精度的里程仪，再通过数据融合分析处理技术，将被里程仪修正后INS的加速度信息、GPS信号的多普勒频移信息以及GPS伪距双差的相对定位信息进行融合分析处理，从而改进了GPS伪距双差的相对定位性能。在本文中，要求车辆GPS接收机可见卫星数不少于5颗，通过仿真对比，证明了使用INS后的伪距双差定位性能得到显著提高。

参考文献

[1] WELLENHOF B H， LICHTENEGGER H， WASLE E. 全球卫星导航系统 GPS， GLONASS， Galileo及其他系统[M].程鹏飞，蔡艳辉，文汉江，等，译.北京：测绘出版社，2009.

[2] KALPLAN E D， HEGARTY C J.GPS原理与应用[M].寇艳红，等，译.北京：电子工业出版社，2007.

[3] 常促宇，向勇，史美林.车载自组网的现状与发展[J].通信学报，2008，28(11)：116-126.

[4] JUNGLAS I A， WATSON R T. Location-based services[J]. Communications of the ACM， 2008， 51(3)：65-69.

[5] DISSANAYAKE G， SUKKARIEH S， NEBOT E， et al. The aiding of a low-cost strapdown inertial measurement unit using vehicle model constraints for land vehicle applications[J]. IEEE Transactions on Robotics and Automation， 2001， 17(5)：731-747.

[6] 任小伟，熊帅，杨文辉.卫星导航相对定位技术研究[J].全球定位系统，2014，39(2)：54-59.

[7] 胡彦亮.基于INS/GPS、INS/里程仪冗余的车载组合导航系统[D].哈尔滨：哈尔滨工程大学，2008.

[8] GREWAL M S， ANDREWS A P. Kalman filtering： theory and practice using Matlab(3rd edition)[M]. Hoboken：Wiley & Sons， Inc， 2008.

[9] AKEILA E， SALCIC Z， SWAIN A. Direct gravity estimation and compensation in strapdown INS applications[C]. Sensing Technology， 2008. IEEE 3rd International Conference on ICST 2008，2008：218-223.

[10] TEASLEY S P， HOOVER W M， JOHNSON C R. Differential GPS navigation[C]. PLANS′80-Po-sition Location and Navigation Symposium， 1980：9-16.