CN111342949B - Synchronous detection method for underwater acoustic mobile communication - Google Patents

Synchronous detection method for underwater acoustic mobile communication Download PDF

Info

Publication number
CN111342949B
CN111342949B CN202010103115.9A CN202010103115A CN111342949B CN 111342949 B CN111342949 B CN 111342949B CN 202010103115 A CN202010103115 A CN 202010103115A CN 111342949 B CN111342949 B CN 111342949B
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.)
Active
Application number
CN202010103115.9A
Other languages
Chinese (zh)
Other versions
CN111342949A (en
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/en
Publication of CN111342949A publication Critical patent/CN111342949A/en
Application granted granted Critical
Publication of CN111342949B publication Critical patent/CN111342949B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L7/00—Arrangements for synchronising receiver with transmitter
    • H04L7/0079—Receiver details
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B11/00—Transmission systems employing sonic, ultrasonic or infrasonic waves
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B13/00—Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
    • H04B13/02—Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202—Channel estimation
    • H04L25/0224—Channel estimation using sounding signals
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE 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/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Abstract

The invention provides a synchronous detection method of underwater acoustic mobile communication, which adopts a receiver of a scale expansion and contraction combined filter, wherein the combined filter of the receiver consists of the sum of signals with different time expansions and time shifts, and the matching degree of the filter and broadband Doppler echo signals with different movement speeds is improved, so that the detection probability and the synchronous time estimation precision during the broadband Doppler echo are improved. In addition, the receiver provided by the invention has the same operation with the matched filter, and is suitable for synchronous detection of mobile underwater acoustic communication with requirements on real time and low power consumption.

Description

Synchronous detection method for underwater acoustic mobile communication
Technical Field
The invention belongs to the field of underwater acoustic communication and networking, and relates to theories of matched filtering, underwater signal processing, underwater acoustic communication, underwater acoustic networking and the like.
Background
In the burst-oriented cooperative communication, the receiver can accurately detect and estimate the arrival time of the communication signal, which is an important problem to be solved first. In order to enable the receiver to quickly solve the problem, a section of synchronization signal agreed by the transmitting and receiving parties is usually added before a data frame of burst cooperative communication to achieve detection and time synchronization of the burst signal.
In order to obtain high detection probability and accurate arrival time of signals under a certain false alarm probability, underwater acoustic communication workers carry out intensive research and do a lot of work, and usually select a linear frequency modulation signal (LFM), a hyperbolic frequency modulation signal (HFM) or a pseudorandom modulation signal with good autocorrelation as a synchronization signal, and a receiver usually adopts a matched filter. The receiver obtains good detection performance and estimation accuracy of signal arrival time when no radial movement exists, but due to relative movement of underwater acoustic communication transceiving ends, moving sea surface and seabed reflection and refraction caused by uneven sound velocity, an underwater acoustic channel is a continuous time-varying extension channel, the performance of the receiver is sharply reduced by the generated radial velocity, and even if the receiver can correctly detect the underwater acoustic channel, the demodulation error rate of a received signal is increased due to large estimation error of the signal arrival time. At present, various researches assume that no Doppler frequency offset exists, so that a filter designed according to a matched filtering theory has a mismatch problem. Therefore, the Doppler frequency offset is considered, the combined filter is constructed according to scale expansion and contraction, the matching degree of the filter and the broadband Doppler echo signals with different motion speeds can be improved, and the detection performance and the synchronous time estimation precision of the receiver are further improved.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides a synchronous detection method of underwater acoustic mobile communication, which adopts a receiver of a scale-expansion combined filter, wherein the combined filter of the receiver consists of the sum of signals with different time expansions and time shifts, and the matching degree of the filter and broadband Doppler echo signals with different movement speeds is improved, so that the detection probability and the synchronous time estimation precision during the broadband Doppler echo are improved. In addition, the receiver provided by the invention has the same operation with the matched filter, and is suitable for synchronous detection of mobile underwater acoustic communication with requirements on real time and low power consumption.
The technical scheme adopted by the invention for solving the technical problem comprises the following steps:
firstly, emitting a detection signal s (t), wherein the detection signal is an LFM signal, an HFM signal or a pseudo-random modulation signal;
second, construct the scale combination filter, which is the sum of the time delay and the different scale expansion of the transmitting signal, and the impulse response
Figure BDA0002387539540000021
In the formula, vthFor the interval of velocity, satisfy vth＝c(apth-1)，apthIs the scale interval, c is the speed of sound in water,
Figure BDA0002387539540000022
time delay for the p-th path, f0Is the chirp start frequency, M is the chirp frequency modulation slope, vmaxThe number of combined filters for maximum relative speed of motion of the system
Figure BDA0002387539540000023
Figure BDA0002387539540000024
In order to round up the operator,
Figure BDA0002387539540000025
the energy normalization coefficient, sigma, and convolution operation symbol;
thirdly, designing a receiver based on a scale combination filter, wherein the maximum output value of the receiver is a test statistic Ts＝ApT，ApThe amplitude corresponding to the p path is T, and the duration of the transmitted signal is T; false alarm probability P given according to detection taskFADetection threshold of receiver by Nelman-Pearson criterion
Figure BDA0002387539540000026
Probability of detection
Figure BDA0002387539540000027
Wherein N to N (0, σ)2) Is additive white Gaussian noise, Q-1(. h) is the inverse of the Q (. h) function,
Figure BDA0002387539540000028
fourthly, the echo signal is matched with the ith order component in the scale combination filter, and the maximum value output by the scale combination filter
Figure BDA0002387539540000029
Taking the maximum of the scaled combined filter and the maximum of the output of the (i-1) th and (i + 1) th order componentsThe sum is the test statistic;
and fifthly, setting a judgment threshold to carry out detection judgment on the test statistic in the fourth step.
The invention has the beneficial effects that: the combined filter of the receiver consists of the sum of signals with different time expansion and time shift, and improves the matching degree of the filter and broadband Doppler echo signals with different movement speeds, thereby improving the detection probability and the synchronous time estimation precision during the broadband Doppler echo. In addition, the receiver provided by the patent has the same operation with a matched filter, and is suitable for synchronous detection of mobile underwater acoustic communication with real-time and low-power consumption requirements. The invention has strong applicability, not only can be used for synchronous detection of LFM signals, but also is suitable for synchronous detection of HFM or pseudo-random signals. The method provides technical support for underwater networking communication and underwater network safety, and lays a solid foundation for building smart oceans, developing oceanic economy and building oceanic strong countries.
Drawings
FIG. 1 is a block diagram of the receiver filter detection scheme of the present invention;
fig. 2 is a block diagram of a combined filter designed by the present invention.
Detailed Description
The technical scheme adopted by the invention for solving the technical problem comprises the following steps:
the first step is as follows: transmitting a probing signal
The transmitted signal is denoted as s (t), and may be an LFM signal, an HFM signal, or a pseudo-random modulation signal.
The second step is that: construction scale combined filter
The impulse response of the combined scale filter is the sum of the time delay and the different scale expansion of the transmission signal,
Figure BDA0002387539540000031
in the formula, vthFor the interval of velocity, satisfy vth＝c(apth-1)，apthIs the scale interval, c is the speed of sound in water,
Figure BDA0002387539540000032
time delay for the p-th path, f0Is the chirp start frequency, M is the chirp frequency modulation slope, vmaxThe maximum relative movement speed of the system is the number of the combined filters
Figure BDA0002387539540000033
Figure BDA0002387539540000034
In order to round up the operator,
Figure BDA0002387539540000035
the energy normalization coefficient, Σ, the symbol of the summation operation, the symbol of the convolution operation, and s (t) are the transmission signals.
The third step: receiver for designing receiving filter based on scale combination
The key of the receiver design is to select proper test statistics under the condition of constant false alarm to obtain the detection probability meeting the detection task. Taking the maximum output value in the combined filter as the test statistic Ts＝ApT，ApIs the amplitude corresponding to the p-th path, T is the duration of the transmitted signal, N-N (0, sigma)2) Is additive white Gaussian noise, N (0, σ)2) Represents a mean of 0 and a variance of σ2The normalized test statistic is as follows,
Figure BDA0002387539540000036
false alarm probability P given according to detection taskFAThe detection threshold is obtained by the Nelman-Pearson criterion
Figure BDA0002387539540000041
Probability of detection
Figure BDA0002387539540000042
In the formula, Q-1(. h) is the inverse of the Q (. h) function,
Figure BDA0002387539540000043
the fourth step: computing test statistics
The echo signal is matched to the ith order component in a combined filter, the combined filter output has a maximum value of
Figure BDA0002387539540000044
To increase the detection probability, the sum of the maximum of the scale combining filter and the maximum of the output of the i-1 th and i +1 th order components is taken as the test statistic.
The fifth step: setting a judgment threshold to carry out detection judgment on the test statistic
And setting a decision threshold to carry out detection decision on the test statistic in the fourth step, wherein in the invention, the maximum value output by the matched filter is considered to be half of the original value (the original value is the maximum value output by the received signal after passing through the traditional matched filter), so that the maximum value at the moment can not be correctly detected.
The present invention will be further described with reference to the following drawings and examples, which include, but are not limited to, the following examples.
The first step is as follows: transmitting a probing signal
The transmission signal is a rectangular envelope LFM signal of duration T:
Figure BDA0002387539540000045
in the formula (f)0The carrier frequency is 9-16kHz, M is linear frequency modulation index, 500-2000 Hz, and the instantaneous frequency is fu(t)＝f0+MT, the frequency modulation width W ═ M | T,
Figure BDA0002387539540000046
is Woodward function.
The transmission may also be an LFM signal, an HFM signal, or a pseudo-random modulated signal.
The time-varying impulse response of an underwater acoustic channel being a continuous time-varying extension channel can be expressed as
Figure BDA0002387539540000047
Wherein N ispFor the number of independent paths separable from the channel at the sampling point, Ap、τp、apRespectively the amplitude, time delay, Doppler scale of each received path, and
Figure BDA0002387539540000051
Figure BDA0002387539540000059
is the radial velocity and c is the velocity of sound in water.
The signal r (t) arriving at the receiver can be expressed as
Figure BDA0002387539540000052
In the formula, τpThe time at which the signal travels along the p-th path to reach the receiver.
The second step is that: construction scale combined filter
The p path is the path with the maximum amplitude, and the received signal caused by Doppler frequency shift and time delay is
Figure BDA0002387539540000053
The impulse response function of the matched filter being designed to match the transmitted signal, i.e.
Figure BDA0002387539540000054
Received signal rp(t) by hrAfter the filter represented by (t), the output signal is
Figure BDA0002387539540000055
The relationship between the maximum value of the output signal and the relative movement speed is as follows
Figure BDA0002387539540000056
In the formula (I), the compound is shown in the specification,
Figure BDA0002387539540000057
Figure BDA0002387539540000058
is the fresnel integral.
At this time
t＝f0(1-ap)/M+apτp
In order to improve the estimation accuracy of the signal arrival time, the scaled signal added in the receiving filter must be a time-shifted signal. The impulse response of the combined scale filter is the sum of the time delay and the different scale expansion of the transmission signal,
Figure BDA0002387539540000061
in the formula, apthIs a scale interval, vthFor the speed interval, satisfy apth＝1+vthC, c is the sound velocity in water,
Figure BDA0002387539540000062
for time delay, vmaxThe maximum relative movement speed of the system is the number of the combined filters
Figure BDA0002387539540000063
Figure BDA0002387539540000064
In order to round up the operator,
Figure BDA0002387539540000065
for the energy normalization coefficient, Σ is the sum operator symbol, the convolution operator symbol, and s (t) is the transmit signal.
The third step: receiver for designing receiving filter based on scale combination
Taking the maximum output value of the matched filter as a test statistic Ts＝ApT，ApIs the amplitude corresponding to the p-th path, T is the duration of the transmitted signal, N-N (0, sigma)2) Is additive white Gaussian noise, N (0, σ)2) Represents a mean of 0 and a variance of σ2The normalized test statistic is as follows,
Figure BDA0002387539540000066
false alarm probability P given according to detection taskFAThe detection threshold is obtained by the Nelman-Pearson criterion
Figure BDA0002387539540000067
Probability of detection
Figure BDA0002387539540000068
In the formula, Q-1(. h) is the inverse of the Q (. h) function,
Figure BDA0002387539540000069
the fourth step: computing test statistics
The echo signal is matched to the ith order component in a combined filter, the combined filter output has a maximum value of
Figure BDA00023875395400000610
To increase the detection probability, the sum of the maximum of the scale combining filter and the maximum of the output of the i-1 th and i +1 th order components is taken as the test statistic.
The fifth step: setting a judgment threshold to carry out detection judgment on the test statistic
And setting a decision threshold to carry out detection decision on the test statistic in the fourth step, wherein in the invention, the maximum value output by the matched filter is considered to be half of the original value (the original value is the maximum value output by the received signal after passing through the traditional matched filter), so that the maximum value at the moment can not be correctly detected.

Claims (1)

1. A synchronous detection method of underwater sound mobile communication is characterized by comprising the following steps:
firstly, emitting a detection signal s (t), wherein the detection signal is an LFM signal, an HFM signal or a pseudo-random modulation signal;
second, construct the scale combination filter, which is the sum of the time delay and the different scale expansion of the transmitting signal, and the impulse response
Figure FDA0002975427790000011
In the formula, vthFor the interval of velocity, satisfy vth＝c(apth-1)，apthIs the scale interval, c is the speed of sound in water,
Figure FDA0002975427790000012
time delay for the p-th path, f0Is the chirp start frequency, M is the chirp frequency modulation slope, vmaxIs the maximum relative speed of movement of the system,number of combined filters
Figure FDA0002975427790000013
Figure FDA0002975427790000014
In order to round up the operator,
Figure FDA0002975427790000015
for energy normalization coefficients, ∑ for sum operator, s*Represents the conjugate of s;
thirdly, designing a receiver based on a scale combination filter, wherein the maximum output value of the receiver is a test statistic Ts＝ApT，ApThe amplitude corresponding to the p path is T, and the duration of the transmitted signal is T; false alarm probability P given according to detection taskFADetection threshold of receiver by Nelman-Pearson criterion
Figure FDA0002975427790000016
Probability of detection
Figure FDA0002975427790000017
Wherein N to N (0, σ)2) Is additive white Gaussian noise, Q-1(. h) is the inverse of the Q (. h) function,
Figure FDA0002975427790000018
fourthly, the echo signal is matched with the ith order component in the scale combination filter, and the maximum value output by the scale combination filter
Figure FDA0002975427790000019
τpThe time delay corresponding to the p path is taken as the test statistic by taking the sum of the maximum value of the scale combination filter, the maximum value output by the (i-1) th order component and the maximum value output by the (i + 1) th order component,
Figure FDA00029754277900000110
represents hr(τp) Convolution with s (t);
and fifthly, setting a judgment threshold to carry out detection judgment on the test statistic in the fourth step.
CN202010103115.9A 2020-02-19 2020-02-19 Synchronous detection method for underwater acoustic mobile communication Active CN111342949B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010103115.9A CN111342949B (en) 2020-02-19 2020-02-19 Synchronous detection method for underwater acoustic mobile communication

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010103115.9A CN111342949B (en) 2020-02-19 2020-02-19 Synchronous detection method for underwater acoustic mobile communication

Publications (2)

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

Family

ID=71185341

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010103115.9A Active CN111342949B (en) 2020-02-19 2020-02-19 Synchronous detection method for underwater acoustic mobile communication

Country Status (1)

Country Link
CN (1) CN111342949B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112152959B (en) * 2020-11-25 2021-02-23 鹏城实验室 Underwater acoustic communication detection signal processing method, device, equipment and storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108243138A (en) * 2018-01-11 2018-07-03 福建星海通信科技有限公司 A kind of combination Doppler estimation suitable for underwater sound communication system
CN109391573A (en) * 2017-08-02 2019-02-26 中国科学院声学研究所 A kind of time-frequency two-dimensional synchronous method based on linear FM signal
CN109474305A (en) * 2018-12-05 2019-03-15 唐山照澜海洋科技有限公司 A kind of 5G multi-carrier spread spectrum underwater acoustic communication method
CN109474295A (en) * 2018-12-24 2019-03-15 南京屹信航天科技有限公司 A kind of baseband processing method for spaceborne measuring and controlling equipment

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8654610B2 (en) * 2009-12-16 2014-02-18 Shb Instruments, Inc. Underwater acoustic navigation systems and methods
CN105323198B (en) * 2014-06-13 2018-08-17 中国科学院声学研究所 A method of carrying out underwater signal transmitting and reception using Hyperbolic Frequency Modulation
GB2555205B (en) * 2016-08-18 2021-07-21 Allen Vanguard Corp System and method for providing adaptive synchronization of LTE communication systems
CN108761384B (en) * 2018-04-28 2022-04-19 西北工业大学 Target positioning method for robust sensor network
CN110247867B (en) * 2019-05-16 2021-03-12 中国科学院声学研究所 Underwater acoustic Doppler estimation method and device, and underwater acoustic communication method and system
CN110535537B (en) * 2019-09-25 2021-04-06 哈尔滨工程大学 Underwater communication and detection integrated method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109391573A (en) * 2017-08-02 2019-02-26 中国科学院声学研究所 A kind of time-frequency two-dimensional synchronous method based on linear FM signal
CN108243138A (en) * 2018-01-11 2018-07-03 福建星海通信科技有限公司 A kind of combination Doppler estimation suitable for underwater sound communication system
CN109474305A (en) * 2018-12-05 2019-03-15 唐山照澜海洋科技有限公司 A kind of 5G multi-carrier spread spectrum underwater acoustic communication method
CN109474295A (en) * 2018-12-24 2019-03-15 南京屹信航天科技有限公司 A kind of baseband processing method for spaceborne measuring and controlling equipment

Also Published As

Publication number Publication date
CN111342949A (en) 2020-06-26

Similar Documents

Publication Publication Date Title
CN106879068B (en) Method for estimating arrival time of signal in strong multipath environment
CN109814094B (en) Multi-target underwater acoustic positioning delay estimation algorithm
CN101594185A (en) The Doppler of mobile water sound communication signal estimates and method for synchronous
CN103176166B (en) Tracking algorithm for time difference of arrival of signals for acoustic passive positioning
CN106019263B (en) Target radial speed measuring method based on more bright spot echo models
CN109765545A (en) A kind of hydrolocation Time Delay Estimation Algorithms based on relevant peaks envelope detected
KR101580427B1 (en) Doppler Frequency Estimation and receiving method for Time-varying Underwater Acoustic Communication Channel
CN109655148A (en) A kind of autonomous extracting method of ship noise non-stationary low frequency spectrum lines
CN110247867B (en) Underwater acoustic Doppler estimation method and device, and underwater acoustic communication method and system
CN104808207A (en) Chaotic underwater sound positioning method
CN110535537A (en) A kind of subsurface communication detection integral method
CN103618575B (en) Automatic real-time frame synchronization method of frequency-modulated hydro-acoustic communication system
CN111342949B (en) Synchronous detection method for underwater acoustic mobile communication
Ma et al. A combined doppler scale estimation scheme for underwater acoustic OFDM system
Mallat et al. UWB based positioning in multipath channels: CRBs for AOA and for hybrid TOA-AOA based methods
CN111580110A (en) Composite code underwater acoustic ranging method based on shallow sea multipath time delay
KR20170029899A (en) Apparatus and method for estimating lfm signal parameter of active sonar system
CN111273294A (en) Underwater acoustic pulse echo ranging method based on time reversal
Diamant Clustering approach for detection and time of arrival estimation of hydrocoustic signals
CN114019453B (en) Ranging method based on underwater acoustic baseline positioning system
CN110673118A (en) Active sonar single-frequency pulse train waveform design and detection algorithm
Yi-Jia et al. Modeling and estimation of the space-time varying channels
CN106646436B (en) A kind of reconnaissance signal method for parameter estimation based on signal wide and narrow strip fuzziness
CN109768812B (en) Underwater acoustic communication Doppler estimation and synchronization method based on chaotic frequency modulation
Liu et al. A positive and negative HFM for speed measurement

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