CN105021182A - 一种重力辅助惯性导航适配区的选择方法 - Google Patents

一种重力辅助惯性导航适配区的选择方法 Download PDF

Info

Publication number
CN105021182A
CN105021182A CN201510300526.6A CN201510300526A CN105021182A CN 105021182 A CN105021182 A CN 105021182A CN 201510300526 A CN201510300526 A CN 201510300526A CN 105021182 A CN105021182 A CN 105021182A
Authority
CN
China
Prior art keywords
coefficient
gravity
theta
gravity anomaly
omega
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
CN201510300526.6A
Other languages
English (en)
Other versions
CN105021182B (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.)
Beijing techlink intelligent Polytron Technologies Inc
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201510300526.6A priority Critical patent/CN105021182B/zh
Publication of CN105021182A publication Critical patent/CN105021182A/zh
Application granted granted Critical
Publication of CN105021182B publication Critical patent/CN105021182B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/005Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 with correlation of navigation data from several sources, e.g. map or contour matching

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明提出一种重力辅助惯性导航适配区的选择方法,解决了传统的适配区选取方法忽略了适配区方向性的技术问题。本发明通过将载体的航行区域分为若干个子区域,在每个子区域上对重力异常数据进行分类统计分析,分别得出各个子区域的重力异常三维曲面等方性系数、重力异常变异系数、重力异常偏态系数、重力异常峰态系数,根据各个参数与重力异常概率分布的关系,构建综合特征参数表达式,从而对适配区进行选择。

Description

一种重力辅助惯性导航适配区的选择方法
技术领域
本发明涉及一种重力辅助惯性导航适配区的选择方法,属于重力辅助惯性导航系统航迹规划领域。
背景技术
水下载体最常用的导航方式是惯性导航,惯导本质上是利用加速度信息不断积分求取速度、位置信息,其定位误差会随时间积累,目前常用重力场来辅助惯导进行位置校正。由于重力场是一个缓变的随机物理场,只有在变化显著的区域才能有效的辅助惯导进行位置校正,在重力变化不显著的区域匹配效果较差,甚至会出现发散的现象,因此需要在载体航行区域内选出重力变化明显的区域作为适配区,在适配区以外区域的航线上只用惯导系统进行导航,在适配区内的航线上则用重力导航系统来校正惯导位置。
目前已经有SPSS分析法、参数阈值法、层次分析法、拐点判据法等适配区选择方法,SPSS分析法通过对重力数据进行回归分析,得到适配区的回归函数关系式,以此作为判断准则对重力匹配区域进行选择;参数阈值法以重力场标准差和经纬度相关系数作为适配区选择的数量指标,给出了重力适配区经验选择准则;层次分析法通过对地球重力场的各种特征参数进行归一化处理并且分类统计分析的基础上运用层次分析法将重力场的局部标准差、经度方向粗糙度、纬度方向粗糙度、经度方向坡度,相关系数进行组合,给出重力场适配区选择准则;拐点判据法通过运用拐点选择和拐点扩充机制,对重力场特征曲线拐点进行选择和扩充,得到重力有效匹配区域。
在实际导航过程中,载体由于起点和终点的不同进入适配区的位置及航行 路线也不同,若载体在适配区内途径的航线无法获得理想的匹配效果,则适配区的选择就失去了意义,目前所给出的适配区选择方案中并没有给出一种结合载体航线的适配区选择方法。
发明内容
本发明的目的是为了克服已有技术存在的不足,提出一种重力辅助惯性导航适配区的选择方法,解决了传统的适配区选取方法忽略了适配区方向性的技术问题。
本发明的基本原理是:通过将载体的航行区域分为若干个子区域,在每个子区域上对重力异常数据进行分类统计分析,分别得出各个子区域的重力异常三维曲面等方性系数、重力异常变异系数、重力异常偏态系数、重力异常峰态系数,根据各个参数与重力异常概率分布的关系,构建综合特征参数表达式,从而对适配区进行选择。
本发明的目的是通过下述技术方案实现的。
一种重力辅助惯性导航适配区的选择方法,其特征在于:包括如下步骤:
步骤一、将航行区域分为若干个子区域,分析每个子区域的等方性;
进一步地,对于重力异常序列三维曲面轮廓的表征,分形维数为D的W-M函数形式如下:
g ( x , y ) = G D - 1 Σ n = n 1 ∞ cos 2 π γ n x γ ( 2 - D ) n - - - ( 1 )
g(x,y)为重力异常数据;D为分形维数;G为特征尺度系数;γ为大于1的常数,γn为随机过程的空间频率;n1与地形表面轮廓的最低截断频率相对应。
进一步地,对于近似服从正态分布的重力异常三维曲面,γ取1.5。
(1)式的连续功率谱函数为
P ( &omega; ) = G 2 ( D - 1 ) 2 ln &gamma; &CenterDot; 1 &omega; ( 5 - 2 D ) , &gamma; n 1 < &omega; < + &infin; - - - ( 2 )
其中ω为频率,重力异常曲面是一个三维均匀的随机表面,γ阶轮廓谱距mr和表面谱距mγ-q,q(q=0,1…γ)定义分别如下:
m &gamma; = &Integral; - &infin; + &infin; &omega; &gamma; P ( &omega; ) d&omega; - - - ( 3 )
m &gamma; - q , q = &Integral; - &infin; + &infin; &Integral; - &infin; + &infin; &omega; 1 &gamma; - q &omega; 2 q Q ( &omega; 1 , &omega; 2 ) d &omega; 1 d &omega; 2 - - - ( 4 )
其中Q(ω12)为重力异常轮廓的二元频谱函数,根据Longuet-Higgins等式,设定一个参考方向x,与参考方向成θ角方向上的γ阶轮廓谱距mr与表面谱距mr-q,q有如下关系:
m &gamma; ( &theta; ) = m &gamma; - q , q C r q cos r - q &theta; sin q &theta; - - - ( 5 )
由(5)式推出在θ方向上轮廓的二阶谱距m2(θ)为
m2(θ)=m20cos2θ+2m11cosθsinθ+m02sin2θ    (6) 
根据(6)式可得,分别在三个方向上测量重力异常数据,求得重力异常序列轮廓表面的二阶谱距,通过(7)式即求得重力异常序列的表面谱距;
m 20 m 11 m 02 = T 2 - 1 m 2 ( &theta; 1 ) m 2 ( &theta; 2 ) m 2 ( &theta; 3 ) - - - ( 7 )
其中 T 2 = cos 2 &theta; 1 2 cos &theta; 1 sin &theta; 1 sin 2 &theta; 1 cos 2 &theta; 2 2 cos &theta; 2 sin &theta; 2 sin 2 &theta; 2 cos 2 &theta; 2 2 cos &theta; 3 sin &theta; 3 sin 2 &theta; 3 , θ1,θ2,θ3为任选的与参考方向成一定夹角的三个方向,表面谱距m20和m02为相互垂直的两个方向上的斜率方差,而m11为两个方向上的协方差;
采用下式评定各向异性表面的等方性系数△:
&Delta; = 2 &Delta; 2 M 2 = 2 m 20 &CenterDot; m 02 - m 11 2 m 20 + m 02 - - - ( 8 )
其中当m11=0,m20=m02,△=1时,表示重力异常三维曲面是完全等方性的,所以△越接近1,等方性越好。
步骤二、根据每个子区域的等方性,分析每个子区域的重力异常变异系数CV
Cv = G . S . D G . mean &times; 100 % - - - ( 9 )
其中G.S.D为重力异常标准差,G.mean为重力异常平均值;
步骤三、分析每个子区域的重力异常偏态系数Cs和峰态系数Ce
首先分析子区域的重力异常序列,得到归中后的重力异常值
&Delta;g i , j * = &Delta; g i , j - 1 NM &Sigma; i = 1 N &Sigma; j = 1 M &Delta;g i , j - - - ( 10 )
其中△gi,j为重力异常值,M,N为在重力异常网格图上所选的M*N的子区域。用归中后的重力异常值作为起始观测数据,得到均方差C0
C 0 = 1 NM &Sigma; i = 1 N &Sigma; j = 1 M &Delta;g i , j * 2 - - - ( 11 )
偏态系数Cs绝对值越大,表明重力异常序列在子区域内的偏斜程度越大;
C s = N ( N - 1 ) ( N - 2 ) &CenterDot; M ( M - 1 ) ( M - 2 ) &CenterDot; &Sigma; i = 1 N &Sigma; j = 1 M &Delta;g i , j * 3 C 0 3 / 2 - - - ( 12 )
峰态系数Ce描述的是重力异常序列分布形态的陡缓程度:
C e = N 2 - 2 N + 3 ( N - 1 ) ( N - 2 ) ( N - 3 ) &CenterDot; M 2 - 2 M + 3 ( M - 1 ) ( M - 2 ) ( M - 3 ) &CenterDot; &Sigma; i = 1 N &Sigma; j = 1 M &Delta;g i , j * 4 C 0 2 - 3 ( 2 N - 3 ) N ( N - 1 ) ( N - 2 ) ( N - 3 ) &CenterDot; 3 ( 2 M - 3 ) M ( M - 1 ) ( M - 2 ) ( M - 3 ) &CenterDot; [ &Sigma; i = 1 N &Sigma; j = 1 M &Delta;g i , j * 2 ] 2 c 0 2 - - - ( 13 )
步骤四、结合等方性系数△、变异系数CV、偏态系数Cs和峰态系数Ce构建综合特征参数,从而对适配区进行选取。
进一步地,选取的原则为:在以上四个参数中等方性是测量坐标系旋转时,该随机表面各条轮廓高度方向的概率分布不变的特征,等方性系数越大的区域越适合做匹配区;重力异常变异系数是衡量重力异常数据离散程度的参数,变异系数越大的区域越适合做匹配区;重力异常偏态系数是衡量重力异常数据不对称性质的参数,偏态系数越大,表明该区域数据的偏斜越严重,只有在某个方向上能获得较好的匹配效果,因此偏态系数大的区域不适合作为适配区;重力异常峰态系数是衡量重力异常聚中性的性质,重力异常峰态系数越大,则重力异常数据分布越陡峭,适合作为匹配区;好的适配区应该具有等方性系数大,变异系数大,峰态系数大,偏态系数小的性质,因此定义综合特征参数T,该参数与等方性系数、变异系数、峰态系数成正比,与偏态系数成反比,定义式如下所示:
T &Proportional; &Delta; &CenterDot; Cv &CenterDot; Ce Cs - - - ( 14 )
根据上述公式分析重力异常数据,结合载体的实际路线设定综合特征参数T的阈值,从而进行适配区的选取。
本发明的有益效果:
与已有的重力辅助惯性导航适配区选取方法相比较,本发明方法将载体所要航行的区域分为若干个子区域,在对每个区域上重力异常数据进行分类统计分析基础上,通过各个子区域的重力异常三维曲面等方性系数、重力异常变异系数、重力异常偏态系数、重力异常峰态系数构建综合特征参数,从而对适配区进行选择,所选的适配区不仅是重力异常变化明显的区域,而且适配区的等方性较好,适配区内各个方向上均可进行匹配,得到较好的匹配效果,载体经过适配区时,行进路线不用做出太大的改变就可以进行匹配校正。
具体实施方式
本实施例中,取载体起点和终点所在的方形区域,将此区域分为1°×1°的若干子区域,对各个子区域进行统计分析。其过程如下:
步骤一、将航行区域分为若干个1°×1°子区域,分析每个子区域的等方性;
大自然中的一切景观包括重力场都具有分形的特征,分形学中Werierstrass-Mandelbrot函数(W-M函数)有类似指数为α的分数布朗运动上的数学特征,我们在此用于重力场重力异常序列三维曲面轮廓的表征,分形维数为D的W-M函数形式如下:
g ( x , y ) = G D - 1 &Sigma; n = n 1 &infin; cos 2 &pi; &gamma; n x &gamma; ( 2 - D ) n - - - ( 1 )
g(x,y)为重力场重力异常数据;D为分形维数;G为特征尺度系数;γ为大于1的常数,对于近似服从正态分布的重力场三维曲面,γ取1.5。
相应于(1)式的连续功率谱为
P ( &omega; ) = G 2 ( D - 1 ) 2 ln &gamma; &CenterDot; 1 &omega; ( 5 - 2 D ) , &gamma; n 1 < &omega; < + &infin; - - - ( 2 )
其中分形维数D可以由随机游走法求得,步骤如下:
(1)首先将重力异常值看作随机游走的结果,设定一个间隔值R(如R=3),计算每个点跟它上下左右相邻R的重力异常值差,即:
P=g(x2,y2)-g(x1,y1)      (3) 
其中||(x2,y2)-(x1,y1)||=R;
(2)计算P的期望E(P);
(3)计算E(P)和R的双对数log-log函数,即log(E)=(3-D)log R+C,取不同R得到一组对应的E(P)和R,进行线性拟合,得到斜率k,由D=3-k求得重力异常三维曲面图的分形维数;
根据(2)式的重力异常序列方差可表示为
m 0 = &Integral; w l w h p ( &omega; ) d&omega; = G 2 ( D - 1 ) 2 ln &gamma; &CenterDot; 1 4 - 2 D &CenterDot; [ 1 &omega; l ( 4 - 2 D ) - 1 &omega; h ( 4 - 2 D ) ] - - - ( 4 )
计算实际重力异常的方差可以得到m0,与(16)式相比即可求得特征尺度系数G,其中△为取样间隔,L为取样长度。
重力异常曲面是一个三维均匀的随机表面,γ阶轮廓谱距mr和表面谱距mr-q,r定义分别如下:
m &gamma; = &Integral; - &infin; + &infin; &omega; &gamma; P ( &omega; ) d&omega; - - - ( 5 )
m r - q , r = &Integral; - &infin; + &infin; &Integral; - &infin; + &infin; &omega; 1 r - q &omega; 2 q P ( &omega; 1 , &omega; 2 ) d &omega; 1 d &omega; 2 - - - ( 6 )
根据Longuet-Higgins等式,θ方向上的γ阶轮廓谱距mr与表面谱距mr-q,r有如下关系:
m r ( &theta; ) = m r - q , r C r q cos r - q &theta; sin q &theta; - - - ( 7 )
由(7)式可以推出在θ方向上轮廓的二阶谱距m2(θ)为
m2(θ)=m20cos2θ+2m11cosθsinθ+m02sin2θ     (8) 
根据(8)式,分别在三个方向上测量重力异常数据,可以求得重力异常序列轮廓表面的二阶谱距,根据(9)式可求得重力异常三维曲面的表面谱距:
m 20 m 11 m 02 = T 2 - 1 m 2 ( &theta; 1 ) m 2 ( &theta; 2 ) m 2 ( &theta; 3 ) - - - ( 9 )
其中 T 2 = cos 2 &theta; 1 2 cos &theta; 1 sin &theta; 1 sin 2 &theta; 1 cos 2 &theta; 2 2 cos &theta; 2 sin &theta; 2 sin 2 &theta; 2 cos 2 &theta; 2 2 cos &theta; 3 sin &theta; 3 sin 2 &theta; 3 , 表面谱距m20和m02为相互垂直的两个方向上的斜率方差,而m11为两个方向上的协方差。
我们用下式评定各向异性表面的等方性:
&Delta; = 2 &Delta; 2 M 2 = 2 m 20 &CenterDot; m 02 - m 11 2 m 20 + m 02 - - - ( 10 )
当m11=0,m20=m02,△=1时,表示重力异常三维曲面是完全等方性的,所以△越接近1,等方性越好。
步骤二、分析每个区域的重力异常变异系数Cv;
在概率论与统计学中,变异系数又称为离散系数,是概率分布离散程度的一个归一化量度,其定义为标准差与平均值之比,也称为标准离差率或单位风险,是衡量数据变异程度的另一个统计量。
变异系数没有量纲,同时又按照其均数大小进行了标准化,这样可以对数据进行客观的比较。因此变异系数可以消除单位和平均数对重力场变异程度比较的影响,比标准差更能衡量重力场变异程度。重力异常变异系数的计算公式如下:
Cv = G . S . D G . mean &times; 100 % - - - ( 11 )
步骤三、分析每个子区域的重力异常偏态系数Cs和峰态系数Ce;
首先分析子区域的重力异常序列,得到归中后的重力异常值
&Delta;g i , j * = &Delta; g i , j - 1 NM &Sigma; i = 1 N &Sigma; j = 1 M &Delta;g i , j - - - ( 12 )
用归中后的重力异常值作为起始观测数据,得到均方差C0
C 0 = 1 NM &Sigma; i = 1 N &Sigma; j = 1 M &Delta;g i , j * 2 - - - ( 13 )
偏态系数又称偏差系数,重力异常偏态系数是说明重力异常序列在子区域内分布不对称程度的统计参数,偏态系数Cs绝对值越大,表明重力异常序列在子区域内的偏斜程度越大。
C s = N ( N - 1 ) ( N - 2 ) &CenterDot; M ( M - 1 ) ( M - 2 ) &CenterDot; &Sigma; i = 1 N &Sigma; j = 1 M &Delta;g i , j * 3 C 0 3 / 2 - - - ( 14 )
重力异常三维峰态系数表征重力异常分布在平均值处峰值高低的特征数, 简单来讲,峰态系数Ce描述的是重力异常序列分布形态的陡缓程度。
C e = N 2 - 2 N + 3 ( N - 1 ) ( N - 2 ) ( N - 3 ) &CenterDot; M 2 - 2 M + 3 ( M - 1 ) ( M - 2 ) ( M - 3 ) &CenterDot; &Sigma; i = 1 N &Sigma; j = 1 M &Delta;g i , j * 4 C 0 2 - 3 ( 2 N - 3 ) N ( N - 1 ) ( N - 2 ) ( N - 3 ) &CenterDot; 3 ( 2 M - 3 ) M ( M - 1 ) ( M - 2 ) ( M - 3 ) &CenterDot; [ &Sigma; i = 1 N &Sigma; j = 1 M &Delta;g i , j * 2 ] 2 c 0 2 - - - ( 15 )
步骤四、结合等方性系数△、变异系数Cv、偏态系数Cs和峰态系数Ce构建综合特征参数T,从而对适配区进行选取;
在对重力异常数据进行分析选择适配区时,好的适配区应该既能保证适配区内的重力异常有明显的变化,又能保证适配区可以在各个方向上都能得到较好的匹配效果,这样无论载体以什么角度进入适配区,不用对航向做出太大的变化就可以对惯导进行匹配校正。
好的适配区应该具有等方性系数大,变异系数大,峰态系数大,偏态系数小的性质,因此定义综合特征参数T,该参数与等方性系数、变异系数、峰态系数成正比,与偏态系数成反比,定义式如下所示:
T &Proportional; &Delta; &CenterDot; Cv &CenterDot; Ce Cs - - - ( 16 )
根据上述公式分析重力异常数据,结合载体的实际路线设定综合特征参数T的阈值,从而进行适配区的选取。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进,或者对其中部分技术特征进行等同替换,这些改进和替换也应视为本发明的保护范围。

Claims (4)

1.一种重力辅助惯性导航适配区的选择方法,其特征在于:包括如下步骤:
步骤一、将航行区域分为若干个子区域,分析每个子区域的等方性系数△;
步骤二、根据每个子区域的等方性,分析每个子区域的重力异常变异系数CV
步骤三、分析每个子区域的重力异常偏态系数Cs和峰态系数Ce
步骤四、结合等方性系数△、变异系数CV、偏态系数Cs和峰态系数Ce构建综合特征参数,从而对适配区进行选取。
2.如权利要求1所述的一种重力辅助惯性导航适配区的选择方法,其特征在于:进一步地,对于重力异常序列三维曲面轮廓的表征,分形维数为D的W-M函数形式如下:
g ( x , y ) = G D - 1 &Sigma; n = n 1 &infin; cos 2 &pi;&gamma; n x &gamma; ( 2 - D ) n - - - ( 1 )
g(x,y)为重力异常数据;D为分形维数;G为特征尺度系数;γ为大于1的常数,γn为随机过程的空间频率;n1与地形表面轮廓的最低截断频率相对应;
(1)式的连续功率谱函数为
P ( &omega; ) = G 2 ( D - 1 ) 2 ln &gamma; &CenterDot; 1 &omega; ( 5 - 2 D ) , &gamma; n 1 < &omega; < + &infin; - - - ( 2 )
其中ω为频率,重力异常曲面是一个三维均匀的随机表面,γ阶轮廓谱距mr和表面谱距mγ-q,q,q=0,1…γ,定义分别如下:
m &gamma; = &Integral; - &infin; + &infin; &omega; &gamma; P ( &omega; ) d&omega; - - - ( 3 )
m &gamma; - q , q = &Integral; - &infin; + &infin; &Integral; - &infin; + &infin; &omega; 1 &gamma; - q &omega; 2 &gamma; Q ( &omega; 1 , &omega; 2 ) d&omega; 1 d&omega; 2 - - - ( 4 )
其中Q(ω12)为重力异常轮廓的二元频谱函数,根据Longuet-Higgins等式,设定一个参考方向x,与参考方向成θ角方向上的γ阶轮廓谱距mr与表面谱距mr-q,q有如下关系:
m r ( &theta; ) = m r - q , q C r q cos r - q &theta; sin q &theta; - - - ( 5 )
由(5)式推出在θ方向上轮廓的二阶谱距m2(θ)为
m2(θ)=m20cos2θ+2m11cosθsinθ+m02sin2θ  (6)
根据(6)式可得,分别在三个方向上测量重力异常数据,求得重力异常序列轮廓表面的二阶谱距,通过(7)式即求得重力异常序列的表面谱距;
m 20 m 11 m 02 = T 2 - 1 m 2 ( &theta; 1 ) m 2 ( &theta; 2 ) m 2 ( &theta; 3 ) - - - ( 7 )
其中 T 2 = cos 2 &theta; 1 2 cos &theta; 1 sin &theta; 1 sin 2 &theta; 1 cos 2 &theta; 2 2 cos &theta; 2 sin &theta; 2 sin 2 &theta; 2 cos 2 &theta; 2 2 cos &theta; 3 sin &theta; 3 sin 2 &theta; 2 , θ1,θ2,θ3为任选的与参考方向成一定夹角的三个方向,表面谱距m20和m02为相互垂直的两个方向上的斜率方差,而m11为两个方向上的协方差;
采用下式评定各向异性表面的等方性系数△:
&Delta; = 2 &Delta; 2 M 2 = 2 m 20 &CenterDot; m 02 - m 11 2 m 20 + m 02 - - - ( 8 )
其中M2=m20+m02,当m11=0,m20=m02,△=1时,表示重力异常三维曲面是完全等方性的,所以△越接近1,等方性越好。
3.如权利要求2所述的一种重力辅助惯性导航适配区的选择方法,其特征在于:进一步地,对于近似服从正态分布的重力异常三维曲面,γ取1.5。
4.如权利要求1或2或3所述的一种重力辅助惯性导航适配区的选择方法,其特征在于:进一步地,选取的原则为:在以上四个参数中等方性是测量坐标系旋转时,该随机表面各条轮廓高度方向的概率分布不变的特征,等方性系数越大的区域越适合做匹配区;重力异常变异系数是衡量重力异常数据离散程度的参数,变异系数越大的区域越适合做匹配区;重力异常偏态系数是衡量重力异常数据不对称性质的参数,偏态系数越大,表明该区域数据的偏斜越严重,只有在某个方向上能获得较好的匹配效果,因此偏态系数大的区域不适合作为适配区;重力异常峰态系数是衡量重力异常聚中性的性质,重力异常峰态系数越大,则重力异常数据分布越陡峭,适合作为匹配区;好的适配区应该具有等方性系数大,变异系数大,峰态系数大,偏态系数小的性质,因此定义综合特征参数T,该参数与等方性系数、变异系数、峰态系数成正比,与偏态系数成反比,定义式如下所示:
T &Proportional; &Delta; &CenterDot; Cv &CenterDot; Ce Cs - - - ( 14 )
根据上述公式分析重力异常数据,结合载体的实际路线设定综合特征参数T的阈值,从而进行适配区的选取。
CN201510300526.6A 2015-06-03 2015-06-03 一种重力辅助惯性导航适配区的选择方法 Active CN105021182B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510300526.6A CN105021182B (zh) 2015-06-03 2015-06-03 一种重力辅助惯性导航适配区的选择方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510300526.6A CN105021182B (zh) 2015-06-03 2015-06-03 一种重力辅助惯性导航适配区的选择方法

Publications (2)

Publication Number Publication Date
CN105021182A true CN105021182A (zh) 2015-11-04
CN105021182B CN105021182B (zh) 2018-04-24

Family

ID=54411331

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510300526.6A Active CN105021182B (zh) 2015-06-03 2015-06-03 一种重力辅助惯性导航适配区的选择方法

Country Status (1)

Country Link
CN (1) CN105021182B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105953799A (zh) * 2016-04-22 2016-09-21 北京理工大学 基于熵值法的重力场适配区内水下运载体航迹规划方法
CN108871342A (zh) * 2018-07-06 2018-11-23 北京理工大学 基于纹理特征的水下重力辅助惯性导航适配区选取方法
CN109141426A (zh) * 2018-08-10 2019-01-04 中国空间技术研究院 一种水下重力匹配导航适配区的优选方法
CN110031000A (zh) * 2019-05-21 2019-07-19 北京理工大学 一种重力辅助惯性导航区域适配性的评价方法
CN110031001A (zh) * 2019-05-21 2019-07-19 北京理工大学 一种重力辅助惯性导航的适配区选取方法
CN110081884A (zh) * 2019-05-21 2019-08-02 北京理工大学 基于虚拟航向的重力辅助惯性导航区域适配性评价方法
CN110906930A (zh) * 2019-12-18 2020-03-24 中国人民解放军61540部队 一种联合auv的水下重力灯塔潜艇导航方法及系统
CN111044041A (zh) * 2019-12-31 2020-04-21 北京理工大学 基于重力场三维特征的重力辅助惯性导航适配区选取方法
CN111504302A (zh) * 2020-05-13 2020-08-07 中国人民解放军61540部队 一种联合海动力信息的重力灯塔导航路径规划方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1508020B1 (en) * 2002-05-30 2007-04-18 Rafael - Armament Development Authority Ltd. Airborne reconnaissance system
CN101358850A (zh) * 2008-09-19 2009-02-04 北京航空航天大学 一种基于地磁熵的匹配区选择方法
CN102128625A (zh) * 2010-12-08 2011-07-20 北京航空航天大学 重力辅助惯性导航系统中重力图匹配的初始匹配方法
CN102788578A (zh) * 2012-07-25 2012-11-21 中国人民解放军海军工程大学 一种基于局部重力场逼近的匹配导航方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1508020B1 (en) * 2002-05-30 2007-04-18 Rafael - Armament Development Authority Ltd. Airborne reconnaissance system
CN101358850A (zh) * 2008-09-19 2009-02-04 北京航空航天大学 一种基于地磁熵的匹配区选择方法
CN102128625A (zh) * 2010-12-08 2011-07-20 北京航空航天大学 重力辅助惯性导航系统中重力图匹配的初始匹配方法
CN102788578A (zh) * 2012-07-25 2012-11-21 中国人民解放军海军工程大学 一种基于局部重力场逼近的匹配导航方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
MAANI A O, AL-DAOUD E: "Improving an approximate string matching algorithm", 《INTERNATIONAL JOURNAL OF COMPUTER SCIENCE ISSUES》 *
张国忠等: "基于等方性系数的地形可导航性的方向性研究", 《国防科技大学学报》 *
李珊珊: "水下重力辅助惯性导航的理论与方法研究", 《中国博士学位论文全文数据库 基础科学辑》 *
王虎彪等: "重力垂直梯度导航值特征分布与可导航性分析", 《测绘学报》 *
蔡体菁,陈鑫巍: "基于层次分析法的重力匹配区域选择准则", 《中国惯性技术学报》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105953799A (zh) * 2016-04-22 2016-09-21 北京理工大学 基于熵值法的重力场适配区内水下运载体航迹规划方法
CN105953799B (zh) * 2016-04-22 2018-12-04 北京理工大学 基于熵值法的重力场适配区内水下运载体航迹规划方法
CN108871342A (zh) * 2018-07-06 2018-11-23 北京理工大学 基于纹理特征的水下重力辅助惯性导航适配区选取方法
CN109141426A (zh) * 2018-08-10 2019-01-04 中国空间技术研究院 一种水下重力匹配导航适配区的优选方法
CN109141426B (zh) * 2018-08-10 2020-11-10 中国空间技术研究院 一种水下重力匹配导航适配区的方法
CN110081884A (zh) * 2019-05-21 2019-08-02 北京理工大学 基于虚拟航向的重力辅助惯性导航区域适配性评价方法
CN110031001A (zh) * 2019-05-21 2019-07-19 北京理工大学 一种重力辅助惯性导航的适配区选取方法
CN110031000A (zh) * 2019-05-21 2019-07-19 北京理工大学 一种重力辅助惯性导航区域适配性的评价方法
CN110031001B (zh) * 2019-05-21 2020-12-11 北京理工大学 一种重力辅助惯性导航的适配区选取方法
CN110906930A (zh) * 2019-12-18 2020-03-24 中国人民解放军61540部队 一种联合auv的水下重力灯塔潜艇导航方法及系统
CN110906930B (zh) * 2019-12-18 2021-08-27 中国人民解放军61540部队 一种联合auv的水下重力灯塔潜艇导航方法及系统
CN111044041A (zh) * 2019-12-31 2020-04-21 北京理工大学 基于重力场三维特征的重力辅助惯性导航适配区选取方法
CN111044041B (zh) * 2019-12-31 2021-04-27 北京理工大学 基于重力场三维特征的重力辅助惯性导航适配区选取方法
CN111504302A (zh) * 2020-05-13 2020-08-07 中国人民解放军61540部队 一种联合海动力信息的重力灯塔导航路径规划方法及系统

Also Published As

Publication number Publication date
CN105021182B (zh) 2018-04-24

Similar Documents

Publication Publication Date Title
CN105021182A (zh) 一种重力辅助惯性导航适配区的选择方法
Setianto et al. Comparison of kriging and inverse distance weighted (IDW) interpolation methods in lineament extraction and analysis
CN102136155B (zh) 基于三维激光扫描数据的物体立面矢量化方法和系统
TWI464443B (zh) 人工智慧地震預警方法與地震偵測系統
CN108459318A (zh) 基于遥感技术的潜在滑坡早期识别方法
US20180158202A1 (en) Omnidirectional Roughness Algorithm for Topographic Signature Analysis of Lunar Craters
US9383325B2 (en) Method for determining geometric relationships of crystal reciprocal vectors on two-dimensional planes obtained from single EBSD pattern
CN104881547B (zh) 一种用于定向井井眼轨迹的误差分析方法
KR101541519B1 (ko) 레이더 관측자료를 3차원 격자 데이터로 구축하여 활용하는 강우량 추정 장치
EP3437396B1 (en) Determining a location of a wireless transmitter
CN105627938A (zh) 一种基于车载激光扫描点云的路面沥青厚度检测方法
Liu et al. A Kriging algorithm for location fingerprinting based on received signal strength
CN105550688A (zh) 点云数据的分类方法及装置
CN103456020A (zh) 基于treelet特征融合的遥感图像变化检测方法
JP2017138124A5 (zh)
Davies et al. Assessment of which scales of orography can be credibly resolved in a numerical model
CN106568412A (zh) 一种不规则形状地块的面积测量方法
CN102508198A (zh) 一种基于最大似然估计的目标被动定位方法
Vaz et al. Tectonic lineament mapping of the Thaumasia Plateau, Mars: Comparing results from photointerpretation and a semi-automatic approach
Smith et al. A gravity wave drag matrix for complex terrain
Sarp et al. Morphotectonic properties of Yenicaga basin area in Turkey
CN108460422B (zh) 基于深度分布特征的海底地貌类型识别方法
CN112986948B (zh) 基于InSAR技术的建筑形变监测方法和装置
CN106910218A (zh) 基于大区域栅格数字高程模型的地表表面面积计算方法
CN102221690B (zh) 高分辨率卫星遥感定位的组合岭估计光束法平差方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20190325

Address after: 100029 3-4 Floor of Crawford Building, 23 Huixin East Street, Chaoyang District, Beijing

Patentee after: Beijing techlink intelligent Polytron Technologies Inc

Address before: 100081 No. 5, Zhongguancun South Street, Haidian District, Beijing

Patentee before: BEIJING INSTITUTE OF TECHNOLOGY

TR01 Transfer of patent right