
Citation: | HUANG He, QIU Kaiyue, LI Wei, LUO Dean. Indoor Geomagnetic Positioning Based on Joint Algorithm of Particle Filter[J]. Journal of Southwest Jiaotong University, 2019, 54(3): 604-610. doi: 10.3969/j.issn.0258-2724.20170927 |
精密单点定位(precise point positioning, PPP)技术自问世以来得到了极大的关注[1-4]. 得益于国际全球导航卫星服务组织(International GNSS Service, IGS)产品精度的提升,PPP静态定位精度与差分GPS相当,动态定位精度可达厘米级. PPP定位基准通常基于国际地球参考框架(international terrestrial reference frame, ITRF),这使得其在地灾监测、智能交通、精准农业、地壳运动和海啸预警等领域具有独特优势[5].
基于PPP技术的优势,国内外先后开发了一系列相关软件,比较著名的有:喷气推进实验室(Jet Propulsion Laboratory, JPL)的GIPSY OASIS (GNSS- inferred positioning system and orbit analysis simulation software)、新不伦瑞克大学(UNB)的GAPS (GNSS analysis and positioning software)、GMV (Grupo Mecánica del Vuelo Sistemas)的magicPPP (magic precise point positioning)、加拿大自然资源部(Natural Resource Canada, NRCan)的CSRS-PPP (canadian spatial reference system precise point positioning)、加泰罗尼亚理工大学(UPC)的gLab (GNSS-lab tool)、东京海洋大学的RTKLIB(real time kinematic library)和武汉大学的PRIDE (positioning racers to image & decipher the earth)等,其中GIPSY、GAPS、CSRS-PPP和magicPPP免费提供在线数据处理服务[6-9]. 西南交通大学也一直致力于GPS数据处理软件的研发,其中,PLAOD (precise location and orbit determination)软件主要目的是精密定位、定轨和大气反演. 相比较于常规软件,PLAOD功能较为丰富,除提供精密静态坐标外,还可提供高时间分辨率的精密动态坐标、对流层天顶延迟和电离层天顶延迟,其中,在进行大气估计时可固定测站坐标以加快其收敛速度. 此外,在提供低轨卫星天线相位中心的前提条件下,PLAOD可以对其进行精密定轨. PLAOD采用C#语言编写,所有代码均是自主研发,在联网的情况下可根据观测文件自动下载对应的精密轨道(sp3)、精密钟差(clk)和码偏差(DCB)等相关文件. 本文将主要介绍PLAOD的基本理论和性能分析.
对于频率为j的卫星s至接收机r间的伪距
{P(s)r,j=ρ(s)r+ctr−ct(s)+T(s)r+μjI(s)r,1+dr,j−d(s)j+ε(s)r,j,ϕ(s)r,j=ρ(s)r+ctr−ct(s)+T(s)r−μjI(s)r,1+λjZ(s)r,j+δr,j−δ(s)j+ξ(s)r,j, |
(1) |
式中:
利用式(1)构造消电离层(ionosphere-free, IF)观测方程,如式(2)所示.
{P(s)r,IF=ρ(s)r+ctr−ct(s)+T(s)r+dr,IF−d(s)IF,ϕ(s)r,IF=ρ(s)r+ctr−ct(s)+T(s)r+λIFZ(s)r,IF+δr,IF−δ(s)IF, |
(2) |
式中:下标IF表示频率j为消电离层组合的频率.
由于相位观测量中存在模糊度,式(2)中的卫星钟差基准实际上来源于伪距,即IGS等机构提供的卫星钟差为
T(s)r=Tr,Dm(s)r,D+Tr,Wm(s)r,W+gr,ncosa(s)+gr,esina(s)sine(s)tane(s)+c0, |
(3) |
式中:
要想获得真实可降水量,必须利用实际气象数据计算得到真实的天顶干延迟
利用式(1)的相位观测方程还可构造无几何距离(geometry free, GF)观测方程,如式(4)所示.
ϕ(s)r,GF=(1−μj)I(s)r,1+λ1Z(s)r,1−λjZ(s)r,j+δr,GF−δ(s)GF, |
(4) |
式中:j≠1;下标GF表示频率j为无几何距离组合频率;等式右边主要有电离层和模糊度参数,电离层参数可以进一步分解,如式(5)所示.
I(s)r,1=m(s)ION(Ir,0+gr,B(B(s)p−Br)+gr,L(L(s)p−Lr)), |
(5) |
式中:
PLAOD主要数据处理流程如图1所示.
周跳探测和粗差剔除对于最终结果至关重要,其在PLAOD中采用逐历元处理模式,以便于后续拓展实时数据处理应用;不同于传统的周跳探测算法,这里采用MW (melbourne wubbena)组合或GF组合乘以
PLAOD的数据处理核心在于扩展卡尔曼滤波[14]. 卡尔曼滤波不仅便于实时估计参数,同时可以提高数据处理效率,以单个测站间隔30 s采样24.0 h观测数据为例,PLAOD在CPU为I7、内存为16 GB的计算机上仅需10 s即可对其执行完毕. PLAOD在进行滤波时相应的参数处理策略见表2.
参数 | 处理方案 |
测站坐标 | 静态时过程噪声为 0,动态时先验值源自伪距单点定位 |
模糊度 | 过程噪声为 0 |
对流层 | 随机游走,过程噪声默认为 4 cm2/h,动态时顾及历元间高差引起的影响 |
对流层梯度 | 随机游走 |
电离层 | 随机游走 |
电离层梯度 | 随机游走 |
接收机钟差 | 先验值源自伪距单点定位 |
卫星钟差 | 20 min弧段二阶多项式拟合 |
卫星坐标 | 10 阶滑动拉格朗日插值 |
值得一提的是表2中关于卫星钟差的处理,早期的IGS钟差产品输出间隔通常为5 min,对于高采样率的观测文件直接采用线性内插方法来获取钟差可能会导致其精度损失较大,PLAOD中利用20 min弧段内的原始钟差数据构造二阶函数则可有效减少钟差拟合精度的损失.
常规软件中常用高度角的正弦来定义观测量的随机模型(测量噪声),PLAOD中的测量噪声矩阵依赖于对流层投影函数或者电离层投影函数,以便于低高度角观测量得到更合理的利用. 以定位测量噪声R为例,其对应的矩阵如式(6)所示.
R=[(0.300m(1)r,D)20⋯000(0.003m(1)r,D)2⋯0000⋱0000⋯(0.300m(s)r,D)2000⋯0(0.003m(s)r,D)2], |
(6) |
式中:0.300表示该卫星伪距观测量在天顶方向的测量噪声为0.300 m; 0.003表示该卫星的相位观测量噪声为对应伪距的1/100.
选择2020年中国区域1月1日—3日期间能下载的IGS站点单天观测文件作为测试数据,其采间隔为30 s,站点分布如表3所示.
站点 | 经度/(°) | 纬度/(°) |
CHAN(长春) | 125.443 | 43.790 |
BJFS(北京房山) | 115.892 | 39.609 |
BJNM(北京) | 116.224 | 40.245 |
HKSL(香港小冷水) | 113.928 | 22.372 |
HKWS(香港黄石) | 114.335 | 22.434 |
JFNG(九峰) | 114.491 | 30.516 |
LHAZ(拉萨) | 91.104 | 29.657 |
TWTF(桃园) | 121.164 | 24.954 |
基于IGS最终精密星历和卫星钟差,利用PLAOD解算上述观测文件后,所得结果是基于ITRF的时间序列. 通常解算坐标序列最终历元的结果能在一定程度上反映软件解算的准确性和稳定性,将最终历元坐标与IGS提供的SINEX (solution independent exchange)文件坐标作差,并将其转换为站心坐标系中北、东和高方向的偏差(dN、dE和dU),其结果如表4所示. 此外,坐标收敛情况作为软件完好性统计的指标之一,表4罗列了每个测站收敛到10 cm准确度所需的观测时长,需要注意的是这里的时长并不一定代表首次收敛到10 cm所需观测时长,而强调的是后续结果不再出现超过10 cm的偏差,同时10 cm表示三维距离总偏差.
站点 | 年积日 | dN /mm | dE /mm | dU /mm | 收敛时间 /min |
CHAN | 001 | 1.9 | −9.2 | −11.3 | 32.0 |
002 | 2.4 | −5.7 | −7.5 | 19.0 | |
003 | 2.0 | −7.1 | −9.2 | 43.0 | |
BJFS | 001 | 2.0 | 1.2 | −1.8 | 30.0 |
002 | 2.5 | 0.7 | −0.5 | 21.0 | |
003 | 2.5 | −0.2 | 2.8 | 23.0 | |
BJNM | 001 | 2.6 | −3.6 | −0.8 | 29.0 |
002 | 4.5 | 0.4 | 1.4 | 50.0 | |
003 | 4.0 | −4.5 | 4.3 | 46.5 | |
HKSL | 001 | 2.0 | 0.5 | −3.9 | 30.5 |
002 | 3.8 | 0.1 | 1.7 | 25.0 | |
003 | 4.4 | −2.5 | 2.0 | 20.5 | |
HKWS | 001 | 1.0 | 1.5 | 2.2 | 30.0 |
002 | 3.3 | 1.9 | −0.1 | 52.5 | |
003 | 3.5 | −0.3 | 1.9 | 104.5 | |
JFNG | 001 | 2.0 | 3.2 | 0.6 | 20.0 |
002 | 2.8 | 4.4 | 9.0 | 13.0 | |
003 | 3.1 | 3.4 | 8.1 | 19.0 | |
LHAZ | 001 | 3.6 | −1.5 | −3.8 | 42.0 |
002 | 4.1 | 1.5 | −0.3 | 57.0 | |
003 | 3.5 | −1.5 | −4.4 | 44.0 | |
TWTF | 001 | 0.8 | 3.9 | −9.3 | 26.5 |
002 | 4.4 | 1.0 | −5.4 | 21.5 | |
003 | 4.0 | 0.1 | −6.3 | 19.0 |
通过表4可以得出:各测站每天最终历元的坐标单方向偏差均在毫米级,除了CHAN站1月1日的高程方向偏差,这表明PLAOD与IGS分析中心的误差模型吻合度较高. 通过对各测站的收敛情况进行检查发现它们的收敛时间均在1.0 h以内,除了HKWS站1月3日的结果,其中,JFNG站三天均保持在20.0 min以内,最短时为13.0 min;而HKWS站1月3日的收敛时间长达104.5 min,经进一步分析发现该站伪距噪声较大.
为进一步展示PLAOD解算测站坐标偏差序列整体情况,选取上述偏差最大的结果和收敛时间最长的结果,分别对应CHAN站 1月1日和HKWS站 1月3日偏差序列,如图2和图3所示,图中:dX、dY和dZ分别为X、Y和Z方向的偏差,历元间隔为采样间隔(如无特殊说明下同).
图2显示:CHAN站1月1日偏差序列在收敛后波动较小,其收敛后的前期结果优于后期结果,这与卫星的几何形状分布等因素有一定的关联,但从整体上看经近3.5 h收敛后的偏差持续保持在毫米级. 图3显示:HKWS站1月3日偏差序列收敛时间较长,但在近14.0 h收敛后波动一直保持较小,不超过5 mm. 上述结果均证实PLAOD软件具备较好的稳健性.
相对于静态定位,动态定位的实践应用更广泛,但在动态定位过程中,待估位置参数个数有所增加. 为检验PLAOD的动态定位性能,首先,将上述CHAN站1月1日的静态数据模拟成动态模式进行解算,所得结果如图4所示.
如图4所示,CHAN站1月1日动态定位偏差序列在收敛后,单方向均在cm级,将偏差序列转换成站心坐标系后显示高程方向误差占主要部分. 尽管每历元均需解算位置参数,但图4显示动态定位收敛速度也较快,仅比静态定位时慢几分钟.
另外利用2011年3月11日日本台站MIZU的观测数据来进行测试,MIZU站距离宫城震中较近,并且由于地震原因导致该站观测数据在地震后不久便发生中断,其观测数据解算结果如图5所示.
30 s分辨率的MIZU站坐标偏差序列显示,X方向最大震时位移将近2.5 m,Y方向最大震时位移超过1.5 m,Z方向最大震时位移超过1.0 m,这与相关文献的结果高度吻合[15]. 此外,尽管数据观测发生中断,从图4仍可显示:X、Y和Z方向的位移在一定时间后趋于稳定,这表明PLAOD不仅可以反演弹性形变,也可以探测永久形变.
虽然地震时期,站点形变的速率不是很高,而生产实践中移动载体的速率千变万化,但PLAOD不受该因素的制约,可以应用到各种精密动态定位场景中.
作为PPP的副产品,GPS反演的ZTD是气象分析和预报的重要参考指标. 相比于其他软件,PLAOD的ZTD过程噪声设置值较大,这主要是用以应对恶劣天气而设定的,以暴雨事件为例,ZTD在1.0 h内变化完全可以达到2 cm. 以2020年8月12日BJFS站反演的ZTD为例,现给出PLAOD、CSRS-PPP和IGS分析中心JPL提供的静态ZTD解算结果:其中PLAOD和CSRS-PPP 30 s提供一次ZTD值,均只采用向前滤波法;JPL则是5 min提供一次ZTD值,而且采用了向后滤波法,即类似批处理.
图6显示:3条曲线吻合度较高,尤其是CSRS-PPP与PLAOD对应的时序变化趋势. 由于JPL的结果采用了向后滤波,无法判断其ZTD收敛速度,不过这可以为PLAOD和CSRS-PPP的ZTD收敛速度提供参照. 易见CSRS-PPP的ZTD收敛速度较PLAOD快,这与不同软件采用的先验ZTD模型有关:CSRS-PPP采用VMF (vienna mapping function)投影函数模型,该模型的优点是精度高,缺点是需要实时更新相应文件,即每6.0 h更新一次[16];PLAOD采用的是GPT (global pressure and temperature)和GMF (global mapping function)模型,GPT模型准确度较VMF略低,但该模型不需要额外提供相应文件[17-18]. 此外,在各项误差模型均准确的前提下,投影函数是影响ZTD解算准确度的主要因素,上述结果验证了GMF与VMF的高度吻合性. 大量数据验算表明PLAOD所得静态ZTD转化成可降水量对应的精度通常在1~2 mm.
PLAOD在估计ZTD时可以将测站坐标固定,以减少参数个数从而加速参数估计的收敛速度,同时也可以利用动态模式估计ZTD. 静止测站只能获取地面有限点位上空的ZTD,而移动平台能反演空间更多点位上空的ZTD,但移动平台反演的ZTD精度如何有待验证. 这里以两个观测数据来进行测试:第一个观测数据即上述TWTF站2020年1月1日对应的观测量,将该静止观测量按动态模式进行解算;第二个观测数据则来自2008年8月1日直升机载观测量,采用的接收机为拓普康,但天线类型未知,故无法对其做接收机天线相位中心改正. 利用PLAOD解算TWTF站的ZTD时,采用了常规静止模式(S模式)、动态模式(K模式)和固定测站坐标模式(F模式),参考值同样为JPL提供的ZTD,具体结果如图7所示;PLAOD在解算动态测站包括机载观测量时会根据相邻历元间的高差判断是否自动增加对流层过程噪声,噪声增加的大小依赖于GPT模型,而CSRS-PPP中动态解算时不估计ZTD,直接利用VMF模型,两软件所得机载GPS ZTD如图8所示.
从图7可以看出:PLAOD解算出的F模式ZTD和S模式ZTD与JPL的结果高度吻合,其中,F模式ZTD的收敛速度更快,而K模式ZTD的解算精度有所下降,其局部波动明显,但在整体范围内也可以达到1~2 cm的准确度,对应可降水量则为2~3 mm的准确度. 从图8可以看出:PLAOD的ZTD估计值与VMF的模型值之间具有一定的吻合性,ZTD的变化与直升机的高度有较大关联,如果加入接收机天线相位中心改正能进一步提高ZTD的估计精度. 算例结果意味着移动平台GPS接收机反演ZTD具有一定的可行性,这将会为ZTD的获取提供新的途径,尤其是弥补其空间分辨率,有利于水汽层析等研究.
前述内容中提及PLAOD具备电离层估计功能,由于电离层估计是建立在相位GF组合观测量的基础上,忽略相位硬件延迟改正,实际上所得的电离层中可能包含部分相位硬件延迟误差. 以上述TWTF站1月1日的第一波段频率天顶电离层延迟估计为例,利用IGS的IONEX (ionosphere exchange)文件值作为参考,所得结果如图9所示.
图9显示:PLAOD解算的电离层天顶延迟与IONEX对应的电离层天顶延迟具有一定的吻合度,虽然在局部会有一定的差异,但整体趋势基本相同,在当地时间PM 1:00—2:00到达峰值(中国区域内地方时比GPS时快8.0 h,地方时PM 1:00对应图中历元为600处),在凌晨降至谷值. 实际上该结果是IGS跟踪站观测量反演电离层的典型结果,实验证明相比较利用非差非组合PPP技术估计电离层,该结果与IONEX的吻合度更高,这意味着利用PLAOD监测电离层的变化具有一定的可行性.
随着社会科技和应用需求的发展,地球周围上空发射了越来越多的低轨卫星,这些卫星的应用是建立在其精密轨道确定的基础之上. PLAOD具备低轨卫星精密轨道确定功能,采用的方式为几何法定轨,其本质为高动态PPP,所得结果也被称为运动学轨道,该方法有利于地球重力场模型反演. 以GRACE(gravity recovery and climate experiment)双星A和B为例,利用PLAOD对2016年1月3日对应的间隔10 s采样观测量进行全历元解算,与JPL提供的简化动力学轨道值作差,其结果如图10和图11所示.
从图10和图11可看出:低轨卫星轨道准确度收敛到厘米级所需时间较短,通常只需几分钟,这得益于低轨卫星和GPS卫星均在高速飞行,使得GPS卫星对应低轨卫星的几何形状分布变化较快,有利于提高待估参数的收敛速度. 从整体上看,图中大多数历元的X、Y和Z方向偏差均在厘米级,但部分历元存在重收敛或者跳变,这主要与所接收的GPS卫星颗数以及卫星几何分布形状有关,这种情况通过向后滤波可一定程度改善,但无法完全消除. 通常在可接收卫星颗数为4且对应几何精度因子值又较大时,就会导致解算的位置参数存在较大的偏差,这在低轨卫星精密运动学轨道确定时不可避免,因为低轨卫星天线不能像地面接收机天线持续朝向天顶方向,这也是为什么国外众多科研机构提供的运动学轨道不连续的原因. 上述GRACE卫星轨道解算结果表明:和地表运动测站一样,只要低轨卫星能够观测足够的GPS卫星,加上严密的数据预处理方法,几何法定轨也能获取厘米级的运动学轨道结果.
1) PLAOD的长时间静态定位精度可达毫米级,可用于高精度控制点坐标解算和形变监测等领域;
2) 在观测条件尚可的情况下,PLAOD的动态定位精度可达厘米级,适用于精密轨迹追踪和动态形变反演等领域;
3) 作为PPP的副产品,PLAOD估计的ZTD可转换为精密水汽,用于气象监测等领域;
4) 同样作为PLAOD非差解算副产品的电离层延迟,与IGS提供的电离层产品具有较高的吻合度,可用于地球上空电子活动事件分析;
5) PLAOD可解算低轨卫星精密轨道,进而为地球重力场反演等提供前期数据.
为适应GNSS应用的发展,后期将会在PLAOD中加入模糊度固定和多系统多频数据处理功能,尤其是重点拓展北斗的数据处理和应用.
致谢:感谢IGS提供观测文件和精密产品;感谢NRCan提供CSRS-PPP服务作为参考;感谢GFZ提供GRACE相关数据;感谢审稿专家提出的宝贵意见;本文受到河北省高等学校科学技术研究项目(QN2020216)资助.
VALVERDE T G, SOLA A G, HAGRAS H, et al. A fuzzy logic-based system for indoor localization using WiFi in ambient intelligent environments[J]. IEEE Transactions on Fuzzy Systems, 2013, 21(4): 702-718.
|
LUOH L. ZigBee-based intelligent indoor positioning system soft computing[J]. Soft Comput., 2013, 18: 443-456.
|
NI L M, LIU Yunhao, LAU Y C, et al. LANDMARC: indoor location sensing using active RFID[C]//Proceedings of the First IEEE International Conference on Pervasive Computing and Communications, 2003, Fort Worth: [s.n.], 2003: 407-415
|
SKVORTZOV V Y, LEE H K, BANG S, et al. Application of electronic compass for mobile robot in an indoor environment[C]//IEEE International Conference on Robotics and Automation. [S.l.]: IEEE, 2007: 2963-2970
|
LI B, GALLAGHER T, DEMPSTER A G, et al, How feasible is the use of magnetic field alone for indoor positioning?[C]//2012 International Conference on Indoor Positioning and Indoor Navigation (IPIN). Sydney: [s.n.], 2012: 1-9
|
SUKSAKULCHAI S, THONGCHAI S, WILKES D M, et al. Mobile robot localization using an electronic compass for corridor environment[C/OL]//2000 IEEE International Conference on Systems, Man, and Cybernetics. Nashville: IEEE, 2000, 5: 3354-3359
|
SIIKSAKULCHAI S, THONGCHAI S, WILKES D M, et al. Mobile robot localization using an electronic compass for corridor environment[C/OL]//IEEE International Conference on Systems, 2000[2017-12-12]. http://www8.cs.umu.se/research/ifor/dl/LOCALIZATION-NAVIGATION/Mobile%20Robot%20Localization%20using%20an%20Electronic%20Compass%20for%20Corridor.pdf
|
Kim S E, Kim Y, Yoon J, et al. Indoor positioning system using geomagnetic anomalies for smartphones[C/OL]//2012 International Conference on Indoor Positioning and Indoor Navigation. Sydney: University of New South Wales, 2012[2017-10-22]. https://ieeexplore.ieee.org/abstract/document/6418947
|
FRASSL M, ANGERMANN M, LICHTENSTERN M, et al. Magnetic maps of indoor environments for precise localization of legged and non-legged locomotion[C]//IEEE International Conference on Intelligent Robots and Systems. [S.l.]: IEEE, 2013: 913-920
|
PUTTA R, MISRA M, KAPOOR D. Smartphone based indoor tracking using magnetic and indoor maps[C]//2015 IEEE Tenth International Conference on Intelligent Sensors, Sensor Networks and Information Processing. Singapore: IEEE, 2015: 7-9
|
LI B, GALLAGHER T, RIZOS C, et al. Using geomagnetic field for indoor positioning[J]. J. Appl. Geod., 2013, 7: 299-308.
|
RIEHLE T H, ANDERSON S M, LICHTER P A, et al. Indoor waypoint navigation via magnetic anomalies[C]//2011 Annual International Conference of the IEEE Engineering in Medicine and Biology Society. Boston: IEEE, 2011: 5315-5318
|
黄鹤,赵焰,王春来,等. 地磁室内定位基准图数据采集系统设计[J]. 测绘通报,2017(2): 54-59.
HUANG He, ZHAO Yan, WANG Chunlai, et al. Design of the acquisition system of indoor positioning reference map based on magnetic field data[J]. Bulletin of Surveying and Mapping, 2017(2): 54-59.
|
[1] | ZHU Jun, ZHANG Tianyi, XIE Yakun, ZHANG Jie, LI Chuangnong, ZHAO Li, LI Weilian. Intelligent Statistic Method for Video Pedestrian Flow Considering Small Object Features[J]. Journal of Southwest Jiaotong University, 2022, 57(4): 705-712, 736. doi: 10.3969/j.issn.0258-2724.20200425 |
[2] | WANG Tao, ZENG Wenhao, YU Qi. Obstacle Detection and Tracking for Driverless Cars Based on Lidar[J]. Journal of Southwest Jiaotong University, 2021, 56(6): 1346-1354. doi: 10.3969/j.issn.0258-2724.20200240 |
[3] | ZHU Jun, SHE Ping, LI Weilian, CAO Yungang, QI Hua, WANG Bo, WANG Yu. Dynamic Planning Method for Indoor-Fire Escape Path Based on Navigation Grid[J]. Journal of Southwest Jiaotong University, 2020, 55(5): 1103-1110. doi: 10.3969/j.issn.0258-2724.20180912 |
[4] | ZHOU Yan, CHEN Hong, ZHANG Yeting, HUANG Yueying, ZHANG Pengcheng, YANG Weijun. Multi-objective Indoor Path Planning Method with Dynamic Environment Awareness[J]. Journal of Southwest Jiaotong University, 2019, 54(3): 611-618, 632. doi: 10.3969/j.issn.0258-2724.20180098 |
[5] | LU Zhiwen, DONG Dawei, HUA Chunrong, YAN Bing, XIE Xiao, CHEN Jun. Localization Method of Seismic Target in Anisotropic Ground[J]. Journal of Southwest Jiaotong University, 2015, 28(3): 442-448. doi: 10.3969/j.issn.0258-2724.2015.03.009 |
[6] | DENG Ping, WU Xiaomei. An Adaptive Localization Algorithm Based on Particle Filter for Wireless Sensor Networks[J]. Journal of Southwest Jiaotong University, 2014, 27(2): 323-329. doi: 10.3969/j.issn.0258-2724.2014.02.021 |
[7] | XIA Chuangwen, XU Jianmin, LU Yijie, WANG Qionghua. Video-Based Vehicle Tracking and Capturing System for Expressway Tollgates[J]. Journal of Southwest Jiaotong University, 2013, 26(2): 350-356. doi: 10.3969/j.issn.0258-2724.2013.02.025 |
[8] | MA Lili, CHENG Wenming, ZHONG Bin. Research on Load's Tracking Control for Crane's Motor Driving System[J]. Journal of Southwest Jiaotong University, 2013, 26(3): 494-499. doi: 10.3969/j.issn.0258-2724.2013.03.016 |
[9] | XU Yantian, CHENG Pengfei, CAI Yanhui. Kalman Filter Algorithm for Medium-Range Real-Time Kinematic Positioning with One Reference Station[J]. Journal of Southwest Jiaotong University, 2013, 26(2): 317-322,356. doi: 10.3969/j.issn.0258-2724.2013.02.020 |
[10] | CUI Bo, ZHANG Jiashu, YANG Yu. Nonlinear Target Tracking Algorithm Based on Block Ensemble Kalman Filter[J]. Journal of Southwest Jiaotong University, 2013, 26(5): 863-869. doi: 10.3969/j.issn.0258-2724.2013.05.013 |
[11] | BAO Haichao, YANG Genqing, ZHANG Yuning. On-Board Software Self-recovery Technique for Defending MBU Effect[J]. Journal of Southwest Jiaotong University, 2012, 25(4): 669-674. doi: 10.3969/j.issn.0258-2724.2012.04.021 |
[12] | MA Lei, SHI Xizhi. Recent Development on Consensus-Based Kalman Filtering in Multi-agent Systems[J]. Journal of Southwest Jiaotong University, 2011, 24(2): 287-293. doi: 10.3969/j.issn.0258-2724.2011.02.019 |
[13] | HU Lu, PENG Bo, JIANG Yangsheng, ZHANG Gai. Multi-objective Optimization of Cross-Docking Logistics System under Stochastic Demand[J]. Journal of Southwest Jiaotong University, 2011, 24(4): 688-694. doi: 10.3969/j.issn.0258-2724.2011.04.026 |
[14] | ZHONG Xuanming, LIAO Cheng. Microwave Imaging of Perfectly Conducting Cylinder Based on Time Reversal and Single Antenna[J]. Journal of Southwest Jiaotong University, 2011, 24(3): 451-455. doi: 10.3969/j.issn.0258-2724.2011.03.016 |
[15] | JIANG Lian-Xiang, YANG Gen-Qiang, LI Hua-Wang. Dynamic Power Management Based on Subtask and Executing Time[J]. Journal of Southwest Jiaotong University, 2010, 23(3): 435-439. doi: 10. 3969/ j. issn. 0258-2724. |
[16] | ZHOU Xiu-Yun, HUANG Jian-Guo, FENG Zhong-Zheng. Polycrystalline Silicon Video Inspection System Based on Small Area Tracking[J]. Journal of Southwest Jiaotong University, 2010, 23(3): 430-434. doi: 10. 3969/ j. issn. 0258-2724. |
[17] | KOU Yimin, MA Guangcheng, WANG Changhong. Solution to Geomagnetic Systems in Areas with Smooth Magnetic Fields[J]. Journal of Southwest Jiaotong University, 2009, 22(6): 951-957. |
[18] | ZHANG Jianyong, LI Jun, GUO Yaohuang. Insertion Heuristic Algorithm for Dynamic Vehicle Routing Problem with Fuzzy Due-Time[J]. Journal of Southwest Jiaotong University, 2008, 21(1): 107-113. |
[19] | ZHAO Dong-mei, TAO Zhang-hua. Method of Noninferior Matrix of Multi-objective Dynamic Programming with Indefinite Phases[J]. Journal of Southwest Jiaotong University, 2003, 16(6): 677-681. |
[20] | HE Rong, FANG Xu-ming. An Integrated Delivery Tracking System for E-Commerce[J]. Journal of Southwest Jiaotong University, 2002, 15(2): 167-171. |
参数 | 处理方案 |
测站坐标 | 静态时过程噪声为 0,动态时先验值源自伪距单点定位 |
模糊度 | 过程噪声为 0 |
对流层 | 随机游走,过程噪声默认为 4 cm2/h,动态时顾及历元间高差引起的影响 |
对流层梯度 | 随机游走 |
电离层 | 随机游走 |
电离层梯度 | 随机游走 |
接收机钟差 | 先验值源自伪距单点定位 |
卫星钟差 | 20 min弧段二阶多项式拟合 |
卫星坐标 | 10 阶滑动拉格朗日插值 |
站点 | 经度/(°) | 纬度/(°) |
CHAN(长春) | 125.443 | 43.790 |
BJFS(北京房山) | 115.892 | 39.609 |
BJNM(北京) | 116.224 | 40.245 |
HKSL(香港小冷水) | 113.928 | 22.372 |
HKWS(香港黄石) | 114.335 | 22.434 |
JFNG(九峰) | 114.491 | 30.516 |
LHAZ(拉萨) | 91.104 | 29.657 |
TWTF(桃园) | 121.164 | 24.954 |
站点 | 年积日 | dN /mm | dE /mm | dU /mm | 收敛时间 /min |
CHAN | 001 | 1.9 | −9.2 | −11.3 | 32.0 |
002 | 2.4 | −5.7 | −7.5 | 19.0 | |
003 | 2.0 | −7.1 | −9.2 | 43.0 | |
BJFS | 001 | 2.0 | 1.2 | −1.8 | 30.0 |
002 | 2.5 | 0.7 | −0.5 | 21.0 | |
003 | 2.5 | −0.2 | 2.8 | 23.0 | |
BJNM | 001 | 2.6 | −3.6 | −0.8 | 29.0 |
002 | 4.5 | 0.4 | 1.4 | 50.0 | |
003 | 4.0 | −4.5 | 4.3 | 46.5 | |
HKSL | 001 | 2.0 | 0.5 | −3.9 | 30.5 |
002 | 3.8 | 0.1 | 1.7 | 25.0 | |
003 | 4.4 | −2.5 | 2.0 | 20.5 | |
HKWS | 001 | 1.0 | 1.5 | 2.2 | 30.0 |
002 | 3.3 | 1.9 | −0.1 | 52.5 | |
003 | 3.5 | −0.3 | 1.9 | 104.5 | |
JFNG | 001 | 2.0 | 3.2 | 0.6 | 20.0 |
002 | 2.8 | 4.4 | 9.0 | 13.0 | |
003 | 3.1 | 3.4 | 8.1 | 19.0 | |
LHAZ | 001 | 3.6 | −1.5 | −3.8 | 42.0 |
002 | 4.1 | 1.5 | −0.3 | 57.0 | |
003 | 3.5 | −1.5 | −4.4 | 44.0 | |
TWTF | 001 | 0.8 | 3.9 | −9.3 | 26.5 |
002 | 4.4 | 1.0 | −5.4 | 21.5 | |
003 | 4.0 | 0.1 | −6.3 | 19.0 |
参数 | 处理方案 |
测站坐标 | 静态时过程噪声为 0,动态时先验值源自伪距单点定位 |
模糊度 | 过程噪声为 0 |
对流层 | 随机游走,过程噪声默认为 4 cm2/h,动态时顾及历元间高差引起的影响 |
对流层梯度 | 随机游走 |
电离层 | 随机游走 |
电离层梯度 | 随机游走 |
接收机钟差 | 先验值源自伪距单点定位 |
卫星钟差 | 20 min弧段二阶多项式拟合 |
卫星坐标 | 10 阶滑动拉格朗日插值 |
站点 | 经度/(°) | 纬度/(°) |
CHAN(长春) | 125.443 | 43.790 |
BJFS(北京房山) | 115.892 | 39.609 |
BJNM(北京) | 116.224 | 40.245 |
HKSL(香港小冷水) | 113.928 | 22.372 |
HKWS(香港黄石) | 114.335 | 22.434 |
JFNG(九峰) | 114.491 | 30.516 |
LHAZ(拉萨) | 91.104 | 29.657 |
TWTF(桃园) | 121.164 | 24.954 |
站点 | 年积日 | dN /mm | dE /mm | dU /mm | 收敛时间 /min |
CHAN | 001 | 1.9 | −9.2 | −11.3 | 32.0 |
002 | 2.4 | −5.7 | −7.5 | 19.0 | |
003 | 2.0 | −7.1 | −9.2 | 43.0 | |
BJFS | 001 | 2.0 | 1.2 | −1.8 | 30.0 |
002 | 2.5 | 0.7 | −0.5 | 21.0 | |
003 | 2.5 | −0.2 | 2.8 | 23.0 | |
BJNM | 001 | 2.6 | −3.6 | −0.8 | 29.0 |
002 | 4.5 | 0.4 | 1.4 | 50.0 | |
003 | 4.0 | −4.5 | 4.3 | 46.5 | |
HKSL | 001 | 2.0 | 0.5 | −3.9 | 30.5 |
002 | 3.8 | 0.1 | 1.7 | 25.0 | |
003 | 4.4 | −2.5 | 2.0 | 20.5 | |
HKWS | 001 | 1.0 | 1.5 | 2.2 | 30.0 |
002 | 3.3 | 1.9 | −0.1 | 52.5 | |
003 | 3.5 | −0.3 | 1.9 | 104.5 | |
JFNG | 001 | 2.0 | 3.2 | 0.6 | 20.0 |
002 | 2.8 | 4.4 | 9.0 | 13.0 | |
003 | 3.1 | 3.4 | 8.1 | 19.0 | |
LHAZ | 001 | 3.6 | −1.5 | −3.8 | 42.0 |
002 | 4.1 | 1.5 | −0.3 | 57.0 | |
003 | 3.5 | −1.5 | −4.4 | 44.0 | |
TWTF | 001 | 0.8 | 3.9 | −9.3 | 26.5 |
002 | 4.4 | 1.0 | −5.4 | 21.5 | |
003 | 4.0 | 0.1 | −6.3 | 19.0 |