CC BY
65
Моделирование пространственных спектров на выходе векторно-скалярной приемной
системы, работающей на фоне шумов моря
А.В. Аверьянов, Г.М. Глебова ЮФУ, НИИ физики, г. Ростов-на-Дону
Введение
Современное развитие технологии создания векторно-скалярных многокомпонентных приемников вызывает интерес к использованию векторно-скалярных антенн (ВСА) в различных прикладных задачах гидроакустики. Для ВСА предлагаются алгоритмы обработки с различным пространственным разрешением: стандартным, высоким [1] и сверхвысоким [2]. Сравнительный анализ показал ряд преимуществ ВСА по сравнению со скалярными антеннами (СА). Так, в работе [3] показано, что отношение сигнал/помеха на выходе ВСА более чем в 3 раза превышает это отношение при использовании СА, имеющей такую же апертуру. В настоящей работе предлагается алгоритм обработки сигналов для ВСА, эффективно подавляющий шум моря. Алгоритм основан на измерении потока мощности принимаемых сигналов и учитывает тот факт, что при достаточном числе осреднений поток мощности в горизонтальной плоскости по шумовой составляющей стремится к нулю [4].
При моделировании алгоритма обработки для ВСА, работающей на фоне шумов моря, требуется полное и точное представление авто- и взаимно ковариационных характеристик векторно-скалярных шумовых полей. Если скалярные характеристики шумов моря исследованы достаточно детально, то векторные характеристики шумов, обусловленных взволнованной поверхностью океана, требуют дополнительных исследований. Результаты, полученные рядом авторов, достаточно полно представлены в монографии [4]. Для безграничного полупространства получены аналитические выражения коэффициентов пространственной корреляции и соотношение мощностей на компонентах векторноскалярного приемного модуля. Однако получить полное описание характеристик векторноскалярных шумовых полей на основе имеющихся данных достаточно проблематично. Для этого необходимо вернуться к интегральным соотношениям, описывающим векторноскалярные поля в волноводе, а далее либо получить аналитические выражения подобно тому, как это, например, сделано для некоторых характеристик в работе [5], либо выполнить численное интегрирование, которое в данной работе реализовано посредством компьютерного моделирования. Моделирование шумов, вызванных взволнованной морской поверхностью, позволяет рассчитать пространственные спектры на выходе приемных систем, для которых основной помехой являются шумы моря.
Модели сигнала и шума
Приемная система, которая установлена в глубоком море, представляет собой горизонтальную линейную эквидистантную антенну, состоящую из M модулей. Каждый модуль состоит из приемника давления и двух ортогональных приемников колебательной скорости, оси которых расположены в горизонтальной плоскости. На рис.1 представлена геометрия модели: центр линейной антенны, расположенной параллельно оси X на глубине
H, находится в точке O. Местоположение шумящего источника находящегося в точке I, задается координатами (г,а, ß), где а - пеленг, ß - угол возвышения. Шумящая поверхность ограничена кругом радиуса R .
Рис.1. Геометрическая модель шумящей поверхности.
Принимаемые сигналы на m -ом модуле ВСА можно представить в виде:
Um =[Pm ,Vmx ,Vmy 1 m = \-M, (1)
где pm и Vmr, (r = x, y) - звуковое давление и проекции колебательной скорости по направлениям x и у, измеряемые m -ым приемным модулем. Если вектор U формировать в
виде U = [p1,..., Pm ,Vxi,.-.,VxM./ ,Vyi,.. . ,VM J, то его можно записать как
U = \p,Vx ,Vy J (2)
где p, Vx ,Vy - векторы сигналов по давлению и по компонентам колебательной скорости на
всей приемной антенне.
Поверхностные шумы моря, обусловленные взволнованной поверхностью, создаются совокупностью независимых источников, равномерно распределенных вблизи поверхности в пределах круга конечного диаметра [4, 5]. Амплитуда излучения одинакова для всех источников и распределена по Рэлею, а начальная фаза - случайная величина, равномерно распределенная в интервале значений \0,2^J. Мнимые и действительные части сигналов таких источников распределены по Гауссу. Характеристика направленности источников шума в вертикальной плоскости имеет вид
g = cosn p. (3)
Звуковое давление, создаваемое единичным источником, записывается в виде
P = A\exp(-jkr)lrJ- g . (4)
Компоненты колебательные скорости для точечного источника с угловыми координатами а и p в эквивалентных единицах звукового давления имеют вид
V = P cosasin p ,V = p sinasin p, (5)
который справедлив практически при любых условиях распространения сигнала в волноводе [4]. При моделировании сигнала полагалось, что источники шума являются диполями, а локальный источник - монополь, т.е. в выражении (3) n = 1 и n = 0 соответственно.
Для гауссовых сигналов и шумов с нулевым математическим ожиданием статистика измерений полностью определяется матрицей ковариаций, которая рассчитывается для заданной модели сигналов и помех как (U -U* 'j. Матрица ковариаций для вектора U,
соответствующего фиксированной частоте, представляет собой сумму матрицы сигналов от локальных источников - K (s) и матрицы помех - K (n)
К = (и ■ и ) = К (^ + К (и). (4)
В отличие от традиционного представления матрицы помех в виде К (и) = ЫЕ (Е -единичная матрица), что соответствует пространственно-некоррелированным шумам, в данной работе для более полного описания шумов моря используется компьютерное моделирование. Тестирование результатов моделирования векторно-скалярной шумовой составляющей показало совпадение с теоретическими характеристиками поля помех [4]. Так, средние значения элементов ковариационной матрицы равны
1, 1 < I < М
Кг,г (И) = Р ■
0.25, М +1 < I < 2М . (5)
0.25, 2М + 1 < I < 3М
Откуда следует, что мощности на горизонтальных компонентах колебательной скорости в 4 раза меньше мощности давления. Ортогональные компоненты колебательной скорости не коррелированны между собой, не коррелированны также сигналы на приемнике давления и горизонтальных компонентах колебательной скорости. Но это не означает, что измеряемый поток акустической мощности шума Ж в горизонтальной плоскости равен нулю. Обнуляется средний уровень шума, а измеряемая величина определяется дисперсией потока [4]
Ж]=^2(р)*2(у) = ст<(р), (6)
1 J 2 ■ (Д/т) 8 ■ (А/т)' ''
где Д/т определяет частотную полосу и время наблюдения.
Рассчитанные коэффициенты пространственной корреляции для скалярного поля на линейной антенне также совпали с известными зависимостями для модели шумового поля Крона-Шермана [6], имеющими вид
р(й, Р) = р( у, Р) = 2 ^ (кё) / кё; (7)
где - функция Бесселя первого порядка, ё - расстояние между приемными элементами. Также получено совпадение с аналитическими зависимостями для векторных компонент поля [5]:
р(й, V) = 8[ (кё) / кё - 3 У2 (кё) /( кё)2], (8)
р(й,Гу) = р(x,Ух) = &/2(кё)/(Ы)2. (9)
Данные соотношения справедливы при горизонтальном расположении приемной антенны вдоль оси Х.
Алгоритм обработки
Пространственный спектр на выходе СА при обработке сигналов методом, имеющим стандартное разрешение (метод Бартлетта) рассчитывается следующим образом:
I 12 *
Яр (р) = \Т„ (р) ■ Р = ТЦр)КТг (р), (10)
здесь Т(р) - сканирующий вектор, р - вектор координат, по которым выполняется фокусировка в пространственной области. При работе в дальней зоне вектор р(а,Р ) является вектором задержек, обеспечивающим обзор антенны в фиксированном
направлении по пеленгу и углу возвышения, К - спектрально-корреляционная матрица принимаемых сигналов (4).
В ряде работ [1, 7] пространственного спектра по методу Бартлетта при использовании ВСА предлагается рассчитывать следующим образом
Яи (р)=\ти (р) ■ и|2 = т*(р)КТи (р). (11)
Сканирующие векторы для СА и ВСА, соответственно, имеют вид:
Тр р) = [Р(Р)],
Ти(р) = [р(р\ ух ОХ V (р)].
<
В данной работе в отличие от традиционно используемых алгоритмов для ВСА (11) предлагается метод обработки, имеющий вид:
К ф) = Кр ф) ■ Ку ф). (12)
здесь Яу (ф) - пространственный спектр с использованием векторных приемников,
измеряющих градиент давления. Электронное вращение на ортогональных векторных приемниках обеспечивается комбинацией у ,У, таким образом, чтобы направление
результирующего вектора V совпадало с направлением, формируемым по скалярным элементам
куф) = |рф)■ Уф\ . (13)
Величина К (ф) ■ К (ф) является потоком акустической мощности, поэтому метод (12) назовем методом, работающим по потоку мощности.
Результаты моделирования
В этом разделе приводятся результаты моделирования пространственных спектров по методу Бартлетта на выходе СА (10), ВСА (11), а также по потоку мощности (12). На рис.2 представлены пространственные спектры, когда на вход приемной антенны поступают сигналы, обусловленные только шумами моря. Нормировка пространственных спектров выполнена на среднюю мощность сигнала на входе одиночного модуля. Антенна расположена горизонтально в глубоком море, число приемных модулей равно 20, межэлементное расстояние равно Х/2, средняя рабочая частота /0 = 1000 Гц.
До)
-----1 -----2 - - - 3
0 45 90 135 180 225 270 315 360
а,0
Рис.2. Пространственные спектры на выходе СА и ВСА, обусловленные шумами моря:
1 - для СА, 2 - для ВСА, 3 - для ВСА при формировании отклика по потоку мощности.
Выходной отклик для СА (кривая 1) имеет ярко выраженные максимумы и минимумы. Минимальный уровень шума принимается по направлениям близким к а = ±180°. Это объясняется тем, что лепесток, сформированный в этом направлении, захватывает достаточно удаленную область морской поверхности, на которой расположены шумящие источники. Максимальный уровень шумов наблюдается по направлениям, близким к а = ±90°. В данном случае в силу пространственной симметрии формируемый главный лепесток захватывает не только траверсные направлении в горизонтальной плоскости, но и поверхность моря «над головой» антенны. Кривая 2 соответствует пространственному спектру при использовании ВСА. Поскольку диаграмма направленности горизонтально
расположенных векторных приемников обращена к морской поверхности нулем, то при групповой обработке сигналов по направлениям а = ±90° уровень принимаемых шумов несколько уменьшается по сравнению с СА. Кривая 3 соответствует алгоритму обработки (12). Видно, что использование такого подхода приводит к эффективному подавлению шумов моря при конечном значении А/ ■т. В данном примере А/ ■т =1000, что вполне реализуемо для многих прикладных задач гидроакустики.
На рис. 3 представлены аналогичные пространственные отклики для ситуации, когда входные сигналы определяются не только шумами моря, но и локальным источником, находящимся под углом а = 67°. Нормировка для каждого из рассчитанных пространственных спектров выполнена на его максимальное значение.
Рис3. Пространственные спектры на выходе СА и ВСА, обусловленные шумами моря и сигналом от локального источника: 1- для СА, 2 - для ВСА, 3 - для ВСА при формировании
отклика по потоку мощности.
Поскольку линейная СА не различает «право-лево», то для неё получен известный результат: в выходном отклике на фоне шумов наблюдается два лепестка единичного уровня. В пространственном спектре ВСА ложный боковой лепесток подавлен, шумовой фон немного ниже, чем в пространственном спектре для СА, особенно по направлениям близким к 90° и 270°. Пространственный спектр той же ВСА, но с использованием алгоритма по потоку мощности (12), также подавляет ложный боковой лепесток от локального источника, а самое главное, эффективно подавляет шумы моря.
Для характеристики помехоустойчивости алгоритмов обработки используется отношение сигнал/помеха на выходе приемной системы:
в_{Р(я + п) - РЩ а[Р(п)\ ’
здесь Р(^) - значение выходного сигнала, когда на вход поступает сигнал только от локального источника, Р(п) - значение выходного сигнала, обусловленное только шумом.
С учетом того, что принимаемые сигналы имеют нормальное распределение, шум аддитивен, а входные сигналы подвергаются линейному преобразованию, данное выражение преобразуется к виду
О = Ш
= <Р(п)> ■
Средние значения, а также и P(s) были найдены путем моделирования и статистического осреднения для каждого из рассматриваемых алгоритмов. Как видно из рис.3, значение P(n) зависит от угла фокусировки - а для всех рассмотренных алгоритмов, следовательно и отношение сигнал/помеха на выходе также зависит от углового положения источника. Для ситуации, представленной на рис.3 (а = 67° ), получены следующие значения отношения сигнал/помеха на выходе приемной антенны: для СА D =52, для ВСА D =68, а для алгоритма, работающего по потоку мощности, D =800. В данной ситуации предложенный алгоритм (по критерию отношения сигнал/помеха на выходе приемной системы) эффективнее традиционного алгоритма для скалярной антенны в 16 раз.
Заключение
Разработан алгоритм, который использует поток мощности сигналов, принимаемых векторно-скалярной антенной. Предложенный алгоритм эффективно подавляет шумы моря за счет того, что среднее значение потока мощности в горизонтальной плоскости по шумовой составляющей стремится к нулю. Компьютерное моделирование показало, что даже при конечном времени наблюдения и ограниченной полосе частот (А/ -т =1000) разработанный алгоритм обработки эффективно подавляет шумовую составляющую на выходе пространственного спектра векторно-скалярной антенны.
Анализ пространственных спектров показал существенные преимущества предложенного алгоритма для векторно-скалярной антенны: отношение сигнал/помеха на выходе приемной системы увеличено в 16 раз по сравнению со скалярной приемной системой.
Как следствие столь эффективного подавления шумов правомерно предположить, что рабочие характеристики векторно-скалярной приемной системы, использующей предложенный алгоритм обработки, будут намного лучше, чем у скалярной антенны, что позволит обнаруживать малошумные источники на больших расстояниях и с большей точностью.
Литература
1. Hawkes M., Nehorai A. Acoustic vector-sensor beamforming and Capon direction estimation // IEEE Transactions on Signal Processing. Sep. 1998. V. 46. № 9. P. 2291-2304.
2. Wong K.T., Zoltowski M.D. Self-initiating MUSIC-based direction finding in underwater acoustic particle velocity-field beamspace // IEEE Journal of Oceanic Engineering. Apr. 2000. V. 25. № 2. P. 262-273.
3. А.В. Аверьянов, Г.М. Глебова, Г.Н. Кузнецов, Точность оценки параметров локального источника векторно-скалярными антеннами // Сборник трудов XXII сессии РАО и Сессии Научного совета РАН по акустике. Т. 2. - М.: Геос, 2010, С. 332-335.
4. Гордиенко В.А. Векторно-фазовые методы в акустике. М.: Физматлит, 2007. 480 с.
5. Захаров Л.Н., Киршов В.А., Рожин Ф.В. Пространственно-корреляционные функции компонент колебательной скорости для двух моделей звукового поля // Акуст. журн. 1972. Т. 18. №1. С. 49-52)
6. B.F. Cron, C.H. Sherman, “Spatial-correlation function for various noise models,” J. Acoust. Soc. Am., vol. 34, 1962, pp. 1732.
7. Glebova G.M., Kuznetsov G.N. “Estimating parameters of signal sources and characteristics of noise field by using spatially separated vector-scalar modules,” The Formation of Acoustical Fields in Oceanic Waveguides, Reconstruction of Inhomogeneities in Shallow Water, Nizhny Novgorod: Inst. Appl. Phys. 1998. V. 1. P. 109.
