留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

一种利用载波相位差值的伪卫星定位方法

徐亚明 孙福余 张鹏 王金岭

徐亚明, 孙福余, 张鹏, 王金岭. 一种利用载波相位差值的伪卫星定位方法[J]. 武汉大学学报 ● 信息科学版, 2018, 43(10): 1445-1450. doi: 10.13203/j.whugis20170033
引用本文: 徐亚明, 孙福余, 张鹏, 王金岭. 一种利用载波相位差值的伪卫星定位方法[J]. 武汉大学学报 ● 信息科学版, 2018, 43(10): 1445-1450. doi: 10.13203/j.whugis20170033
XU Yaming, SUN Fuyu, ZHANG Peng, WANG Jinling. A Pseudolite Positioning Approach Utilizing Carrier Phase Difference[J]. Geomatics and Information Science of Wuhan University, 2018, 43(10): 1445-1450. doi: 10.13203/j.whugis20170033
Citation: XU Yaming, SUN Fuyu, ZHANG Peng, WANG Jinling. A Pseudolite Positioning Approach Utilizing Carrier Phase Difference[J]. Geomatics and Information Science of Wuhan University, 2018, 43(10): 1445-1450. doi: 10.13203/j.whugis20170033

一种利用载波相位差值的伪卫星定位方法

doi: 10.13203/j.whugis20170033
基金项目: 

国家自然科学基金 41374011

国家自然科学基金 41474005

湖北省自然科学基金 2014CFB717

详细信息
    作者简介:

    徐亚明, 博士, 教授, 主要从事精密工程测量和工业测量等领域的研究。ymxu@sgg.whu.edu.cn

    通讯作者: 张鹏, 博士, 讲师。pzhang@sgg.whu.edu.cn
  • 中图分类号: P228

A Pseudolite Positioning Approach Utilizing Carrier Phase Difference

Funds: 

The National Natural Science Foundation of China 41374011

The National Natural Science Foundation of China 41474005

the Natural Science Foundation of Hubei Province 2014CFB717

More Information
    Author Bio:

    XU Yaming, PhD, professor, specializes in the theories and methods of precise engineering surveying, industrial measurement. E-mail: ymxu@sgg.whu.edu.cn

    Corresponding author: ZHANG Peng, PhD, lecturer. E-mail: pzhang@sgg.whu.edu.cn
  • 摘要: 针对卫星定位技术无法应用于室内环境的情况,提出了一种新颖的利用载波相位差值的伪卫星定位方法,该方法能够实现亚米级的定位精度,而且无需基准站的支持,无需伪卫星之间的时钟同步,也不需要求解整周模糊度。构建了伪卫星定位的系统模型,并阐述了利用载波相位差值进行定位的基本原理。首先将双天线接收机输出的两组载波相位测量值进行单次差分操作,消除共有误差带来的影响,然后通过非线性最小二乘方法迭代解算出双天线连线中点的空间位置。仿真结果和基于双通道软件接收机的实测数据均证明了该方法的可行性,能够作为现有室内定位技术的有效补充。
  • 图  1  伪卫星和双天线的几何关系

    Figure  1.  Geometric Relations Between Pseudolite and Dual-Antenna

    图  2  天线方位角取不同值时的定位结果

    Figure  2.  Positioning Results when Antenna Azimuth Takes Different Values

    图  3  伪卫星天线布置

    Figure  3.  Distribution of Pseudolite Antennas

    图  4  4个测试点的定位结果

    Figure  4.  Positioning Results for Four Test Points

    表  1  5颗伪卫星的三维坐标

    Table  1.   Three-dimensional Coordinates of Five Pseudolites

    伪卫星号 伪卫星三维坐标/m
    x y z
    1 2.506 -3.186 3.204
    2 0.064 2.033 2.952
    3 5.124 1.930 3.128
    4 -2.358 -0.438 2.991
    5 -3.444 -2.871 3.255
    下载: 导出CSV

    表  2  天线方位角取不同值时的定位结果统计

    Table  2.   Statistical Positioning Results when Antenna Azimuth Takes Different Values

    测试点号 θ=0° θ=45° θ=90°
    x/m y/m z/m θ/(°) x/m y/m z/m θ/(°) x/m y/m z/m θ/(°)
    1 0.13 0.33 0.34 1.85 0.75 0.60 0.91 4.89 0.31 0.08 1.58 1.23
    6 0.12 0.34 0.39 1.87 0.12 0.34 0.39 1.87 0.30 0.07 1.56 1.08
    11 0.11 0.25 0.34 1.44 0.29 0.26 0.79 1.24 0.30 0.07 1.53 1.50
    16 0.11 0.59 0.36 2.37 0.11 0.59 0.36 2.37 0.41 0.10 1.62 1.34
    下载: 导出CSV

    表  3  4个测试点的定位结果统计

    Table  3.   Statistical Positioning Results for Four Test Points

    测试点
    颜色
    平均误差 最大误差
    x/m y/m x/m y/m
    0.15 0.29 0.34 1.16
    0.21 0.35 0.47 1.48
    绿 0.11 0.20 0.33 1.02
    0.23 0.22 1.33 0.39
    下载: 导出CSV
  • [1] Gu Y Y, Anthony L, Niemegeers I. A Survey of Indoor Positioning Systems for Wireless Personal Networks[J]. IEEE Communications Surveys & Tutorials, 200911(1):13-32 http://cn.bing.com/academic/profile?id=ac62f913c9565afdeb83c0f3768d99b3&encoded=0&v=paper_preview&mkt=zh-cn
    [2] Wang J. Pseudolite Applications in Positioning and Navigation Progress and Problems[J]. Journal of Global Positioning Systems, 2002, 1(1):48-56 doi:  10.5081/jgps
    [3] Rizos C, Barnes J, Gambale N. Experimental Results of Locata: A High Accuracy Indoor Positioning System[C]. International Conference on Indoor Positioning and Indoor Navigation, Zurich, Switzerland, 2010
    [4] Xu R, Chen W, Xu Y, et al. A New Indoor Positioning System Architecture Using GPS Signals[J]. Sensors, 2015, 15:10074-10087 doi:  10.3390/s150510074
    [5] Chongwon K, Hyoungmin S, Taikjin L, et al. A Pseudolite-Based Positioning System for Legacy GNSS Receivers[J]. Sensors, 2014, 14:6104-6123 doi:  10.3390/s140406104
    [6] Naohiko K, Shusuke M, Satoshi K, et al. Indoor and Outdoor Seamless Positioning Using Indoor Messaging System and GPS[C]. International Conference on Indoor Positioning and Indoor Navigation, Guimarães, Portugal, 2011
    [7] Daniele B, Ciro G. Asynchronous Pseudolite Navigation Using C/N0 Measurements[J]. Journal of Navigation, 2015, 69(3):639-658 http://cn.bing.com/academic/profile?id=fae2398b5e40f1eb17c3c88edbf2bcf4&encoded=0&v=paper_preview&mkt=zh-cn
    [8] Daniele B, Ciro G. Indoor Navigation Using Asynchronous Pseudolites[C]. 6th European Workshop on GNSS Signals and Signal Processing, Munich, Germany, 2013
    [9] Yoshihiro S, Hiroaki A, Takuji E, et al. Hyperbo-lic Positioning with Proximate Multi-channel Pseu-dolite for Indoor Localization[C]. Proceedings of the IGNSS Symposium 2013, Gold Coast, Australia, 2013
    [10] Samama N, Vervisch-Picois A, Taillandier-Loize T. A GNSS-like Indoor Positioning System Implementing an Inverted Radar Approach[C]. International Conference on Indoor Positioning and Indoor Navigation, Alcalá de Henares, Spain, 2016
  • [1] 柳景斌, 赵智博, 胡宁松, 黄格格, 龚晓东, 杨盛.  室内高精度定位技术总结与展望 . 武汉大学学报 ● 信息科学版, 2022, 47(7): 997-1008. doi: 10.13203/j.whugis20220029
    [2] 张星, 刘涛, 孙龙培, 李清泉, 方志祥.  一种视觉与惯性协同的室内多行人目标定位方法 . 武汉大学学报 ● 信息科学版, 2021, 46(5): 672-680. doi: 10.13203/j.whugis20200454
    [3] 甄杰, 吴建新, 刘纪平, 徐胜华, 周振发, 辛海强.  一种单基站高精度室内应急定位方法 . 武汉大学学报 ● 信息科学版, 2020, 45(8): 1146-1154. doi: 10.13203/j.whugis20200128
    [4] 程鹏飞, 文汉江, 刘焕玲, 董杰.  卫星大地测量学的研究现状及发展趋势 . 武汉大学学报 ● 信息科学版, 2019, 44(1): 48-54. doi: 10.13203/j.whugis20180356
    [5] 李团, 章红平, 牛小骥, 张全.  RTK/INS紧组合算法在卫星数不足情况下的性能分析 . 武汉大学学报 ● 信息科学版, 2018, 43(3): 478-484. doi: 10.13203/j.whugis20150305
    [6] 陈锐志, 叶锋.  基于Wi-Fi信道状态信息的室内定位技术现状综述 . 武汉大学学报 ● 信息科学版, 2018, 43(12): 2064-2070. doi: 10.13203/j.whugis20180176
    [7] 蔡成林, 王亮亮, 刘昌盛, 李刚, 秦懿, 邓洪高.  利用改进型TurboEdit算法进行BDS载波相位周跳检测与修复 . 武汉大学学报 ● 信息科学版, 2016, 41(12): 1632-1637. doi: 10.13203/j.whugis20140602
    [8] 施闯, 章红平, 辜声峰, 楼益栋, 唐卫明.  云定位技术及云定位服务平台 . 武汉大学学报 ● 信息科学版, 2015, 40(8): 995-999. doi: 10.13203/j.whugis20150118
    [9] 刘春燕, 王坚.  基于几何聚类指纹库的约束KNN室内定位模型 . 武汉大学学报 ● 信息科学版, 2014, 39(11): 1287-1292.
    [10] 温永智, 吴杰.  高轨飞行器精确导航的载波相位时间差分/捷联惯导紧组合算法 . 武汉大学学报 ● 信息科学版, 2011, 36(10): 1195-1199.
    [11] 韩天柱, 曹建平, 陆明泉.  基于抗干扰天线的远近效应抑制技术 . 武汉大学学报 ● 信息科学版, 2011, 36(10): 1222-1225.
    [12] 王劲, 宋茂忠.  GPS软件接收机弱信号捕获系统设计 . 武汉大学学报 ● 信息科学版, 2010, 35(7): 846-849.
    [13] 田辉, 夏林元, 莫志明, 吴东金.  泛在无线信号辅助的室内外无缝定位方法与关键技术 . 武汉大学学报 ● 信息科学版, 2009, 34(11): 1372-1376.
    [14] 蔡华, 赵齐乐, 楼益栋.  精密卫星钟差确定系统的实现与精度分析 . 武汉大学学报 ● 信息科学版, 2009, 34(11): 1293-1296.
    [15] 寇艳红, 周兴云.  利用数字中频信号模拟器和软件接收机对MBOC信号进行精度测试评估 . 武汉大学学报 ● 信息科学版, 2009, 34(11): 1377-1380.
    [16] 李伟强, 张其善, 杨东凯.  GNSS软件接收机中信号捕获新算法 . 武汉大学学报 ● 信息科学版, 2009, 34(9): 1130-1133.
    [17] 张飞舟, 杨东凯, 陈嘉, 杨伯钢.  GNSS软件接收机的实现与测试分析 . 武汉大学学报 ● 信息科学版, 2009, 34(5): 577-580.
    [18] 杜自钦.  GPS用于海上定位的模拟试验 . 武汉大学学报 ● 信息科学版, 1989, 14(1): 20-26.
    [19] 柳景斌, 王泽民, 吕轩凡, 李维, 尹斐, 仇宏煜.  一种低成本毫米波雷达的室内自定位方法 . 武汉大学学报 ● 信息科学版, 0, 0(0): -. doi: 10.13203/j.whugis20210593
    [20] 毕京学, 甄杰, 姚国标, 桑文刚, 宁一鹏, 郭秋英.  面向智能手机的改进有限状态机步态探测算法 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20200186
  • 加载中
图(4) / 表(3)
计量
  • 文章访问数:  1239
  • HTML全文浏览量:  95
  • PDF下载量:  390
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-10-09
  • 刊出日期:  2018-10-05

一种利用载波相位差值的伪卫星定位方法

doi: 10.13203/j.whugis20170033
    基金项目:

    国家自然科学基金 41374011

    国家自然科学基金 41474005

    湖北省自然科学基金 2014CFB717

    作者简介:

    徐亚明, 博士, 教授, 主要从事精密工程测量和工业测量等领域的研究。ymxu@sgg.whu.edu.cn

    通讯作者: 张鹏, 博士, 讲师。pzhang@sgg.whu.edu.cn
  • 中图分类号: P228

摘要: 针对卫星定位技术无法应用于室内环境的情况,提出了一种新颖的利用载波相位差值的伪卫星定位方法,该方法能够实现亚米级的定位精度,而且无需基准站的支持,无需伪卫星之间的时钟同步,也不需要求解整周模糊度。构建了伪卫星定位的系统模型,并阐述了利用载波相位差值进行定位的基本原理。首先将双天线接收机输出的两组载波相位测量值进行单次差分操作,消除共有误差带来的影响,然后通过非线性最小二乘方法迭代解算出双天线连线中点的空间位置。仿真结果和基于双通道软件接收机的实测数据均证明了该方法的可行性,能够作为现有室内定位技术的有效补充。

English Abstract

徐亚明, 孙福余, 张鹏, 王金岭. 一种利用载波相位差值的伪卫星定位方法[J]. 武汉大学学报 ● 信息科学版, 2018, 43(10): 1445-1450. doi: 10.13203/j.whugis20170033
引用本文: 徐亚明, 孙福余, 张鹏, 王金岭. 一种利用载波相位差值的伪卫星定位方法[J]. 武汉大学学报 ● 信息科学版, 2018, 43(10): 1445-1450. doi: 10.13203/j.whugis20170033
XU Yaming, SUN Fuyu, ZHANG Peng, WANG Jinling. A Pseudolite Positioning Approach Utilizing Carrier Phase Difference[J]. Geomatics and Information Science of Wuhan University, 2018, 43(10): 1445-1450. doi: 10.13203/j.whugis20170033
Citation: XU Yaming, SUN Fuyu, ZHANG Peng, WANG Jinling. A Pseudolite Positioning Approach Utilizing Carrier Phase Difference[J]. Geomatics and Information Science of Wuhan University, 2018, 43(10): 1445-1450. doi: 10.13203/j.whugis20170033
  • 卫星定位技术能够提供可靠和高精度的位置信息,但在室内环境中,卫星定位技术的使用受到了很大的限制,以至于无法持续输出精确和可靠的位置解。虽然WiFi、WSN、RFID和蓝牙等技术[1]都能够被用于室内定位,但是这些定位技术均无法持续地提供亚米级的定位精度。UWB技术尽管定位精度更高,但其占用频谱广、传送距离短的特点限制了UWB技术在定位领域的广泛应用。

    伪卫星定位系统为室内定位领域提供了一种合适的解决方案,其定位原理与传统的卫星系统基本相同,接收机的软硬件结构也与卫星接收机保持一致[2]。最有代表性的伪卫星系统是Locata系统[3],该系统充分利用了接收机输出的载波相位测量值,能够持续输出厘米级精度的位置解,但是Locata系统需要伪卫星保持严格的时间同步。文献[4-5]中的伪卫星定位系统可以利用传统的卫星接收机实现米级定位精度。但是伪距测量值容易受多路径效应的影响,在反射物较多的室内环境中,定位精度会明显降低。IMES系统[6]在伪卫星导航电文中包含了伪卫星的坐标和楼层信息,其定位精度与伪卫星的布设密度有着密切关系,定位精度一般只能达到10 m左右。文献[7-8]选择利用接收机输出的伪卫星信号载噪比,结合无线信号的衰落模型,计算出接收机的空间位置, 只能够达到10 m级的定位精度。

    本研究介绍了一种基于双天线接收机的伪卫星定位技术,双天线之间的距离固定且小于伪卫星信号载波波长的一半。首先将双天线接收机输出的两组载波相位测量值进行单次差分操作;然后通过非线性最小二乘方法解算出双天线连线中点的空间坐标, 并对定位结果进行了仿真分析,讨论了影响定位结果的主要因素;最后在由5颗商用伪卫星构成的真实伪卫星试验场景中,利用双通道软件接收机输出的实测数据,验证了这种利用载波相位差值的伪卫星定位方法能够达到亚米级的定位精度。

    与以Locata为代表的伪卫星定位系统相比,这种利用双天线接收机输出的载波相位差值的伪卫星定位方法无需基准站的支持,也不需要进行伪卫星间的时钟同步,既能够充分利用载波相位测量值高精度特点,又无需进行整周模糊度解算,因此可作为现有伪卫星定位技术的有效补充。

    文献[9-10]与本研究有相似之处。文献[9]提出了一种利用载波相位测量值的双曲线定位方法,采用3个近距离布设的伪卫星发射天线, 并且每个天线发射不同的伪随机码。随着接收机和伪卫星距离的增加,定位结果会显著恶化; 文献[10]进一步利用6个近距离布设的伪卫星发射信号的载波相位测量值来计算接收机的位置,实现更为复杂。而本研究适用于更加通用的分布式伪卫星场景,并且采用的双天线接收机已经在定向和抗干扰等领域得到了应用。

    • 因为研究的重点主要在室内区域,所以本研究选择了独立的本地坐标系。在进行室内定位时,通常将接收天线固定在待定位设备的顶部,而且室内环境的地面相对比较平整,因此认为两个距离接近的接收天线始终位于同一个水平面。

      图 1给出了伪卫星i、接收机天线ab的空间几何关系。天线ab处于同一个水平面,且天线之间的距离为d(d需要小于伪卫星信号载波波长的一半)。由于d的值固定不变,并且可以进行精确标定,因此d被视为已知量。

      图  1  伪卫星和双天线的几何关系

      Figure 1.  Geometric Relations Between Pseudolite and Dual-Antenna

      两个天线几何中心O的三维坐标表示为x, y, zxyz均为未知参数。天线连线与x坐标轴之间的角度表示为θ(以下简称天线方向角)。由于θ的值会随着待定位设备的转向而发生变化,因此也将θ作为未知参数进行求解。

      根据图 1的几何关系,可以推算出天线ab的三维坐标分别为(x-Δx, y-Δy, z)和(xx, yy, z),Δx和Δy分别表示为:

      $$ \Delta x = 0.5d\cos \theta $$ (1)
      $$ \Delta y = 0.5d\sin \theta $$ (2)

      伪卫星i的天线坐标已知,记为(xi, yi, zi),那么伪卫星i和天线a之间的距离可以表示为

      $$ \begin{array}{*{20}{c}} {r_a^i = }\\ {\sqrt {{{\left( {x - \Delta x - {x^i}} \right)}^2} + {{\left( {y - \Delta y - {y^i}} \right)}^2} + {{\left( {z - {z^i}} \right)}^2}} } \end{array} $$ (3)

      类似地,伪卫星i和天线b之间的距离也可以记为:

      $$ \begin{array}{*{20}{c}} {r_b^i = }\\ {\sqrt {{{\left( {x + \Delta x - {x^i}} \right)}^2} + {{\left( {y + \Delta y - {y^i}} \right)}^2} + {{\left( {z - {z^i}} \right)}^2}} } \end{array} $$ (4)

      rairbi的差值可以通过式(5)计算得到:

      $$ \begin{array}{*{20}{c}} {r_{ab}^i = r_a^i - r_b^i = }\\ {\sqrt {{{\left( {x - \Delta x - {x^i}} \right)}^2} + {{\left( {y - \Delta y - {y^i}} \right)}^2} + {{\left( {z - {z^i}} \right)}^2}} - }\\ {\sqrt {{{\left( {x + \Delta x - {x^i}} \right)}^2} + {{\left( {y + \Delta y - {y^i}} \right)}^2} + {{\left( {z - {z^i}} \right)}^2}} } \end{array} $$ (5)

      由式(1)和式(2)可知,Δx和Δy中的未知参数只有θ,因此在式(5)中,有xyzθ共4个未知参数,而rabi的值可以通过下文§2中式(9)计算得到。

    • 一旦天线a和天线b所连接的双天线接收机捕获到了伪卫星i的信号,就会借助锁相环技术,对伪卫星i的信号保持持续跟踪,并输出相应的载波相位测量值

      $$ \varphi _a^i = {\lambda ^{ - 1}}r_a^i + \frac{c}{\lambda }\left( {{\rm{ \mathsf{ δ} }}t - {\rm{ \mathsf{ δ} }}{t^i}} \right) + N_a^i + \varepsilon _{\varphi ,a}^i $$ (6)
      $$ \varphi _b^i = {\lambda ^{ - 1}}r_b^i + \frac{c}{\lambda }\left( {{\rm{ \mathsf{ δ} }}t - {\rm{ \mathsf{ δ} }}{t^i}} \right) + N_b^i + \varepsilon _{\varphi ,b}^i $$ (7)

      式中,λ表示伪卫星载波信号的波长;δtδti分别表示双天线软件接收机的钟差和伪卫星i的钟差, 由于天线a和天线b共享一个双天线接收机,共用同一个晶体振荡器,因此对应同一个δt; NaiNbi分别表示天线a和天线b捕获到的伪卫星i信号的整周模糊度;ε$φ$, aiε$φ$, bi可以被视为载波相位测量值的测量噪声。通过式(8)可以进一步计算$φ$ai和$φ$bi的差值,以移除未知的接收机钟差δt和伪卫星i的钟差δti

      $$ \begin{array}{*{20}{c}} {\varphi _{ab}^i \approx \varphi _a^i - \varphi _b^i = {\lambda ^{ - 1}}\left( {r_a^i - r_b^i} \right) + }\\ {\left( {N_a^i - N_b^i} \right) + \left( {\varepsilon _{\varphi ,a}^i - \varepsilon _{\varphi ,b}^i} \right)} \end{array} $$ (8)

      式(8)两端的单位均为周,因此将式(8)两端同时取余数运算。因为天线ab之间的距离被限制为小于伪卫星信号载波波长λ的一半,根据三角形原理,有(rairbi) < 0.5λ,即λ-1(rairbi)总是小于0.5周,所以取余数之后的结果仍然为λ-1(rairbi)。而(NaiNbi)始终为整数,取余数后的结果为0。又由于载波相位测量值的测量噪声ε$φ$, aiε$φ$, bi一般小于0.01周,因此(ε$φ$, aiε$φ$, bi)取余数的结果也仍然为(ε$φ$, aiε$φ$, bi)。则式(8)两端取余数的结果可以表示为式(9),其中mod表示取余数运算符。

      $$ \begin{array}{*{20}{c}} {\bmod \left( {\varphi _{ab}^i,1} \right) = {\lambda ^{ - 1}}\left( {r_a^i - r_b^i} \right) + }\\ {\left( {\varepsilon _{\varphi ,a}^i - \varepsilon _{\varphi ,b}^i} \right) = {\lambda ^{ - 1}}r_{ab}^i + \varepsilon _{\varphi ,ab}^i} \end{array} $$ (9)

      式(9)给出了$φ$abi的小数部分mod($φ$abi, 1)和rabi之间的关系,而mod($φ$abi, 1)可以根据载波相位测量值$φ$ai和$φ$bi直接计算得到。

      由于式(5)中有xyzθ共4个未知参数,记为u =[x y z θ]T,理论上最少需要4颗伪卫星的测量值mod($φ$abi, 1),i=1, 2, 3, 4。考虑到冗余量,在接下来的仿真和实测场景中,均选择了5颗伪卫星的载波相位测量值的数据,并通过高斯牛顿迭代算法进行定位计算。

      由5颗伪卫星的测量值rabi(i=1, 2…5)组成5组非线性方程。首先设定u =[x y z θ]T的初始值分别为x0y0z0θ0,对式(5)进行线性化:

      $$ \begin{array}{*{20}{c}} {r_{ab}^i\left( {x,y,z,\theta } \right) = r_{ab}^i\left( {{x_0},{y_0},{\theta _0}} \right) + }\\ {\frac{{\partial r_{ab}^i}}{{\partial x}}\left| {_{x = {x_0}}} \right.\Delta {x_0} + \frac{{\partial r_{ab}^i}}{{\partial y}}\left| {_{y = {y_0}}} \right.\Delta {y_0} + }\\ {\frac{{\partial r_{ab}^i}}{{\partial z}}\left| {_{z = {y_0}}} \right.\Delta {z_0} + \frac{{\partial r_{ab}^i}}{{\partial \theta }}\left| {_{\theta = {\theta _0}}} \right.\Delta {\theta _0}} \end{array} $$ (10)

      然后由5组线性化后的方程构成如下矩阵方程:

      $$ \mathit{\boldsymbol{G}} \cdot \Delta {\mathit{\boldsymbol{u}}_0} = \mathit{\boldsymbol{b}} $$ (11)
      $$ \begin{array}{*{20}{c}} {\mathit{\boldsymbol{G}} = }\\ {\left[ {\begin{array}{*{20}{c}} {\frac{{\partial r_{ab}^1}}{{\partial x}}\left| {_{x = {x_0}}} \right.}&{\frac{{\partial r_{ab}^1}}{{\partial y}}\left| {_{y = {y_0}}} \right.}&{\frac{{\partial r_{ab}^1}}{{\partial z}}\left| {_{z = {z_0}}} \right.}&{\frac{{\partial r_{ab}^1}}{{\partial \theta }}\left| {_{\theta = {\theta _0}}} \right.}\\ \vdots&\vdots&\vdots&\vdots \\ {\frac{{\partial r_{ab}^5}}{{\partial x}}\left| {_{x = {x_0}}} \right.}&{\frac{{\partial r_{ab}^5}}{{\partial y}}\left| {_{y = {y_0}}} \right.}&{\frac{{\partial r_{ab}^5}}{{\partial z}}\left| {_{z = {z_0}}} \right.}&{\frac{{\partial r_{ab}^5}}{{\partial \theta }}\left| {_{\theta = {\theta _0}}} \right.} \end{array}} \right]} \end{array} $$ (12)
      $$ \Delta {\mathit{\boldsymbol{u}}_0} = \left[ {\begin{array}{*{20}{c}} {\Delta {x_0}}\\ {\Delta {y_0}}\\ {\Delta {z_0}}\\ {\Delta {\theta _0}} \end{array}} \right] $$ (13)
      $$ \mathit{\boldsymbol{b}} = \left[ {\begin{array}{*{20}{c}} {\lambda \bmod \left( {\varphi _{ab}^1,1} \right) - r_{ab}^1\left( {{\mathit{\boldsymbol{u}}_0}} \right)}\\ \vdots \\ {\lambda \bmod \left( {\varphi _{ab}^5,1} \right) - r_{ab}^5\left( {{\mathit{\boldsymbol{u}}_0}} \right)} \end{array}} \right] $$ (14)

      直接利用最小二乘方法求解式(11),得到:

      $$ \Delta {\mathit{\boldsymbol{\hat u}}_0} = {\left( {{\mathit{\boldsymbol{G}}^{\rm{T}}}\mathit{\boldsymbol{G}}} \right)^{ - 1}}{\mathit{\boldsymbol{G}}^{\rm{T}}}\mathit{\boldsymbol{b}} $$ (15)

      然后进行迭代计算:

      $$ {{\mathit{\boldsymbol{\hat u}}}_1} = {\mathit{\boldsymbol{u}}_0} + \Delta {{\mathit{\boldsymbol{\hat u}}}_0} $$ (16)

      经过几次迭代后,就能够求解出未知参数的估计值$\mathit{\boldsymbol{\hat u = }}{\left[{\hat x\;\hat y\;\hat z\;\hat \theta } \right]^{\rm{T}}}$,从而完成对双天线接收机的定位计算。

    • 为了和真实的试验场景对应,仿真过程中选择了与真实伪卫星试验场景完全一致的伪卫星三维坐标。共包括5颗伪卫星,伪卫星的三维坐标见表 1

      表 1  5颗伪卫星的三维坐标

      Table 1.  Three-dimensional Coordinates of Five Pseudolites

      伪卫星号 伪卫星三维坐标/m
      x y z
      1 2.506 -3.186 3.204
      2 0.064 2.033 2.952
      3 5.124 1.930 3.128
      4 -2.358 -0.438 2.991
      5 -3.444 -2.871 3.255

      在本研究中,伪卫星采用的载波波长和GPS的L1频段的波长一致,均为0.19 m,双天线之间的间距也相应地设为0.095 m。载波相位测量值的精度可以达到载波波长的1%,因此理论上,载波相位测量值能够达到优于2 mm的测量精度。在仿真过程中,给载波相位测量值叠加了最大2 mm的随机测量噪声。仿真过程直接采用了MATLAB中的fsolve函数来进行未知参数的迭代计算,未知参数xyz的初始值均设置为0,而未知天线方位角θ的初始值则设置为在真实值的基础上叠加最大20°的随机误差。在实际的定位过程中,该方位角的初始估计值可由电子罗盘获取。

      在定位范围内选取了16个测试点进行仿真试验。测试点的真实平面位置在图 2中用黑色圆点表示,代表双天线连线的几何中心,16个测试点真实的z坐标均为0.2 m。5颗伪卫星的空间坐标用5个实心小方块表示。在每个测试点均利用叠加了随机测量噪声的载波相位测量值进行60次定位计算,计算出的平面定位结果用绿色点表示。

      图  2  天线方位角取不同值时的定位结果

      Figure 2.  Positioning Results when Antenna Azimuth Takes Different Values

      另外,因为天线方位角θ也需要作为未知参数进行估计,所以在仿真验证阶段,将θ的真实值分别设定为0°、45°和90°等3种情况。在3种情况下分别计算出的16个测试点的平面坐标值分别如图 2所示。

      选择1、6、11和16号4个典型的测试点进行误差分析,这4个测试点的平面坐标分别为(-0.6,-2.4)m、(0.6,-1.2)m、(1.8,0)m和(3.0,1.2)m。针对天线方位角分别为0°、45°、90°时3种情况,表 2给出了在4个测试点位置计算出的xyz方向的最大误差以及方位角θ的最大误差。

      表 2  天线方位角取不同值时的定位结果统计

      Table 2.  Statistical Positioning Results when Antenna Azimuth Takes Different Values

      测试点号 θ=0° θ=45° θ=90°
      x/m y/m z/m θ/(°) x/m y/m z/m θ/(°) x/m y/m z/m θ/(°)
      1 0.13 0.33 0.34 1.85 0.75 0.60 0.91 4.89 0.31 0.08 1.58 1.23
      6 0.12 0.34 0.39 1.87 0.12 0.34 0.39 1.87 0.30 0.07 1.56 1.08
      11 0.11 0.25 0.34 1.44 0.29 0.26 0.79 1.24 0.30 0.07 1.53 1.50
      16 0.11 0.59 0.36 2.37 0.11 0.59 0.36 2.37 0.41 0.10 1.62 1.34

      根据定位结果和误差分析可以得出:

      1) 采用载波相位差值的定位方法在平面方向的最大定位误差均小于1 m,而在z方向的定位误差相对较差。这主要是因为所有伪卫星均分布在接收机天线的上方,因此会造成在垂直方向上的定位误差变大。对于室内定位的场景,重点关注的是平面定位精度,可以考虑预先标定天线高度z,并将z值作为已知数代入定位计算,达到减小计算量的目的,并得到更高的平面定位精度。

      2) 计算出天线方位角的最大误差不超过5°,多数情况在2°以内,说明将天线方位角作为未知参数进行估计在理论上是可行的。但是由于估计方位角涉及到三角函数的泰勒级数展开,实现较为复杂。而且在进行最小二乘方法求解时,会让迭代计算次数明显增加,解的质量变差。由于主流的电子罗盘等定向设备均能够提供2°以内的角度测量精度,因此在实际定位场景中,选择电子罗盘等角度测量设备来估计天线方位角,可以显著降低求解过程的复杂性。

      3) 靠近定位区域中心的测试点的平面定位效果较好,而靠近定位区域边缘的测试点的定位结果相对较差。这也可以解释为中心区域具有较好的伪卫星几何结构,对应着较小的DOP值,而边缘区域的伪卫星几何结构变差,DOP值会明显增大。进一步观察图 2可以看出,水平方向的定位误差与天线方位角的取值有着密切关系。天线方位角为0°,意味着双天线连线平行于x轴,那么此时x方向的定位精度很高,而与双天线连线垂直的方向上的定位精度较低;天线方位角为90°,意味着双天线连线与y轴平行,在y方向上能够达到相对更高的定位精度,而在x方向上的定位精度会有明显的降低。后续的真实伪卫星场景下的实测结果也有类似现象。这种现象可以解释为:当入射的伪卫星信号垂直于双天线的连线,两个天线捕获到的信号的载波相位测量值的分辨率较低;而当入射的伪卫星信号平行于双天线连线方向时,两个天线捕获到伪卫星信号的载波相位测量值有较大的差别,带来了更高的载波相位测量值的分辨率。

    • 在实测试验中,采用的5颗商用伪卫星的型号均为GSG-L1,它可以输出包含简化的导航电文和C/A码的L1频段载波。伪卫星的发射天线选择了线性极化的全向天线,可以支持较大的覆盖范围。USRP是专用的软件无线电的硬件平台,通过搭配合适的子板,可以捕获全系列的导航卫星信号,在本试验中被选择作为双天线软件接收机的射频前端。两块型号为DBSRX的子板被装配在USRP中,可以捕获800 MHz至2.4 GHz的信号,能够完整地捕获到工作在L1频段的伪卫星信号。两块DBSRX的子板共用同一个晶振,保证了两个天线输入数据的完全同步。两个接收机天线选择了常规的贴片天线。USRP将处理后的两路数字中频信号通过USB接口传递到计算机。计算机上的双通道软件接收机程序能够实时地对两路输入信号进行处理,并输出RINEX格式的数据文件,载波相位测量值的输出频率为10 Hz。

      试验在一间面积为10 m×7 m的房间内进行,5颗伪卫星的天线被布设在了屋顶,如图 3所示。试验采用了独立的本地坐标系,5颗伪卫星天线的坐标预先采用全站仪进行了精确测定,三维坐标见表 1

      图  3  伪卫星天线布置

      Figure 3.  Distribution of Pseudolite Antennas

      双天线被依次放置在4个测试点,测试点与双天线几何中心严格对准,双天线间的距离固定为0.095 m,天线高度为0.2 m,天线方位角由型号为DCM260的电子罗盘测量得到,其方位角的标称测量精度为0.8°。在每个测试点处,双天线软件接收机均稳定测量60个历元的数据,然后按照§2的载波相位差值定位原理计算出双天线几何中心的平面坐标。

      定位结果如图 4所示,伪卫星的位置用黑方块表示,4个测试点的真实位置用黑色的点来表示。而在4个测试点逐个历元的定位结果分别用红、蓝、绿和青等颜色的点来表示。

      图  4  4个测试点的定位结果

      Figure 4.  Positioning Results for Four Test Points

      4个测试点的x方向和y方向的定位结果统计见表 3。从表 3可知,平面方向的平均定位误差基本可以保持在0.35 m以内,而最大的定位误差可以控制在1 m左右。这与仿真得到的结果基本一致,证明了利用载波相位差值的伪卫星定位方法的可行性。

      表 3  4个测试点的定位结果统计

      Table 3.  Statistical Positioning Results for Four Test Points

      测试点
      颜色
      平均误差 最大误差
      x/m y/m x/m y/m
      0.15 0.29 0.34 1.16
      0.21 0.35 0.47 1.48
      绿 0.11 0.20 0.33 1.02
      0.23 0.22 1.33 0.39

      图 4中的红色、蓝色、绿色测试点,双天线的连线方向平行于x轴,相应地在x方向能够达到更高的定位精度;而在青色测试点,双天线的连线方向平行于y轴,在y方向上的定位精度更好,这也与之前的仿真结果保持一致。

    • 本文介绍了一种利用载波相位差值的伪卫星定位方法,能够在无需基准站的协助,无需伪卫星之间的时钟同步,也无需求解整周模糊的情况下实现亚米级的定位精度。包含5颗伪卫星的仿真结果和基于商用伪卫星及双通道软件接收机的实测数据均验证了该定位方法的有效性,能够作为现有的室内定位技术的有效补充。由于室内环境中多路径现象严重,在实测的过程中,当双天线被放置在靠近墙壁和玻璃的位置时,双天线软件接收机无法持续捕获到全部伪卫星信号,导致无法完成定位计算,因此在下一阶段,会针对多路径对定位带来的影响和抑制方法进行更深入的研究。

参考文献 (10)

目录

    /

    返回文章
    返回