CN106950586A - 用于农机作业的gnss/ins/车辆组合导航方法 - Google Patents
用于农机作业的gnss/ins/车辆组合导航方法 Download PDFInfo
- Publication number
- CN106950586A CN106950586A CN201710046370.2A CN201710046370A CN106950586A CN 106950586 A CN106950586 A CN 106950586A CN 201710046370 A CN201710046370 A CN 201710046370A CN 106950586 A CN106950586 A CN 106950586A
- Authority
- CN
- China
- Prior art keywords
- gnss
- information
- speed
- ins
- attitude
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 239000011159 matrix material Substances 0.000 claims abstract description 111
- 238000001914 filtration Methods 0.000 claims abstract description 56
- 238000004364 calculation method Methods 0.000 claims abstract description 30
- 238000005259 measurement Methods 0.000 claims abstract description 16
- 230000007704 transition Effects 0.000 claims description 53
- 238000012937 correction Methods 0.000 claims description 36
- 239000013598 vector Substances 0.000 claims description 12
- 239000013078 crystal Substances 0.000 claims description 9
- 230000001360 synchronised effect Effects 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000013213 extrapolation Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000033001 locomotion Effects 0.000 claims description 3
- 230000036544 posture Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 7
- 239000003607 modifier Substances 0.000 abstract 2
- 230000003068 static effect Effects 0.000 description 11
- 238000013461 design Methods 0.000 description 7
- 230000001133 acceleration Effects 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000003306 harvesting Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining 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/42—Determining position
- G01S19/48—Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system
- G01S19/49—Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system whereby the further system is an inertial position system, e.g. loosely-coupled
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明提供了导航技术领域内的一种用于农机作业的GNSS/INS/车辆组合导航方法,步骤1:对惯性测量单元进行零偏修正后进行惯导解算,经卡尔曼滤波计算状态转移矩阵和系统协方差矩阵;步骤2:判断GNSS信息是否已更新,若GNSS信息更新,GNSS信息外推,惯导与GNSS的位置速度差值为观测量进行滤波估计,判断载体的动态,反馈修正平台误差和IMU零偏误差返回至步骤1中循环执行;否则,当前信息为载体的位置速度姿态信息;步骤3:判断车辆信息是否已更新,若车辆信息更新,车辆信息外推,惯导与车辆的速度差值为观测量进行滤波估计,反馈修正平台误差和IMU零偏误差;否则,当前信息为载体的位置速度姿态信息;本发明导航精度高,低成本IMU可以获得较高导航精度。
Description
技术领域
本发明涉及一种自动掉头路径规划及其控制方法,特别涉及一种用于农机作业的GNSS/INS/车辆组合导航方法。
背景技术
随着农业自动化的发展、国家惠农政策的普及以及农民老龄化的比例上升,精准农业逐渐成为国际研究热点,精准农业的普及可以解放双手,提高工作效率、提供土地利用率、及时播种收割、节约种子等投入,达到减少劳动力、减少投入、增加产量的目标。
在现有的导航系统中,最常用的技术为基于GNSS/INS 组合导航方法设计的卫星组合导航系统,该系统包括惯性测量单元和GNSS接收单元,导航方法具体的为,惯性测量单元测量载体的加速度和角速度,通过INS计算单元解算得到速度、位置和姿态信息,GNSS输出位置和速度信息,对GNSS输出的位置速度信息和INS计算单元输出的位置速度姿态信息进行卡尔曼滤波,估计出INS 的位置、速度和姿态误差,以及惯性测量单元的误差,将这些误差反馈给INS计算单元,修正INS的导航解和传感器误差,从而抑制INS误差随时间的增长,获得稳定的导航解,但是在信号被遮挡等复杂环境下时,GNSS信号会瞬间丢失,此时IMU发散,导航精度差,同时,这样的导航系统中需要提供高精度的IMU才能获得较高精度的姿态信息,对传感器自身的要求高,成本高。
发明内容
针对现有技术中的缺陷,本发明的目的在于克服上述现有技术中的不足之处,解决现有技术中信号遮挡等复杂环境下导航精度差的技术问题,提供一种用于农机作业的GNSS/INS/车辆组合导航方法,本发明的导航精度高,对IMU自身的要求低,成本低。
本发明的目的是这样实现的:一种用于农机作业的GNSS/INS/车辆组合导航方法,具体包括以下步骤,
步骤1(INS解算):对惯性测量单元进行零偏修正后,通过INS计算单元进行惯导解算,经过卡尔曼滤波计算出状态转移矩阵,根据状态转移矩阵进行卡尔曼滤波一步预测得到系统协方差矩阵;
步骤2(GNSS与INS的组合导航):判断GNSS信息是否已经更新,若GNSS信息更新,GNSS信息外推与IMU(惯导测量单元)时间同步,惯导与GNSS的位置、速度差值作为观测量进行kalman滤波估计,kalman滤波估计更新步骤1中的状态量和状态协方差矩阵,构成全维状态观测量,判断载体的动静状态,反馈修正平台误差和IMU零偏误差返回至步骤1中循环执行;否则,当前信息就为载体的位置、速度和姿态信息;
步骤3(车辆与INS的组合导航):判断车辆信息是否已经更新,若车辆信息更新,车辆信息外推与IMU更新时间同步,惯导与车辆的速度差值作为观测量进行kalman滤波估计,得到状态误差修正量,用于反馈修正平台误差和IMU零偏误差;否则,当前信息就为载体的位置、速度和姿态信息;
其中,修正平台误差指位置、速度和姿态误差,IMU零偏误差指惯性测量单元的零偏误差;惯性测量单元包括加速度计和陀螺仪。
本发明工作时,当惯导解算得到新的位置、速度、姿态信息,经过卡尔曼滤波计算出状态转移矩阵,根据状态转移矩阵进行卡尔曼滤波一步预测得到状态协方差矩阵,此时,若GNSS信息得到新的位置和速度信息时,惯导与GNSS的位置、速度差值作为观测量进行kalman滤波估计,得到全维状态量的误差修正值,判断载体动静态,反馈修正平台误差和IMU零偏误差至INS解算步骤中的INS计算单元,同时,输出经过kalman滤波估计后解算出来的载体速度、位置和姿态信息;若GNSS信息没有更新,则当前信息就是载体的位置、速度和姿态信息;当惯导解算得到新的位置、速度和姿态信息,且里程计更新得到新的车辆信息,则惯导与车辆的速度差值作为观测量进行kalman滤波估计,得到状态误差修正量,反馈修正平台误差和IMU零偏误差,反馈IMU零偏误差至IMU的零偏修正,修正平台误差反馈至INS计算单元重新进行惯导解算以循环执行,同时,输出当前载体速度、位置和姿态信息;若车辆信息的时间没有更新,则直接输出当前载体速度、位置和姿态信息;本发明提供了GNSS和INS组合导航以及INS和车辆组合导航方法,GNSS/INS组合导航方法,动态初始对准收敛速度快、精度高,当GNSS信号丢失时,该INS和车辆组合导航技术可以很好地抑制IMU发散,提供高精度的位置、速度和姿态信息;可应用于农机在田间作业时自动驾驶的导航工作中。
作为本发明的进一步改进,所述步骤1中,INS计算单元进行惯导解算,解算过程如下,
a)速度解算:
(1-1)
其中,
为k时刻导航坐标系下的载体速度;
为单位矩阵;
为导航系相对于惯性系在k-1到k时刻的角度变化量;
表示叉乘;
为k-1时刻姿态旋转矩阵;
为k-1到k时刻机体系下载体速度增量;
为划桨效应速度补偿量;
为旋转效应速度补偿量;
为由重力和科里奥利引起的速度增量;
初始速度由GNSS提供;
b)位置解算:
位置四元数解算:
(1-2)
其中:
(1-3)
为导航系k-1到k时刻等效旋转矢量;
为k-1到k时刻地球自转引起的等效旋转矢量;
为k时刻位置四元数;
利用更新后的四元数直接转换得到更新后的位置,四元数与位置之间关系为:
(1-4)
其中,
L为当前位置纬度;
λ为当前位置经度;
INS的初始位置由GNSS提供;
c)姿态解算:
姿态四元数解算:
(1-5)
其中:
(1-6)
为导航系k-1到k时刻等效旋转矢量;
为机体系下等效旋转矢量;
为k时刻姿态四元数;
利用更新后的姿态四元数直接转换得到姿态矩阵:
(1-7)
其中,
为姿态矩阵方向余弦形式;
为姿态四元数的第x维;
利用姿态矩阵与方向角之间的关系可以得到姿态角:
(1-8)
(1-9)
其中,
r、p、h分别为翻转角、俯仰角、航向角;
为姿态矩阵的第x/y维;
初始航向由GNSS东向和北向速度解算,初始翻转角和俯仰角设定为0,根据上述姿态角转换矩阵初始化初始姿态阵;
此设计中,通过惯导解算计算出载体的速度、位置和姿态信息。
作为本发明的进一步改进,所述步骤1中,惯导解算后的经过卡尔曼滤波计算出状态转移矩阵和一步预测,具体的为,
d)状态转移矩阵的计算:
选取位置误差、速度误差、平台姿态角以及陀螺和加表的零偏误差共15维作为滤波估计的状态量,即
(2-1)
其中,为陀螺仪的三个轴向漂移,为加表三个轴向的偏置误差,并不是对应IMU的XYZ三轴零偏误差,δL、δλ和δh为INS的三个位置误差由地球坐标系描述,δVN、δVE和δVD为东北天三个方向上的速度误差;
根据所选状态量之间的关系可以得到状态转移矩阵如下,
(2-2)
其中,
F1为位置与位置之间状态转移矩阵;
F2为位置与速度之间状态转移矩阵;
F3为速度与位置之间状态转移矩阵;
F4为速度与速度之间状态转移矩阵;
F5为速度与姿态之间状态转移矩阵;
F6为速度与加表零偏误差之间状态转移矩阵;
F7为姿态与姿态之间状态转移矩阵;
F8为姿态与陀螺零偏误差之间状态转移矩阵;
F9为陀螺零偏误差与陀螺零偏误差之间状态转移矩阵;
F10为加表零偏误差与加表零偏误差之间状态转移矩阵;
e)一步预测:
一步预测n-1时刻到n时刻状态量状态协方差矩阵,
(2-3)
(2-4)
其中,
为n-1时刻状态量X的估计值;
为n-1时刻状态协方差矩阵;
为n-1时刻到n时刻状态转移矩阵;
Q为系统噪声矩阵;
为IMU更新周期;
此设计中,通过以上步骤解算出状态转移矩阵和状态协方差矩阵。
为了获得高精度的同步时间,所述步骤2中,实现GNSS信息与IMU信息同步更新的方法为,GNSS的PPS(脉冲)信号的周期为1s,每个脉冲处系统晶振均重新计时,相当于每1s校准一次晶振,当GNSS信息更新时,将GNSS信息外推到最近时刻的IMU处;所述车辆与INS的组合导航中,实现里程计信息与IMU信息同步更新的方法与GNSS信息与IMU信息同步更新的方法相同;此设计中,通过软校准的时间同步方法,实现GNSS信息与INS信息的同步或者INS信息与里程计信息的同步,实现简单,时间同步精度高,即使使用低精度的晶振也可以获得高精度的同步时间。
为了进一步提高载体静止时载体姿态信息的精度,所述步骤2中,判断载体的动静状态, 反馈修正平台误差和IMU零偏误差具体的为,若载体静止,将GNSS输出的速度约束为零并与IMU测得的当前速度作差,将其差值作为kalman滤波估计的观测量,同时,约束载体航向,经过kalman滤波估计得到状态修正量,反馈经过kalman滤波估计的修正平台误差和IMU零偏误差至INS计算单元进行解算,同时,输出当前载体的速度、位置和姿态信息;若载体为动态,则直接反馈修正平台误差和IMU零偏误差至INS计算单元循环执行,同时,输出当前载体的速度、位置和姿态信息;此设计中,零速修正和航向锁定技术,在载体静止时可以很好地抑制航向漂移,提供高精度的姿态信息。
为了进一步修正观测误差,所述Kalman滤波估计更新状态量和状态协方差矩阵,具体的为,
f)解算滤波增益:
(3-1)
其中,
为n时刻滤波增益;
为n时刻的观测矩阵,由观测量与状态量之间的关系得到;
为n时刻观测噪声协方差矩阵,由GNSS信息提供;
g)状态量的更新:
(3-2)
其中,
为n时刻的状态估计量;
为n时刻的观测量即IMU与GNSS的位置速度差值;
h)状态协方差矩阵的更新:
(3-3)
其中,
为n时刻状态协方差矩阵;
为单位矩阵;
此设计中,在滤波估计时加入了质量控制,即当实际观测误差量与估计的观测误差相差较大时,对Kalman滤波的观测误差进行进一步的修正。
为了进一步提高姿态角的精度,所述步骤3中,获取车辆信息时,具体的为,将农机运动模型简化为简单的四轮模型,以后轴中心为原点的前右下坐标系,将车辆后轴中心的右和下两个轴的速度约束为零,前向轴的速度由里程计测得;将车辆三轴的速度与INS的速度差值作为Kalman滤波估计的观测量,经过Kalman滤波估计后输出修正平台误差和IMU零偏误差并反馈至INS计算单元解算循环执行;此设计,可有效辅助航向信息,提供更高精度的姿态角,当GNSS信号丢失时,该技术可以很好地抑制IMU发散,提供更高精度的位置、速度和姿态信息。
附图说明
图1为本发明中导航方法的总体流程图。
图2为本发明中惯性导航系统的解算流程图。
图3为本发明中GNSS和INS组合导航技术流程。
图4为本发明中INS和车辆组合导航技术流程。
具体实施方式
下面结合附图对本发明作进一步的说明。
如图1~4所示的一种用于农机作业的GNSS/INS/车辆组合导航方法,具体包括以下步骤,
步骤1(INS解算):对惯性测量单元进行零偏修正后,通过INS计算单元进行惯导解算,经过卡尔曼滤波计算出状态转移矩阵,根据状态转移矩阵进行卡尔曼滤波一步预测得到状态协方差矩阵;
步骤2(GNSS与INS的组合导航):判断GNSS信息是否已经更新,若GNSS信息更新,GNSS信息外推与IMU(惯导测量单元)时间同步,惯导与GNSS的位置、速度差值作为观测量进行kalman滤波估计,kalman滤波估计更新步骤1中的状态量和状态协方差矩阵,构成全维状态观测量,判断载体的动静状态,反馈修正平台误差和IMU零偏误差返回至步骤1中循环执行;否则,当前信息就为载体的位置、速度和姿态信息;
步骤3(车辆与INS的组合导航):判断车辆信息是否已经更新,若车辆信息更新,车辆信息外推与IMU更新时间同步,惯导与车辆的速度差值作为观测量进行kalman滤波估计,得到状态误差修正量,用于反馈修正平台误差和IMU零偏误差;否则,当前信息就为载体的位置、速度和姿态信息;
其中,修正平台误差指位置、速度和姿态误差,IMU零偏误差指惯性测量单元的零偏误差;惯性测量单元包括加速度计和陀螺仪;
步骤1中,INS计算单元进行惯导解算,解算过程如下,
a)速度解算:
(1-1)
其中,
为k时刻导航坐标系下的载体速度;
为单位矩阵;
为导航系相对于惯性系在k-1到k时刻的角度变化量;
表示叉乘;
为k-1时刻姿态旋转矩阵;
为k-1到k时刻机体系下载体速度增量;
为划桨效应速度补偿量;
为旋转效应速度补偿量;
为由重力和科里奥利引起的速度增量;
初始速度由GNSS提供;
b)位置解算:
位置四元数解算:
(1-2)
其中:
(1-3)
为导航系k-1到k时刻等效旋转矢量;
为k-1到k时刻地球自转引起的等效旋转矢量;
为k时刻位置四元数;
利用更新后的四元数直接转换得到更新后的位置,四元数与位置之间关系为:
(1-4)
其中,
L为当前位置纬度;
λ为当前位置经度;
INS的初始位置由GNSS提供;
c)姿态解算:
姿态四元数解算:
(1-5)
其中:
(1-6)
为导航系k-1到k时刻等效旋转矢量;
为机体系下等效旋转矢量;
为k时刻姿态四元数;
利用更新后的姿态四元数直接转换得到姿态矩阵:
(1-7)
其中,
为姿态矩阵方向余弦形式;
为姿态四元数的第x维;
利用姿态矩阵与方向角之间的关系可以得到姿态角:
(1-8)
(1-9)
其中,
r、p、h分别为翻转角、俯仰角、航向角;
为姿态矩阵的第x/y维;
初始航向由GNSS东向和北向速度解算,初始翻转角和俯仰角设定为0,根据上述姿态角转换矩阵初始化初始姿态阵;
步骤1中,惯导解算后的经过卡尔曼滤波计算出状态转移矩阵和一步预测,具体的为,
d)状态转移矩阵的计算:
选取位置误差、速度误差、平台姿态角以及陀螺和加表的零偏误差共15维作为滤波估计的状态量,即
(2-1)
其中,为陀螺仪的三个轴向漂移,为加表三个轴向的偏置误差,并不是对应IMU的XYZ三轴零偏误差,δL、δλ和δh为INS的三个位置误差由地球坐标系描述,δVN、δVE和δVD为东北天三个方向上的速度误差;
根据所选状态量之间的关系可以得到状态转移矩阵如下,
(2-2)
其中,
F1为位置与位置之间状态转移矩阵;
F2为位置与速度之间状态转移矩阵;
F3为速度与位置之间状态转移矩阵;
F4为速度与速度之间状态转移矩阵;
F5为速度与姿态之间状态转移矩阵;
F6为速度与加表零偏误差之间状态转移矩阵;
F7为姿态与姿态之间状态转移矩阵;
F8为姿态与陀螺零偏误差之间状态转移矩阵;
F9为陀螺零偏误差与陀螺零偏误差之间状态转移矩阵;
F10为加表零偏误差与加表零偏误差之间状态转移矩阵;
e)一步预测:
一步预测n-1时刻到n时刻状态量状态协方差矩阵,
(2-3)
(2-4)
其中,
为n-1时刻状态量X的估计值;
为n-1时刻状态协方差矩阵;
为n-1时刻到n时刻状态转移矩阵;
Q为系统噪声矩阵;
为IMU更新周期;
GNSS信息相较于IMU信息更新较慢,当GNSS信息更新时,IMU不一定同步更新,所以要严格同步IMU与GNSS时间;实现GNSS信息与IMU信息同步更新的方法为,GNSS的PPS(脉冲)信号的周期为1s,每个脉冲处系统晶振均重新计时,相当于每1s校准一次晶振,当GNSS信息更新时,将GNSS信息外推到最近时刻的IMU处;车辆信息相较于IMU信息更新较慢,当车辆信息更新时,IMU不一定同步更新,所以要严格同步IMU与车辆信息更新的时间,实现车辆信息与IMU信息同步更新的方法与GNSS信息与IMU信息同步更新的方法相同;此设计中,通过软校准的时间同步方法,实现GNSS信息与INS信息的同步或者INS信息与里程计信息的同步,实现简单,时间同步精度高,即使使用低精度的晶振也可以获得高精度的同步时间。
为了进一步提高载体静止时载体姿态信息的精度,步骤3中,判断载体的动静状态, 反馈修正平台误差和IMU零偏误差具体的为,若载体静止,将GNSS输出的速度约束为零并与IMU测得的当前速度作差,将其差值作为kalman滤波估计的观测量,同时,约束载体航向,经过kalman滤波估计得到状态修正量,反馈经过kalman滤波估计的修正平台误差和IMU零偏误差至INS计算单元进行解算,同时,输出当前载体的速度、位置和姿态信息;若载体为动态,则直接反馈修正平台误差和IMU零偏误差至INS计算单元循环执行,同时,输出当前载体的速度、位置和姿态信息;
为了进一步修正观测误差, Kalman滤波估计更新状态量和状态协方差矩阵,具体的为,
f)解算滤波增益:
(3-1)
其中,
为n时刻滤波增益;
为n时刻的观测矩阵,由观测量与状态量之间的关系得到;
为n时刻观测噪声协方差矩阵,由GNSS信息提供;
g)状态量的更新:
(3-2)
其中,
为n时刻的状态估计量;
为n时刻的观测量即IMU与GNSS的位置速度差值;
h)状态协方差矩阵的更新:
(3-3)
其中,
为n时刻状态协方差矩阵;
为单位矩阵;
为了进一步提高姿态角的精度,步骤3中,将农机运动模型简化为简单的四轮模型,以后轴中心为原点的前右下坐标系,将车辆后轴中心的右和下两个轴的速度约束为零,前向轴的速度由里程计测得,将车辆三轴的速度与INS的速度差值作为Kalman滤波估计的观测量,经过Kalman滤波估计后输出修正平台误差和IMU零偏误差并反馈至INS计算单元解算循环执行。
本发明工作时,当惯导解算得到新的位置、速度、姿态信息,经过卡尔曼滤波计算出状态转移矩阵,根据状态转移矩阵进行卡尔曼滤波一步预测得到状态协方差矩阵,此时,若GNSS信息得到新的位置和速度信息时,惯导与GNSS的位置、速度差值作为观测量进行kalman滤波估计,得到15维状态量的误差修正值,判断载体动静态,若此时载体静止,则利用将GNSS输出的速度约束为零并与惯导元件测得的当前速度作差,将其差值作为观测量,进一步进行kalman滤波估计得到状态修正量,然后反馈修正平台误差和IMU零偏误差;若载体为动态,则直接反馈IMU零偏误差进行零偏修正,以及修正平台误差至步骤1中的INS计算单元进行位置、速度和姿态的解算,同时,输出当前经过kalman滤波估计后解算出来的载体速度、位置和姿态信息;若GNSS信息没有更新,则当前信息就是载体的位置、速度和姿态信息;当惯导解算得到新的位置、速度和姿态信息,且里程计更新得到新的车辆信息,则惯导与车辆的速度差值作为观测量进行kalman滤波估计,得到状态误差修正量,反馈修正平台误差和IMU零偏误差,反馈IMU零偏误差至IMU的零偏修正,修正平台误差反馈至步骤1中的INS计算单元重新进行惯导解算以循环执行,同时,输出当前载体速度、位置和姿态信息;若车辆信息没有更新,则直接输出当前载体速度、位置和姿态信息;本发明提供了GNSS和INS组合导航以及INS和车辆组合导航方法,GNSS/INS组合导航方法,动态初始对准收敛速度快、精度高,当GNSS信号丢失时,该INS和车辆组合导航技术可以很好地抑制IMU发散,提供高精度的位置、速度和姿态信息;采用的零速修正和航向锁定技术,在载体静止时可以很好地抑制航向漂移,提供高精度的姿态信息;使用软校准时间同步方法,实现简单,时间同步精度高,利用低精度的晶振可以获得高精度的同步时间;可应用于农机在田间作业时自动驾驶的导航工作中。
本发明并不局限于上述实施例,在本发明公开的技术方案的基础上,本领域的技术人员根据所公开的技术内容,不需要创造性的劳动就可以对其中的一些技术特征作出一些替换和变形,这些替换和变形均在本发明保护范围内。
Claims (7)
1.一种用于农机作业地形感知的GNSS/INS/车辆组合导航方法,其特征在于,
步骤1(INS解算):对惯性测量单元进行零偏修正后,通过INS计算单元进行惯导解算,经过卡尔曼滤波计算出状态转移矩阵,根据状态转移矩阵进行卡尔曼滤波一步预测得到系统协方差矩阵;
步骤2(GNSS与INS的组合导航):判断GNSS信息是否已经更新,若GNSS信息更新,GNSS信息外推与IMU(惯导测量单元)时间同步,惯导与GNSS的位置、速度差值作为观测量进行kalman滤波估计,kalman滤波估计更新步骤1中的状态量和状态协方差矩阵,构成全维状态观测量,判断载体的动态,反馈修正平台误差和IMU零偏误差返回至步骤1中循环执行;否则,当前信息就为载体的位置、速度和姿态信息;
步骤3(车辆与INS的组合导航):判断车辆信息是否已经更新,若车辆信息更新,车辆信息外推与IMU更新时间同步,惯导与车辆的速度差值作为观测量进行kalman滤波估计,得到状态误差修正量,用于反馈修正平台误差和IMU零偏误差;否则,当前信息就为载体的位置、速度和姿态信息;
其中,修正平台误差指位置、速度和姿态误差,IMU零偏误差指惯性测量单元的零偏误差;惯性测量单元包括加速度计和陀螺仪。
2.根据权利要求1所述的用于农机作业地形感知的GNSS/INS/车辆组合导航方法,其特征在于,所述步骤1,INS计算单元进行惯导解算,解算过程如下,
a)速度解算:
其中,
为k时刻导航坐标系下的载体速度;
I为单位矩阵;
ξk为导航系相对于惯性系在k-1到k时刻的角度变化量;
X表示叉乘;
为k-1时刻姿态旋转矩阵;
为k-1到k时刻机体系下载体速度增量;
Δvscul为划桨效应速度补偿量;
Δvcot为旋转效应速度补偿量;
为由重力和科里奥利引起的速度增量;
初始速度由GNSS提供;
b)位置解算:
位置四元数解算:
其中:
ζk为导航系k-1到k时刻等效旋转矢量;
ξk为k-1到k时刻地球自转引起的等效旋转矢量;
为k时刻位置四元数;
利用更新后的四元数直接转换得到更新后的位置,四元数与位置之间关系为:
其中,
L为当前位置纬度;
λ为当前位置经度;
INS的初始位置由GNSS提供;
c)姿态解算:
姿态四元数解算:
其中:
ζk为导航系k-1到k时刻等效旋转矢量;
φk为机体系下等效旋转矢量;
为k时刻姿态四元数;
利用更新后的姿态四元数直接转换得到姿态矩阵
其中,
为姿态矩阵方向余弦形式;
qx为姿态四元数的第x维;
利用姿态矩阵与方向角之间的关系可以得到姿态角:
其中,
r、p、h分别为翻转角、俯仰角、航向角;
为姿态矩阵的第x/y维;
初始航向由GNSS东向和北向速度解算,初始翻转角和俯仰角设定为0,根据上述姿态角转换矩阵初始化初始姿态阵。
3.根据权利要求1或2所述的用于农机作业地形感知的GNSS/INS/车辆组合导航方法,其特征在于,步骤1中,惯导解算后的经过卡尔曼滤波计算出状态转移矩阵和一步预测,具体的为,
d)状态转移矩阵的计算:
选取位置误差、速度误差、平台姿态角以及陀螺和加表的零偏误差共15维作为滤波估计的状态量,即
其中,εbfεbrεbd为陀螺仪的三个轴向漂移,▽bf▽br▽bd为加表三个轴向的偏置误差,并不是对应IMU的XYZ三轴零偏误差,δL、δλ和δh为INS的三个位置误差由地球坐标系描述,δVN、δVE和δVD为东北天三个方向上的速度误差;
根据所选状态量之间的关系可以得到状态转移矩阵如下:
其中,
F1为位置与位置之间状态转移矩阵;
F2为位置与速度之间状态转移矩阵;
F3为速度与位置之间状态转移矩阵;
F4为速度与速度之间状态转移矩阵;
F5为速度与姿态之间状态转移矩阵;
F6为速度与加表零偏误差之间状态转移矩阵;
F7为姿态与姿态之间状态转移矩阵;
F8为姿态与陀螺零偏误差之间状态转移矩阵;
F9为陀螺零偏误差与陀螺零偏误差之间状态转移矩阵;
F10为加表零偏误差与加表零偏误差之间状态转移矩阵;
e)一步预测:
一步预测n-1时刻到n时刻状态量状态协方差矩阵Pn|n-1,
Pn|n-1=φn|n-1Pn-1φn|n-1 T+(φn|n-1Q+Qφn|n-1 T)δt/2 (2-4)
其中,
为n-1时刻状态量X的估计值;
Pn-1为n-1时刻状态协方差矩阵;
φn|n-1为n-1时刻到n时刻状态转移矩阵;
Q为系统噪声矩阵;
δt为IMU更新周期。
4.根据权利要求3所述的用于农机作业地形感知的GNSS/INS/车辆组合导航方法,其特征在于,所述步骤2中,实现GNSS信息与IMU信息同步更新的方法为,GNSS的PPS(脉冲)信号的周期为1s,每个脉冲处系统晶振均重新计时,相当于每1s校准一次晶振,当GNSS信息更新时,将GNSS信息外推到最近时刻的IMU处;所述车辆与INS的组合导航中,实现里程计信息与IMU信息同步更新的方法与GNSS信息与IMU信息同步更新的方法相同。
5.根据权利要求4所述的用于农机作业地形感知的GNSS/INS/车辆组合导航方法,其特征在于,所述步骤2中,判断载体的动态,反馈修正平台误差和IMU零偏误差具体的为,若载体静止,将GNSS输出的速度约束为零并与IMU测得的当前速度作差,将其差值作为kalman滤波估计的观测量,同时,约束载体航向,经过kalman滤波估计得到状态修正量,反馈经过kalman滤波估计的修正平台误差和IMU零偏误差至INS计算单元进行解算,同时,输出当前载体的速度、位置和姿态信息;若载体为动态,则直接反馈修正平台误差和IMU零偏误差至INS计算单元循环执行,同时,输出当前载体的速度、位置和姿态信息。
6.根据权利要求5所述的用于农机作业地形感知的GNSS/INS/车辆组合导航方法,其特征在于,所述Kalman滤波估计更新状态量和状态协方差矩阵,具体的为,
f)解算滤波增益:
Kn=Pn|n-1Hn T(HnPn|n-1Hn T+Rn)-1 (3-1)
其中,
Kn为n时刻滤波增益;
Hn为n时刻的观测矩阵,由观测量与状态量之间的关系得到;
Rn为n时刻观测噪声协方差矩阵,由GNSS信息提供;
g)状态量的更新:
其中,
为n时刻的状态估计量;
Zn为n时刻的观测量即IMU与GNSS的位置速度差值;
h)状态协方差矩阵的更新:
Pn=(I-KnHn)Pn|n-1(I-KnHn)T+KnRnKn T (3-3)
其中,
Pn为n时刻状态协方差矩阵;
I为单位矩阵。
7.根据权利要求1、2或4~6任一项所述的用于农机作业地形感知的GNSS/INS/车辆组合导航方法,其特征在于,所述步骤3中,将农机运动模型简化为简单的四轮模型,以后轴中心为原点的前右下坐标系,将车辆后轴中心的右和下两个轴的速度约束为零,前向轴的速度由里程计测得,将车辆三轴的速度与INS的速度差值作为Kalman滤波估计的观测量,经过Kalman滤波估计后输出修正平台误差和IMU零偏误差并反馈至INS计算单元解算循环执行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710046370.2A CN106950586A (zh) | 2017-01-22 | 2017-01-22 | 用于农机作业的gnss/ins/车辆组合导航方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710046370.2A CN106950586A (zh) | 2017-01-22 | 2017-01-22 | 用于农机作业的gnss/ins/车辆组合导航方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106950586A true CN106950586A (zh) | 2017-07-14 |
Family
ID=59465555
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710046370.2A Pending CN106950586A (zh) | 2017-01-22 | 2017-01-22 | 用于农机作业的gnss/ins/车辆组合导航方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106950586A (zh) |
Cited By (41)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108036797A (zh) * | 2017-11-30 | 2018-05-15 | 深圳市隐湖科技有限公司 | 基于四轮独立驱动且结合imu的里程推算方法 |
CN108955851A (zh) * | 2018-07-12 | 2018-12-07 | 北京交通大学 | 利用ins和dtm确定gnss误差的方法 |
CN109059904A (zh) * | 2018-06-01 | 2018-12-21 | 浙江亚特电器有限公司 | 用于移动载具的组合导航方法 |
CN109143304A (zh) * | 2018-09-30 | 2019-01-04 | 百度在线网络技术(北京)有限公司 | 用于确定无人驾驶车辆位姿的方法和装置 |
CN109443347A (zh) * | 2017-07-31 | 2019-03-08 | 意法半导体股份有限公司 | 用于陆地车辆导航的系统以及对应的方法 |
CN109737985A (zh) * | 2018-12-06 | 2019-05-10 | 成都路行通信息技术有限公司 | 一种基于gnss角度的初始对准优化方法 |
CN109781096A (zh) * | 2017-11-15 | 2019-05-21 | 洛阳中科晶上智能装备科技有限公司 | 一种用于智能农机的组合导航定位系统和方法 |
CN109931926A (zh) * | 2019-04-04 | 2019-06-25 | 山东智翼航空科技有限公司 | 一种基于站心坐标系的小型无人机无缝自主式导航算法 |
CN109991640A (zh) * | 2017-12-29 | 2019-07-09 | 上海司南卫星导航技术股份有限公司 | 一种组合导航系统及其定位方法 |
CN110057381A (zh) * | 2019-03-20 | 2019-07-26 | 深圳市元征科技股份有限公司 | 一种导航系统的零速修正方法及系统 |
CN110231031A (zh) * | 2018-03-05 | 2019-09-13 | 高德信息技术有限公司 | 一种姿态角确定方法、装置及系统 |
CN110988950A (zh) * | 2018-10-03 | 2020-04-10 | 古野电气株式会社 | 导航装置、航行辅助信息的生成方法及存储介质 |
CN111006675A (zh) * | 2019-12-27 | 2020-04-14 | 西安理工大学 | 基于高精度重力模型的车载激光惯导系统自标定方法 |
CN111158356A (zh) * | 2018-11-08 | 2020-05-15 | 苏州宝时得电动工具有限公司 | 自动割草机及其控制方法 |
CN111578925A (zh) * | 2019-12-25 | 2020-08-25 | 重庆自行者科技有限公司 | 一种基于imu和码盘融合的车辆定位方法 |
CN111982106A (zh) * | 2020-08-28 | 2020-11-24 | 北京信息科技大学 | 导航方法、装置、存储介质及电子装置 |
CN112050807A (zh) * | 2020-08-03 | 2020-12-08 | 河北汉光重工有限责任公司 | 一种基于时间同步补偿的sins_gnss组合导航方法 |
CN112146561A (zh) * | 2020-09-09 | 2020-12-29 | 无锡卡尔曼导航技术有限公司 | 一种霍尔角度传感器安装角度偏置的估计方法 |
CN112326990A (zh) * | 2020-10-22 | 2021-02-05 | 北京农业智能装备技术研究中心 | 一种作业车辆测速方法及系统 |
CN112363249A (zh) * | 2020-09-02 | 2021-02-12 | 广东工业大学 | 一种移动式气象测量方法及装置 |
CN112577512A (zh) * | 2019-09-27 | 2021-03-30 | 北京初速度科技有限公司 | 一种基于轮速融合的状态量误差确定方法及车载终端 |
CN112577513A (zh) * | 2019-09-27 | 2021-03-30 | 北京初速度科技有限公司 | 一种状态量误差确定方法及车载终端 |
CN112710315A (zh) * | 2020-12-15 | 2021-04-27 | 广州导远电子科技有限公司 | 一种基于智能车辆的车辆定位方法及装置 |
CN112729281A (zh) * | 2020-12-18 | 2021-04-30 | 无锡卡尔曼导航技术有限公司 | 一种约束惯性/卫星组合导航静止时航向漂移的方法 |
CN112945275A (zh) * | 2021-02-07 | 2021-06-11 | 无锡卡尔曼导航技术有限公司 | 一种判断载体前进后退及组合导航航向初始对准的方法 |
CN113029139A (zh) * | 2021-04-07 | 2021-06-25 | 中国电子科技集团公司第二十八研究所 | 基于运动检测的机场飞行区车辆差分北斗/sins组合导航方法 |
CN113188539A (zh) * | 2021-04-27 | 2021-07-30 | 深圳亿嘉和科技研发有限公司 | 一种巡检机器人的组合定位方法 |
CN113256147A (zh) * | 2021-06-10 | 2021-08-13 | 农业农村部南京农业机械化研究所 | 基于联合收获机相对位移的作物产量确定方法及系统 |
CN113252048A (zh) * | 2021-05-12 | 2021-08-13 | 高新兴物联科技有限公司 | 一种导航定位方法、导航定位系统及计算机可读存储介质 |
CN113325454A (zh) * | 2021-05-18 | 2021-08-31 | 武汉大学 | 一种基于ArduRover无人车的组合定位方法 |
CN113432612A (zh) * | 2021-06-21 | 2021-09-24 | 北京信息科技大学 | 飞行体的导航方法、装置及系统 |
CN113503882A (zh) * | 2021-06-03 | 2021-10-15 | 北京自动化控制设备研究所 | 一种车载惯性/地磁组合导航方法及装置 |
CN113566850A (zh) * | 2021-07-29 | 2021-10-29 | 深圳元戎启行科技有限公司 | 惯性测量单元的安装角度标定方法、装置和计算机设备 |
CN113703019A (zh) * | 2021-08-26 | 2021-11-26 | 北京宇系航通科技有限公司 | 一种导航系统的故障处理方法、电子设备及存储介质 |
CN113984044A (zh) * | 2021-10-08 | 2022-01-28 | 杭州鸿泉物联网技术股份有限公司 | 一种基于车载多感知融合的车辆位姿获取方法及装置 |
CN114001730A (zh) * | 2021-09-24 | 2022-02-01 | 深圳元戎启行科技有限公司 | 融合定位方法、装置、计算机设备和存储介质 |
CN114199239A (zh) * | 2022-01-05 | 2022-03-18 | 中国兵器工业计算机应用技术研究所 | 结合北斗导航的双视觉辅助惯性差分座舱内头部测姿系统 |
CN114295122A (zh) * | 2021-12-02 | 2022-04-08 | 河北汉光重工有限责任公司 | 用于嵌入式系统的sins_gnss时间同步方法及系统 |
CN114413934A (zh) * | 2022-01-20 | 2022-04-29 | 北京经纬恒润科技股份有限公司 | 一种车辆定位系统校正方法和装置 |
CN114485641A (zh) * | 2022-01-24 | 2022-05-13 | 武汉梦芯科技有限公司 | 一种基于惯导卫导方位融合的姿态解算方法及装置 |
CN116718153A (zh) * | 2023-08-07 | 2023-09-08 | 成都云智北斗科技有限公司 | 一种基于gnss和ins的形变监测方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1948910A (zh) * | 2006-11-09 | 2007-04-18 | 复旦大学 | 利用gps与陀螺仪、里程计的组合定位方法与装置 |
CN101867868A (zh) * | 2010-03-26 | 2010-10-20 | 东南大学 | 一种组合导航装置和实施方法 |
US20120050095A1 (en) * | 2010-08-27 | 2012-03-01 | Bruno Scherzinger | Systems and methods for computing vertical position |
CN104697520A (zh) * | 2015-02-05 | 2015-06-10 | 南京航空航天大学 | 一体化无陀螺捷联惯导系统与gps系统组合导航方法 |
CN105021192A (zh) * | 2015-07-30 | 2015-11-04 | 华南理工大学 | 一种基于零速校正的组合导航系统的实现方法 |
CN205067758U (zh) * | 2015-09-07 | 2016-03-02 | 广西南宁至简至凡科技咨询有限公司 | 一种基于gps或ins组合的导航系统 |
CN105824039A (zh) * | 2016-03-17 | 2016-08-03 | 孙红星 | 基于里程计的克服卫星失锁时的gnss/ins车载组合定位定向算法 |
CN106154299A (zh) * | 2016-06-22 | 2016-11-23 | 陕西宝成航空仪表有限责任公司 | 一种gps/sins组合导航系统时间同步方法 |
-
2017
- 2017-01-22 CN CN201710046370.2A patent/CN106950586A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1948910A (zh) * | 2006-11-09 | 2007-04-18 | 复旦大学 | 利用gps与陀螺仪、里程计的组合定位方法与装置 |
CN101867868A (zh) * | 2010-03-26 | 2010-10-20 | 东南大学 | 一种组合导航装置和实施方法 |
US20120050095A1 (en) * | 2010-08-27 | 2012-03-01 | Bruno Scherzinger | Systems and methods for computing vertical position |
CN104697520A (zh) * | 2015-02-05 | 2015-06-10 | 南京航空航天大学 | 一体化无陀螺捷联惯导系统与gps系统组合导航方法 |
CN105021192A (zh) * | 2015-07-30 | 2015-11-04 | 华南理工大学 | 一种基于零速校正的组合导航系统的实现方法 |
CN205067758U (zh) * | 2015-09-07 | 2016-03-02 | 广西南宁至简至凡科技咨询有限公司 | 一种基于gps或ins组合的导航系统 |
CN105824039A (zh) * | 2016-03-17 | 2016-08-03 | 孙红星 | 基于里程计的克服卫星失锁时的gnss/ins车载组合定位定向算法 |
CN106154299A (zh) * | 2016-06-22 | 2016-11-23 | 陕西宝成航空仪表有限责任公司 | 一种gps/sins组合导航系统时间同步方法 |
Non-Patent Citations (6)
Title |
---|
EUN-HWAN SHIN: "Estimation Techniques for Low-Cost Inertial Navigation", 《道客巴巴》 * |
丁贝: "基于1pps的GPS/INS组合导航系统数据同步方法", 《海峡科技与产业》 * |
中国卫星导航定位协会编著: "《卫星导航定位与北斗系统应用 壮大北斗产业 创新位置服务 2014》", 31 December 2014, 北京:测绘出版社 * |
常乐等: "零速修正在GNSS/INS紧组合导航中的应用研究", 《大地测量与地球动力学》 * |
徐海鑫等: "零角度修正在 GNSS/INS 组合导航中的应用", 《测绘科学技术学报》 * |
李兵等: "基于GPS/SINS/里程计的车载组合导航研究", 《测控技术》 * |
Cited By (58)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109443347A (zh) * | 2017-07-31 | 2019-03-08 | 意法半导体股份有限公司 | 用于陆地车辆导航的系统以及对应的方法 |
CN109781096A (zh) * | 2017-11-15 | 2019-05-21 | 洛阳中科晶上智能装备科技有限公司 | 一种用于智能农机的组合导航定位系统和方法 |
CN108036797A (zh) * | 2017-11-30 | 2018-05-15 | 深圳市隐湖科技有限公司 | 基于四轮独立驱动且结合imu的里程推算方法 |
CN109991640A (zh) * | 2017-12-29 | 2019-07-09 | 上海司南卫星导航技术股份有限公司 | 一种组合导航系统及其定位方法 |
CN110231031A (zh) * | 2018-03-05 | 2019-09-13 | 高德信息技术有限公司 | 一种姿态角确定方法、装置及系统 |
US11566901B2 (en) * | 2018-06-01 | 2023-01-31 | Zhejiang Yat Electrical Appliance Co., Ltd | Integrated navigation method for mobile vehicle |
CN109059904A (zh) * | 2018-06-01 | 2018-12-21 | 浙江亚特电器有限公司 | 用于移动载具的组合导航方法 |
WO2019228437A1 (zh) * | 2018-06-01 | 2019-12-05 | 浙江亚特电器有限公司 | 用于移动载具的组合导航方法 |
CN108955851A (zh) * | 2018-07-12 | 2018-12-07 | 北京交通大学 | 利用ins和dtm确定gnss误差的方法 |
CN109143304A (zh) * | 2018-09-30 | 2019-01-04 | 百度在线网络技术(北京)有限公司 | 用于确定无人驾驶车辆位姿的方法和装置 |
CN109143304B (zh) * | 2018-09-30 | 2020-12-29 | 百度在线网络技术(北京)有限公司 | 用于确定无人驾驶车辆位姿的方法和装置 |
CN110988950A (zh) * | 2018-10-03 | 2020-04-10 | 古野电气株式会社 | 导航装置、航行辅助信息的生成方法及存储介质 |
CN111158356A (zh) * | 2018-11-08 | 2020-05-15 | 苏州宝时得电动工具有限公司 | 自动割草机及其控制方法 |
CN109737985A (zh) * | 2018-12-06 | 2019-05-10 | 成都路行通信息技术有限公司 | 一种基于gnss角度的初始对准优化方法 |
CN110057381A (zh) * | 2019-03-20 | 2019-07-26 | 深圳市元征科技股份有限公司 | 一种导航系统的零速修正方法及系统 |
CN109931926A (zh) * | 2019-04-04 | 2019-06-25 | 山东智翼航空科技有限公司 | 一种基于站心坐标系的小型无人机无缝自主式导航算法 |
CN112577513A (zh) * | 2019-09-27 | 2021-03-30 | 北京初速度科技有限公司 | 一种状态量误差确定方法及车载终端 |
CN112577512A (zh) * | 2019-09-27 | 2021-03-30 | 北京初速度科技有限公司 | 一种基于轮速融合的状态量误差确定方法及车载终端 |
CN111578925A (zh) * | 2019-12-25 | 2020-08-25 | 重庆自行者科技有限公司 | 一种基于imu和码盘融合的车辆定位方法 |
CN111006675A (zh) * | 2019-12-27 | 2020-04-14 | 西安理工大学 | 基于高精度重力模型的车载激光惯导系统自标定方法 |
CN111006675B (zh) * | 2019-12-27 | 2022-10-18 | 西安理工大学 | 基于高精度重力模型的车载激光惯导系统自标定方法 |
CN112050807A (zh) * | 2020-08-03 | 2020-12-08 | 河北汉光重工有限责任公司 | 一种基于时间同步补偿的sins_gnss组合导航方法 |
CN112050807B (zh) * | 2020-08-03 | 2023-08-18 | 河北汉光重工有限责任公司 | 一种基于时间同步补偿的sins_gnss组合导航方法 |
CN111982106A (zh) * | 2020-08-28 | 2020-11-24 | 北京信息科技大学 | 导航方法、装置、存储介质及电子装置 |
CN112363249A (zh) * | 2020-09-02 | 2021-02-12 | 广东工业大学 | 一种移动式气象测量方法及装置 |
CN112146561A (zh) * | 2020-09-09 | 2020-12-29 | 无锡卡尔曼导航技术有限公司 | 一种霍尔角度传感器安装角度偏置的估计方法 |
CN112326990A (zh) * | 2020-10-22 | 2021-02-05 | 北京农业智能装备技术研究中心 | 一种作业车辆测速方法及系统 |
CN112710315A (zh) * | 2020-12-15 | 2021-04-27 | 广州导远电子科技有限公司 | 一种基于智能车辆的车辆定位方法及装置 |
CN112729281A (zh) * | 2020-12-18 | 2021-04-30 | 无锡卡尔曼导航技术有限公司 | 一种约束惯性/卫星组合导航静止时航向漂移的方法 |
CN112945275A (zh) * | 2021-02-07 | 2021-06-11 | 无锡卡尔曼导航技术有限公司 | 一种判断载体前进后退及组合导航航向初始对准的方法 |
CN113029139A (zh) * | 2021-04-07 | 2021-06-25 | 中国电子科技集团公司第二十八研究所 | 基于运动检测的机场飞行区车辆差分北斗/sins组合导航方法 |
CN113029139B (zh) * | 2021-04-07 | 2023-07-28 | 中国电子科技集团公司第二十八研究所 | 基于运动检测的机场飞行区车辆差分北斗/sins组合导航方法 |
CN113188539A (zh) * | 2021-04-27 | 2021-07-30 | 深圳亿嘉和科技研发有限公司 | 一种巡检机器人的组合定位方法 |
CN113252048A (zh) * | 2021-05-12 | 2021-08-13 | 高新兴物联科技有限公司 | 一种导航定位方法、导航定位系统及计算机可读存储介质 |
CN113252048B (zh) * | 2021-05-12 | 2023-02-28 | 高新兴物联科技股份有限公司 | 一种导航定位方法、导航定位系统及计算机可读存储介质 |
CN113325454A (zh) * | 2021-05-18 | 2021-08-31 | 武汉大学 | 一种基于ArduRover无人车的组合定位方法 |
CN113325454B (zh) * | 2021-05-18 | 2022-06-14 | 武汉大学 | 一种基于ArduRover无人车的组合定位方法 |
CN113503882A (zh) * | 2021-06-03 | 2021-10-15 | 北京自动化控制设备研究所 | 一种车载惯性/地磁组合导航方法及装置 |
CN113503882B (zh) * | 2021-06-03 | 2023-09-12 | 北京自动化控制设备研究所 | 一种车载惯性/地磁组合导航方法及装置 |
CN113256147A (zh) * | 2021-06-10 | 2021-08-13 | 农业农村部南京农业机械化研究所 | 基于联合收获机相对位移的作物产量确定方法及系统 |
CN113432612A (zh) * | 2021-06-21 | 2021-09-24 | 北京信息科技大学 | 飞行体的导航方法、装置及系统 |
CN113566850A (zh) * | 2021-07-29 | 2021-10-29 | 深圳元戎启行科技有限公司 | 惯性测量单元的安装角度标定方法、装置和计算机设备 |
CN113566850B (zh) * | 2021-07-29 | 2024-03-08 | 深圳元戎启行科技有限公司 | 惯性测量单元的安装角度标定方法、装置和计算机设备 |
CN113703019B (zh) * | 2021-08-26 | 2023-12-22 | 北京北航天宇长鹰无人机科技有限公司 | 一种导航系统的故障处理方法、电子设备及存储介质 |
CN113703019A (zh) * | 2021-08-26 | 2021-11-26 | 北京宇系航通科技有限公司 | 一种导航系统的故障处理方法、电子设备及存储介质 |
CN114001730B (zh) * | 2021-09-24 | 2024-03-08 | 深圳元戎启行科技有限公司 | 融合定位方法、装置、计算机设备和存储介质 |
CN114001730A (zh) * | 2021-09-24 | 2022-02-01 | 深圳元戎启行科技有限公司 | 融合定位方法、装置、计算机设备和存储介质 |
CN113984044A (zh) * | 2021-10-08 | 2022-01-28 | 杭州鸿泉物联网技术股份有限公司 | 一种基于车载多感知融合的车辆位姿获取方法及装置 |
CN114295122B (zh) * | 2021-12-02 | 2024-06-11 | 河北汉光重工有限责任公司 | 用于嵌入式系统的sins_gnss时间同步方法及系统 |
CN114295122A (zh) * | 2021-12-02 | 2022-04-08 | 河北汉光重工有限责任公司 | 用于嵌入式系统的sins_gnss时间同步方法及系统 |
CN114199239A (zh) * | 2022-01-05 | 2022-03-18 | 中国兵器工业计算机应用技术研究所 | 结合北斗导航的双视觉辅助惯性差分座舱内头部测姿系统 |
CN114199239B (zh) * | 2022-01-05 | 2024-04-16 | 中国兵器工业计算机应用技术研究所 | 结合北斗导航的双视觉辅助惯性差分座舱内头部测姿系统 |
CN114413934B (zh) * | 2022-01-20 | 2024-01-26 | 北京经纬恒润科技股份有限公司 | 一种车辆定位系统校正方法和装置 |
CN114413934A (zh) * | 2022-01-20 | 2022-04-29 | 北京经纬恒润科技股份有限公司 | 一种车辆定位系统校正方法和装置 |
CN114485641A (zh) * | 2022-01-24 | 2022-05-13 | 武汉梦芯科技有限公司 | 一种基于惯导卫导方位融合的姿态解算方法及装置 |
CN114485641B (zh) * | 2022-01-24 | 2024-03-26 | 武汉梦芯科技有限公司 | 一种基于惯导卫导方位融合的姿态解算方法及装置 |
CN116718153B (zh) * | 2023-08-07 | 2023-10-27 | 成都云智北斗科技有限公司 | 一种基于gnss和ins的形变监测方法及系统 |
CN116718153A (zh) * | 2023-08-07 | 2023-09-08 | 成都云智北斗科技有限公司 | 一种基于gnss和ins的形变监测方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106950586A (zh) | 用于农机作业的gnss/ins/车辆组合导航方法 | |
CN109059904A (zh) | 用于移动载具的组合导航方法 | |
KR102018450B1 (ko) | 단일 안테나에 기초한 gnss―ins 차량 자세 결정 방법 | |
CN108535755B (zh) | 基于mems的gnss/imu车载实时组合导航方法 | |
CN108180925B (zh) | 一种里程计辅助车载动态对准方法 | |
CN108594283B (zh) | Gnss/mems惯性组合导航系统的自由安装方法 | |
CN108051866B (zh) | 基于捷联惯性/gps组合辅助水平角运动隔离的重力测量方法 | |
CN112378400A (zh) | 一种双天线gnss辅助的捷联惯导组合导航方法 | |
CN108680942B (zh) | 一种惯性/多天线gnss组合导航方法及装置 | |
CA3003298A1 (en) | Gnss and inertial navigation system utilizing relative yaw as an observable for an ins filter | |
CN108088443B (zh) | 一种定位定向设备速度补偿方法 | |
CN109855617A (zh) | 一种车辆定位方法、车辆定位装置及终端设备 | |
CN110133694B (zh) | 基于双天线gnss航向和轮速辅助的车辆定位方法及系统 | |
CN107063241A (zh) | 基于双gnss天线及单轴mems陀螺的前轮测角系统 | |
CN104880189B (zh) | 一种动中通天线低成本跟踪抗干扰方法 | |
CN114076610B (zh) | Gnss/mems车载组合导航系统的误差标定、导航方法及其装置 | |
JP2012173190A (ja) | 測位システム、測位方法 | |
CN110133695A (zh) | 一种双天线gnss位置延迟时间动态估计系统及方法 | |
CN103644911A (zh) | 陀螺仪辅助定位方法 | |
CN114526731A (zh) | 一种基于助力车的惯性组合导航方向定位方法 | |
CN116255998A (zh) | 一种基于多惯导互观测的惯导振荡误差精确估计方法 | |
CN115727843A (zh) | 用于航位推算的轮速确定方法、装置及设备 | |
CN111307114A (zh) | 基于运动参考单元的水面舰船水平姿态测量方法 | |
CN113008229B (zh) | 一种基于低成本车载传感器的分布式自主组合导航方法 | |
CN112484731B (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170714 |
|
RJ01 | Rejection of invention patent application after publication |