CN111342949A - 一种水声移动通信的同步检测方法 - Google Patents

一种水声移动通信的同步检测方法 Download PDF

Info

Publication number
CN111342949A
CN111342949A CN202010103115.9A CN202010103115A CN111342949A CN 111342949 A CN111342949 A CN 111342949A CN 202010103115 A CN202010103115 A CN 202010103115A CN 111342949 A CN111342949 A CN 111342949A
Authority
CN
China
Prior art keywords
signal
detection
receiver
filter
scale
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.)
Granted
Application number
CN202010103115.9A
Other languages
English (en)
Other versions
CN111342949B (zh
Inventor
申晓红
张之琛
王领
花飞
王海燕
马石磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202010103115.9A priority Critical patent/CN111342949B/zh
Publication of CN111342949A publication Critical patent/CN111342949A/zh
Application granted granted Critical
Publication of CN111342949B publication Critical patent/CN111342949B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L7/00Arrangements for synchronising receiver with transmitter
    • H04L7/0079Receiver details
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B11/00Transmission systems employing sonic, ultrasonic or infrasonic waves
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B13/00Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
    • H04B13/02Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/0224Channel estimation using sounding signals
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明提供了一种水声移动通信的同步检测方法,采用尺度展缩组合滤波器的接收机,该接收机的组合滤波器由不同时间展缩及时移的信号之和组成,提高了滤波器与不同运动速度的宽带多普勒回波信号的匹配度,从而提高宽带多普勒回波时的检测概率及同步时间估计精度。此外,本发明提出的接收机与匹配滤波器的运算相同,适合于实时、低功耗要求的移动移动水声通信的同步检测。

Description

一种水声移动通信的同步检测方法
技术领域
本发明属于水声通信与组网领域,涉及匹配滤波、水下信号处理、水声通信和水声组网等理论。
背景技术
水声通信大多为突发协作通信,在这种突发协作通信中,接收机能正确地检测并估计通信信号的到达时刻是突发协作通信首先要解决的重要问题。为了使接收机快速解决该问题,通常在突发协作通信的数据帧前增加一段收发双方约定好的同步信号,以实现突发信号的检测和时间同步。
为了在一定的虚警概率下获得高的检测概率及信号精准的到达时间,水声通信工作者进行了深入的研究,做了大量的工作,通常选取自相关性好的线性调频信号(LFM)、双曲调频信号(HFM)或伪随机调制信号作为同步信号,接收机通常采用匹配滤波器。这种接收机在无径向运动时获得了很好的检测性能及信号到达时刻的估计精度,但由于水声通信收发端的相对运动,运动的海面、海底反射以及声速不均匀引起的折射,使得水声信道是连续的时变扩展信道,此时产生的径向速度使得接收机的性能急剧下降,此时即使能够正确检测到,但是由于信号到达时刻的估计误差较大,使得接收信号的解调误码率增大。目前各类的研究都是假设没有多普勒频偏,这就使得根据匹配滤波理论设计的滤波器出现失配问题。因此考虑多普勒频偏,根据尺度展缩构造组合滤波器,可以提高滤波器与不同运动速度的宽带多普勒回波信号的匹配度,进而提高接收机的检测性能和同步时间估计精度。
发明内容
为了克服现有技术的不足,本发明提供一种水声移动通信的同步检测方法,采用尺度展缩组合滤波器的接收机,该接收机的组合滤波器由不同时间展缩及时移的信号之和组成,提高了滤波器与不同运动速度的宽带多普勒回波信号的匹配度,从而提高宽带多普勒回波时的检测概率及同步时间估计精度。此外,本发明提出的接收机与匹配滤波器的运算相同,适合于实时、低功耗要求的移动移动水声通信的同步检测。
本发明解决其技术问题所采用的技术方案包括以下步骤:
第一步,发射探测信号s(t),探测信号为LFM信号、HFM信号或伪随机调制信号;
第二步,构造尺度组合滤波器,尺度组合滤波器为发射信号经过不同的尺度展缩及时间延迟之和,其冲激响应
Figure BDA0002387539540000021
式中,vth为速度间隔,满足vth=c(apth-1),apth为尺度间隔,c为水中声速,
Figure BDA0002387539540000022
为第p条路径的时间延迟,f0为线性调频信号起始频率,M为线性调频信号频率调制斜率,vmax为系统的最大相对运动速度,组合滤波器的个数
Figure BDA0002387539540000023
Figure BDA0002387539540000024
为向上取整运算符,
Figure BDA0002387539540000025
为能量归一化系数,∑为求和运算符号,*为卷积运算符号;
第三步,设计基于尺度组合滤波器的接收机,接收机的输出最大值为检验统计量Ts=ApT,Ap为第p条路径对应的幅值,T为发射信号持续时间;根据检测任务给定的虚警概率PFA,由奈曼-皮尔逊准则得接收机的检测门限
Figure BDA0002387539540000026
检测概率
Figure BDA0002387539540000027
式中,n~N(0,σ2)为加性高斯白噪声,Q-1(·)为Q(·)函数的反函数,
Figure BDA0002387539540000028
第四步,回波信号与尺度组合滤波器中的第i阶组件相匹配,尺度组合滤波器输出的最大值
Figure BDA0002387539540000029
取尺度组合滤波器的最大值及第i-1阶和i+1阶组件的输出的极大值之和为检验统计量;
第五步,设置判决门限对第四步中的检验统计量进行检测判决。
本发明的有益效果是:提出了一种尺度展缩组合滤波器的接收机,该接收机的组合滤波器由不同时间展缩及时移的信号之和组成,提高了滤波器与不同运动速度的宽带多普勒回波信号的匹配度,从而提高宽带多普勒回波时的检测概率及同步时间估计精度。此外,本专利提出的接收机与匹配滤波器的运算相同,适合于实时、低功耗要求的移动移动水声通信的同步检测。本发明具有较强的适用性,不仅可用于对LFM信号的同步检测,同样适用于对HFM或伪随机信号的同步检测。为水下组网通信及水下网络安全提供技术支撑,为建设智慧海洋、发展海洋经济、建立海上强国打下坚实基础。
附图说明
图1是本发明设计的接收机滤波检测原理框图;
图2是本发明设计的组合滤波器结构图。
具体实施方式
本发明解决其技术问题所采用的技术方案包括以下步骤:
第一步:发射探测信号
发射信号记为s(t),可以为LFM信号、HFM信号或伪随机调制信号。
第二步:构造尺度组合滤波器
尺度组合滤波器为发射信号经过不同的尺度展缩及时间延迟之和,其冲激响应如下式所示,
Figure BDA0002387539540000031
式中,vth为速度间隔,满足vth=c(apth-1),apth为尺度间隔,c为水中声速,
Figure BDA0002387539540000032
为第p条路径的时间延迟,f0为线性调频信号起始频率,M为线性调频信号频率调制斜率,vmax为系统的最大相对运动速度为,组合滤波器的个数
Figure BDA0002387539540000033
Figure BDA0002387539540000034
为向上取整运算符,
Figure BDA0002387539540000035
为能量归一化系数,∑为求和运算符号,*为卷积运算符号,s(t)为发射信号。
第三步:设计基于尺度组合接收滤波器的接收机
接收机设计的关键是在恒虚警条件下选取合适的检验统计量获取满足检测任务的检测概率。取组合滤波器中输出最大值为检验统计量Ts=ApT,Ap为第p条路径对应的幅值,T为发射信号持续时间,n~N(0,σ2)为加性高斯白噪声,N(0,σ2)代表均值为0,方差为σ2的正态分布,归一化检验统计量如下,
Figure BDA0002387539540000036
根据检测任务给定的虚警概率PFA,由奈曼-皮尔逊准则得检测门限
Figure BDA0002387539540000041
检测概率
Figure BDA0002387539540000042
式中,Q-1(·)为Q(·)函数的反函数,
Figure BDA0002387539540000043
第四步:计算检验统计量
回波信号与组合滤波器中的第i阶组件相匹配,组合滤波器输出的最大值为
Figure BDA0002387539540000044
为了提高检测概率,取尺度组合滤波器的最大值及第i-1阶和i+1阶组件的输出的极大值之和为检验统计量。
第五步:设置判决门限对检验统计量进行检测判决
设置判决门限对第四步中的检验统计量进行检测判决,本发明中认为匹配滤波器输出最大值变为原始值(这里的原始值为接收信号经过传统匹配滤波器后输出的最大值)的一半时无法正确检测,因此选择此时的最大值作为判决门限。
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
第一步:发射探测信号
发射信号为持续时间T的矩形包络LFM信号:
Figure BDA0002387539540000045
式中,f0为载波频率,一般取9-16kHz,M为线性调频指数,一般取500Hz~2000Hz,瞬时频率为fu(t)=f0+Mt,频率调制宽度W=|M|T,
Figure BDA0002387539540000046
为Woodward函数。
此外,发射还可以为LFM信号、HFM信号或伪随机调制信号。
水声信道是连续的时变扩展信道其时变冲击响应可表示为
Figure BDA0002387539540000047
其中,Np为信道在采样点上可分离出的独立路径数,Ap、τp、ap分别为接收到的每一条路径的幅值、时间延迟、多普勒尺度,且
Figure BDA0002387539540000051
Figure BDA0002387539540000059
为径向速度,c为水中声速。
到达接收机的信号r(t)可表示为
Figure BDA0002387539540000052
式中,τp为信号沿第p条路径传播到达接收机的时刻。
第二步:构造尺度组合滤波器
第p条路径为幅值最大的路径,由多普勒频移和时延引起的接收信号为
Figure BDA0002387539540000053
将匹配滤波器的冲激响应函数设计为与发射信号相匹配,即
Figure BDA0002387539540000054
接收信号rp(t)经过hr(t)表示的滤波器后,其输出信号为
Figure BDA0002387539540000055
输出信号最大值与相对运动速度的关系如下
Figure BDA0002387539540000056
式中,
Figure BDA0002387539540000057
Figure BDA0002387539540000058
为菲涅尔积分。
此时
t=f0(1-ap)/M+apτp
为了提高信号到达时刻的估计精度,接收滤波器中增加的尺度展缩信号必须为时移信号。尺度组合滤波器为发射信号经过不同的尺度展缩及时间延迟之和,其冲激响应如下式所示,
Figure BDA0002387539540000061
式中,apth为尺度间隔,vth为速度间隔,满足apth=1+vth/c,c为水中声速,
Figure BDA0002387539540000062
为时间延迟,vmax为系统的最大相对运动速度为,组合滤波器的个数
Figure BDA0002387539540000063
Figure BDA0002387539540000064
为向上取整运算符,
Figure BDA0002387539540000065
为能量归一化系数,Σ为求和运算符号,*为卷积运算符号,s(t)为发射信号。
第三步:设计基于尺度组合接收滤波器的接收机
取匹配滤波器输出最大值为检验统计量Ts=ApT,Ap为第p条路径对应的幅值,T为发射信号持续时间,n~N(0,σ2)为加性高斯白噪声,N(0,σ2)代表均值为0,方差为σ2的正态分布,归一化检验统计量如下,
Figure BDA0002387539540000066
根据检测任务给定的虚警概率PFA,由奈曼-皮尔逊准则得检测门限
Figure BDA0002387539540000067
检测概率
Figure BDA0002387539540000068
式中,Q-1(·)为Q(·)函数的反函数,
Figure BDA0002387539540000069
第四步:计算检验统计量
回波信号与组合滤波器中的第i阶组件相匹配,组合滤波器输出的最大值为
Figure BDA00023875395400000610
为了提高检测概率,取尺度组合滤波器的最大值及第i-1阶和i+1阶组件的输出的极大值之和为检验统计量。
第五步:设置判决门限对检验统计量进行检测判决
设置判决门限对第四步中的检验统计量进行检测判决,本发明中认为匹配滤波器输出最大值变为原始值(这里的原始值为接收信号经过传统匹配滤波器后输出的最大值)的一半时无法正确检测,因此选择此时的最大值作为判决门限。

Claims (1)

1.一种水声移动通信的同步检测方法,其特征在于包括下述步骤:
第一步,发射探测信号s(t),探测信号为LFM信号、HFM信号或伪随机调制信号;
第二步,构造尺度组合滤波器,尺度组合滤波器为发射信号经过不同的尺度展缩及时间延迟之和,其冲激响应
Figure FDA0002387539530000011
式中,vth为速度间隔,满足vth=c(apth-1),apth为尺度间隔,c为水中声速,
Figure FDA0002387539530000012
为第p条路径的时间延迟,f0为线性调频信号起始频率,M为线性调频信号频率调制斜率,vmax为系统的最大相对运动速度,组合滤波器的个数
Figure FDA0002387539530000013
Figure FDA0002387539530000014
为向上取整运算符,
Figure FDA0002387539530000015
为能量归一化系数,∑为求和运算符号,*为卷积运算符号;
第三步,设计基于尺度组合滤波器的接收机,接收机的输出最大值为检验统计量Ts=ApT,Ap为第p条路径对应的幅值,T为发射信号持续时间;根据检测任务给定的虚警概率PFA,由奈曼-皮尔逊准则得接收机的检测门限
Figure FDA0002387539530000016
检测概率
Figure FDA0002387539530000017
式中,n~N(0,σ2)为加性高斯白噪声,Q-1(·)为Q(·)函数的反函数,
Figure FDA0002387539530000018
第四步,回波信号与尺度组合滤波器中的第i阶组件相匹配,尺度组合滤波器输出的最大值
Figure FDA0002387539530000019
取尺度组合滤波器的最大值及第i-1阶和i+1阶组件的输出的极大值之和为检验统计量;
第五步,设置判决门限对第四步中的检验统计量进行检测判决。
CN202010103115.9A 2020-02-19 2020-02-19 一种水声移动通信的同步检测方法 Active CN111342949B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010103115.9A CN111342949B (zh) 2020-02-19 2020-02-19 一种水声移动通信的同步检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010103115.9A CN111342949B (zh) 2020-02-19 2020-02-19 一种水声移动通信的同步检测方法

Publications (2)

Publication Number Publication Date
CN111342949A true CN111342949A (zh) 2020-06-26
CN111342949B CN111342949B (zh) 2021-06-11

Family

ID=71185341

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010103115.9A Active CN111342949B (zh) 2020-02-19 2020-02-19 一种水声移动通信的同步检测方法

Country Status (1)

Country Link
CN (1) CN111342949B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112152959A (zh) * 2020-11-25 2020-12-29 鹏城实验室 水声通信探测信号处理方法、装置、设备及存储介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140126334A1 (en) * 2009-12-16 2014-05-08 Shb Instruments, Inc. Underwater acoustic navigation systems and methods
CN105323198A (zh) * 2014-06-13 2016-02-10 中国科学院声学研究所 一种利用双曲调频进行水下信号发射和接收的方法
US20180054787A1 (en) * 2016-08-18 2018-02-22 Allen-Vanguard Corporation System and method for providing adaptive synchronization of lte communication systems
CN108243138A (zh) * 2018-01-11 2018-07-03 福建星海通信科技有限公司 一种适用于水声通信系统的组合多普勒估计方法
CN108761384A (zh) * 2018-04-28 2018-11-06 西北工业大学 一种抗差的传感器网络目标定位方法
CN109391573A (zh) * 2017-08-02 2019-02-26 中国科学院声学研究所 一种基于线性调频信号的时频二维同步方法
CN109474295A (zh) * 2018-12-24 2019-03-15 南京屹信航天科技有限公司 一种用于星载测控设备的基带处理方法
CN109474305A (zh) * 2018-12-05 2019-03-15 唐山照澜海洋科技有限公司 一种5g多载波扩频水声通信方法
CN110247867A (zh) * 2019-05-16 2019-09-17 中国科学院声学研究所 水声多普勒估计方法及装置、水声通信方法及系统
CN110535537A (zh) * 2019-09-25 2019-12-03 哈尔滨工程大学 一种水下通信探测一体化方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140126334A1 (en) * 2009-12-16 2014-05-08 Shb Instruments, Inc. Underwater acoustic navigation systems and methods
CN105323198A (zh) * 2014-06-13 2016-02-10 中国科学院声学研究所 一种利用双曲调频进行水下信号发射和接收的方法
US20180054787A1 (en) * 2016-08-18 2018-02-22 Allen-Vanguard Corporation System and method for providing adaptive synchronization of lte communication systems
CN109391573A (zh) * 2017-08-02 2019-02-26 中国科学院声学研究所 一种基于线性调频信号的时频二维同步方法
CN108243138A (zh) * 2018-01-11 2018-07-03 福建星海通信科技有限公司 一种适用于水声通信系统的组合多普勒估计方法
CN108761384A (zh) * 2018-04-28 2018-11-06 西北工业大学 一种抗差的传感器网络目标定位方法
CN109474305A (zh) * 2018-12-05 2019-03-15 唐山照澜海洋科技有限公司 一种5g多载波扩频水声通信方法
CN109474295A (zh) * 2018-12-24 2019-03-15 南京屹信航天科技有限公司 一种用于星载测控设备的基带处理方法
CN110247867A (zh) * 2019-05-16 2019-09-17 中国科学院声学研究所 水声多普勒估计方法及装置、水声通信方法及系统
CN110535537A (zh) * 2019-09-25 2019-12-03 哈尔滨工程大学 一种水下通信探测一体化方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112152959A (zh) * 2020-11-25 2020-12-29 鹏城实验室 水声通信探测信号处理方法、装置、设备及存储介质
CN112152959B (zh) * 2020-11-25 2021-02-23 鹏城实验室 水声通信探测信号处理方法、装置、设备及存储介质
WO2022110409A1 (zh) * 2020-11-25 2022-06-02 鹏城实验室 水声通信探测信号处理方法、装置、设备及存储介质

Also Published As

Publication number Publication date
CN111342949B (zh) 2021-06-11

Similar Documents

Publication Publication Date Title
CN101594185B (zh) 移动水声通信信号的多普勒估计与同步方法
CN109814094B (zh) 一种多目标水声定位时延估计算法
CN106019263B (zh) 基于多亮点回波模型的目标径向速度测量方法
CN103197318B (zh) 一种基于Pattern时延编码水下声定位的时延估计方法
CN111580110B (zh) 一种基于浅海多径时延的复合码水声测距方法
CN103176166B (zh) 一种用于水声被动定位的信号到达时延差跟踪算法
CN106879068B (zh) 一种强多径环境下信号的到达时间估计方法
CN109765545A (zh) 一种基于相关峰包络检测的水声定位时延估计算法
CN110247867B (zh) 水声多普勒估计方法及装置、水声通信方法及系统
CN103618575B (zh) 一种调频水声通信系统的自动实时帧同步方法
Mallat et al. UWB based positioning in multipath channels: CRBs for AOA and for hybrid TOA-AOA based methods
CN111342949B (zh) 一种水声移动通信的同步检测方法
KR20170029899A (ko) 능동 소나 시스템의 lfm신호 파라미터 추정 장치 및 방법
CN111263321A (zh) 一种提高toa室内超宽带定位精度的方法
Cui et al. Timing estimation of multiple hyperbolic frequency‐modulated signals based on multicarrier underwater acoustic communication
Liu et al. Sch: a speed measurement method of combined hyperbolic frequency modulation signals
WO2004111674A2 (en) Estimation of background noise and its effect on sonar range estimation
Yi-Jia et al. Modeling and estimation of the space-time varying channels
CN114019453B (zh) 一种基于水声基线定位系统的测距方法
Diamant Clustering approach for detection and time of arrival estimation of hydrocoustic signals
CN106646436B (zh) 一种基于信号宽窄带模糊度的侦察信号参数估计方法
CN109768812B (zh) 一种基于混沌调频的水声通信多普勒估计与同步方法
CN112466330A (zh) 一种多途信道下噪声源的声源级估计方法
McGaughey et al. Improved torpedo range estimation using the fast orthogonal search
Zhang et al. Experimental studies of underwater acoustic communications over multipath channels

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
GR01 Patent grant
GR01 Patent grant