CN107247275A - 基于公交车的城市gnss脆弱性监测系统及其方法 - Google Patents

基于公交车的城市gnss脆弱性监测系统及其方法 Download PDF

Info

Publication number
CN107247275A
CN107247275A CN201710390589.4A CN201710390589A CN107247275A CN 107247275 A CN107247275 A CN 107247275A CN 201710390589 A CN201710390589 A CN 201710390589A CN 107247275 A CN107247275 A CN 107247275A
Authority
CN
China
Prior art keywords
mrow
msub
msubsup
mtr
mtd
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710390589.4A
Other languages
English (en)
Other versions
CN107247275B (zh
Inventor
战兴群
刘宝玉
刘金玉
朱兵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201710390589.4A priority Critical patent/CN107247275B/zh
Publication of CN107247275A publication Critical patent/CN107247275A/zh
Application granted granted Critical
Publication of CN107247275B publication Critical patent/CN107247275B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/23Testing, monitoring, correcting or calibrating of receiver elements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/18Stabilised platforms, e.g. by gyroscope

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Navigation (AREA)

Abstract

一种基于公交车的城市GNSS脆弱性监测系统及其方法,包括:惯性导航模块、GNSS导航模块、地图匹配导航模块、组合导航分析模块、GNSS脆弱性分析模块和网络通信模块,本发明能充分利用城市公交车这一大众交通工具及其车载设备,实施成本低,采集的数据实时性好、样本数量多、覆盖范围广,能得到可靠、准确、实时的统计分析结果,充分结合公交车的行车特点,利用里程计和电子地图进行特征点匹配和航位推测导航,即使GNSS导航失效、受到干扰或欺骗时,也能保证系统有效地进行导航,从而实现连续的监测,更加准确地计算出公交车经受的GNSS脆弱性。

Description

基于公交车的城市GNSS脆弱性监测系统及其方法
技术领域
本发明涉及的是一种交通运输领域的技术,具体是一种基于公交车的城市GNSS脆弱性 监测系统及其方法。
背景技术
全球卫星导航系统(GNSS)能全球、全天候、实时地提供定位、测速和授时等服务。在高 楼林立、树木茂密、电磁辐射污染较严重的城市区域中,GNSS信号很容易受到遮挡、多径、 干扰等影响,使城市GNSS导航的脆弱性大大提高。GNSS脆弱性监测是监测GNSS系统内外因 素对GNSS服务性能的不利影响,其结果可指示GNSS服务性能受不利影响的程度。为了提高 城市GNSS应用的性能和可靠性,有必要监测、掌握城市GNSS脆弱性。而现在对城市GNSS 导航的监测主要是基于设立的少数固定监测站点,相关科学研究也是通过专用实验车辆对城市 少数特定的街道进行GNSS导航数据的采样与分析,都存在覆盖面小、采样样本数有限、实时 性不足等缺点,无法对城市GNSS脆弱性进行实时、全面的监测。
发明内容
本发明针对现有技术或没有考虑和分析GNSS脆弱性的影响、或没有结合地图信息和里 程计对惯性导航进行反馈校正,导致其在卫星导航及其差分定位被欺骗或不可用时,将导致系 统输出的定位结果出现严重偏差或随惯性导航误差的累积而发散等缺陷,提出一种基于公交车 的城市GNSS脆弱性监测系统及其方法。
本发明是通过以下技术方案实现的:
本发明涉及一种基于公交车的城市GNSS脆弱性监测系统,包括:惯性导航模块、GNSS 导航模块、地图匹配导航模块、组合导航分析模块、GNSS脆弱性分析模块和网络通信模块, 其中:GNSS导航模块测得公交车位置GNSS导航坐标[Lg,λg,hg],当地东北天地理坐标系 OXEYNZU中GNSS导航速度和伪距伪距率惯性导航模块采集公交车的加速度信息和角速度信息通过捷联惯性导航算法得到惯性导航坐 标[Lins,λins,hins]和惯性导航速度地图匹配导航模块采集行车信息后通过地图匹 配算法得到地图匹配导航坐标[Lodon,λodon,hodon]和公交车本体坐标系OXBYBZB中地图匹配导航速 度组合导航分析模块接收惯性导航模块、GNSS导航模块和地图匹配导航模 块的坐标和速度等得到公交车的估计坐标[Lfinal,λfinal,hfinal]和估计速度GNSS 脆弱性分析模块根据估计坐标[Lfinal,λfinal,hfinal]和估计速度计算GNSS脆弱性 参数并通过网络通信模块传送到监测中心。
所述的地图匹配算法是指采集公交车与设置于公交车运行路线中识别点k之间的距离 和速度比照公交车在电子地图中运行路线获得地图匹配导航坐标[Lodon,λodon,hodon]和地 图匹配导航速度
所述的识别点分别设置于公交车站台、道路转弯处和道路陡坡处。
所述的距离通过公式进行校正,速度通过公式进行校正。
所述的组合导航分析模块首先对惯性导航坐标[Lins,λins,hins]、惯性导航速度地图匹配导航坐标[Lodon,λodon,hodon]和地图匹配导航速度进行松 组合导航滤波,对GNSS伪距伪距率惯性导航坐标[Lins,λins,hins]和惯性导航速度进行紧组合导航滤波,再经联邦滤波得到估计坐标 [Lfinal,λfinal,hfinal]和估计速度
所述的松组合导航滤波的观测方程为:
分别为惯性导航航向角和根据磁强计输出估计的航向角,δL,δλ,δh(纬度误差、 经度误差、高度误差)、δVE,δVN,δVU(坐标系OXEYNZU下的速度误差)、φE,φN,φU(姿态误差)为 组合导航滤波状态量中包含的惯性导航误差状态。
所述的GNSS脆弱性参数包括伪距误差方差、GNSS信号功率、载噪比、DOP值、定位误差方差。
所述的伪距误差方差其中:ρfinal为通过估计坐标[Lfinal,λfinal,hfinal]计算的第i颗GNSS卫星的伪距值。
本发明涉及一种基于上述系统的公交车的城市GNSS脆弱性监测方法,包括以下步骤:
1)测量公交车位置GNSS导航坐标[Lg,λg,hg]和坐标系OXEYNZU中GNSS导航速度和伪距伪距率
2)采集公交车的加速度信息和角速度信息通过捷联惯性导航算法得到惯性导航坐标 [Lins,λins,hins]和惯性导航速度
3)通过地图匹配算法得到公交车的地图匹配导航坐标[Lodon,λodon,hodon]和地图匹配导航 速度
4)对惯性导航坐标[Lins,λins,hins]、惯性导航速度地图匹配导航坐标[Lodon,λodon,hodon]和地图匹配导航速度进行松组合导航滤波,对GNSS伪距伪距率惯性导航坐标[Lins,λins,hins]和惯性导航速度进行紧组合导航滤波,再经联邦滤波得到估计坐标[Lfinal,λfinal,hfinal]和估计速度
5)计算GNSS脆弱性参数并通过网络通信模块传送到监测中心。
技术效果
与现有技术相比,本发明能充分利用城市公交车这一大众交通工具及其车载设备,实施 成本低,采集的数据实时性好、样本数量多、覆盖范围广,能得到可靠、准确、实时的统计分 析结果,充分结合公交车的行车特点,利用里程计和电子地图进行特征点匹配和航位推测导航, 即使在GNSS导航不可用、受到干扰或欺骗时,也能保证系统有效地进行导航,从而实现连续 的监测,更加准确地计算出公交车经受的GNSS脆弱性。
附图说明
图1为本发明系统示意图;
图2为实施例坐标示意图;
图3为实施例实现场景示意图。
具体实施方式
如图1所示,本实施例中的基于公交车的城市GNSS脆弱性监测系统,包括:惯性导航 模块、GNSS导航模块、地图匹配导航模块、组合导航分析模块、GNSS脆弱性分析模块和网络 通信模块,其中:GNSS导航模块测得公交车位置GNSS导航坐标[Lg,λg,hg]和坐标系OXEYNZU中GNSS导航速度和伪距伪距率惯性导 航模块采集公交车的加速度信息和角速度信息通过捷联惯性导航算法得到惯性导航坐标 [Lins,λins,hins]和惯性导航速度地图匹配导航模块采集行车信息后通过地图匹配 算法得到地图匹配导航坐标[Lodon,λodon,hodon]和地图匹配导航速度组合导航 分析模块接收惯性导航模块、GNSS导航模块和地图匹配导航模块的坐标和速度等得到公交车 的估计坐标[Lfinal,λfinal,hfinal]和估计速度GNSS脆弱性分析模块根据估计坐标 [Lfinal,λfinal,hfinal]和估计速度计算GNSS脆弱性参数并通过网络通信模块传送 到监测中心。
所述的公交车上设有MEMS传感器,加速度计输出陀螺仪输出磁强计输出磁北角ψmag。MEMS传感器三轴方向与公交车本体坐标系OXBYBZB三轴方向平行,OXB轴指向公交车正前方,OZB轴指向公交车上方,OYB轴与OXB轴、OZB轴形 成右手直角坐标系。根据公交车的行车信息即里程计和车门开关状态量,可以得到车辆航向角的变化量车辆行驶距离的变化量实时车速车辆航向角为在当地东北天地理坐标系OXEYNZU中的XEYN平面内,公交车轴线OXB的投影与OYN轴的夹角,北偏东为正,坐 标系原点O为MEMS传感器中心位置。
所述的GNSS导航模块接收GNSS信号,解算得到GNSS导航坐标[Lg,λg,hg]和坐标系OXEYNZU中GNSS导航速度和伪距伪距率
所述的惯性导航模块接收MEMS传感器的加速度信息和角速度信息运行捷联惯性导航算法获得公交车惯性导航坐标[Lins,λins,hins]和坐标系OXEYNZU中的惯性导航速度以及姿态信息。姿态信息包括俯仰角滚动角和车辆 航向角
所述的地图匹配导航模块采集行车信息后通过地图匹配算法得到地图匹配导航坐标 [Lodon,λodon,hodon]和地图匹配导航速度所述的地图匹配算法是指采集公交车 与设置于公交车运行路线中识别点k之间的距离和速度比照公交车在电子地图中运行 路线获得地图匹配导航坐标[Lodon,λodon,hodon]和地图匹配导航速度
在公交车的运行线路上依次设置识别点,并编号为识别点设置于公交站 台、道路转弯处或道路陡坡处。公交车行驶在特征点Pk-1与Pk之间时,因为公交车的行车线路 固定,且线路上各点的坐标值已知,根据可知公交车相对于特征点Pk-1在其行车线路上已行 驶的距离,结合公交车电子地图可得到公交车在其行车线路上的位置,得到地图匹配导航坐标 [Lodon,λodon,hodon]和地图匹配导航速度 为公交车从识别点Pk-1开始的里程计 的累加之和。
所述的距离通过公式进行校正,速度通过公式进行校正。kodo为里程计输出的误差因子,为噪声量,主要为道路路面平整度、车轮滑行和打滑等对的影响,公交车里程计车速,Sk为对应的无误差理想值。
所述的组合导航分析模块首先对惯性导航坐标[Lins,λins,hins]、惯性导航速度地图匹配导航坐标[Lodon,λodon,hodon]和地图匹配导航速度进行松 组合导航滤波,对GNSS伪距伪距率惯性导航坐标[Lins,λins,hins]和惯性导航速度进行紧组合导航滤波,再经联邦滤波得到估计坐标 [Lfinal,λfinal,λfinal]和估计速度
所述的松组合导航滤波观测方程为:
,其中:RM为当地子午圈曲率半径,RN为当地卯酉圈曲率半径,ωloose为测量噪声, C={Cij},i,j=1,2,3为公交车本体坐标系OXBYBZB至坐标系OXEYNZU的转换矩阵。再经Kalman 滤波器进行滤波,得到惯性导航模块和地图匹配导航模块的松组合导航滤波结果。
如图2所示,所述的MEMS传感器中心位置取为坐标系OXBYBZB与坐标系OXEYNZU的坐标原点,其中为公交车所在道路的前行方向在坐标平面XEYN的投影与OYN轴的夹角,坐标轴OYN绕OZU轴旋转后与该投影重合,为公交车所在道路 与坐标平面XEYN平面的夹角。当公交车行驶在识别点之间时,可认为公交车行的行车方向与行车道路保持一致,此时,地图匹配导航模块得到的公交车速度估计值在坐标系OXBYBZB中可表 示为其估计的公交车位置误差及里程计输出的车速误差体现在公交车行车线路方 向上分别为对于里程计,kodo一般是一个较小的量,当kodo被不断估计和 校正后,其残差δkodo可被认为是一个小量,里程计输出的位置误差和车速误差可分别表示为
所述的紧组合导航滤波的观测值为未接收到各GNSS导航卫星的伪距校正量时,观测值为:其中:为基于惯性导航结果计算的第i颗GNSS卫星的伪距值,为基于惯性导航结果计算的第i颗GNSS卫星的伪距率值。假定第i颗GNSS卫星的高度角为且相对于坐标系OXBYBZB中OXB轴的方位角,指向OXB轴为0°,顺时钟为正,为观测值的测量方差设定为其中为公交车运营线路的电子地图中对应于该位置处第i颗GNSS卫星的脆弱性信息中的伪距误差方差。
经过松组合导航滤波和紧组合导航滤波,再通过联邦滤波最终得到估计坐标[Lfinal,λfinal,hfinal]和估计速度
所述的GNSS脆弱性分析模块根据估计坐标[Lfinal,λfinal,hfinal]和估计速度对GNSS脆弱性参数进行分析。GNSS脆弱性参数包括伪距误差方差、GNSS 信号功率、载噪比、DOP值、定位误差方差、定位异常率和可用性。
所述的伪距误差方差其中:ρfinal为通过估计坐标[Lfinal,λfinal,hfinal]计算的第i颗GNSS卫星的伪距值。定位误差为 ELλh=[Lg,λg,hg]T-[Lfinal,λfinal,hfinal]T。对于信号功率、载噪比和DOP值可通过相应公式得到。 用新得到的GNSS脆弱性参数测量值和α滤波器对GNSS脆弱性参数进行更新。GNSS定位的可 用性定义为该位置处GNSS定位误差小于一定阀值的定位次数占总定位次数的比率,定位异常 率定义为该位置处GNSS定位误差大于一定阀值的定位次数占总定位次数的比率。
本实施例涉及一种基于公交车的城市GNSS脆弱性监测方法,包括以下步骤:
1)测量公交车位置GNSS导航坐标[Lg,λg,hg]和坐标系OXEYNZU中GNSS导航速度和伪距伪距率
2)采集公交车的加速度信息和角速度信息通过捷联惯性导航算法得到惯性导航坐标 [Lins,λins,hins]和惯性导航速度
3)通过地图匹配算法得到公交车的地图匹配导航坐标[Lodon,λodon,hodon]和地图匹配导航 速度
4)对惯性导航坐标[Lins,λins,hins]、惯性导航速度地图匹配导航坐标[Lodon,λodon,hodon]和地图匹配导航速度进行松组合导航滤波,对GNSS伪距伪距率惯性导航坐标[Lins,λins,hins]和惯性导航速度进行紧组合导航滤波,再经联邦滤波得到估计坐标[Lfinal,λfinal,hfinal]和估计速度
5)计算GNSS脆弱性参数并通过网络通信模块传送到分析监测中心。
如图3所示,分析监测中心根据接收到各公交车的组合导航定位结果、GNSS脆弱性监 测信息和存储数据后综合分析得到整个城市的GNSS脆弱性,并向用户发布GNSS脆弱性监测 信息。
与现有技术相比,本发明能充分利用城市公交车这一大众交通工具及其车载设备,实施 成本低,采集的数据实时性好、样本数量多、覆盖范围广,能得到可靠、准确、实时的统计分 析结果,充分结合公交车的行车特点,利用里程计和电子地图进行特征点匹配和航位推测导航, 即使在GNSS导航不可用、受到干扰或欺骗时,也能保证系统有效地进行导航,从而实现连续 的监测,更加准确地计算出公交车经受的GNSS脆弱性。
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式 对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围 内的各个实现方案均受本发明之约束。

Claims (8)

1.一种基于公交车的城市GNSS脆弱性监测系统,其特征在于,包括:惯性导航模块、GNSS导航模块、地图匹配导航模块、组合导航分析模块、GNSS脆弱性分析模块和网络通信模块,其中:GNSS导航模块测得公交车位置GNSS导航坐标[Lg,λg,hg]和坐标系OXEYNZU中GNSS导航速度和伪距伪距率惯性导航模块采集公交车的加速度信息和角速度信息通过捷联惯性导航算法得到惯性导航坐标[Lins,λins,hins]和惯性导航速度地图匹配导航模块采集行车信息后通过地图匹配算法得到地图匹配导航坐标和地图匹配导航速度组合导航分析模块接收惯性导航模块、GNSS导航模块和地图匹配导航模块的坐标和速度等得到公交车的估计坐标[Lfinal,λfinal,hfinal]和估计速度GNSS脆弱性分析模块根据估计坐标[Lfinal,λfinal,Lfinal]和估计速度计算GNSS脆弱性参数并通过网络通信模块传送到监测中心。
2.根据权利要求1所述的基于公交车的城市GNSS脆弱性监测系统,其特征是,所述的地图匹配算法是指采集公交车与设置于公交车运行路线中识别点k之间的距离和速度比照公交车在电子地图中运行路线获得地图匹配导航坐标[Lodon,λodon,hodon]和地图匹配导航速度
3.根据权利要求2所述的基于公交车的城市GNSS脆弱性监测系统,其特征是,所述的距离通过公式进行校正,速度通过公式进行校正。
4.根据权利要求3所述的基于公交车的城市GNSS脆弱性监测系统,其特征是,所述的组合导航分析模块首先对惯性导航坐标[Lins,λins,hins]、惯性导航速度地图匹配导航坐标[Lodon,λodon,hodon]和地图匹配导航速度进行松组合导航滤波,对GNSS伪距伪距率惯性导航坐标[Lins,λins,hins]和惯性导航速度进行紧组合导航滤波,再经联邦滤波得到估计坐标[Lfinal,λfinal,hfinal]和估计速度
5.根据权利要求4所述的基于公交车的城市GNSS脆弱性监测系统,其特征是,所述的松组合导航滤波的观测方程为:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>L</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>s</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;lambda;</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>s</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>h</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>s</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>V</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>s</mi> </mrow> <mi>E</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>V</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>s</mi> </mrow> <mi>N</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>V</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>s</mi> </mrow> <mi>U</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>s</mi> </mrow> <mrow> <mi>y</mi> <mi>a</mi> <mi>w</mi> </mrow> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>s</mi> </mrow> <mrow> <mi>y</mi> <mi>a</mi> <mi>w</mi> </mrow> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>L</mi> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> <mi>n</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;lambda;</mi> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> <mi>n</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>h</mi> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> <mi>n</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <msub> <mover> <mi>V</mi> <mo>^</mo> </mover> <mi>k</mi> </msub> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> </mrow> </msup> <mi>cos</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>p</mi> <mi>i</mi> <mi>t</mi> <mi>c</mi> <mi>h</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mi>sin</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>y</mi> <mi>a</mi> <mi>w</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <msub> <mover> <mi>V</mi> <mo>^</mo> </mover> <mi>k</mi> </msub> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> </mrow> </msup> <mi>cos</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>p</mi> <mi>i</mi> <mi>t</mi> <mi>c</mi> <mi>h</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mi>cos</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>y</mi> <mi>a</mi> <mi>w</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <msub> <mover> <mi>V</mi> <mo>^</mo> </mover> <mi>k</mi> </msub> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> </mrow> </msup> <mi>sin</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>p</mi> <mi>i</mi> <mi>t</mi> <mi>c</mi> <mi>h</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>y</mi> <mi>a</mi> <mi>w</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>&amp;psi;</mi> <mo>^</mo> </mover> <mrow> <mi>m</mi> <mi>a</mi> <mi>g</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>&amp;delta;</mi> <mi>L</mi> <mo>-</mo> <msup> <mi>&amp;delta;k</mi> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> </mrow> </msup> <msub> <mover> <mi>S</mi> <mo>^</mo> </mover> <mi>k</mi> </msub> <mi>cos</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>p</mi> <mi>i</mi> <mi>t</mi> <mi>c</mi> <mi>h</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mi>cos</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>y</mi> <mi>a</mi> <mi>w</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>/</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>R</mi> <mi>M</mi> </msub> <mo>+</mo> <mi>h</mi> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;delta;</mi> <mi>&amp;lambda;</mi> <mo>-</mo> <msup> <mi>&amp;delta;k</mi> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> </mrow> </msup> <msub> <mover> <mi>S</mi> <mo>^</mo> </mover> <mi>k</mi> </msub> <mi>cos</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>p</mi> <mi>i</mi> <mi>t</mi> <mi>c</mi> <mi>h</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mi>sin</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>y</mi> <mi>a</mi> <mi>w</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mi>sec</mi> <mrow> <mo>(</mo> <msub> <mi>L</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>s</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>R</mi> <mi>N</mi> </msub> <mo>+</mo> <mi>h</mi> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;delta;</mi> <mi>h</mi> <mo>-</mo> <msup> <mi>&amp;delta;k</mi> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> </mrow> </msup> <msub> <mover> <mi>S</mi> <mo>^</mo> </mover> <mi>k</mi> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>p</mi> <mi>i</mi> <mi>t</mi> <mi>c</mi> <mi>h</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;delta;V</mi> <mi>E</mi> </msub> <mo>-</mo> <msup> <mi>&amp;delta;k</mi> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> </mrow> </msup> <msup> <msub> <mover> <mi>V</mi> <mo>^</mo> </mover> <mi>k</mi> </msub> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> </mrow> </msup> <mi>cos</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>p</mi> <mi>i</mi> <mi>t</mi> <mi>c</mi> <mi>h</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mi>sin</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>y</mi> <mi>a</mi> <mi>w</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;delta;V</mi> <mi>N</mi> </msub> <mo>-</mo> <msup> <mi>&amp;delta;k</mi> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> </mrow> </msup> <msup> <msub> <mover> <mi>V</mi> <mo>^</mo> </mover> <mi>k</mi> </msub> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> </mrow> </msup> <mi>cos</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>p</mi> <mi>i</mi> <mi>t</mi> <mi>c</mi> <mi>h</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mi>cos</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>y</mi> <mi>a</mi> <mi>w</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;delta;V</mi> <mi>U</mi> </msub> <mo>-</mo> <msup> <mi>&amp;delta;k</mi> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> </mrow> </msup> <msup> <msub> <mover> <mi>V</mi> <mo>^</mo> </mover> <mi>k</mi> </msub> <mrow> <mi>o</mi> <mi>d</mi> <mi>o</mi> </mrow> </msup> <mi>sin</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>p</mi> <mi>i</mi> <mi>t</mi> <mi>c</mi> <mi>h</mi> </mrow> <mrow> <mi>l</mi> <mi>o</mi> <mi>a</mi> <mi>d</mi> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mfrac> <mrow> <msub> <mi>C</mi> <mn>12</mn> </msub> <msub> <mi>C</mi> <mn>32</mn> </msub> </mrow> <mrow> <msubsup> <mi>C</mi> <mn>12</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>C</mi> <mn>22</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> <msub> <mi>&amp;phi;</mi> <mi>E</mi> </msub> <mo>-</mo> <mfrac> <mrow> <msub> <mi>C</mi> <mn>22</mn> </msub> <msub> <mi>C</mi> <mn>32</mn> </msub> </mrow> <mrow> <msubsup> <mi>C</mi> <mn>12</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>C</mi> <mn>22</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> <msub> <mi>&amp;phi;</mi> <mi>N</mi> </msub> <mo>+</mo> <msub> <mi>&amp;phi;</mi> <mi>U</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mfrac> <mrow> <msub> <mi>C</mi> <mn>12</mn> </msub> <msub> <mi>C</mi> <mn>32</mn> </msub> </mrow> <mrow> <msubsup> <mi>C</mi> <mn>12</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>C</mi> <mn>22</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> <msub> <mi>&amp;phi;</mi> <mi>E</mi> </msub> <mo>-</mo> <mfrac> <mrow> <msub> <mi>C</mi> <mn>22</mn> </msub> <msub> <mi>C</mi> <mn>32</mn> </msub> </mrow> <mrow> <msubsup> <mi>C</mi> <mn>12</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>C</mi> <mn>22</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> <msub> <mi>&amp;phi;</mi> <mi>N</mi> </msub> <mo>+</mo> <msub> <mi>&amp;phi;</mi> <mi>U</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <msub> <mi>&amp;omega;</mi> <mrow> <mi>l</mi> <mi>o</mi> <mi>o</mi> <mi>s</mi> <mi>e</mi> </mrow> </msub> </mrow>
分别为惯性导航航向角和根据磁强计输出估计的航向角。
6.根据权利要求5所述的基于公交车的城市GNSS脆弱性监测系统,其特征是,所述的GNSS脆弱性参数包括伪距误差方差、GNSS信号功率、载噪比、DOP值、定位误差方差。
7.根据权利要求6所述的基于公交车的城市GNSS脆弱性监测系统,其特征是,所述的伪距误差方差其中:ρfinal为通过估计坐标[Lfinal,λfinal,hfinal]计算的第i颗GNSS卫星的伪距值。
8.一种基于公交车的城市GNSS脆弱性监测方法,其特征在于,包括以下步骤:
1)测量公交车位置GNSS导航坐标[Lg,λg,hg]和坐标系OXEYNZU中GNSS导航速度和伪距伪距率
2)采集公交车的加速度信息和角速度信息通过捷联惯性导航算法得到惯性导航坐标[Lins,λins,hins]和惯性导航速度
3)通过地图匹配算法得到公交车的地图匹配导航坐标[Lodon,hodon,hodon]和地图匹配导航速度
4)对惯性导航坐标[Lins,λins,hins]、惯性导航速度地图匹配导航坐标[Lodon,λodon,hodon]和地图匹配导航速度进行松组合导航滤波,对GNSS伪距伪距率惯性导航坐标[Lins,λins,hins]和惯性导航速度进行紧组合导航滤波,再经联邦滤波得到估计坐标[Lfinal,λfinal,hfinal]和估计速度
5)计算GNSS脆弱性参数并通过网络通信模块传送到监测中心。
CN201710390589.4A 2017-05-27 2017-05-27 基于公交车的城市gnss脆弱性监测系统及其方法 Active CN107247275B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710390589.4A CN107247275B (zh) 2017-05-27 2017-05-27 基于公交车的城市gnss脆弱性监测系统及其方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710390589.4A CN107247275B (zh) 2017-05-27 2017-05-27 基于公交车的城市gnss脆弱性监测系统及其方法

Publications (2)

Publication Number Publication Date
CN107247275A true CN107247275A (zh) 2017-10-13
CN107247275B CN107247275B (zh) 2020-12-08

Family

ID=60017662

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710390589.4A Active CN107247275B (zh) 2017-05-27 2017-05-27 基于公交车的城市gnss脆弱性监测系统及其方法

Country Status (1)

Country Link
CN (1) CN107247275B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108051839A (zh) * 2017-10-27 2018-05-18 成都天合世纪科技有限责任公司 一种车载三维定位装置及三维定位的方法
CN108846598A (zh) * 2018-03-29 2018-11-20 宏图物流股份有限公司 一种车辆定位的方法及装置
CN108920892A (zh) * 2018-09-30 2018-11-30 南京地铁集团有限公司 一种城市轨道交通站点脆弱性测量方法
CN108931260A (zh) * 2018-07-04 2018-12-04 山东省科学院自动化研究所 Gps欺骗时无人驾驶车辆定位安全测试系统及方法
CN109631886A (zh) * 2018-12-28 2019-04-16 江苏满运软件科技有限公司 车辆定位方法、装置、电子设备、存储介质
CN110795518A (zh) * 2019-10-24 2020-02-14 中国人民解放军63653部队 一种基于vc++的gnss定位导航软件系统
CN111323802A (zh) * 2020-03-20 2020-06-23 北京百度网讯科技有限公司 车辆定位方法、装置及设备
CN111965671A (zh) * 2020-09-28 2020-11-20 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种gnss信号质量监测与干扰监测定位系统及方法
CN117092665A (zh) * 2023-08-03 2023-11-21 广州海格晶维信息产业有限公司 一种组合导航设备抗多径干扰的方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1530635A (zh) * 2003-03-14 2004-09-22 ���ǵ�����ʽ���� 导航系统中检测可移动物体位置的设备及其方法
US20080082265A1 (en) * 2006-09-15 2008-04-03 Sin Etke Technology Co., Ltd. Accurate positioning system for a vehicle and its positioning method
CN101334288A (zh) * 2008-08-07 2008-12-31 北京工业大学 基于标准线路匹配的公交准确定位方法
CN101688910A (zh) * 2007-05-24 2010-03-31 电子地图有限公司 使用绝对定位系统及相对定位系统来确定位置的定位装置及方法、计算机程序及数据载体
CN105424041A (zh) * 2016-01-18 2016-03-23 重庆邮电大学 一种基于bd/ins紧耦合的行人定位算法
CN106710204A (zh) * 2017-01-13 2017-05-24 浙江大学 公交车辆的交叉定位装置和方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1530635A (zh) * 2003-03-14 2004-09-22 ���ǵ�����ʽ���� 导航系统中检测可移动物体位置的设备及其方法
US20080082265A1 (en) * 2006-09-15 2008-04-03 Sin Etke Technology Co., Ltd. Accurate positioning system for a vehicle and its positioning method
CN101688910A (zh) * 2007-05-24 2010-03-31 电子地图有限公司 使用绝对定位系统及相对定位系统来确定位置的定位装置及方法、计算机程序及数据载体
CN101334288A (zh) * 2008-08-07 2008-12-31 北京工业大学 基于标准线路匹配的公交准确定位方法
CN105424041A (zh) * 2016-01-18 2016-03-23 重庆邮电大学 一种基于bd/ins紧耦合的行人定位算法
CN106710204A (zh) * 2017-01-13 2017-05-24 浙江大学 公交车辆的交叉定位装置和方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张再勇: "车载GPS_MIMU_DM组合导航系统研究", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》 *
梅浩等: "城市道路GNSS脆弱性评估技术", 《哈尔滨工业大学学报》 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108051839A (zh) * 2017-10-27 2018-05-18 成都天合世纪科技有限责任公司 一种车载三维定位装置及三维定位的方法
CN108051839B (zh) * 2017-10-27 2021-11-05 成都天合世纪科技有限责任公司 一种车载三维定位装置及三维定位的方法
CN108846598A (zh) * 2018-03-29 2018-11-20 宏图物流股份有限公司 一种车辆定位的方法及装置
CN108931260A (zh) * 2018-07-04 2018-12-04 山东省科学院自动化研究所 Gps欺骗时无人驾驶车辆定位安全测试系统及方法
CN108931260B (zh) * 2018-07-04 2022-06-07 山东省科学院自动化研究所 Gps欺骗时无人驾驶车辆定位安全测试系统及方法
CN108920892A (zh) * 2018-09-30 2018-11-30 南京地铁集团有限公司 一种城市轨道交通站点脆弱性测量方法
CN109631886A (zh) * 2018-12-28 2019-04-16 江苏满运软件科技有限公司 车辆定位方法、装置、电子设备、存储介质
CN109631886B (zh) * 2018-12-28 2021-07-27 江苏满运物流信息有限公司 车辆定位方法、装置、电子设备、存储介质
CN110795518A (zh) * 2019-10-24 2020-02-14 中国人民解放军63653部队 一种基于vc++的gnss定位导航软件系统
CN110795518B (zh) * 2019-10-24 2022-10-11 中国人民解放军63653部队 一种基于vc++的gnss定位导航软件系统
CN111323802A (zh) * 2020-03-20 2020-06-23 北京百度网讯科技有限公司 车辆定位方法、装置及设备
CN111323802B (zh) * 2020-03-20 2023-02-28 阿波罗智能技术(北京)有限公司 智能驾驶车辆定位方法、装置及设备
CN111965671A (zh) * 2020-09-28 2020-11-20 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种gnss信号质量监测与干扰监测定位系统及方法
CN117092665A (zh) * 2023-08-03 2023-11-21 广州海格晶维信息产业有限公司 一种组合导航设备抗多径干扰的方法及系统
CN117092665B (zh) * 2023-08-03 2024-04-19 广州海格晶维信息产业有限公司 一种组合导航设备抗多径干扰的方法及系统

Also Published As

Publication number Publication date
CN107247275B (zh) 2020-12-08

Similar Documents

Publication Publication Date Title
CN107247275A (zh) 基于公交车的城市gnss脆弱性监测系统及其方法
CN108196289B (zh) 一种卫星信号受限条件下的列车组合定位方法
CN101360971B (zh) 定位装置及导航系统
CN103453913B (zh) 停车场中车辆的组合定位方法
CN104061899B (zh) 一种基于卡尔曼滤波的车辆侧倾角与俯仰角估计方法
JP4090852B2 (ja) Gps測位による列車走行情報検出装置及びその列車走行情報検出方法
CN109059909A (zh) 基于神经网络辅助的卫星/惯导列车定位方法与系统
CN101476894A (zh) 车载sins/gps组合导航系统性能增强方法
CN101907714A (zh) 基于多传感器数据融合的gps辅助定位系统及其定位方法
CN104713555A (zh) 应用全天域中性点辅助定向的车辆自主导航方法
CN104990554B (zh) Gnss盲区中基于vanet车辆间协作的惯性导航定位方法
CN101900573B (zh) 一种实现陆用惯性导航系统运动对准的方法
CN109813306A (zh) 一种无人车规划轨迹卫星定位数据可信度计算方法
Li et al. Constrained unscented Kalman filter based fusion of GPS/INS/digital map for vehicle localization
CN106093992A (zh) 一种基于cors的亚米级组合定位导航系统及导航方法
CN102538790A (zh) 惯性导航中陀螺仪参数的差异性解决方法
CN113758483A (zh) 一种自适应fkf地图匹配方法及系统
Iqbal et al. Experimental results on an integrated GPS and multisensor system for land vehicle positioning
Jiménez et al. Limitations of positioning systems for developing digital maps and locating vehicles according to the specifications of future driver assistance systems
Wang et al. Performance analysis of GNSS/MIMU tight fusion positioning model with complex scene feature constraints
CN113074757B (zh) 车载惯导安装误差角的标定方法
Wang et al. Land vehicle navigation using odometry/INS/vision integrated system
CN104422459A (zh) 一种确定车辆位置的方法及装置、导航方法和导航系统
Song et al. Enhanced map-aided gps/3d riss combined positioning strategy in urban canyons
Moussa et al. Investigation of different low-cost land vehicle navigation systems based on CPD sensors and vehicle information

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant