CC BY
20Journal of Siberian Federal University. Engineering & Technologies, 2019, 12(2), 138-145
yflK 534.01
Method for Assessing the Distance to the Underwater Object of the Curvature of the Wave Front in the Conditions of Exposure to Correlated Noise
Sergei V. Shostak, Paul A. Starodubtsev* and Roman N. Aliafanov
Far Eastern Federal University 8 Sukhanova Str., Vladivostok, 690090, Russia
Received 29.10.2018, received in revised form 15.01.2019, accepted 10.02.2019
The article presents the results of the substantiation of a new methodfor estimating the distance to any underwater object or physical phenomenon based on the curvature of the wave front in terms of the impact on the measuring sonar system of correlated noise. At the same time, the main task of space-time signal processing is the analysis of the resulting wave field to determine the position of a spatial object, which is represented by a point model of a reflecting (radiating) object, which creates a spherical wave in a homogeneous unlimited medium. Such a model is fundamental for analyzing signals from real sea objects, since in many cases these objects are well described by the model in the form of a certain set of "bright" points. In this case, the curvature of the wave front can be used to determine the distance to the object when it is in the Fresnel zone relative to the receiving antenna array, which allows for independent processing along the time and spatial coordinates.
Keywords: wavefront curvature, correlated noise, bright points, Fresnel zone, temporal and spatial coordinates.
Citation: Shostak S.V., Starodubtsev P. A., Aliafanov R.N. Method for assessing the distance to the underwater object of the curvature of the wave front in the conditions of exposure to correlated noise, J. Sib. Fed. Univ. Eng. technol., 2019, 12(2), 138145. DOI: 10.17516/1999-494X-0122.
© Siberian Federal University. All rights reserved
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0). Corresponding author E-mail address: spa1958@mail.ru
Способ оценки дальности до подводного объекта по кривизне волнового фронта в условиях воздействия коррелированного шума
СВ. Шостак, П.А. Стародубцев, Р.Н. Алифанов
Дальневосточный федеральный университет
Россия, 690090, Владивосток, ул. Суханова, 8
В статье представлены результаты обоснования нового способа оценки дальности до любого подводного объекта или физического явления по кривизне волнового фронта в условиях воздействия на измерительную гидроакустическую систему коррелированного шума. При этом основной задачей пространственно-временной обработки сигналов такой системой является анализ результирующего волнового поля для определения положения пространственного объекта, который представлен точечной моделью отражающего (излучающего) объекта, который создает в однородной безграничной среде сферическую волну. Такая модель основная для анализа сигналов реальных морских объектов, так как эти объекты во многих случаях хорошо описываются моделью в виде некоторого набора «блестящих» точек. Кривизну волнового фронта можно использовать для определения дальности до объекта, когда он находится в зоне Френеля относительно приемной антенной решетки, что позволяет проводить независимую обработку по временной и пространственной координате.
Ключевые слова: кривизна волнового фронта, коррелированный шум, блестящие точки, зона Френеля, временные и пространственные координаты.
Введение
Измерительные гидроакустические системы (ИГС) получают информацию об удаленных подводных объектах (ПО) путем анализа волновых полей (ВП), которые создаются ими в результате собственного излучения или отражения от них зондирующих сигналов. Структура и параметры принятого ВП в области, где осуществляется такой анализ, зависят от положения ПО относительно этой области и его характеристик. Такое ВП несет практически всю информацию об его источнике. В то же время для получения информации о ПО используется не все волновое поле, излучаемое им, а только та его часть, которая попадает в раскрыв приемной антенны ИГС. Воздействуя на элементы антенны, ВП образует пространственно-временной сигнал (ПВС), который обрабатывается ИГС.
Основной задачей пространственно-временной обработки (ПВО) сигналов ИГС является анализ результирующего ВП для определения положения наблюдаемого ВП дальности и направления на ПО. Принципы построения ИГС и методы анализа качества их работы определяются хорошо известными теориями пространственно-временной обработки сигналов [1, 2]. Основные положения этих теорий основаны на анализе пространственно-временных соотношений между габаритными размерами области наблюдения ИГС, т.е. раскрыва приемной антенны Ь, и расстоянием до наблюдаемого объекта. При этом известные в теории и на практике методы ПВО сигналов получены в основном для случая, когда принимаемый волновой фронт считается практически плоским в пределах раскрыва антенны [1].
Вместес темстремлениеразработчиков ИГС кувеличению дальности их действия и повышению разрешающей способности по угловым координатам при наблюдении за ПО и связанная с ними тенденция к увеличению размеров антенн приводят к тому, что в ряде случаев пренебрегать кривизной волнового фронта в пределах раскрыва приемной антенны практически невозможно.
С точнизрееся криввмныволновогофроита, обребачываемого СВЧ, вьюобласть пространства, где могут находиться источники сигналов, можно разделить на три зоны, отсчитываемые относительно приемной антенны [2]:
1) ближнюю,кьгда г|, <О,38£— ,
сивг -дальнсслсо бъе ксь;Ч-еабаритныЧ (елстет анте нны;Т-длинаволны;
21}
2) зощоФренеля^огдс тбз < т3ф <-;
21}
3) дальнююзону,сьгда гдз >-.
X
В зависимости от нахождения ПО в одной из указанных зон в ИГС применяется соответствующая модель формирования ПВС.
Рассмотрим новое математическое решение для оценки дальности до ПО по кривизне волнового фронта (ВФ), находящегося в зоне Френеля, на основе анализа пространственной частоты.
Примем в качестве модели ПО точечную модель отражающего (излучающего) объекта, который создает в однородной безграничной среде сферическую волну. Такая модель является основной для анализа ПВО сигналов реальных морских объектов, так как эти объекты во многих случаях хорошо описываются моделью в виде некоторого набора «блестящих» точек [3].
Соотношения, справедливые для точечного объекта, применимы также и для малоразмерных объектов при условии, 1^тоих величины много меньше элемента разрешен ия системы ПВО сигналов [1-3].
Для теоретического описания процесса формирования ПВС в зоне Френеля антенной решетки построим графическую модель (рис. 1).
На рис. 1 г - дальность ,г0 - лальностьобъекла ,т- номер элемента АР,-0,5 (М -1)<т<0,5(М -1); М - число элементов, d - расстояние между элементами АР, п - номер отсчета во временной области, п = 0...Ы -1, N - числоотсчетов.
Так как зондитующий еигнал предполагается гармоиическим, то ПВС в т-м элементе АР определяетсявыраже нием
Основная часть
(1)
f '
где a - амплитуда сигнал. /0 = — - нормированная частота; At _ частота сианала. fà - ^г^стота
Ле
дискретизации; m - номер элемента в антенной и>еше тке (АР)-
Д-- даАьнейших рассуждсний значение нмплитуды не играет существенной роли, поэтому будем рассматривать ПНС вида
s (m, n ) = exp
У2п/о
n-f+
dm r
(2)
МодельПВСв зоне Френелямежшт 6a>iTi> лреднтавлена слеаьующмм об]эгл;?(5м [2;]:
s (m,n) = exp
J 2ж Л
( (
r c
1 | dm 1 + -1-
2 l r^
2 ^ ^
: eXp(]2жЛ • n) • eXp(-ш'2пм0 • m2 + jме),
(3)
f -2 r 2ее еде u0 = —-; p = 2л/0- — = к0-г0 - начальная фаза; k0 =--волновое число; Xg - длина вол-
ны.
2г0с
Как видим из вырожиHho(h), ПВСфнАиоризуеесяв видепроизведенив временшПипро-странственнойсостаилеме щпй a имее евид04]
s(a,w) = se (и) • s2 (а),
(4)
юи
_ m-1 _ t*.
/"
//ф и ¡)
) в )
, I I
^ /
//н\\ч
п\\
\
w\
W
/
<M-1 _*
/
Рис. 1. Модель формировения ПВСезонеФренеля>штеннойрешетке Fig. 1. Model of the formation of space-time signal in the Fresnel zone of the antenna array
- 141 -
0
m
d
n
где ^ (и) = ехр(д'2тг/0 -и); .)2 (т) = ехср:)ф-рфжи^ • т2 + ).
Причем, как видим из (3) и (4), информацию о дальности несет пространственная составляющая s2 (т).
Полагая, чтоПО малоподвижен, и учите1вая, что (п) = ехр(}2ж/0 ■ п), для выделения со-ставляюще й s2 (да)умиожим и (т, п) на ^ (и), где е<*» - знио кхмплексного с опряжения;
щ (т) = Л]*(п)•.?(т,п) = ехр(-_/2;гм0-т2 + апк). (5)
Выражение дляо2 (ж)включкет квндратичную пространственную переменную т2, что усложняео ьбраоотну таного сичнало для вьгделеянм из номомеличины г0. Для упрощения обработки воспользуемся следующим приемом: сдвинем у2 ВП идр ооптетев и перемножим с исхорнзду квмпиекнно сопряженныо
д(о) = (о - р] • в; (о) = ехрКт2о- м2г/0['' о - _/2пыо • _р2) =
= ехр('2пу0 • л + рЧ"^, (6)
где; у0 = 2р• и0 = рр^—; ¥ = -2пи0 • р2 - ночаввпоя Ьааза.
Лого
В результате поделанных преобразований получим гармонический сигнал частоты V0 и задача определения дальности до ПО свелась к определению пространственной частоты V,,, откуда [5]
(7)
В настоящее время известно большое количество методов оценки частоты гармонического сигнала [6]. Для примера рас смотрим метод дискретного преобразования Фурье (ДПФ) для оценки V,
0'5(М-Д ( 2ЯТ
2:(о)я 2° о - е>°;Е51 -д'—V — 1 =
-о,5(м-1)+р V 21 )
0-Им-Д и С 2я Л
= 22 ехр ((2^0 • м + -То ехр I-У—V • м 1= (8)
-0,5(2--1)+ р V МП )
0,5(М-1)
= exp(/;F)- Z exp
-0,5(М-1)+р
-j—(v-M • v0)
Модуль выражения (8) имеет максимальное значение, когда V = М • V,. И тогда согласно выражению (7)[7]
^ = 40^ (9)
Приведенные выше рассуждения не учитывают влияние окружающей среды на формирование ПВС.
Однако в реальной обстановке ПВС s (т, п) испытывает искажающее воздействие окружающей среды в виде аддитивной составляющей ю (т, п), куда отнесем шум элементов АР и помехигидроакустическогоканала
X (м,пТ = $(м,пО + Си(м,п), (10)
где в(да,тг)- стационарныйгаус совский шумснулевымсреднимирорреляционной матрицей Сии = е||\к(т,п) • (т,п«мл- гильбераовс сопряжение, ОУ (да, п) - вектор шума. Предетавим сформированыый ПВС АР в лекторно-матричном виде [8]
" x(-0,5(M -1),0) " ■ ,,(0) 0 • - 0 " " Ц-0,5(М-1),0) "
x(-0,5(M - 1), (N-1)) со (- 0,5(M - 1),N+-1)
S1 (N-1) 0 - 0 "s2 (-0,5 (M -1))"
■x (0,0) = 0 "1 (0) • 0 (0) + ю(0,0)
x(0, NN -1) 0 "1 ((N 1) • 0 m(0,Ne 1 )
0 0 • )(2) ■s2(0,5(M-
x (0,5(М-1),0) 0 0 • • (-З е Р (0,5(М-1),0)
д:((,5(Л4-1), N-1) и о> ( 0 , 5 (M - 1 ) , N - 1.)
(11)
Перепишемвыражение (11) is более компактномвиде
x = h-o+w,
(12)
где x - (МЫ х 1) вектор ПВС; н - (МЫ х М) известная матрица, формируется в приемнике; 0 -(М х 1) вектнр оценивеемлк пкраметров; \Т - (МЫ х 1) веотор шума.
ТреХметья получикн несмещк-ную хуцттс^пестм вектора 0 с минималрннй дикперсией. Для: получения такой оценки всспо льзуомся теоремой Гаусса)Марковн [0>, К-)], котормя утве рждает,что если данные предстсвнеются йинийной модеусю вора (ИТ, а которой гсумноой СЛ спрн-
извольной плотностью функции распредвления, нулевым средним и известной ковариационной матртцев С[ии=1° Щ\)- • УЮН - «о» - опирацие вильрерто во сопряжение, по луршая некмесцем-ная оценкавклассе линвйкьох оценокимвет ]вид
Ц^С-н^-НС. х.
Приэтомминимальнаядисперсияоценокдаетсявыражением
уаг(т) = Г(ннс!.н)-
(13)
(14)
Выражение [А]„,в означает выделение диагональных илбментов -вадратной метрицы А. Если шумовая составляющая XV представляется белым гауссовским шумом с ковариационной матрицей с2 • I, где I - единичная матрица, с2 - дисперсия, то выражение (13) и (14) преобразуется соответстве нно [5,6]в
0 = (HH xH) •HH X
var(вт ) = a2 T(HH • H
(15)
(16)
Выполняя вычисления в соответствии (13) или (15), фактически мы получим оценки s2 (m) с минимальной дисперсией. Дальнейшие расчеты для получения окончательного результата,
и
т.е. определения дальности до объекта, необходимо проводить в соответствии с выражением (6), (8), (9).
Заключение
Как показано в изложенных материалах, кривизну волнового фронта ПВС можно использовать для определения дальности до ПО, когда он находится в зоне Френеля относительно приемной АР. Модель образования ПВС (3) позволяет факторизовать его на временную и пространственную составляющие (4), что позволяет проводить независимую обработку по временной и пространственной координате. Кроме того, это дает возможность представить сформированный на АР ПВС реальный сигнал в виде линейной векторно-матричной модели (11) и применить хорошо известные теоремы минимизации дисперсии оценок (13), (15). Однако обработка по пространственной координате серьезно осложняется наличием квадратичной переменной (5). Введенный прием сдвига и перемножения позволяет представить сигнал по пространственной координате в виде гармонического, пространственная частота которого зависит от дальности до ПО (6). В результате проделанных операций задача определения дальности ПО в зоне Френеля свелась к оценке пространственной частоты (8).
Список литературы
[1] Гусев В.Г.. Системы пространственно-временной обработки гидроакустической информации. Л.: Судостроение, 1988. 262 c. [Gusev V.G. Systems of space-time processing of hydroacoustic information. L.: Shipbuilding, 1988 262 p. (in Russian)]
[2] Кремер И.Я. Пространственно-временная обработка сигналов. М.: Радиосвязь, 1984. 224 с. [Kremer I.Y. Spatio-temporal signal processing: M.: Radio communication, 1984. 224 p. (in Russian)]
[3] Бурдик В.С. Анализ гидроакустических систем. Л.: Судостроение, 1988 391 с. [Burdick V.S. Analysis of sonar systems. L.: Shipbuilding, 1988. 391 p. (in Russian)]
[4] Марпл С.Л. Цифровой спектральный анализ и его приложения. М.: Мир, 1990. 265 с. [Marpl S.L. Digital spectral analysis and its applications. - M.: Mir, 1990. 265 p. (in Russian)]
[5] Graybill F.A. Theory and Application of the Linear Model, Duxbury Press, North Scituate, Mass., 1976. 704 p.
[6] Рао С.Р. Линейные статистические методы и их применения.. Пер. с англ. под ред. Лин-ника. М.: Наука, 1968 547 p. [Rao S.R. Linear statistical methods and their applications. Translate from English ed. Linnik. - M.: Science, 1968. 547 с. (in Russian)]
[7] Больбасова Л.А. Адаптивная коррекция атмосферных искажений оптических изображений на основе искусственного опорного источника. М.: ФИЗМАТЛИТ, 2012. 126 c. [Bolbasova L.A. Adaptive correction of atmospheric distortions of optical images based on an artificial reference source. M.: FIZMATLIT, 2012. 126 p. (in Russian)]
[8] Канев Ф.Ю., Лукин В.П. Адаптивная оптика. Численные и экспериментальные исследования. М.: Изд. Института оптики атмосферы СО РАН, 2005, 250 с. [Kanev F.Yu., Lukin V.P. Adaptive optics. Numerical and experimental studies. М.: Ed. Institute of Atmosp heric Optics SB RAS, 2005, 250 p. (in Russian)]
[9] Воронцов М.А., Кудряшов А.В., Шмальгаузен В.И. Гибкие зеркала для адаптивных систем атмосферной оптики. Теоретический анализ. Изв. вузов: Радиофизика, 27(11), 1419-1430. [Vorontsov M.A., Kudryashov A.V., Schmalgauzen V.I. Flexible mirrors for adaptive atmospheric optics systems. Theoretical analysis. Izv. universities: Radiophysics, 27(11), 1419-1430. (in Russian)]
[10] Шелдакова Ю.В., Кудряшов А.В., Рукосуев А.Л., Черезова Т.Ю. Использование гибридного алгоритма управления биморфным зеркалом для фокусировки светового излучения. Оптика атмосферы и океана, 2007, 20(4), 380-383. [Yu.V. Sheldakova, A.V. Kudryashov, A.L.Rukosuev, T.Yu. Cherezova, Using a hybrid algorithm for controlling a bimorph mirror to focus light radiation. Optics of the Atmosphere and the Ocean, 2007, 20 (4), 383-383 (in Russian)]
