CN103743402B - 一种基于地形信息量的水下智能自适应地形匹配方法 - Google Patents
一种基于地形信息量的水下智能自适应地形匹配方法 Download PDFInfo
- Publication number
- CN103743402B CN103743402B CN201410002945.7A CN201410002945A CN103743402B CN 103743402 B CN103743402 B CN 103743402B CN 201410002945 A CN201410002945 A CN 201410002945A CN 103743402 B CN103743402 B CN 103743402B
- Authority
- CN
- China
- Prior art keywords
- sigma
- fuzzy
- equal
- terrain
- landform
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments 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)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于地形信息量的水下智能自适应地形匹配方法,其主要目的在于解决水下航行区由于地形信息量不充分或不完整导致地形辅助导航***的误差增大,以至不能满足航行器长时间精确导航定位需求的问题。本发明的主要步骤包括:航行区内不同地形特征参数值的计算、基于智能方法的地形信息量综合指标计算,基于地形信息量综合指标值的采样时间计算、自适应地形匹配算法。本发明可以有效改善航行区域地形信息量不足引起的匹配误差增大以至不能匹配的问题,以及使用单一地形特征指标会导致对地形信息量的评价不全面的问题,利用智能方法综合考虑各个指标的影响,基于地形信息量自适应调节采样时间增大匹配范围并提高匹配精度。
Description
技术领域
本发明涉及地形导航技术领域,尤其设计一种能够满足长时间高精度航行要求的基于地形信息量的水下智能自适应地形匹配方法。
背景技术
地形导航技术在导航领域具有巨大的应用前景,主要是因为它是一个自主导航***,几乎不受外界干扰。匹配算法的设计是地形导航技术的核心技术之一,算法设计好坏直接影响了导航精度。地形导航信息量是地形的一种统计特性,是以地形高程提供地形平面位置信息为基础的,对区域内所有地形高程提供地形平面位置信息能力的一种平均。在实际应用中,对于不同的地形其包含的导航信息量是完全不一样的,由于使用单一地形导航信息参数会导致对地形信息的评价不全面,需综合考虑不同参数的影响。对地形导航信息量进行综合评价一直是地形导航技术研究的一个重点,也是影响地形匹配技术广泛应用的重要因素,因此通过最大化地形导航信息利用率来提高导航精度是地形导航技术的研究目标之一。
发明内容
发明目的:为解决现有技术存在的不足,本发明提出一种基于地形信息量的水下智能自适应地形匹配方法,本发明方法针对不同地形包含的导航信息量不同,其对地形匹配算法精度的影响也不同这一特点,设计了一种基于模糊控制技术综合评价地形信息量的方法,并利用地形信息量综合评价结果设计自适应地形匹配算法。
本发明的技术方案具体如下:
一种基于地形信息量的水下智能自适应地形匹配方法,包括如下步骤:
步骤S1:适配区内地形特征参数值的计算:
通过地形高程数据库给定地形高度样本,设定航行区经度、纬度跨度为m×n网格,其中m、n为大于1的整数,h(i,j)为经度方向网格点i(i=1,2,......m)及纬度方向网格点j(j=1,2,......n)的交叉点(i,j)处的高度值,则各地形统计特征参数值的具体计算方法如下:
(1)地形粗糙度:
r=(rλ+rφ)/2
其中rλ和rφ分别为经度方向和纬度方向的粗糙度:
(2)地形标准差:
其中为地形高度均值;
(3)地形高度熵
其中pij为地形点坐标(i,j)处的归一化高程值;
(4)相关系数:
R=(Rλ+Rφ)/2
其中Rλ和Rφ分别表示经度方向和纬度方向的相关系数,
步骤S2:基于步骤S1中对航行区内导航信息量的各种特征值计算结果,建立一个特征向量X,包括上述所有地形特征参数,即X=[rσHfR],应用模糊控制方法,将该边界不清、不易定量的特征参数定量化,具体地,包括如下步骤:
(a)将航行器航行区域等范围划分成L块子航行区,其中L大于1的整数,则航行区m×n网格区变为L个ml×nl网格区,且l=1,2,...,L,其中L值越大,模拟精度也高,但计算也会越复杂,成本对应特征向量Xl=[rlσlHflRl],构建4个模糊控制器,分别计算第l个网格区内4个特征参数包含的信息量,并计算不同网格区地形信息量的权值系数;
其中,对于计算第l个网格区内地形粗糙度r包含信息量的模糊控制器构造过程如下:
模糊控制器的输入变量为rl,输出为其所包含信息量的权值系数该模糊控制器输入变量rl的语言值为Erl,其模糊子集定义为Erl={equal0,less1,equal1},其中模糊语言变量equal0,less1,equal1分别表示基本等于0,基本小于1,基本等于1;该模糊控制器输出变量的语言值为其模糊子集定义为其中模糊语言变量equal0,less1,equal1分别表示基本等于0,基本小于1,基本等于1;
根据地形匹配精度随地形粗糙度的增大而增大这一规律,选择模糊控制规则Rr如下:
该控制规则表明如果地形起伏变化不明显,则地形粗糙度rl值很小,包含地形信息量少,此时应减小其权值系数降低该地形的利用率甚至不利用;如果地形起伏变化不明显,则形粗糙度r值不大,进而包含地形信息量不大,此时其权值系数也不大,应适度利用该地形进行匹配;如果地形起伏变化很大,则地形粗糙度r很大,进而包含地形信息量增大,此时应增大其权值系数提高该地形的利用率。
输入变量为rl的隶属度函数选取正态函数: 输出变量的隶属度函数选取正态函数
利用模糊合成运算得到模糊控制器的输出模糊量为该模糊量不能直接用于计算,需将模糊量转化为清晰的数字量,采用重心法将其清晰化:
对于计算第l个网格区内地形标准差σ包含信息量的模糊控制器构造过程如下:
模糊控制器的输入变量为σl,输出为其所包含信息量的权值系数该模糊控制器输入变量σl的语言值为其模糊子集定义为其中模糊语言变量equal0,less1,equal1分别表示基本等于0,基本小于1,基本等于1;该模糊控制器输出变量的语言值为其模糊子集定义为其中模糊语言变量equal0,less1,equal1分别表示基本等于0,基本小于1,基本等于1;
根据地形匹配精度随地形标准差的增大而增大这一规律,选择模糊控制规则Rσ如下:
该控制规则表明如果地形起伏程度不大,则地形标准差σl值很小,包含地形信息量少,此时应减小其权值系数降低该地形的利用率甚至不利用;如果地形起伏程度不明显,则地形标准差σl值不大,包含地形信息量不大,此时其权值系数也不大,应适度利用该地形进行匹配;如果地形起伏程度很大,则地形标准差σl很大,包含地形信息量增大,此时应增大其权值系数提高该地形的利用率。
输入变量为σl的隶属度函数选取正态函数: 输出变量的隶属度函数选取正态函数
利用模糊合成运算得到模糊控制器的输出模糊量为该模糊量不能直接用于计算,需将模糊量转化为清晰的数字量,采用重心法将其清晰化:
对于计算第l个网格区内地形高度熵Hf包含信息量的模糊控制器构造过程如下:
模糊控制器的输入变量为Hfl,输出为其所包含信息量的权值系数该模糊控制器输入变量Hfl的语言值为其模糊子集定义为其中模糊语言变量equal0,more0,less1,equal1分别表示基本等于0,基本大于0,基本小于1,基本等于1;该模糊控制器输出变量的语言值为其模糊子集定义为 其中模糊语言变量equal0,more0,less1,equal1分别表示基本等于0,基本大于0,基本小于1,基本等于1;
根据地形匹配精度随地形高度熵的增大而减小这一规律,选择模糊控制规则如下:
该控制规则表明如果地形高度变化越剧烈程度很大,则地形高度熵Hfl值很小,包含地形信息量很大,此时应增大其权值系数最大限度的提高该地形的利用率;如果地形高度变化越剧烈程度相对比较大,则形高度熵Hfl的值相对较小,包含地形信息量相对比较多,此时其权值系数也相对较大,应适度利用该地形进行匹配;如果地形高度变化越剧烈程度相对不太大,则地形高度熵Hfl的值相对较大,包含地形信息量相对比较少,此时其权值系数也相对较小,应降低该地形的利用率;如果地形高度变化越剧烈程度很小或几乎没有,则地形高度熵Hfl的值很大,包含地形信息量相对比很少,此时其权值系数也很小,应尽量减少利用该地形进行匹配。
输入变量为Hfl的隶属度函数选取正态函数: 输出变量的隶属度函数选取正态函数
利用模糊合成运算得到模糊控制器的输出模糊量为该模糊量不能直接用于计算,需将模糊量转化为清晰的数字量,采用重心法将其清晰化:
对于计算第l个网格区内地形相关系数R包含信息量的模糊控制器构造过程如下:
模糊控制器的输入变量为Rl,输出为其所包含信息量的权值系数该模糊控制器输入变量Rl的语言值为其模糊子集定义为其中模糊语言变量equal0,more0,less1,equal1分别表示基本等于0,基本大于0,基本小于1,基本等于1;该模糊控制器输出变量的语言值为其模糊子集定义为其中模糊语言变量equal0,more0,less1,equal1分别表示基本等于0,基本大于0,基本小于1,基本等于1;
根据地形匹配精度随地形相关系数值的增大而减小这一规律,选择模糊控制规则RR如下:
该控制规则表明如果地形相关程度很大,则地形相关系数Rl值很小,包含地形信息量很大,此时应增大其权值系数最大限度的提高该地形的利用率;如果地形相关程度相对比较大,则形相关系数Rl的值相对较小,包含地形信息量相对比较多,此时其权值系数也相对较大,应适度利用该地形进行匹配;如果地形相关程度相对不太大,则地形相关系数Rl的值相对较大,包含地形信息量相对比较少,此时其权值系数也相对较小,应降低该地形的利用率;如果地形高度变化越剧烈程度很小或几乎没有,则地形相关系数Rl的值很大,包含地形信息量相对比很少,此时其权值系数也很小,应尽量减少利用该地形进行匹配。
输入变量为Rl的隶属度函数选取正态函数: 输出变量的隶属度函数选取正态函数
利用模糊合成运算得到模糊控制器的输出模糊量为该模糊量不能直接用于计算,需将模糊量转化为清晰的数字量,采用重心法将其清晰化:
(b)根据最大熵原理计算每个特征参量的权值Wr、Wσ、和WR,即各个特征参量在特征向量集中所占比重:
令:
则:
(c)然后计算表示每个网格区地形信息含量的权系数Wl:
步骤S3:计算l个网格区的匹配采样时间:根据划分的每个子网格区内地形信息含量的权系数Wl计算采样时间tl,其中T为常用典型采样时间间隔,一般在1~50s的范围内,更为常见的为20s或30s,根据航行距离和航行速度来确定;
步骤S4:水下航行器驶入地形数据库包含区域时,进行自适应地形匹配:在主惯性导航***INS的导航下航行一段距离,与此同时每隔一个时间段tl,通过航行器的高程测量装置获得实测高程Dsp,sp=1,2…10,并获得与此对应的惯导***指示航迹序列 sp为一次匹配需要的采样点个数,采样点个数主要是根据导航精度和试试性的需要进行选择,采样点数越多,则精度越高越可靠,但计算量大、实时性小,而采样越小则可能产生误匹配或匹配失效,根据多次试验选可选为10;提取高程为实测高程Dsp的等高线Csp;通过位置序列Psp向等高线Csp作垂线得到一个匹配垂足序列Ysp,通过随机旋转和平移位置序列Psp逐渐逼近匹配垂足序列Ysp,获得匹配位置序列P′sp。
有益效果:目前常用的地形匹配算法ICCP算法采样时间都是等长的,忽略了航行区导航信息量不同这一特点,本发明对传统地形匹配算法进行了改进,利用航行区的综合地形导航信息含量对采样时间进行自适应调节,克服了原有方法采样时间单一,不能充分高效地利用航行区内的地形导航信息量的缺点,大大提高了导航效率和导航定位精度,同时,本发明利用模糊控制方法对各地形参量进行了综合评定,应用模糊数学的原理,将边界不清、不易定量的地形导航信息参数定量化,给出综合信息参量值作为地形匹配采样时间权系数,克服了使用单一地形导航信息参数会导致对信息的评价不全面的缺陷。通过将定性指标定量化,对各单一信息参量进行综合评价,给出一个综合信息参量值作为地形匹配采样时间权系数,当航行区域包含的导航信息量比较大时,缩短采样时间,尽可能多的利用该处地形信息;当航行区域包含的导航信息量不丰富时,则增加采样时间间隔,减少该地形的利用率以保证匹配精度。
附图说明
图1为本发明所描述的智能自适应地形匹配算法的结构图。
图2为地形粗糙度信息量模糊计算结构框图。
图3为实施例中所采用的地形图。
图4为实施例中传统ICCP算法和本发明自适应匹配算法匹配结果图。
具体实施方式
下面结合具体实施例和附图,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域普通技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
实施例1一种基于地形信息量的水下智能自适应地形匹配方法。
图1为本发明所描述的智能自适应地形匹配算法结构图。
仿真实验中地形高程数据库为10000m×14000m范围内的矩形区域,数据库地形的分辨率为100m×100m,将该区域划分为9个网格区,如图3所示,计算9个网格区的地形信息量权系数,表1为9个网格区各单一地形信息量值及综合信息量权系数及采样时间间隔计算结果。
其中,地形粗糙度r的计算公式为:
r=(rλ+rφ)/2
其中rλ和rφ分别为经度方向和纬度方向的粗糙度:
地形高度标准差σ的计算公式为:
其中为地形高度均值;
地形高度熵Hf的计算公式为:
其中pij为地形点坐标(i,j)处的归一化高程值;
地形的相关系数R的计算公式为R=(Rλ+Rφ)/2
其中Rλ和Rφ分别表示经度方向和纬度方向的相关系数,
建立一个特征向量X,包括上述所有地形特征参数,即X=[rσHfR],应用模糊控制方法,将该边界不清、不易定量的特征参数定量化,具体地,利用如下方法
(a)将航行器航行区域等范围划分成L块子航行区,则航行区m×n网格区变为L个ml×nl网格区,且l=1,2,...,L,对应特征向量Xl=[rlσlHflRl],构建4个模糊控制器,分别计算第l个网格区内4个特征参数包含的信息量,并计算不同网格区地形信息量的权值系数;
其中,对于计算第l个网格区内地形粗糙度r包含信息量的模糊控制器构造过程如下:
模糊控制器的输入变量为rl,输出为其所包含信息量的权值系数该模糊控制器输入变量rl的语言值为Erl,其模糊子集定义为Erl={equal0,less1,equal1},其中模糊语言变量equal0,less1,equal1分别表示基本等于0,基本小于1,基本等于1;该模糊控制器输出变量的语言值为其模糊子集定义为其中模糊语言变量equal0,less1,equal1分别表示基本等于0,基本小于1,基本等于1;
根据地形匹配精度随地形粗糙度的增大而增大这一规律,选择模糊控制规则Rr如下:
输入变量为rl的隶属度函数选取正态函数: 输出变量的隶属度函数选取正态函数
利用模糊合成运算得到模糊控制器的输出模糊量为该模糊量不能直接用于计算,需将模糊量转化为清晰的数字量,采用重心法将其清晰化:
对于计算第l个网格区内地形标准差σ包含信息量的模糊控制器构造过程如下:
模糊控制器的输入变量为σl,输出为其所包含信息量的权值系数该模糊控制器输入变量σl的语言值为其模糊子集定义为其中模糊语言变量equal0,less1,equal1分别表示基本等于0,基本小于1,基本等于1;该模糊控制器输出变量的语言值为其模糊子集定义为其中模糊语言变量equal0,less1,equal1分别表示基本等于0,基本小于1,基本等于1;
根据地形匹配精度随地形标准差的增大而增大这一规律,选择模糊控制规则Rσ如下:
输入变量为σl的隶属度函数选取正态函数: 输出变量的隶属度函数选取正态函数
利用模糊合成运算得到模糊控制器的输出模糊量为该模糊量不能直接用于计算,需将模糊量转化为清晰的数字量,采用重心法将其清晰化:
对于计算第l个网格区内地形高度熵Hf包含信息量的模糊控制器构造过程如下:
模糊控制器的输入变量为Hfl,输出为其所包含信息量的权值系数该模糊控制器输入变量Hfl的语言值为其模糊子集定义为其中模糊语言变量equal0,more0,less1,equal1分别表示基本等于0,基本大于0,基本小于1,基本等于1;该模糊控制器输出变量的语言值为其模糊子集定义为其中模糊语言变量equal0,more0,less1,equal1分别表示基本等于0,基本大于0,基本小于1,基本等于1;
根据地形匹配精度随地形高度熵的增大而减小这一规律,选择模糊控制规则如下:
输入变量为Hfl的隶属度函数选取正态函数: 输出变量的隶属度函数选取正态函数
利用模糊合成运算得到模糊控制器的输出模糊量为该模糊量不能直接用于计算,需将模糊量转化为清晰的数字量,采用重心法将其清晰化:
对于计算第l个网格区内地形相关系数R包含信息量的模糊控制器构造过程如下:
模糊控制器的输入变量为Rl,输出为其所包含信息量的权值系数该模糊控制器输入变量Rl的语言值为其模糊子集定义为其中模糊语言变量equal0,more0,less1,equal1分别表示基本等于0,基本大于0,基本小于1,基本等于1;该模糊控制器输出变量的语言值为其模糊子集定义为其中模糊语言变量equal0,more0,less1,equal1分别表示基本等于0,基本大于0,基本小于1,基本等于1;
根据地形匹配精度随地形相关系数值的增大而减小这一规律,选择模糊控制规则RR如下:
输入变量为Rl的隶属度函数选取正态函数: 输出变量的隶属度函数选取正态函数
利用模糊合成运算得到模糊控制器的输出模糊量为该模糊量不能直接用于计算,需将模糊量转化为清晰的数字量,采用重心法将其清晰化:
(b)根据最大熵原理计算每个特征参量的权值Wr、Wσ、和WR,即各个特征参量在特征向量集中所占比重:
令:
则:
(c)然后计算表示每个网格区地形信息含量的权系数Wl:
然后根据划分的每个子网格区内地形信息含量的权系数Wl计算采样时间tl,其中T为常用典型采样时间间隔20s。
整个过程如图2所示,即先计算控制变量,然后进行模糊量化处理,之后构建模糊控制规则,进行模糊推理,最后进行非模糊量化处理。
初始航向为38.7°,航向变化范围为±7°,初始航速为11.5节,航速变换范围为±1.3节。水下航行器惯导***的陀螺仪漂移为0.02°/h,加速度计零偏偏置为0.0005m/s2,测深装置的测量误差取为幅值为0.2的白噪声,仿真结果如图4所示,其中实线条为真实航行轨迹,虚线条为传统相关匹配算法匹配结果,点划线为本发明采用的智能自适应匹配算法匹配结果。
表1
从仿真结果可以看出,与传统的ICCP地形匹配算法相比本发明采用的智能自适应匹配算法的匹配结果精度得到很大提高,这是由于该方法在地形变化明显的区域缩短采用时间间隔,增加该区域的地形利用率;在地形变化较小的区域增大采样时间间隔,在保证精度的同时充分利用该区域的地形信息,克服了原有方法采样时间单一,不能充分高效地利用航行区内的地形导航信息量的缺点,大大提高了相同硬件环境下的匹配效率,因此,本发明可以解决不同地形条件下单一采样时间带来匹配效率低的问题,利用智能自适应的方法提高了导航定位的精度。
Claims (1)
1.一种基于地形信息量的水下智能自适应地形匹配方法,其特征在于,包括如下步骤:
步骤S1:适配区内地形特征参数值的计算:
通过地形高程数据库给定地形高度样本,设定航行区经度、纬度跨度为m×n网格,其中m、n为大于1的整数,h(i,j)为经度方向网格点i(i=1,2,......m)及纬度方向网格点j(j=1,2,......n)的交叉点(i,j)处的高度值,则各地形统计特征参数值的具体计算方法如下:
(1)地形粗糙度:
r=(rλ+rφ)/2
其中rλ和rφ分别为经度方向和纬度方向的粗糙度:
(2)地形标准差:
其中为地形高度均值;
(3)地形高度熵
其中pij为地形点坐标(i,j)处的归一化高程值;
(4)相关系数:
R=(Rλ+Rφ)/2
其中Rλ和Rφ分别表示经度方向和纬度方向的相关系数,
步骤S2:基于步骤S1中对航行区内导航信息量的各种特征值计算结果,建立一个特征向量X,包括上述所有地形特征参数,即X=[rσHfR],应用模糊控制方法,将该边界不清、不易定量的特征参数定量化,具体地,包括如下步骤:
(a)将航行器航行区域等范围划分成L块子航行区,其中L大于1的整数,则航行区m×n网格区变为L个ml×nl网格区,且l=1,2,...,L,对应特征向量Xl=[rlσlHflRl],构建4个模糊控制器,分别计算第l个网格区内4个特征参数包含的信息量,并计算不同网格区地形信息量的权值系数;
其中,对于计算第l个网格区内地形粗糙度r包含信息量的模糊控制器构造过程如下:
模糊控制器的输入变量为rl,输出为其所包含信息量的权值系数该模糊控制器输入变量rl的语言值为Erl,其模糊子集定义为Erl={equal0,less1,equal1},其中模糊语言变量equal0,less1,equal1分别表示等于0,小于1,等于1;该模糊控制器输出变量的语言值为其模糊子集定义为其中模糊语言变量equal0,less1,equal1分别表示等于0,小于1,等于1;
根据地形匹配精度随地形粗糙度的增大而增大这一规律,选择模糊控制规则Rr如下:
输入变量为rl的隶属度函数选取正态函数:输出变量的隶属度函数选取正态函数
利用模糊合成运算得到模糊控制器的输出模糊量为该模糊量不能直接用于计算,需将模糊量转化为清晰的数字量,采用重心法将其清晰化:
对于计算第l个网格区内地形标准差σ包含信息量的模糊控制器构造过程如下:
模糊控制器的输入变量为σl,输出为其所包含信息量的权值系数该模糊控制器输入变量σl的语言值为其模糊子集定义为其中模糊语言变量equal0,less1,equal1分别表示等于0,小于1,等于1;该模糊控制器输出变量的语言值为其模糊子集定义为其中模糊语言变量equal0,less1,equal1分别表示等于0,小于1,等于1;
根据地形匹配精度随地形标准差的增大而增大这一规律,选择模糊控制规则Rσ如下:
输入变量为σl的隶属度函数选取正态函数:输出变量的隶属度函数选取正态函数
利用模糊合成运算得到模糊控制器的输出模糊量为该模糊量不能直接用于计算,需将模糊量转化为清晰的数字量,采用重心法将其清晰化:
对于计算第l个网格区内地形高度熵Hf包含信息量的模糊控制器构造过程如下:
模糊控制器的输入变量为Hfl,输出为其所包含信息量的权值系数该模糊控制器输入变量Hfl的语言值为其模糊子集定义为其中模糊语言变量equal0,more0,less1,equal1分别表示等于0,大于0,小于1,等于1;该模糊控制器输出变量的语言值为其模糊子集定义为其中模糊语言变量equal0,more0,less1,equal1分别表示等于0,大于0,小于1,等于1;
根据地形匹配精度随地形高度熵的增大而减小这一规律,选择模糊控制规则如下:
输入变量为Hfl的隶属度函数选取正态函数:输出变量的隶属度函数选取正态函数
利用模糊合成运算得到模糊控制器的输出模糊量为该模糊量不能直接用于计算,需将模糊量转化为清晰的数字量,采用重心法将其清晰化:
对于计算第l个网格区内地形相关系数R包含信息量的模糊控制器构造过程如下:
模糊控制器的输入变量为Rl,输出为其所包含信息量的权值系数该模糊控制器输入变量Rl的语言值为其模糊子集定义为其中模糊语言变量equal0,more0,less1,equal1分别表示等于0,大于0,小于1,等于1;该模糊控制器输出变量的语言值为其模糊子集定义为其中模糊语言变量equal0,more0,less1,equal1分别表示等于0,大于0,小于1,等于1;
根据地形匹配精度随地形相关系数值的增大而减小这一规律,选择模糊控制规则RR如下:
输入变量为Rl的隶属度函数选取正态函数:输出变量的隶属度函数选取正态函数
利用模糊合成运算得到模糊控制器的输出模糊量为该模糊量不能直接用于计算,需将模糊量转化为清晰的数字量,采用重心法将其清晰化:
(b)根据最大熵原理计算每个特征参量的权值Wr、Wσ、和WR,即各个特征参量在特征向量集中所占比重:
令:
则:
(c)然后计算表示每个网格区地形信息含量的权系数Wl:
步骤S3:计算l个网格区的匹配采样时间:根据划分的每个子网格区内地形信息含量的权系数Wl计算采样时间tl,其中T为常用典型采样时间间隔;
步骤S4:水下航行器驶入地形数据库包含区域时,进行自适应地形匹配:在主惯性导航***INS的导航下航行一段距离,与此同时每隔一个时间段tl,通过航行器的高程测量装置获得实测高程Dsp,sp=1,2…10,并获得与此对应的惯导***指示航迹序列sp=1,2…10,sp为一次匹配需要的采样点个数;提取高程为实测高程Dsp的等高线Csp;通过位置序列Psp向等高线Csp作垂线得到一个匹配垂足序列Ysp,通过随机旋转和平移位置序列Psp逐渐逼近匹配垂足序列Ysp,获得匹配位置序列P′sp。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410002945.7A CN103743402B (zh) | 2014-01-03 | 2014-01-03 | 一种基于地形信息量的水下智能自适应地形匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410002945.7A CN103743402B (zh) | 2014-01-03 | 2014-01-03 | 一种基于地形信息量的水下智能自适应地形匹配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103743402A CN103743402A (zh) | 2014-04-23 |
CN103743402B true CN103743402B (zh) | 2016-06-08 |
Family
ID=50500445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410002945.7A Expired - Fee Related CN103743402B (zh) | 2014-01-03 | 2014-01-03 | 一种基于地形信息量的水下智能自适应地形匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103743402B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104484891A (zh) * | 2014-10-22 | 2015-04-01 | 哈尔滨工程大学 | 一种基于纹理特征和地形特征参数的水下地形匹配方法 |
CN105243276B (zh) * | 2015-10-14 | 2017-03-22 | 中国地震局地壳应力研究所 | 一种建筑物震害分析方法 |
CN105953799B (zh) * | 2016-04-22 | 2018-12-04 | 北京理工大学 | 基于熵值法的重力场适配区内水下运载体航迹规划方法 |
CN106815397A (zh) * | 2016-12-08 | 2017-06-09 | 中国科学院新疆生态与地理研究所 | 利用等水位线进行水文地质参数分区的方法及装置 |
CN106767834B (zh) * | 2017-01-24 | 2019-11-08 | 哈尔滨工程大学 | 一种基于模糊熵值的auv水下地形匹配适配区划分方法 |
CN107727096B (zh) * | 2017-09-15 | 2020-10-27 | 哈尔滨工程大学 | 基于有效节点筛选的auv地形匹配定位方法 |
CN108871376A (zh) * | 2018-04-28 | 2018-11-23 | 南京信息工程大学 | 一种基于模糊综合评价的水下导航特征适配性分析方法 |
CN108985327B (zh) * | 2018-06-05 | 2022-02-11 | 东南大学 | 一种基于因子分析的地形匹配区自组织优化分类方法 |
CN109522903B (zh) * | 2018-08-31 | 2022-10-25 | 中国人民解放军海军大连舰艇学院 | 基于双向滚动球变换的tin_ddm地形特征线提取方法 |
CN110617819B (zh) * | 2019-10-17 | 2022-09-30 | 国营芜湖机械厂 | 一种基于蚁群算法路径规划的无人机地形辅助导航方法 |
CN111649763B (zh) * | 2020-04-26 | 2021-08-27 | 中国人民解放军61540部队 | 一种基于重力灯塔建立的潜艇导航方法及*** |
CN114002683B (zh) * | 2021-11-01 | 2024-05-28 | 长沙有色冶金设计研究院有限公司 | 防结晶水下地形探测方法及*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101339036A (zh) * | 2008-08-20 | 2009-01-07 | 北京航空航天大学 | 地形辅助导航方法和设备 |
CN102221363A (zh) * | 2011-04-12 | 2011-10-19 | 东南大学 | 一种水下潜器用捷联惯性组合导航***容错组合方法 |
-
2014
- 2014-01-03 CN CN201410002945.7A patent/CN103743402B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101339036A (zh) * | 2008-08-20 | 2009-01-07 | 北京航空航天大学 | 地形辅助导航方法和设备 |
CN102221363A (zh) * | 2011-04-12 | 2011-10-19 | 东南大学 | 一种水下潜器用捷联惯性组合导航***容错组合方法 |
Non-Patent Citations (4)
Title |
---|
Goal-Enhancement Technique for Fuzzy Terrain-Based Navigation;Wong Kit Yee etal;《Proceedings of 2012 4th International Conference on Computer Engineering and Technology(ICCET 2012)》;20120512;第159-163 * |
基于模糊决策的地形辅助导航区域选择准则;张涛等;《大连海事大学学报》;20090228;第35卷(第1期);第5-8页 * |
海底地形匹配导航技术研究;刘勇;《中国优秀硕士学位论文全文数据库 信息科技辑》;20090615;全文 * |
粒子群优化算法研究及其在海底地形辅助导航中的应用;谭佳琳;《中国博士学位论文全文数据库 工程科技II辑》;20110715;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103743402A (zh) | 2014-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103743402B (zh) | 一种基于地形信息量的水下智能自适应地形匹配方法 | |
Masseran et al. | Fitting a mixture of von Mises distributions in order to model data on wind direction in Peninsular Malaysia | |
CN104197945B (zh) | 一种基于低采样率浮动车数据的全局投票地图匹配方法 | |
US10439594B2 (en) | Actually-measured marine environment data assimilation method based on sequence recursive filtering three-dimensional variation | |
CN111695299B (zh) | 一种中尺度涡轨迹预测方法 | |
CN110608738B (zh) | 一种无人船全局气象航线动态规划方法和*** | |
CN111199270B (zh) | 一种基于深度学习的区域波高预报方法及终端 | |
CN110232471B (zh) | 一种降水传感网节点布局优化方法及装置 | |
CN106643714A (zh) | 一种自主实时机载地形辅助惯性导航方法和*** | |
CN107193060B (zh) | 一种多路径台风风暴潮快速预测方法及*** | |
CN104915534A (zh) | 基于序列学习的电力铁塔变形分析与决策方法 | |
CN114445634A (zh) | 一种基于深度学习模型的海浪波高预测方法及*** | |
CN112946784B (zh) | 一种基于深度学习的SuperDARN雷达对流图短期预报方法 | |
CN106156267A (zh) | 一种面向高速公路的车道级增强型矢量数字地图制作方法 | |
CN102999667A (zh) | 一种基于模糊回归的机场噪声预测方法及其*** | |
CN114254767A (zh) | 基于Stacking集成学习的气象水文特征预测方法及*** | |
CN105953799A (zh) | 基于熵值法的重力场适配区内水下运载体航迹规划方法 | |
CN110457795A (zh) | 基于神经网络的高层建筑物风致位移状态评估方法 | |
CN114139566A (zh) | 基于机器学***均融合特征提取提高测高精度的方法 | |
CN114742206A (zh) | 综合多时空尺度多普勒雷达数据的降雨强度估计方法 | |
CN107292039A (zh) | 一种基于小波聚类的uuv巡岸轮廓构建方法 | |
CN113962426A (zh) | 一种近海水上通航安全智能预报方法和装置 | |
CN116523152B (zh) | 漂浮物的漂移预测方法、装置、设备和介质 | |
CN111896002B (zh) | 地形辅助导航匹配区在线规划与适配性分析方法及*** | |
Yang et al. | The research on real-time map-matching algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160608 Termination date: 20190103 |
|
CF01 | Termination of patent right due to non-payment of annual fee |