相比较于室外广域无线定位，室内定位有自己的特点，其范围一般较小，范围通常在几百米以内；另一方面由于室内环境的复杂，干扰因素较多，容易对定位信号产生影响。而卫星导航定位在室内由于受到建筑物、家具等的遮挡，卫星导航信号衰减很快，而且即使有一些信号，由于环境因素复杂，导致多径现象严重，使其无法在室内定位应用。
自20世纪90年代末起，室内无线定位技术开始受到人们的广泛关注。经过近20年的研究，室内定位技术已经取得了一些成果[2]，有一些比较成熟的技术得到了初步应用。室内定位技术分为局部室内定位和大范围室内定位，目前最常见的局域定位技术是WiFi室内定位，其基于信号能量的指纹匹配定位技术，定位精度能够达到3 m，但是信号受环境影响较大，而且需要前期采集大量的数据，定位时需要与数据库的数据进行比对。另一种室内定位技术UWB(Ultra Wideband)能够达到很高的定位精度，典型值为0.3 m，但缺点也很明显，即覆盖范围小。除此之外，还有一些新的局域室内定位技术，比如Zig Bee、蓝牙、射频标签等[3-4]，这些局域室内定位技术几乎都不是专门用作定位应用的，主要作用还是通信应用。
目前能够实现大范围室内定位的技术主要是架设地面伪卫星系统和利用通信基站的室内定位系统。北京邮电大学的时分&码分正交频分复用(TC-OFDM)[5]广域室内定位技术利用现有的通信基站，通过将定位信号承载到手机通信的信号上，能够实现大范围室内定位，系统精度虽然较高，但是大范围室内定位的技术复杂，而且需要较长的布设周期和较大的成本。地面伪卫星系统实现室内定位的代表是澳大利亚的Locata系统[6]，通过增大地面伪卫星发射信号的功率，在室内也可达到较高的定位精度，但其覆盖范围有限或需架设很多伪卫星基站。
目前，针对大型、复杂、多功能建筑，其内部结构复杂、信号环境恶劣、导航信号衰减幅度大，而且多径影响严重，无法直接利用卫星导航信号进行定位；并且在大型火灾、暴力恐怖袭击等灾害现场或事件发生地，电力供应常常被切断或无法使用，停电范围可能波及多个街区，导致目前已有的室内定位手段无法展开应用。针对突发事件现场，对定位系统的基本要求是独立、快速、易用及高精度，目前还缺乏成熟的技术或系统可以应用[7]。
当前室内定位方案难以成熟应用的问题，主要是覆盖范围与定位精度不可兼得[8]。如UWB的精度较高，但穿透能力有限；Locata和TC-OFDM系统为了保证覆盖，需要布设较多基站[9]；惯导模块或设备没有覆盖范围的限制，但存在积累误差[10]。
针对以上问题，本文提出了一种基于地面基站的大区域室内定位方案，不同于发射卫星信号的伪卫星，地面基站发送的是专门设计的导航信号，在频率、信号体制、定位方式等方面进行改进。考虑到频率与信号穿越性能的关系，选择了较低频段作为信号载波，以及采用码伪距和载波伪距联合定位的方式。本文利用所提出的新型定位方式，搭建了一个基本测试系统，通过信号的产生、发射，以及空间无线传播，进行信号的捕获、跟踪和伪距求差解算，初步验证了本文方法的可行性及有效性。
1 系统模型 系统的设计采用覆盖与精度统筹考虑的策略，选用了76 MHz的载波频率，具有良好的信号穿透能力以及障碍翻越能力。采用伪随机码进行初步定位，获取大致的用户位置。采用载波进行精细测距，最终实现高精度的定位。
根据以上设计思路，专门研制了实际系统。系统包括多个发射基站及接收终端，下面分别进行介绍。发射端的组成框图如图 1所示。接收端的组成框图如图 2所示。
 |
| 图 1 发射端组成 Fig. 1 Transmitter structure |
| 图选项 |
 |
| 图 2 接收端组成 Fig. 2 Receiver structure |
| 图选项 |
图中各部分的功能介绍如下：
现场可编程门阵列(FPGA)：完成数字基带信号的生成与调制，生成的数字信号包括导航电文、伪随机码和信号载波，采用低数字中频方式，将导航电文、伪随机码对数字载波进行调制。
AD9361射频：完成数字信号向射频信号的转换，数字基带信号输入给AD9361模块后，将进行滤波、数模转换、上变频等一系列处理，产生所需射频信号。对AD9361射频本振频率进行设置，产生中心频率为76 MHz的射频信号。
47 dB功放：将AD9361输出的模拟信号功率放大到1 W。
自动增益+固定增益功放：将天线接收后的信号进行功率放大，使信号的适应性增加，能够正常进行采样。
中频采样：将天线接收到的模拟信号使用62 MHz进行采样。采样后数据使用软件接收机进行捕获、跟踪，输出码相位和载波相位观测值。
2 定位方案 系统研制完成后，通过在地面合适位置布设4~6个地面基站，使得在所服务的区域内都可以收到至少4个基站的定位信号，接收机对接收到的信号进行捕获跟踪，得到到达时间(TDOA)即距离的测量值，并进行定位解算，就可以获得用户的实际位置。定位过程如下：
1) 按布局设计搭建多个基站，使得服务区域内的GDOP值尽量小。
2) 通过地理测绘或全球导航卫星系统(GNSS)差分技术等，对基站的真实位置进行标定。
3) 通过光纤技术、时间双向比对技术等，实现对基站之间的高精度时间同步。
4) 各基站同步发射导航信号，利用低频载波增加覆盖范围。
5) 接收机接收4个质量较好、GDOP值较小的基站信号进行距离测量并实时定位，并且通过载波相位测量提高定位精度。
本系统的定位算法原理如下：利用载波波长较小，载波跟踪精度高和受多径影响小的特点，使用载波来进行高精度定位，缺点是载波整周是位置的。本定位方案使用码伪距来确定未知的载波整周，将载波伪距与码伪距确定的载波整周相结合，得到较精确的目标到信号源距离差，使用下面的TDOA方法完成定位结算。
基于TDOA的无线电定位是一种重要的目标三维定位方法，被广泛应用于导航、雷达等定位系统，定位方法主要分为2种：一是迭代算法[11]，二是解析算法[12-13]。迭代算法的优越性体现在能充分利用系统的冗余信息，计算精度高，但计算复杂，实时性差。解析算法精度较低，但实时性好。考虑到室内运动速度不高，本文采用迭代算法。
当已知信号源的位置坐标为Si(xi, yi, zi)T，目标位置坐标为Tr(xr, yr, zr)T，得到目标到信号源的距离差为
 | (1) |
对目标到信号源的距离取平方：
 | (2) |
通过对目标接收机码相位及载波相位观测值分别两两做差，可以得到距离差dij，在实际测试中，距离差含有各种因素引起的误差εij，假设距离差的理论值为dijt，所以距离差测量值可以表示为
 | (3) |
将式(1)移项平方，得到ri2=(dij+rj)2，将式(2)代入，可得
 | (4) |
式中：
由于测试误差的存在，式(4)的等号左边并不等于0，可用微小量δij代替。实际情况下，可以将式(4)表示成矩阵形式，通过最小二乘法就可以求解目标位置。
这种实现方案，优点是可以避免卫星导航中的主要误差因素，如星历、星钟、对流层、电离层等误差，导航信号的完好性也大大增强，并且利用低频载波可以增加覆盖范围和穿透能力，缺点是在地球表面的多径信号较多，需要通过抗多径技术减弱多径的影响。
与其他室内定位技术相比[14]，比如WiFi定位，本文方法利用伪随机码的相关性进行测距和定位，精度明显高于利用WiFi信号的功率测距方法，所以WiFi定位目前大多通过指纹匹配技术来提高精度；而UWB定位，其测距精度最高，但UWB工作频率也很高，受周围环境的影响很大，覆盖范围十分受限；Locata系统是一种地基伪卫星定位系统，目前有一定应用，但其在每个室内都需要布设多个伪卫星基站并实现同步，成本及布局问题限制了其进一步应用[15]。
3 测试结果处理与分析 3.1 在室外空旷条件下测试 首先测试信号源间码相位和载波相位同步情况，接着测试76 MHz发射频率时中垂线处和偏移处接收码相位和载波相位能否体现距离差的变化。
以操场为测试环境(见图 3)，在主席台右侧看台最高处台阶放置两台信号发射源，发射功率为30 dBm，信号源1、2号星间距为45 m。在对侧跑道上取3个等间隔的测试点(接收机位置)，分别为1号点、2号点、3号点。信号源与测试点构成长方形。2号点(取1号星、2号星中垂线与长方形的交点)距两信号源中点87 m。1号点、3号点与两信号源的距离差约为11 m。信号源与接收机的位置如图 3所示。
 |
| 图 3 操场测试环境 Fig. 3 Playground test environment |
| 图选项 |
1) 首先进行2个信号源的同步测试，启动1号星与2号星，经过同步模块进行时钟校准后，使用电脑同时配置信号源，使两台信号源分别产生给定的定位信号。
将两台信号源通过导线同时连接接收机，通过接收机显示界面观察码相位，并采集码相位及载波相位观测量，其码相位差和载波相位差的测试曲线分别如图 4和图 5所示(测试2次)。其对应的测试数据如表 1所示。
 |
| 图 4 码相位差值(同步测试) Fig. 4 Code phase difference value (synchronous test) |
| 图选项 |
 |
| 图 5 载波相位差值(同步测试) Fig. 5 Carrier phase difference value (synchronous test) |
| 图选项 |
表 1 测试数据(同步测试) Table 1 Test data (synchronous test)
| 观测量 | 测试分组 | 均值/m | 标准差/cm |
| 码相位差 | 第1次 | 1.63 | 57.7 |
| 第2次 | 1.68 | 53.7 | |
| 载波相位差 | 第1次 | 0.991 5 | 0.18 |
| 第2次 | 2.943 1 | 0.11 |
表选项
测试结果表明，系统工作正常，但不同组测试时可能会有半周差，可在后期处理时校正。可以看出，在导线连接条件下测试的精度较高，原因是由于没有引入空间传播方面的误差。
2) 下面进行运动状态下的测试，分别连接好发射、同步、接收设备，并启动其正常工作，在操场携带接收端设备首先从1号点沿直线走到2号点，再走到3号点，在该过程中在3个测试点各停留约200 s左右，实时采集并处理2个信号源发出的定位信号，其码相位差和载波相位差的测试曲线分别如图 6和图 7所示。其对应的测试数据如表 2所示。
 |
| 图 6 码相位差值(运动状态下测试，1号点~3号点) Fig. 6 Code phase difference value (test under motion, point 1 to 3) |
| 图选项 |
 |
| 图 7 载波相位差值(运动状态下测试，1号点~3号点) Fig. 7 Carrier phase difference value (test under motion, point 1 to 3) |
| 图选项 |
表 2 测试数据(运动状态下测试，1号点~3号点) Table 2 Test data (test under motion, point 1 to 3)
| m | |||||
| 测试点 | 码相位差 | 载波相位差 | |||
| 均值 | 标准差 | 均值 | 标准差 | ||
| 1 | 3.66 | 66.8 | 13.2 | 8.397 | |
| 2 | 2.83 | 5.6 | 31.1 | 0.19 | |
| 3 | -16.6 | 5.4 | 39.5 | 0.13 | |
表选项
测试结果表明：码相位差值在1号点有较大的波动，在2号点、3号点比较平稳，标准差符合预期。同样，载波相位差在1号点时出现了几次跳跃，与静止状态不符，但2号点、3号点基本平稳，标准差也小，符合预期，并且比较理想地反映了从2号点到3号点的运动过程。
3) 最后再从3号点沿直线返回2号点，再返回1号点，同样在3个测试点各停留约200 s左右。实时采集并处理2个信号源发出的定位信号，其码相位差和载波相位差的测试曲线分别如图 8和图 9所示。其对应的测试数据如表 3所示。
 |
| 图 8 码相位差值(运动状态下测试，3号点~1号点) Fig. 8 Code phase difference value (test under motion, point 3 to 1) |
| 图选项 |
 |
| 图 9 载波相位差值(运动状态下测试，3号点~1号点) Fig. 9 Carrier phase difference value (test under motion, point 3 to 1) |
| 图选项 |
表 3 测试数据(运动状态下测试，3号点~1号点) Table 3 Test data (test under motion, point 3 to 1)
| m | |||||
| 测试点 | 码相位差 | 载波相位差 | |||
| 均值 | 标准差 | 均值 | 标准差 | ||
| 3 | -26.5 | 7.1 | -1.38 | 0.138 | |
| 2 | 4.4 | 4 | -9.38 | 0.125 | |
| 1 | 24.5 | 53.8 | -8.67 | 6.4 | |
表选项
测试结果表明：在2号点、3号点及其连接路径上的定位精度比较好，同样在1号点出现了大的误差和波动，从3号点走到2号点的路径刻画十分理想，说明了测试数据的可信性。虽然行走的方向与1号点~3号点相反，但2个方向的测试结果是一致的，2号点和3号点载波相位差对应的两次位置差值分别为8.4、8 m，比较接近，说明了测试的可重复性及有效性。
3.2 在复杂室内条件下测试 在前期室外测试验证了系统工作原理及测试结果正确性的基础上，开展了实际复杂室内条件下的测试，测试场景为北京航空航天大学新主楼，如图 10所示。
 |
| 图 10 新主楼室内测试环境 Fig. 10 New main building indoor test environment |
| 图选项 |
图中发射信号源放在位于大楼中间露天的地面绿点处，发射功率4 W，测试点为2楼室内的红点处，可以看到环境比较复杂，尤其外圈的一些红点，信号要穿过多重建筑实体才能够到达。
通过在红点处逐个测试，采集信号源穿过楼体的到达信号，通过事后的相应处理，得到以下结果：
1) 在所有测试点(红点处)都可以接收、捕获到信号源发出的导航信号，说明本文的低频策略是有效的，该方法的覆盖能力较强。如在测试点2对3个基站信号的捕获情况如图 11所示。
 |
| 图 11 测试点2的捕获相关峰 Fig. 11 Capture correlation peak at test Point 2 |
| 图选项 |
2) 由于环境复杂、多径严重等原因，测试数据的跟踪结果表明码相位和载波相位的误差都比较大，尤其信号路径长的测试点还出现了大的偏差，严重的已无法用于定位解算，需要进一步进行抗多径等处理。同样在测试点2，对3个基站的码跟踪曲线及基站间的码伪距差曲线如图 12所示。图 12(d)为PRN1-2码相位示意图，最大值和最小值分别为83.071 m和72.570 m，标准差为1.847 9 m，均值为78.605 4 m; 图 12(e)为PRN1-3码相位示意图，最大值和最小值分别为-21.904 7 m和-121.623 9 m，标准差为34.529 4 m，均值为-63.764 4 m; 图 12(f)为PRN1-3码相位示意图，最大值和最小值分别为57.223 8 m和-44.515 9 m，标准差为34.770 9 m，均值为14.840 9 m。从图中可以看出，各基站间的码伪距差波动较小，但其均值与事先标定的真值偏差较大，分析是由多径因素造成的固定偏差。
 |
| 图 12 测试点2跟踪曲线及基站间距离差曲线 Fig. 12 Tracking curve and distance difference curves between base stations at test point 2 |
| 图选项 |
4 结论 本文经过不同环境的测试并对结果进行分析，初步验证了所提方法的可行性，证明了系统定位的可行性和良好的覆盖能力，但也发现了系统进一步应用需解决的问题。综合分析本文测试结果，可以得出如下结论：
1) 在室外开阔环境(操场2号点、3号点)的测试，位置差的精度较高，动态测试符合预期，说明了方法的可行性和有效性。
2) 如果接收天线被行人遮挡(操场1号点，有跑步人员)，会出现码相位误差变大和载波相位跳变的情况，在使用中应该尽量避免。
3) 本文信号的穿透能力较强，对复杂室内的覆盖范围较大，但需要进一步解决多径严重，以及信号衰减大、信号质量下降等问题，才可以实现大区域的室内定位，这也是下一步的工作重点。
参考文献
| [1] | 刘蔚, 康永. 室内定位技术应用分析研究[J]. 现代导航, 2016, 7(2): 86-93. LIU W, KANG Y. Application analysis of indoor positioning technology[J]. Modern Navigation, 2016, 7(2): 86-93. (in Chinese) |
| [2] | 阮陵, 张翎, 许越, 等. 室内定位:分类、方法与应用综述[J]. 地理信息世界, 2015, 22(2): 8-14. RUAN L, ZHANG L, XU Y, et al. Indoor positioning:A review of classification, methods and applications[J]. Geomatics World, 2015, 22(2): 8-14. DOI:10.3969/j.issn.1672-1586.2015.02.002 (in Chinese) |
| [3] | LIU H, DARABI H, BANERJEE P, et al. Survey of wireless indoor positioning techniques and systems[J]. IEEE Transactions on Systems Man & Cybernetics Part C, 2007, 37(6): 1067-1080. |
| [4] | MAINETTI L, PATRONO L, SERGI I.A survey on indoor positioning systems[C]//International Conference on Software, Telecommunications and Computer Networks.Piscataway, NJ: IEEE Press, 2015: 111-120. |
| [5] | 邓中亮, 尹露. 基于TC-OFDM体制的室内定位系统[J]. 电信网技术, 2015(3): 32-35. DENG Z L, YIN L. Indoor positioning system based on TC-OFDM system[J]. Information and Communications Technology and Policy, 2015(3): 32-35. (in Chinese) |
| [6] | BARNES J, RIZOS C, WANG J, et al. High precision indoor and outdoor positioning using LocataNet[J]. Positioning, 2003, 2(2): 73-82. |
| [7] | 赵锐, 钟榜, 朱祖礼. 室内定位技术及应用综述[J]. 电子科技, 2014, 27(3): 154-157. ZHAO R, ZHONG B, ZHU Z L. Overview of indoor positioning technology and application[J]. Electronic Science and Technology, 2014, 27(3): 154-157. DOI:10.3969/j.issn.1007-7820.2014.03.046 (in Chinese) |
| [8] | 许华.连续相位调制技术研究及其实现[D].西安: 西安电子科技大学, 2009. XU H.Research and implementation of continuous phase modulation technology[D].Xi'an: Xidian University of Electronic Technology, 2009(in Chinese). http://cdmd.cnki.com.cn/Article/CDMD-10701-2009064701.htm |
| [9] | 冯思泉.连续相位调制(CPM)技术研究[D].昆明: 昆明理工大学, 2007. FENG S Q.Research on continuous phase modulation (CPM) technology[D].Kunming: Kunming University of Science and Technology, 2007(in Chinese). http://cdmd.cnki.com.cn/Article/CDMD-10674-2007073076.htm |
| [10] | 黄嵩.CPFSK调制解调的研究与实现[D].绵阳: 西南科技大学, 2011. HUANG S.Research and implementation of CPFSK modulation and demodulation[D].Mianyang: Southwest University of Science and Technology, 2011(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y1998080 |
| [11] | FOY W H. Position-location solutions by Taylor-series estimation[J]. IEEE Transactions on Aerospace and Electronic Systems, 1976, 12(2): 187-194. |
| [12] | 王成, 李少洪, 王鑫全, 等. 测时差被动定位算法的研究[J]. 系统工程与电子技术, 2001, 23(11): 9-12. WANG C, LI S H, WANG X Q, et al. Research on time difference differential positioning algorithm[J]. Systems Engineering and Electronics, 2001, 23(11): 9-12. (in Chinese) |
| [13] | HUANG Y T, BENESTY J, ELKO G W, et al. Real-time passive source localization:A practical linear-correction least-squares approach[J]. IEEE Transoctions on Speech and Audio Processing, 2001, 9(8): 943-956. DOI:10.1109/89.966097 |
| [14] | 齐恒, 彭林宁, 姜禹, 等. 基于临近听域超声波TDOA室内定位的实现[J]. 数据采集与处理, 2018, 33(6): 1094-1100. QI H, PENG L N, JIANG Y, et al. Realization of indoor positioning based on near-listening ultrasonic TDOA[J]. Journal of Data Acquisition and Processing, 2018, 33(6): 1094-1100. (in Chinese) |
| [15] | 张令文, 杨刚. 超宽带室内定位关键技术[J]. 数据采集与处理, 2013, 28(6): 706-713. ZHANG L W, YANG G. Ultra-wideband indoor positioning key technology[J]. Journal of Data Acquisition and Processing, 2013, 28(6): 706-713. DOI:10.3969/j.issn.1004-9037.2013.06.002 (in Chinese) |
