本发明属于电子、通信与自动控制领域，涉及到基于导航卫星的用户定位，特别涉及到一种多天线gnss差分定位方法。
背景技术：
传统室外gps单点定位精度一般为5-10m，已经不能够满足各种实际应用需求。鉴于测量误差在时间空间上的相关性，gnss(全球卫星导航系统)差分技术是通过建立基准站使得接收机获得相关信息对伪距或相位测量值进行差分。当前所使用的差分技术都需要设立差分站，需要较高的成本，并且接收机到差分站的距离远近限制了差分站的有效定位范围，在较大的区域内定位需要多个差分站，在人口稀少的地区内的用户数目也限制了差分站的建立。同时在一些应用场景比如海上、丘陵地带或是偏远地区，使用地面固定基准站进行差分定位会使基线向量(基准站与接收机相对位置)过长，同时某些地区由于其地理条件的限制可能不适合建造差分站。另外，可通过建立动态差分基准站来避免上述地面固定基准站存在的问题，但是这种方法仍要建设差分站，成本较高，同时应用范围较为狭窄。
技术实现要素：
为解决上述问题，本发明提供一种利用多天线进行差分定位的gnss导航定位方法，在没有地面固定基准站的情况下，实现高精度定位。
本发明的技术方案：
一种多天线gnss差分定位方法，na个天线同时对相同的卫星进行跟踪，接收机利用从na个天线接收到的信号进行载波相位测量，并对测量值进行差分运算；利用天线之间几何关系构建定位方程，通过相位差分值求解天线位置坐标，然后接收机再利用接收机中心点与天线的几何关系求解接收机位置坐标；具体步骤如下：
(1)构建系统模型
gnss导航定位系统包含一个接收机和na个天线a1，a2，…，an，其中，n＝1，2，…，na，na≥4，接收机与天线通过电缆相连；安装时，天线和接收机的位置不固定，根据应用场合改变天线与接收机之间的位置关系，构建一机多天线应用场景的系统模型；当接收机与天线安装成功后，得到接收机与天线之间的几何关系，包括各个天线到接收机的距离、各个天线之间的距离以及各个位置坐标之间的方向向量；其中，na个天线到接收机的距离分别为l1，l2，…，ln，单位m。
(2)测量值的差分处理
接收机的天线跟踪相同的ns颗卫星，同时对卫星信号载波相位进行测量，在接收机处对测量值进行差分处理以消除测量误差；根据na个天线的几何关系和差分后的载波相位测量值构建方程，具体过程如下：
对于天线a1，a2，…，an，接收机测得的关于卫星si信号的载波相位如公式(1)所示：
其中，
在同一时刻，对同一颗卫星si，以天线ar为参考，将经天线an和ar转发的载波相位测量值在接收机处进行差分，得到单差观测方程式，见公式(2)：
其中，(·)nr＝(·)n-(·)r，
由于天线之间的距离小且安装在同一个载体上，在同一时刻和同一地区的电离层与对流层误差的相关性高，
整理公式(3)，令
由于天线到卫星之间的距离远大于各天线之间的基线距离，所以在同一时刻，处于同一直线上的两天线与卫星之间连线是平行的；此时卫星到两个天线之间的距离差如公式(5)所示：
其中，bnr表示天线an与ar之间的距离，
将公式(5)代入公式(4)，得到公式(4)的等价形式，见公式(6)：
其中，n＝1，2，…，r-1,r+1，…，na；公式(6)表示经na-1个天线转发测得的相位值分别与经天线ar转发的相位值进行差分操作，经天线ar转发的相位值作为参考值；i＝1，2，…，ns，表示不同的卫星。
(3)构建定位方程组并计算天线位置坐标
通过周整模糊度确定算法来获取周整模糊度
(na-1)ns+na-1≥3na(7)
求解公式(7)，得到公式(8)：
ns≥2+1/(na-1)(8)
由于na≥4，则0＜1/(na-1)≤1/3，所以无论天线个数多少，使用3颗卫星就可以满足定位条件，将差分信息与天线几何关系所构建的方程进行联立，构建定位方程组，当r＝1时，定位方程组见公式(9)：
其中，
求解公式(9)，得到na个天线位置坐标。
(4)获得接收机位置坐标
根据设定的一机多天线应用场景的系统模型中的接收机与天线位置关系，通过天线位置坐标计算接收机位置坐标(xu,yu,zu)；根据天线与接收机几何位置关系，在以天线a1为原点的站心坐标系中，天线a1到接收机的方向向量为[δeδnδu]t，则在地心地固坐标系中对应的向量见公式(10)：
其中，坐标变换矩阵s的表达式见公式(11)：
其中，λ和
则接收机位置坐标见公式(12)：
其中，[x1y1z1]t是天线a1的位置坐标，[δxδyδz]t是天线a1到接收机的方向向量。
对于除a1外的其它天线，即当r取除1以外的值时，采用与r＝1同样的方式进行计算以获得接收机位置坐标。
本发明的方法适用于不同天线个数及天线与接收机位置关系。当天线个数发生变化和/或天线与接收机之间的几何位置关系发生变化时，可采用类同方法对天线处获得的载波相位测量值进行处理及构造定位方程。在天线个数大于4的情况下，在定位过程中，如果检测到某个天线的载噪比小于设定阈值，可以将此天线去除，利用剩余的天线进行定位，增加定位系统的鲁棒性。
本发明的有益效果：本发明提出的一种多天线gnss差分定位方法，可以在不借助基准站的情况下使用差分方法进行定位。相较于传统的差分定位法来说，节省了建立差分站的成本，减少了差分定位的复杂度，同时提供了较高的定位精度。本发明成本较低且适用范围广，并可获得较高的定位精度。系统结构和大小可根据实际情况进行调整，灵活性较高。
具体实施方式
以下结合技术方案，进一步说明本发明的具体实施方式。
一种多天线gnss差分定位方法，具体步骤如下：
(1)构建系统模型
gnss导航定位系统包含一个接收机和na个天线a1，a2，…，an，其中，n＝1，2，…，na，na≥4，接收机与天线通过电缆相连；安装时，天线和接收机的位置不固定，根据应用场合改变天线与接收机之间的位置关系，构建一机多天线应用场景的系统模型；当接收机与天线安装成功后，得到接收机与天线之间的几何关系，包括各个天线到接收机的距离、各个天线之间的距离以及各个位置坐标之间的方向向量；其中，na个天线到接收机的距离分别为l1，l2，…，ln，单位m。
(2)测量值的差分处理
接收机的天线跟踪相同的ns颗卫星，同时对卫星信号载波相位进行测量，在接收机处对测量值进行差分处理以消除测量误差；根据na个天线的几何关系和差分后的载波相位测量值构建方程，具体过程如下：
对于天线a1，a2，…，an，接收机测得的关于卫星si信号的载波相位如公式(1)所示：
其中，
在同一时刻，对同一颗卫星si，以天线ar为参考，将经天线an和ar转发的载波相位测量值在接收机处进行差分，得到单差观测方程式，见公式(2)：
其中，(·)nr＝(·)n-(·)r，
由于天线之间的距离小且安装在同一个载体上，在同一时刻和同一地区的电离层与对流层误差的相关性高，
整理公式(3)，令
由于天线到卫星之间的距离远大于各天线之间的基线距离，所以在同一时刻，处于同一直线上的两天线与卫星之间连线是平行的；此时卫星到两个天线之间的距离差如公式(5)所示：
其中，bnr表示天线an与ar之间的距离，
将公式(5)代入公式(4)，得到公式(4)的等价形式，见公式(6)：
其中，n＝1，2，…，r-1,r+1，…，na；公式(6)表示经na-1个天线转发测得的相位值分别与经天线ar转发的相位值进行差分操作，经天线ar转发的相位值作为参考值；i＝1，2，…，ns，表示不同的卫星。
(3)构建定位方程组并计算天线位置坐标
通过周整模糊度确定算法来获取周整模糊度
(na-1)ns+na-1≥3na(7)
求解公式(7)，得到公式(8)：
ns≥2+1/(na-1)(8)
由于na≥4，则0＜1/(na-1)≤1/3，所以无论天线个数多少，使用3颗卫星就可以满足定位条件，将差分信息与天线几何关系所构建的方程进行联立，构建定位方程组，当r＝1时，定位方程组见公式(9)(其它情况类似)：
其中，
求解公式(9)，得到na个天线位置坐标；
(4)获得接收机位置坐标
根据设定的一机多天线应用场景的系统模型中的接收机与天线位置关系，通过天线位置坐标计算接收机位置坐标(xu,yu,zu)；根据天线与接收机几何位置关系，可知在以天线a1为原点的站心坐标系中，天线a1到接收机的方向向量为[δeδnδu]t，则在地心地固坐标系中对应的向量见公式(10)：
其中，坐标变换矩阵s的表达式见公式(11)：
其中，λ和
则接收机位置坐标见公式(12)：
其中，[x1y1z1]t是天线a1的位置坐标，[δxδyδz]t是天线a1到接收机的方向向量。
对于除a1外的其它天线，即当r取除1以外的值时，采用与r＝1同样的方式进行计算以获得接收机位置坐标。