RU2700005C1 - Method of estimating channel parameters in broadband hydroacoustic communication and a device for realizing said channel - Google Patents

Method of estimating channel parameters in broadband hydroacoustic communication and a device for realizing said channel Download PDF

Info

Publication number
RU2700005C1
RU2700005C1 RU2018139703A RU2018139703A RU2700005C1 RU 2700005 C1 RU2700005 C1 RU 2700005C1 RU 2018139703 A RU2018139703 A RU 2018139703A RU 2018139703 A RU2018139703 A RU 2018139703A RU 2700005 C1 RU2700005 C1 RU 2700005C1
Authority
RU
Russia
Prior art keywords
reference sequences
likelihood
extended
spectrum
reference sequence
Prior art date
Application number
RU2018139703A
Other languages
Russian (ru)
Inventor
Алексей Сергеевич Шустов
Тимур Георгиевич Келин
Original Assignee
федеральное государственное бюджетное образовательное учреждение высшего образования "Санкт-Петербургский государственный морской технический университет" (СПбГМТУ)
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 федеральное государственное бюджетное образовательное учреждение высшего образования "Санкт-Петербургский государственный морской технический университет" (СПбГМТУ) filed Critical федеральное государственное бюджетное образовательное учреждение высшего образования "Санкт-Петербургский государственный морской технический университет" (СПбГМТУ)
Priority to RU2018139703A priority Critical patent/RU2700005C1/en
Application granted granted Critical
Publication of RU2700005C1 publication Critical patent/RU2700005C1/en

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems

Landscapes

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

Abstract

FIELD: information technology.
SUBSTANCE: invention relates to information transmission systems. Technical result of the proposed technical solution is achieved due to that in the transmitter of the information transmission system generate the spread spectrum reference sequences in the frequency subbands of the broadband channel, and then formed spreading reference sequences are combined with information part of transmitted symbol. At the receiving side of the information transmission system, reference sequences are formed for the expected range of Doppler scaling of the received reference sequence with a spread spectrum, coherent detection of reference sequences in the received signal is carried out.
EFFECT: high accuracy of estimating time of arrival and a coefficient of Doppler zooming for information symbols in a data transmission system in broadband hydroacoustic channels with a short coherence interval and/or high frequency selectivity, which increases reliability of communication, as well as wider range of application conditions for hydroacoustic communication systems.
25 cl, 3 dwg

Description

Изобретение относится к области систем передачи информации и предназначено для оценки времени прихода и коэффициента доплеровского масштабирования информационного символа в широкополосном гидроакустическом канале.The invention relates to the field of information transmission systems and is intended to assess the time of arrival and Doppler scaling factor of an information symbol in a broadband sonar channel.

Задача оценки параметров канала возникает, например, при обмене данными между объектами, находящимися в водной среде, когда в качестве канала передачи информации выступает широкополосный гидроакустический канал. При передаче информационных символов в широкополосном гидроакустическом канале возникают их искажения, обусловленные доплеровским масштабированием и задержкой распространения. Оценка параметров этих искажений на приемной стороне необходима для успешного приема и демодуляции информационных символов, и, как следствие, для декодирования переданного сообщения.The task of estimating the channel parameters arises, for example, in the exchange of data between objects located in the aquatic environment, when a broadband sonar channel acts as an information transmission channel. When transmitting information symbols in a broadband sonar channel, their distortions arise due to Doppler scaling and propagation delay. The estimation of the parameters of these distortions at the receiving side is necessary for the successful reception and demodulation of information symbols, and, as a result, for decoding the transmitted message.

Известны способ и устройство, обеспечивающие оценку времени прихода и коэффициента доплеровского масштабирования для группы информационных символов, описанные в патенте США 7859944 от 28.12.2010 (Н04В 11/00), где в начале и в конце группы информационных символов передают опорные последовательности "Preamble" и "Postamble" соответственно. При приеме за начало группы символов принимают момент обнаружения опорной последовательности "Preamble", а коэффициент доплеровского масштабирования определяют как функцию от временного интервала между моментом обнаружения опорной последовательности "Preamble" и моментом обнаружения опорной последовательности "Postamble". Однако данное техническое решение не может быть в полной мере использовано в гидроакустических каналах с малым временем когерентности, поскольку в таких каналах единые оценки времени прихода и коэффициента доплеровского масштабирования, полученные для группы символов, будут обладать значительной погрешностью для каждого символа в группе. Наиболее близкими к заявленному решению являются способ и устройство по патенту США 8023596 от 20.09.2011 (H04L 27/00), в котором в информационных символах осуществляется передача опорных несущих с фиксированными частотами, мощности которых существенно превышают мощности несущих, используемых для передачи информации. На приемной стороне осуществляется оценка коэффициента доплеровского масштабирования за счет фильтрации набором полосовых фильтров каждой из несущих с повышенной мощностью, выделение огибающих за счет некогерентного детектирования сигнала на выходе каждого фильтра и последующей комбинации значений с выходов некогерентных детекторов. Однако данное техническое решение не может обеспечить требуемую функциональность, т.к. не позволяет работать в гидроакустических каналах с высокой частотной избирательностью, поскольку оценки энергетического параметра, полученные при некогерентном детектировании, искажаются влиянием неравномерной Амплитудно-частотной Характеристики (АЧХ) канала. К другим недостаткам данного технического решения можно отнести невозможность его работы в гидроакустических каналах с временем когерентности, соизмеримым с длительностью символа, поскольку даже незначительная потеря ортогональности, вызванная погрешностью оценки коэффициента доплеровского масштабирования, будет приводить к проникновению внеполосных спектральных составляющих от несущих с повышенным уровнем мощности на частоты, используемые для передачи информации, и, как следствие, к искажению передаваемой информации. Также данное техническое решение имеет ограничение на применение в многолучевых гидроакустических каналах, в которых пути распространения сигнала имеют разный коэффициент доплеровского масштабирования, так называемый, Doppler spread, поскольку, даже будучи ослабленными на пути распространения, смещенные по частоте несущие с повышенным уровнем мощности будут вносить помехи на частотах, используемых для передачи информации.A known method and apparatus for estimating the arrival time and Doppler scaling factor for a group of information symbols described in US patent 7859944 dated 12/28/2010 (HB04/00), where at the beginning and at the end of the group of information symbols transmit reference sequences "Preamble" and "Postamble" respectively. When a group of symbols is taken as the beginning, the moment of detection of the reference sequence "Preamble" is taken, and the Doppler scaling factor is determined as a function of the time interval between the moment of detection of the reference sequence "Preamble" and the moment of detection of the reference sequence "Postamble". However, this technical solution cannot be fully used in hydroacoustic channels with a short coherence time, since in such channels the unified estimates of the arrival time and Doppler scaling coefficient obtained for a group of symbols will have a significant error for each symbol in the group. Closest to the claimed solution are the method and device according to US patent 8023596 from 09.20.2011 (H04L 27/00), in which information carriers transmit reference carriers with fixed frequencies, the powers of which significantly exceed the powers of the carriers used to transmit information. On the receiving side, the Doppler scaling coefficient is estimated by filtering with a set of band-pass filters of each of the carriers with increased power, envelopes are selected due to incoherent signal detection at the output of each filter and the subsequent combination of values from the outputs of incoherent detectors. However, this technical solution cannot provide the required functionality, because it does not allow working in sonar channels with high frequency selectivity, since estimates of the energy parameter obtained with incoherent detection are distorted by the influence of the non-uniform amplitude-frequency characteristics (AFC) of the channel. Other disadvantages of this technical solution include the impossibility of its operation in sonar channels with a coherence time commensurate with the symbol duration, since even a slight loss of orthogonality caused by an error in the estimation of the Doppler scaling coefficient will lead to the penetration of out-of-band spectral components from carriers with an increased power level at frequencies used to transmit information, and, as a result, to distortion of the transmitted information. Also, this technical solution has a restriction on the use in multipath hydroacoustic channels in which the signal propagation paths have different Doppler scaling factors, the so-called Doppler spread, because even if weakened along the propagation path, frequency-shifted carriers with an increased power level will introduce interference at frequencies used to transmit information.

Техническим результатом заявляемого технического решения является повышение точности оценки времени прихода и коэффициента доплеровского масштабирования для информационных символов в системе передачи данных в широкополосных гидроакустических каналах с малым интервалом когерентности и/или с высокой частотной избирательностью, повышение надежности коммуникации, а также расширение диапазона применимости для систем гидроакустической связи.The technical result of the proposed technical solution is to increase the accuracy of estimating the time of arrival and Doppler scaling factor for information symbols in a data transmission system in broadband sonar channels with a small coherence interval and / or with high frequency selectivity, increasing the reliability of communication, as well as expanding the range of applicability for sonar systems communication.

Технический результат достигается тем, что в передатчике системы передачи информации формируют опорные последовательности с расширенным спектром в частотных поддиапазонах широкополосного канала, а затем сформированные опорные последовательности с расширенным спектром комбинируют с информационной частью передаваемого символа. На приемной стороне системы передачи информации формируют эталонные последовательности для ожидаемого диапазона доплеровского масштабирования принимаемой опорной последовательности с расширенным спектром, производят когерентное детектирование эталонных последовательностей в принимаемом сигнале, а в случае, когда для опорной последовательности с расширенным спектром формируют эталонный ряд отсчетов, производят формирование исходного ряда отсчетов для принимаемого сигнала, формирование скорректированных рядов отсчетов из исходного ряда отсчетов для ожидаемого диапазона доплеровского масштабирования принимаемых опорных последовательностей с расширенным спектром, после этого осуществляют когерентное детектирование эталонных рядов отсчетов в каждом из скорректированных рядов отсчетов, затем, как и в случае после когерентного детектирования эталонных последовательностей в принятом сигнале, производят формирование оценок правдоподобия, комбинирование оценок правдоподобия и определяют оценку времени прихода и величины доплеровского масштабирования символа, а также значений правдоподобия оценок.The technical result is achieved by the fact that in the transmitter of the information transmission system, reference sequences with a spread spectrum are formed in the frequency subbands of the broadband channel, and then the formed reference sequences with a spread spectrum are combined with the information part of the transmitted symbol. On the receiving side of the information transmission system, reference sequences are formed for the expected range of Doppler scaling of the received spread spectrum reference sequence, the reference sequences are coherently detected in the received signal, and in the case when a reference series of samples are formed for the spread spectrum reference sequence, the initial series is generated samples for the received signal, the formation of the corrected series of samples from the outcome a series of samples for the expected range of Doppler scaling of the received extended-spectrum reference sequences, then coherent detection of the reference series of samples in each of the corrected series of samples is carried out, then, as in the case after coherent detection of the reference sequences in the received signal, likelihood estimates are generated, combining likelihood estimates and determining the estimated time of arrival and magnitude of Doppler scaling symbols la, as well as the likelihood values of the estimates.

В отличие от известного способа оценки коэффициента доплеровского масштабирования в широкополосных гидроакустических каналах включающего формирование и передачу сигнала с использованием несущих, на которых передается информационная посылка, а также набор опорных несущих с фиксированными частотами, уровни которых существенно превышают уровни несущих, используемых для передачи информационной посылки, фильтрацию набором фильтров каждой из опорных несущих с фиксированной частотой и повышенным уровнем, формирование значений огибающих для сигналов на выходе каждого фильтра, комбинирование значений огибающих, в заявляемом способе оценки времени прихода сигнала и коэффициента доплеровского масштабирования в широкополосных гидроакустических каналахIn contrast to the known method for estimating the Doppler scaling coefficient in broadband sonar channels, which includes generating and transmitting a signal using carriers on which the information packet is transmitted, as well as a set of reference carriers with fixed frequencies, the levels of which significantly exceed the levels of carriers used to transmit the information packet, filtering by a set of filters of each of the supporting carriers with a fixed frequency and a high level, the formation of envelope values x for the signals at the output of each filter, combining the envelope values, in the inventive method for estimating the signal arrival time and the Doppler scaling factor in broadband sonar channels

- формируют опорные последовательности с расширенным спектром в частотных поддиапазонах широкополосного канала,- form reference sequences with an extended spectrum in the frequency subbands of the broadband channel,

- комбинируют опорные последовательности с расширенным спектром с информационной частью передаваемого символа,- combine reference sequences with an extended spectrum with the information part of the transmitted symbol,

- формируют эталонные последовательности на приемной стороне,- form reference sequences on the receiving side,

- осуществляют когерентное детектирование эталонных последовательностей в принимаемом сигнале,- carry out coherent detection of reference sequences in the received signal,

- формируют оценки правдоподобия,- form the likelihood estimates,

- комбинируют оценки правдоподобия.- combine likelihood estimates.

При этом эталонные последовательности формируют для ожидаемого диапазона коэффициентов доплеровского масштабирования принимаемых опорных последовательностей с расширенным спектром.In this case, reference sequences are formed for the expected range of Doppler scaling coefficients of the received extended-spectrum reference sequences.

Когерентное детектирование эталонных последовательностей в принимаемом сигнале производят как во временной, так и в частотной областях.Coherent detection of reference sequences in the received signal is performed both in the time and in the frequency domains.

В случае использования отсчетов при приеме для вычисления оценок правдоподобия и коэффициента доплеровского масштабирования в широкополосных гидроакустических каналах в принимаемом сигнале для каждой опорной последовательности с расширенным спектромIn the case of using reception samples to calculate the likelihood estimates and the Doppler scaling coefficient in broadband sonar channels in the received signal for each reference sequence with an extended spectrum

- формируют эталонный ряд отсчетов,- form a reference series of readings,

- формируют исходный ряд отсчетов, из которого- form the initial series of samples from which

- формируют скорректированные ряды отсчетов, затем- form the corrected series of readings, then

- когерентно детектируют эталонный ряд отсчетов в каждом из скорректированных рядов отсчетов,- coherently detect a reference series of samples in each of the adjusted series of samples,

- формируют оценки правдоподобия,- form the likelihood estimates,

- комбинируют оценки правдоподобия.- combine likelihood estimates.

При этом скорректированные ряды отсчетов формируют для ожидаемого диапазона доплеровского масштабирования принимаемой опорной последовательностей с расширенным спектром.In this case, the corrected series of samples form for the expected range of Doppler scaling of the received reference sequences with extended spectrum.

Чертежи, поясняющие заявляемый способ, представлены наThe drawings explaining the claimed method are presented on

Фиг. 1 - структура передаваемого Символа 100.FIG. 1 - structure of the transmitted Symbol 100.

Фиг. 2 - техническая реализация приемной подсистемы для оценивания времени прихода и доплеровского масштабирования информационного символа.FIG. 2 - technical implementation of the receiving subsystem for estimating the arrival time and Doppler scaling of the information symbol.

Фиг. 3 - техническая реализация блока детектирования.FIG. 3 is a technical implementation of the detection unit.

Заявляемый способ оценки времени прихода и коэффициента доплеровского масштабирования в широкополосном канале передачи информации работает следующим образом:The inventive method of estimating the time of arrival and the Doppler scaling factor in a broadband information channel works as follows:

- на передающей стороне формируют N Опорных Последовательностей с Расширенным Спектром 131…13N в частотных поддиапазонах широкополосного канала, затем комбинируют Опорные Последовательности с Расширенным Спектром 131…13N с Информационной Частью 111…11М передаваемого Символа 100 (фиг. 1). На приемной стороне:- N Reference Sequences with Extended Spectrum 131 ... 13N in the frequency subbands of the broadband channel are formed on the transmitting side, then Reference Sequences with Extended Spectrum 131 ... 13N with Information Part 111 ... 11M of transmitted Symbol 100 are combined (Fig. 1). On the receiving side:

- формируют эталонные последовательности для ожидаемого диапазона доплеровского масштабирования принимаемых Опорных Последовательностей с Расширенным Спектром;- form reference sequences for the expected range of Doppler scaling of the received Reference Sequences with Extended Spectrum;

- осуществляют когерентное детектирование эталонных последовательностей в принимаемом сигнале;- carry out coherent detection of reference sequences in the received signal;

- формируют оценки правдоподобия;- form a likelihood assessment;

- комбинируют оценки правдоподобия.- combine likelihood estimates.

При выборе Опорных Последовательностей с Расширенным Спектром 131…13N учитывают ширину их полосы для обеспечения возможно меньшего уровня внеполосных спектральных составляющих, а также форму их функции неопределенности, обеспечивающей необходимое разрешение по времени прихода и доплеровскому масштабированию. Длительность Опорных Последовательностей с Расширенным Спектром 131…13N может быть ограничена сверху временем когерентности канала передачи, а ширина полосы каждой последовательности может быть ограничена сверху шириной полосы когерентности канала передачи. Вариантами таких последовательностей могут быть последовательности с угловой модуляцией (фазовой или частотной), псевдослучайные последовательности, а также последовательности, сформированные с использованием мультиплексирования с ортогональным частотным разделением каналов (см. Nadav Levanon, Eli Mozeson; Radar Signals, ISBN: 978-0-471-47378-7, Jul 2004, Wiley-IEEE Press, 432 pages).When choosing Reference Sequences with an Extended Range of 131 ... 13N, their bandwidth is taken into account to ensure the lowest possible level of out-of-band spectral components, as well as the shape of their uncertainty function, which provides the necessary resolution for arrival time and Doppler scaling. The duration of the Reference Sequences with the Extended Spectrum 131 ... 13N may be limited from above by the coherence time of the transmission channel, and the bandwidth of each sequence may be limited from above by the bandwidth of the coherence of the transmission channel. Options for such sequences may include angular modulation (phase or frequency) sequences, pseudo-random sequences, and sequences generated using orthogonal frequency division multiplexing (see Nadav Levanon, Eli Mozeson; Radar Signals, ISBN: 978-0-471 -47378-7, Jul 2004, Wiley-IEEE Press, 432 pages).

Для исключения влияния нулей диаграммы направленности при использовании передающей системы с несколькими излучателями декоррелируют сигналы, передаваемые разными излучателями, используя разные циклические сдвиги Опорных Последовательностей с Расширенным Спектром 131…13N одновременно передаваемых разными излучателями.To exclude the influence of zeros of the radiation pattern when using a transmitting system with several emitters, signals transmitted by different emitters are decorrelated using different cyclic shifts of the Reference Sequences with Extended Spectrum 131 ... 13N simultaneously transmitted by different emitters.

Для уменьшения влияния межсимвольных искажений при передаче последовательности из нескольких Символов 100 в канале с многолучевым распространением используют повторение части Опорных Последовательностей с Расширенным Спектром 131…13N для создания циклического префикса.To reduce the influence of intersymbol distortions when transmitting a sequence of several Symbols 100 in a multipath channel, use the repetition of part of the Reference Sequences with Extended Spectrum 131 ... 13N to create a cyclic prefix.

Гидроакустические широкополосные излучатели и гидрофоны отличает высокая неравномерность АЧХ. Для компенсации влияния априорно известных параметров передающего и/или приемного блока на качество оценки времени прихода и коэффициента доплеровского масштабирования возможно изменение количества Опорных Последовательностей с Расширенным Спектром 131…13N, значений их центральных частот, а также их спектральных плотностей мощности.Hydroacoustic broadband emitters and hydrophones are characterized by high frequency response unevenness. To compensate for the influence of a priori known parameters of the transmitting and / or receiving unit on the quality of estimation of the arrival time and Doppler scaling coefficient, it is possible to change the number of Reference Sequences with Extended Spectrum 131 ... 13N, the values of their center frequencies, as well as their power spectral densities.

Возможным методом комбинирования Опорных Последовательностей 131…13N с Информационной Частью 111…11М является метод, при котором для передачи Символа 100 используют мультиплексирование с ортогональным частотным разделением каналов (OFDM), при этом часть несущих выделяют для передачи Опорных Последовательностей 131…13N, а другую часть несущих - для передачи Информационной Части 111…11М. На выходе Блока Комбинирования формируется Символ 100, который поступает в широкополосный гидроакустический канал.A possible method of combining the Reference Sequences 131 ... 13N with the Information Part 111 ... 11M is the method in which orthogonal frequency division multiplexing (OFDM) is used to transmit Symbol 100, while some of the carriers are allocated for transmitting the Reference Sequences 131 ... 13N, and the other part carriers - for transmission of the Information Part 111 ... 11M. At the output of the Combination Unit, Symbol 100 is formed, which enters the broadband sonar channel.

Опорные Последовательности с Расширенным Спектром 131…13N также используют в качестве пилот-сигналов при демодуляции Информационной Части 111…11М Символа 100 для чего осуществляют когерентное комбинирование Опорных Последовательностей с Расширенным Спектром 131…13N с Информационной Частью 111…11М передаваемого Символа 100.Reference Sequences with Extended Spectrum 131 ... 13N are also used as pilot signals for demodulation of the Information Part 111 ... 11M of Symbol 100 for which coherent combination of Reference Sequences with Extended Spectrum 131 ... 13N with Information Part 111 ... 11M of transmitted Symbol 100 is performed.

Для уменьшения влияния внеполосных составляющих от Информационной Части 111…11М на Опорные Последовательности с Расширенным Спектром 131…13N и внеполосных составляющих от Опорных Последовательностей с Расширенным Спектром 131…13N на Информационную Часть 111…11М, которые возникают при потере ортогональности между отдельными несущими при приеме Символа 100 за счет влияния флуктуаций канала и ошибок синхронизации, и, как следствие, для улучшения качества оценки и/или демодуляции, возможно формирование Защитных Частотных Интервалов 121…12N и 141…14N между Опорными Последовательностями с Расширенным Спектром 131…13N и Информационной Частью 111…11М.To reduce the influence of out-of-band components from the Information Part 111 ... 11M on the Reference Sequences with the Extended Spectrum 131 ... 13N and out-of-band components from the Reference Sequences with the Extended Spectrum 131 ... 13N on the Information Part 111 ... 11M, which occur when the orthogonality between the individual carriers is lost when receiving the Symbol 100 due to the influence of channel fluctuations and synchronization errors, and, as a result, to improve the quality of estimation and / or demodulation, it is possible to form Protective Frequency Intervals 121 ... 12N and 141 ... 14N between Reference Sequences with Extended Spectrum 131 ... 13N and Information Part 111 ... 11M.

На приемной стороне для каждой из Опорных Последовательностей с Расширенным Спектром 131…13N на выходах Формирователей Эталонных Последовательностей 321…32K формируют набор эталонных последовательностей, каждая из которых в наборе может представлять собой соответствующую Опорную Последовательность с Расширенным Спектром, формируемую в передатчике, и подвергнутую доплеровскому масштабированию, а именно, если коэффициент доплеровского масштабирования равен (1+α), то каждая частотная компонента с частотой ƒ передаваемой Опорной Последовательности с Расширенным Спектром будет преобразована в (1+α)ƒ, а длительность Т переданной Опорной Последовательности с Расширенным Спектром будет преобразована в (1-α)Т. Для Опорной Последовательности с Расширенным Спектром с номером n=1…N количество Kn эталонных последовательностей Unk), k=1…Kn в наборе может определяться необходимой разрешающей способностью по коэффициенту доплеровского масштабирования Δnk+1k и ожидаемым диапазоном значений коэффициента доплеровского масштабирования (1+αmin)≤(1+αk)≤(1+αmах).On the receiving side, for each of the Reference Sequences with the Extended Spectrum 131 ... 13N, at the outputs of the Reference Sequencers 321 ... 32K, a set of reference sequences is formed, each of which in the set can be a corresponding Reference Sequence with the Extended Spectrum generated in the transmitter and subjected to Doppler scaling namely, if the Doppler scaling factor is (1 + α), then each frequency component with a frequency ƒ of transmitted Reference Sequence The extended spectrum will be converted to (1 + α) ƒ, and the duration T of the transmitted extended spectrum will be converted to (1-α) T. For the Reference Sequence with the Extended Spectrum with number n = 1 ... N, the number K n of reference sequences U nk ), k = 1 ... K n in the set can be determined by the necessary resolution by the Doppler scaling coefficient Δ n = α k + 1 - α k and the expected range of values of the Doppler scaling coefficient (1 + α min ) ≤ (1 + α k ) ≤ (1 + α max ).

При наличии априорно известной информации о параметрах канала на приемной стороне ее можно учитывать при детектировании эталонных последовательностей в принимаемом сигнале. Так, например, известность уточненных диапазонов времени прихода и/или коэффициента доплеровского масштабирования сокращает пространство поиска за счет сокращения количества эталонных последовательностей и, как следствие, уменьшает вероятность ложного обнаружения Опорной Последовательности, а также сложность приемной аппаратуры, что в свою очередь сокращает расход энергии и продливает время автономной работы системы. Диапазон значений коэффициента доплеровского масштабирования может зависеть от взаимной радиальной скорости передающей и приемной антенн, а также от скорости движения среды в радиальном направлении. Возможными значениями, формируемыми в моменты времени τ на выходе Блоков Когерентного Детектирования 311…31K, являются значения взаимно-корреляционной функции Rnk(τ) принимаемого сигнала и k-й эталонной последовательности Unk), сформированной для n-й Опорной Последовательности с Расширенным Спектром 131…13N. Таким образом для одного Символа 100 количество значений Rnk(τ), вычисляемых в Блоках Когерентного Детектирования 311..31K для N Опорных Последовательностей с Расширенным Спектром 131..13N, составляет

Figure 00000001
. Когерентное детектирование эталонных последовательностей в принимаемом сигнале может производиться как во временной, так и в частотной областях.If there is a priori known information about the channel parameters on the receiving side, it can be taken into account when detecting reference sequences in the received signal. For example, the popularity of the specified arrival time ranges and / or the Doppler scaling factor reduces the search space by reducing the number of reference sequences and, as a result, reduces the probability of false detection of the Reference Sequence, as well as the complexity of the receiving equipment, which in turn reduces energy consumption and extends the battery life of the system. The range of values of the Doppler scaling coefficient may depend on the mutual radial velocity of the transmitting and receiving antennas, as well as on the speed of the medium in the radial direction. Possible values generated at time instants τ at the output of Coherent Detection Units 311 ... 31K are the values of the cross-correlation function R nk (τ) of the received signal and the k-th reference sequence U nk ) generated for the nth Reference Sequence with Extended Spectrum 131 ... 13N. Thus, for one Symbol 100, the number of R nk (τ) values calculated in the Coherent Detection Blocks 311..31K for N Reference Sequences with the Extended Spectrum 131..13N, is
Figure 00000001
. Coherent detection of reference sequences in the received signal can be performed both in the time and in the frequency domains.

В случае использования отсчетов при приеме для вычисления значений Rnk(τ) для каждой n-й Опорной Последовательности с Расширенным Спектром 131..13N формируют эталонный ряд отсчетов с частотой отсчетов Fs;In the case of using samples at reception to calculate the values of R nk (τ) for each n-th Reference Sequence with Extended Spectrum 131..13N, a reference series of samples with a sampling frequency F s is formed ;

- для принимаемого сигнала формируют исходный ряд отсчетов с частотой отсчетов Fs, используя, например, аналогово-цифровой преобразователь (АЦП);- for the received signal form the initial series of samples with a sampling frequency F s using, for example, an analog-to-digital converter (ADC);

- из исходного ряда отсчетов формируют Kn скорректированных рядов отсчетов;- from the original series of samples form K n adjusted series of samples;

- когерентно детектируют эталонный ряд отсчетов в каждом из скорректированных рядов отсчетов, вычисляя в моменты времени τ значения взаимно-корреляционной функции Rnk(τ), где k=1…Kn;- coherently detect the reference series of samples in each of the corrected series of samples, calculating at times τ the values of the cross-correlation function R nk (τ), where k = 1 ... K n ;

при этом скорректированные ряды отсчетов формируют для ожидаемого диапазона доплеровского масштабирования принимаемой опорной последовательностей с расширенным спектром, а частота отсчетов в k-м скорректированном ряде отсчетов составляет (1+αk)Fs.in this case, the corrected series of samples form for the expected range of Doppler scaling of the received reference sequences with an extended spectrum, and the frequency of samples in the k-th corrected series of samples is (1 + α k ) F s .

Для формирования из исходного ряда отсчетов k-го скорректированного ряда отсчетов с измененной в (1+αk) раз частотой отсчетов можно использовать полиномиальную интерполяцию отсчетов из исходного ряда отсчетов.To form the kth corrected series of samples from the original series of samples with the sample frequency changed by (1 + α k ) times, you can use polynomial interpolation of samples from the original series of samples.

Оценкой правдоподобия, вычисляемой в каждом из формирователей Оценки Правдоподобия 331…33K, для входных значений взаимно-корреляционной функции Rnk(τ) при низких отношениях сигнал-шум в канале является значение

Figure 00000002
. Для высоких отношений сигнал-шум в канале оценкой правдоподобия является значение
Figure 00000003
, где Еk - энергия k-й эталонной последовательности Unk), которая может быть вычислена заранее, а
Figure 00000004
- оценка энергии n-й принятой опорной последовательности в моменты времени τ с учетом влияния коэффициента доплеровского масштабирования (1+αk), т.е. с учетом изменения длительности в (1-αk) раз. Большее значение Qnk(τ) показывает, что n-я опорная последовательность в принятом сигнале с большей достоверностью соответствует эталонной последовательности Unk).The likelihood estimate calculated in each of the formers of the Likelihood Assessment 331 ... 33K for the input values of the cross-correlation function R nk (τ) at low signal-to-noise ratios in the channel is the value
Figure 00000002
. For high signal-to-noise ratios in a channel, the likelihood estimate is
Figure 00000003
where E k is the energy of the kth reference sequence U nk ), which can be calculated in advance, and
Figure 00000004
- estimation of the energy of the nth received reference sequence at time instants τ taking into account the influence of the Doppler scaling coefficient (1 + α k ), i.e. taking into account the change in duration by (1-α k ) times. A larger value of Q nk (τ) indicates that the nth reference sequence in the received signal corresponds with greater certainty to the reference sequence U nk ).

Для уменьшения вероятности ложного обнаружения Опорной Последовательности с Расширенным Спектром при приеме и, как следствие, для повышения точности оценивания времени прихода и коэффициента доплеровского масштабирования значения Оценок Правдоподобия сравниваются с порогом, при этом значение Оценки Правдоподобия поступает в Блок Комбинирования Оценок Правдоподобия 270 в том случае, если порог превышен. Такая пороговая обработка обеспечивает выделение и последующее комбинирование оценок только с высокой достоверностью в Блоке Комбинирования Оценок Правдоподобия 270.To reduce the likelihood of false detection of the Reference Sequence with the Extended Spectrum during reception and, as a result, to increase the accuracy of estimating the time of arrival and the Doppler scaling coefficient, the likelihood estimates are compared with a threshold, and the likelihood estimates are sent to the Combination of Likelihood Assessment Block 270 if if the threshold is exceeded. Such threshold processing ensures the selection and subsequent combination of estimates with only high reliability in the Combination Block of Credibility Estimates 270.

В Блоке Комбинирования Оценок Правдоподобия 270 для комбинирования оценок правдоподобия, получаемых с Формирователей Оценки Правдоподобия 331…33K для каждой из принятых Опорных Последовательностей с Расширенным Спектром 131…13N и набора эталонных последовательностей, и формирования оценки времени прихода

Figure 00000005
, оценки коэффициента доплеровского масштабирования
Figure 00000006
, а также значения правдоподобия оценок
Figure 00000007
для принятого информационного символа вычисляется взвешенное среднее значение. Для его вычисления множество значений времени прихода и коэффициента доплеровского масштабирования усредняется с весами, равными значениям их оценок правдоподобия т.е решается задача о нахождении координат центра масс
Figure 00000008
и эквивалентной массы
Figure 00000007
для множества точек Qnk(τ) на плоскости {τ, α}.In the Combination of Likelihood Assessment Unit 270 for combining the likelihood estimates obtained from the Likelihood Estimator formers 331 ... 33K for each of the accepted Reference Sequences with Extended Spectrum 131 ... 13N and a set of reference sequences, and forming an estimate of the time of arrival
Figure 00000005
Doppler scaling factor estimates
Figure 00000006
and likelihood estimates
Figure 00000007
for the received information symbol, a weighted average value is calculated. To calculate it, the set of values of the arrival time and the Doppler scaling coefficient is averaged with weights equal to the values of their likelihood estimates, i.e., the problem of finding the coordinates of the center of mass is solved
Figure 00000008
and equivalent mass
Figure 00000007
for the set of points Q nk (τ) on the plane {τ, α}.

Другим возможным способом комбинирования в Блоке Комбинирования Оценок Правдоподобия 270 является нахождение максимального значения среди оценок правдоподобия Qnk(τ). В этом случае оценки времени прихода

Figure 00000005
и коэффициента доплеровского масштабирования
Figure 00000006
будут соответствовать найденному максимальному значению оценки правдоподобия
Figure 00000007
. Уточнение оценки времени прихода
Figure 00000005
и оценки коэффициента доплеровского масштабирования
Figure 00000006
осуществляют за счет интерполяции значений оценок правдоподобия в окрестности максимума и последующего нахождения максимума среди интерполированных значений.Another possible combination method in the Combination Block of Credibility Estimates 270 is to find the maximum value among the likelihood estimates Q nk (τ). In this case, the estimated time of arrival
Figure 00000005
and Doppler scaling factor
Figure 00000006
will match the maximum likelihood score found
Figure 00000007
. Refinement of the estimated time of arrival
Figure 00000005
and estimates of the Doppler scaling factor
Figure 00000006
carried out by interpolating the values of the likelihood estimates in the vicinity of the maximum and the subsequent finding of the maximum among the interpolated values.

Также в Блоке Комбинирования Оценок Правдоподобия 270 сначала могут формироваться оценки времени прихода и коэффициента доплеровского масштабирования, а также значение правдоподобия этих оценок отдельно для каждой из принятых Опорных Последовательностей с Расширенным Спектром 131…13N, а затем формироваться окончательная общая оценка за счет комбинирования значений из всех N поддиапазонов широкополосного канала передачи информации. При этом для комбинирования оценок правдоподобия для каждой из принятых Опорных Последовательностей с Расширенным Спектром 131…13N может использоваться один метод, а для комбинирования значений из всех N поддиапазонов широкополосного канала передачи информации - другой.Also, in the Combination Block of Likelihood Estimates 270, first, estimates of the time of arrival and Doppler scaling coefficient can be formed, as well as the likelihood of these estimates separately for each of the accepted Reference Sequences with Extended Spectrum 131 ... 13N, and then the final general estimate can be formed by combining the values from N subbands of a broadband information channel. At the same time, one method can be used to combine the likelihood estimates for each of the accepted Reference Sequences with the Extended Spectrum 131 ... 13N, and the other to combine the values from all N subbands of the broadband information channel.

Устройство для оценки коэффициента доплеровского масштабирования в широкополосных гидроакустических каналах по патенту-прототипу (см. патент США 8023596 от 20.09.2011 (H04L 27/00)) содержит:A device for evaluating the coefficient of Doppler scaling in broadband sonar channels according to the patent prototype (see US patent 8023596 from 09.20.2011 (H04L 27/00)) contains:

- формирователь и передатчик. сигнала, включающего несущие, на которых передается информационная посылка, а также набор опорных несущих с фиксированными частотами, уровни которых существенно превышают уровни несущих, используемых для передачи информационной посылки;- shaper and transmitter. a signal including carriers on which the information package is transmitted, as well as a set of reference carriers with fixed frequencies, the levels of which significantly exceed the levels of carriers used to transmit the information package;

- наборы с одинаковым количеством полосовых фильтров для каждой из опорных несущих с фиксированной частотой и повышенным уровнем;- sets with the same number of band-pass filters for each of the reference carriers with a fixed frequency and an increased level;

- квадратичные детекторы огибающей;- quadratic envelope detectors;

- устройство комбинирования значений огибающих,- a device for combining envelope values,

причем принимаемый сигнал поступает на входы фильтров, выходы фильтров подключены ко входам квадратичных детекторов огибающей, а выходы квадратичных детекторов огибающей подключены к устройству комбинирования, при этом общая полоса фильтров в каждом наборе перекрывает возможный диапазон изменения частоты соответствующей опорной несущей, вызванного доплеровским масштабированием.moreover, the received signal is fed to the inputs of the filters, the outputs of the filters are connected to the inputs of the quadratic envelope detectors, and the outputs of the quadratic envelope detectors are connected to the combiner, while the common filter band in each set covers the possible frequency range of the corresponding reference carrier caused by Doppler scaling.

Техническим результатом предлагаемого устройства является повышение точности при реализации заявляемого способа оценки времени прихода и коэффициента доплеровского масштабирования для информационных символов в системе передачи данных в широкополосных гидроакустических каналах с малым интервалом когерентности и/или с высокой частотной избирательностью, повышение надежности коммуникации, а также расширение диапазона применимости для систем гидроакустической связи.The technical result of the proposed device is to increase the accuracy when implementing the proposed method for estimating the arrival time and Doppler scaling factor for information symbols in a data transmission system in broadband sonar channels with a small coherence interval and / or with high frequency selectivity, increasing the reliability of communication, as well as expanding the range of applicability for sonar communication systems.

Технический результат достигается тем, что в устройство введены:The technical result is achieved by the fact that the device introduced:

- формирователи опорных последовательностей с расширенным спектром в частотных поддиапазонах широкополосного канала;- shapers of reference sequences with extended spectrum in the frequency subbands of the broadband channel;

- блок комбинирования опорных последовательностей с расширенным спектром с информационной частью передаваемого символа;- block combining reference sequences with an extended spectrum with the information part of the transmitted symbol;

- формирователи эталонных последовательностей на приемной стороне;- shapers of reference sequences on the receiving side;

- блоки когерентного детектирование эталонных последовательностей в принимаемом сигнале;- blocks of coherent detection of reference sequences in the received signal;

- формирователи оценок правдоподобия;- likelihood estimators;

- блок комбинирования оценок правдоподобия;- a block combining likelihood ratings;

причем выходы формирователей опорных последовательностей с расширенным спектром подключены к входам блока комбинирования опорных последовательностей с расширенным спектром с информационной частью передаваемого символа, а принимаемый сигнал поступает на первые входы блоков когерентного детектирования, на вторые входы блоков когерентного детектирования поступают эталонные последовательности с выходов формирователей эталонных последовательностей, значения с выходов блоков когерентного детектирования поступают на входы формирователей оценки правдоподобия, выходы которых подключены к входам блока комбинирования оценок правдоподобия, на выходе которого формируются оценки времени прихода и доплеровского масштабирования символа, а также значение правдоподобия этих оценок.moreover, the outputs of the extended-spectrum reference sequence shapers are connected to the inputs of the extended-spectrum reference sequence combining unit with the information part of the transmitted symbol, and the received signal is fed to the first inputs of the coherent detection blocks, the reference sequences from the outputs of the reference sequence shapers are sent to the second inputs of the coherent detection blocks, the values from the outputs of the coherent detection blocks are fed to the inputs of the teley likelihood estimator, whose outputs are connected to inputs of the combination unit likelihood estimators, the output of which generates an estimate arrival-time and the Doppler scaling the symbol and the value of the likelihood estimates.

Заявляемое устройство для оценивания времени прихода и доплеровского масштабирования символа в широкополосном канале передачи информации (фиг. 2) включает:The inventive device for estimating the time of arrival and Doppler scaling of a symbol in a broadband information channel (Fig. 2) includes:

- формирователи Опорных Последовательностей с Расширенным Спектром 211…21N в частотных поддиапазонах широкополосного канала;- Shapers of Reference Sequences with Extended Spectrum 211 ... 21N in the frequency subbands of the broadband channel;

- блок Комбинирования Опорных Последовательностей с Расширенным Спектром 220 с Информационной Частью 111…11М передаваемого Символа 100;- Combination block of Reference Sequences with Extended Spectrum 220 with Information Part 111 ... 11M of transmitted Symbol 100;

- передающий Блок 230;- transmitting unit 230;

- приемный Блок 240;- receiving unit 240;

- формирователи Эталонных Последовательностей 321…32K на приемной стороне;- shapers of Reference Sequences 321 ... 32K on the receiving side;

- блоки Когерентного Детектирования 311…31K эталонных последовательностей в принимаемом сигнале;- Coherent Detection blocks 311 ... 31K reference sequences in the received signal;

- формирователи Оценок Правдоподобия 331…33K;- Likelihood Estimators formers 331 ... 33K;

- блок Комбинирования Оценок Правдоподобия 270,- Combination of Likelihood Assessment Block 270,

причем выходы формирователей Опорных Последовательностей с Расширенным Спектром 211…21N подключены к входам блока Комбинирования Опорных Последовательностей с Расширенным Спектром 220 с Информационной Частью 111…11М передаваемого Символа 100, выход Блока Комбинирования Опорных Последовательностей с Расширенным Спектром 220 подключен к входу передающего блока 230, а принимаемый сигнал поступает на приемный Блок 240 и далее с выхода приемного Блока 240 на первые входы блоков Когерентного Детектирования 311…31K, на вторые входы блоков Когерентного Детектирования 311…31K поступают эталонные последовательности с выходов Формирователей Эталонных Последовательностей 321…32K, значения с выходов блоков Когерентного Детектирования 311…31K поступают на входы формирователей Оценок Правдоподобия 331…33K, выходы которых подключены к входам Блока Комбинирования Оценок Правдоподобия 270, на выходе которого получают оценку времени прихода

Figure 00000005
и оценку доплеровского масштабирования
Figure 00000006
принимаемого Символа 100, а также значение правдоподобия оценок
Figure 00000007
.moreover, the outputs of the Reference Sequencers with the Extended Spectrum 211 ... 21N are connected to the inputs of the Combination Block of the Reference Sequences with the Extended Spectrum 220 with the Information Part 111 ... 11M of the transmitted Symbol 100, the output of the Combination Block of the Reference Sequences with the Extended Spectrum 220 is connected to the input of the transmitting unit 230, and the received the signal enters the receiving Block 240 and then from the output of the receiving Block 240 to the first inputs of the Coherent Detection blocks 311 ... 31K, to the second inputs of the Coherent D blocks 311 ... 31K detectors receive reference sequences from the outputs of the Shapers of Reference Sequences 321 ... 32K, values from the outputs of the Coherent Detection blocks 311 ... 31K go to the inputs of the likelihood Assessment generators 331 ... 33K, the outputs of which are connected to the inputs of the Combination of the Block of Likelihood Assessment 270, estimated time of arrival
Figure 00000005
and evaluation of Doppler scaling
Figure 00000006
received Symbol 100, as well as the likelihood of estimates
Figure 00000007
.

Вариантом реализации заявляемого устройства является устройство, в котором выходы формирователя Опорных Последовательностей с Расширенным Спектром 211…21N подключены к входам Блоков Управления Мощностью Опорных Последовательностей, а выходы Блоков Управления Мощностью Опорных Последовательностей подключены к входам Блока Комбинирования Опорных Последовательностей с Расширенным Спектром 220 с Информационной Частью 111…11М передаваемого Символа 100, что компенсирует априорно известные неравномерности АЧХ гидроакустического канала и тем самым улучшает качество оценок времени прихода и коэффициента доплеровского масштабирования.An embodiment of the inventive device is a device in which the outputs of the Reference Sequence generator with the Extended Spectrum 211 ... 21N are connected to the inputs of the Power Control Units of the Reference Sequences, and the outputs of the Power Control Units of the Reference Sequences are connected to the inputs of the Combination of Reference Sequences with the Extended Spectrum 220 with Information Part 111 ... 11M of transmitted Symbol 100, which compensates for the a priori known irregularities in the frequency response of the sonar channel and those This improves the quality of estimates of the arrival time and the Doppler scaling factor.

Заявляемое устройство работает следующим образом. На передающей стороне на выходах Формирователей Опорных Последовательностей с Расширенным Спектром 211…21N в частотных поддиапазонах широкополосного канала формируются Опорные Последовательности с Расширенным Спектром 131…13N, которые затем поступают на входы Блока Комбинирования Опорных Последовательностей с Расширенным Спектром 220 с Информационной Частью 111…11М передаваемого Символа 100. Возможным методом комбинирования Опорных Последовательностей 131…13N с Информационной Частью 111…11М является метод, при котором для передачи Символа 100 используется мультиплексирование с ортогональным частотным разделением каналов, при этом часть несущих выделяется для передачи Опорных Последовательностей 131…13N, а другая часть несущих - для передачи Информационной Части 111…11М. На выходе Блока Комбинирования Опорных Последовательностей 220 с Информационной Частью формируется Символ 100, который через Передающий Блок 230 передается в широкополосный гидроакустический канал.The inventive device operates as follows. On the transmitting side, at the outputs of the Reference Sequencers of the Extended Spectrum 211 ... 21N in the frequency subbands of the broadband channel, the Reference Sequences with the Extended Spectrum 131 ... 13N are formed, which are then fed to the inputs of the Combination Block of the Reference Sequences with the Extended Spectrum 220 with the Information Part 111 ... 11M of the transmitted Symbol 100. A possible method of combining the Reference Sequences 131 ... 13N with the Information Part 111 ... 11M is the method in which to transmit the Symbol 100 orthogonal frequency division multiplexing is used, with part of the carriers allocated for transmitting the Reference Sequences 131 ... 13N, and the other part of the carriers for transmitting the Information Part 111 ... 11M. At the output of the Block of Combination of Reference Sequences 220 with the Information Part, Symbol 100 is formed, which is transmitted through the Transmitting Unit 230 to a broadband sonar channel.

При распространении от передающей части к приемной части Символ 100 в широкополосном гидроакустическом канале подвергается, помимо прочего, задержке, которая обусловлена конечной скоростью распространения акустических волн в водной среде, а также доплеровскому масштабированию, которое обусловлено влиянием эффекта Доплера на каждую частотную компоненту Символа 100. Влияние доплеровского масштабирования характеризуется коэффициентом доплеровского масштабирования (1+α). При приеме каждая частотная компонента, переданная с частотой ƒ, будет преобразована в (1+α)ƒ, а сигнал, переданный с длительностью T, будет преобразован в сигнал длительностью (1-α)Т. При работе на мелководье при высокой подвижности объектов, между которыми производится обмен информационными символами, гидроакустические каналы характеризуются малым временем когерентности, поэтому для успешной демодуляции и декодирования Символа 100 на приемной стороне необходимо надежно оценивать как задержку, так и коэффициент доплеровского масштабирования для каждого Символа 100.When propagating from the transmitting part to the receiving part, Symbol 100 in the broadband sonar channel undergoes, among other things, a delay due to the finite propagation velocity of acoustic waves in the aquatic environment, as well as to Doppler scaling due to the influence of the Doppler effect on each frequency component of Symbol 100. Effect Doppler scaling is characterized by a Doppler scaling factor (1 + α). Upon reception, each frequency component transmitted with a frequency ƒ will be converted to (1 + α) ƒ, and a signal transmitted with a duration T will be converted to a signal of (1-α) T duration. When working in shallow water with high mobility of objects between which information symbols are exchanged, sonar channels are characterized by a short coherence time, therefore, for successful demodulation and decoding of Symbol 100 on the receiving side, it is necessary to reliably estimate both the delay and the Doppler scaling coefficient for each Symbol 100.

На приемной стороне принимаемый из широкополосного гидроакустического канала сигнал поступает в Приемный Блок 240 и далее с выхода Приемного Блока 240 поступает на первые входы Блоков Когерентного Детектирования 311…31K (фиг. 3), при этом на вторые входы Блоков Когерентного Детектирования 311…31K поступают эталонные последовательности с выходов Формирователей Эталонных Последовательностей 321..32K. На выходах Формирователей Эталонных Последовательностей 321…32K формируют набор эталонных последовательностей для каждой из Опорных Последовательностей с Расширенным Спектром 131…13N. Каждая из эталонных последовательностей в наборе может представлять собой соответствующую Опорную Последовательность с Расширенным Спектром, формируемую в передатчике, и подвергнутую доплеровскому масштабированию. Для Опорной Последовательности с Расширенным Спектром с номером n=1…N количество Kn эталонных последовательностей Unk), k=1…Kn в наборе может определяться необходимой разрешающей способностью по коэффициенту доплеровского масштабирования Δnk+lk и ожидаемым диапазоном значений коэффициента доплеровского масштабирования (1+αmin)≤(1+αk)≤(1+αmax).On the receiving side, the signal received from the broadband sonar channel enters the Receiving Unit 240 and then from the output of the Receiving Block 240 it enters the first inputs of the Coherent Detection Units 311 ... 31K (Fig. 3), while the reference inputs of the Coherent Detection Units 311 ... 31K receive the reference sequences from the outputs of the Shapers of Reference Sequences 321..32K. At the outputs of the Shapers of Reference Sequences 321 ... 32K, a set of reference sequences is formed for each of the Reference Sequences with an Extended Range of 131 ... 13N. Each of the reference sequences in the set may be a corresponding Spread Spectrum Reference Pattern formed in the transmitter and subjected to Doppler scaling. For the Reference Sequence with the Extended Spectrum with the number n = 1 ... N, the number K n of reference sequences U nk ), k = 1 ... K n in the set can be determined by the necessary resolution by the Doppler scaling coefficient Δ n = α k + l - α k and the expected range of values of the Doppler scaling coefficient (1 + α min ) ≤ (1 + α k ) ≤ (1 + α max ).

При приеме, для исключения влияния неравномерности АЧХ канала частотные составляющие, находящиеся вне рабочего частотного диапазона некоторого Блока Когерентного Детектирования 311…31K, могут значительно превышать уровень сигнала внутри рабочего частотного диапазона Блока Когерентного Детектирования 311…31K, т.е. представлять для некоторого Блока Когерентного Детектирования 311…31K внеполосную помеху, воздействие которой снижает точность оценки времени прихода и коэффициента доплеровского масштабирования в устройстве устанавливают фильтры, через которые принимаемый сигнал поступает на первые входы Блоков Когерентного Детектирования 311…31K. Полоса каждого Фильтра 251…25N выбрана так, что фильтр пропускает соответствующую принимаемую Опорную Последовательность с Расширенным Спектром 131…13N с учетом возможных изменений ее центральной частоты и ширины полосы, вызванных доплеровским масштабированием, и подавляет частотные составляющие принимаемого Символа 100, находящиеся вне рабочего частотного диапазона соответствующего Блока Когерентного Детектирования 311…31K. Это повышает точность оценки времени прихода и коэффициента доплеровского масштабирования для информационных символов в системе передачи данных в широкополосных гидроакустических каналах с высокой частотной избирательностью. Наличие Защитных Частотных Интервалов 121…12N и 141…14N может ослабить требования к ширине переходной полосы АЧХ Фильтров 211…21N.When receiving, in order to eliminate the influence of the channel frequency response unevenness, frequency components that are outside the working frequency range of a Coherent Detection Unit 311 ... 31K can significantly exceed the signal level inside the working frequency range of a Coherent Detection Unit 311 ... 31K, i.e. to present out-of-band interference for some Coherent Detection Unit 311 ... 31K, the effect of which reduces the accuracy of estimating the arrival time and Doppler scaling coefficient, the filters are installed in the device through which the received signal is fed to the first inputs of the Coherent Detection Unit 311 ... 31K. The band of each Filter 251 ... 25N is selected so that the filter passes the corresponding received Reference Sequence with the Extended Spectrum 131 ... 13N taking into account possible changes in its center frequency and bandwidth caused by Doppler scaling, and suppresses the frequency components of the received Symbol 100 that are outside the operating frequency range corresponding Coherent Detection Unit 311 ... 31K. This increases the accuracy of estimating the arrival time and Doppler scaling factor for information symbols in a data transmission system in broadband sonar channels with high frequency selectivity. The presence of Protective Frequency Intervals 121 ... 12N and 141 ... 14N can weaken the requirements for the transition bandwidth of the frequency response of the Filters 211 ... 21N.

Блок Когерентного Детектирования 311…31K в моменты времени τ может вычислять значения взаимно-корреляционной функции Rnk(τ) принимаемого сигнала и k-й эталонной последовательности Unk), сформированной для n-й Опорной Последовательности с Расширенным Спектром 131…13N. Таким образом для одного Символа 100, включающего N Опорных Последовательностей с Расширенным Спектром 131…13N, в Блоках Когерентного Детектирования 311…31K вычисляется

Figure 00000009
значений Rnk(τ). Когерентное детектирование эталонных последовательностей в принимаемом сигнале может производиться как во временной, так и в частотной областях.The Coherent Detection Unit 311 ... 31K at times τ can calculate the values of the cross-correlation function R nk (τ) of the received signal and the k-th reference sequence U nk ) generated for the n-th Reference Sequence with Extended Spectrum 131 ... 13N . Thus, for one Symbol 100, including N Reference Sequences with the Extended Range 131 ... 13N, in Coherent Detection Units 311 ... 31K
Figure 00000009
values of R nk (τ). Coherent detection of reference sequences in the received signal can be performed both in the time and in the frequency domains.

Значения с выходов Блоков Когерентного Детектирования 311…31K поступают на входы Формирователей Оценок Правдоподобия 331…33K, которые вычисляют оценки правдоподобия Qnk(τ). Оценкой правдоподобия может являться значение

Figure 00000010
, где Еk - энергия k-й эталонной последовательности Unk), которая может быть вычислена заранее,
Figure 00000011
- оценка энергии n-й принятой опорной последовательности в моменты времени τ с учетом влияния коэффициента доплеровского масштабирования (1+αk), т.е. с учетом изменения длительности. Такая оценка правдоподобия Qnk(τ) может принимать значения в диапазоне [0, 1]. Большее значение Qnk(τ) показывает, что n-я опорная последовательность в принятом сигнале с большей достоверностью соответствует эталонной последовательности Unk).The values from the outputs of the Coherent Detection Units 311 ... 31K are fed to the inputs of the Likelihood Estimators 331 ... 33K, which calculate the likelihood estimates Q nk (τ). The likelihood score may be
Figure 00000010
where E k is the energy of the kth reference sequence U nk ), which can be calculated in advance,
Figure 00000011
- estimation of the energy of the nth received reference sequence at time instants τ taking into account the influence of the Doppler scaling coefficient (1 + α k ), i.e. subject to changes in duration. Such a likelihood estimate Q nk (τ) can take values in the range [0, 1]. A larger value of Q nk (τ) indicates that the nth reference sequence in the received signal corresponds with greater certainty to the reference sequence U nk ).

С выходов Формирователей Оценок Правдоподобия 331…33K значения поступают на вход Блока Комбинирования Оценок Правдоподобия 270, где происходит комбинирование оценок правдоподобия и формирования оценки времени прихода

Figure 00000005
и оценки коэффициента доплеровского масштабирования
Figure 00000006
, а также значения правдоподобия оценок
Figure 00000007
для принятого информационного символа.From the outputs of the Likelihood Estimators, 331… 33K, the values are input to the Block of Combining Likelihood Estimates 270, where the likelihood estimates are combined with the formation of the estimate of the time of arrival
Figure 00000005
and estimates of the Doppler scaling factor
Figure 00000006
and likelihood estimates
Figure 00000007
for the received information symbol.

Список цитируемых документовList of cited documents

US Pat. 7859944 Apparatus, systems and methods for enhanced multi-carrier based underwater acoustic communications.US Pat. 7859944 Apparatus, systems and methods for enhanced multi-carrier based underwater acoustic communications.

US Pat. 8023596 Doppler tracking method and device for a wide band modem. US Pat. 8023596 Doppler tracking method and device for a wide band modem.

Nadav Levanon, Eli Mozeson; Radar Signals, ISBN: 978-0-471-47378-7, Jul 2004, Wiley-IEEE Press, 432 pages.Nadav Levanon, Eli Mozeson; Radar Signals, ISBN: 978-0-471-47378-7, Jul 2004, Wiley-IEEE Press, 432 pages.

Claims (25)

1. Способ оценки параметров канала в широкополосной гидроакустической связи, состоящий в определении времени прихода и величины доплеровского масштабирования символа, включающий при передаче формирование опорных последовательностей с расширенным спектром в частотных поддиапазонах широкополосного канала, комбинирование опорных последовательностей с расширенным спектром с информационной частью символа, комбинация которых поступает на вход передающего блока, а при приеме - для опорной последовательности с расширенным спектром, принятой с выхода приемного блока, формирование эталонных последовательностей для ожидаемого диапазона доплеровского масштабирования, затем когерентное детектирование эталонных последовательностей в принятом сигнале, а в случае, когда для опорной последовательности с расширенным спектром формируют эталонный ряд отсчетов, производят формирование исходного ряда отсчетов для принимаемого сигнала, формирование скорректированных рядов отсчетов из исходного ряда отсчетов для ожидаемого диапазона доплеровского масштабирования принимаемых опорных последовательностей с расширенным спектром, после этого осуществляют когерентное детектирование эталонных рядов отсчетов в каждом из скорректированных рядов отсчетов, затем, как и в случае после когерентного детектирования эталонных последовательностей в принятом сигнале, производят формирование оценок правдоподобия, комбинирование оценок правдоподобия и определяют оценку времени прихода и величины доплеровского масштабирования символа, а также значений правдоподобия оценок.1. A method for estimating channel parameters in broadband sonar communication, which consists in determining the arrival time and magnitude of the Doppler scaling of a symbol, including the transmission of the formation of reference sequences with extended spectrum in the frequency subbands of the broadband channel, the combination of reference sequences with extended spectrum and the information part of the symbol, a combination of which arrives at the input of the transmitting unit, and upon reception, for a reference sequence with an extended spectrum, from the output of the receiving unit, the formation of reference sequences for the expected range of Doppler scaling, then coherent detection of the reference sequences in the received signal, and in the case when a reference series of samples is formed for the reference sequence with an extended spectrum, the initial series of samples for the received signal are generated, formation the adjusted series of samples from the original series of samples for the expected range of Doppler scaling are accepted reference sequences with an extended spectrum, then coherent detection of reference series of samples in each of the adjusted series of samples is carried out, then, as in the case after coherent detection of reference sequences in a received signal, likelihood estimates are generated, likelihood estimates are combined and the estimated time of arrival is determined and the magnitude of the Doppler scaling of the symbol, as well as the likelihood of estimates. 2. Способ по п. 1, характеризующийся тем, что формирование опорных последовательностей с расширенным спектром в частотных поддиапазонах широкополосного канала включает формирование сигнала с угловой модуляцией.2. The method according to p. 1, characterized in that the formation of reference sequences with spread spectrum in the frequency subbands of the broadband channel includes generating a signal with angular modulation. 3. Способ по п. 1, характеризующийся тем, что формирование опорных последовательностей с расширенным спектром в частотных поддиапазонах широкополосного канала включает формирование псевдослучайной последовательности.3. The method according to p. 1, characterized in that the formation of reference sequences with spread spectrum in the frequency subbands of the broadband channel includes forming a pseudo-random sequence. 4. Способ по п. 1, характеризующийся тем, что формирование опорных последовательностей с расширенным спектром в частотных поддиапазонах широкополосного канала включает формирование модулированных несущих.4. The method according to p. 1, characterized in that the formation of reference sequences with spread spectrum in the frequency subbands of the broadband channel includes the formation of modulated carriers. 5. Способ по п. 1, характеризующийся тем, что формирование опорных последовательностей с расширенным спектром в частотных поддиапазонах широкополосного канала включает циклический сдвиг опорной последовательности.5. The method according to p. 1, characterized in that the formation of reference sequences with spread spectrum in the frequency subbands of the broadband channel includes a cyclic shift of the reference sequence. 6. Способ по п. 1, характеризующийся тем, что формирование опорных последовательностей с расширенным спектром в частотных поддиапазонах широкополосного канала включает повторение части опорной последовательности.6. The method according to p. 1, characterized in that the formation of reference sequences with spread spectrum in the frequency subbands of the broadband channel includes repeating part of the reference sequence. 7. Способ по п. 1, характеризующийся тем, что количество опорных последовательностей с расширенным спектром устанавливают в зависимости от амплитудно-частотной характеристики канала передачи информации.7. The method according to p. 1, characterized in that the number of reference sequences with an extended spectrum is set depending on the amplitude-frequency characteristics of the information transmission channel. 8. Способ по п. 1, характеризующийся тем, что центральную частоту опорной последовательности с расширенным спектром устанавливают в зависимости от амплитудно-частотной характеристики канала передачи информации.8. The method according to p. 1, characterized in that the center frequency of the reference sequence with an extended spectrum is set depending on the amplitude-frequency characteristics of the information transmission channel. 9. Способ по п. 1, характеризующийся тем, что спектральную плотность мощности опорной последовательности с расширенным спектром устанавливают в зависимости от амплитудно-частотной характеристики канала передачи информации.9. The method according to p. 1, characterized in that the power spectral density of the reference sequence with an extended spectrum is set depending on the amplitude-frequency characteristics of the information transmission channel. 10. Способ по п. 1, характеризующийся тем, что комбинирование опорных последовательностей с расширенным спектром с информационной частью передаваемого символа производят когерентно.10. The method according to p. 1, characterized in that the combination of reference sequences with an extended spectrum with the information part of the transmitted symbol is produced coherently. 11. Способ по п. 1, характеризующийся тем, что комбинирование опорных последовательностей с расширенным спектром с информационной частью передаваемого символа включает формирование защитного частотного интервала между опорной последовательностью и информационной частью символа.11. The method according to p. 1, characterized in that the combination of reference sequences with spread spectrum with the information part of the transmitted symbol includes forming a protective frequency interval between the reference sequence and the information part of the symbol. 12. Способ по п. 1, характеризующийся тем, что набор эталонных последовательностей формируют на приемной стороне с учетом взаимной радиальной скорости передающей и приемной антенн и скорости движения среды в радиальном направлении.12. The method according to p. 1, characterized in that a set of reference sequences is formed on the receiving side, taking into account the mutual radial speed of the transmitting and receiving antennas and the speed of the medium in the radial direction. 13. Способ по п. 1, характеризующийся тем, что когерентное детектирование эталонной последовательности в принимаемом сигнале производят с учетом параметров канала передачи информации.13. The method according to p. 1, characterized in that the coherent detection of the reference sequence in the received signal is performed taking into account the parameters of the information transmission channel. 14. Способ по п. 1, характеризующийся тем, что когерентное детектирование эталонной последовательности в принимаемом сигнале включает вычисление преобразования Фурье.14. The method according to p. 1, characterized in that the coherent detection of the reference sequence in the received signal includes the calculation of the Fourier transform. 15. Способ по п. 1, характеризующийся тем, что формирование скорректированных рядов отсчетов из исходного ряда отсчетов включает полиномиальную интерполяцию отсчетов из исходного ряда отсчетов.15. The method according to p. 1, characterized in that the formation of the corrected series of samples from the original series of samples includes polynomial interpolation of samples from the original series of samples. 16. Способ по п. 1, характеризующийся тем, что формирование оценки правдоподобия включает получение оценки отношения сигнал-шум.16. The method according to p. 1, characterized in that the formation of a likelihood assessment includes obtaining an estimate of the signal-to-noise ratio. 17. Способ по п. 1, характеризующийся тем, что формирование оценки правдоподобия включает вычисление энергетического параметра опорной последовательности.17. The method according to p. 1, characterized in that the formation of a likelihood assessment includes calculating the energy parameter of the reference sequence. 18. Способ по п. 1, характеризующийся тем, что сформированную оценку правдоподобия сравнивают с порогом оценки правдоподобия.18. The method according to p. 1, characterized in that the generated likelihood assessment is compared with the likelihood assessment threshold. 19. Способ по п. 1, характеризующийся тем, что при комбинировании оценок правдоподобия вычисляют взвешенное среднее значение оценок правдоподобия.19. The method according to p. 1, characterized in that when combining the likelihood estimates calculate the weighted average value of the likelihood estimates. 20. Способ по п. 1, характеризующийся тем, что при комбинировании оценок правдоподобия определяют максимум оценок правдоподобия.20. The method according to p. 1, characterized in that when combining the likelihood estimates determine the maximum likelihood estimates. 21. Способ по п. 1, характеризующийся тем, что при комбинировании оценок правдоподобия первоначально выполняют комбинирование оценок правдоподобия в каждом поддиапазоне широкополосного канала передачи информации.21. The method according to claim 1, characterized in that when combining the likelihood estimates, the combination of likelihood estimates is initially performed in each subband of the broadband information channel. 22. Способ по п. 1, характеризующийся тем, что при комбинировании оценок правдоподобия выполняют интерполяцию результатов.22. The method according to claim 1, characterized in that when combining the likelihood estimates, the results are interpolated. 23. Устройство для реализации способа оценки параметров канала при широкополосной гидроакустической связи, включающее: формирователи опорных последовательностей с расширенным спектром в частотных поддиапазонах широкополосного канала; блок комбинирования опорных последовательностей с расширенным спектром с информационной частью передаваемого символа; передающий блок; приемный блок; формирователи эталонных последовательностей на приемной стороне; блоки когерентного детектирования эталонных последовательностей в принимаемом сигнале; формирователи оценок правдоподобия; блок комбинирования оценок правдоподобия, при этом в передающей части - выходы формирователей опорных последовательностей с расширенным спектром подключены к входам блока комбинирования опорных последовательностей с расширенным спектром с информационной частью передаваемого символа, выход которого подключен к передающему блоку, а в приемной части - выход приемного блока и выходы формирователей эталонных последовательностей подключены к входам блоков когерентного детектирования эталонных последовательностей, выходы которых подключены к входам формирователей оценок правдоподобия, выходы которых подключены к входам блока комбинирования оценок правдоподобия, на выходе которого формируются оценки времени прихода и доплеровского масштабирования принимаемого символа, а также значения правдоподобия оценок.23. A device for implementing a method for estimating channel parameters for broadband sonar communication, including: extended-range reference sequence formers in the frequency subbands of the broadband channel; block combining reference sequences with an extended spectrum with the information part of the transmitted symbol; transmitting unit; receiving unit; shapers of reference sequences on the receiving side; blocks of coherent detection of reference sequences in the received signal; likelihood estimators; a combination block of likelihood estimates, while in the transmitting part, the outputs of the extended-sequence reference sequence formers are connected to the inputs of the extended-sequence reference sequence combining unit with the information part of the transmitted symbol, the output of which is connected to the transmitting unit, and in the receiving part, the output of the receiving unit and the outputs of the reference sequence generators are connected to the inputs of the coherent detection blocks of the reference sequences, the outputs of which are connected cheny to the input of the likelihood estimates whose outputs are connected to inputs of the combination unit likelihood estimators, the output of which generates an estimate arrival-time and the Doppler scaling the received symbols, as well as estimates of the likelihood values. 24. Устройство по п. 23, характеризующееся тем, что выходы формирователей опорных последовательностей с расширенным спектром подключены к входам блоков управления мощностью опорных последовательностей с расширенным спектром, а выходы блоков управления мощностью опорных последовательностей с расширенным спектром подключены к входам блока комбинирования опорных последовательностей с расширенным спектром с информационной частью передаваемого символа.24. The device according to claim 23, characterized in that the outputs of the extended-range reference sequence formers are connected to the inputs of the extended-spectrum reference power control units, and the outputs of the extended-spectrum reference sequence power control units are connected to the inputs of the extended reference sequence combining unit spectrum with the information part of the transmitted symbol. 25. Устройство по п. 23, характеризующееся тем, что в приемной части перед блоками когерентного детектирования установлены фильтры.25. The device according to p. 23, characterized in that in the receiving part in front of the coherent detection units, filters are installed.
RU2018139703A 2018-11-08 2018-11-08 Method of estimating channel parameters in broadband hydroacoustic communication and a device for realizing said channel RU2700005C1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2018139703A RU2700005C1 (en) 2018-11-08 2018-11-08 Method of estimating channel parameters in broadband hydroacoustic communication and a device for realizing said channel

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2018139703A RU2700005C1 (en) 2018-11-08 2018-11-08 Method of estimating channel parameters in broadband hydroacoustic communication and a device for realizing said channel

Publications (1)

Publication Number Publication Date
RU2700005C1 true RU2700005C1 (en) 2019-09-12

Family

ID=67989871

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2018139703A RU2700005C1 (en) 2018-11-08 2018-11-08 Method of estimating channel parameters in broadband hydroacoustic communication and a device for realizing said channel

Country Status (1)

Country Link
RU (1) RU2700005C1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111366897A (en) * 2020-04-02 2020-07-03 中国科学院深海科学与工程研究所 Signal design and processing method for high-precision underwater sound positioning system
CN115987742A (en) * 2023-03-20 2023-04-18 四川创智联恒科技有限公司 Frequency domain Doppler frequency offset resistant extension method based on OFDM communication system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101505291A (en) * 2009-03-16 2009-08-12 哈尔滨工程大学 Hydroacoustic communication differential decoding method based on OFDM coding
CN103095639A (en) * 2013-01-15 2013-05-08 哈尔滨工程大学 Orthogonal frequency division multiplexing (OFDM) underwater acoustic communication parallel iterative inter-carrier interference (ICI) elimination method
CN103618686A (en) * 2013-11-22 2014-03-05 江苏科技大学 Method for accurately estimating underwater sound OFDM Doppler factor
CN107426120A (en) * 2017-07-24 2017-12-01 哈尔滨工程大学 A kind of underwater sound OFDM MFSK channel equalization methods based on least mean-square error
RU183781U1 (en) * 2018-06-14 2018-10-02 Публичное акционерное общество "Российский институт мощного радиостроения", Патентное бюро DEVICE FOR DOPPLER FREQUENCY DEFINITION DETERMINATION DETERMINATION BY THE INFORMATION PHASOMANIPULATED SIGNAL BY APPROXIMATION OF PHASE DEFLECTION

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101505291A (en) * 2009-03-16 2009-08-12 哈尔滨工程大学 Hydroacoustic communication differential decoding method based on OFDM coding
CN103095639A (en) * 2013-01-15 2013-05-08 哈尔滨工程大学 Orthogonal frequency division multiplexing (OFDM) underwater acoustic communication parallel iterative inter-carrier interference (ICI) elimination method
CN103618686A (en) * 2013-11-22 2014-03-05 江苏科技大学 Method for accurately estimating underwater sound OFDM Doppler factor
CN107426120A (en) * 2017-07-24 2017-12-01 哈尔滨工程大学 A kind of underwater sound OFDM MFSK channel equalization methods based on least mean-square error
RU183781U1 (en) * 2018-06-14 2018-10-02 Публичное акционерное общество "Российский институт мощного радиостроения", Патентное бюро DEVICE FOR DOPPLER FREQUENCY DEFINITION DETERMINATION DETERMINATION BY THE INFORMATION PHASOMANIPULATED SIGNAL BY APPROXIMATION OF PHASE DEFLECTION

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111366897A (en) * 2020-04-02 2020-07-03 中国科学院深海科学与工程研究所 Signal design and processing method for high-precision underwater sound positioning system
CN111366897B (en) * 2020-04-02 2023-03-21 中国科学院深海科学与工程研究所 Signal design and processing method for high-precision underwater sound positioning system
CN115987742A (en) * 2023-03-20 2023-04-18 四川创智联恒科技有限公司 Frequency domain Doppler frequency offset resistant extension method based on OFDM communication system
CN115987742B (en) * 2023-03-20 2023-06-16 四川创智联恒科技有限公司 OFDM communication system-based frequency domain Doppler frequency offset expansion resisting method

Similar Documents

Publication Publication Date Title
Zheng et al. Super-resolution delay-Doppler estimation for OFDM passive radar
JP4551714B2 (en) Orthogonal frequency division multiplexed signal processing method
Baczyk et al. Reconstruction of the reference signal in DVB-T-based passive radar
CN101005475A (en) Method and system for synchronizing time and frequency in orthogonal frequency division multiplex communication
Perrine et al. Doppler estimation and correction for shallow underwater acoustic communications
WO2015135295A1 (en) Method, device and system for detecting random access signal in interference environment
CN110736977A (en) integration signal generation method suitable for sonar detection
RU2700005C1 (en) Method of estimating channel parameters in broadband hydroacoustic communication and a device for realizing said channel
Aparicio et al. Asynchronous detection and identification of multiple users by multi-carrier modulated complementary set of sequences
JP2012088279A (en) Radar device and mobile target detecting method to be applied for radar device
Chi et al. Utilizing orthogonal coprime signals for improving broadband acoustic Doppler current profilers
Li et al. Cross power spectral density based beamforming for underwater acoustic communications
Baldone et al. Doppler estimation and correction for JANUS underwater communications
Dokhanchi et al. A half-duplex joint communications and sensing system using ZP-OFDM
US7003415B2 (en) Delay profile estimation apparatus and a correlating unit
Kochańska et al. Underwater acoustic communications in time-varying dispersive channels
US8385438B1 (en) System and method for adaptive synchronization
CN111539323B (en) Frequency estimation method and device for cyclic prefix linear frequency modulation signal
CN110191079B (en) Non-coherent combined capturing method and device
Babich et al. A simple method for TOA estimation in OFDM systems
Socheleau Cyclostationarity of communication signals in underwater acoustic channels
CN115333912B (en) Communication interference integrated signal design and processing method based on parameter modulation
Rafique et al. A novel frame design for integrated communication and sensing based on position modulation
WO2019142150A1 (en) Method for vehicle location estimation using orthogonal frequency-division multiplexing
Unal et al. Field experiments with doppler compensation in High-Frequency underwater acoustic communication system