Например, Бобцов

ИМИТАТОР СИГНАЛОВ НА ВЫХОДЕ ПРИЕМНЫХ ЭЛЕМЕНТОВ ПАССИВНОЙ ГИДРОАКУСТИЧЕСКОЙ СТАНЦИИ С ГИБКОЙ ПРОТЯЖЕННОЙ БУКСИРУЕМОЙ АНТЕННОЙ

ИМИТАТОР СИГНАЛОВ НА ВЫХОДЕ ПРИЕМНЫХ ЭЛЕМЕНТОВ ПАССИВНОЙ …
УДК 681.88
ИМИТАТОР СИГНАЛОВ НА ВЫХОДЕ ПРИЕМНЫХ ЭЛЕМЕНТОВ ПАССИВНОЙ ГИДРОАКУСТИЧЕСКОЙ СТАНЦИИ С ГИБКОЙ ПРОТЯЖЕННОЙ БУКСИРУЕМОЙ АНТЕННОЙ
Ю.В. Шафранюк
Предложена математическая модель имитатора сигналов на выходе приемных элементов гибкой протяженной буксируемой антенны, учитывающая особенности излучения, распространения и приема сигналов и помех в реальных морских условиях. Приводятся требования к имитатору. Предлагается общая схема и принцип построения. На примере имитационных данных иллюстрируется его работоспособность. Ключевые слова: гидроакустика, пассивные гидроакустические станции, антенны, многолучевое распространение, имитация.
Введение Пассивные гидроакустические станции с гибкими протяженными буксируемыми антеннами (ГАС с ГПБА) являются одними из основных информационных систем освещения подводной обстановки и позволяют на достаточно больших дальностях осуществлять обнаружение подводных объектов [1]. Наиболее сложными задачами, решаемыми в ГАС с ГПБА, являются обнаружение шумящих объектов, их классификация и определение координат и параметров движения. Создание эффективных алгоритмов, решающих эти задачи, представляет большой научный и практический интерес. Поскольку такие алгоритмы отличаются высокой сложностью, одним из основных этапов их проектирования является их отработка на модельных данных, учитывающих особенности излучения, распространения и приема сигналов и помех в реальных морских условиях. Целью работы являлась разработка имитатора, обеспечивающего моделирование (имитацию) смеси сигналов и помех на выходе приемных элементов ГПБА в интересах отработки алгоритмов классификации и определения координат и параметров движения обнаруженных объектов. Научная новизна разработанного имитатора заключается в учете всех значимых факторов, влияющих на эффективность функционирования ГАС с ГПБА, в первую очередь, особенностей распространения сигнала цели в реальной морской среде.
Требования к имитатору Для обеспечения адекватной отработки алгоритмов обнаружения, классификации и определения координат и параметров движения шумящих объектов в реальных морских условиях имитатор должен учитывать:  конструкцию ГПБА (ее длину, количество приемных элементов, расстояние между ними, рабочий диапазон частот) [1];  особенности маневрирования и шумоизлучения морских объектов (скорость, курс, глубина, уровень и спектр шумоизлучения в рабочем диапазоне частот в функции скорости объекта, параметры дискрет-
64 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)

Ю.В. Шафранюк
ных составляющих в спектре шумоизлучения, обусловленные работой корабельных механизмов, параметры амплитудной вально-лопастной модуляции шума цели, а также модуляции вследствие качки на волнении) [2–4];  особенности распространения сигнала цели в реальной морской среде (пространственное затухание, многолучевость и аномалию распространения) [5];  четырехкомпонентную модель распределенной помехи (шумы моря, шумы носителя, гидродинамические шумы обтекания ГПБА, шумы удаленного судоходства) [3, 5];  особенности преобразования сигнала в приемных каналах антенны [6, 7].
Структура имитатора
С учетом перечисленных требований была предложена структура имитатора, изображенная в виде блок-схемы на рис. 1. Ниже приводится описание блоков имитатора.

Рис. 1. Структура имитатора
Блок 1. Имитатор внешней обстановки. Блок позволяет задать произвольный тестовый эпизод, описываемый:
 конструктивными параметрами ГПБА (см. выше), спектром гидродинамической помехи в функции скорости буксировки;
 начальным расположением носителя ГАС с ГПБА и произвольного числа объектов произвольных классов и их движением по произвольным траекториям;
 спектром сигнала каждого моделируемого объекта (включая носитель ГАС с ГПБА) и спектром его амплитудной огибающей;
 гидроакустическими условиями, характеризуемыми глубиной района, вертикальным распределением скорости звука, волнением поверхности моря, частотной и угловой зависимостью коэффициента отражения акустического сигнала от дна, частотной зависимостью пространственного затухания сигнала, спектром шумов моря, спектром шумов удаленного судоходства. Имитатор внешней обстановки с заданным шагом по времени (как правило,
1 с) выдает текущее относительное (относительно носителя ГАС с ГПБА) положение всех объектов, их параметры движения. Он также синхронизирует работу всех остальных блоков.
Блок 2. Гидроакустические расчеты. В блоке на моменты времени, генерируемые в блоке 1, рассчитываются параметры сигнала каждого объекта и всех видов помех на входе ГПБА. Расчет осуществляется по специальной программе гидроакустических расчетов, учитывающей распространение акустических сигналов и помех в реальном морском волноводе [5].
Блок 3. Формирование распределенных помех. В блоке производится расчет воздействия на приемные элементы антенны помеховой составляющей акустического поля и формирование распределенных помех на выходе антенны.
Блок 4. Формирование сигналов локальных источников. В блоке производится вычисление параметров сигнала от локальных источников на входе антенны и формирование сигналов от локальных источников, приведенных к выходу каждого приемного элемента ГПБА (с учетом временного сдвига, зависящего от направления прихода сигнала на антенну).
Полученный сигнал от локального источника суммируется с сигналами от других источников и помехой на выходе текущего приемного элемента, затем подвергается дискретизации и усиливается, со-

Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)

65

ИМИТАТОР СИГНАЛОВ НА ВЫХОДЕ ПРИЕМНЫХ ЭЛЕМЕНТОВ ПАССИВНОЙ …

гласно амплитудно-частотной характеристике усиления приемного тракта антенны. В итоге формируется цифровой сигнал с выхода элементарного приемного канала во временной области. Все программное обеспечение имитатора написано на языке программирования С++ в стандарте ANSI.
Ниже приведены математические модели, предложенные автором, для реализации блока 3 и блока 4 имитатора.

Метод формирования распределенных помех

Особенность моделирования сигналов локальных источников и распределенных помех состоит в том, что сигналы локальных источников можно считать плосковолновыми и, как следствие, полностью коррелированными на входе приемных элементов ГПБА (с учетом временного сдвига, зависящего от направления прихода сигнала на антенну), а коэффициент корреляции распределенных помех зависит от
частоты f , расстояния между элементами антенны d и скорости звука в воде с по формуле

sin 2fd

rрп 

c 2fd

.

c
В связи с этим предлагается следующая математическая модель, описывающая шумовую помеху в частотной области:

 S(if )  norm 0,1, iSR  1

   Scltr (if

, iSR )



 S (if 

 )  norm 
 Scltr

rpn i f , iSR

Scltr

i f , iSR 1 ,

1  rp2n

 , 

iSR

1,

   S(if ) 



Psmo

 

f

if f0

2  





Pngb

 

f

if f0

4 , 

 где if  1, N f – индекс частотного отсчета; f if – рабочие частоты ГПБА, Гц; Psmo – приведенный к

частоте 1 кГц уровень шумов моря, мПа·мПа; Pngb – приведенный к частоте 1 кГц уровень гидродинами-
ческой помехи, мПа·мПа; f0 – частота приведения, Гц; iSR – индекс приемного элемента ГПБА;
norm(,) – функция плотности распределения вероятности для случая гауссового распределения (пер-

вый аргумент – математическое ожидание, второй – среднеквадратическое отклонение). Для получения шумовой помехи во временной области используется обратное преобразование Фурье.

Метод формирования сигналов от локальных источников

Сигнал от каждого локального источника, приведенный к входу антенны, моделируется по следующей формуле [3]:

Nray

            u(ti

)



iray

1

 

Pnoise

iray

unoise

ti  t

iray

1  uM ti  t iray

 uDS

ti  t

iray

, iray

 

,

где ti – временной отсчет; u(ti ) – ti -й отсчет сигнала; iray – индекс луча; Nray – число лучей;

  unoise ti  t iray – ti -й отсчет широкополосной составляющей сигнала объекта с заданным наклоном

спектра, моделируемый датчиком нормально распределенных случайных чисел с нулевым математиче-

 ским ожиданием и единичным среднеквадратическим отклонением; Pnoise iray – уровень широкополос-

 ной составляющей сигнала объекта, Па; t iray – временная разница в набеге iray -го луча, с;

      uM ti  t iray – модулирующая часть сигнала; uDS ti  t iray ,iray – полигармоническая часть сигнала

(дискретная часть спектра). Модулирующая составляющая сигнала моделируется по формуле [4]

NM    uM ti   kM iM  cos 2fM iM ti  M  , M  norm 0, 2fM iM  fM iM  ,
iM 1
где iM – индекс гармоники; NM – число гармоник; kM iM  – коэффициент модуляции для iM -й гармо-

ники модулирующего сигнала; fM iM  – частота iM -й гармоники модулирующего сигнала, Гц; M –

66 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)

Ю.В. Шафранюк

случайная фаза iM -й гармоники; fM iM  – ширина iM -й гармоники. Полигармоническая составляющая
сигнала моделируется по формуле

N DS

      uDS ti ,iray 

 P iDS , iray cos 2fDS iDS ti  DS  ,

iDS 1

     PDS iDS ,iray

 PD*S

iDS

100,05R6 A iray R

,

 DS  norm 0, 2fDS iDS   fDS iDS ,

 где iDS – индекс гармоники; NDS – число гармоник; P iDS ,iray – уровень iDS -й гармоники полигармониче-

ского сигнала, приведенный ко входу антенны, Па;  PD*S iDS – уровень iDS -й гармоники полигармониче-

 ской составляющей сигнала в точке излучения, Па; A iray – аномалия распространения для iray -го луча от

 объекта; R – расстояние до объекта в горизонтальной плоскости, км;   a  b favg f0 c – коэффициент

пространственного затухания; a , b , c – эмпирические коэффициенты модели пространственного затуха-

ния, взятые для среднегеометрической частоты диапазона спектрального анализа – favg  fhigh  flow , где

flow – нижняя частота спектрального диапазона, а fhigh – верхняя; fDS iDS  – частота iDS -й гармоники по-

лигармонического сигнала, Гц; DS – случайная фаза iDS -й гармоники; где fDS iDS  – ширина iDS -й

гармоники. Уровень широкополосной составляющей сигнала для каждого объекта на входе антенны вы-

числяется на основе данных помехо-сигнальной обстановки по следующей формуле:

   P inoise ray



P* noise

iray

f0 f avg

 fhigh 

100,05R6 A flow R

iray

,

 где

P i* noise ray

– приведенный уровень широкополосной составляющей сигнала в точке излучения,

Па/√Гц.

Результаты тестирования имитатора
С целью проверки работоспособности имитатора были проведены проверки по следующим параметрам:  приведенная мощность помехи на входе антенны;  наклон сплошной части спектра помехи;  частоты формирования дискретных составляющих (ДС) в спектре сигнала;  частоты формирования амплитудных огибающих (АО) в спектре сигнала.
Для проведения проверок были составлены эпизоды, имитирующие взаимное положение носителя ГАС и объектов. При этом взаимное положение объектов описывалось на плоскости с помощью галсов, где галс характеризовался равномерным прямолинейным движением. Длительность каждого эпизода составляла 10 мин.
В табл. 1 приведены описания тестовых эпизодов. Эпизод 1 позволяет провести первые две проверки – проверки приведенной мощности помехи и наклона сплошной части помехи. В свою очередь, эпизод 2, имитирующий один объект, позволяет проверить частоты формирования ДС в спектре сигнала и его АО.

Номер эпизода
1 2

Число объек-
тов
0 1

Описание поведения объекта

Координаты Х, км Y, км
–– –10 17,5

Курс, град
– 90

Скорость, узлов
– 10

Глубина, м
– 5

Описание поведения носителя ГАС

Координаты Х, км Y, км
00 00

Курс, град
0 –90

Скорость, Глубина,

узлов

м

0 100

0 100

Таблица 1. Описание тестовых эпизодов

В табл. 2 приведены полученное и заданное значения давления на входе антенны. Как видно из табл. 2, разность полученного и заданного значения давления не превышает 0,01 дБ, что вполне допустимо для практических приложений. Это позволяет сделать вывод, что проверка приведенной мощности помехи на входе антенны прошла успешно.

Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)

67

ИМИТАТОР СИГНАЛОВ НА ВЫХОДЕ ПРИЕМНЫХ ЭЛЕМЕНТОВ ПАССИВНОЙ …

Моделируемое значение давления, дБ Измеренное значение давления, дБ Разность, дБ

35,56 35,57 0,01

Таблица 2. Результаты проверки давления на входе антенны

Ниже представлены результаты проверки величины наклона сплошной части спектра помехи. На рис. 2 изображен график зависимости наклона сплошной части спектра помехи от частоты в октавах (единица измерения частотного диапазона; 1 октава соответствует интервалу между частотами, логарифм отношения которых при основании 2 равен 1).
В табл. 3 приведены полученное и заданное значения наклона спектра на входе антенны. Нетрудно видеть, что разность заданного и полученного значений наклона спектра помехи не превышает 0,04 дБ, что вполне допустимо для практических приложений. Следовательно, проверка прошла успешно.

Моделируемое значение наклона спектра, дБ/октава Измеренное значение наклона спектра, дБ/октава Разность, дБ

–6,0 –6,04 –0,04

Таблица. 3. Результаты проверки наклона спектра на входе антенны

Рис. 2. Оценка наклона спектра: 1 – спектр сигнала; 2 – оценка наклона спектра

Измеренная частота ДС, Гц 50,05 68,0 101,1 313,15

Моделируемая частота ДС, Гц 51,0 69,0 102,0 314,0

Разность, Гц –0,95 –1,0 –0,9 –1,15

Таблица 4. Результаты проверки частот формирования ДС в спектре сигнала

Рис. 3. Зависимость частоты ДС от времени
68 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)

Ю.В. Шафранюк

В табл. 4 приведены результаты проверки частот формирования ДС в спектре сигнала. Как видно из табл. 4, разность частот формирования ДС и истинных значений частот ДС не превышает 1,15 Гц, что для практических приложений вполне допустимо. Следовательно, проверка прошла успешно.
На рис. 3 приведен график зависимости частот ДС от времени, иллюстрирующий стабильность формирования ДС во времени.
В табл. 5 приводятся результаты проверки частот формирования ДС в спектре АО. Нетрудно видеть, что полученное значение разности не превышает 0,01 Гц, что для практических приложений вполне допустимо. Проверка прошла успешно.

Измеренная частота ДС, Гц 0,19

Моделируемая частота ДС, Гц 0,2

Разность, Гц –0,01

Таблица 5. Результаты проверки частот формирования ДС в спектре АО

Заключение

Разработан и протестирован имитатор сигналов на выходе приемных элементов с гибкими протяженными буксируемыми антеннами, учитывающий особенности излучения, распространения и приема сигналов и помех в реальных морских условиях, влияющие на эффективность работы гидроакустических станций с гибкими протяженными буксируемыми антеннами. Предложены модели блока формирования распределенных помех и блока формирования сигналов локальных источников имитатора. Разработанный имитатор позволяет проводить отработку алгоритмов обнаружения, классификации и определения координат морских шумящих объектов.
Работа выполнена при поддержке Российского фонда фундаментальных исследований (проект 1108-01097-а и 12-08-00511-а).

Литература

1. Корякин Ю.А., Смирнов С.А., Яковлев Г.В. Корабельная гидроакустическая техника. Состояние и актуальные проблемы. – СПб: Наука, 2004. – 143 с.
2. Справочник штурмана / Под ред. В.Д. Шандабылова. – М.: Воениздат, 1968. – 544 с. 3. Урик Р. Дж. Основы гидроакустики. – Л.: Судостроение, 1978. – 445 с. 4. Кудрявцев А.А., Лугинец К.П., Машошин А.И. Oб амплитудной модуляции подводного шумоизлуче-
ния гражданских судов // Акустический журнал. – 2003. – Т. 49. – № 2. – С. 224–228. 5. Акустика океана / Под ред. акад. Л.М. Бреховских. – М.: Наука, 1974. – 693 с. 6. Малышкин Г.С. Оптимальные и адаптивные методы обработки гидроакустических сигналов. Т.1. Оп-
тимальные методы. Т. 2. Адаптивные методы. – СПб: ОАО «Концерн «ЦНИИ «Электроприбор», 2009; 2011. – 400 с.; 374 с. 7. Сергиенко А.Б. Цифровая обработка сигналов. – СПб: Питер, 2002. – 607 с.

Шафранюк Юлия Валентиновна – ОАО «Концерн «ЦНИИ «Электроприбор», аспирант, Julis85@yandex.ru

Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)

69