CN110333526B - 一种基于gps的仿真航迹推算传感器方法 - Google Patents

一种基于gps的仿真航迹推算传感器方法 Download PDF

Info

Publication number
CN110333526B
CN110333526B CN201910438926.1A CN201910438926A CN110333526B CN 110333526 B CN110333526 B CN 110333526B CN 201910438926 A CN201910438926 A CN 201910438926A CN 110333526 B CN110333526 B CN 110333526B
Authority
CN
China
Prior art keywords
gps
time
psi
reckoning sensor
longitude
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.)
Expired - Fee Related
Application number
CN201910438926.1A
Other languages
English (en)
Other versions
CN110333526A (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.)
University of Electronic Science and Technology of China Zhongshan Institute
Original Assignee
University of Electronic Science and Technology of China Zhongshan Institute
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 University of Electronic Science and Technology of China Zhongshan Institute filed Critical University of Electronic Science and Technology of China Zhongshan Institute
Priority to CN201910438926.1A priority Critical patent/CN110333526B/zh
Publication of CN110333526A publication Critical patent/CN110333526A/zh
Application granted granted Critical
Publication of CN110333526B publication Critical patent/CN110333526B/zh
Expired - Fee Related 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/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/52Determining velocity

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)

Abstract

本发明涉及一种基于GPS的仿真航迹推算传感器方法,根据GPS提供的测量仿真实现航迹推算传感器,其特征在于在GPS信号不能覆盖的区域利用GPS采集的历史数据估计线速度和角速度。

Description

一种基于GPS的仿真航迹推算传感器方法
技术领域
本发明涉及一种基于GPS的仿真航迹推算传感器方法,根据GPS提供的测量仿真实现航迹推算传感器,其特征在于在GPS信号不能覆盖的区域利用GPS采集的历史数据估计线速度和角速度。
背景技术
GPS提供基于通行卫星的精确导航和定位信息,广泛应用于汽车、飞机、船舶、手机等设备,在智能交通、智慧城市、机器人导航等领域被广泛使用。
航迹推算是根据本体的运行线速度和角速度来估计本体自身位姿的过程,通常利用编码器和陀螺仪等传感器来测量线速度和角速度。当没有这类传感器时,如何利用GPS信息来仿真实现编码器和陀螺仪等航迹推算传感器的功能,具有非常重要的实际价值。
在现有技术中,李昇智等人提出了一种基于GPS 轨迹数据的混合多步Markov 位置预测方法(参见,李昇智,乔建忠,林树宽,杨迪;《基于GPS 轨迹数据的混合多步Markov位置预测》,东北大学学报(自然科学版),2017年第12期1686-1690页)。该方法提出了基于混合多步Markov模型的位置预测方法,在将原始GPS 轨迹转化为区域轨迹的基础上,对各多步模型进行融合,利用Adaboost算法确定多步模型影响系数。
其不足在于没有利用GPS测量提供的精度信息。
发明内容
本发明的目的在于克服上述技术存在的缺陷,提供一种基于GPS的仿真航迹推算传感器方法,根据GPS提供的测量仿真实现航迹推算传感器,其特征在于利用了GPS测量精度信息,在GPS信号不能覆盖的区域利用GPS采集的历史数据估计线速度和角速度。
其具体技术方案包括以下步骤:
输入:GPS测量数据
Figure 923191DEST_PATH_IMAGE001
,
Figure 146362DEST_PATH_IMAGE002
,
Figure 617795DEST_PATH_IMAGE003
Figure 406759DEST_PATH_IMAGE004
,N;
其中,
Figure 39866DEST_PATH_IMAGE001
表示t时刻的纬度(单位为弧度),
Figure 750333DEST_PATH_IMAGE002
表示t时刻的经度(单位为弧度),
Figure 494298DEST_PATH_IMAGE005
表示t时刻的GPS测量精度(误差范围半径,单位为米),
Figure 341031DEST_PATH_IMAGE003
表示t时刻的时间(单位为 秒),
Figure 207356DEST_PATH_IMAGE004
为设定的精度阈值,默认值为60m,N为设定的循环次数;
输出:
Figure 641005DEST_PATH_IMAGE006
Figure 454240DEST_PATH_IMAGE007
,其中
Figure 889901DEST_PATH_IMAGE006
表示t时刻的线速度(m/s),
Figure 130389DEST_PATH_IMAGE007
表示t时刻的偏航率(rad/ s);
步骤1:r=6371,
Figure 346607DEST_PATH_IMAGE008
Figure 697954DEST_PATH_IMAGE009
分别设置为GPS最初时刻测量的纬度和经度,t=0,
Figure 253700DEST_PATH_IMAGE010
Figure 665090DEST_PATH_IMAGE011
步骤2:t=t+1;
步骤3:获取时间步t的纬度
Figure 368604DEST_PATH_IMAGE001
,经度
Figure 523641DEST_PATH_IMAGE002
,精度
Figure 432429DEST_PATH_IMAGE005
和时间
Figure 14720DEST_PATH_IMAGE003
步骤4:如果
Figure 408793DEST_PATH_IMAGE005
小于
Figure 101942DEST_PATH_IMAGE004
Figure 897860DEST_PATH_IMAGE004
为设定的精度阈值,缺省值为60m,转步骤;否者转步 骤5;
步骤5:
Figure 651052DEST_PATH_IMAGE012
步骤6:
Figure 532420DEST_PATH_IMAGE013
步骤7:
Figure 763682DEST_PATH_IMAGE014
步骤8:
Figure 903852DEST_PATH_IMAGE015
步骤9:
Figure 624683DEST_PATH_IMAGE016
步骤10:
Figure 727768DEST_PATH_IMAGE017
,其中
Figure 497141DEST_PATH_IMAGE018
Figure 267651DEST_PATH_IMAGE019
步骤11:
Figure 159384DEST_PATH_IMAGE020
步骤12:如果
Figure 749765DEST_PATH_IMAGE021
,则
Figure 322829DEST_PATH_IMAGE022
,如果
Figure 947845DEST_PATH_IMAGE023
,则
Figure 213742DEST_PATH_IMAGE024
步骤13:
Figure 789954DEST_PATH_IMAGE025
步骤14:如果t<N,转步骤2;如果t>=N,结束。
与现有技术相比,本发明的有益效果为:在GPS信号不能覆盖的区域也能利用GPS采集的历史数据估计线速度和角速度。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合附图和具体实例进一步阐述本发明。
一种基于GPS的仿真航迹推算传感器方法,根据GPS提供的测量仿真实现航迹推算传感器,其特征在于在GPS信号不能覆盖的区域利用GPS采集的历史数据估计线速度和角速度,包括以下步骤:
输入:GPS测量数据
Figure 166709DEST_PATH_IMAGE001
,
Figure 646232DEST_PATH_IMAGE002
,
Figure 879767DEST_PATH_IMAGE003
Figure 444741DEST_PATH_IMAGE004
,N;
其中,
Figure 828449DEST_PATH_IMAGE001
表示t时刻的纬度(单位为弧度),
Figure 162478DEST_PATH_IMAGE002
表示t时刻的经度(单位为弧度),
Figure 566914DEST_PATH_IMAGE005
表示t时刻的GPS测量精度(误差范围半径,单位为米),
Figure 619184DEST_PATH_IMAGE003
表示t时刻的时间(单位为 秒),
Figure 104785DEST_PATH_IMAGE004
为设定的精度阈值,默认值为60m,N为设定的循环次数;
输出:
Figure 293321DEST_PATH_IMAGE006
Figure 806342DEST_PATH_IMAGE007
,其中
Figure 80329DEST_PATH_IMAGE006
表示t时刻的线速度(m/s),
Figure 602577DEST_PATH_IMAGE007
表示t时刻的偏航率(rad/ s);
步骤1:r=6371,
Figure 911199DEST_PATH_IMAGE008
Figure 595121DEST_PATH_IMAGE009
分别设置为GPS最初时刻测量的纬度和经度,t=0,
Figure 120518DEST_PATH_IMAGE010
Figure 977615DEST_PATH_IMAGE011
步骤2:t=t+1;
步骤3:获取时间步t的纬度
Figure 875164DEST_PATH_IMAGE001
,经度
Figure 261146DEST_PATH_IMAGE002
,精度
Figure 509725DEST_PATH_IMAGE005
和时间
Figure 170513DEST_PATH_IMAGE003
步骤4:如果
Figure 188148DEST_PATH_IMAGE005
小于
Figure 479452DEST_PATH_IMAGE004
Figure 12065DEST_PATH_IMAGE004
为设定的精度阈值,缺省值为60m,转步骤;否者转步 骤5;
步骤5:
Figure 679806DEST_PATH_IMAGE012
步骤6:
Figure 551947DEST_PATH_IMAGE013
步骤7:
Figure 515617DEST_PATH_IMAGE014
步骤8:
Figure 535526DEST_PATH_IMAGE015
步骤9:
Figure 741379DEST_PATH_IMAGE016
步骤10:
Figure 468027DEST_PATH_IMAGE017
,其中
Figure 366713DEST_PATH_IMAGE018
Figure 873918DEST_PATH_IMAGE019
步骤11:
Figure 617883DEST_PATH_IMAGE020
步骤12:如果
Figure 464616DEST_PATH_IMAGE021
,则
Figure 268624DEST_PATH_IMAGE022
,如果
Figure 263125DEST_PATH_IMAGE023
,则
Figure 810781DEST_PATH_IMAGE024
步骤13:
Figure 10556DEST_PATH_IMAGE025
步骤14:如果t<N,转步骤2;如果t>=N,结束。
以上所述,仅为本发明最佳实施方式,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可显而易见地得到的技术方案的简单变化或等效替换均落入本发明的保护范围内。

Claims (1)

1.一种基于GPS的仿真航迹推算传感器方法,根据GPS提供的测量仿真实现航迹推算传感器,其特征在于利用GPS测量精度信息,在GPS信号不能覆盖的区域利用GPS采集的历史数据估计线速度和角速度,包括以下步骤:
输入:GPS测量数据
Figure FDA0003215430470000011
λt,τt,Ψ,N,ωt
其中,
Figure FDA0003215430470000012
表示t时刻的纬度,单位为弧度,λt表示t时刻的经度,单位为弧度,ωt表示t时刻的GPS测量精度,误差范围半径,单位为米,τt表示t时刻的时间,单位为秒,Ψ为设定的精度阈值,默认值为60m,N为设定的循环次数;
输出:vt
Figure FDA0003215430470000013
其中vt表示t时刻的线速度,单位为m/s,
Figure FDA0003215430470000014
表示t时刻的偏航率,单位为rad/s;
步骤1:r=6371,
Figure FDA0003215430470000015
λ0分别设置为GPS最初时刻测量的纬度和经度,t=0,v0=0,
Figure FDA0003215430470000016
步骤2:t=t+1;
步骤3:获取时刻t的纬度
Figure FDA0003215430470000017
经度λt,精度ωt和时间τt
步骤4:如果ωt小于Ψ,Ψ为设定的精度阈值,缺省值为60m,转步骤8;否则转步骤5;
步骤5:δ=vt-1*(τtt-1)*0.001/r;
步骤6:
Figure FDA0003215430470000018
步骤7:
Figure FDA0003215430470000019
步骤8:
Figure FDA00032154304700000110
步骤9:
Figure FDA00032154304700000111
步骤10:bt=2*atan2(yt,xt),其中
Figure FDA00032154304700000112
Figure FDA00032154304700000113
步骤11:vt=dt/(τtt-1);
步骤12:如果bt-bt-1>π,则
Figure FDA00032154304700000114
如果bt-bt-1≤-π,则
Figure FDA00032154304700000115
步骤13:
Figure FDA00032154304700000116
步骤14:如果t<N,转步骤2;如果t≥N,结束。
CN201910438926.1A 2019-05-22 2019-05-22 一种基于gps的仿真航迹推算传感器方法 Expired - Fee Related CN110333526B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910438926.1A CN110333526B (zh) 2019-05-22 2019-05-22 一种基于gps的仿真航迹推算传感器方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910438926.1A CN110333526B (zh) 2019-05-22 2019-05-22 一种基于gps的仿真航迹推算传感器方法

Publications (2)

Publication Number Publication Date
CN110333526A CN110333526A (zh) 2019-10-15
CN110333526B true CN110333526B (zh) 2021-11-02

Family

ID=68139476

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910438926.1A Expired - Fee Related CN110333526B (zh) 2019-05-22 2019-05-22 一种基于gps的仿真航迹推算传感器方法

Country Status (1)

Country Link
CN (1) CN110333526B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012163509A (ja) * 2011-02-09 2012-08-30 Sumitomo Rubber Ind Ltd 車両用推定航法装置、車両用推定航法及び車両用推定航法のプログラム
CN202614240U (zh) * 2011-08-26 2012-12-19 康讯科技股份有限公司 结合车辆诊断系统的航迹推算装置
CN103344963A (zh) * 2013-07-12 2013-10-09 电子科技大学中山学院 一种基于激光雷达的鲁棒航迹推算方法
CN109765569A (zh) * 2017-11-09 2019-05-17 电子科技大学中山学院 一种基于激光雷达实现虚拟航迹推算传感器的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012163509A (ja) * 2011-02-09 2012-08-30 Sumitomo Rubber Ind Ltd 車両用推定航法装置、車両用推定航法及び車両用推定航法のプログラム
CN202614240U (zh) * 2011-08-26 2012-12-19 康讯科技股份有限公司 结合车辆诊断系统的航迹推算装置
CN103344963A (zh) * 2013-07-12 2013-10-09 电子科技大学中山学院 一种基于激光雷达的鲁棒航迹推算方法
CN109765569A (zh) * 2017-11-09 2019-05-17 电子科技大学中山学院 一种基于激光雷达实现虚拟航迹推算传感器的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Automotive Dead-Reckoning Navigation System Based on Vehicle Speed and Yaw Rate";Nicolae Dicu et al.;《SACI 2018.IEEE 12th International Symposium on Applied Computational Intelligence and Informatics》;20180519;第000225-000228页 *
"Positioning of vehicle on undulating ground using GPS and dead reckoning";Toshihiro Aono et al.;《Proceedings of the I998 IEEE International Conference on Robotics & Automation》;19980531;第3443-3448页 *

Also Published As

Publication number Publication date
CN110333526A (zh) 2019-10-15

Similar Documents

Publication Publication Date Title
CN104061899B (zh) 一种基于卡尔曼滤波的车辆侧倾角与俯仰角估计方法
CN107990891B (zh) 基于长基线和信标在线标定的水下机器人组合导航方法
CN104165641B (zh) 一种基于捷联惯导/激光测速仪组合导航系统的里程计标定方法
CN106595715B (zh) 基于捷联惯导与卫星组合导航系统里程计标定方法及装置
CN109343095B (zh) 一种车载导航车辆组合定位装置及其组合定位方法
CN104729506A (zh) 一种视觉信息辅助的无人机自主导航定位方法
CN104197958B (zh) 一种基于激光测速仪航位推算系统的里程计标定方法
CN109507706B (zh) 一种gps信号丢失的预测定位方法
CN108387236B (zh) 一种基于扩展卡尔曼滤波的偏振光slam方法
CN106767894A (zh) 一种用于捷联惯导的北斗/里程计组合标定方法
US20200182627A1 (en) Segmented path coordinate system
CN103604430A (zh) 一种基于边缘化ckf重力辅助导航的方法
CN102608643A (zh) 一种组合式车辆位置测量方法
WO2018214226A1 (zh) 一种无人车实时姿态测量方法
Gao et al. An integrated land vehicle navigation system based on context awareness
CN110333526B (zh) 一种基于gps的仿真航迹推算传感器方法
CN103245793B (zh) 一种组合导航洋流速度测量方法
CN108801238B (zh) 一种基于gpx的路径优化方法
Farroni et al. A comparison among different methods to estimate vehicle sideslip angle
CN102865880B (zh) 一种预计水下机器人导航系统圆概率偏差的方法
CN115046547A (zh) 一种复杂城市环境下的自动驾驶汽车定位系统及方法
CN112902950B (zh) 一种面向低速运动载体中mems级imu的初始对准方法
CN101793529B (zh) 一种双伪卫星辅助惯性导航系统位置校准方法
Liu et al. A tightly-coupled GNSS/INS/MM integrated system based on binary search algorithm for train localization applications
CN103389090A (zh) 一种惯性导航系统极区导航模式初始速度的测量方法

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20211102

CF01 Termination of patent right due to non-payment of annual fee