CN113932811A - 一种新的地形匹配导航系统、方法 - Google Patents

一种新的地形匹配导航系统、方法 Download PDF

Info

Publication number
CN113932811A
CN113932811A CN202110981651.3A CN202110981651A CN113932811A CN 113932811 A CN113932811 A CN 113932811A CN 202110981651 A CN202110981651 A CN 202110981651A CN 113932811 A CN113932811 A CN 113932811A
Authority
CN
China
Prior art keywords
positioning
point
navigation
underwater vehicle
terrain
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
CN202110981651.3A
Other languages
English (en)
Other versions
CN113932811B (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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN202110981651.3A priority Critical patent/CN113932811B/zh
Publication of CN113932811A publication Critical patent/CN113932811A/zh
Application granted granted Critical
Publication of CN113932811B publication Critical patent/CN113932811B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/20Instruments for performing navigational calculations

Landscapes

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

Abstract

本发明公开了一种新的地形匹配导航系统、方法,包括初始定位模块、滤波器初始化模块、初始对准滤波过程模块、滤波器状态监测模块一、跟踪导航滤波过程模块、滤波器状态监测模块二和定位输出模块。在地形匹配定位似然函数存在多个伪波峰的情况下,本发明所述地形匹配导航系统、方法可获得较传统方法更加稳定的定位信息和更加高效的滤波收敛速度。

Description

一种新的地形匹配导航系统、方法
技术领域
本发明涉及一种地形匹配导航系统,尤其涉及一种新的地形匹配导航系统、方法。
背景技术
地形匹配导航采用环境特征匹配与跟踪方式进行定位解算,其本质是一种组合导航系统,由主导航系统(一般为惯性导航)和地形匹配定位系统通过信息融合模块组合而成,而地形匹配定位信息仅仅作为主导航系统的辅助信息。由于主导航系统往往存在随时间无限累积的误差,在地形匹配导航系统启动时刻,主导航系统往往会输出较大的定位误差,导致地形匹配导航起始时刻的匹配定位误差和融合定位误差均较大,由此引起较大的地形匹配导航系统滤波器初始化误差。在滤波器迭代过程中,初始化误差随迭代过程向前传递导致滤波过程及其不稳定。
发明内容
本发明的目的是提供一种新的地形匹配导航系统、方法,可获得较传统方法更加稳定的定位信息和更加高效的滤波收敛速度。
本发明为实现上述发明目的采用如下技术方案:
本发明提供了一种新的地形匹配导航系统,包括初始定位模块、滤波器初始化模块、初始对准滤波过程模块、滤波器状态监测模块一、跟踪导航滤波过程模块、滤波器状态监测模块二和定位输出模块;
所述初始定位模块获取系统在导航起点的初始定位信息,将该初始定位信息输入到所述滤波器初始化模块;
所述滤波器初始化模块获取到初始定位信息后执行滤波器初始化过程,并获得初始化信息,将该初始化信息输入到所述初始对准滤波过程模块;
所述初始对准滤波过程模块中的滤波器通过获取所述滤波器初始化模块输入的初始化信息完成滤波器启动,并实时地输出地形匹配系统的初始对准信息,将该信息输入到所述滤波器状态监测模块一;
所述滤波器状态监测模块一判断所述初始对准滤波过程模块中的初始对准信息是否收敛,如果初始对准信息收敛,则系统完成初始对准,进入所述跟踪导航滤波过程模块,如果初始对准信息未收敛,则继续执行所述初始对准滤波过程模块;
所述跟踪导航滤波过程模块的滤波器输出地形匹配导航信息至所述滤波器状态监测模块二,如果所述滤波器状态监测模块二判断所述跟踪导航滤波过程模块输出的地形匹配导航信息未发散,则所述跟踪导航滤波过程的滤波器输出的地形匹配导航信息为有效的地形匹配导航信息,否则,则系统重新执行所述初始定位模块。
进一步地,所述初始定位信息包括定位点、定位误差及定位概率分布函数。
本发明提供了一种新的地形匹配导航系统、方法,采用上述一种新的地形匹配导航系统,所述方法包括:
进行地形匹配导航系统的初始匹配定位模块;
在完成初始匹配定位后进行地形匹配导航系统滤波器初始化模块;
在完成地形匹配导航系统滤波器初始化模块之后,开始运行系统初始对准滤波过程模块;
在运行系统初始对准滤波过程模块的过程中,根据系统初始对准滤波过程模块的输出信息实时监测输出信息的收敛性;
在监测出系统初始对准滤波过程模块的输出信息收敛后,执行地形匹配导航跟踪导航滤波模块;
在执行地形匹配导航跟踪导航滤波模块的过程中,根据地形匹配导航跟踪导航滤波模块的输出信息实时监测输出信息的收敛性,如果输出信息始终保持收敛,则系统输出定位信息,如果监测到输出信息发散,则进行重新初始匹配定位。
进一步地,进行地形匹配导航系统的初始匹配定位模块的方法包括:
在地形匹配导航规划点的起点wpi,i=1,基于主导航系统给出的导航定位点
Figure BDA0003229219220000021
i=1和导航误差
Figure BDA0003229219220000022
i=1计算主导航系统的误差椭圆及矩形搜索区间;
以矩形搜索区间为匹配搜索区间,获得地形匹配导航规划点的起点wpi,i=1的地形匹配定位点
Figure BDA0003229219220000023
i=1和地形匹配定位点的定位误差
Figure BDA0003229219220000024
i=1:
其中,ki表示地形匹配定位点的索引号,ki=1,2…Ki,Ki表示地形匹配导航规划点wpi获得的地形匹配定位点总数。
进一步地,在完成初始匹配定位模块后进行地形匹配导航系统滤波器初始化模块的方法包括:
在地形匹配导航规划点的起点wpi;i=1,根据主导航系统给出的主导航定位点
Figure BDA0003229219220000025
i=1、主导航定位点的定位误差
Figure BDA0003229219220000026
i=1、地形匹配定位点
Figure BDA0003229219220000027
i=1和地形匹配定位点的定位误差
Figure BDA0003229219220000028
i=1,计算水下航行器路径定位点
Figure BDA0003229219220000029
i=1和水下航行器路径定位点的定位误差
Figure BDA00032292192200000210
i=1,计算公式如下:
Figure BDA00032292192200000211
式中,ni∈{1,2,3…Ni}表示地形匹配导航规划点wpi点的水下航行器路径定位点的索引,Ni表示有效路径的总数,其中N0=1,n0=1;
在地形匹配导航规划点的起点wpi,i=1,计算水下航行器路径定位点的估计代价
Figure BDA0003229219220000031
Ni=ni-1×ki表示wpi点获得的有效路径总数,公式如下:
Figure BDA0003229219220000032
Figure BDA0003229219220000033
按照由小到大的顺序排序,得到序列
Figure BDA0003229219220000034
保存
Figure BDA0003229219220000035
对应在
Figure BDA0003229219220000036
中的索引号;
根据
Figure BDA0003229219220000037
的索引号对地形匹配导航规划点wpi获得的水下航行器路径定位点
Figure BDA0003229219220000038
水下航行器路径定位点的定位误差
Figure BDA0003229219220000039
水下航行器路径定位点的定位偏差
Figure BDA00032292192200000310
水下航行器路径定位点的定位误差的方差下确界
Figure BDA00032292192200000311
进行排序,得到序列
Figure BDA00032292192200000312
Figure BDA00032292192200000313
进一步地,在完成地形匹配导航系统初始化模块之后,开始执行系统初始对准滤波过程模块的方法包括:
在i>1的时,水下航行器经过水下航行器路径定位点
Figure BDA00032292192200000314
到达wpi+1点的过程中,主导航系统给出的航行距离Di,i+1,主导航系统给出的主导航误差增值Dei,i+1,则水下航行器经过wpi点的Ni条路径定位点序列
Figure BDA00032292192200000315
到达wpi+1点后,主导航系统的定位点为:
Figure BDA00032292192200000316
主导航定位误差为:
Figure BDA00032292192200000317
根据wpi+1点获得的地形匹配定位点
Figure BDA00032292192200000318
地形匹配定位误差
Figure BDA00032292192200000319
计算水下航行器路径定位点序列
Figure BDA00032292192200000320
和水下航行器路径定位点的定位误差
Figure BDA00032292192200000321
其中ni+1∈{1,2,3…Ni+1}表示wpi+1点的水下航行器路径定位点索引,ki+1∈{1,2,3…Ki+1}为wpi+1点的地形匹配定位点索引号,公式如下:
Figure BDA0003229219220000041
在获得水下航行器路径定位点
Figure BDA0003229219220000042
和水下航行器路径定位点的定位误差
Figure BDA0003229219220000043
后,计算④中的Ni+1条水下航行器路径定位点的估计代价
Figure BDA0003229219220000044
公式如下:
Figure BDA0003229219220000045
式中,k∈{1,2,3…Ki+1}为wpi+1点的地形匹配定位点索引号,ni∈{1,2,3…Ni}为wpi点的水下航行器路径定位点索引号;
Figure BDA0003229219220000046
按照由小到大的顺序排序,得到序列
Figure BDA0003229219220000047
保存
Figure BDA0003229219220000048
对应在
Figure BDA0003229219220000049
中的索引号;
根据
Figure BDA00032292192200000410
的索引号对地形匹配导航规划点wpi+1获取的水下航行器路径定位点
Figure BDA00032292192200000411
水下航行器路径定位点的定位误差
Figure BDA00032292192200000412
水下航行器路径定位点的定位偏
Figure BDA00032292192200000413
水下航行器路径定位点的定位误差的方差下确界
Figure BDA00032292192200000414
进行排序,得到序列
Figure BDA00032292192200000415
Figure BDA00032292192200000416
根据水下航行器路径定位点的定位偏
Figure BDA00032292192200000417
和水下航行器路径定位点的定位误差
Figure BDA00032292192200000418
对水下航行器路径定位点的有效性进行筛选;
获得筛选后的水下航行器路径定位点、水下航行器路径定位点的定位误差分别为:Ddte(m)、Ddt(m),并将此序列重新赋值给
Figure BDA00032292192200000419
从而获得新的水下航行器路径定位点、水下航行器路径定位点的定位误差;
根据新的水下航行器路径定位点
Figure BDA00032292192200000420
和水下航行器路径定位点的定位误差
Figure BDA00032292192200000421
及水下航行器路径定位点的估计代价
Figure BDA0003229219220000051
选出水下航行器路径定位点的估计代价中的极小值
Figure BDA0003229219220000052
并提取其在序列
Figure BDA0003229219220000053
中对应的索引号
Figure BDA0003229219220000054
根据索引号
Figure BDA0003229219220000055
选择水下航行器路径定位点
Figure BDA0003229219220000056
中与索引号
Figure BDA0003229219220000057
对应的最优水下航行器路径定位点
Figure BDA0003229219220000058
随后执行水下航行器路径定位过程的收敛性判断。
进一步地,进行系统初始对准滤波过程模块的输出信息的状态检测和收敛性判断的方法包括:
地形匹配导航系统状态检测模块一根据获得的水下航行器路径定位点
Figure BDA0003229219220000059
和水下航行器路径定位点的定位误差
Figure BDA00032292192200000510
及水下航行器路径定位点的估计代价
Figure BDA00032292192200000511
计算水下航行器路径定位点的融合定位点
Figure BDA00032292192200000512
α是为了确保
Figure BDA00032292192200000513
成立的归一化矩阵,I为单位矩阵;
计算最优水下航行器路径定位点
Figure BDA00032292192200000514
Figure BDA00032292192200000515
之间的偏差
Figure BDA00032292192200000516
Figure BDA00032292192200000517
水下航行器路径定位收敛判断关系式:D_Ddtlimlnli+1min<γ·Ddtei+1ni+1,1<γ<3,是否成立,如果成立,则表明水下航行器路径定位满足收敛要求,地形匹配导航系统初始对准滤波过程模块的输出信息收敛,完成系统初始对准过程模块的运行,i+1时刻的最优水下航行器路径定位点
Figure BDA00032292192200000518
即为地形匹配导航输出,同时地形匹配导航系统进入地形匹配导航跟踪导航滤波模块。
本发明的有益效果如下:
在地形匹配定位似然函数存在多个伪波峰的情况下,本发明所述地形匹配导航系统、方法可获得较传统方法更加稳定的定位信息和更加高效的滤波收敛速度;
本发明所述新的地形匹配导航系统包含系统状态监测模块一和系统状态监测模块二,该模块通过监测和反馈系统滤波器运行状态,可实现地形匹配导航系统运行状态的自主感知和运行状态发散情况下的自主重新初始对准;
新的地形匹配导航系统可自主感知地形匹配导航系统滤波器发散运行状态,并在地形匹配导航系统滤波器发散状态下自主控制地形匹配导航系统完成系统重新初始定位,新的地形匹配导航系统运行架构可使地形匹配导航系统具有更高的运行可靠性;
新的地形匹配导航系统可自主感知地形匹配导航系统滤波器收敛运行状态,并在地形匹配导航系统滤波器收敛状态下自主完成导航模式的切换,将系统导航模式由初始对准模式切换至跟踪导航模式,并输出有效的地形匹配导航信息。
附图说明
图1为根据本发明实施例提供的一种新的地形匹配导航系统的框架图;
图2为根据本发明实施例提供的一种新的地形匹配导航方法中初始定位误差椭圆及初始地形匹配搜索区间图;
图3为根据本发明实施例提供的一种新的地形匹配导航方法中地形匹配定位点及其置信区间图;
图4为根据本发明实施例提供的一种新的地形匹配导航方法中初始地形匹配定位及其到达下一wpi的航迹图;
图5为根据本发明实施例提供的一种新的地形匹配导航方法中wp1点的有效航迹和wp2点的地形匹配定位点图;
图6为根据本发明实施例提供的一种新的地形匹配导航方法中地形匹配导航初始对准过程中数据结果和流程图。
具体实施方式
下面结合具体实施例对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如图1所示,本发明公开的一种新的地形匹配导航系统,①包括地形匹配导航系统初始定位模块1、地形匹配导航系统滤波初始化模块2、地形匹配导航系统初始对准滤波模块3、地形匹配导航系统滤波器状态监测模块一4、地形匹配导航系统收敛性判断5、地形匹配导航系统跟踪导航滤波6、地形匹配导航系统滤波器状态监测模块二7、地形匹配导航系统滤波发散判定8、地形匹配导航系统输出9。其中地形匹配导航系统初始定位1、地形匹配导航系统滤波初始化2、地形匹配导航系统初始对准滤波3、地形匹配导航系统滤波器状态监测一4构成了地形匹配导航系统的初始对准模块10,该模块也是本算法运行框架与其他地形匹配导航系统不同之处。
②在新的地形匹配导航算法运行框架中的地形匹配导航系统初始对准滤波过程模块3和地形匹配导航跟踪滤波过程模块6后面均加入了地形匹配导航滤波器状态检测模块一4和地形匹配导航滤波器状态检测模块二7;在地形匹配导航滤波器状态检测模块一4和地形匹配导航滤波器状态检测模块二7后又分别加入了系统滤波器状态收敛检测模块5和滤波器发散检测模块8,这也是本算法运行框架与其他算法运行框架不同之处。
③地形匹配导航系统运行框架中的各个模块及其内部算法的执行流程概述如下:
a)地形匹配导航系统初始定位模块1获取系统在导航起点的定位点、定位误差及定位概率分布信息,将该信息输入到地形匹配导航系统滤波器初始化模块2;
b)地形匹配导航系统滤波器初始化模块2获取到初始定位信息后执行滤波器初始化过程,并获得地形匹配导航系统初始对准过程滤波器的初始化信息,包括初始定位点、定位误差、定位概率分布等滤波器初始化信息,该信息直接输入到地形匹配导航系统初始对准滤波过程模块3。
c)地形匹配导航系统初始对准滤波过程模块3执行地形匹配导航系统的初始对准过程,在该过程中模块3通过获取地形匹配导航系统滤波器初始化模块2的初始化信息完成初始对准滤波器启动,并实时输出地形匹配导航的定位信息、定位误差信息、实时测量地形、地形测量误差等信息,这些信息被输入到地形匹配导航滤波器状态监测模块一4。
d)地形匹配导航滤波器状态监测模块一4判断地形匹配导航系统初始对准滤波过程模块3输出的导航信息是否收敛,如果模块3输出的导航信息是收敛的,则系统完成初始对准,进入跟踪导航阶段。如果尚未收敛,则继续执行地形匹配导航系统初始对准滤波过程模块3。
e)在地形匹配导航滤波器状态监测模块一4判断地形匹配导航系统初始对准滤波过程模块3输出的导航信息是收敛的情况下,系统将运行地形匹配导航跟踪导航滤波模块6,进入跟踪导航阶段。跟踪导航阶段同样是执行连续的定位和定位误差估计,跟踪导航阶段的滤波器初始值即在图1地形匹配导航滤波器状态监测模块一判断初始对准过程滤波收敛后,初始对准滤波器输出的定位信息作为地形匹配导航跟踪导航滤波模块6的初始化信息。
f)地形匹配导航跟踪导航滤波模块6的输出信息同样的被输入到地形匹配导航滤波器状态监测模块二7中,如果地形匹配导航滤波器状态监测模块二判断地形匹配导航跟踪导航滤波模块6的输出信息是稳定和收敛的,也即是跟踪导航滤波器状态没有发散,则地形匹配导航跟踪导航滤波模块6输出的定位信息即是有效的地形匹配导航定位信息。否则,跟踪导航滤波器发散,则系统进入重新初始对准过程,此时地形匹配导航跟踪导航滤波模块6输出定位点、定位误差等信息,并执行地形匹配导航系统初始定位模块1。
新的地形匹配导航方法的执行过程为:
①地形匹配导航系统的初始匹配定位模块执行过程。初始匹配定位模块也即是图1中模块1,模块1执行后将获取地形匹配导航系统初始定位信息。假设地形匹配导航规划点wpi,其中i表示地形匹配导航规划点索引号,主导航系统在地形匹配导航规划点wpi输出导航定位点Di和导航误差Dei。当i=1时表示第一个地形匹配导航规划点,也就是地形匹配导航的起点。因此,主导航定位点位于wpi,i=1时,将触发地形匹配导航系统的初始匹配定位模块执行。初始匹配定位的执行过程如下:
a)基于主导航系统给出的导航定位点D111和导航误差De1计算主导航系统的误差椭圆12及矩形搜索区间13。
b)以矩形搜索区间13为匹配搜索区间,计算该区间内的地形匹配定位概率分布,并获得地形匹配定位点15、16和地形匹配定位置信区间17、18。
c)地形匹配定位可能存在多个地形匹配定位点和地形匹配定位置信区间。
Figure BDA0003229219220000071
表示地形匹配定位点,
Figure BDA0003229219220000072
表示地形匹配定位误差,其中i表示地形匹配导航规划点的索引号,
ki∈{1,2…Ki}表示地形匹配定位获得的定位点索引号,K1表示wp1点获得的地形匹配定位点总数(即i=1,地形匹配到起点)。
②在完成初始匹配定位模块执行过程①后进行进入地形匹配导航系统滤波器初始化模块执行过程,其滤波器初始化模块的执行过程的流程如下:
a)在地形匹配导航起点wp1(即i=1时),根据主导航系统给出的主导航定位点D1,主导航误差De1,地形匹配定位点
Figure BDA0003229219220000081
地形匹配定位误差
Figure BDA0003229219220000082
计算水下航行器路径定位点
Figure BDA0003229219220000083
i=1和水下航行器路径定位点的定位误差
Figure BDA0003229219220000084
i=1。
Figure BDA0003229219220000085
式中,ni∈{1,2,3…Ni}表示地形匹配导航规划点wpi点的水下航行器路径定位点索引,Ni表示水下航行器路径定位点的总数,其中N0=1,ni=1。
b)然后,计算i=1时水下航行器路径定位点估计的代价函数值
Figure BDA0003229219220000086
ki=1,2…K1,K1表示在地形匹配导航规划点wp1获得的地形匹配定位点总数,ki表示地形匹配定位点的索引号,N1=ni-1×ki表示在地形匹配导航规划点wp1获得的水下航行器路径定位点总数。
Figure BDA0003229219220000087
c)然后,将
Figure BDA0003229219220000088
按照由小到大的顺序排序,得到序列
Figure BDA0003229219220000089
ni∈{1,2,3…Ni},
Figure BDA00032292192200000810
保存
Figure BDA00032292192200000811
对应在
Figure BDA00032292192200000812
中的索引号。根据
Figure BDA00032292192200000813
的索引号对地形匹配导航规划点wpi获得的水下航行器路径定位点
Figure BDA00032292192200000814
水下航行器路径定位点定位误差
Figure BDA00032292192200000815
水下航行器路径定位点定位偏差
Figure BDA00032292192200000816
水下航行器路径定位点定位误差的方差下确界
Figure BDA00032292192200000817
进行排序,得到序列
Figure BDA00032292192200000818
ni∈{1,2,3…Ni}。
③在完成地形匹配导航系统初始化模块的执行过程②之后,进入地形匹配导航系统初始对准滤波模块的执行过程,该过程的具体流程如下。
a)水下航行器根据wpi点获得的水下航行器路径定位点继续向前航行到达下一个地形匹配导航规划点wpi+1。如图4所示,11和21分别表示在wpi点获得的水下航行器路径定位点对应的航迹,水下航行器可通过19和22到达wpi+1
获得水下航行器由wpi点的Ni个水下航行器路径定位点到达wpi+1点的主导航距离:Di,i+1,主导航误差为Dei,i+1,则水下航行器经过地形匹配导航路径规划点wpi的水下航行器路径定位点由wpi到达地形匹配导航规划点wpi+1后的主导航定位点为:
Figure BDA0003229219220000091
主导航误差为:
Figure BDA0003229219220000092
根据地形匹配导航系统在wpi+1点获得的地形匹配定位点
Figure BDA0003229219220000093
地形匹配定位误差
Figure BDA0003229219220000094
计算水下航行器路径定位点
Figure BDA0003229219220000095
和水下航行器路径定位点的定位误差
Figure BDA0003229219220000096
其中ni+1∈{1,2,3…Ni+1}表示wpi+1点的水下航行器路径定位点索引,ki+1∈{1,2,3…Ki+1}为wpi+1点的地形匹配定位点索引号。如图5所示,22表示wp2的地形匹配定位搜索区间,23、24、25为wp2点的地形匹配定位点,26、27、28分别表示其定位置信区间,20和23表示由wp1点的水下航行器路径定位点到达wp2点的主导航点,19和22为wp1的水下航行器路径定位点对应的主导航航迹,则wp2点的有效航迹数为N2=2×3=6。
Figure BDA0003229219220000097
在获得水下航行器路径定位点
Figure BDA0003229219220000098
后,计算④中的Ni+1个水下航行器路径定位点的估计代价值
Figure BDA0003229219220000099
Figure BDA00032292192200000910
式中,k∈{1,2,3…Ki+1}为wpi+1点的地形匹配定位点索引号,ni∈{1,2,3…Ni}为wpi点的有效航迹的索引号。
b)将
Figure BDA00032292192200000911
按照由小到大的顺序排序,得到序列
Figure BDA00032292192200000912
ni+1∈{1,2,3…Ni+1},
Figure BDA00032292192200000913
保存
Figure BDA00032292192200000914
对应在
Figure BDA00032292192200000915
中的索引号。
其次,根据
Figure BDA00032292192200000916
的索引号对地形匹配导航规划点wpi+1的水下航行器路径定位点
Figure BDA00032292192200000917
水下航行器路径定位点的定位误差
Figure BDA00032292192200000918
水下航行器路径定位点的定位偏差
Figure BDA00032292192200000919
水下航行器路径定位误差的方差下确界
Figure BDA0003229219220000101
进行排序,得到序列
Figure BDA0003229219220000102
Figure BDA0003229219220000103
Figure BDA0003229219220000104
至此,在地形匹配导航路径规划点wpi,i≥2获得经过筛选的
Figure BDA0003229219220000105
个水下航行器路径定位点、同时获得筛选后的水下航行器路径定位点的定位误差、水下航行器路径定位点的估计代价分别为:
Figure BDA0003229219220000106
c)根据b)中的序列
Figure BDA0003229219220000107
选出其中的极小值
Figure BDA0003229219220000108
并提取其在序列
Figure BDA0003229219220000109
中对应的索引号
Figure BDA00032292192200001010
根据索引号
Figure BDA00032292192200001011
选择水下航行器路径定位点
Figure BDA00032292192200001012
中与索引号
Figure BDA00032292192200001013
对应的最优水下航行器路径定位点
Figure BDA00032292192200001014
d)在每一个地形匹配导航规划点wpi进行定位算法的收敛性判断。具体如下:首先,地形匹配导航系统状态检测模块一根据获得的水下航行器路径定位点
Figure BDA00032292192200001015
和水下航行器路径定位点的定位误差
Figure BDA00032292192200001016
及水下航行器路径定位点的估计代价
Figure BDA00032292192200001017
计算水下航行器路径定位点的融合定位点
Figure BDA0003229219220000111
α是为了确保
Figure BDA0003229219220000112
成立的归一化矩阵,I为单位矩阵。然后,计算最优水下航行器路径定位点
Figure BDA0003229219220000113
之间的偏差
Figure BDA0003229219220000114
其次,水下航行器路径定位收敛判断关系式:
Figure BDA0003229219220000115
1<γ<3,是否成立,如果成立,则表明水下航行器路径定位满足收敛要求,地形匹配导航系统初始对准滤波过程模块的输出信息收敛,完成系统初始对准过程模块的运行,i+1时刻的最优水下航行器路径定位点
Figure BDA0003229219220000116
即为地形匹配导航输出,同时地形匹配导航系统进入地形匹配导航跟踪导航滤波模块。如果不满足要求,水下航行器继续航行,到达下一个地形匹配导航规划点wpi+2,i≥2,在wpi+2,i≥2点重复步骤③,直到水下航行器路径定位收敛判断关系式满足为止。随着i值(29)的增加,形成如图6所示的数据流程图。
④在完成系统初始对准滤波过程模块执行过程③后,地形匹配导航系统初始对准滤波过程模块的输出信息收敛,地形匹配导航系统进入地形匹配导航跟踪导航滤波模块的执行阶段。地形匹配导航跟踪导航滤波模块在运行过程中向地形匹配导航滤波器状态检测模块二输入地形匹配导航信息、测量地形和地形测量误差信息等,地形匹配导航滤波器状态检测模块二根据这些信息判断地形匹配导航跟踪导航滤波模块的输出信息是否发散,以此判断该模块的运行状态是否稳定,若地形匹配导航跟踪导航滤波模块的输出信息发散,则地形匹配导航跟踪导航滤波模块的运行状态不稳定,则地形匹配导航系统需要重新回到地形匹配导航初始定位模块,并重新完成上述程序。若地形匹配导航跟踪导航滤波模块的输出信息没有发散,则地形匹配导航跟踪导航滤波模块输出有效的导航定位信息。
以上对本发明的较佳实施进行了具体说明,当然,本发明还可以采用与上述实施方式不同的形式,熟悉本领域的技术人员在不违背本发明精神的前提下所作的等同的变换或相应的改动,都应该属于本发明的保护范围内。

Claims (7)

1.一种新的地形匹配导航系统,其特征在于,包括初始定位模块、滤波器初始化模块、初始对准滤波过程模块、滤波器状态监测模块一、跟踪导航滤波过程模块、滤波器状态监测模块二和定位输出模块;
所述初始定位模块获取系统在导航起点的初始定位信息,将该初始定位信息输入到所述滤波器初始化模块;
所述滤波器初始化模块获取到初始定位信息后执行滤波器初始化过程,并获得初始化信息,将该初始化信息输入到所述初始对准滤波过程模块;
所述初始对准滤波过程模块中的滤波器通过获取所述滤波器初始化模块输入的初始化信息完成滤波器启动,并实时地输出地形匹配系统的初始对准信息,将该信息输入到所述滤波器状态监测模块一;
所述滤波器状态监测模块一判断所述初始对准滤波过程模块中的初始对准信息是否收敛,如果初始对准信息收敛,则系统完成初始对准,进入所述跟踪导航滤波过程模块,如果初始对准信息未收敛,则继续执行所述初始对准滤波过程模块;
所述跟踪导航滤波过程模块的滤波器输出地形匹配导航信息至所述滤波器状态监测模块二,如果所述滤波器状态监测模块二判断所述跟踪导航滤波过程模块输出的地形匹配导航信息未发散,则所述跟踪导航滤波过程的滤波器输出的地形匹配导航信息为有效的地形匹配导航信息,否则,则系统重新执行所述初始定位模块。
2.根据权利要求1所述的一种新的地形匹配导航系统,其特征在于,所述初始定位信息包括定位点、定位误差及定位概率分布函数。
3.一种新的地形匹配导航方法,其特征在于,采用如权利要求1或2所述的一种新的地形匹配导航系统,所述方法包括:
进行地形匹配导航系统的初始匹配定位模块;
在完成初始匹配定位后进行地形匹配导航系统滤波器初始化模块;
在完成地形匹配导航系统滤波器初始化模块之后,开始运行系统初始对准滤波过程模块;
在运行系统初始对准滤波过程模块的过程中,根据系统初始对准滤波过程模块的输出信息实时监测输出信息的收敛性;
在监测出系统初始对准滤波过程模块的输出信息收敛后,执行地形匹配导航跟踪导航滤波模块;
在执行地形匹配导航跟踪导航滤波模块的过程中,根据地形匹配导航跟踪导航滤波模块的输出信息实时监测输出信息的是否发散,如果输出信息没有发散,则系统输出定位信息,如果监测到输出信息发散,则进行回到初始匹配定位模块,进行重新进行初始对准。
4.根据权利要求3所述的一种新的地形匹配导航方法,其特征在于,进行地形匹配导航系统的初始匹配定位模块的方法包括:
在地形匹配导航规划点的起点wpi,i=1,基于主导航系统给出的导航定位点
Figure FDA0003229219210000021
和导航误差
Figure FDA0003229219210000022
计算主导航系统的误差椭圆及矩形搜索区间;
以矩形搜索区间为匹配搜索区间,获得地形匹配导航规划点的起点wpi,i=1的地形匹配定位点
Figure FDA0003229219210000023
和地形匹配定位点的定位误差
Figure FDA0003229219210000024
其中,ki表示地形匹配定位点的索引号,ki=1,2...Ki,Ki表示地形匹配导航规划点wpi获得的地形匹配定位点总数。
5.根据权利要求4所述的一种新的地形匹配导航方法,其特征在于,在完成初始匹配定位模块后进行地形匹配导航系统滤波器初始化模块的方法包括:
在地形匹配导航规划点的起点wpi;i=1,根据主导航系统给出的主导航定位点
Figure FDA0003229219210000025
主导航定位点的定位误差
Figure FDA0003229219210000026
地形匹配定位点
Figure FDA0003229219210000027
和地形匹配定位点的定位误差
Figure FDA0003229219210000028
计算水下航行器路径定位点
Figure FDA0003229219210000029
和水下航行器路径定位点的定位误差
Figure FDA00032292192100000210
计算公式如下:
Figure FDA00032292192100000211
式中,ni∈{1,2,3...Ni}表示地形匹配导航规划点wpi点的水下航行器路径定位点的索引,Ni表示有效路径的总数,其中N0=1,n0=1;
在地形匹配导航规划点的起点wpi,i=1,计算水下航行器路径定位点的估计代价
Figure FDA00032292192100000212
Ni=ni-1×ki表示wpi点获得的有效路径总数,公式如下:
Figure FDA00032292192100000213
Figure FDA0003229219210000031
按照由小到大的顺序排序,得到序列
Figure FDA0003229219210000032
Figure FDA0003229219210000033
保存
Figure FDA0003229219210000034
对应在
Figure FDA0003229219210000035
中的索引号;
根据
Figure FDA0003229219210000036
的索引号对地形匹配导航规划点wpi获得的水下航行器路径定位点
Figure FDA0003229219210000037
水下航行器路径定位点的定位误差
Figure FDA0003229219210000038
水下航行器路径定位点的定位偏差
Figure FDA0003229219210000039
水下航行器路径定位点的定位误差的方差下确界
Figure FDA00032292192100000310
进行排序,得到序列
Figure FDA00032292192100000311
Figure FDA00032292192100000312
6.根据权利要求5所述的一种新的地形匹配导航方法,其特征在于,在完成地形匹配导航系统初始化模块之后,开始执行系统初始对准滤波过程模块的方法包括:
在i>1的时,水下航行器经过水下航行器路径定位点
Figure FDA00032292192100000313
到达wpi+1点的过程中,主导航系统给出的航行距离Di,i+1,主导航系统给出的主导航误差增值Dei,i+1,则水下航行器经过wpi点的Ni条路径定位点序列
Figure FDA00032292192100000314
到达wpi+1点后,主导航系统的定位点为:
Figure FDA00032292192100000315
主导航定位误差为:
Figure FDA00032292192100000316
根据wpi+1点获得的地形匹配定位点
Figure FDA00032292192100000317
地形匹配定位误差
Figure FDA00032292192100000318
计算水下航行器路径定位点序列
Figure FDA00032292192100000319
和水下航行器路径定位点的定位误差
Figure FDA00032292192100000320
其中ni+1∈{1,2,3...Ni+1}表示wpi+1点的水下航行器路径定位点索引,ki+1∈{1,2,3...Ki+1}为wpi+1点的地形匹配定位点索引号,公式如下:
Figure FDA00032292192100000321
在获得水下航行器路径定位点
Figure FDA00032292192100000322
和水下航行器路径定位点的定位误差
Figure FDA00032292192100000323
后,计算④中的Ni+1条水下航行器路径定位点的估计代价
Figure FDA00032292192100000324
公式如下:
Figure FDA0003229219210000041
式中,k∈{1,2,3...Ki+1}为wpi+1点的地形匹配定位点索引号,ni∈{1,2,3...Ni}为wpi点的水下航行器路径定位点索引号;
Figure FDA0003229219210000042
按照由小到大的顺序排序,得到序列
Figure FDA0003229219210000043
Figure FDA0003229219210000044
保存
Figure FDA0003229219210000045
对应在
Figure FDA0003229219210000046
中的索引号;
根据
Figure FDA0003229219210000047
的索引号对地形匹配导航规划点wpi+1获取的水下航行器路径定位点
Figure FDA0003229219210000048
水下航行器路径定位点的定位误差
Figure FDA0003229219210000049
水下航行器路径定位点的定位偏
Figure FDA00032292192100000410
水下航行器路径定位点的定位误差的方差下确界
Figure FDA00032292192100000411
进行排序,得到序列
Figure FDA00032292192100000412
Figure FDA00032292192100000413
根据水下航行器路径定位点的定位偏
Figure FDA00032292192100000414
和水下航行器路径定位点的定位误差
Figure FDA00032292192100000415
对水下航行器路径定位点的有效性进行筛选;
获得筛选后的水下航行器路径定位点、水下航行器路径定位点的定位误差分别为:Ddte(m)、Ddt(m),并将此序列重新赋值给
Figure FDA00032292192100000416
从而获得新的水下航行器路径定位点、水下航行器路径定位点的定位误差;
根据新的水下航行器路径定位点
Figure FDA00032292192100000417
和水下航行器路径定位点的定位误差
Figure FDA00032292192100000418
及水下航行器路径定位点的估计代价
Figure FDA00032292192100000419
选出水下航行器路径定位点的估计代价中的极小值
Figure FDA00032292192100000420
并提取其在序列
Figure FDA00032292192100000421
中对应的索引号
Figure FDA00032292192100000422
根据索引号
Figure FDA00032292192100000423
选择水下航行器路径定位点
Figure FDA00032292192100000424
中与索引号
Figure FDA00032292192100000425
对应的最优水下航行器路径定位点
Figure FDA00032292192100000426
随后执行水下航行器路径定位过程的收敛性判断。
7.根据权利要求6所述的一种新的地形匹配导航系统初始化方法,其特征在于,进行系统初始对准滤波过程模块的输出信息的状态检测和收敛性判断的方法包括:
地形匹配导航系统状态检测模块一根据获得的水下航行器路径定位点
Figure FDA00032292192100000427
和水下航行器路径定位点的定位误差
Figure FDA0003229219210000051
及水下航行器路径定位点的估计代价
Figure FDA0003229219210000052
计算水下航行器路径定位点的融合定位点
Figure FDA0003229219210000053
α是为了确保
Figure FDA0003229219210000054
成立的归一化矩阵,I为单位矩阵;
计算最优水下航行器路径定位点
Figure FDA0003229219210000055
Figure FDA0003229219210000056
之间的偏差
Figure FDA0003229219210000057
水下航行器路径定位收敛判断关系式:
Figure FDA0003229219210000058
是否成立,如果成立,则表明水下航行器路径定位满足收敛要求,地形匹配导航系统初始对准滤波过程模块的输出信息收敛,完成系统初始对准过程模块的运行,i+1时刻的最优水下航行器路径定位点
Figure FDA0003229219210000059
即为地形匹配导航输出,同时地形匹配导航系统进入地形匹配导航跟踪导航滤波模块。
CN202110981651.3A 2021-08-25 2021-08-25 一种新的地形匹配导航系统、方法 Active CN113932811B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110981651.3A CN113932811B (zh) 2021-08-25 2021-08-25 一种新的地形匹配导航系统、方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110981651.3A CN113932811B (zh) 2021-08-25 2021-08-25 一种新的地形匹配导航系统、方法

Publications (2)

Publication Number Publication Date
CN113932811A true CN113932811A (zh) 2022-01-14
CN113932811B CN113932811B (zh) 2023-10-27

Family

ID=79274432

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110981651.3A Active CN113932811B (zh) 2021-08-25 2021-08-25 一种新的地形匹配导航系统、方法

Country Status (1)

Country Link
CN (1) CN113932811B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070297289A1 (en) * 2006-06-22 2007-12-27 Northrop Grumman Corporation Underwater navigation by aided light sensor
EP3133413A2 (en) * 2015-08-17 2017-02-22 The Boeing Company Navigation system for a self-guided aerial vehicle utilizing multiple optical sensors
CN106767836A (zh) * 2017-02-17 2017-05-31 哈尔滨工程大学 一种auv地形匹配导航滤波方法
US20180080787A1 (en) * 2016-09-16 2018-03-22 The Charles Stark Draper Laboratory, Inc. Star Tracker-Aided Airborne or Spacecraft Terrestrial Landmark Navigation System
US20180313651A1 (en) * 2016-09-16 2018-11-01 The Charles Stark Draper Laboratory, Inc. Star Tracker-Aided Airborne or Spacecraft Terrestrial Landmark Navigation System
CN109900271A (zh) * 2019-02-25 2019-06-18 河海大学 地形匹配定位中有效伪波峰估计与多点融合初始定位方法
CA3094969A1 (en) * 2018-03-28 2019-10-03 Psionic Llc Navigation apparatus and method
CN110441760A (zh) * 2019-09-16 2019-11-12 河海大学 一种基于先验地形图的大范围海底地形图拓展构图方法
CN111121778A (zh) * 2019-12-03 2020-05-08 河海大学 一种导航系统初始化方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070297289A1 (en) * 2006-06-22 2007-12-27 Northrop Grumman Corporation Underwater navigation by aided light sensor
EP3133413A2 (en) * 2015-08-17 2017-02-22 The Boeing Company Navigation system for a self-guided aerial vehicle utilizing multiple optical sensors
US20180080787A1 (en) * 2016-09-16 2018-03-22 The Charles Stark Draper Laboratory, Inc. Star Tracker-Aided Airborne or Spacecraft Terrestrial Landmark Navigation System
US20180313651A1 (en) * 2016-09-16 2018-11-01 The Charles Stark Draper Laboratory, Inc. Star Tracker-Aided Airborne or Spacecraft Terrestrial Landmark Navigation System
CN106767836A (zh) * 2017-02-17 2017-05-31 哈尔滨工程大学 一种auv地形匹配导航滤波方法
CA3094969A1 (en) * 2018-03-28 2019-10-03 Psionic Llc Navigation apparatus and method
CN109900271A (zh) * 2019-02-25 2019-06-18 河海大学 地形匹配定位中有效伪波峰估计与多点融合初始定位方法
CN110441760A (zh) * 2019-09-16 2019-11-12 河海大学 一种基于先验地形图的大范围海底地形图拓展构图方法
CN111121778A (zh) * 2019-12-03 2020-05-08 河海大学 一种导航系统初始化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴康;赵龙;: "适用航空的地形匹配导航算法研究", 压电与声光, no. 05 *
王汝鹏;李晔;马腾;丛正;: "AUV地形匹配导航快速收敛滤波", 华中科技大学学报(自然科学版), no. 07 *

Also Published As

Publication number Publication date
CN113932811B (zh) 2023-10-27

Similar Documents

Publication Publication Date Title
JP6427908B2 (ja) 地図情報生成システム、方法およびプログラム
CA2870381C (en) Adaptive mapping with spatial summaries of sensor data
WO2017162036A1 (zh) 偏航识别方法、终端和存储介质
Lacaze et al. Path planning for autonomous vehicles driving over rough terrain
García et al. Augmented state Kalman filtering for AUV navigation
US20190391583A1 (en) Method and device for detecting road boundary
Luetteke et al. Implementation of the hungarian method for object tracking on a camera monitored transportation system
CN111736586B (zh) 用于路径规划的自动驾驶车辆位置的方法及其装置
JP2020123346A (ja) 各領域において最適化された自律走行を遂行できるように位置基盤アルゴリズムの選択によってシームレスパラメータ変更を遂行する方法及び装置
CN111354022B (zh) 基于核相关滤波的目标跟踪方法及系统
CN114636414A (zh) 高清晰度城市地图绘制
Wang et al. Chemical plume tracing using an auv based on pomdp source mapping and a-star path planning
Girija et al. Tracking filter and multi-sensor data fusion
Wyffels et al. Joint tracking and non-parametric shape estimation of arbitrary extended objects
CN109387198A (zh) 一种基于序贯检测的惯性/视觉里程计组合导航方法
CN112904901B (zh) 一种基于双目视觉slam与融合算法的路径规划方法
CN113932811A (zh) 一种新的地形匹配导航系统、方法
Aufrère et al. A state exchange approach in real conditions for multi-robot cooperative localization
CN111507151A (zh) 使得能够评价自动行驶客体检测过程的支援方法及装置
WO2022085368A1 (ja) 情報処理装置、情報処理システム、および方法、並びにプログラム
Tiar et al. ICP-SLAM methods implementation on a bi-steerable mobile robot
Kleiner et al. Mapping for the support of first responders in critical domains
Li et al. An enhanced direction calibration based on reinforcement learning for indoor localization system
Baltzakis et al. Closing multiple loops while mapping features in cyclic environments
Su et al. 3D Multi-Object Tracking based on Two-Stage Data Association for Collaborative Perception Scenarios

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