Радиолокация и радионавигация
УДК 621.396.06 Оригинальная статья
https://doi.org/10.32603/1993-8985-2021 -24-4-48-56
Усечение падающего поля в задаче рассеяния электромагнитных волн на случайных поверхностях конечной длины
В. В. Леонтьев, М. А. Бородин^
Санкт-Петербургский государственный электротехнический университет "ЛЭТИ" им. В. И. Ульянова (Ленина), Санкт-Петербург, Россия
н boroda84@gmail.com
Аннотация
Введение. При анализе рассеяния радиоволн на случайных поверхностях часто используют интегральные уравнения, решаемые численными методами, работоспособными лишь при ограничении размеров поверхности. Указанное ограничение приводит к появлению краевых токов, которые, в свою очередь, вызывают существенные ошибки в расчетах эффективной площади рассеяния (ЭПР), особенно при скользящих углах облучения. Для снижения влияния краевых токов обычно используют функцию усечения падающего поля, которая должна удовлетворять совокупности требований - обеспечивать малое подавление поля по всей площади поверхности конечного размера между ее краями и одновременно снижение амплитуды поля до пренебрежимо малых значений при приближении к этим краям. Падающее поле с используемой функцией усечения должно удовлетворять волновому уравнению с минимальной погрешностью. Известны различные функции усечения падающей волны (гауссовская, Торсоса, интегральная), однако ни одна из них полностью не удовлетворяет предъявляемым требованиям. Цель работы. Предложить новую функцию усечения амплитуды электромагнитной волны, падающей на возмущенную поверхность ограниченного размера, для расчета ее ЭПР. По сравнению с известными функциями усечения новая функция должна удовлетворять всей совокупности предъявляемых требований. Материалы и методы. Выполнено сравнение предложенной функции усечения амплитуды падающего поля с известными функциями усечения, в том числе проведена оценка погрешности удовлетворения волновому уравнению. Для доказательства применимости предложенной функции усечения в среде Ма^аЬ выполнено математическое моделирование бистатической диаграммы рассеяния двумерной морской поверхности конечной длины с пространственным спектром Эльфохейли методом Монте-Карло. Результаты. Установлено, что предложенная функция усечения амплитуды падающего поля удовлетворяет предъявляемым требованиям лучше известных функций. Результаты математического моделирования показали, что новая функция обеспечивает приемлемую точность оценки ЭПР случайных поверхностей конечной длины.
Заключение. Получена новая функция усечения амплитуды падающего поля для уменьшения влияния краевых токов на точность оценки ЭПР двумерных случайных поверхностей конечной длины, которая рекомендуется к использованию при решении задачи рассеяния.
Ключевые слова: функция усечения амплитуды падающего поля, эффективная площадь рассеяния, морская поверхность, рассеяние радиоволн, скользящий угол облучения
Для цитирования: Леонтьев В. В., Бородин М. А. Усечение падающего поля в задаче рассеяния электромагнитных волн на случайных поверхностях конечной длины // Изв. вузов России. Радиоэлектроника. 2021. Т. 24, № 4. С. 48-56. doi: 10.32603/1993-8985-2021 -24-4-48-56
Конфликт интересов. Авторы заявляют об отсутствии конфликта интересов.
Статья поступила в редакцию 20.04.2021; принята к публикации после рецензирования 25.05.2021; опубликована онлайн 29.09.2021
Контент доступен по лицензии Creative Commons Attribution 4.0 License This work is licensed under a Creative Commons Attribution 4.0 License
© Леонтьев В. В., Бородин М. А., 2021
Radar and Navigation
Original article
Tapering the Incident Field When Solving Problems of Electromagnetic Wave Scattering Over Finite-Size Random Surfaces
Viktor V. Leonfev, Mikhail A. BorodinH
Saint Petersburg Electrotechnical University, St Petersburg, Russia
H boroda84@gmail.com
Abstract
Introduction. An analysis of radio wave scattering over random surfaces frequently involves integral equations, which are solved by numerical methods. These methods are feasible only provided limited dimensions of the surface. The requirement of surface limitation leads to the appearance of edge currents, resulting in significant errors when calculating the radar cross section (RCS), particularly for grazing incident angles. The influence of edge currents is reduced by a function tapering the incident field amplitude. This function should satisfy the following requirements: to provide a low suppression of the field along the entire finite-size surface between its edges at the same time as decreasing the incident field amplitude to negligible values when approaching the edges. The incident field under the application of the tampering function should satisfy the wave equation with a minimum error. Although various tapering functions are applied for incident field amplitude (i.e. Gaussian, Thorsos, integral), none of them satisfies the aforementioned requirements.
Aim. To suggest a novel function for tapering the amplitude of an electromagnetic wave incident on a perturbed finite-size surface when calculating RCS. In comparison with the known functions, the proposed function must satisfy the entire set of requirements.
Materials and methods. A comparison of the proposed tapering function for incident field amplitude with the known tapering functions was performed, including the estimation of the error of satisfying the wave equation. To prove the applicability of the proposed tapering function, a mathematical modeling of the bistatic scatter diagram of a two-dimensional sea-like finite surface with a spatial Elfouhaily spectrum was carried out using Monte Carlo calculations in the Matlab environment.
Results. Compared to the known tapering functions, the proposed tapering function satisfies the entire set of requirements. The results of mathematical modeling showed that the proposed function for tapering the incident field amplitude provides acceptable accuracy of estimating the RCS of finite-size random surfaces. Conclusion. A novel function for tapering the incident field amplitude was derived. This function reduces the influence of edge currents on the accuracy of RCS estimation of two-dimensional finite-size random surfaces, thus being instrumental for solving scattering problems.
Keywords: tapering function for incident field amplitude, radar cross section, sea surface, radio wave scattering, grazing incident angle
For citation: Leont'ev V. V., Borodin M. A. Tapering the Incident Field When Solving Problems of Electromagnetic Wave Scattering Over Finite-Size Random Surfaces. Radioelectronics. 2021, vol. 24, no. 4, pp. 48-56. doi: 10.32603/1993-8985-2021 -24-4-48-56
Conflict of interest. The authors declare no conflicts of interest. Submitted 20.04.2021; accepted 25.05.2021; published online 29.09.2021
Введение. Одной из основных задач анализа распространения электромагнитных волн является задача их рассеяния на различных поверхностях. Рассмотрим такую задачу в плоскости, проходящей через точки расположения источника излучения и отражения от поверхности, считая последнюю шероховатой, т. е. имеющей различные
отклонения г точек поверхности для различных координат поверхности х: г = г (х) (рис. 1). Будем считать поверхность идеально проводящей шероховатой, представляющей собой реализацию случайной морской волны. Начало 0 декартовой прямоугольной системы координат 0хг размещено в точке зеркального отражения, имеющей г = 0.
49
Усечение падающего поля в задаче рассеяния электромагнитных волн на случайных поверхностях конечной длины
Tapering the Incident Field When Solving Problems of Electromagnetic Wave Scattering Over Finite-Size Random Surfaces
а б
Рис. 1. Геометрия задачи рассеяния: а - при горизонтальной поляризации (TE); б - при вертикальной поляризации (TM)
Fig. 1. Geometry of scattering problem: a - with horizontal polarization (TE); б - with vertical polarization (TM)
В эту точку электромагнитная волна падает под углом 0, отражение происходит под углом 9р. Углы
отсчитываются от осей 0х и 0z соответственно в положительном направлении оси x по часовой стрелке. В частном случае зеркального отражения (или рассеяния вперед) 9р = 9отр = л/ 2 — 9.
Направление обратного отражения, важное для од-нопозиционных радиолокационных систем, определит угол 9р = 9 — л/ 2.
Задачи рассеяния решаются с помощью интегральных уравнений. В случае бесконечной поверхности в качестве падающего поля рассматривают плоскую волну
Уп (x, z, t) = = У о exp (—/'raí) exp [ik (x cos 9 — z sin 9)], (1)
где x и z — координаты точек поверхности; t — время; Уо — амплитуда; ra — круговая частота; k = 2 л/Л — волновое число (Л — длина электромагнитной волны); 9 — угол скольжения. Далее множитель exp (—irat), учитывающий зависимость поля от времени, опускается.
В рассматриваемом двумерном случае в зависимости от поляризации удобнее работать с напряженностью электрического поля E либо с напряженностью магнитного поля H. При горизонтальной (TE) поляризации (рис. 1, а) выбирают напряженность электрического поля [1] и в (1) Уп = Еп, при вертикальной (TM) — работают с напряженностью магнитного поля и в (1) Уп = Нп (рис. 1, б).
На практике интегральные уравнения приходится решать численными методами, в силу чего любую поверхность необходимо ограничивать. Так, на рис. 1 введено ограничение длины поверхности в пределах ± L¡ 2 относительно начала координат. Ограничение поверхности ведет к появлению краевых токов, которые, в свою очередь,
приводят к существенным ошибкам в расчетах эффективной площади рассеяния (ЭПР). Особенно эта проблема актуальна при скользящих углах облучения. Для снижения влияния краевых токов обычно используют усечение падающей волны. Цель усечения - посредством взвешивания функцией усечения обеспечить уменьшение ее амплитуды при приближении к краям поверхности конечного размера до пренебрежимо малых значений, оптимально - до нуля. В этом случае в качестве падающего поля рассматривают волну [1]:
Уп ( x, z ) = = УоТ ( x, z ) exp[ik (x cos 9 — z sin 9)],
(2)
где Т (х, г) - функция усечения.
Из (2) следует, что падающее поле больше не является плоской волной, а его энергия распределена в узком луче, сосредоточенном около среднего угла скольжения 9.
Требования к функции усечения предъявляются не только с учетом изменения амплитуды падающей волны (2) на краях поверхности. Необходимо, чтобы она обеспечивала малое подавление поля по всей длине поверхности между ее краями. Кроме того, поле (2) с используемой функцией усечения должно удовлетворять волновому уравнению (по крайней мере, с малой погрешностью).
Цель настоящей статьи - предложить новую функцию усечения амплитуды электромагнитной волны, падающей на возмущенную поверхность ограниченного размера при расчете ее ЭПР. По сравнению с известными функциями усечения новая функция должна удовлетворять всей совокупности предъявляемых требований.
Первой и часто используемой функцией усечения является гауссовская функция [2-5]:
\2
Tr (x, z) = exp
( x + z ctg 9)
(3)
Усечение падающего поля в задаче рассеяния электромагнитных волн
на случайных поверхностях конечной длины Tapering the Incident Field When Solving Problems of Electromagnetic Wave
Scattering Over Finite-Size Random Surfaces
где g - параметр, контролирующий ширину луча. Торсос показал [2], что поле (2) с усечением (3) представляет собой грубое приближение решения волнового уравнения Гельмгольца:
(V2 + к2 (х, г) = 0, (4)
где V2 = б2/йг2 + б2/&2 - оператор Лапласа.
В силу того, что поле (2) с усечением (3) является грубой аппроксимацией реального падающего поля, Торсос предложил функцию (3) модернизировать [2]:
\2
Тт (x, z) = exp
(x+z ctg e)2
2
g
.k (x cos e-z sin e)rr/ ч - г---U (x, z)
где
U (x, z ) =
(kg sin e)
2 (x + z ctg e)2
(5)
-1
Функцию (5) принято называть функцией усечения Торсоса [1]. Из сравнения (3) и (5) следует, что гауссовская функция и функция Торсоса совпадают по амплитуде, а суть усечения Торсоса состоит в коррекции фазового сдвига облучающего поля.
Падающее поле (2) с функцией усечения Торсоса (5) удовлетворяет волновому уравнению (4) более точно, чем поле (2) с гауссовской функцией усечения (3).
Поле (2) с функцией усечения Торсоса (5) является аппроксимацией интегрального представления падающего поля [2]:
( x, z ) =
1
VVAe
V 2
í exP
-V 2
(e -e)2
(Ae)2
+ ik (x cos e - z sin ei)] de-
(6)
где A9 = 2/ (kg sin 9). Будучи суперпозицией плоских волн, поле (6) удовлетворяет волновому уравнению (4), однако, как будет показано далее, обеспечивает при прочих равных условиях существенно меньшее уменьшение амплитуды на краях, чем усечение (5) Торсоса.
Оценке значения контролирующего ширину луча параметра g для представленных моделей в литературе уделено значительное внимание. Обычно рекомендуют выбирать его из диапазона от i/10 до L/4 [5].
Сравнив рассмотренные модели при одних и тех же условиях задачи, можно заключить, что ни
одна из них полностью не удовлетворяет предъявляемым требованиям. Модель (6) удовлетворяет волновому уравнению, но меньше двух других моделей уменьшает амплитуду падающей волны на краях и сильно уменьшает ее по координате х между краями ("внутри") поверхности. Модель поля с гауссовской функцией усечения очень плохо удовлетворяет волновому уравнению. По сравнению с ней модель с функцией усечения Торсоса удовлетворяет уравнению лучше, вызывает такое же уменьшение амплитуды волны на краях, но сильно уменьшает ее по координате х "внутри" поверхности.
Синтез новой функции усечения. Введем новую функцию усечения, которая будет лучше удовлетворять сформулированным требованиям. Расчет ЭПР морской поверхности базируется на методе Монте-Карло, при котором ее случайные реализации генерируют математическим моделированием. Так как профиль поверхности априори неизвестен, сравнивать функции усечения принято на плоской поверхности (при г = 0). В связи с этим сначала запишем выражение для новой функции усечения в зависимости только от координаты х:
( x ) = sin ( xg ) .
(7)
Для гарантированного уменьшения амплитуды падающей волны до нуля на краях поверхности примем g = L¡ (2 я). Тогда (7) примет следующий вид:
, ч sin (2rcx/L) . , , ч , N
T (x) = —i-= Sinc (2roc/L). (8)
s 2nx/L v ' '
Рис. 2 иллюстрирует влияние типа модели падающей волны на характер изменения ее амплитуды вдоль координаты x для плоской поверхности (z = 0). Линия 1 соответствует падающему полю (2) с предложенной функцией усечения (8), линия 2 -с функцией усечения Торсоса (5), линия 3 - интегральному представлению (6). Для всех трех моделей поля параметр, контролирующий ширину луча, одинаков и определен как g = L/(2я). Значения остальных параметров задачи: L = 10 м; 9 = 2°; А, = 0.03 м.
Убывание значений рассмотренных функций усечения (рис. 2) от центра к краям поверхности различно и носит нелинейный характер. Значения предложенной функции усечения (8) в большинстве точек "внутри" рассматриваемой поверхности превышают значения функций усечения (5) и (6).
Усечение падающего поля в задаче рассеяния электромагнитных волн на случайных поверхностях конечной длины
Tapering the Incident Field When Solving Problems of Electromagnetic Wave Scattering Over Finite-Size Random Surfaces
+
I* п
1 4,/i
/ /
//0.6
.>/ 0.4 "" /
2 / 0.2 V
w " \\
V \
\
ч
-5 -4 -3 -2 -1 0 1 2 3 4 x, м Рис. 2. Модуль падающего поля c функцией усечения: 1 - (8); 2 - (5); 3 - (6)
Fig. 2. Absolute value incident wave of the tapering function: 1 - (8); 2 - (5); 3 - (6)
На краях поверхности при |x| = 5 м амплитуда волны принимает следующие значения: для интегрального представления (6) |*п| = 0.244, при
функции усечения Торсоса (5) |*п| = 5 -10-5. Согласно определению (8) значения предложенной функции усечения на краях рассматриваемой поверхности равны нулю, что, в отличие от известных функций, обеспечивает значения *п| = 0 в этих точках.
Таким образом, по критериям уменьшения амплитуды волны на краях до нуля и меньшего ее подавления "внутри" поверхности функция усечения (8) является наилучшей из рассмотренных.
Попутно отметим, что, хотя волна с функцией усечения Торсоса (5) и хуже, чем волна (6), удовлетворяет волновому уравнению, по сравнению с последней она обеспечивает лучшее подавление амплитуды на краю. По этой причине ранее во многих работах (см., например, [5-9]) в качестве функции усечения использовали функцию Торсоса.
Следует обратить внимание на еще один важный факт. При аппроксимации интегрального представления падающего поля, в результате которой получена функция усечения Торсоса, введен ряд ограничений [1]:
g >
(3...4)У2 .
k0 sin 0 ' kg sin 0 > 1;
|x + zctg9| «с 7ig; \z - xctg Q kg sin 0
g
2
(9)
(10) (11)
(12)
Эти ограничения не позволяют выбирать параметры задачи произвольно. Они устанавливают
жесткие связи между параметром g, углом скольжения 9, длиной волны электромагнитного поля Л, параметрами поверхности.
Современные радиолокационные станции (РЛС) имеют высокое разрешение по дальности. В этом случае необходим расчет ЭПР поверхностей малой длины. Особенно актуален такой расчет при разработке алгоритмов обработки сигналов в РЛС, предназначенных для экологического мониторинга морской поверхности с целью обнаружения разливов нефти [10]. Морские судовые навигационные РЛС сантиметрового диапазона длин волн облучают поверхность моря под малыми углами скольжения. Анализ показывает, что для таких РЛС удовлетворить требованиям (9)-(12) сложно, что ограничивает область применения функции усечения Торсоса.
Проверка на удовлетворение волновому уравнению. Проверим представление падающего поля (2) с функцией усечения (8) на соответствие волновому уравнению. Подставив (2) с (8) в (4), получим:
(У2 + к2 (х, г ) = к \, (13)
где 58 - некий параметр. Поясним его физический смысл и приведем для него аналитическое
выражение. Очевидно, что к > 0 и правая часть выражения (13) стремится к нулю только при условии стремления к нулю 58. По этой причине параметр 58 можно назвать погрешностью, с которой
волна не удовлетворяет волновому уравнению. Для волны (2) с функцией усечения (8) погрешность определяется следующей формулой:
где
A =-
B =
Ss =|(A + iBi)*п(x, z),
4л2 4 л cos ( x/g ) ( xk )2 ( kL)2 xk 2L sin ( x/g )' 4лcos(0)cos(x/g) 2cos(0)
(14)
kL sin (x/g )
xk
Аналогичная погрешность для волны (2) с функцией усечения Торсоса (5) имеет вид [4]:
5Т =|(^2 + В )^п (х, г )|, (15)
где A2 = -W2 -
4 2 2 16x cos 0
k 4 g 8sin6 0;
52 Усечение падающего поля в задаче рассеяния электромагнитных волн
на случайных поверхностях конечной длины Tapering the Incident Field When Solving Problems of Electromagnetic Wave
Scattering Over Finite-Size Random Surfaces
-5 -4 -3 -2 -1 0 1 2 3 4 x, м Рис. 3. Оценка погрешности удовлетворения волновому уравнению
Fig. 3. Estimation of the error of satisfying the wave equation
B =
4kx cos 9
(kg sin 9)
1 -
4x
W =
2 (x/g )2 -1 (kg sin 9)2
Сравним погрешности (14) и (15). При этом будем полагать, что параметр, контролирующий ширину луча, для указанных моделей поля одинаков и определяется как g = Ь/(2л). Другие параметры определены так же, как при расчете кривых на рис. 2.
Изменение погрешностей по оси х показано на рис. 3. Штриховая линия соответствует погрешности для волны с функцией усечения Торсоса, сплошная - умноженной на 100 погрешности для волны с предложенной функцией усечения (8). Масштабирование второй зависимости выполнено для отображения обеих зависимостей на одном рисунке. Максимальное значение погрешности для волны с функцией усечения Торсоса 5ттах = 0.4, для волны с функцией усечения (8) -5 зтах = 0.0026.
Сравнение показывает, что волна с предложенной функцией усечения (8) удовлетворяет волновому уравнению существенно лучше, чем волна с функцией усечения Торсоса.
Обобщим введенную функцию усечения с учетом координаты г поверхности. Каждая реализация поверхности при моделировании случайна, и значения ее ординат на краях априори не известны. В этих условиях функция усечения должна обеспечивать на краях равную нулю амплитуду падающей волны, поэтому приобретает сложный вид:
/ 4_ísinc^i (x, z), -L/ 2 < x < 0; s (x z) = { sincU2 (x, z), 0 < x < L2,
(16)
где U1 2 =
x±z(x)ctg9 +Z/2 ± z (+Z/2) ctg 9'
-5
-0.04 -0.06
а
F п1
Л sí A
/л '!
( >Ш 0.8 —
f / * v
1/ \ \
0.6 — \1
/ 1
Л f /
f\f /
0.4 — V
oí f V Л\ /л
XV /
jT / 1 К Л
0.2 — 1 \ 1 1 / м
\ /
V lili 1 1 1 1
-5 -4 -3
-2 -1
0
б
1
2 3
4 x, м
Рис. 4. Реализация морской поверхности (а) и функции усечения (б)
Fig. 4. Realization of the sea surface (a) and tapering function (б) причем z(+Z/2) - значения функции z — z{x) при x = + L¡2.
На рис. 4, а представлена реализация волнового профиля случайной морской волны, полученная с помощью спектральной модели [11-13]. В качестве пространственного спектра моря использован спектр Эльфохейли [14-15], лучше других спектров описывающий как гравитационные, так и гравитационно-капиллярные волны. Моделирование выполнено при следующих значениях параметров: L = 10 м; скорость ветра на высоте 19.5 м ув = 2 м/с; среднеквадратическое отклонение ординат волнового профиля стz = 0.02 м.
На рис. 4, б сплошной линией изображена реализация функции усечения (16), соответствующая реализации волнового профиля морской волны, представленной на рис. 4, а. Штриховой линии на рис. 4, б соответствует функция усечения (8).
Расчет бистатической диаграммы рассеяния. Определение плотности поверхностного тока (ППТ), наведенного падающим полем (2) с функцией усечения (16), базируется на решении интегрального уравнения для электрического поля (EFIE) [1], [5]. По ППТ вычисляется напряженность электрического поля рассеянной электромагнитной волны и бистатическая диаграмма рассеяния (ДР) - зависимость модуля коэффициента
Усечение падающего поля в задаче рассеяния электромагнитных волн 53
на случайных поверхностях конечной длины
Tapering the Incident Field When Solving Problems of Electromagnetic Wave Scattering Over Finite-Size Random Surfaces
1Г1, дБ
-20 -10
-90 -80 -70 -60 -50 -40 -30
|Г|, дБ -20 -10
80
Г:
-100 б
Рис. 5. Бистатическая диаграмма рассеяния реализации морской поверхности: а - по одной реализации; б - усредненная по 100 реализациям
Fig. 5. A bistatic scatter diagram of the realization of the sea surface: а - one realizations at a time; б - averaged over 100 realizations
рассеяния Г от угла рассеяния 9р при фиксированном угле скольжения 9:
1Г(0р )!=№, Ч^. (17)
Бистатическая ДР для реализации волнового профиля морской волны, представленной на рис. 4, а, изображена на рис. 5, а.
Численный расчет ДР (17) является приближенным. В [7] предложен критерий оценки точности приближенных численных решений дифракционных задач:
V 2
i
-V 2
Г(0р )| d% = 1.
(18)
Степень отклонения от 1 в (18) служит мерой погрешности всех вычислений. Для представленной на рис. 5, а случайной реализации бистатиче-ской ДР значение интеграла в (18) равно 1.1, что соответствует погрешности расчета 10 %.
На рис. 5, б представлена бистатическая ДР, полученная усреднением 100 ДР для различных реализаций волнового профиля случайной морской волны.
Моделирование волнового профиля случайной морской волны выполнено при тех же значениях па-
раметров, что и для рис. 4, а. Средняя (по 100 реализациям) погрешность расчета ЭПР равна 11.8 %.
Критерий (18) является необходимым, но не достаточным. В качестве достаточного критерия справедливости расчета ДР может служить ее соответствие физическим представлениям. По критерию Рэлея моделируемая поверхность близка к гладкой, и при скользящем облучении следует ожидать сильного отражения в зеркальном направлении.
Из рис. 5, б следует, что при угле скольжения 9 = 2° большая часть рассеянной энергии сосредоточена в области зеркального отражения в лепестке ДР шириной 2.6° вокруг угла зеркального отражения 9отр =л/ 2 -9 = 88°. Максимум ДР
|г(9р )| = 11.5 дБ наблюдается при 9р = 86.3°.
Для угла зеркального отражения 9отр = 88° модуль коэффициента рассеяния равен 8.35 дБ, в то время как для направления обратного отражения (при угле 9р = 9 - V 2 = -88°) |Г(9р )| = -84.9 дБ.
Заключение. Предложенная функция усечения (16) не только лучше других удовлетворяет предъявляемым к ним требованиям, но и обеспечивает приемлемую точность оценки ЭПР случайных поверхностей ограниченной длины.
Усечение падающего поля в задаче рассеяния электромагнитных волн
на случайных поверхностях конечной длины Tapering the Incident Field When Solving Problems of Electromagnetic Wave
Scattering Over Finite-Size Random Surfaces
а
Список литературы
1. Toporkov J. V., Awadallah R. S., Brown G. S. Issues related to the Use of a Gaussian-like Incident Field for Low-grazing-angle Scattering // J. Optical Society of America A. 1999. Vol. 16, № 1. P. 176-187. doi: 10.1364/JOSAA.16.000176
2. Thorsos E. The validity of the Kirchhhoff Approximation for Rough Surface Scattering using Gaussian Roughness Spectrum // J. Acoustical Society of America. 1988. Vol. 83, № 1. P. 78-92. doi: 10.1121/1.396188
3. Pan G., Zhang L. Closed Form Solution to the Incident Power of Gaussian-Like Beam for Scattering Problems // IEEE Trans. on antennas and propagation. 2019. Vol. 67, № 2. P. 1364-1367. doi: 10.1109/TAP.2018.2884851
4. Zhang Y., Wang Y., Zheng H. EM Scattering from a Simple Water Surface composed of Two Time-varying Sinusoidal Waves // Proc. of IEEE Intern. Conf. on Computational Electromagnetics (ICCEM). 2020. Vol. 8. P. 200684-200694. doi: 10.1109/C0MPEM.2019.8779021
5. Scattering of Electromagnetic Waves: Numerical Simulations / L. Tsang, J. A. Kong, K.-H. Ding, C. O. Ao. New York: John Wiley & Sons, 2001. 736 p. doi: 10.1002/0471224278
6. Ye H., Jin Y.-Q. Parameterization of the Tapered Incident Wave for Numerical Simulation of Electromagnetic Scattering from Rough Surface // IEEE Trans. on Antennas and Propagation. 2005. Vol. 53, № 3. P. 1234-1237. doi: 10.1109/TAP.2004.842586
7. Borodin M. A., Leont'ev V. V. Analysis of the Accuracy of an Iterative Algorithm for Calculating the Field Scattered by a Rough Surface // J. of Communications Technology and Electronics. 2009. Vol. 54, № 9. P. 989-994. doi: 10.1134/S1064226909090034
8. Numerical Analysis of Microwave Scattering from Layered Sea Ice Based on the Finite Element Method / X. Xu,
C. Brekke, A. P. Doulgeris, F. Melandso // Remote sensing. 2018. Vol. 10, № 9. P. 1-16. doi: 10.3390/rs10091332
9. Jun M., Guo L-X., Zeng H. Study on 1D Large-scale Rough Surface EM scattering at Low Grazing Incident Angle by Parallel MOM based on PC Clusters // Wave in Random and Complex Media. 2009. Vol. 19, № 4. P. 585-599. doi: 10.1080/17455030903033190
10. Леонтьев В. В., Бородин М. А., Игнатьева О. А. Бистатические диаграммы рассеяния морской поверхности, покрытой мономолекулярной пленкой нефти // Радиотехника. 2012. № 7. C. 39-44.
11. Toporkov J. V., Brown G. S. Numerical Simulations of Scattering from Time-varing Randomly Rough Surfaces // IEEE Trans. on Geoscience and Remote Sensing. 2000. Vol. 38, № 4. P. 1616-1624. doi: 10.1109/36.851961
12. Johnson J. T., Toporkov J., Brown G. A Numerical Study of Backscattering from Time-evolving Sea Surfaces: Comparison of Hydrodynamic Models // IEEE Trans. on Geoscience and Remote Sensing. 2001. Vol. 39, № 11. P. 2411-2420. doi: 10.1109/36.964977
13. Numerical Simulation of Backscatter from Linear and Nonlinear Ocean Surface Realization / C. L. Rino, T. L. Crystal, A. K. Koide, H. D. Ngo, H. Guthart // Radio Science. 1991. Vol. 26, № 1. P. 51 -71. doi: 10.1029/90RS01687
14. A Unified Directional Spectrum for Long and Shot Wind-driven Waves / T. Elfouhaily, B. Chapron, K. Katsa-ros, D. Vandemark // J. of Geophysical Research. Oceans. 1997. Vol. 102. № C7. P. 15781-15796. doi: 10.1029/97JC00467
15. Bourlier C., Saillard J., Berginc G. Intrinsic Infrared Radiation of the Sea Surface // Progress in Electromagnetics Research. 2000. Vol. 27. P. 185-335. doi: 10.2528/PIER99080103
Информация об авторах
Леонтьев Виктор Валентинович - доктор технических наук (2000), профессор (2003) кафедры радиотехнических систем Санкт-Петербургского государственного электротехнического университета "ЛЭТИ" им. В. И. Ульянова (Ленина). Почетный работник высшего профессионального образования РФ. Автор более 140 научных и учебно-методических работ. Сфера интересов - радиолокация, распространение и рассеяние радиоволн, радиотехнические системы мониторинга окружающей среды.
Адрес: Санкт-Петербургский государственный электротехнический университет "ЛЭТИ" им. В. И. Ульянова (Ленина), ул. Профессора Попова, д. 5, Санкт-Петербург, 197376, Россия E-mail: vvleontyev@mail.ru http://orcid.org/0000-0002-6800-7124
Бородин Михаил Анатольевич - кандидат технических наук (2011), доцент кафедры радиотехнических систем Санкт-Петербургского государственного электротехнического университета "ЛЭТИ" им. В. И. Ульянова (Ленина). Автор 29 научных работ. Сфера интересов - радиолокация, распространение и рассеяние радиоволн, радиотехнические системы мониторинга окружающей среды.
Адрес: Санкт-Петербургский государственный электротехнический университет "ЛЭТИ" им. В. И. Ульянова (Ленина), ул. Профессора Попова, д. 5, Санкт-Петербург, 197376, Россия E-mail: boroda84@gmail.com http://orcid.org/0000-0002-5237-9118
Усечение падающего поля в задаче рассеяния электромагнитных волн 55
на случайных поверхностях конечной длины
Tapering the Incident Field When Solving Problems of Electromagnetic Wave Scattering Over Finite-Size Random Surfaces
References
1. Toporkov J. V., Awadallah R. S., Brown G. S. Issues related to the Use of a Gaussian-like Incident Field for Low-grazing-angle Scattering. J. Optical Society of America A. 1999, vol. 16, no. 1, pp. 176-187. doi: 10.1364/JOSAA.16.000176
2. Thorsos E. The validity of the Kirchhhoff Approximation for Rough Surface Scattering using Gaussian Roughness Spectrum. J. Acoustical Society of America. 1988, vol. 83, no. 1, pp. 78-92. doi: 10.1121/1.396188
3. Pan G., Zhang L. Closed Form Solution to the Incident Power of Gaussian-Like Beam for Scattering Problems. IEEE Trans. on antennas and propagation. 2019, vol. 67, no. 2, pp. 1364-1367. doi: 10.1109/TAP.2018.2884851
4. Zhang Y., Wang Y., Zheng H. EM Scattering from a Simple Water Surface composed of Two Time-varying Sinusoidal Waves. Proc. of IEEE Intern. Conf. on Computational Electromagnetics (ICCEM). 2020, vol. 8, pp. 200684-200694. doi: 10.1109/C0MPEM.2019.8779021
5. Tsang L., Kong J. A., Ding K.-H., Ao C. O. Scattering of Electromagnetic Waves: Numerical Simulations. New York, John Wiley & Sons, 2001, 736 p. doi: 10.1002/0471224278
6. Ye H., Jin Y.-Q. Parameterization of the Tapered Incident Wave for Numerical Simulation of Electromagnetic Scattering from Rough Surface. IEEE Trans. on Antennas and Propagation. 2005, vol. 53, no. 3, pp. 1234-1237. doi: 10.1109/TAP.2004.842586
7. Borodin M. A., Leont'ev V. V. Analysis of the Accuracy of an Iterative Algorithm for Calculating the Field Scattered by a Rough Surface. J. of Communications Technology and Electronics. 2009, vol. 54, no. 9, pp. 989-994. doi: 10.1134/S1064226909090034
8. Xu X., Brekke C., Doulgeris A. P., Melandso F. Numerical Analysis of Microwave Scattering from Layered
Sea Ice Based on the Finite Element Method. Remote sensing. 2018, vol. 10, no. 9, pp. 1-16. doi: 10.3390/rs10091332
9. Jun M., Guo L-X., Zeng H. Study on 1D Large-scale Rough Surface EM scattering at Low Grazing Incident Angle by Parallel MOM based on PC Clusters. Wave in Random and Complex Media. 2009, vol. 19, no. 4, pp. 585-599. doi: 10.1080/17455030903033190
10. Leont'ev V. V., Borodin M. A., Ignatieva O. A. Bi-static scattering diagrams of the sea surface covered with a monomolecular film of oil. Radio Engineering. 2012, no. 7, pp. 39-44. (In Russ.)
11. Toporkov J. V., Brown G. S. Numerical Simulations of Scattering from Time-varing Randomly Rough Surfaces. IEEE Trans. on Geoscience and Remote Sensing. 2000, vol. 38, no. 4, pp. 1616-1624. doi: 10.1109/36.851961
12. Johnson J. T., Toporkov J., Brown G. A Numerical Study of Backscattering from Time-evolving Sea Surfaces: Comparison of Hydrodynamic Models. IEEE Trans. on Geo-science and Remote Sensing. 2001, vol. 39, no. 11, pp. 2411 -2420. doi: 10.1109/36.964977
13. Rino C. L., Crystal T. L., Koide A. K., Ngo H. D., Guthart H. Numerical Simulation of Backscatter from Linear and Nonlinear Ocean Surface Realization. Radio Science. 1991, vol. 26, no. 1, pp. 51 -71. doi: 10.1029/90RS01687
14. Elfouhaily T., Chapron B., Katsaros K., Vandemark D. A Unified Directional Spectrum for Long and Shot Wind-driven Waves. J. of Geophysical Research. Oceans. 1997, vol. 102, no. C7, pp. 15781-15796. doi: 10.1029/97JC00467
15. Bourlier C., Saillard J., Berginc G. Intrinsic Infrared Radiation of the Sea Surface. Progress in Electromagnetics Research. 2000, vol. 27, pp. 185-335. doi: 10.2528/PIER99080103
Information about the authors
Viktor V. Leont'ev, Dr. Sci. (Eng.) (2000), Professor (2003) of the Department of Radio Engineering System of Saint Petersburg Electrotechnical University. Honorable Worker of High Professional Education of the Russian Federation. The author of more than 140 scientific and teaching publications. Area of expertise: radiolocation, propagation and scattering of radio waves, radio engineering system for ecology monitoring.
Address: Saint Petersburg Electrotechnical University, 5 Professor Popov St., St Petersburg 197376, Russia
E-mail: vvleontyev@mail.ru
http://orcid.org/0000-0002-6800-7124
Mikhail A. Borodin, Cand. Sci. (Eng.) (2011), assistant professor of the Department of Radio Engineering System of Saint-Petersburg Electrotechnical University. The author of 29 scientific publications. Area of expertise: radiolocation, propagation and scattering of radio waves, radio engineering system for ecology monitoring. Address: Saint Petersburg Electrotechnical University, 5 Professor Popov St., St Petersburg 197376, Russia E-mail: boroda84@gmail.com http://orcid.org/0000-0002-5237-9118
56 Усечение падающего поля в задаче рассеяния электромагнитных волн
на случайных поверхностях конечной длины Tapering the Incident Field When Solving Problems of Electromagnetic Wave
Scattering Over Finite-Size Random Surfaces