本发明专利技术提供的是一种AUV海底地形匹配导航的滤波方法。包括对地形匹配定位的误差估计,地形匹配定位的似然函数修正,地形匹配定位的置信区间估计,参考导航与地形匹配结果的融合滤波。本发明专利技术的主要部分是地形匹配定位的误差估计方法、地形匹配定位数据和参考导航数据的融合方法,地形匹配导航计算机接收参考导航数据和实时地形测量数据,通过匹配定位,融合滤波两个步骤得到定位位置估计,然后将位置信息送入参考导航计算机进行导航修正,实现连续的递推导航。
【技术实现步骤摘要】
一种AUV地形匹配导航滤波方法
本专利技术涉及的是一种水下地形匹配导航方法,具体地说是一种AUV的地形匹配导航滤波方法。
技术介绍
水下地形匹配导航具有无累计导航偏差的优点可作为推算导航的偏差修正,修正过程实际上是地形匹配定位结果与推算导航结果的融合滤波过程,由于地形匹配定位通常不具有连续性以及地形变化的具有很强的随机性容易导致滤波结果的发散。目前地形匹配导航方法主要有TERCOM方法,SIATN方法。TERCOM法属于批处理算法,通过获得航线下方的部分地形与载体携带的先验地形配准获得定位位置,TERCOM算法不涉及滤波和递推导航的问题,由于没有充分利用推算导航的信息,导致算法本身不具备较好的可靠性和连续导航能力;SITAN方法属于递推形式的导航方法,但递推过程需要进行地形的线性化,由于地形具有很强的非线性特征和较大的测量误差,这容易导致滤波发散。
技术实现思路
本专利技术的目的在于提供一种能实现连续的地形匹配导航和非连续的地形匹配导航,提高导航方式的灵活性的AUV地形匹配导航滤波方法。本专利技术的目的是这样实现的:AUV在水下航行lt|t距离之后到达地形匹配导航的修正点,推算导航给出了推算导航位置推算导航的误差并开启地形匹配导航程序,此时按照如下的步骤完成地形匹配导航的滤波;步骤一、向地形匹配导航系统输入初始数据,所述初始数据包括:推算导航位置推算导航误差AUV航程当前时刻测量得到的地形Z以及先验地形图h;步骤二、计算地形匹配定位的搜索区间估计,所述搜索区间由下式确定:其中:λx,λy表示导航偏差占载体的运行里程的百分比,由参导航系统的统计结果事先给出;α≥1表示放大倍数,取(1,1.5];以为中心,和为x方向和y方向的长度构成的矩形区域就是搜索区间XS;步骤三、地形匹配定位,根据步骤二中计算得到的搜索范围XS开始进行地形搜索匹配,假设潮差为tij,搜索点ij的先验地形插值序列与测量地形的高度偏差序列:假设地形高度偏差序列满足白噪声分布而潮差是常量,得到搜索点ij的潮差和地形测量误差估计:得到了潮差的估计值之后,计算测量地形和先验地形的真实相似程度:其中:h(xij,yij)表示测量序列点(xij,yij)在先验地形图h中的插值结果;zij∈Z,Z表示测量地形,zij表示测量地形内的点;查找似然函数Lij值最大的位置得到地形匹配定位的位置得到定位位置后将定位位置的潮差估计和测量地形的误差估计作为当前定位位置的最优潮差估计和测量误差估计记为:步骤四、地形匹配地位结果的有效性检验根据地形匹配定位位置检验该定位点的有效性,建立如下的有效性检验模型:其中Ct|t为有效区间,由下面的方法得到:(1)以步骤三中得到的似然函数Lij为输入,计算定位点似然函数的下界Llow;其中:σp表示地形匹配地位点的测量误差估计;mn表示地形匹配地位点的个数;Zα/2表示以σp为正态分布的上分位点对应的值;(2)以Llow为等高平面截取步骤三中的似然函数Lij,得到的等高线其中i表示等高线索引号,假设共有n条等高线;(3)查找包围地形匹配定位点的等高线在等高线所包围区域的内部,就是地形匹配定位的置信区间Ct|t;步骤五、地形匹配定位结果的定位误差估计,其估计方法如下,先计算地形匹配定位点的信息量,I+,I-表示正向和负向的信息量,先计算I+:式中:m,n表示测量地形的节点数目,k,l表示地形节点的索引号,σp表示定位点的测量误差估计,表示测量地形中的(k,l)号节点在先验地形中的插值结果;由于定位误差等于信息矩阵的逆P+=(I+)-1进一步简化:其中:由此得到了定位点在x+和y+方向的定位误差将I+计算式中的+d换成-d同理计算P-:由此得到了定位点在x-和y-方向的定位误差以及x+和y+方向的定位误差将定位误差标记为:步骤六、地形匹配定位与推算导航的融合滤波利用推算导航输出误差Pt|t和推算导航的位置地形匹配定位位置地形匹配定位误差和求得两个估计值的滤波结果:再计算滤波误差滤波误差为:由于地形的强非线性和随机性特征导致地形匹配导航的滤波难以采用常规的方法(如Kalman滤波)容易滤波发散和收敛太慢,通常滤波过程需要计算观测预测值,这就需要进行地形的线性化操作,不仅实现困难还会引入更大的误差。本专利技术就是为了简化滤波过程,实现地形匹配导航的连续或非连续递推而设计的。本专利技术提出了一种新的地形匹配导航的滤波方法。该方法可以实现连续的地形匹配导航和非连续的地形匹配导航,并且导航过程中不需要对地形做额外的处理,避免地形的线性化等其他复杂的操作,导航过程可以对地形匹配定位点进行有效性判定和推算导航与地形匹配定位结果的递推融合滤波。本专利技术滤波过程简单不需要进行额外的地形线性化处理,由于滤波过程包含了地形匹配定位结果的有效性评估具有较好的鲁棒性,不存在滤波结果发散的问题。利用本方法不仅可以实现连续的地形匹配定位与推算导航的融合滤波,而且还可以进行间断性的融合滤波,也就是说两次滤波之间可以有一段时间只使用推算导航,当获得有效的地形匹配定位位置后再进行滤波融合,提高了导航方式的灵活性。附图说明图1是地形匹配导航的滤波过程;图2是地形匹配定位过程流程图;图3是地形匹配定位点有效性检测流程图;图4是地形匹配定位与推算导航定位结果的相对位置关系。具体实施方式下面举例对本专利技术作更详细的描述。AUV在水下航行lt|t距离之后到达地形匹配导航的修正点,推算导航给出了推算导航位置推算导航的误差并开启地形匹配导航程序,此时按照如下的步骤完成地形匹配导航的滤波。步骤100:初始数据输入,向地形匹配导航系统输入初始数据:推算导航位置推算导航误差AUV航程当前时刻测量得到的地形Z;先验地形图h。1.计算地形匹配定位的搜索区间估计(结合图中101)该搜索区间即使AUV可能的位置区间,搜索区间由下式确定:式中:λx,λy表示导航偏差占载体的运行里程的百分比,由参导航系统的统计结果事先给出;α≥1表示放大倍数,一般取(1,1.5];以为中心,和为x方向和y方向的长度构成的矩形区域就是搜索区间XS。2.地形匹配定位(结合图中102)根据专利技术步骤1中计算得到的搜索范围XS(结合图中101)开始进行地形搜索匹配,假设潮差为tij,搜索点ij的先验地形插值序列(结合图中201)与测量地形(结合图中103)的高度偏差序列:假设地形高度偏差序列满足白噪声分布而潮差是常量,得到搜索点ij的潮差和地形测量误差估计(结合图中202):得到了潮差的估计值之后,计算测量地形和先验地形的真实相似程度(结合图中203):式中:h(xij,yij)表示测量序列点(xij,yij)在先验地形图h中的插值结果;zij∈Z,Z表示测量地形,zij表示测量地形内的点。查找似然函数Lij值最大的位置(结合图中206)可得到地形匹配定位的位置得到定位为之后将定位位置的潮差估计和测量地形的误差估计作为当前定位位置的最优潮差估计和测量误差估计记为:(结合图中207)。3.地形匹配地位结果的有效性检验(结合图中107)根据步骤2中获得的地形匹配定位位置下面即使检验该定位点的有效性,建立如下的有效性检验模型:(结合图中305,结合图中306)式中Ct|t为有效区间,由下面的方法得到:(1)以步骤3中得到的似然函数Lij为输本文档来自技高网...
【技术保护点】
一种AUV地形匹配导航滤波方法,其特征是:AUV在水下航行l
【技术特征摘要】
1.一种AUV地形匹配导航滤波方法,其特征是:AUV在水下航行lt|t距离之后到达地形匹配导航的修正点,推算导航给出了推算导航位置推算导航的误差并开启地形匹配导航程序,此时按照如下的步骤完成地形匹配导航的滤波;步骤一、向地形匹配导航系统输入初始数据;步骤二、计算地形匹配定位的搜索区间估计;步骤三、地形匹配定位;步骤四、地形匹配地位结果的有效性检验;步骤五、地形匹配定位结果的定位误差估计;步骤六、地形匹配定位与推算导航的融合滤波。2.根据权利要求1所述的AUV地形匹配导航滤波方法,其特征是:所述初始数据包括:推算导航位置推算导航误差AUV航程当前时刻测量得到的地形Z以及先验地形图h。3.根据权利要求3所述的AUV地形匹配导航滤波方法,其特征是:所述搜索区间由下式确定:其中:λx,λy表示导航偏差占载体的运行里程的百分比,由参导航系统的统计结果事先给出;α≥1表示放大倍数,取(1,1.5];以为中心,和为x方向和y方向的长度构成的矩形区域就是搜索区间XS。4.根据权利要求3所述的AUV地形匹配导航滤波方法,其特征是所述地形匹配定位具体包括:根据步骤二中计算得到的搜索范围XS开始进行地形搜索匹配,假设潮差为tij,搜索点ij的先验地形插值序列与测量地形的高度偏差序列:假设地形高度偏差序列满足白噪声分布而潮差是常量,得到搜索点ij的潮差和地形测量误差估计:得到了潮差的估计值之后,计算测量地形和先验地形的真实相似程度:其中:h(xij,yij)表示测量序列点(xij,yij)在先验地形图h中的插值结果;zij∈Z,Z表示测量地形,zij表示测量地形内的点;查找似然函数Lij值最大的位置得到地形匹配定位的位置得到定位位置后将定位位置的潮差估计和测量地形的误差估计作为当前定位位置的最优潮差估计和测量误差估计记为:5.根据权利要求4所述的AUV地形匹配导航滤波方法,其特征是所述地形匹配地位结果的有效性检验具体包括:根据地形匹配定位位置检验该定位点的有效性,建立如下的有效性检验模型:其中Ct|t为有效区间,由下面的方法...
【专利技术属性】
技术研发人员:李晔,王汝鹏,姜言清,马腾,丛正,贡雨森,郭宏达,安力,龚昊,何佳雨,
申请(专利权)人:哈尔滨工程大学,
类型:发明
国别省市:黑龙江,23
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。