CN102096064B - 一种用于短基线水声定位系统中精确测量时延差的方法及其系统 - Google Patents
一种用于短基线水声定位系统中精确测量时延差的方法及其系统 Download PDFInfo
- Publication number
- CN102096064B CN102096064B CN 201010540329 CN201010540329A CN102096064B CN 102096064 B CN102096064 B CN 102096064B CN 201010540329 CN201010540329 CN 201010540329 CN 201010540329 A CN201010540329 A CN 201010540329A CN 102096064 B CN102096064 B CN 102096064B
- Authority
- CN
- China
- Prior art keywords
- signal
- correlator
- positioning system
- calculation module
- delay inequality
- 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
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明所述的用于短基线水声定位系统中精确测量时延差的方法为信号产生在信号处理单元1中,通过发射信号样本存储器、发射信号控制模块和发射信号生成模块实现;信号解算在信号处理单元2中,通过参考信号样本存储器、相关器1、相关器2、包络计算模块1和包络计算模块2、峰值检测器1、峰值检测器2、时延差初测值计算模块、相位计算模块1、相位计算模块2、减法器1、时延差修正值计算模块和减法器2实现。本发明所述的测量并计算多普勒偏移方法的系统包括接收端和发射端,发射端包括信号处理单元1、功放、发射换能器,接收端包括接收换能器1、接收换能器2、滤波器1、滤波器2和信号处理单元2。
Description
技术领域:
本发明提出了一种用于短基线水声定位系统中精确测量时延差的方法和系统,以应用于对定位精度要求较高的短基线水声定位系统,为定位解算提供高精度的时延差测量值。
背景技术:
短基线水声定位系统广泛用于导航,水下测量,海洋勘探和开发等众多领域。随着对海洋探索的不断深入,人们对定位精度的要求日益提高。
根据短基线定位的原理:短基线水声定位系统的定位精度与定位基阵的基线上两接收换能器之间的距离(即基线的长度)成正比;定位精度与对定位基阵的基线上两接收换能器接收到声波时延差的测量精度成正比。通常,人们希望短基线水声定位系统基阵的尽可能的小,以便节省空间;而一些设备对定位基阵的尺寸有着严格的要求。因而,提高时延差测量的精度值,就成了提高短基线定位系统定位精度的必然要求。
发明内容:
本发明的目的是研究出一种用于短基线水声定位系统中精确测量时延差的方法及其系统。
本发明针对短基线水声定位系统中希望在基线长度较短条件下有较高定位精度的要求,给出了一种利用相位差计算时延差补偿值的方法和系统,以提高测量时延差的精度。
本发明所述的用于短基线水声定位系统中精确测量时延差的方法,至少包括以下步骤:
1)短基线水声定位系统的发射端发射线性调频脉冲信号作为定位信号,形式如下:
其中,TS1是短基线水声定位系统发射信号的采样间隔,T是线性调频信号由最低频扫到最高频所经过的时间,k0是信号的调频斜率,fL是信号s(n)的最低频率;
2)短基线水声定位系统的接收端有相关器1和相关器2,相关器1接收信号x1(n)并进行相关处理,相关器2接收信号x2(n)并进行相关处理;相关器1和相关器2进行相关处理时,使用的参考信号相同,参考信号形式如下:
其中,TS2是短基线水声定位系统接收信号的采样间隔;
设处理后,从相关器1输出的输出信号R1(m)和从相关器2输出的输出信号R2(m)的包络出现峰值对应的位置为M1和M2,那么时延差初测值为:
根据峰值位置M1和M2,由相位计算模块1生成的生成信号p1(m)和相位计算模块2生成的生成信号p2(m)分别为:
p1(m)=R1(m-M1-L)0≤m≤2L-1
p2(m)=R2(m-M2-L)0≤m≤2L-1
设存在W个非负整数满足:
则时延差修正值为:
那么,时延差初测值和时延差修正值的差值,即为时延差值
Δτ=τR-τM
参数Δτ的物理意义是,在发射端发出的信号经过信道传播后,接收端两个不同的通道接收到信号的时刻相差了Δτ。
本发明所述的用于短基线水声定位系统中精确测量时延差方法的系统包括:接收端和发射端,发射端包括信号处理单元1,功放,发射换能器,接收端包括接收换能器1,接收换能器2,滤波器1,滤波器2,信号处理单元2;信号处理单元1与功放之间通过数据线连接;功放和发射换能器之间通过数据线连接;发射换能器和接收换能器1之间是无线水声信道;发射换能器和接收换能器2之间是无线水声信道;接收换能器1和滤波器1之间通过数据线连接;接收换能器2和滤波器2之间通过数据线连接;滤波器1和信号处理单元之间通过数据线连接;滤波器2和信号处理单元之间通过数据线连接;
信号处理单元1用于生成待发送的信号;功放用于放大待发送的信号;发射换能器用于将电信号转化为声信号发送出去;接收换能器1和接收换能器2用于接收声信号,并转化为电信号;滤波器1和滤波器2用于选择特定频段的信号,以提高信噪比;信号处理单元2用于对接收到的信号进行处理,计算出两路的时延差。
信号处理单元1由一块DSP实现,它包括发射信号样本存储器,发射信号控制模块和发射信号生成模块,发射信号样本存储器和发射信号生成模块连接,发射信号控制模块和发射信号生成模块连接。信号处理单元2由一块FPGA和一块DSP实现,FPGA与DSP间由数据线连接;FPGA包括参考信号样本存储器、相关器1、相关器2、包络计算模块1和包络计算模块2;参考信号样本存储器和相关器1、相关器2连接,相关器1、相关器2分别是与包络计算模块1和包络计算模块2连接;DSP包括峰值检测器1、峰值检测器2、时延差初测值计算模块、相位计算模块1、相位计算模块2、减法器1、时延差修正值计算模块和减法器2;包络计算模块1和包络计算模块分别与峰值检测器1、峰值检测器2连接,相位计算模块1、相位计算模块2分别与峰值检测器1、峰值检测器2连接,时延差初测值计算模块上连接有峰值检测器1、峰值检测器2和减法器2,减法器1与相位计算模块1、相位计算模块2和时延差修正值计算模连接。
本发明的优点是:本发明的方法和系统应用于对定位精度要求较高的短基线水声定位系统,为定位解算提供高精度的时延差测量值,提高测量时延差的精度。
附图说明:
图1为本发明系统的组成框图。
图2为信号处理单元1组成框图。
图3为信号处理单元2组成框图。
具体实施方式:
下面结合附图及实施例对本发明作进一步的说明。
如图所示,本发明所述的用于短基线水声定位系统中精确测量时延差的方法,其发射信号产生在信号处理单元1(1)中,通过发射信号样本存储器(9),发射信号控制模块(10)和发射信号生成模块(11)实现,信号时延差的测定在信号处理单元2(8)中通过参考信号样本存储器(12),相关器1(13),相关器,2(14),包络计算模块1(15),包络计算模块2(16),峰值检测器1(17),峰值检测器2(18),时延差初测值计算模块(19),相位计算模块1(20),相位计算模块2(21),减法器1(22),时延差修正值计算模块(23)和减法器2(24)实现,包括如下步骤:
步骤1、根据实际测量需要,选择线性调频脉冲信号的幅度A,频率下限fL,频率上限fH,扫频时间T。令:
获得线性调频脉冲信号的解析表达式如下:
根据信号发射的采样率FS1,对信号s(t)做离散化处理生成数字信号,表达式如下:
其中,TS1=1/FS1表示信号发射时的采样间隔。将s1(n)的数值存入发射信号样本存储器(9)。
根据信号接收的采样率FS2,对信号s(t)做离散化处理生成数字信号,表达式如下:
其中,TS2=1/FS2表示信号接收时的采样间隔。将s2(n)的数值存入参考信号样本存储器(12)。
步骤2、根据实际测量需要,设定发射信号控制模块(10)的控制逻辑。在发射信号控制模块(10)的控制下,发射信号生成模块(11)从发射信号样本存储器(9)中取出发射信号样本,经功放(2)放大后,由发射换能器(3)发射出去。接收换能器1(4)接收到的信号经滤波器1(6)处理后,以采样率FS2进行采样,得到信号x1(n)送入信号处理单元单元2(8);接收换能器2(5)接收到的信号经滤波器1(7)处理后,以采样率FS2进行采样,得到信号x2(n)送入信号处理单元单元2(8)。
步骤3、采集到的信号x1(n)和参考信号样本存储器(12)中存储的信号s2(n)在相关器1(13)进行相关运算,相关器1(13)的输出信号表达式如下:
采集到的信号x2(n)和参考信号样本存储器(12)中存储的信号s2(n)在相关器2(14)进行相关运算,相关器2(14)的输出信号表达式如下:
步骤4、R1(m)送入包络计算模块1(15),计算出相关峰包络B1(m);R2(m)送入包络计算模块2(16),计算出相关峰包络B2(m)。B1(m)送入峰值检测器1(17),计算出相关峰包络峰值位置M1;B2(m)送入峰值检测器2(18),计算出相关峰包络峰值位置M2。
峰值检测器1(17)的输出M1和峰值检测器2(18)的输出M2送入时延差初测值计算模块(19),时延差初测值计算模块(19)的输出为:
步骤5、R1(m)送入相位计算模块1(20),根据峰值检测器1(17)计算出的B1(m)的峰值位置M1,生成信号p1(m);R2(m)送入相位计算模块2(21),根据峰值检测器2(18)计算出的B2(m)的峰值位置M2,生成信号p2(m)。p1(m)和p2(m)的表达式如下:
p1(m)=R1(m-M1-L)0≤m≤2L-1
p2(m)=R2(m-M2-L)0≤m≤2L-1
其中,Im[·]表示取表达式的虚部,Re[·]表示取表达式的虚部。
步骤7、减法器1(22)的输出送入时延差修正值计算模块(23)计算时延差修正值τM。设存在W个非负整数满足
则时延差修正值计算模块(23)的输出为:
步骤8、时延差初测值计算模块(19)的输出τR和时延差修正值计算模块(23)的输出τM送入减法器2(24),减法器2(24)的输出
Δτ=τR-τM即为信号x1(n)和信号x2(n)的时延差。
上述测量并计算多普勒偏移的方法采用的系统为短基线水声定位系统的发射和接收系统。该系统包括:接收端和发射端。
发射端包括信号处理单元1(1),功放(2),发射换能器(3),接收端包括接收换能器1(4),接收换能器2(5),滤波器1(6),滤波器2(7),信号处理单元2(8)。信号处理单元1(1)与功放(2)之间通过数据线连接;功放(2)和发射换能器(3)之间通过数据线连接;发射换能器(3)和接收换能器1(4)之间是无线水声信道;发射换能器(3)和接收换能器2(5)之间是无线水声信道;接收换能器1(4)和滤波器1(6)之间通过数据线连接;接收换能器2(5)和滤波器2(7)之间通过数据线连接;滤波器1(6)和信号处理单元(8)之间通过数据线连接;滤波器2(7)和信号处理单元(8)之间通过数据线连接。
信号处理单元1(1)用于生成待发送的信号,可由一块DSP实现。它包括发射信号样本存储器(9),发射信号控制模块(10)和发射信号生成模块(11)。发射信号样本存储器(9)和发射信号生成模块(11)连接,发射信号控制模块(10)和发射信号生成模块(11)连接。
功放(2)用于放大待发送的信号。
发射换能器(3)用于将电信号转化为声信号发送出去。
接收换能器1(4)和接收换能器2(5)用于接收声信号,并转化为电信号。
滤波器1(6)和滤波器2(7)用于选择特定频段的信号,以提高信噪比。
信号处理单元2(8)用于对接收到的信号进行处理,计算出两路的时延差。
信号处理单元2由一块FPGA和一块DSP实现,FPGA与DSP间由数据线连接;FPGA包括参考信号样本存储器(12)、相关器1(13)、相关器2(14)、包络计算模块1(15)和包络计算模块2(16);参考信号样本存储器(12)和相关器1(13)、相关器2(14)连接,相关器1(13)、相关器2(14)分别是与包络计算模块1(15)和包络计算模块2(16)连接;DSP包括DSP实现峰值检测器1(17)、峰值检测器2(18)、时延差初测值计算模块(19)、相位计算模块1(20)、相位计算模块2(21)、减法器1(22)、时延差修正值计算模块(23)和减法器2(24);包络计算模块1(15)和包络计算模块2(16)分别与峰值检测器1(17)、峰值检测器2(18)连接,相位计算模块1(20)、相位计算模块2(21)分别与峰值检测器1(17)、峰值检测器2(18)连接,时延差初测值计算模块(19)上连接有峰值检测器1(17)、峰值检测器2(18)和减法器2(24),减法器1(22)与相位计算模块1(20)、相位计算模块2(21)和时延差修正值计算模块(23)连接。
Claims (1)
1.一种用于短基线水声定位系统中精确测量时延差的方法,其特征在于至少包括以下步骤:
1)短基线水声定位系统的发射端发射线性调频脉冲信号作为定位信号,形式如下:
其中,TS1是短基线水声定位系统发射信号的采样间隔,T是线性调频信号由最低频扫到最高频所经过的时间,k0是信号的调频斜率,fL是信号s1(n)的最低频率;
2)短基线水声定位系统的接收端有相关器1和相关器2,相关器1接收信号x1(n)并进行相关处理,相关器2接收信号x2(n)并进行相关处理;相关器1和相关器2进行相关处理时,使用的参考信号相同,参考信号形式如下:
其中,TS2是短基线水声定位系统接收信号的采样间隔;
设处理后,从相关器1输出的输出信号R1(m)和从相关器2输出的输出信号R2(m)的包络出现峰值对应的位置为M1和M2,那么时延差初测值为:
根据峰值位置M1和M2,由相位计算模块1生成的生成信号p1(m)和相位计算模块2生成的生成信号p2(m)分别为:
p1(m)=R1(m-M1-L) 0≤m≤2L-1
p2(m)=R2(m-M2-L) 0≤m≤2L-1
设存在W个非负整数满足:
则时延差修正值为:
这里的FS2是信号接收的采样率,fH是频率上限;
那么,时延差初测值和时延差修正值的差值,即为时延差值
Δτ=τR-τM
参数Δτ的物理意义是,在发射端发出的信号经过信道传播后,接收端两个不同的通道接收到信号的时刻相差了Δτ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010540329 CN102096064B (zh) | 2010-11-12 | 2010-11-12 | 一种用于短基线水声定位系统中精确测量时延差的方法及其系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010540329 CN102096064B (zh) | 2010-11-12 | 2010-11-12 | 一种用于短基线水声定位系统中精确测量时延差的方法及其系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102096064A CN102096064A (zh) | 2011-06-15 |
CN102096064B true CN102096064B (zh) | 2013-10-30 |
Family
ID=44129234
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010540329 Active CN102096064B (zh) | 2010-11-12 | 2010-11-12 | 一种用于短基线水声定位系统中精确测量时延差的方法及其系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102096064B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102901839B (zh) * | 2012-11-05 | 2014-08-13 | 杭州电子科技大学 | 一种低复杂度的声学多普勒测流系统速度估算方法 |
CN103200677B (zh) * | 2013-04-02 | 2015-05-13 | 东南大学 | 一种基于lte定位参考信号特征的精确时延计算方法 |
CN103777177B (zh) * | 2013-11-08 | 2018-02-02 | 中国船舶重工集团公司七五○试验场 | 一种基于宽带信号时延检测的超短基线水下目标定位方法 |
US9594062B2 (en) * | 2014-09-29 | 2017-03-14 | General Electric Company | System and method for testing transducer elements of an acoustic probe |
CN105721077A (zh) * | 2016-02-16 | 2016-06-29 | 武汉滨湖电子有限责任公司 | 一种射频通道间时延差测量装置及测量方法 |
CN108123902B (zh) * | 2016-11-29 | 2021-05-25 | 联发科技股份有限公司 | 用来估计一影音信号的信道状态的估计方法及相关的估计电路与接收器 |
CN106603117A (zh) * | 2016-12-09 | 2017-04-26 | 江苏理工学院 | 一种水下测量传播时延的方法 |
CN109782661B (zh) * | 2019-01-04 | 2020-10-16 | 中国科学院声学研究所东海研究站 | 基于fpga实现可重配置和多输出的实时处理系统及方法 |
CN109814094B (zh) * | 2019-02-01 | 2021-01-26 | 中国科学院声学研究所东海研究站 | 一种多目标水声定位时延估计算法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101021562A (zh) * | 2006-11-20 | 2007-08-22 | 重庆大学 | 一种低信噪比下超声波测距的方法及其测距仪 |
CN101692629A (zh) * | 2009-05-07 | 2010-04-07 | 嘉兴中科声学科技有限公司 | 一种测量并计算多普勒偏移的方法及其系统 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1110709C (zh) * | 2001-09-13 | 2003-06-04 | 中国科学院声学研究所 | 用于测量海底微地貌的高分辨测深侧扫声纳系统和测量方法 |
CN100492051C (zh) * | 2006-11-10 | 2009-05-27 | 哈尔滨工程大学 | 多功能水声测距仪 |
CN101872020B (zh) * | 2010-06-18 | 2012-10-31 | 华南理工大学 | 基于频谱变换的水下运动目标定位导航方法和装置 |
CN202110277U (zh) * | 2010-11-12 | 2012-01-11 | 嘉兴中科声学科技有限公司 | 一种用于短基线水声定位系统中精确测量时延差的装置 |
-
2010
- 2010-11-12 CN CN 201010540329 patent/CN102096064B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101021562A (zh) * | 2006-11-20 | 2007-08-22 | 重庆大学 | 一种低信噪比下超声波测距的方法及其测距仪 |
CN101692629A (zh) * | 2009-05-07 | 2010-04-07 | 嘉兴中科声学科技有限公司 | 一种测量并计算多普勒偏移的方法及其系统 |
Non-Patent Citations (4)
Title |
---|
张倩等.一种变步长的自适应特征值分解时延估计方法.《声学技术》.2009,第28卷(第2期),137-141. * |
惠娟等.时延/时延差组合定位模型抗距离模糊技术研究.《声学技术》.2007,第26卷(第3期),385-388. * |
李国栋等.一种水声脉冲测距系统测时方法.《声学技术》.2010,第29卷(第1期),28-31. * |
殷敬伟等.Pattern时延差编码四信道水声通信技术研究.《应用声学》.2006,第25卷(第3期),180-186. * |
Also Published As
Publication number | Publication date |
---|---|
CN102096064A (zh) | 2011-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102096064B (zh) | 一种用于短基线水声定位系统中精确测量时延差的方法及其系统 | |
CN104237871B (zh) | 一种基于相位补偿的时延差估计方法 | |
CN106772393B (zh) | 一种改进的基于飞行时间检测的超声波测距方法 | |
CN101975949B (zh) | 一种多用户水下超短基线定位设备与方法 | |
CN100395547C (zh) | 混凝土灌注桩质量检测系统 | |
CN103941259B (zh) | 一种具备高抗干扰性的超声波测距方法与测距装置 | |
CN102901949B (zh) | 一种二维空间分布式声音相对定位方法及装置 | |
CN101813673B (zh) | 检测微量二元混合气体浓度的声学信号处理装置及方法 | |
CN102749107A (zh) | 高精度时差式单脉冲超声波流量计系统及其测量流量方法 | |
CN104374532A (zh) | 航天器在轨泄漏定向方法 | |
CN102841343A (zh) | 一种基于工控机的回声测深仪校准系统及其校准方法 | |
CN103308915A (zh) | 一种声波测距装置与方法 | |
CN101692629B (zh) | 一种测量并计算多普勒偏移的方法及其系统 | |
CN103728464A (zh) | 一种用于声学多普勒流速剖面仪的组合脉冲测速方法 | |
CN111796259A (zh) | 一种基于lfm互相关模特性逼近的子样本时延估计方法 | |
CN101266294B (zh) | 一种抑制雷达回波信号包络移动的补偿方法 | |
CN101576618A (zh) | 基于小波变换的水声定位测量方法和测量系统 | |
CN203037849U (zh) | 一种超声波测厚仪 | |
CN101504459A (zh) | 一种测量水声阵阵元相幅一致性的新方法 | |
CN105116371A (zh) | 一种基于连续发射调频信号的目标定位方法与装置 | |
CN110440896B (zh) | 一种超声波测量系统及测量方法 | |
CN105467395B (zh) | 一种超远距离超声波测量仪 | |
CN102645267A (zh) | 检测岩石类构件弹性波速度的方法和系统 | |
CN202110277U (zh) | 一种用于短基线水声定位系统中精确测量时延差的装置 | |
CN102707288B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |