一种重力辅助惯性导航适配区可导航性的评价方法
【技术领域】
[0001] 本发明涉及一种重力辅助惯性导航适配区可导航性的评价方法，属于重力辅助惯 性导航系统航迹规划领域。
【背景技术】
[0002] 水下载体在航行过程中常用重力来辅助惯导进行位置校正，用重力来辅助惯性导 航就是利用重力异常值或者重力梯度值等辅助信息来修正惯性导航随时间积累的定位误 差，以实现载体的精确导航。重力辅助惯性导航的核心是匹配算法，其中序列匹配算法的匹 配精度与局部重力场特征有很大的关系，局部区域重力异常数据离散程度越高，校正精度 越高；反之校正精度较差，甚至会发散，因此保证载体可以长航时、高精度航行的关键是提 取丰富的适配区使得惯导的定位误差可以得到及时的校正。
[0003]目前有很多文献已给出了提取适配区的方法，但其所选的适配区究竟能否有校辅 助惯导进行校正，并没有给出具体的评价方法。在实际航行过程中，一片海域内的适配区是 固定的，但载体进入海域的方向是不定的，当载体沿不同方向进入适配区时，参与匹配的重 力场曲面剖面也不同，得到的匹配效果也不同，在适配区内不一定在所有方向上都能得到 好的匹配效果，因此为了全面评价重力场适配区的可导航性，本发明引入重力异常三维曲 面的等方性系数来对适配区进行评价，从而有效的利用重力辅助惯导进行位置校正。
【发明内容】
[0004] 本发明的目的是为了克服已有技术存在的不足，提出一种重力辅助惯性导航适配 区可导航性的评价方法，解决了目前无法评价所选适配区是否有效的技术问题。
[0005] 本发明的基本原理是：以分形几何理论为基础，通过分析重力异常曲面的轮廓谱 距和表面谱距的特性，提出一种重力场特征度量值一一等方性系数，分别得出等方性好与 差时，适配区的匹配概率与方向之间的定性关系，有助于提高重力场适配区可导航性分析 的效率。
[0006] 本发明的目的是通过下述技术方案实现的。
[0007] -种重力辅助惯性导航适配区可导航性的评价方法，包括如下步骤：
[0008] 步骤一、根据重力场重力异常值确定适配区重力场三维曲面图，利用随机游走法 计算其三维曲面图的分形维数；
[0009] 进一步地，采用如下方法计算三维分形维数，步骤如下：
[0010] (1)首先将重力异常值看作随机游走的结果，设定一个间隔值R，计算每个点跟它 上下左右相邻R的重力异常值差，即：
[0011]
[0012] 其中P为两坐标点的重力异常值差，g(x, y)为在重力异常网格图上坐标为（X，y) 的重力异常值，I I (x2, y2)_ (X1, I I = R ;
[0013] (2)计算P的期望E (P);
[0014] (3)计算E (P)和R的双对数log-log函数，取不同R得到一组对应的E (P)和R， 进行线性拟合，得到斜率k，由D = 3-k求得重力异常三维曲面图的分形维数D。
[0015] 步骤二、计算实际重力场重力异常方差，反推求得重力异常序列特征尺度系数；
[0016] 用于重力场重力异常序列轮廓的表征，分形维数为D的W-M函数形式如下：
[0017]
(2)
[0018] g(x，y)为重力场重力异常数据；D为分形维数；G为特征尺度系数；γ为大于1的 常数，7"为随机过程的空间频率与重力异常面轮廓的最低阶段频率相对应。
[0019] 进一步地，对于近似服从正态分布的重力异常三维曲面，γ取L 5。
[0020] 相应于（2)式的连续功率谱为
[0021]
(3)
[0022] 其中ω为频率；
[0023] 根据（3)式得重力异常序列方差表示为
[0024]
(4；
[0025] 计算实际重力异常的方差得到m。，根据（4)式求得特征尺度系数G，其中
:，Λ为取样间隔，L为取样长度。
[0026] 步骤三、计算重力异常三维曲面图的表面二阶谱距；
[0027] 重力异常曲面是一个三维均匀的随机表面，γ阶轮廓谱距πν和表面谱距 义分别如下：
[0028]
[0029] (6)
[0030] 其中Ww1, ω2)为重力异常轮廓的二元频谱函数，根据Longuet-Higgins等式，设 定一个参考方向X，与参考方向成Q角方向上的γ阶轮廓谱距πν与表面谱距有如下 关系：
[0031]
[0032] 由（7)式推出在Θ方向上轮廓的二阶谱距％( Θ )为
[0033] m2( Θ ) = m20cos2 Θ +2mncos Θ sin Θ +m02sin2 θ (8)
[0034] 根据（8)式，分别在三个方向上测量重力异常数据，可以求得重力异常序列轮廓 表面的二阶谱距，根据（9)式可以求得重力异常三维曲面的表面谱距：
[0035] υ?Ν 丄 υυιυ? ?υο λ j < 〇/ υ
[0036] ;为任选的与参考方向成 一定夹角的三个方向，表面谱距m2。和m。2为相互垂直的两个方向上的斜率方差，而m n为两 个方向上的协方差；
[0037] 步骤四、根据重力异常三维曲面的表面谱距计算其等方性系数，进行适配区可导 航性的评价，
[0038] 测量坐标系旋转时，该随机表面各条轮廓高度方向的概率分布不变。轮廓的二阶 谱矩，既可以反映轮廓幅度分布，又可以反映频率特性，而且一般情况下很强地依赖于被测 轮廓的方向，因此可以用其来表征各向异性表面的等方性。如果三维表面为等方性的，则必 定m 2 ( Θ )与Θ无关，此时IH2 ( Θ )的轨迹为圆。
[0039] 由 m2 ( Θ )= m2(]COS2 Θ +2mncos Θ sin Θ +m(]2Sin2 Θ 可知轮廓的二阶谱距不仅依赖于 Θ。对于均一的随机表面，轮廓斜率方差最大值和最小值所在的方向是相互垂直的，这两个 方向称为主方向，而斜率方向最大的那条轮廓曲线与X夹角称为主方向Θ p。
[0040] 对（8)式求导求取最大值和最小值可以求得
[0043] 因此若三维表面是完全等方性的，则^ = 4Δ2，我们用丁式评定各向异性表面的 等方性
[0044]
[0045] 当mn= 0，m2Q=mQ2, Λ= 1时，表示重力异常三维曲面是完全等方性的，所以Λ越 接近1，等方性越好。实际评价时，对于同一重力异常曲面，测量的网格数不仅仅代表三条轮 廓线，这时Λ不为单一数值，等方性参数可取其平均值
[0046] 根据等方性系数较低和较高的情况下特定区域上的适配性度量值和方向之间的 关系，在进行导航任务的预处理时，可以先计算某方向上适配区的等方性系数，如果该适配 区等方性较好，则它的导航性能对方向不敏感，各个方向上的匹配概率相差很小，不需要再 经过复杂的分析就可以全面得出该区域的可导航性；如果该区域等方性较差，则它的导航 性能对方向敏感，适配区在各个方向上的可导航性能可能不一样，并且各个方向上的匹配 概率起伏很大，只有选择较小的方向粒度作为步调再计算分析该区域在各个方向上的匹配 性能指标，才能全面分析该适配区的可导航性，按照上述策略可以大大提高导航规划系统 预处理的效率。
[0047] 本发明的有益效果：
[0048] 本发明方法给出了一种重力辅助惯性导航适配区可导航性的评价方法，利用分形 几何的原理，通过分析重力异常曲面的轮廓谱距和表面谱距的特性，提出用等方性系数来 衡量适配区的可导航性，分别得出等方性好与差时，适配区的匹配概率与方向之间的定性 关系，有助于提高重力场适配区可导航性分析的效率。
【具体实施方式】
[0049] 本实施例中，取10° XKT的海洋范围，分辨率为10' X10'，局部区域的网格取 m = n = 5,也即5X5的局部区域网格，重力场背景数据为重力异常值。其过程如下：
[0050]
[0051] 步骤一、根据重力场重力异常值确定适配区重力场三