CC BY
31
Vol. 19, No. 05, 2016
Ovil Aviation High TECHNOLOGIES
УДК 621.396
ВОССТАНОВЛЕНИЕ ФОРМЫ ПОВЕРХНОСТИ ОБЪЕКТА ПО ПОЛЯРИЗАЦИОННОЙ СТРУКТУРЕ ПОЛЯ ОТРАЖЕННОЙ ЭЛЕКТРОМАГНИТНОЙ ВОЛНЫ
А.И. КОЗЛОВ1, В.Ю. МАСЛОВ2
1МГТУ ГА, г. Москва, Россия 2МИРЭА, г. Москва, Россия
Рассматривается задача обратного рассеяния электромагнитной волны на объекте при однопозиционной радиолокации, когда комплексная матрицы рассеяния симметрична. Исследуется, как будет изменяться величина поля, регистрируемого приемной антенной, при плавном изменении параметров матрицы рассеяния. Решается задача обратного рассеяния электромагнитной волны на произвольном трехмерном объекте. Приводятся графики изменения величины модуля поляризационного коэффициента отраженной волны от угла поворота объекта и соответствующие им диаграммы обратного рассеяния. Получено дифференциальное уравнение, которое связывает изменение поляризационного коэффициента отраженной волны с изменением элементов матрицы рассеяния объекта. Полученные соотношения дают возможность разработать алгоритмы быстрого численного решения обратной задачи рассеяния на этом объекте, т. е. определение комплексной функции рассеяния поверхности объекта и восстановления формы неизвестного объекта по распределению фазы отраженной волны. Метод использует лучевые представления полей рассеяния, базирующиеся на принципе Гюйгенса - Френеля. Алгоритм восстановления формы объекта по фазе отраженной волны позволяет восстанавливать не только гладкие поверхности, но и поверхности объектов, имеющие неровности меньшие, чем длина волны.
Ключевые слова: поляризация радиоволн, матрица рассеяния радиоволн, обратная задача радио-поляриметрии.
Рассмотрим задачу обратного рассеяния электромагнитной волны на объекте при одно-позиционной радиолокации, когда комплексная матрицы рассеяния S симметрична (512 = ^21). Тогда в приближении плоских электромагнитных волн связь отраженной объектом электромагнитной волны с падающей на этот объект волной определяется соотношением [ 1]
Eo = SEn
(1)
или в развернутом виде
ГEoi Л С
VEo2 J
ill ¿12
V-12 - 22 J
CE Л
Ь п1
V E п2 J
(2)
Исследуем, как будет изменяться величина поля Е о, регистрируемого приемной антенной, при плавном изменении параметров матрицы рассеяния. Например, это может быть поворот на небольшой угол исследуемого объекта вокруг начала координат или можно зафиксировать исследуемый объект и произвести поворот на небольшой угол вектор облучения объекта. Предполагается, что в обоих случаях величина поля облучения объекта одинакова и поддерживается постоянной.
Рассмотрим матрицу рассеяния S (г) в качестве оператора, у которого элементы матрицы являются функциями некоторого параметра г. Связь отраженной объектом волны с падающей на этот объект волной определяется уравнением
Eo = S (г) Eп. (3)
Научный Вестник МГТУ ГА_Том 19, № 05, 2016
Civil Aviation High TECHNOLOGIES Vol. 19, No. 05, 2016
В этом случае изменение отраженной объектом волны будет связано с изменениями элементов матрицы рассеяния S (t) уравнением
Eo = S (t) Еп . (4)
Из уравнения (3) следует, что
Еп = S-1 (t) Eo( *), (5)
где 8 ) - обратная матрица.
Подставляя соотношение (5) в уравнение (4), получаем дифференциальное уравнение для вектора Е0 (V) отраженной волны
Ё о (V) = 8 (V )8-1 (V )Ео (V) = Ж(/ )Ео (V), (6)
где матрица Ж (*) = )8-1 (г) •
Так как все элементы, входящие в уравнение (6), зависят от параметра V, то в дальнейшем указание на эту зависимость опускаем. Матрица Ж является матрицей дифференциальных форм
W = S S-1 =■
1
С К.. К-- — ¿..я. - — .¿.-.¿..Л f
det S
¿11 ¿22 ¿12 ¿12 ¿11 ¿12 ¿12 ¿11 V¿ 12 ¿22 _ ¿22 ¿12 ¿ 12 ¿12 — ¿22 ¿11J
W11 W12
V W21 W22 J
(7)
В произвольном поляризационном базисе поляризация электромагнитной волны может быть задана отношением ортогонально поляризованных компонент этой волны E01 и E02 и
р = = ре", (8)
Е 01
которое называют поляризационным коэффициентом (поляризационным фазором), а р и ф -модулем и аргументом поляризационного коэффициента. Матричное уравнение (6) можно представить в виде
E01 =W11 E01 +W12E02; E02 = w21 E01 +w22E02 • (9)
Если в системе уравнений (9) перейти к отношению ортогонально поляризованных компонент волны E02 и Eoi, то с учетом соотношения (8) можно получить уравнение для поляризационного коэффициента
p = E02 E01 - E02 E01 = (W21E01 + W22 E02 ) E01 - E02 (W11E01 + W12 E02 ) (10)
— "С2 "С2
E01 E01
После упрощения соотношения (10) получается искомое дифференциальное уравнение для поляризационного коэффициента
Vol. 19, No. 05, 2016
Oivil Aviation High TECHNOLOGIES
Р = -.12 Р2 - .21 + (.22 - .11) Р . (11)
Выражение (11) является комплексным уравнением Риккати, которое связывает изменение поляризационного коэффициента отраженной волны с изменением элементов матрицы рассеяния объекта. Изменения элементов матрицы рассеяния объекта описываются с помощью дифференциальных форм
.11 = _ .22 — _ .12, .12 = _ .12 — _ .11, (12) . 21 = .12 .22 — .22 .12' .22 = .12 .12 — .22 .11.
Для полного описания поведения поляризационного коэффициента отраженной волны необходимо задать три комплексных элемента матрицы рассеяния S (г) и три комплексных величины приращения этих элементов. Варианты видов решения уравнения (11) отличаются большим разнообразием. Кроме того, полный анализ поведения решений уравнения (11) представляет собой отдельную сложную математическую задачу. Поэтому в настоящей статье приходится ограничиться несколькими примерами.
Рассмотрим некоторые характерные решения уравнения (11), в которых скорость изменения элементов матрицы рассеяния .ц, ._12'22 связывается с изменением ракурса наблюдения объекта, т. е. предполагается, что объект совершает в течение времени вращение вокруг одной из своих собственных осей.
На рис. 1-3 приводятся графики изменения величины модуля поляризационного коэффициента отраженной волны от угла поворота объекта и соответствующие им диаграммы обратного рассеяния.
Рассмотренные выше примеры показывают, что, задавая соответствующие коэффициенты в уравнении (11), можно моделировать диаграммы обратного рассеяния реальных объектов. Для сложных и протяженных объектов необходимо будет суммировать определенное количество решений уравнения (11) имеющих различные коэффициенты.
Рис. 1. Изменения модуля поляризационного коэффициента отраженной волны в зависимости от угла поворота объекта и соответствующие им диаграммы обратного рассеяния
Рис. 2. Изменения модуля поляризационного коэффициента отраженной волны в зависимости от угла поворота объекта и соответствующие им диаграммы обратного рассеяния
Civil Aviation High TECHNOLOGIES
Vol. 19, No. 05, 2016
Полученные соотношения дают возможность перейти к решению задачи обратного рассеяния электромагнитной волны на объекте при однопозиционной радиолокации, когда комплексная матрицы рассеяния 8 симметрична (512 = 5 21). Тогда в приближении плоских электромагнитных волн связь отраженной объектом электромагнитной волны с падающей на этот объект волной определяется соотношением [ 1].
Решаемая задача заключается в следующем. Пусть на неизвестный трехмерный объект с гладкой поверхностью падает плоская электромагнитная волна. Поверхность объекта имеет комплексную функцию рассеяния. На некотором известном расстоянии от объекта измеряется электромагнитное поле отраженной волны. Необходимо восстановить форму поверхности.
Метод численного решения этой задачи использует лучевые представления полей рассеяния, базирующиеся на принципе Гюйгенса - Френеля. Метод включает в себя специально разработанные для данного случая алгоритмы быстрого численного прямого и обратного преобразования Гюйгенса - Френеля [1] и алгоритм восстановления формы неизвестного объекта по распределению фазы отраженной волны.
Рассмотрим разработанный в работе [ 1] алгоритм быстрого численного прямого и обратного преобразования Гюйгенса - Френеля, обобщив его на случай восстановления формы объекта. При этом предполагается, что ни для падающей волны, ни для рассеянной нет затенений каких-либо элементов поверхности объекта.
Если на поверхность объекта падает плоская электромагнитная волна, то в точке наблюдения (рис. 3) комплексная амплитуда напряженности электрического поля отраженной волны находится по формуле
где §_Ху - функция рассеяния объекта, р - расстояние между рассеивающей точкой на поверхности объекта и точкой, где определяется величина отраженного поля, к = 2п/Л - волновое число, Л - длина волны, Б(х, у) - контур объекта. Начало координат совпадает с первыми элементами соответствующих двумерных массивов §_Х у и ет п.
Поэтому если поверхность трехмерного объекта описывается уравнением г = /Х у и облучается плоской волной, нормаль к которой совпадает с осью 2, то расстояние
Точки, в которых задается функция рассеяния объекта, и точки, где определяются величины отраженного поля, располагаются в параллельных плоскостях на расстоянии с1 друг от друга. В области, в которой отраженные волны имеют сферический фронт, в сомножителе 1/ р
расстояние р полагается р & г, а расстояние р в показателе экспоненты можно представить в виде двух членов разложения в ряд Тейлора:
ikp
(13)
D (х, ^)
(14)
2r
Vol. 19, No. 05, 2016
Oivil Aviation High TECHNOLOGIES
Рис. 3. Схема взаимного расположения объекта и плоскости наблюдения отраженного от объекта электромагнитного поля
С учетом сделанных предположений комплексная амплитуда вектора напряженности электрического поля отраженной волны определяется как
N N Qikp N N eik (р0 +pf)
em, n = ZISx, y— = ZZS.
y=1 x=1
y=1 x=1
N N eikPo
em, n = ZZ-
x, y
Jkpf
____x, y
1 x=1 r У r
(16) (17)
где размерность двумерных массивов §_х у и етп определяется произведением N х N. Выражение (17) можно представить в виде произведения матриц
G = ФР.
При этом двумерный массив р = S_x y exp
x, y
ik
^ 1 Л f + — f2 2r x,у
(18)
в уравнении (18)
преобразуется в матрицу Р размером N х 1:
P
1( ...
Д£ц £1,2 £1, N £2,1 £
2,2 £N, N
(19)
Матрица Ф размером N2 х N2 имеет элементы
фх' y = exp J— ikd'
-m,n r I 2r
( x-m )2+( y-n )2
(20)
Матрица Ф зависит только от двух величин г и ё, определяющих взаимное положение объекта и точек наблюдения отраженного поля. Поэтому при неизменных величинах г и ё требуется лишь однократный расчет этой матрицы. Этим достигается высокое быстродействие предложенного алгоритма прямого преобразования Гюйгенса - Френеля.
Civil Aviation High TECHNOLOGIES Vol. 19, No. 05, 2016
Полученную в результате расчета по формуле (7) матрицу G размером N х 1
g = ( eu
-1,2
-1, N
-2,1
-2,2
-N, N ,
(21)
следует затем преобразовать в двумерную матрицу размером N х N
G о =
-1,1
VeN,1
-1,N
eN,N )
(22)
На рис. 4 в качестве отражающего объекта изображена поверхность гиперболического параболоида с постоянной функцией рассеяния §_х у. Результаты расчета по формуле (7) комплексной амплитуды напряженности электрического поля отраженной волны ет п на
расстоянии г от поверхности гиперболического параболоида представлены на рис. 5.
Алгоритм обратного преобразования Гюйгенса - Френеля заключается в нахождении матрицы Р размером N2 х 1, элементы которой характеризуются функцией рассеяния объекта §_х у . Эта задача сводится к решению системы линейных алгебраических уравнений вида
ФР = G. (23)
В уравнении (23) матрица Ф определяется из соотношения (19) при заданных геометрических величинах г и ё, определяющих взаимное положение объекта и точек наблюдения отраженного поля. Матрица О размером N х 1 образуется из массива комплексных значений отраженного поля етп (21). Для решения системы (23) можно использовать, например, алгоритм
¿^-разложения. При постоянных величинах г и ё, определяющих взаимное положение объекта и точек наблюдения отраженного поля, требуется однократное нахождение разложения матрицы
Ф = ии, (24)
где Ь - нижняя треугольная, а и - верхняя треугольная матрицы. Этим обстоятельством определяется высокое быстродействие предложенного алгоритма обратного преобразования Гюйгенса - Френеля.
Рис. 4. Отражающая поверхность гиперболического параболоида (N = 20)
Рис. 5. Массив значений модуля комплексной амплитуды напряженности
электрического поля отраженной волны етп N = 20)
Vol. 19, No. 05, 2016
Ovil Aviation High TECHNOLOGIES
Рис. 8. Массив значений аргументов комплексных элементов матрицы Р о гауссовой поверхности (Ы = 30)
Рис. 9. Пример восстановления формы гауссовой поверхности, имеющей неровности (Ы = 30)
Рис. 6 дает представление обо всем массиве значений аргументов комплексных элементов матрицы Ро при N = 20, с опорой на который проводилось восстановление формы поверхности гиперболического параболоида, для чего использовались значения фазы обратного преобразования Гюйгенса - Френеля. Полученная поверхность иллюстрируется рис. 7.
После нахождения искомой матрицы Р ее надо преобразовать в квадратную матрицу Ро размером N х N. На рис. 8 представлен массив фазы комплексных элементов матрицы Ро.
Для нахождения формы поверхности неизвестного объекта используется алгоритм восстановления формы объекта по фазе отраженной волны, основанный на нахождении по аргументам элементов матрицы Ро зависимости абсолютной фазы, которая пропорциональна расстоянию до соответствующей точки объекта крт,„ (к = 2п\к ).
Алгоритм восстановления формы объекта по фазе отраженной волны позволяет восстанавливать не только гладкие поверхности, но и поверхности объектов, имеющие неровности меньшие чем длина волны X. На рис. 9 изображен пример восстановления формы поверхности объекта, имеющего неровности. Восстановление произведено по массиву значений аргументов комплексных элементов матрицы Ро гауссовой поверхности (рис. 9).
Таким образом, разработанный алгоритм дает возможность проводить восстановление формы объекта по фазе отраженной волны. При этом речь может идти не только об объектах с гладкой поверхностью, но и об объектах с неровной поверхностью, размер шероховатостей на которых меньше длины волны.
СПИСОК ЛИТЕРАТУРЫ
1. Козлов А.И., Логвин А.И., Сарычев В.А. Поляризация радиоволн. Кн. 2. Радиолокационная поляриметрия. М.: Радиотехника, 2007.
2. Козлов А.И., Маслов В.Ю. Численный метод решения трехмерной обратной задачи рассеяния электромагнитных волн на препятствии // Научный Вестник МГТУ ГА. 2012. № 179 (5).
Civil Aviation High TECHNOLOGIES
Vol. 19, No. 05, 2016
3. Козлов А.И., Маслов В.Ю. Численный метод определения неоднородной комплексной диэлектрической проницаемости плоской поверхности объектов по поляризационной структуре поля отраженной электромагнитной волны // Научный Вестник МГТУ ГА. 2012. № 179 (5).
4. Козлов А.И., Маслов В.Ю. Дистанционное определение диэлектрической проницаемости поверхности в оптическом диапазоне // Научный Вестник МГТУ ГА. 2014. № 210.
5. Козлов А.И., Маслов В.Ю. Дифференциальные уравнения эволюции матрицы рассеяния // Научный вестник МГТУ ГА. 2014. № 210 (12). С. 43-46.
6. Козлов А.И., Логвин А.И., Сарычев В.А. Поляризация радиоволн. Кн. 3. Радиополя-риметрия сложных по структуре сигналов. М.: Радиотехника, 2008.
7. Козлов А.И., Маслов В.Ю. Дифференциальные свойства матрицы рассеяния // Научный Вестник МГТУ ГА. 2004. № 79. С. 19-25.
8. Маслов В.Ю. Дифференциальные свойства поляризационного коэффициента электромагнитной волны // Научный Вестник МГТУ ГА. 2004. № 79. С. 26-30.
9. Маслов В.Ю. Дифференциальные свойства матрицы когерентности отраженной от объекта электромагнитной волны // Научный Вестник МГТУ ГА. 2004. № 79. С. 31-34.
10. Маслов В.Ю. Матрица когерентности в дифференциальной радиополяриметрии и ее основные свойства // Научный Вестник МГТУ ГА. 2005. № 93. С. 13-18.
СВЕДЕНИЯ ОБ АВТОРАХ
Козлов Анатолий Иванович, Заслуженный деятель науки и техники РФ, профессор, доктор физико-математических наук, профессор кафедры ТЭРЭО ВТ МГТУ ГА, электронный адрес: vilandes@yandex.ru.
Маслов Виктор Юрьевич, доктор технических наук, профессор МТУ(МИРЭА), электронный адрес: vilandes@yandex.ru.
SHAPE RESTORATIONS OF OBJECT SURFACE ON POLARIZATION STRUCTURE OF REFLECTED ELECTROMAGNETIC WAVE FIELD
Anatoliy I. Kozlov
Moscow State Technical University of Civil Aviation, Moscow, Russia, vilandes@yandex.ru
Viktor Yu. Maslov
Moscow Technological University, Moscow, Russia, vilandes@yandex.ru ABSTRACT
A problem of electromagnetic wave backscattering on a chosen 3D object is solved. A differential equation which is linked change of polarization coefficient of reflected wave with variation of matrix elements of object scattering is obtained. Obtained relation enables to develop algorithms of fast numerical solution of inverse problem of scattering on this object that is determination of complex function of object surface scattering and restoration of unknown object shape on phase distribution of reflected wave. The method uses ray representation of scattering fields based on principle Huygens-Fresnel. The algorithm of object shape restoration on phase of reflected wave allows to restore not only smooth surfaces, but also object surfaces with smaller roughness than a wave length.
Key words: radiowave polarization, matrix of radiowave scattering, inverse problem of radiopolarization.
REFERENCE
1. Kozlov A.I., Logvin A.I., Sarychev V.A. Polyarizaciya radiovoln. Kn. 2. Radiol okacion-naya polyarimetriya. - M.: Radiotexnika, 2007. (in Russian)
Vol. 19, No. 05, 2016
Civil Aviation High TECHNOLOGIES
2. Kozlov A.I., Maslov V.Yu. Chislennyj metod resheniya trexmernoj obratnoj zadachi rasseyaniya elektromagnitnyx voln na prepyatstvii. Nauchnyj Vestnik MGTU GA, 2012, № 179 (5). (in Russian)
3. Kozlov A.I., Maslov V.Yu. Chislennyj metod opredeleniya neodnorodnoj kompleksnoj diel-ektricheskoj pronicaemosti ploskoj poverxnosti obektov po polyarizacionnoj strukture polya otrazhen-noj elektromagnitnoj volny. Nauchnyj Vestnik MGTU GA, 2012, № 179 (5). (in Russian)
4. Kozlov A.I., Maslov V.Yu. Distancionnoe opredelenie dielektricheskoj pronicaemosti poverxnosti v opticheskom diapazone. Nauchnyj Vestnik MGTU GA, 2014, № 210. (in Russian)
5. Kozlov A.I., Maslov V.Yu. Differencialnye uravneniya evolyucii matricy rasseyaniya. Nauchnyj Vestnik MGTU GA, 2014, № 210 (12), s. 43-46. (in Russian)
6. Kozlov A.I., Logvin A.I., Sarychev V.A. Polyarizaciya slosznich sygnalov. Kn. 3. Radiop-olyarimetriya slozhnyh po strukture signalov. M.: Radiotehnika, 2008. (in Russian)
7. Kozlov A.I., Maslov V.Yu. Differencialnie svoistva matrici rasseiania. Nauchnyj Vestnik MGTU GA, 2004, № 79, s. 19-25. (in Russian)
8. Maslov V.Yu. Differencialnie svoistva polariacionnogo koefficienta EMW. Nauchnyj Vestnik MGTU GA, 2004, № 79, s. 26-30. (in Russian)
9. Maslov V.Yu. Differencialnie svoistva matrici kogerentnosti otrazheennoj EMW. Nauchnyj Vestnik MGTU GA, 2004, № 79, s. 31-34. (in Russian)
10. Maslov V.Yu. Matritsa kogerentnosti v differentsialnoy radiopolyarimetrii I eyo osnovnye svoystva. Nauchnyj Vestnik MGTU GA, 2005, № 93, s. 13-18. (in Russian)
