CN111273294A - 一种基于时间反转的水声脉冲回波测距方法 - Google Patents
一种基于时间反转的水声脉冲回波测距方法 Download PDFInfo
- Publication number
- CN111273294A CN111273294A CN202010103114.4A CN202010103114A CN111273294A CN 111273294 A CN111273294 A CN 111273294A CN 202010103114 A CN202010103114 A CN 202010103114A CN 111273294 A CN111273294 A CN 111273294A
- Authority
- CN
- China
- Prior art keywords
- signal
- time
- path
- channel
- emax
- 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
- 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/06—Systems determining the position data of a target
- G01S15/08—Systems for measuring distance only
- G01S15/10—Systems for measuring distance only using transmission of interrupted, pulse-modulated waves
- G01S15/102—Systems for measuring distance only using transmission of interrupted, pulse-modulated waves using transmission of pulses having some particular characteristics
-
- 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/06—Systems determining the position data of a target
- G01S15/08—Systems for measuring distance only
- G01S15/10—Systems for measuring distance only using transmission of interrupted, pulse-modulated waves
- G01S15/102—Systems for measuring distance only using transmission of interrupted, pulse-modulated waves using transmission of pulses having some particular characteristics
- G01S15/104—Systems for measuring distance only using transmission of interrupted, pulse-modulated waves using transmission of pulses having some particular characteristics wherein the transmitted pulses use a frequency- or phase-modulated carrier wave
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/523—Details of pulse systems
- G01S7/526—Receivers
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/523—Details of pulse systems
- G01S7/526—Receivers
- G01S7/527—Extracting wanted echo signals
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明提供了一种基于时间反转的水声脉冲回波测距方法,通过对多径信道各路径时延进行估计,获取并设计虚拟信道;将接收信号通过虚拟信道,消除由多径效应引起的接收信号时间扩展现象;通过自适应能量加窗算法截断时反输出信号的旁峰成分,并提取信号的聚焦峰;对加窗信号与发射信号进行相关处理寻找接收信号的准确到达时间,信号到达时间与信号发射时间之差为信号在水下传播时长,最终利用声传播公式对测距系统与目标之间的距离进行测量。本发明消除了多径效应引起的信号畸变,提高了提取目标回波到达时间的准确性,进而提高多径环境下测距精度。
Description
技术领域
本发明属于信息信号处理领域,涉及一种脉冲测距方法。
背景技术
脉冲测距法是利用发射脉冲与接收信号之间的时间差来进行测距的方法。由于该方法理论简单易于实现,因此广泛应用于主动声呐系统。在水声脉冲测距系统中,测时精度决定了测距精度,如何提高测时精度是脉冲测距方法的关键问题之一。
测距系统存在系统时延,采用多次标定测量系统时延的方法可以有效降低系统时延带来的测时误差;水声接收机具有一定带宽使得一些边频不能通过,接收信号频谱的完整性受到影响,接收机的输出波形会产生一定失真,在时域上表现为脉冲前后沿出现起伏不再陡峭,利用三次样条插值法重构接收信号,精确求出每个载波周期的峰值提取出信号包络,在对信号包络进行平滑后通过求解斜率最大值来估计时延从而提高测时精度;复杂水声环境下,由界面反射或散射体散射引起的多径效应将会产生强混响,使得测距系统不能准确判断出反射信号到达的真实时刻,进而影响系统的测时精度。这种由多径效应产生的混响有两类常采用的处理方式,一类是针对混响的统计特性对其进行预白化处理,另一类是根据目标反射信号与混响相关性的差异,采用子空间分解理论从回波信号中将反射信号和混响信号分离开来,这两类方法具有一定提高信混比的能力,但是这些算法都是对多径效应带来的混响进行抑制,算法复杂实时性差。此外,由于不能对所有的混响进行抑制,干扰伪峰依然存在不利于对到达时间的精确估计。时间反转是一种空时聚焦技术,操作简单且利用其聚焦特性,可以消除多径效应引起的信号畸变,适用于提高多径环境下测距精度的问题。
发明内容
为了克服现有技术的不足,本发明提供一种基于时间反转的水声脉冲回波测距方法,通过对多径信道各路径时延进行估计,获取并设计虚拟信道;将接收信号通过虚拟信道,消除由多径效应引起的接收信号时间扩展现象;通过自适应能量加窗算法截断时反输出信号的旁峰成分,并提取信号的聚焦峰;对加窗信号与发射信号进行相关处理寻找接收信号的准确到达时间,信号到达时间与信号发射时间之差为信号在水下传播时长;最终利用声传播公式对测距系统与目标之间的距离进行精确测量。
本发明解决其技术问题所采用的技术方案包括以下步骤:
第一步,发射探测脉冲信号s(t),发射时间记为τ1;
式中,ζmn=ζm-ζn为第m路径与第n条路径时延之差;
第四步,通过自适应能量加窗算法截断旁峰成分并提取信号的聚焦峰,加窗后的接收信号TEmax≤t≤TEmax+Ts-1,式中,Ts为探测脉冲信号s(t)的时宽,TEmax为能量函数最大值所在位置,arg为当|z(t)|2取最大值时TEmax的取值,Ans(-t)为截断的虚拟信道输出信号,为截断噪声;
本发明的有益效果是:通过对多径信道各路径时延进行估计,获取并设计虚拟信道;将接收信号通过虚拟信道,消除由多径效应引起的接收信号时间扩展现象;通过自适应能量加窗算法截断时反输出信号的旁峰成分,并提取信号的聚焦峰;对加窗信号与发射信号进行相关处理寻找接收信号的准确到达时间,信号到达时间与信号发射时间之差为信号在水下传播时长。最终利用声传播公式对测距系统与目标之间的距离进行测量。本发明利用时间反转技术的聚焦特性,消除了多径效应引起的信号畸变,从而提高了提取目标回波到达时间的准确性,进而提高多径环境下测距精度。本发明为水下目标探测定位跟踪等技术提供理论支撑,为我国实施远海防卫与近海防御战略打下坚实基础。
附图说明
图1是本发明的方法流程图。
具体实施方式
针对水下多径效应影响测时精度进而影响测距精度的问题,本发明提出基于时间反转的脉冲回波测距方法。
本发明的主步骤如下:
第一步:发射探测脉冲信号
测距系统探测脉冲信号记为s(t),发射时间记为τ1。
第二步:接收目标反射信号
测距系统接收回波信号,
式中,为射线模型下水声多径信道,s(t)为第一步发射的探测脉冲信号,t为时间变量,*为卷积运算,∑为求和符号,n为信道多径路径变量,N为信道多径数量,An为第n条路径对应的幅值,ξn为第n条路径对应的时延,d(t)为接收端加性噪声。
第三步:估计水声信道各路径时延
第四步:对接收信号进虚拟时间反转处理
式中,yAVTR(t)为输出信号,ζmn=ζm-ζn为第m路径与第n条路径时延之差。
第四步:对虚拟时间反转处理后的信号进行自适应能量加窗处理
通过自适应能量加窗算法截断旁峰成分并提取信号的聚焦峰,
式中,z(t)为加窗后的接收信号为窗函数,Ts为探测脉冲信号s(t)的时宽,TEmax为能量函数最大值所在位置,其获取方法为arg为当|z(t)|2取最大值时TEmax的取值,Ans(-t)为截断的虚拟信道输出信号,为截断噪声。
第五步:获取相关函数
把上加窗后输出结果与发射信号s(-t)(s(t)取反)进行相关处理,
式中,Rz(τ)为相关函数,τ为时延变量,E[]为求期望,R-s(τ)=E[s(-t)s(-t)]为s(-t)的自相关函数。
第六步:估计接收信号到达时刻
相关函数中峰值所在位置为回波信号到达时刻,
第七步:估计测距系统与目标之间的距离
水声信号从发射到接收,传播的距离是测距系统与目标之间距离的两倍,
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
第一步:发射探测脉冲信号
为了便于对信道估计,采用设计波形作为探测脉冲信号,如线性调频信号,双曲调频信号或伪随机调制信号等。测距系统探测脉冲信号选用线性调频信号记为s(t),发射频率为10kHz~12kHz,发射时宽为100ms。
第二步:接收目标反射信号
测距系统接收回波信号,
式中,yAS(t)为接收的目标回波信号,为射线模型下水声多径信道,s(t)为,t为时间变量,*为卷积运算,∑为求和符号,n为信道多径路径变量,N为信道多径数量,An为第n条路径对应的幅值,ξn为第n条路径对应的时延,d(t)为接收端加性噪声。为了收集完整的时延扩展信号,增大信号接收时长,这里接收时长为200ms
第三步:估计水声信道各路径时延
根据相关法对水声信道各路径时延进行估计,
第四步:对接收信号进虚拟时间反转处理
式中,yAVTR(t)为输出信号,ζmn=ζm-ζn为第m路径与第n条路径时延之差。
第四步:对虚拟时间反转处理后的信号进行自适应能量加窗处理
通过自适应能量加窗算法截断旁峰成分并提取信号的聚焦峰,其实现方法如下所述。
yAVTR(t)的数字信号为yAVTR(n),n=1,2,3,…,Nz,Nz代表时反接收序列的长度满足Tz=Nzfs,fs为采样率取为100kHzs,Tz为yAVTR(t)的时宽;s(t)对应的数字信号为s(n),n=1,2,3,…,Ns,Ns代表时反接收序列的长度满足Ts=Nsfs且Ns<Nz,Ts为探测脉冲信号s(t)的时宽;由于时反具有恢复原始信号的特点,因此窗函数w(t)的数字信号w(n)的长度与s(n)同为Ns,
式中,NEmax为初值时刻,满足TEmax=NEmaxfs,TEmax为窗函数初始时间,窗函数w(n)中NEmax的选取遵从能量最大原则,其获取方法参考自适应能量加窗算法。加窗后的信号记为,
zw(n)=yAVTR(n)w(n),n=1,2,3,…,Ns
初始化设置为NEmax=1;
以信号能量作为衡量准则,其公式如下
式中,E(NEmax)为初始时刻为NEmax时zw(n)的能量;
NEmax=NEmax+L,L代表步长,重复步骤2,当L=1时至NEmax=Nz-Ns+1停止循环,当L>1时至NEmax=[(Nz-Ns)/L]L停止循环,[·]代表取整函数;
将数字信号形式转换成连续时间形式
第五步:获取相关函数
把上加窗后输出结果与发射信号s(-t)(s(t)取反)进行相关处理,
式中,Rz(τ)为相关函数,τ为时延变量,E[]为求期望,R-s(τ)=E[s(-t)s(-t)]为s(-t)的自相关函数。
第六步:估计接收信号到达时刻
相关函数中峰值所在位置为回波信号到达时刻,
第七步:估计测距系统与目标之间的距离
水声信号从发射到接收,传播的距离是测距系统与目标之间距离的两倍,
Claims (1)
1.一种基于时间反转的水声脉冲回波测距方法,其特征在于包括下述步骤:
第一步,发射探测脉冲信号s(t),发射时间记为τ1;
式中,ζmn=ζm-ζn为第m路径与第n条路径时延之差;
第四步,通过自适应能量加窗算法截断旁峰成分并提取信号的聚焦峰,加窗后的接收信号TEmax≤t≤TEmax+Ts-1,式中,Ts为探测脉冲信号s(t)的时宽,TEmax为能量函数最大值所在位置,arg为当|z(t)|2取最大值时TEmax的取值,Ans(-t)为截断的虚拟信道输出信号,为截断噪声;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010103114.4A CN111273294A (zh) | 2020-02-19 | 2020-02-19 | 一种基于时间反转的水声脉冲回波测距方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010103114.4A CN111273294A (zh) | 2020-02-19 | 2020-02-19 | 一种基于时间反转的水声脉冲回波测距方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111273294A true CN111273294A (zh) | 2020-06-12 |
Family
ID=71002159
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010103114.4A Pending CN111273294A (zh) | 2020-02-19 | 2020-02-19 | 一种基于时间反转的水声脉冲回波测距方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111273294A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111800202A (zh) * | 2020-06-28 | 2020-10-20 | 西北工业大学 | 一种基于Labview平台的水声网络节点测距方法 |
CN113253280A (zh) * | 2021-04-23 | 2021-08-13 | 青岛农业大学 | 一种联合双曲调频与线性调频的测距测速方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101183150A (zh) * | 2007-11-30 | 2008-05-21 | 浙江大学 | 基于高频时反阵的时反发射聚焦反波束形成方法 |
US20090323473A1 (en) * | 2008-06-27 | 2009-12-31 | Yoshiaki Tsurugaya | Target searching device, target searching program, and target searching method |
CN103091677A (zh) * | 2012-07-08 | 2013-05-08 | 西北工业大学 | 一种基于时反的均匀线列阵波束形成方法 |
CN107607937A (zh) * | 2017-08-21 | 2018-01-19 | 西安电子科技大学 | 基于时间反演的雷达目标测距方法 |
CN107942322A (zh) * | 2017-11-03 | 2018-04-20 | 西安电子科技大学 | 多径利用雷达时间反演目标距离估计方法 |
US20180149476A1 (en) * | 2016-11-29 | 2018-05-31 | The Trustees Of The Stevens Institute Of Technology | Method and apparauts for quantum measurement via mode matched photon conversion |
CN108845309A (zh) * | 2018-08-28 | 2018-11-20 | 西北工业大学 | 一种基于主动虚拟时反方法的水下目标doa估计方法 |
-
2020
- 2020-02-19 CN CN202010103114.4A patent/CN111273294A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101183150A (zh) * | 2007-11-30 | 2008-05-21 | 浙江大学 | 基于高频时反阵的时反发射聚焦反波束形成方法 |
US20090323473A1 (en) * | 2008-06-27 | 2009-12-31 | Yoshiaki Tsurugaya | Target searching device, target searching program, and target searching method |
CN103091677A (zh) * | 2012-07-08 | 2013-05-08 | 西北工业大学 | 一种基于时反的均匀线列阵波束形成方法 |
US20180149476A1 (en) * | 2016-11-29 | 2018-05-31 | The Trustees Of The Stevens Institute Of Technology | Method and apparauts for quantum measurement via mode matched photon conversion |
CN107607937A (zh) * | 2017-08-21 | 2018-01-19 | 西安电子科技大学 | 基于时间反演的雷达目标测距方法 |
CN107942322A (zh) * | 2017-11-03 | 2018-04-20 | 西安电子科技大学 | 多径利用雷达时间反演目标距离估计方法 |
CN108845309A (zh) * | 2018-08-28 | 2018-11-20 | 西北工业大学 | 一种基于主动虚拟时反方法的水下目标doa估计方法 |
Non-Patent Citations (6)
Title |
---|
RUIQIN ZHAO ET AL.: "Medium Access Control Mechanism for Multi-hop Underwater Acoustic Networks Utilizing Time Reversal", 《OCEANS 2015 - MTS/IEEE WASHINGTON》, 11 February 2016 (2016-02-11), pages 1 * |
徐复等: "多途条件下聚焦波束近程定位", 《声学技术》, vol. 26, no. 06, 15 December 2007 (2007-12-15), pages 1101 * |
惠俊英等: "二元阵被动时间反转镜定位技术研究", 《哈尔滨工程大学学报》, vol. 28, no. 11, 30 November 2007 (2007-11-30), pages 1247 * |
王铭灏等: "水声信号频率对时间反转探测性能的影响", 《水下无人系统学报》, vol. 26, no. 05, 15 October 2018 (2018-10-15), pages 480 * |
荆海霞等: "基于主动时间反转的目标探测性能研究", 《电视技术》, vol. 40, no. 08, 17 August 2016 (2016-08-17), pages 103 * |
赵发勇等: "《数字信号处理》", 30 April 2018, 北京师范大学出版集团安徽大学出版社, pages: 172 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111800202A (zh) * | 2020-06-28 | 2020-10-20 | 西北工业大学 | 一种基于Labview平台的水声网络节点测距方法 |
CN113253280A (zh) * | 2021-04-23 | 2021-08-13 | 青岛农业大学 | 一种联合双曲调频与线性调频的测距测速方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Van Walree et al. | A discrete-time channel simulator driven by measured scattering functions | |
CN109814094B (zh) | 一种多目标水声定位时延估计算法 | |
Walree | Channel sounding for acoustic communications: techniques and shallow-water examples | |
CN107561486B (zh) | 一种基于主动时反的浅海目标波达方向估计方法 | |
CN111273294A (zh) | 一种基于时间反转的水声脉冲回波测距方法 | |
CN109765545A (zh) | 一种基于相关峰包络检测的水声定位时延估计算法 | |
CN103618575B (zh) | 一种调频水声通信系统的自动实时帧同步方法 | |
Su et al. | Fast estimation of underwater acoustic multipath channel based on LFM signal | |
Bates et al. | Signal-to-reverberation ratio comparison of linear frequency modulated continuous active sonar and pulsed active sonar | |
CN113704685B (zh) | 一种基于垂直线列阵的深海盲解卷积方法 | |
EP1631840A2 (en) | Estimation of background noise and its effect on sonar range estimation | |
WO2005006009A2 (en) | Noise adaptive sonar signal processor | |
CN110673118A (zh) | 一种主动声纳单频脉冲串波形设计及检测算法 | |
Ling et al. | Probing waveforms and adaptive receivers for active sonar | |
CN103197282B (zh) | 基于幅度补偿的mvdr时反聚焦定位方法 | |
CN115685169A (zh) | 基于宽带keystone变换的水声弱运动目标检测方法 | |
CN111342949B (zh) | 一种水声移动通信的同步检测方法 | |
Wang et al. | Measurement-based analysis of characteristics of fast moving underwater acoustic communication channel | |
Das et al. | Analysis of effective signal design for active sensing of undersea objects/bottoms in tropical shallow waters | |
Chen et al. | On probing waveforms and adaptive receivers for active sonar | |
Malarkodi et al. | Characterization of underwater acoustic communication channel | |
CN111273303B (zh) | 双基地声纳强直达波干扰下目标微弱散射信号检测的方法 | |
Abtahi et al. | Energy allocation for parameter estimation in block cs-based distributed mimo systems | |
CN117970298B (zh) | 一种基于对称hfm信号的目标检测方法及系统 | |
Pelekanakis et al. | ADAPTIVE EQUALISATION FOR CONTINUOUS ACTIVE SONAR? |
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: 20200612 |