CN112986913A - 一种基于差分多普勒与到达时延差的水下目标定位方法 - Google Patents
一种基于差分多普勒与到达时延差的水下目标定位方法 Download PDFInfo
- Publication number
- CN112986913A CN112986913A CN202110161121.4A CN202110161121A CN112986913A CN 112986913 A CN112986913 A CN 112986913A CN 202110161121 A CN202110161121 A CN 202110161121A CN 112986913 A CN112986913 A CN 112986913A
- Authority
- CN
- China
- Prior art keywords
- receiving node
- target
- underwater
- vector
- node
- 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
Images
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
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/50—Systems of measurement, based on relative movement of the target
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明属于水下目标定位技术领域,涉及一种利用差分多普勒与到达时延差的水下目标定位方法,包括以下步骤:第一步,实际测量获得各个接收节点与随机编号的第1个接收节点的水下目标DDS和TDOA;第二步,利用TDOA计算得到第i个接收节点与第1个接收节点之间的距离差ri1;第三步,将第一步获得的DDS和第二步获得的距离差表示为矩阵方程的形式,然后估计得到目标位置坐标和速度向量;第四步,获得精度更高的目标位置坐标和速度向量估计值,以及时间同步因子的估计值。本发明的有益效果:建立更加实际的定位节点时钟模型,通过算法实现测距节点之间时钟偏离的有效补偿,降低了定位节点硬件设计的复杂度,获得了相比于传统的TDOA方法更高的水下目标定位精度。
Description
技术领域
本发明属于水下目标定位技术领域,涉及一种时钟不同步情况下利用差分多普勒与到达时延差的水下目标定位方法。
背景技术
水下目标定位技术是一个一直以来备受关注的研究内容。通常来说,水下目标定位是通过若干个已知位置的接收节点来接收水下目标发出的噪声信号,然后利用接收到的信息综合推断目标在水下的位置。目前较为常用的一种水下目标定位方法是基于各节点的到达时延(time-of-arrival,TOA)估计方法,该方法利用各节点接收到目标信号的时延来对目标位置进行解算并定位(具体步骤详见文献:N.Patwari,et.al.,“Locating thenodes:cooperative localization in wireless sensor networks,”IEEE SignalProcess.Mag.,vol.22,no.4,pp.54-69,Jul.2005.)。但是该方法要求各节点准确地知道目标信号的发出时刻(这种条件称为“时钟同步”),否则无法准确测量目标信号传输到各节点的时延,严重影响目标定位精度(通常采用作为目标位置估计值与真实值的均方根误差作为指标来衡量定位精度)。事实上,绝大多数水下目标定位的实际情况中都存在时钟不同步的问题。Liang等人指出采用到达时延差(time-difference-of-arrival,TDOA)估计方法可以较为有效地解决这一问题,因为该方法计算的是每个节点与某一个参考节点接收到目标信号的时延差,从而不会受到时钟同步的约束(具体步骤详见文献:Q.Liang,et.al.,“TDoAfor passive localization:Underwater versus terrestrial environment,”IEEETrans.Parallel Distrib.Syst.,vol.24,no.10,pp.2100-2108,Oct.2013.)。但是,该方法的劣势在于仅仅利用到达时延差观测量获得的目标定位精度不尽如人意;此外,该方法跳过了对时钟不同步程度的估计过程,因而无法获得对时钟不同步程度的估计值。由于接收节点接收到的目标信号不仅仅含有时延信息,还含有多普勒、相位等信息,充分利用这些信息进行信号处理有望获得更高的水下目标定位精度。
发明内容
为了进一步提高水下目标的定位精度,本发明在传统的TDOA方法基础上,引入水下目标的差分多普勒(differential Doppler scale,DDS,即各个接收节点之间的多普勒差值)这一参数,提出一种时钟不同步情况下利用差分多普勒与到达时延差的水下目标定位方法。
本发明的技术方案包括以下步骤:
第一步,实际测量获得各个接收节点与随机编号的第1个接收节点的水下目标DDS和TDOA:
水下目标发射已知时宽和带宽的合作信号由各接收节点接收,各接收节点不重复地随机编号为1~N,N表示接收节点数目;
其中i=2,3,...,N,和分别表示第i个接收节点接收到的目标多普勒和时延,是第i个接收节点直接接收信号计算获得的观测量,和εi1分别表示叠加在和τi1上的独立高斯白噪声,是第i个接收节点直接接收背景噪声计算获得的观测量,和εi1的单位分别为“m/s”和“s”。ω表示时间同步因子,反映了时钟不同步程度;ω=1表示各接收节点本地时钟与目标节点本地时钟频率完全一致,即时间同步;时间不同步时ω会在1附近正或负偏离,正或负偏离越大则时间越不同步,由于ω的值在实际中较难通过接收节点测量获得,第四步中将给出ω的计算方法;和为第一个接收节点接收到的目标多普勒和时延。
第二步,利用TDOA计算得到第i个接收节点与第1个接收节点之间的距离差ri1:
其中,ri o表示水下目标在探测场景中与第i个节点的距离,r1 o表示水下目标在探测场景中与第1个接收节点的距离;
第三步,将第一步获得的DDS和第二步获得的距离差表示为矩阵方程的形式,然后估计得到目标位置坐标和速度向量:
令ri o=||u-si||,其中u表示水下目标在探测场景中的三维坐标,si表示第i个节点在探测场景中的三维坐标(下标为1则表示第1个节点在探测场景中的三维坐标,其他类似,不再单独定义)。和分别为u和si的微分,分别表示水下目标和第i个节点的三维速度向量。为ri o的微分,表示第i个接收节点接收到的目标多普勒,两者的关系式为然后将(1)和(2)式进行两边取平方,并遍历所有i值后将(1)和(2)式合并表示为如下矩阵形式:
e1=h1-G1φ1 (4)
公式(4)中的φ1为一个具有10个元素的列向量,记为
其中上标T表示向量转置。
对公式(4)利用最小二乘法(具体步骤详见文献:S.M.Kay,Fundamentals ofStatistical Signal Processing,Estimation Theory.Upper Saddle River,NJ,USA:Prentice-Hall,1993.),得到φ1的估计向量
第四步,对第三步获得的目标位置坐标和速度向量估计值进行进一步迭代计算,获得精度更高的目标位置坐标和速度向量估计值,以及时间同步因子的估计值:
从前面φ1的定义中可以看出,φ1向量的各个元素并非独立的,部分元素之间存在关联性。考虑φ1向量中元素之间的关联性可构建如下矩阵方程:
e2=h2-G2φ2 (6)
公式(6)中的e2为e1通过迭代计算后的高斯白噪声组合矩阵,记为
e2=B2Δφ1
同样,对公式(6)利用最小二乘法(具体步骤详见文献:S.M.Kay,Fundamentals ofStatistical Signal Processing,Estimation Theory.Upper Saddle River,NJ,USA:Prentice-Hall,1993.),得到φ2的估计向量
至此,精度更高的目标位置坐标和速度向量估计值,以及时间同步因子的估计值分别表示为和(水下目标定位的目的是获得目标位置坐标而和两个估计结果是求解目标位置的过程中获得的附加结果。但是,代表目标的速度向量,则反映了目标节点和接收节点之间的时钟同步情况,为我们提供了更加丰富的目标参数信息,对后续目标处理(如跟踪)是有益的。)
本发明的有益效果是:
1、相比传统的TOA和TDOA方法,本发明方法计入了时间同步因子,建立更加实际的定位节点时钟模型,通过算法实现测距节点之间时钟偏离的有效补偿。
2、本发明方法具备时钟修正能力,可实现测距节点之间时钟不同步情况下的目标定位,降低了定位节点硬件设计的复杂度,使本发明方法对探测环境的适应性更强。
3、利用水下目标的DDS这一信息,获得了相比于传统的TDOA方法更高的水下目标定位精度,使得定位结果更加可信。
附图说明
图1是本发明的实现流程图;
图2是利用本发明的具体实施方式进行实测数据试验的结果:(a)为时间同步因子RMSE(root mean squared error,RMSE,均方根误差)随TDOA上叠加的独立高斯白噪声标准差σ的变化情况,(b)为本发明方法与时间同步的加权最小二乘方法的水下目标位置坐标RMSE随σ的变化情况的对比结果,(c)为本发明方法与时间同步的加权最小二乘方法的水下目标速度向量RMSE随σ的变化情况的对比结果。
具体实施方式
图1是本发明的实现流程图。具体实施方式包括以下步骤:
第一步,实际测量获得各个接收节点与随机编号的第1个接收节点的水下目标DDS和TDOA。
第二步,利用TDOA计算得到第i个接收节点与第1个接收节点之间的距离差ri1。
第三步,将第一步获得的DDS和第二步获得的距离差表示为矩阵方程的形式,然后估计得到目标位置坐标和速度向量。
第四步,对第三步获得的目标位置坐标和速度向量估计值进行进一步迭代计算,获得精度更高的目标位置坐标和速度向量估计值,以及时间同步因子的估计值。
图2是利用本发明的具体实施方式进行数值仿真实验的结果。仿真实验场景设计如下:场景中存在1个水下目标和10个接收节点,其位置坐标均匀随机地分布于[0~1000,0~1000,0~1000]m这一区域;水中有效声速c设为1500m/s;水下目标和10个接收节点的运动速度为-3~3m/s之间的随机值;时间同步因子ω为0.995~1.005之间的随机值;DDS和TDOA上叠加的独立高斯白噪声和εi1的方差分别为0.01σ2和σ2,其单位分别为“(m/s)2”和“s2”。通过2000次蒙特卡洛仿真,统计得到不同σ情况下时间同步因子RMSE、水下目标位置坐标和速度向量RMSE如图2所示。图中横坐标为σ,单位为“1”;纵坐标从图2(a)至图2(c)分别为时间同步因子RMSE,单位为“1”;水下目标位置坐标RMSE,单位为“m”;水下目标速度向量RMSE,单位为“m/s”。用于对比的方法为基于TDOA的时间同步(即ω=1)的加权最小二乘方法(具体步骤详见文献:K.C.Ho and W.Xu,“An accurate algebraic solution formoving source location using TDOA and FDOA measurements,”IEEE Trans.SignalProcess.,vol.52,no.9,pp.2453-2463,Sep.2004.)。对比结果发现,本发明所提方法可以较为准确地获得时间同步因子的估计值,而现有方法无法对时间同步因子进行估计。另一方面,本发明所提方法的水下目标位置坐标和速度向量RMSE值在不同σ情况下均比时间同步的加权最小二乘方法更小,即本发明方法估计得到的水下目标位置坐标和速度向量更加接近目标的真实情况,表明本发明方法具有相比于现有方法更高的水下目标定位精度。
Claims (1)
1.一种基于差分多普勒与到达时延差的水下目标定位方法,其特征在于,该方法包括以下步骤:
第一步,实际测量获得各个接收节点与随机编号的第1个接收节点的水下目标DDS和TDOA:
水下目标发射已知时宽和带宽的合作信号由各接收节点接收,各接收节点不重复地随机编号为1~N,N表示接收节点数目;
其中i=2,3,...,N,和分别表示第i个接收节点接收到的目标多普勒和时延,是第i个接收节点直接接收信号计算获得的观测量,和εi1分别表示叠加在和τi1上的独立高斯白噪声,是第i个接收节点直接接收背景噪声计算获得的观测量,和εi1的单位分别为“m/s”和“s”;ω表示时间同步因子,反映了时钟不同步程度;和为第一个接收节点接收到的目标多普勒和时延;
第二步,利用TDOA计算得到第i个接收节点与第1个接收节点之间的距离差ri1:
其中,ri o表示水下目标在探测场景中与第i个节点的距离,r1 o表示水下目标在探测场景中与第1个接收节点的距离;
第三步,将第一步获得的DDS和第二步获得的距离差表示为矩阵方程的形式,然后估计得到目标位置坐标和速度向量:
令ri o=||u-si||,其中u表示水下目标在探测场景中的三维坐标,si表示第i个节点在探测场景中的三维坐标;和分别为u和si的微分,分别表示水下目标和第i个节点的三维速度向量;为ri o的微分,表示第i个接收节点接收到的目标多普勒,两者的关系式为然后将(1)和(2)式进行两边取平方,并遍历所有i值后将(1)和(2)式合并表示为如下矩阵形式:
e1=h1-G1φ1 (4)
φ1为一个具有10个元素的列向量,记为
其中上标T表示向量转置;
第四步,对第三步获得的目标位置坐标和速度向量估计值进行进一步迭代计算,获得精度更高的目标位置坐标和速度向量估计值,以及时间同步因子的估计值:
考虑φ1向量中元素之间的关联性可构建如下矩阵方程:
e2=h2-G2φ2 (6)
e2为e1通过迭代计算后的高斯白噪声组合矩阵,记为
e2=B2Δφ1
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110161121.4A CN112986913A (zh) | 2021-02-05 | 2021-02-05 | 一种基于差分多普勒与到达时延差的水下目标定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110161121.4A CN112986913A (zh) | 2021-02-05 | 2021-02-05 | 一种基于差分多普勒与到达时延差的水下目标定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112986913A true CN112986913A (zh) | 2021-06-18 |
Family
ID=76348017
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110161121.4A Pending CN112986913A (zh) | 2021-02-05 | 2021-02-05 | 一种基于差分多普勒与到达时延差的水下目标定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112986913A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114924229A (zh) * | 2022-05-18 | 2022-08-19 | 中国人民解放军国防科技大学 | 一种时钟不同步的时延-多普勒域水下目标差分定位方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103969622A (zh) * | 2014-04-25 | 2014-08-06 | 西安电子科技大学 | 一种基于多运动接收站的时差定位方法 |
CN106707234A (zh) * | 2016-12-16 | 2017-05-24 | 西北工业大学 | 一种联合时延差与角度测量的传感器网络目标定位方法 |
CN108051779A (zh) * | 2017-10-13 | 2018-05-18 | 西安电子科技大学 | 一种面向tdoa的定位节点优选方法 |
CN109581281A (zh) * | 2018-11-08 | 2019-04-05 | 宁波大学 | 基于到达时间差和到达频率差的移动目标定位方法 |
-
2021
- 2021-02-05 CN CN202110161121.4A patent/CN112986913A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103969622A (zh) * | 2014-04-25 | 2014-08-06 | 西安电子科技大学 | 一种基于多运动接收站的时差定位方法 |
CN106707234A (zh) * | 2016-12-16 | 2017-05-24 | 西北工业大学 | 一种联合时延差与角度测量的传感器网络目标定位方法 |
CN108051779A (zh) * | 2017-10-13 | 2018-05-18 | 西安电子科技大学 | 一种面向tdoa的定位节点优选方法 |
CN109581281A (zh) * | 2018-11-08 | 2019-04-05 | 宁波大学 | 基于到达时间差和到达频率差的移动目标定位方法 |
Non-Patent Citations (1)
Title |
---|
张兵兵: "海洋环境下的目标定位技术研究", 《中国优秀博硕士学位论文全文数据库(博士) 基础科学辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114924229A (zh) * | 2022-05-18 | 2022-08-19 | 中国人民解放军国防科技大学 | 一种时钟不同步的时延-多普勒域水下目标差分定位方法 |
CN114924229B (zh) * | 2022-05-18 | 2024-04-09 | 中国人民解放军国防科技大学 | 一种时钟不同步的时延-多普勒域水下目标差分定位方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106842128B (zh) | 运动目标的声学跟踪方法及装置 | |
CN106405533B (zh) | 基于约束加权最小二乘的雷达目标联合同步与定位方法 | |
CN114924229B (zh) | 一种时钟不同步的时延-多普勒域水下目标差分定位方法 | |
CN106093849B (zh) | 一种基于测距和神经网络算法的水下定位方法 | |
CN111157943B (zh) | 在异步网络中基于toa的传感器位置误差抑制方法 | |
CN105353351A (zh) | 一种基于多信标到达时间差改进型定位方法 | |
Poursheikhali et al. | TDOA based target localization in inhomogenous underwater wireless sensor network | |
WO2016112758A1 (zh) | 终端的定位方法及装置 | |
CN112986907B (zh) | 一种时钟偏差和时钟漂移条件下的运动目标定位方法 | |
Wang et al. | TOA-based joint synchronization and source localization with random errors in sensor positions and sensor clock biases | |
KR101627419B1 (ko) | 이동 노드의 위치 추정 방법 및 그 장치 | |
CN107765216B (zh) | 非同步无线网络中的目标位置和时钟参数联合估计方法 | |
CN112986913A (zh) | 一种基于差分多普勒与到达时延差的水下目标定位方法 | |
White et al. | Localisation of sperm whales using bottom-mounted sensors | |
CN111770528B (zh) | 基于信道参数萃取方法的视距与非视距识别方法及装置 | |
Stefanski | Asynchronous wide area multilateration system | |
KR101480834B1 (ko) | 다중 경로 음파 전달 모델 및 표적 식별을 이용한 표적 기동분석 방법 | |
Jespersen et al. | An indoor multipath-assisted single-anchor UWB localization method | |
Spiesberger | Geometry of locating sounds from differences in travel time: Isodiachrons | |
CN108111440A (zh) | 一种ofdm的子信道估计测距方法 | |
JP2570110B2 (ja) | 水中音源位置推定システム | |
Nikoo et al. | A disambiguation technique for passive localization using trajectory analysis | |
CN111337881B (zh) | 一种利用螺旋桨噪声的水下目标探测方法 | |
Ricks et al. | Passive acoustic tracking for cueing a multistatic active acoustic tracking system | |
Khalaf-Allah | An extended closed-form least-squares solution for three-dimensional hyperbolic geolocation |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20210618 |