CN108120955A - 一种用于长基线水声定位系统性能评估的方法 - Google Patents
一种用于长基线水声定位系统性能评估的方法 Download PDFInfo
- Publication number
- CN108120955A CN108120955A CN201711260848.8A CN201711260848A CN108120955A CN 108120955 A CN108120955 A CN 108120955A CN 201711260848 A CN201711260848 A CN 201711260848A CN 108120955 A CN108120955 A CN 108120955A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtd
- mfrac
- mtr
- 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.)
- Pending
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)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明提供了一种用于长基线水声定位系统性能评估的方法,选取GDOP、HDOP和VDOP中的一项或者几项作为评估参数,对长基线水声定位系统性能进行评估,评估方法如下:确定长基线声信标几何布局方案;对预投放三维空间区域划分网格点;计算声信标与指定网格点的空间三维距离;计算声信标与指定网格点的观测矩阵H;计算指定网格点的与声信标的权系数矩阵Q;计算指定网格点的评估参数值;计算其他网格点的评估参数值;对长基线声信标几何布局进行评估;根据评估结果,选择是否需要更改声信标布局。本发明可以从系统层面上评估长基线水声定位系统声信标的空间几何分布是否满足用户需求,从而可以指导长基线水声定位系统声信标布放的空间位置。
Description
技术领域
本发明涉及水下定位导航领域,具体涉及一种用于长基线水声定位系统性能评估的方法。
背景技术
地表目标的定位通常采用卫星定位系统作为基本的技术手段,惯性定位导航技术为辅助的方式。但当目标处于水下时,由于水对无线电波的强烈吸收,导致目标无法利用卫星导航系统等地表成熟的导航方法;而传统常用的惯性定位导航系统的误差随时间增加而大量累积。因此,深海航行器在水下长时自主运行时,需要使用声波作为信息载体来进行水声定位技术,从而作为惯性定位导航技术的有效辅助校准手段。水声定位技术现已经成功应用于军事、商用、民用工程。它能够为水下潜水器,如AUV(Autonomous UnderwaterVehicle)等重要的定位、导航工作进行支撑。通过在作业海区的水面工作母船或者其他水下平台加装和布放声学定位设备,实现水面对水下目标位置的实时监控从而完成海洋科学考察、海洋资源勘探及开发、深海空间站建设等大工程。
根据定位系统基线长度的差别,一般将其划分为长基线水声定位系统、短基线水声定位系统和超短基水声定位系统。本发明针对长基线水声定位系统,因此重点介绍长基线水声定位系统。长基线水声定位系统的基线一般在几百米到数千米的量级上,具有作用距离远、定位精度高等优点,主要应用于水下载体定位、海底电缆铺设、海上石油勘探等方面。通过测量载体到相对/绝对位置已知的海底应答器的相对距离,采用距离交汇的方法解算目标坐标。
现有针对长基线水声定位系统的研究主要集中在抗多径、抗距离模糊、声线修正、信号设计优化等如何开发出传感器误差更小的产品上,鲜有研究关注如何定量分析如何评估长基线水声定位系统的性能。因此,本发明将长基线水声定位系统的性能与GDOP、HDOP和VDOP等指标联系起来,对长基线水声定位系统的性能评估,从而优化长基线水声定位系统的声信标布局方案。
发明内容
本发明的目的在于提供一种用于长基线水声定位系统性能评估的方法,以解决上述问题。
本发明是通过以下技术方案实现的:一种用于长基线水声定位系统性能评估的方法,其特征在于,选取GDOP、HDOP和VDOP中的一项或者几项作为评估参数,对长基线水声定位系统性能进行评估,确定长基线水声定位系统声信标布局是否合理,其评估方法步骤如下:
S1.根据声信标数目N,确定长基线声信标几何布局方案;
S2.根据空间采样分辨率,对预投放的三维空间区域进行网格点划分;
S3.计算所述声信标与指定网格点的空间三维距离;
S4.计算所述声信标与指定网格点之间的观测矩阵H;
S5.计算所述指定网格点的与声信标之间的权系数矩阵Q;
S6.计算所述指定网格点的评估参数值;
S7.重复步骤S3~S6,计算其他网格点的评估参数值;
S8.选择合适的评估参数,对步骤S1中所确定的长基线声信标几何布局进行评估;
S9.根据评估结果,选择是否需要更改声信标布局。
进一步的,步骤S1中,所述声信标数目N必须大于或等于3,所述长基线声信标几何布局方案为同一平面的正N边形布局。
进一步的,步骤S2中,通过下列方法计算网格点数目M:
当所述预投放的三维空间区域为规则几何图形时,采用下列公式计算所划分的网格点数目M:
式中,a为预投放区域的长度值,b为预投放区域的宽度值,c为预投放区域的深度值,D1为预投放区域在长度上的空间采样分辨率,D2为预投放区域在宽度上的空间采样分辨率,D3为预投放区域在深度上的空间采样分辨率。
当预投放的三维空间区域为不规则几何图形时,同样采用式(1),按能将所述不规则图形包含在内的最小长方体计算所划分的网格点数目M,如果某一网格点正好处于声信标坐标,此时需要使用其他N-1个声信标坐标对此网格点进行评估。
进一步的,步骤S3中,根据下列公式计算所述声信标与指定网格点的空间三维距离:
式中,xi、yi、zi为第i个声信标的三维坐标(xi,yi,zi),di为第i个声信标与指定网格点的空间三维距离,其中,i取1、2、3、…、n,x、y、z为指定网格点的三维坐标(x,y,z)。
进一步的,步骤S4中,采用下列方法计算所述声信标与指定网格点之间的观测矩阵H:
所述声信标空间位置与指定网格点空间位置不处于同一平面时,采用下式计算观测矩阵H:
式中,xi、yi、zi为第i个声信标的三维坐标(xi,yi,zi),di为第i个声信标与指定网格点的空间三维距离,其中,i取1、2、3、…、n,x、y、z为指定网格点的三维坐标(x,y,z)。
进一步的,所述声信标空间位置与指定网格点空间位置处于同一平面,即zi等于z时,采用下式计算观测矩阵H:
xi、yi、zi为第i个声信标的三维坐标(xi,yi,zi),di为第i个声信标与指定网格点的空间三维距离,其中,i取1、2、3、…、n,x、y、z为指定网格点的三维坐标(x,y,z)。
进一步的,步骤S5中,根据下列公式计算所述指定网格点的与声信标之间的权系数矩阵Q:
进一步的,步骤S6中,采用下列公式计算所述指定网格点的评估参数GDOP、HDOP和VDOP值:
进一步的,所述声信标三维位置与指定网格点三维位置不处于同一平面时,选择评估参数GDOP、HDOP及VDOP对长基线水声定位系统声信标布局进行评估。
进一步的,所述声信标三维位置与指定网格点三维位置处于同一平面时,只选择评估参数HDOP对长基线水声定位系统声信标布局进行评估。
与现有技术相比,本发明达到的有益效果如下:
本发明提供的一种用于长基线水声定位系统性能评估的方法,在预设计及投放海底声信标前采用GDOP(Geometric Dilution Of Precision,几何精度因子)、HDOP(Horizontal Dilution Of Precision,水平精度因子)和VDOP(Vertical Dilution OfPrecision,垂直精度因子)作为评估参数,对长基线水声定位系统性能进行评估,从而确定长基线水声定位系统声信标布局是否合理,该评估方法填补了我国对长基线水声定位领域的一项研究空白,其实用性强,具有十分重要的科研与实用价值。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的优选实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例1提供的一种用于长基线水声定位系统性能评估的方法流程图。
图2为本发明实施例1提供的,c取350m时,指定区域的GDOP值及其等高线。
图3为为本发明实施例1提供的,c取400m时,指定区域的GDOP值及其等高线。
图4为本发明实施例1提供的,c取450m时,指定区域的GDOP值及其等高线。
图5为本发明实施例1提供的,c取350m时,指定区域的HDOP值及其等高线。
图6为本发明实施例1提供的,c取400m时,指定区域的HDOP值及其等高线。
图7为本发明实施例1提供的,c取450m时,指定区域的HDOP值及其等高线。
图8为本发明实施例1提供的,c取350m时,指定区域的VDOP值及其等高线。
图9为本发明实施例1提供的,c取400m时,指定区域的VDOP值及其等高线。
图10为本发明实施例1提供的,c取450m时,指定区域的VDOP值及其等高线。
图11为本发明实施例2提供的,c取350m时,指定区域的GDOP值及其等高线。
图12为本发明实施例2提供的,c取400m时,指定区域的GDOP值及其等高线。
图13为本发明实施例2提供的,c取450m时,指定区域的GDOP值及其等高线。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅只是本发明的一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的其他所有实施例,都属于本发明保护的范围。
实施例1
如图1所示,一种用于长基线水声定位系统性能评估的方法,其特征在于,选取GDOP(Geometric Dilution Of Precision,几何精度因子)、HDOP(Horizontal DilutionOf Precision,水平精度因子)和VDOP(Vertical Dilution Of Precision,垂直精度因子)中的一项或者几项作为评估参数,对长基线水声定位系统性能进行评估,确定长基线水声定位系统声信标布局是否合理,其评估方法步骤如下:
S1.根据声信标数目N,确定长基线声信标几何布局方案;
具体的,在本实施例中,采用四声信标对预投放区域进行几何布局,所采用的布局方案为正四边形布局,四声信标坐标位置分别为p1=[400,400,0]m,p2=[-400,400,0]m,p3=[-400,-400,0]m,p4=[400,-400,0]m。
S2.根据空间采样分辨率,对预投放的三维空间区域进行网格点划分;
具体的,设定所述预投放的三维空间区域范围为[a,b,c],其中a∈[-400,400]m,b∈[-400,400]m,c∈[350,450]m,设定所述空间采样分辨率为1m×1m×50m,当所述预投放的三维空间区域为规则几何图形时,采用下列公式计算所划分的网格点数目:
式中,a为预投放区域的长度值,b为预投放区域的宽度值,c为预投放区域的深度值,D1为预投放区域在长度上的空间采样分辨率,D2为预投放区域在宽度上的空间采样分辨率,D3为预投放区域在深度上的空间采样分辨率。优选的,的值均为不大于的整数。
进一步的,当预投放的三维空间区域为不规则几何图形时,同样采用式(1),按能将所述不规则图形包含在内的最小长方体计算所划分的网格点数目M,如果某一网格点正好处于声信标坐标,此时需要使用其他N-1个声信标坐标对此网格点进行评估。
S3.计算所述声信标与指定网格点的空间三维距离;
具体的,根据下列公式计算所述声信标与指定网格点的空间三维距离:
式中,xi、yi、zi为第i个声信标的三维坐标(xi,yi,zi),di为第i个声信标与指定网格点的空间三维距离,其中,i取1、2、3、…、n,x、y、z为指定网格点的三维坐标(x,y,z)。
S4.计算所述声信标与指定网格点之间的观测矩阵H;
具体的,采用下列方法计算所述声信标与指定网格点之间的观测矩阵H:
当所述声信标空间位置与指定网格点空间位置不处于同一平面时,采用下式计算观测矩阵H1:
式中,xi、yi、zi为第i个声信标的三维坐标(xi,yi,zi),di为第i个声信标与指定网格点的空间三维距离,其中,i取1、2、3、…、n,x、y、z为指定网格点的三维坐标(x,y,z)。
进一步的,所述声信标空间位置与指定网格点空间位置处于同一平面,即zi等于z时,采用下式计算观测矩阵H2:
式中,xi、yi、zi为第i个声信标的三维坐标(xi,yi,zi),di为第i个声信标与指定网格点的空间三维距离,其中,i取1、2、3、…、n,x、y、z为指定网格点的三维坐标(x,y,z)。
S5.计算所述指定网格点的与声信标之间的权系数矩阵Q;
具体的,根据下列公式计算所述指定网格点的与声信标之间的权系数矩阵Q:
式中,H可取H1或H2。
S6.计算所述指定网格点的评估参数值;
具体的,采用下列公式计算所述指定网格点的评估参数GDOP、HDOP和VDOP值:
进一步的,当所述声信标空间位置与指定网格点空间位置处于不同一平面时,需根据式(6)、式(7)、式(8)分别计算评估参数GDOP、HDOP和VDOP值;当所述声信标空间位置与指定网格点空间位置处于同一平面时,只需根据式(7)计算评估参数HDOP的值。优选的,当用户对三维定位精度比较关注时,可以重点选取GDOP作为评估参数;优选的,当用户对水平面定位的精度比较关注时,可以重点选取HDOP作为评估参数;优选的,当用户对垂直定位精度比较关注,可以重点选择VDOP作为评估参数。
S7.重复步骤S3~S6,计算其他网格点的评估参数GDOP、HDOP和VDOP值;
S8.选择合适的评估参数,对步骤S1中所确定的长基线声信标几何布局进行评估;
具体的,设定所述预投放的三维空间区域范围为[a,b,c],其中a∈[-400,400]m,b∈[-400,400]m,c∈[350,450]m,设定所述空间采样分辨率为1m×1m×50m,根据所述空间采样分辨率。c分别取350m、400m、450m,可获得GDOP与其等高线值图,如图2、图3、图4所示;可获得HDOP与其等高线值图,如图5、图6、图7所示;可获得VDOP与其等高线值图,如图8、图9、图10所示。
进一步的,对比图2、图3、图4,所述GDOP值随着感兴趣目标的离底高度增加而增加,且性能比较佳的位置位于中心位置。性能较佳区域近乎矩形,靠近声信标的位置,性能反而不佳。
进一步的,对比图5、图6、图7,所述HDOP值随着感兴趣目标的离底高度增加而增加,且性能比较佳的区域位于中心位置,形状近乎圆形。
进一步的,对比图8、图9、图10,所述VDOP值随着感兴趣目标的离底高度增加而减少,且性能比较佳的区域位于中心位置。
S9.根据评估结果,选择是否需要更改声信标布局。
实施例2
本实施例与实施例1的区别在于,所述四声信标坐标位置分别为p1=[400,400,100]m,p2=[-400,400,200]m,p3=[-400,-400,300]m,p4=[400,-400,400]m。
进一步的,设定所述预投放的三维空间区域范围为[a,b,c],其中a∈[-400,400]m,b∈[-400,400]m,c∈[350,450]m,设定所述空间采样分辨率为1m×1m×50m,c分别取350m、400m、450m,可获得GDOP与其等高线值图,如图11、图12、图13所示;
进一步的,对比图11、图12、图13,可以发现从GDOP值的角度来说,同一平面的正N边形声信标布局方案最优。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。
Claims (9)
1.一种用于长基线水声定位系统性能评估的方法,其特征在于,选取GDOP、HDOP和VDOP中的一项或者几项作为评估参数,对长基线水声定位系统性能进行评估,确定长基线水声定位系统声信标布局是否合理,其评估方法步骤如下:
S1.根据声信标数目N,确定长基线声信标几何布局方案;
S2.根据空间采样分辨率,对预投放的三维空间区域进行网格点划分;
S3.计算所述声信标与指定网格点的空间三维距离;
S4.计算所述声信标与指定网格点之间的观测矩阵H;
S5.计算所述指定网格点的与声信标之间的权系数矩阵Q;
S6.计算所述指定网格点的评估参数值;
S7.重复步骤S3~S6,计算其他网格点的评估参数值;
S8.选择合适的评估参数,对步骤S1中所确定的长基线声信标几何布局进行评估;
S9.根据评估结果,选择是否需要更改声信标布局。
2.根据权利要求1所述的一种用于长基线水声定位系统性能评估的方法,其特征在于,步骤S1中,所述声信标数目N必须大于或等于3,所述长基线声信标几何布局方案为同一平面的正N边形布局。
3.根据权利要求1所述的一种用于长基线水声定位系统性能评估的方法,其特征在于,步骤S2中,通过下列方法计算所划分的网格点数目M;
当所述预投放的三维空间区域为规则几何图形时,采用下列公式计算所划分的网格点数目M:
<mrow>
<mi>M</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mo>&lsqb;</mo>
<mfrac>
<mi>a</mi>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mo>&lsqb;</mo>
<mfrac>
<mi>b</mi>
<msub>
<mi>D</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mo>&lsqb;</mo>
<mfrac>
<mi>c</mi>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,a为预投放区域的长度值,b为预投放区域的宽度值,c为预投放区域的深度值,D1为预投放区域在长度上的空间采样分辨率,D2为预投放区域在宽度上的空间采样分辨率,D3为预投放区域在深度上的空间采样分辨率。
当预投放的三维空间区域为不规则几何图形时,同样采用式(1),按能将所述不规则图形包含在内的最小长方体计算所划分的网格点数目M,如果某一网格点正好处于声信标坐标,此时需要使用其他N-1个声信标坐标对此网格点进行评估。
4.根据权利要求1所述的一种用于长基线水声定位系统性能评估的方法,其特征在于,步骤S3中,根据下列公式计算所述声信标与指定网格点的空间三维距离:
<mrow>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,xi、yi、zi为第i个声信标的三维坐标(xi,yi,zi),di为第i个声信标与指定网格点的空间三维距离,其中,i取1、2、3、…、n,x、y、z为指定网格点的三维坐标(x,y,z)。
5.根据权利要求1所述的一种用于长基线水声定位系统性能评估的方法,其特征在于,步骤S4中,采用下列方法计算所述声信标与指定网格点之间的观测矩阵H:
所述声信标空间位置与指定网格点空间位置不处于同一平面时,采用下式计算观测矩阵H:
<mrow>
<mi>H</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mi>y</mi>
</mrow>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mi>z</mi>
</mrow>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mi>y</mi>
</mrow>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mi>z</mi>
</mrow>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>y</mi>
</mrow>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>z</mi>
</mrow>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,xi、yi、zi为第i个声信标的三维坐标(xi,yi,zi),di为第i个声信标与指定网格点的空间三维距离,其中,i取1、2、3、…、n,x、y、z为指定网格点的三维坐标(x,y,z)。
6.根据权利要求5所述的一种用于长基线水声定位系统性能评估的方法,其特征在于,所述声信标空间位置与指定网格点空间位置处于同一平面,即zi等于z时,采用下式计算观测矩阵H:
<mrow>
<mi>H</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mi>y</mi>
</mrow>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mi>y</mi>
</mrow>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>y</mi>
</mrow>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
xi、yi、zi为第i个声信标的三维坐标(xi,yi,zi),di为第i个声信标与指定网格点的空间三维距离,其中,i取1、2、3、…、n,x、y、z为指定网格点的三维坐标(x,y,z)。
7.根据权利要求5或6所述的一种用于长基线水声定位系统性能评估的方法,其特征在于,步骤S5中,根据下列公式计算所述指定网格点的与声信标之间的权系数矩阵Q:
<mrow>
<mi>Q</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>H</mi>
<mi>T</mi>
</msup>
<mi>H</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>Q</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Q</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Q</mi>
<mn>13</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Q</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Q</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Q</mi>
<mn>23</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Q</mi>
<mn>31</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Q</mi>
<mn>32</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Q</mi>
<mn>33</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
8.根据权利要求1所述的一种用于长基线水声定位系统性能评估的方法,其特征在于,步骤S6中,采用下列公式计算所述指定网格点的评估参数GDOP、HDOP和VDOP值:
<mrow>
<mi>G</mi>
<mi>D</mi>
<mi>O</mi>
<mi>P</mi>
<mo>=</mo>
<msqrt>
<mrow>
<msub>
<mi>Q</mi>
<mn>11</mn>
</msub>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mn>22</mn>
</msub>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mn>33</mn>
</msub>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>H</mi>
<mi>D</mi>
<mi>O</mi>
<mi>P</mi>
<mo>=</mo>
<msqrt>
<mrow>
<msub>
<mi>Q</mi>
<mn>11</mn>
</msub>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mn>22</mn>
</msub>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>V</mi>
<mi>D</mi>
<mi>O</mi>
<mi>P</mi>
<mo>=</mo>
<msqrt>
<msub>
<mi>Q</mi>
<mn>33</mn>
</msub>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
9.根据权利要求1所述的一种用于长基线水声定位系统性能评估的方法,其特征在于,所述声信标三维位置与指定网格点三维位置不处于同一平面时,选择评估参数GDOP、HDOP及VDOP对长基线水声定位系统声信标布局进行评估;所述声信标三维位置与指定网格点三维位置处于同一平面时,只选择评估参数HDOP对长基线水声定位系统声信标布局进行评估。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711260848.8A CN108120955A (zh) | 2017-12-04 | 2017-12-04 | 一种用于长基线水声定位系统性能评估的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711260848.8A CN108120955A (zh) | 2017-12-04 | 2017-12-04 | 一种用于长基线水声定位系统性能评估的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108120955A true CN108120955A (zh) | 2018-06-05 |
Family
ID=62229663
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711260848.8A Pending CN108120955A (zh) | 2017-12-04 | 2017-12-04 | 一种用于长基线水声定位系统性能评估的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108120955A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109116361A (zh) * | 2018-08-03 | 2019-01-01 | 中国人民解放军91388部队 | 一种抗距离模糊水声合作定位信号 |
CN109814110A (zh) * | 2019-02-21 | 2019-05-28 | 哈尔滨工程大学 | 深海长基线定位阵形拓扑结构的布阵方法 |
CN111443329A (zh) * | 2020-03-25 | 2020-07-24 | 北京东方振动和噪声技术研究所 | 声源定位方法及装置、计算机存储介质和电子设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106199519A (zh) * | 2016-07-26 | 2016-12-07 | 山东省科学院海洋仪器仪表研究所 | 一种超短基线五基元立体空间基阵及其水声定位方法 |
CN106483501A (zh) * | 2015-09-01 | 2017-03-08 | 北京自动化控制设备研究所 | 一种基于dop值分析声学定位系统多应答器最优布阵方法 |
-
2017
- 2017-12-04 CN CN201711260848.8A patent/CN108120955A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106483501A (zh) * | 2015-09-01 | 2017-03-08 | 北京自动化控制设备研究所 | 一种基于dop值分析声学定位系统多应答器最优布阵方法 |
CN106199519A (zh) * | 2016-07-26 | 2016-12-07 | 山东省科学院海洋仪器仪表研究所 | 一种超短基线五基元立体空间基阵及其水声定位方法 |
Non-Patent Citations (3)
Title |
---|
孙万卿: "浅海水声定位技术及应用研究", 《中国博士学位论文全文数据库 工程科技II辑》 * |
张文照 等: "基于遗传算法的潜艇自航模水声定位系统定位站布置研究", 《船舶力学》 * |
杨日杰 等: "主动全向声纳浮标跟踪潜艇优化布放方法", 《系统工程与电子技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109116361A (zh) * | 2018-08-03 | 2019-01-01 | 中国人民解放军91388部队 | 一种抗距离模糊水声合作定位信号 |
CN109814110A (zh) * | 2019-02-21 | 2019-05-28 | 哈尔滨工程大学 | 深海长基线定位阵形拓扑结构的布阵方法 |
CN111443329A (zh) * | 2020-03-25 | 2020-07-24 | 北京东方振动和噪声技术研究所 | 声源定位方法及装置、计算机存储介质和电子设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Seafloor geodetic network establishment and key technologies | |
CN102495420B (zh) | 一种水下目标精确定位系统及方法 | |
CN103412198B (zh) | 船舶防护电场的三维空间分布特性测量装置及测量方法 | |
CN105157701B (zh) | 一种水下位置基准布设的方法 | |
CN105547290B (zh) | 一种基于超短基线定位系统的从潜器导航方法 | |
CN107390177B (zh) | 一种基于纯测向的被动水下声学定位方法 | |
CN107664758B (zh) | 基于长基线或超短基线组网的深海导航定位系统及方法 | |
CN103868493B (zh) | 一种基于ppp技术的深度基准面大地高测量方法 | |
CN110941017A (zh) | 基于磁矢量数据的海缆三维路由测量方法及测量仪 | |
CN105182390B (zh) | 一种运载体水下定位的方法 | |
CN108120955A (zh) | 一种用于长基线水声定位系统性能评估的方法 | |
CN105258684A (zh) | 基于激光点云为约束的多波束低掠射角波束归位方法 | |
CN102692217A (zh) | 运用地效应飞行器进行河道测量的方法 | |
CN107576939A (zh) | 一种基于虚拟测距信标的单信标测距定位方法 | |
CN101441266B (zh) | 水下多应答器组合导航方法 | |
CN104049241B (zh) | 目标位置坐标未知的双基地合成孔径雷达的空间同步方法 | |
CN103792595B (zh) | 用于定位水下运动磁目标的磁传感器阵列校准方法 | |
CN105738869B (zh) | 一种适用于单水听器的深水信标搜索定位方法 | |
CN104484891A (zh) | 一种基于纹理特征和地形特征参数的水下地形匹配方法 | |
CN109490927A (zh) | 一种水下整平架定位系统及其定位方法 | |
CN110441736B (zh) | 多关节水下无人潜行器变基线三维空间定位方法 | |
CN109059746A (zh) | 一种基于精密pos的海底地形测量方法 | |
CN107462891A (zh) | 一种三点式深海潜标定位方法 | |
CN101937106B (zh) | 一种海底大起伏测线磁测数据的处理方法 | |
Grządziel et al. | Estimation of effective swath width for dual-head multibeam echosounder |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180605 |
|
RJ01 | Rejection of invention patent application after publication |