CN106483501B - 一种基于dop值分析声学定位系统多应答器最优布阵方法 - Google Patents
一种基于dop值分析声学定位系统多应答器最优布阵方法 Download PDFInfo
- Publication number
- CN106483501B CN106483501B CN201510551195.3A CN201510551195A CN106483501B CN 106483501 B CN106483501 B CN 106483501B CN 201510551195 A CN201510551195 A CN 201510551195A CN 106483501 B CN106483501 B CN 106483501B
- Authority
- CN
- China
- Prior art keywords
- transponders
- optimal
- transponder
- circular area
- positioning system
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于DOP值分析声学定位系统多应答器最优布阵方法。本发明属于声学定位系统多应答器布阵方法技术领域,具体涉及一种适用于多应答器声学定位系统的应答器空间几何最优布阵方法。本方法建立了一种从空间几何角度定量分析多应答器布局优劣的思路方法;该理论可直接用于应答器阵的布放和应答器位置的标校。并以三应答器为例,对空间几何最优布局思路进行详细分析,利用本方法定量的分析得到了平面圆形区域内三应答器阵的最优布局方案。
Description
技术领域
本发明属于声学定位系统多应答器布阵方法技术领域,具体涉及一种适用于多应答器声学定位系统的应答器空间几何最优布阵方法。
背景技术
对于声学定位系统而言,当有三个应答器甚至更多应答器能够使用时,应答器的空间几何阵型布置以及应答器和基阵间的相对几何位置关系都和定位精度有一定的关系。尤其对于以球面交汇为定位原理的声学定位系统,多应答器空间几何最优布阵方案设计对提高整体精度有重要意义。
随着声学定位技术在水下导航领域的发展和应用,国内外研制了包括长基线、短基线和超短基线在内的多种声学定位系统典型产品,声学定位系统研制单位专注于开发出高性能(传感器误差小)的产品,但仍缺少一种能够定量分析多应答器空间几何布局优劣的方法,未能在传感器误差一定的情况下,通过多应答器的阵型优化来减小声学定位系统系统误差的方法。
因此,亟需研制一种适用于多应答器声学定位系统的应答器空间几何最优布阵方法,从而能够定量分析多应答器空间几何布局优劣,
发明内容
本发明要解决的技术问题是提供一种声学定位系统多应答器空间几何最优布阵方法,从空间几何角度定量的分析多应答器布局的优劣,进而得到一种多应答器阵型布放的最优化设计,一定程度上减小声学定位系统系统误差。
为了实现这一目的,本发明采取的技术方案是:
一种基于DOP值分析声学定位系统多应答器最优布阵方法,包括以下步骤:
(1)利用超短基线定位系统USBL输出的方位角α和高度角β,建立DOP观测矩阵A:
其中,αi和βi分别表示第i个应答器的方位角和高度角,i=1~n,n=3;
不考虑USBL测量的时间差,建立3应答器的位置误差最小二乘解:
其中Δρ为观测噪声向量,其方差为σ2;
简化后的位置误差方差为:
(2)建立3应答器DOP观测方程的权系数矩阵Q:
(3)根据式0定义多应答器USBL定位的三种DOP值的计算公式如下:
水平几何精度因子(Horizontal Dilution of Precision,HDOP):
垂直几何精度因子(Vertical Dilution of Precision,VDOP):
三维几何精度因子(Positional Dilution of Precision,PDOP):
(4)确定平面圆形区域内三应答器的相对位置关系
(4.1)确定平面圆形区域内观测矩阵中高度角和方位角
半径为R的平面圆形区域内,观测矩阵中高度角β1,β2,β3均设定为0,设三个应答器相对于基阵的方位角依次为α1,α2,α3;
DOP值观测矩阵为:
权系数矩阵为:
其中:
(4.2)分析平面圆形区域内水平几何精度因子HDOP值,将0中非零部分取出,构成平面内简化HDOP权系数矩阵QH:
其中:
则HDOP为
最优布局需要设定A=Amax,α1+α2+α3=360°,得到当三个方位角相差均为120°时,可获得半径为R的平面圆形区域内的最优布局,从而确定三个应答器分别在边长为的等边三角形三个顶点上;
结合0和0确定,位于同一水平面上的三个应答器的方位角依次为:α,α+120°,α+240°,在保证三个应答器两两之间方位角差值为120°的情况下,α为任一角度,HDOP值不会受到影响;
(5)确定圆形区域内半径和深度关系的最优解
在步骤(4)平面圆形区域内最优布阵方案的前提下,求解USBL绝对定位三维几何精度因子最优点;
取半径为R的圆的内切等边三角形的三个顶点为应答器的布放位置,三个方位角依次为α1=0°,α2=120°,α3=240°,从数学几何关系计算确定,最优解在圆的轴向上,其到各平面圆形区域内个应答器的高度h都相等;令高度角β1=β2=β3=β;
其中,
DOP值观测矩阵为:
权系数矩阵为:
平面定位精度因子为:
垂直几何精度因子为:
三维定位精度因子为:
最优布局需要设定PDOP值最小,化简上式:
可得当且仅当即时,PDOP值取得最小值,此时为确定圆形区域内半径和深度关系的最优解。
进一步的,如上所述的一种基于DOP值分析声学定位系统多应答器最优布阵方法,应答器数量多于3个时,取n>3,采用步骤(1)~(5)的方法求解最优布阵。
本发明的有益效果在于,建立了一种从空间几何角度定量分析多应答器布局优劣的思路方法;该理论可直接用于应答器阵的布放和应答器位置的标校。并以三应答器为例,对空间几何最优布局思路进行详细分析,利用本方法定量的分析得到了平面圆形区域内三应答器阵的最优布局方案。此外,该方法还有助于提高应答器的绝对位置标定精度。
具体实施方式
下面结合具体实施例对本发明一种基于DOP值分析声学定位系统多应答器最优布阵方法进行详细说明。该方法包括以下步骤:
(1)利用超短基线定位系统USBL输出的方位角α和高度角β,建立DOP观测矩阵A:
其中,αi和βi分别表示第i个应答器的方位角和高度角,i=1~n,n=3;
不考虑USBL测量的时间差,建立3应答器的位置误差最小二乘解:
其中Δρ为观测噪声向量,其方差为σ2;
简化后的位置误差方差为:
(2)建立3应答器DOP观测方程的权系数矩阵Q:
(3)根据式0定义多应答器USBL定位的三种DOP值的计算公式如下:
水平几何精度因子(Horizontal Dilution of Precision,HDOP):
垂直几何精度因子(Vertical Dilution of Precision,VDOP):
三维几何精度因子(Positional Dilution of Precision,PDOP):
(4)确定平面圆形区域内三应答器的相对位置关系
(4.1)确定平面圆形区域内观测矩阵中高度角和方位角
半径为R的平面圆形区域内,观测矩阵中高度角β1,β2,β3均设定为0,设三个应答器相对于基阵的方位角依次为α1,α2,α3;
DOP值观测矩阵为:
权系数矩阵为:
其中:
(4.2)分析平面圆形区域内水平几何精度因子HDOP值,将0中非零部分取出,构成平面内简化HDOP权系数矩阵QH:
其中:
则HDOP为
最优布局需要设定A=Amax,α1+α2+α3=360°,得到当三个方位角相差均为120°时,可获得半径为R的平面圆形区域内的最优布局,从而确定三个应答器分别在边长为的等边三角形三个顶点上;
结合0和0确定,位于同一水平面上的三个应答器的方位角依次为:α,α+120°,α+240°,在保证三个应答器两两之间方位角差值为120°的情况下,α为任一角度,HDOP值不会受到影响;
(5)确定圆形区域内半径和深度关系的最优解,在步骤(4)平面圆形区域内最优布阵方案的前提下,求解USBL绝对定位三维几何精度因子最优点;
即分析基阵位置精度最优解所在的位置,有助于提高各应答器的绝对位置标定精度。在平面圆形区域内最优布阵方案的前提下,求解USBL绝对定位三维几何精度因子最优点。这个最优点的分析可以对USBL应答器的布放使用以及标定提供理论基础,比如应答器初始位置的高精度校准,或者在固定半径为R的圆形区域内应答器应该布放多深(h)才能获得相对较小的PDOP值。
取半径为R的圆的内切等边三角形的三个顶点为应答器的布放位置,三个方位角依次为α1=0°,α2=120°,α3=240°,从数学几何关系计算确定,最优解在圆的轴向上,其到各平面圆形区域内个应答器的高度h都相等;令高度角β1=β2=β3=β;
其中,
DOP值观测矩阵为:
权系数矩阵为:
平面定位精度因子为:
垂直几何精度因子为:
三维定位精度因子为:
最优布局需要设定PDOP值最小,化简上式:
可得当且仅当即时,PDOP值取得最小值,此时为确定圆形区域内半径和深度关系的最优解。
当应答器数量多于3个时,取n>3,采用步骤(1)~(5)的方法求解最优布阵从而扩展本方法的应用范围,该理论可直接用于应答器阵的布放和应答器位置的标校。
Claims (2)
1.一种基于DOP值分析声学定位系统多应答器最优布阵方法,其特征在于,包括以下步骤:
(1)利用超短基线定位系统USBL输出的方位角α和高度角β,建立DOP观测矩阵A:
其中,αi和βi分别表示第i个应答器的方位角和高度角,i=1~n,n=3;不考虑USBL测量的时间差,建立3应答器的位置误差最小二乘解:
其中Δρ为观测噪声向量,其方差为σ2;
简化后的位置误差方差为:
(2)建立3应答器DOP观测方程的权系数矩阵Q:
(3)根据式(4)定义多应答器USBL定位的三种DOP值的计算公式如下:
水平几何精度因子HDOP:
垂直几何精度因子VDOP:
三维几何精度因子PDOP:
(4)确定平面圆形区域内三应答器的相对位置关系
(4.1)确定平面圆形区域内观测矩阵中高度角和方位角
半径为R的平面圆形区域内,观测矩阵中高度角β1,β2,β3均设定为0,设三个应答器相对于基阵的方位角依次为α1,α2,α3;
DOP值观测矩阵为:
权系数矩阵为:
其中:
(4.2)分析平面圆形区域内水平几何精度因子HDOP值,将(9)中非零部分取出,构成平面内简化HDOP权系数矩阵QH:
其中:
则HDOP为
最优布局需要设定A=Amax,α1+α2+α3=360°,得到当三个方位角相差均为120°时,可获得半径为R的平面圆形区域内的最优布局,从而确定三个应答器分别在边长为的等边三角形三个顶点上;
结合(11)和(12)确定,位于同一水平面上的三个应答器的方位角依次为:α,α+120°,α+240°,在保证三个应答器两两之间方位角差值为120°的情况下,α为任一角度,HDOP值不会受到影响;
(5)确定圆形区域内半径和深度关系的最优解
在步骤(4)平面圆形区域内最优布阵方案的前提下,求解USBL绝对定位三维几何精度因子最优点;
取半径为R的圆的内切等边三角形的三个顶点为应答器的布放位置,三个方位角依次为α1=0°,α2=120°,α3=240°,从数学几何关系计算确定,最优解在圆的轴向上,其到各平面圆形区域内个应答器的高度h都相等;令高度角β1=β2=β3=β;
其中,
DOP值观测矩阵为:
权系数矩阵为:
水平几何精度因子为:
垂直几何精度因子为:
三维几何精度因子为:
最优布局需要设定PDOP值最小,化简上式:
当时,PDOP值取得最小值,此时为确定圆形区域内半径和深度关系的最优解。
2.如权利要求1所述的一种基于DOP值分析声学定位系统多应答器最优布阵方法,其特征在于,应答器数量多于3个时,取n>3,采用步骤(1)~(5)的方法求解最优布阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510551195.3A CN106483501B (zh) | 2015-09-01 | 2015-09-01 | 一种基于dop值分析声学定位系统多应答器最优布阵方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510551195.3A CN106483501B (zh) | 2015-09-01 | 2015-09-01 | 一种基于dop值分析声学定位系统多应答器最优布阵方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106483501A CN106483501A (zh) | 2017-03-08 |
CN106483501B true CN106483501B (zh) | 2019-04-23 |
Family
ID=58235652
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510551195.3A Active CN106483501B (zh) | 2015-09-01 | 2015-09-01 | 一种基于dop值分析声学定位系统多应答器最优布阵方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106483501B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108120955A (zh) * | 2017-12-04 | 2018-06-05 | 中国科学院深海科学与工程研究所 | 一种用于长基线水声定位系统性能评估的方法 |
CN110703196B (zh) * | 2019-10-16 | 2021-07-13 | 沈阳建筑大学 | 基于超声波传感器三位置测量的环境轮廓重构方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1908682A1 (en) * | 2000-03-20 | 2008-04-09 | Francois Bernard | Arrangement of an accoustic array with a sound velocity meter |
CN101799533A (zh) * | 2010-04-09 | 2010-08-11 | 华北电力大学 | 平面四元阵列电站锅炉承压管泄漏定位方法 |
CN103454616A (zh) * | 2013-08-27 | 2013-12-18 | 西北工业大学 | 一种十字型振速梯度水听器的方位估计方法 |
CN103728642A (zh) * | 2013-12-30 | 2014-04-16 | 深圳先进技术研究院 | 定位方法和系统 |
GB2508314A (en) * | 2012-07-17 | 2014-05-28 | Silixa Ltd | Detecting vibration and monitoring position of a structure using an optical fibre distributed acoustic sensor. |
-
2015
- 2015-09-01 CN CN201510551195.3A patent/CN106483501B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1908682A1 (en) * | 2000-03-20 | 2008-04-09 | Francois Bernard | Arrangement of an accoustic array with a sound velocity meter |
CN101799533A (zh) * | 2010-04-09 | 2010-08-11 | 华北电力大学 | 平面四元阵列电站锅炉承压管泄漏定位方法 |
GB2508314A (en) * | 2012-07-17 | 2014-05-28 | Silixa Ltd | Detecting vibration and monitoring position of a structure using an optical fibre distributed acoustic sensor. |
CN103454616A (zh) * | 2013-08-27 | 2013-12-18 | 西北工业大学 | 一种十字型振速梯度水听器的方位估计方法 |
CN103728642A (zh) * | 2013-12-30 | 2014-04-16 | 深圳先进技术研究院 | 定位方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106483501A (zh) | 2017-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110440692B (zh) | 激光跟踪仪与结构光3d扫描仪组合式测量标定方法 | |
CN108490473B (zh) | 一种融合gnss和uwb的无人机增强定位方法与系统 | |
CN108317993B (zh) | 一种集成gnss和激光跟踪仪的垂线偏差测量装置及方法 | |
CN109883443B (zh) | 一种线结构光传感器空间姿态标定方法 | |
US10670689B2 (en) | System and method for determining geo location of a target using a cone coordinate system | |
WO2018049945A1 (zh) | 一种高精度实时卫星定位方法 | |
CN106407560B (zh) | 表征大气各向异性的对流层映射函数模型的构建方法 | |
CN107894241A (zh) | 一种基于椭球拟合的无人机磁传感器校准方法、无人机 | |
CN110017812B (zh) | 一种垂线偏差、重力加速度和重力梯度的测量装置及方法 | |
CN103791868B (zh) | 一种基于虚拟球的空间标定体及其标定方法 | |
CN103323855A (zh) | 一种基线动态测量系统的精度获取方法 | |
CN108413988A (zh) | 机器人末端经纬仪坐标系快速标定方法 | |
CN103926058A (zh) | 非球面检测中采用自准平面镜测量光轴的方法 | |
CN110297220B (zh) | 一种双天线InSAR系统基线矢量的测量方法 | |
CN103364803A (zh) | 选星方法及应用该选星方法的卫星导航定位方法 | |
CN109959898B (zh) | 一种座底式水声被动定位基阵自校准方法 | |
US9383429B2 (en) | System and method for determining geo location of a target using locus of emitter positions (LEP) | |
CN108562872B (zh) | 一种超短基线水下声学定位校准时异常值探测的方法 | |
CN106483501B (zh) | 一种基于dop值分析声学定位系统多应答器最优布阵方法 | |
CN108413983A (zh) | 一种sins/usbl一体化定位系统安装误差标定的机械方法 | |
CN106568365A (zh) | 一种球面孔系复合位置度误差的检测与评定方法 | |
CN108917698B (zh) | 一种方位角计算方法 | |
CN111046584A (zh) | 星载仪器设备的精密调整方法 | |
CN106093987B (zh) | 一种应用于无人机的低成本差分gps系统及其实现方法 | |
CN102564383A (zh) | 基于轴面角确定三轴磁传感器轴间角的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |