CN105487046A - 大入射角声线跟踪定位方法 - Google Patents
大入射角声线跟踪定位方法 Download PDFInfo
- Publication number
- CN105487046A CN105487046A CN201510801428.0A CN201510801428A CN105487046A CN 105487046 A CN105487046 A CN 105487046A CN 201510801428 A CN201510801428 A CN 201510801428A CN 105487046 A CN105487046 A CN 105487046A
- Authority
- CN
- China
- Prior art keywords
- layer
- rsqb
- lsqb
- sound
- sound ray
- 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
- G01S1/00—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith
- G01S1/72—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith using ultrasonic, sonic or infrasonic waves
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种水下声学定位方法,具体为基于分层等梯度声线跟踪的大入射角声线跟踪定位方法,根据分层等梯度假设中第一层曲率半径的公式用第一层曲率半径R来表示Snell常数p;将声波在水层中水平位移公式中的Snell常数p用曲率半径R表示;对新建立的公式用弦截法对R进行迭代;将迭代结果转换为入射角;用所求得的入射角进行声线跟踪。本发明提供的大入射角声线跟踪定位方法,将传统的分层等梯度声线跟踪进行了改进,分析了大入射角声线情况下定位失败的原因并进行了改进,既达到了扩大初始入射角迭代的最大范围,又提高了定位的效率。
Description
技术领域
本发明涉及一种水下声学定位方法,具体为基于分层等梯度声线跟踪的大入射角声线跟踪定位方法。
背景技术
为获取水下高精度位置信息,现常采用声线跟踪的方法,即跟踪声波在换能器与应答器之间的传播路径,利用最小二乘法解算水下应答器的坐标。其中,分层等梯度声线跟踪方法采用逐层计算获取的声线轨迹与真实声线轨迹最为接近,精度较高。但当换能器与应答器组成的观测几何结构不好时,例如当水面的换能器距离应答器较远时,声线入射角会增大,此时若采用分层等梯度声线跟踪方法计算,会导致入射角迭代发散,无法得到正确结果。
发明内容
针对上述技术问题,本发明提供一种基于分层等梯度声线跟踪的大入射角声线跟踪定位方法,方法能够有效扩大入射角解算范围,并提高定位效率。
本发明的目的可以通过以下技术方案来实现:
大入射角声线跟踪定位方法,包括以下过程:
假设声波经历的水柱分为N个等梯度层,即声波在每层的传播速度是等梯度的,在水深为zi的i层,分别用Ci和θi表示声波的传播速度和入射角度,第i层内的声速梯度gi为:
gi=(Ci+1-Ci)/Δzi
式中,Δzi为第i层的水层厚度,Ci+1为第i+1层声速;
声波波束在第i层内的实际传播轨迹为一连续的、对应一定曲率半径Ri的弧段,Ri的表达式为:
第i层内声线的水平位移为:
式中,θi为第i层声线的入射角,θi+1为第i+1层声线的入射角;
波束在该层经历的弧段长度为:
Si=Ri(θi-θi+1);
经历该段的时间为:
p为Snell常数,转换成曲率半径R进行迭代计算,p与R的关系为:
用第一层的曲率半径R表示p:
声线水平位移为:
用弦截法迭代求解曲率半径R的方程:
将迭代计算得到的R换算成初始入射角,计算声波传播时间,与实测传播时间进行对比,列立误差方程并用最小二乘方法进行vcx迭代求解。
本发明提供的大入射角声线跟踪定位方法,将传统的分层等梯度声线跟踪进行了改进,分析了大入射角声线情况下定位失败的原因并进行了改进,既达到了扩大初始入射角迭代的最大范围,又提高了定位的效率。
附图说明
图1为本发明的方法流程图;
图2为实施例的海面浮标位置数据与海底应答器位置示意图;
图3为实施例采用的中国南海声速剖面数据。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
根据模拟浮标位置对比原始定位方法作为对比例,与本方法的大入射角声线跟踪定位方法同时进行对比:
(1)定位数据
模拟四个浮标布设成正方形形状,浮标间距分别模拟为2200m、2400m和4000m。给浮标的三维坐标分别加上2米的余弦波动。如图2,海底模拟4个应答器A、B、C、D,坐标分别为(0,0,-500),(0,40,-500),(40,0,-500),(40,40,-500),单位均为m。总共模拟100个历元。水域模拟采用如图3所示的中国南海声速剖面数据,采样间隔为1m。误差模拟如下:换能器定位中误差x、y、z方向各为10cm,应答器时延偏差为10cm,声速测量误差为0.15m/s,浮标位置随机误差为10cm。
(2)浮标间距为2200m时两种方法定位结果,见表1:
表1浮标间距为2200m时两种方法定位结果
在入射角不太大的情况下,大入射角声线跟踪方法与原方法定位精度相当,但效率较原方法有所提升。
3、浮标间距为2400m时两种方法定位结果见表2:
表2浮标间距为2400m时两种方法定位结果
随着浮标间距的扩大,入射角的增大,原始定位方法失效,大入射角声线跟踪方法仍能有效定位。
4、浮标间距为4000m时两种方法定位结果见表3:
表3浮标间距为4000m时两种方法定位结果
针对此声速剖面数据,大入射角声线跟踪方法将最大入射角的解算范围扩大到88.8°。
Claims (1)
1.大入射角声线跟踪定位方法,其特征在于,包括以下步骤:
假设声波经历的水柱分为N个等梯度层,声波在每层的传播速度是等梯度的,在水深为zi的i层,分别用Ci和θi表示声波的传播速度和入射角度,第i层内的声速梯度gi为:
gi=(Ci+1-Ci)/Δzi
式中,Δzi为第i层的水层厚度,Ci+1为第i+1层声速;
声波波束在第i层内的实际传播轨迹为一连续的、对应一定曲率半径Ri的弧段,Ri的表达式为:
第i层内声线的水平位移为:
式中,θi为第i层声线的入射角,θi+1为第i+1层声线的入射角;
波束在该层经历的弧段长度为:
Si=Ri(θi-θi+1);
经历该段的时间为:
p为Snell常数,转换成曲率半径R进行迭代计算,p与R的关系为:
用第一层的曲率半径R表示p:
声线水平位移为:
用弦截法迭代求解曲率半径R的方程:
将迭代计算得到的R换算成初始入射角,计算声波传播时间,与实测传播时间进行对比,列立误差方程并用最小二乘方法进行vcx迭代求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510801428.0A CN105487046A (zh) | 2015-11-19 | 2015-11-19 | 大入射角声线跟踪定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510801428.0A CN105487046A (zh) | 2015-11-19 | 2015-11-19 | 大入射角声线跟踪定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105487046A true CN105487046A (zh) | 2016-04-13 |
Family
ID=55674155
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510801428.0A Pending CN105487046A (zh) | 2015-11-19 | 2015-11-19 | 大入射角声线跟踪定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105487046A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110297250A (zh) * | 2019-06-17 | 2019-10-01 | 东南大学 | 基于泰勒展开的初始掠射角求解方法、声线弯曲修正方法和设备 |
WO2020228547A1 (zh) * | 2019-05-16 | 2020-11-19 | 浙江大学 | 基于倒置式多波束回声仪的声速剖面反演方法 |
CN113406645A (zh) * | 2021-05-10 | 2021-09-17 | 山东科技大学 | 一种新的平均声速水下声呐定位方法 |
CN113761462A (zh) * | 2021-09-10 | 2021-12-07 | 山东大学 | 基于弦截法的初始入射角迭代计算改进方法 |
CN117252035A (zh) * | 2023-11-14 | 2023-12-19 | 天津大学 | 水面航行器入射角的确定方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101806884A (zh) * | 2010-04-23 | 2010-08-18 | 哈尔滨工程大学 | 基于超短基线的深海信标绝对位置精确定位方法 |
CN104613906A (zh) * | 2015-02-06 | 2015-05-13 | 长江水利委员会水文局 | 基于声线跟踪的库区深水水深测量方法 |
-
2015
- 2015-11-19 CN CN201510801428.0A patent/CN105487046A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101806884A (zh) * | 2010-04-23 | 2010-08-18 | 哈尔滨工程大学 | 基于超短基线的深海信标绝对位置精确定位方法 |
CN104613906A (zh) * | 2015-02-06 | 2015-05-13 | 长江水利委员会水文局 | 基于声线跟踪的库区深水水深测量方法 |
Non-Patent Citations (2)
Title |
---|
张红梅 等编著: "《水下导航定位技术》", 30 November 2010, 武汉大学出版社 * |
李圣雪 等: "一种适用于深海长基线定位的自适应分层声线跟踪法", 《海洋通报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020228547A1 (zh) * | 2019-05-16 | 2020-11-19 | 浙江大学 | 基于倒置式多波束回声仪的声速剖面反演方法 |
US11733381B2 (en) | 2019-05-16 | 2023-08-22 | Zhejiang University | Sound velocity profile inversion method based on inverted multi-beam echo sounder |
CN110297250A (zh) * | 2019-06-17 | 2019-10-01 | 东南大学 | 基于泰勒展开的初始掠射角求解方法、声线弯曲修正方法和设备 |
CN113406645A (zh) * | 2021-05-10 | 2021-09-17 | 山东科技大学 | 一种新的平均声速水下声呐定位方法 |
CN113761462A (zh) * | 2021-09-10 | 2021-12-07 | 山东大学 | 基于弦截法的初始入射角迭代计算改进方法 |
CN117252035A (zh) * | 2023-11-14 | 2023-12-19 | 天津大学 | 水面航行器入射角的确定方法 |
CN117252035B (zh) * | 2023-11-14 | 2024-02-13 | 天津大学 | 水面航行器入射角的确定方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106886024B (zh) | 深海多波束声线精确跟踪方法 | |
CN105487046A (zh) | 大入射角声线跟踪定位方法 | |
CN105388457B (zh) | 一种基于等效声速梯度的长基线水声定位方法 | |
CN104237849B (zh) | 一种双五元十字阵被动声定位融合方法 | |
CN105629220B (zh) | 一种基于单水听器的深海水声被动测距方法 | |
CN101644759B (zh) | 一种基于二元结构的水下目标定位***及定位方法 | |
CN110399680B (zh) | 一种浅海弹性结构辐射声场计算方法 | |
CN102749622B (zh) | 基于多波束测深的声速剖面及海底地形的联合反演方法 | |
CN104678384B (zh) | 一种波束域的声压差互相关谱分析水下目标速度估计方法 | |
CN112505666B (zh) | 一种二维多波束合成孔径声呐目标回波仿真方法 | |
CN101806884A (zh) | 基于超短基线的深海信标绝对位置精确定位方法 | |
CN105547290B (zh) | 一种基于超短基线定位***的从潜器导航方法 | |
CN101086534A (zh) | 检波点二次定位方法 | |
CN102262226A (zh) | 基于差分技术的水下定位方法 | |
CN112540348A (zh) | 一种基于空间尺度上的声线修正算法在长基线水声定位***上的应用 | |
CN105445694A (zh) | 等梯度声线跟踪的声速剖面自适应分层方法 | |
CN105182347B (zh) | 一种基于小孔径基阵的高精度水下目标被动定位方法 | |
US10677919B2 (en) | Self-positioning system and self-positioning identification method of deepwater underwater robot of irregular dam surface of reservoir | |
CN113671443A (zh) | 基于掠射角声线修正的水声传感器网络深海目标定位方法 | |
CN109738902A (zh) | 一种基于同步信标模式的水下高速目标高精度自主声学导航方法 | |
CN104678359A (zh) | 一种声场识别的多孔声全息方法 | |
CN110780263A (zh) | 一种基于卡西尼卵形线的多基地声纳***定位精度分析方法 | |
JP6207817B2 (ja) | 水中位置関係情報取得システム | |
CN110441736A (zh) | 多关节水下无人潜行器变基线三维空间定位方法 | |
CN110309581A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160413 |