Метод решения задачи множественного рассеяния электромагнитных волн в однородной безграничной среде
Автор: Шарфарец Б.П.
Журнал: Научное приборостроение @nauchnoe-priborostroenie
Рубрика: Оригинальные статьи
Статья в выпуске: 4 т.14, 2004 года.
Бесплатный доступ
В статье предлагается подход к решению задачи множественного электромагнитного рассеяния на совокупности тел в однородном безграничном пространстве. Для этого рассматривается задача о множественном рассеянии двух тел, находящихся в первичном поле плоской волны. Предполагаются известными исходные невозмущенные матрицы рассеяния каждого рассеивателя. Решение строится путем учета многократного перерассеяния плоских волн между рассеивателями. Получены интегральные уравнения, позволяющие рассчитать результирующие матрицы рассеяния каждого из них и совокупную матрицу рассеяния системы, состоящей из двух рассеивателей. Показано, что решение этой задачи позволяет решить задачу о поле рассеяния произвольного числа рассеивателей. Приведены выражения для поля рассеяния в случае произвольного первичного поля. Продемонстрирована связь интегральных уравнений для задач о множественном рассеянии в однородном пространстве и о множественном рассеянии одиночного тела вблизи границы раздела. Приведены приближенные выражения для расчета поля рассеяния при множественном рассеянии.
Короткий адрес: https://sciup.org/14264361
IDR: 14264361
Текст научной статьи Метод решения задачи множественного рассеяния электромагнитных волн в однородной безграничной среде
Задаче множественного рассеяния волн, в частности электромагнитных волн, когда в некотором первичном поле находится не менее двух рассеивателей, а их взаимным влиянием друг на друга (вторичным перерассеянием) пренебречь нельзя, посвящено большое число публикаций. Упомянем лишь некоторые из них — [1–15].
Рассмотренная впервые в работе [ 1 ] для взаимодействия двух рассеивающих тел в задаче дифракции на щели проблема множественного рассеяния получила свое значительное развитие. Укажем лишь некоторые методы ее решения.
Так, в работе [ 2 ] дан обзор метода, когда решение ищется в виде суммы первичной волны и результирующих полей всех рассеивателей. На границах рассеивателей накладываются краевые условия и условия излучения на бесконечности для рассеянного поля. В работах [ 3-6 ] рассматривается метод решения задачи множественного рассеяния на объемных рассеивателях путем сведения ее к решению интегральных уравнений. В работе [ 7 ] метод, развитый в работах [ 3-6 ] , применяется для решения соответствующей задачи в плоских волноводах. В работах [ 8-12 ] решение ищется с привлечением теорем сложения для различных решений векторного уравнения Гельмгольца, позволяющих аналитически решать задачу множественного рассеяния. Наконец в работах [ 14, 15 ] рассматриваемая задача решается методом Т-матриц, развитым в работе [ 13 ] для случая одиночного рассеивателя.
В настоящей работе для решения задачи множественного электромагнитного рассеяния на нескольких телах предлагается альтернативный подход, заключающийся в учете многократно переот-раженных (перерассеянных) между телами плоских волн. Учет перерассеянных плоских волн осуществляется с помощью использования невозмущенных матриц рассеяния каждого из тел, что приводит к интегральным уравнениям, позволяющим вычислять возмущенные матрицы рассеяния через невозмущенные. Ранее сходная методика была применена автором при решении аналогичной акустической задачи [16], а также задач расчета результирующего акустического рассеяния протяженных излучателей и рассеивателей в условиях влияния границ раздела и слоистой неоднородности среды [17–19].
Единственное ограничение предлагаемого подхода заключается в том, что должны быть известны невозмущенные матрицы рассеяния тел, однако это требование не представляется чрезмерным, т. к. "паспортом" рассеивателя является именно его матрица рассеяния в условиях однородной безграничной среды при отсутствии посторонних тел.
ПОСТАНОВКА ЗАДАЧИ
Без ограничения общности рассмотрим задачу для вакуума в случае трехмерного неограниченного пространства. Вначале сформулируем математически прямую задачу рассеяния для случая единичного рассеивателя, занимающего ограничен- ную область D с R 3, находящегося в поле первичной электромагнитной плоской волны. Суммарное поле, как обычно, рассматривается в виде суммы падающей ui и рассеянной волны us (i в индексах — признак принадлежности к падающей волне (от "incident")). Обычно рассматриваются две классические постановки прямой задачи рассеяния: а) рассеиватель представляет собой прозрачную ограниченную неоднородную область D с R3; б) рассеиватель представляет собой непрозрачное неоднородное включение, ограниченное поверхностью dD, где dD — граница области D.
Пусть падающая волна имеет вид (временнóй фактор полагаем равным exp( -itot ))
E i( x ) = E 0i ex P ( i k i x X
H i( x ) = —Vx E i( x ) = k i x E 0i e xp ( i k i x ) .
ik
Здесь и везде i — мнимая единица.
Следуя в основном [20], запишем выражения для падающей (2–4) и рассеянной (5–8) волн:
Vx E - ik H = 0,
Vx H + ikn ( x ) E = 0, x e R 3 ;
E ( x ) = E 0i exp ( i k i x ) + E s( x ),
H ( x ) = ki x E 0i exp ( i k i x ) + H s ( x );
lim( H s x x - R E s) = 0,
R ^~ lim(E s x x + RH s) = 0 .
R ^~
Здесь Ei , Hi и Es , Hs — компоненты падающей и рассеянной электромагнитных волн соответственно; n = e / e0 — показатель преломления; e(x) (xe D) и e0— диэлектрические посто- янные неоднородной среды и окружающего вакуума соответственно (магнитная проницаемость полагается постоянной во всем пространстве x = = (x,y,z) e R3); R = |x|; ki = kki — волновой вектор падающей волны; k = to /c — волновое число; kˆ i — единичный вектор, совпадающий по на- правлению с вектором ki (далее везде вектор с "крышечкой" обозначает единичный вектор); c — скорость света в вакууме; E0i — постоянный вектор поляризации электрического поля падающей плоской волны. Условия излучения Сильвера—Мюллера (4) служат для выделения единственного решения прямой задачи рассеяния, каждое из них является следствием другого (обычно ограничиваются одним из них), и для векторного уравнения Гельмгольца они переходят в условия излучения Зоммерфельда для компонентов векторного поля (см., напр., [21]).
Для рассеянной волны:
Vx E - ikH = 0,(5)
Vx H + ik E = 0, x e R 3\ D ;
E(x) = E0i exP(ikix) + Es (x),
0i i s
H ( x ) = k i x E 0i exp ( i k i x ) + H s( x ) ;
r(E(x), H(x)) = 0, x edD;(7)
lim( H s x x - R E s) = 0,
R ^~ lim(Es xx + RHs) = 0.
Здесь условие (7) символически определяет краевые условия на границе непрозрачного рассеивателя D .
В дальнейшем везде будем полагать, что решения прямых задач рассеяния (2–4) или (5–8) для одиночного рассеивателя существуют и единственны. Расчеты будут вестись для электрической составляющей поля, магнитная составляющая вне рассеивателя может быть найдена из уравнений Максвелла (2)
H = — VxE.
ik
ПРЕДСТАВЛЕНИЯ ПОЛЯ РАССЕЯНИЯ
Пусть для определенности начало координат находится внутри области D . Как известно ([22] для скалярного случая и [23] для случая электромагнитного поля), поле излучателя (рассеивателя), компоненты которого удовлетворяют условиям излучения Зоммерфельда, можно единственным образом представить вне полосы ^ D = { x , У e R 2, z 2 ^ z ^ z 1 } , включающей этот излучатель, в виде двумерного интеграла Фурье. Здесь z 1 и z 2 соответственно максимальная и минимальная аппликаты области D . Используя, к примеру, результаты работы [22], справедливые для решения скалярного уравнения Гельмгольца, легко выписать выражение для векторного потенциала электромагнитного поля рассеяния, являющегося решением векторного уравнения Гельмгольца при произвольном первичном поле
A s ( x ) = Ja ( k ^ 2 exp ( i k m x ) d k x d k ,
2n a
R 2
m = 1,2; x e R 3 \ Q D .
Здесь k 1 = (kx, ky, a) и k 2 = (kx, ky ,-a) — сопря- женные волновые векторы; а = (к2 - kx2 - ky2)1/2; к — волновое число; m = 1 соответствует полю в верхнем полупространстве z > z 1, m = 2 — в нижнем полупространстве z < z2; a(k m) — спектр векторного потенциала.
Если применить к каждой плоской волне в (9) тождество из [ 24 ]
E = ik (A + ^ W * A), то для рассеянного электрического поля получим:
E 5 = JTk m ) exp ( i k m x ) d kx d ky ,
2 п RR 2 а (10)
m = 1,2, xeR 3\QD, где векторы
T ( k m ) = ik ( a ( k m ) + W * a ( k m ))
k 2
представляют собой векторы, лежащие в плоскости, нормальной к векторам k m (см., напр., [ 23 ] ). Таким образом, как видно из (10), падающая волна вызывает сплошной спектр плоских рассеянных волн не только однородных k x 2 + k y 2 < k , но и неоднородных k x 2 + k y 2 > k .
Асимптотическая оценка интегралов (10) при R ^^ в направлении ( 0 , ф ) ( k x 2 + ky 2 < k 2 , что соответствует так называемой области видимости), где R = |х| и ( R,0,ф ) — сферические координаты, дает [ 22-24 ]
дальнего поля рассеивателя) [25, 26]. Для определения матрицы рассеяния введем следующие системы координат, связанные с векторами падающей k i и рассеянной k s плоских волн. Пусть k i и k s — соответствующие им единичные вектора, тогда остальные два вектора этих систем координат вычисляются следующим образом [ 26 ] :
/X /X /X /X h is = z X k is, v is = h is X k is, где z — единичный орт вдоль оси OZ ; h — соответствующий горизонтальный орт, лежащий в плоскости, нормальной к k; V — вертикальный орт, лежащий в той же нормальной плоскости. Отметим, что орт h лежит в плоскости, параллельной плоскости k Ok , и действительно явля-xy ется горизонтальным вектором. Впредь подразумеваем, что вектора h и V зависят от вектора k.
Пусть на рассеиватель падает плоская электромагнитная волна вида (1). Разложим вектор E 0i по системе векторов ( h i , v i ):
E 0i = (E 0 h h i + E 0 v Vi), а амплитуду рассеяния по направлению ks по системе векторов (h s, vs)
/X
T ( k s ) = ( T h ( k s ) h 5 + T v ( k s ) V s ).
Тогда, как показано в [ 25, 26 ] , для данного конкретного рассеивателя справедливо следующее матричное равенство:
Г T h ( k • > ]
T (k ) | v v s /у
(T T T 11 T 12
T T
I T 21 T 22

E ' s ( R , 0 , ф ) = T ' s ( 0 , ф)exp(i k Rl R
Здесь
T 's ( 0 , ф ) = T ( k 1 ),
T 's ( п - 0 , ф ) = T ( k 2 ),
1 + O | — I kR
. (11)
Матрица T ( k i , k s)
Г T„ T ,2 ]
T T
I T 21 T 22 у
называется мат-
рицей рассеяния (мр). Она не зависит от свойств
падающей плоской волны, а зависит только от
/X /X свойств рассеивателя и векторов k i и k s (отме-
0 e [ 0, п /2 ] , ф е [ 0,2 п ] , а e [ k ,0 ] ; а = k cos 0 ; k x = £ cos ф ; k y = ^ sin ф ; ^ = k sin 0 ; ^ = = ( k x , k y ) = ( 5 , ф ) — горизонтальная проекция волнового вектора. Функция T 's ( 0 , ф ) носит название амплитуды рассеяния, а ее значение совпадает с весовым множителем в (10) плоской волны, распространяющейся в данном направлении.
Фундаментальную роль в дальнейшем изложении играет матрица рассеяния (характеристика
тим, что речь идет о матрице рассеяния по электрической составляющей поля рассеяния, по магнитной составляющей мр вводится аналогично, и они аналитически связаны). В работе [ 25 ] приведена техника вычисления матриц рассеяния для конкретных рассеивателей.
Таким образом, мр напрямую ассоциируется с откликом рассеивателя на падение плоских волн. Поскольку выражение (10) двузначно зависит от переменных интегрирования k x , k y , обозначим матрицу рассеяния следующим образом: t m ( k i , k s ), n , m = 1,2 . Нижний индекс m = 1,2 ха-
рактеризует направление движения падающей плоской волны: m = 1 — волна движется в сторону возрастания z , (т.е. снизу вверх), m = 2 — наоборот. Индекс n характеризует направление рассеяния n = 1 — рассматривается рассеяние вверх (значение + а ), n = 2 — рассеяние вниз ( - а ). Рассматриваются области определения как внутри, так и вне области видимости.
Для введения необходимого в дальнейшем векторного оператора образуем следующий аффинор Am n ( k i , k s ) , ассоциированный с мр:
A m ( k i , k s ) =
ˆ nˆ n h s\(^11) m h i + (T12) m v i +
+ v s ( ( t 21 ) m • h i + ( t 22 ) m • V i ). (12)
Отметим, что все слагаемые зависят от векторов k i , k s ; ( Tlj ) n m — соответствующие элементы матрицы T m ( k i , k s), 1 , j = 1,2. После этого легко найти поле рассеяния при падении плоской волны с произвольным вектором поляризации E 0i :
E = f Am (ki,ks)Eoi exp(ik )dk dk s 2n R2 as s “ sy (13)
x e R 3 \ Q D .
Следует обратить внимание, что в (13) и далее имеет значение порядок следования аффинора и вектора вследствие невыполнения коммутативности их произведения в общем случае.
Когда центр рассеивателя смещается из начала координат в точку x 0 , вследствие линейности задачи отклик на падение плоской волны (1) имеет вид
E s = -i- exp(k i x o) x 2n xJ
R 2
A mn ( k i , k s ) E 0i
a s
exp( i k s( x - x 0))d k s x d k s y , (14)
x e R 3 \ Q D .
Пусть центр источника первичного падающего поля находится в начале координат. Тогда электрическое поле источника может быть представлено в виде, аналогичном полю рассеяния (10) [ 22, 23 ] :
E i = f D ( k m ) exP ( i k m x )d kx d k y ,
2 n R 2 a (15)
m = 1,2; x e R 3\ Q D .
Здесь D ( k m ) — вектор диаграммы направленности электрического поля, лежащий в нормальной плоскости к вектору k m . В случае когда рассеиватель с центром в точке x 0 находится в первичном поле (15), результирующий вектор поляризации рассеянного поля может быть найден с помощью техники, изложенной в [ 16, 17 ] , и имеет вид
T n ( k s , x o ) =
= г J An(ki,ks)Dm(ki) exp(ikixo)dkixdkiy, (16) 2n R 2 a i n = 1,2.
Здесь для устранения избыточных индексов принято обозначение D m ( k ) = D ( k m ) ( D ( k m ) фигурирует в (15)).
Знание результирующего вектора поляризации рассеянного поля позволяет находить поле рассеивателя в виде
E s = г J T n ( k s, x o) exp ( i k s ( x - x o))d k s x d k s y , 2 П R 2 a s (17)
n = 1,2.
Величину T n ( k s, x 0) из (16) при переменном x 0 можно трактовать как некое поле рассеяния, т. к. оно по форме совпадает с (13–15) и удовлетворяет векторному уравнению Гельмгольца. В случае электромагнитных волн, как и в скалярном случае, имеет место полезное разложение из [ 21 ] , которое справедливо в случае выполнения условий (4), (8). А именно любое решение векторного уравнения Гельмгольца вне наименьшей сферы, включающей источник, может быть представлено в виде абсолютно сходящегося ряда
E = exp( ikR ) у E i ( 0 , Ф )
R l ^o ( kR ) 1 ,
E i + 1 ( 0 , ф ) = L [ E i ( 0 , ф )]( 0 , ф, 1 ). (19)
Здесь L — оператор, определяемый выражением
дтж ф, 1 )
Л 9 , ф + 1 ( 1 + 1) 2 i ( 1 + 1)
W, ф ),
V
позволяет рекурсивно найти все члены ряда (18) по члену с индексом 0, который в свою очередь легко находится с помощью асимптотической оценки интегралов типа (10) и представляющих собой двумерные преобразования Фурье. Как видно из (11), они представляют собой соответствующее значение спектра, представленное в сферической системе координат, поделенное на exp( iaz I)
функцию ------!—. Тогда применительно, на- a пример, к (10) имеем
E о ( 9 , ф ) = T *s ( 9 , ф ). (21)
Окончательно для (16) имеем
T '( R о , 9 о , ф о , 9 s , ф s ) =
= ex p ( ikR о ) у A i ( 9 о , ф о , 9 s , Ф s) (22)
R о 1 = о ( kR о ) 1
A 1 + 1( 9 о , ф о , 9 s , ф s ) =
= L[ A z ( 9 о , ф о , 9 s , ф s )]( 9 о , ф о , 1 ),
A о( 9 о , ф о , 9 s , ф s ) =
= • (24)
= A m ( 9 о , ф о 9 , ф s ) D '( 9 о , ф о ).
Здесь приняты обозначения: ( R о, 9 о, ф о) — сферические координаты точки x о = ( x о, у о, z о), центра рассеивателя; ( R s, 9 s, ф s) — сферические координаты произвольной точки относительно точки-центра ( x о, у о, z о); Л 9ф — оператор Бельтрами на сфере (угловая часть лапласиана); штрих вверху функции говорит о том, что она рассматривается в соответствующей сферической системе координат.
Представление (22-24) позволяет легко проанализировать ошибки при вычислениях результирующего вектора поляризации в виде конечной суммы. Так, в зоне Фраунгофера можно ограни -читься нулевым приближением
T'(R о,9 о, фо,9 в, ф s) = n т т т \ eXp(ikR о) ,
= A m ( 9 о , ф о , 9 s , ф s) D ( 9 о , фо) D + (25)
R о
+ O
IT
что соответствует представлению о том, что поле источника в точке рассеивателя представляет собой плоскую волну с весом n т О т Cxp(ikR о)
Am ( 9 о , ф о , 9 s , ф s ) D ( 9 о , ф о ) D •
R о
Если подставить T из (25) в интеграл для поля рассеяния (17), то получаем следующую асимптотику:
E s ( R s , 9 s , ф s ) =
= exp(ikRо) exp(ikRs) у A1 (9s,Фs)
R о R s у ( kR s ) 1 ’
A1+1 (9s, ф s) = L [A 1 (9s, ф s )](9 s, Фs, 1),(27)
'
Aо(9s,фs) = Amn (9о,фо,9s,фs)D'(9о,фо).(28)
Если теперь уже в (26-28) ограничиться нулевым приближением, то получим лучевую версию рассеяния:
E s ( R s , 9 s , ф s ) =
'
= A mn ( 9 о , ф о , 9 s , ф s ) D '( 9 о , ф о ) x , (29)
exp(ikRо) exp(ikRs ) x, R о Rs а именно, луч, вышедший из направленного источника в сторону рассеивателя под углами 9о, фо в системе координат источника, взвешивается его диаграммой направленности D' (9о, фо), достигает рассеивателя в точке (Rо,9о,фо), отражается под углами 9s, фs относительно рассеивателя, одновременно взвешиваясь мр, и, достигая точки (Rs,9s,фs) в системе координат рассеивателя, приобретает амплитуду (29).
ДВА РАССЕИВАТЕЛЯ В ПОЛЕ ПЛОСКОЙ ВОЛНЫ
Рассмотрим далее случай, когда в однородном пространстве расположено два рассеивателя S 1 и S2 с соответствующими исходными мр (невозмущенными, когда их влияние друг на друга исключено), из которых строятся соответствующие аффиноры 1 Amn (k i, k s) и 2 Amn(k i, k s) (назовем их аффиноры рассеяния (ар)). Пусть центр первого z рассеивателя расположен в точке (о,о,——), вто- z рого в точке (о,о,-^-). Будем рассматривать сум- марное поле рассеяния в случае падения на указанную пару рассеивателей плоской волны E0i exp(ikix) , распространяющейся для определенности в сторону возрастания величины z , т. е. снизу вверх. Поле ищется суммированием многократно рассеянных плоских волн (принцип суперпозиции уместен, т. к. каждой падающей плоской волне соответствует своя совокупность рассеянных волн, которые вместе с падающей обеспечивают выполнение граничных условий на рассеивателе, и так для всех пар (падающая плоская вол-на—соответствующие ей рассеянные волны)). Отметим, кроме того, что построенное таким образом решение является единственным в силу единственности составляющих его слагаемых, начиная с первичных полей рассеяния, что определяется выбором соответствующих условий излучения Сильвера—Мюллера (а следовательно, и Зоммер-фельда).
Пусть результирующее поле каждого рассеивателя после всех переотражений определяется выражением (17), где под x 0 понимаются координаты центра соответствующего рассеивателя. Обозначим в (17) результирующие спектры векторов поляризации (рсвп) 1 T m n ( k i, k s) и 2 T m n ( k i, k s) соответственно. Рассмотрим, например, процесс формирования рсвп S 1 в верхней полусфере, т. е. 1 T 11( k i , k s) . Происходит следующее. Плоская волна E 0i exp( i k i x ) , распространившись мимо рассеивателей S 1 и S 2 , вызовет на них первичные поля рассеяния соответственно со спектрами векторов поляризации (свп)
1 A 1 n ( k i , k s ) E 0i exP ( - i a i 4 ° ) и
2Ain(ki,ks)Eoi exp(iai 40), которые и окажут влияние на формирование интересующей нас рсвп 1 T11(ki,ks) (множители exp(±iai Z°) отвечают фазе падающей плоской волны соответственно на S2 и S1 ). Поясним это на примере первичного рассеянного поля с свп 1 A11(ki,ks)E0i exp(-iai Z°). Это поле, порожденное плоской волной на S1 , распространяется вверх, достигает S2 , переотражается от него, идет далее вниз, переотражается от S1 , порождая первую добавку к рсвп, снова идет вверх и т. д. до бесконечности. В итоге формируется составляю- щая 11T11(ki,ks) рсвп 1 T11 (k i , k s ) , обусловленная первичным рассеянным полем на S1 с свп 1 A11(ki,ks)E0i exp(-iai 4°). Аналогично формируется составляющая 12T11(ki,ks) рсвп 1 T11(ki,ks), вызванная первичным рассеянным на S2 полем с свп 2 A12 (k i, k s) E 0i exp(iai Z°): первичное поле 2A12(ki,ks)E0i exp(iai -20-) на S2 распространится вниз, переотразится от S1 , вызвав свою первую добавку к 12T11(ki ,ks) , затем пойдет вверх, снова переотразится от S2 и т.д. до бесконечности. Таким образом, необходимо вычислить сумму
-
1 T/ ( k i , k s ) = 1 T 1( k i , k s ) + 2 T i1 ( k i , k s ). (30)
Используя выражение (16) для расчетов вторичных свп и действуя совершенно идентично тому, как это проделывалось для скалярного случая [16] или когда роль второго рассеивателя играла граница раздела [17], получим следующие ставления
1 1 T 1 1( k i , k s )
exp( - ia i z 0 /2)
= £ ( A1 ) n [ 1 A E 01 ] ( k i , k s ) = n = 0
= (l - A11 )—1 [1 AEoi ](ki,ks), пред-
-
2 T/( k i , k s ) =
= :£(A11)n [1 T11 ](k i, k s) = n=0
= ( l - A 1 1 )- 1 [ 1 T 1 1 ] ( k i , k s ), (32)
где 1 T ˆ 11( k i, k s) определяется выражением
-
1 T 11 ( k i , k s ) =
zi
= exp( ia i —) — x 2 2n
[1 A 2 1( k n , k s )][2 A 1 2( k i , k n ) E 0i ]
------------------------x xj
R 2
an x exp(ianzo)dknxdkny ,
а оператор A 11 определяется следующим образом:
^ 1 1[ T ]( R i , k s ) =
( i ) 2Г1 A 21 ( k m , k s ) f f 2 A i 2( k n , k m ) T ( k i , k n ) I
V П J R 2 a m V R 2 a n
exP ( i a n z 0)d k nx d k ny
X exp ( i«mZ 0)d k mx d k my .
Как видно из (34), оператор A 11 вычисляет добавку к рсвп 1 T 11( k i, k s) за один цикл: рассеянная на S 1 волна достигает S 2 , переотражается обратно, достигает S 1 и вызывает очередное приращение рсвп 1 T 11 ( k i , k s ) . Заметим, что учет конечного числа членов рядов (31), (32) может быть использовано для приближенного вычисления составляющих рсвп.
Формально ряды Неймана (31), (32) сходятся, когда норма оператора (34) меньше единицы. Однако из физических соображений ясно, что если соответствующие рсвп существуют, то по построению единственного решения представляют собой именно ряды типа (31), (32), т. к. эти ряды лишь отражают физику происходящих процессов и представляют собой решение, альтернативное описанному в работах [ 1-15 ] .
Для сходящихся рядов (31), (32) можно обратить оператор ( l - A 1 1 ) стандартным образом, после чего получаем интегральные уравнения для определения составляющих рсвп 1 T 11( k i , k s) :
( l — 41 ^T ] ( k i , k s ) =
= 1 A 1 1 ( k i , k s) E 0i exP ( - i a i y). (35)
( l - 41 )[ 2 T 11 ] ( k i , k s ) = 1 T/ ( k i , k s ), (36)
где оператор A 11 определяется из (34).
Отметим, что после определения величины рсвп 1 T 11( k i, k s) рассеивателя S 1, легко получить рсвп 2 T 11( k i, k s) рассеивателя S 2 . Для этого, как легко заметить, к первичной составляющей 2 A 11 ( k i , k s) E 0i exp( ia i “ у), вызванной падением плоской первичной волны, необходимо добавить составляющую, вызванную рассеянием на S 2 суммарного поля рассеяния рассеивателя S 1 с рсвп 1 T 11 ( k i , k s ) . Или в окончательном виде
2 T 11 ( k i , k s ) = 2 A /( k i , k s ) E oi exp( ia i Z 0) +
i 2 A 1 1( k n , k s )1 T 1 1( k i , k n )
+X
2 П R 2 a n
X exp ( i a n z o)d k nx d k ny . (37)
Поступая совершенно аналогично, найдем соответствующие интегральные уравнения для 21 T 12( k i, k s) и 22 T 12( k i, k s) , являющихся составляющими рсвп 2 T 12( k i, k s) рассеивателя S 2 в нижнем полупространстве:
2 T12(k i, k s)=21T12(k i, k s)+^(k i, ks),(38)
( l - A 22 )[ 2 1 T 1 2 ] ( k i , k s ) =
=2A12(ki, ks)E0i exp(iai Z0),(39)
(i - A22 )[22T12 ](ki,ks)=2T12(ki,ks),(40)
где 2 T ˆ 12 ( k i, k s) определяется выражением
2 T12(k i, ks) = zi
= exp( - ia i0) X
2 2 n
Xj
R 2
[2 A 1 2( k n , k s )][ 1 A 11 ( k i , k n ) E 0i ] X
a n
X exp ( i a n z o)d k nx d k ny ,
а оператор A 22[ T ]( k i, k s) равен:
A 2 2[ T ]( k i , k s ) =
(
Xj
V R 2
2 A 1 2( k m , k s )
X
T J R 2 a m
1 A 2 1( k n , k m ) T ( k i , k n )
a n
^ e xp ( i a n z 0)d k nx d k ny
J
X
X exp( i a m z 0 )d kmx d k my .
Рсвп рассеивателя S 1 в нижнем полупространстве определяется из выражения
1 T i 2( k i , k s ) =
=1 A12(ki, ks)Eoi exP(-iai zO) + i 1A12 (kn ,ks)2 T12(ki,kn)
+ X
2 n R 2 a n
X exp ( KZ^ o)d k nx d k ny . (42)
Пусть теперь первичная волна распространяется сверху вниз. Тогда останутся верными все выражения (30–42), если в них произвести следующие изменения. Первичные свп будут соответственно
1 An ( k i , k s ) E 0i exp( i a i z0-) и
-
2 A 2 n ( k i , k s ) E oi exp( - i a i z 2).
Далее в выражениях (30–42) во всех функциях свп, содержащих аргумент k i , необходимо правый нижний индекс в обозначении функции и аффинора поменять с 1 на 2, а во всех экспонентах, содержащих a i , поменять знак на противоположный. Все остальные обозначения и выражения останутся без изменений. Например, выражения (39), (40) примут соответственно вид
( l - A 22 ) [ 21 T 2 2 ] ( k i , k s ) =
= 2 A 22 ( k i , k s ) E 01 exP ( — i a i ^ Z 0) ,
( l — A 22 )[ 22 T 2 2 ] ( k i , k s ) = 2 T 22 ( k i , k s ).
Таким образом, получены все выра ж ения для определения рсвп двух рассеивателей 1 T m n ( k i, k s) и 2 T mn ( k i , k s), m , n = 1,2 . Такая информация является, безусловно, полезной хотя бы для того, чтобы судить о степени возмущения исходных свп 1 A mn ( k i , k s ) E 0i и 2 A mn ( k i , k s ) E 0i .
ОБОБЩЕНИЕ НА ПРОИЗВОЛЬНОЕ ЧИСЛО РАССЕИВАТЕЛЕЙ И ПРОИЗВОЛЬНОЕ
ПЕРВИЧНОЕ ПОЛЕ
В тех случаях, когда возмущением свп пренебречь нельзя, практически более значимой является, по-видимому, информация о совокупной свп системы, состоящей из рассматриваемых двух рассеивателей. Последнюю легко получить, пользуясь выражением (17) и принципом суперпозиции. Окончательно имеем для совокупной свп T mn ( k i , k s ):
T> i , k s ) =
= 1 T 1 ( k i , k s )exp( ia s Z 0 ) +
+ 2 Tm(k i, k s)exp(-ias yX Tm2(k i, k s) = m = 1,2 .
= 1 T m 2 ( k i , k s )exp( - ia s Z 0 ) + + 2 T m 2 ( k i , k s )exp( ia s Z 0 )
После вычисления свп уже можно проделать три важные вещи.
Во-первых, вычислить результирующую мр. В скалярном случае рассматриваемой задачи [16] вычисление скалярного спектра, аналогичного свп (43), уже означало вычисление совокупной амплитуды рассеяния, аналогичной мр для двух тел при падающей плоской волне единичной амплитуды и нулевой фазы. В векторном случае это не так. Анализ показывает, что в этом случае необходимо для определения четырех элементов результирующей мр дважды находить свп (43): отдельно для падающих плоских волн единичной амплитуды и нулевой фазы, поляризованных по hi и vˆ i векторам соответственно. Тогда элементы резуль- n тирующей мр T m (k i, k s) = nn
T m ,и T m ,12
г n t"
^ T m ,21 T m ,22
равны
n n n n
T m,11 = hihs T m (ki,ks), T m,21 = hivs T m (ki,ks), nn nn
T m ,12 = v i h s T m ( k i , k s ), T m ,22 = v i v s T m ( k i , k s ),
где, например, h i v s T nm ( k i, k s) означает v ˆs -составляющую свп (43), вычисленную при падении h i -поляризованной плоской волны единичной амплитуды и нулевой фазы, и т.д. Как показывает (17), центр совокупного рассеивателя, состоящего из описанных двух рассеивателей, находится посредине между центрами обоих рассеивателей, т. е. в данном случае в начале координат.
Во-вторых, умение вычислять результирую- щую мр Tm(ki ,ks) системы, состоящей из двух рассеивателей, позволяет последовательно решить задачу о N рассеивателях. А именно: сначала рассматривается система из двух рассеивателей, находится их результирующая мр, и оба рассеивателя заменяются одним, фиктивным, который обладает их результирующей мр с центром, лежа- щим посредине между центрами исходных рассеивателей. Затем рассматривается система рассеивателей, состоящая из третьего рассеивателя и описанного фиктивного, находится их результирующая мр и т.д., пока не будут исчерпаны все рассеиватели и не найдена единая результирующая мр всех N рассеивателей. Правомерность описанной итерационной процедуры следует из принципа суперпозиции.
В-третьих, после вычисления результирующей мр необходимого числа рассеивателей можно вычислить результирующий свп (16) всей системы рассеивателей при облучении системы рассеивателей направленным источником с полем (15). В (16) теперь уже будет фигурировать аффинор, ассоциированный с результирующей мр. Само поле рассеяния можно вычислить из выражения (17). При этом очевидно, что центр системы рассеивателей находится в их геометрическом центре.
МНОЖЕСТВЕННОЕ РАССЕЯНИЕ В СИСТЕМЕ: (ПЛОСКАЯ ГРАНИЦА РАЗДЕЛА— РАССЕИВАТЕЛЬ)
Представляет интерес получить из вышеприведенных выражений для системы двух рассеивателей результаты для случая, когда одним из рассеивателей является отражающая плоскость (плоская граница раздела). Для этого вначале найдем матрицу рассеяния границы раздела.
Пусть плоскость z = 0 служит границей раздела однородного верхнего полупространства z > 0 и нижнего слоисто-неоднородного полупространства z < 0 (идеальные границы являются частным случаем этой модели). При падении плоской волны E oi exp ( i ( k ix x + к у у - a l z )) из верхнего полупространства на границу раздела в нем образуется отраженная от границы плоская волна E s exp( i ( klxx + k ly y + a l z )), где E s — вектор поляризации отраженной волны. Для вычисления мр отражающей плоскости формально воспользуемся выражением (13). В этом случае имеем
E 5 exP ( i ( kix x + ку У + a i z )) =
= —Vn i ( a s ) a s S ( k i x - k sx ) S ( k i у - k s у ), i (46)
-
n , l = 1,2.
Здесь использовано фильтрующее свойство δ -функций. Функции Vnl (a i ) = Vnl (a s) имеют смысл соответствующих коэффициентов отражения для волн различной поляризации, рассмотренных во многих работах, например в [ 27 ] . Действительно, пусть единичный вектор падающей плоской волны горизонтально поляризован, тогда формально вектор E s из (45) равен:
E s = E H( a i ) h s + V 2i ( a i ) vs. (47)
В случае вертикальной поляризации имеем
E s = V ,2 ( a i ) h s + V 22 ( a i )vs. (48)
Впрочем, в случае плоскослоистой среды из физических соображений следует, что недиагональные члены мр должны быть равны нулю V 12 = V 21 = 0 и мр — диагональная.
Перейдем к постановке задачи. Пусть центр рассеивателя S 2 по-прежнему находится в точке z
(0,0, 20 ) , рассеиватель S1 отсутствует, однако присутствует граница раздела z = -z0 /2 с элементами мр вида (46). Первичное поле имеет вид E 01 exp( i ( k ixx + k 1уу - ai z)). Найдем поле рассеяния выше и ниже рассеивателя. Для этого в выражениях (30–42) во всех функциях свп, содержащих аргумент ki необходимо правый нижний индекс в обозначении функции и аффинора поменять с 1 на 2, а во всех экспонентах, содержащих ai, поменять знак на противоположный. Все остальные обозначения и выражения останутся без изменений. Роль S1 , разумеется, играет отражающая плоскость. С учетом этих замечаний и выражения (46) выпишем выражения для операторов (34) и (41). Для оператора A11 имеем i f A> i, k s)E oi x
2n j as x exp(i(ksxx + ksуУ + asz))dksxdksу .
A l [ T ]( k i , k s ) =
i
= — exp( ia s z 0 ) x
2 n
Анализ последнего выражения показывает, что элементы мр отражающей плоскости равны:
Tnl ,2 ( k l , k s ) =
= '^V ( aa ( k i x - k s x ) S ( k t у - k s у ) =
xj
R 2
(1 a 2 1( k s ))2 A 1 2( k n , k s ) T ( k i , k n )
x
an x exp(ianz0)dk„xdk„y ,
где 1 a 21( k s) — аффинор, равный
-
1 a> s) = hs Vn(as)h s + Vs V„(as)Vs.(50)
Для оператора A 22 имеем
A 2 [ T ]( k i , k s ) =
= f 2 A 1 2( k m , k s V a 2 ( k m ) T ( k i , k m ) x
2П R2
X exP(i2amz0)dkmxdkmy ,(51)
где 1 a 2 1 ( k m ) также определяется из (50) при k = = k m .
Теперь для определения свп поля рассеяния ниже рассеивателя S 2 необходимо воспользоваться видоизмененными интегральными уравнениями (39), (40):
( l - A 2 ) [ 21 T 22 ] ( k i , k s ) =
= 2 A 2( k i , k s ) E 0i exp( ia i z 0),
(l - A 2 )[22 T22 ](k i, k s) = 2 T22(k i, k s), где
-
2 T> i , k s ) =
= [2 A 2 ( k i , k s )][ 1 a 2 ( k i ) E 0i ] exp( ia i 32), а также (35), (36)
( l — A i 1 ’I T' ] ( k i , k s ) =
= A 2 1 ( k i , k s) E 0i e xp ( - i a i y), (52)
( l - A )[ 2 T 2 1 ] ( k i , k s ) = 1 T 2 * ( k i , k s ), (53)
где
-
1 rT 2 1( k i , k s ) =
= [ 1 a 2 1 ( k s )][2 A 2 ( k i , k s ) E 0i ] exp( ia i z 0).
Кроме того, необходимо видоизменить выражение (37), а именно
-
2 T 21 ( k i , k s ) = 2 A 2 1 ( k i , k s ) E 0i exp( ia i z 0) +
i 2 A 1 1( k n , k s )1 T 2 1( k i , k n )
+x
2 n R 2 a
-
x exp ( i a n z 0)d k nx d k ny .
Эти выражения позволяют вычислить рсвп 1 T 2 1( k i , k s ), 2 T 2 1( k i , k s ) и 2 T 2 2( k i , k s ). Однако в данном случае объединять рассеиватели в виде совокупной свп (43) нецелесообразно, а следует представить поле рассеяния в виде полей реального и мнимого рассеивателей. Центр мнимого рассеивателя, очевидно, находится в зеркальной от-3
носительно границы точке (0,0, - —z 0). Воспользовавшись выражением (17), получаем для поля рассеивателя
Es = —J s 2n 12
1 T 21( k i , k s)ex P ( ia s z 0/2)
+
-
2 T2 ( k i , k s )exp( - ia s z 0 /2) +
a
X exp(i(ksxx + ksуУ + as z))dksxdksy выше рассеивателя и
E s
=f
2 n J,
R
1 T 2 1 ( k i , k s)e xp ( ia s z 0/2)
x as
X exp(i(ksxx + ksуУ + asz))dksxdksу + f T2 (k i, ks)exp(ias z0/2) +X
1 2 a s
X e xp ( i ( k s x x + k s у У + a s z ))d k s x d k s у
ниже рассеивателя.
После этого выделятся составляющие поля реального и мнимого рассеивателя, как в скалярном случае [17]. Кроме того, использование выражений (18–20) позволяет получить геометрооптическое представление поля рассеяния реального и мнимого рассеивателей.
ГЕОМЕТРООПТИЧЕСКОЕ ПРЕДСТАВЛЕНИЕ РЕШЕНИЯ ЗАДАЧИ МНОЖЕСТВЕННОГО
РАССЕЯНИЯ
Рассмотрим далее случай, когда в поле первичной волны, образованной некоторым излучателем, по-прежнему находятся два рассеивателя. При этом ограничимся ситуацией, когда можно принять допущение, что рассеиватели находятся друг по отношению к другу и к излучателю в зоне Фраунгофера. Тогда можно воспользоваться выражениями (25), (29) для расчета свп рассеивателей и далее вычислить результирующие свп рассеивателей. Выражения (25), (29) при фиксированных расстояниях между объектами имеют порядок точности O(к -2), т. е. точность этих приближений растет с частотой.
Как всегда, отклик вначале будем искать на воздействие первичной плоской волны. Рассеяние при произвольной падающей волне может быть найдено по вышеописанной методике. Пусть плоская первичная волна E 0i exp( i k i x ), идущая для определенности снизу вверх, вызовет на рассеивателях S 1 и S 2 два невозмущенных рассеянных поля с свп соответственно
1 A '( 6 i , Ф i , 6 1s , Ф s ) E 0i ex P ( - i a i у) и
2 A '( 6 i , Ф i , 6 2s , Ф s ) E 01 ex P ( i a i y)
(свп представлены в сферических координатах; 6 1s, фs,6 2s, ф s — сферические углы относительно центров). Как показано выше, каждое из этих первичных рассеянных полей будет принимать участие в формировании результирующих рассеянных полей обоих рассеивателей. Получим приближение, например, для рсвп 1 T '( 6 i, ф i, 61s, ф s),
6 1s е [0,-]:
-
1 T '( 6 i , ф i , 6 1s , Ф s ) =
= 1 T '( 6 i , Ф i , 6 1s , Ф s ) + 2 T '( 6 i , Ф i , 6 1s , Ф s ).
Как нетрудно увидеть, вычисление последней рсвп распадается в этом случае на вычисление двух геометрических прогрессий:
mT'(6i,Фi,61s,Фs) ^(q)la0m , m = 1,2, l=0
где a01 =1 A'(6i,Фi,61s,Фs)E0i exp(-iai Z°);(58)
a 02 = P A '(61i = П 61s, Фs)]x x2A'E0iexp(iai z0) exp^ikzo);(59)
q =
2 e xp ( ikz 0 ) ) 2
V
x
z 0

[ 1 A '( 6 H = п , 6 1s , Ф s ) ] x 0, 6 2s = п ) ] .
Здесь q — аффинор, ассоциированный с матрицей, равной произведению матриц, ассоциированных с аффинорами 1 A ' и 2 A ' из (60) (произведение осуществляется в том же порядке, что и в (60)); в (59-60) учтено также, что при 6 = 0 или 6 = п зависимость от ф пропадает.
Полагая справедливым условие ||T=||< 1, где Г= — матрица, ассоциированная с аффинором q , q получим для рсвп
1 T '( 6 i , ф i , 6 1s , ф s ) = Q ' ( a 01 + a 02 ) , 6 1s е [0, П -].
Здесь Q ' — аффинор, ассоциированный с матрицей ( I - T q ) - 1, где I — единичная матрица. Равенство (61) следует из того факта, что ряду (57) соответствует матричный операторный ряд
Y = 2 TqX • где X , Y — вектор-столбцы. Тогда имеем
Y = ( I - T =) - 1 X .
q
Учитывая, что
X ' = (( a 01 + a 02 ) h , ,( a 01 + a 02 ) v/ ), Y ' =( * T h s , * T . ), имеем (61) (штрих означает здесь транспонирование).
П
Используя (37), получаем для 6 2s е [0,— ]
-
2 T '( 6 1 , ф i , 6 2s , ф s ) =
= 2 A ' ( 6 1 , Ф i , 6 2s , Ф s ) E 0i exp( i a i Z ° ) +
+ ^x^ [ 2 A '( 6 i = 0, 6 „, ф . )] x z 0
x [ * T '( 6 i , Ф i , 6 1s = 0)]. (62)
В выражениях (61), (62) фs е [0,2п]. Аналогично находятся результирующие свп в оставшейся п области определения 6s е [—, п], фs е [0,2п].
А именно:
-
2 T '( 6 1 , Ф 1 ,6 2s , Ф s ) = P '(а 01 + a 02 ), (63)
1 T '( 6 1 , ф i , 6 1s , Ф s ) =
= 1 A '( 6 i , ф i , 6 1s , ф s) E 0i exp ( - i a i 4 ° ) +
+ p z ) [ 1 1че,,= п, 6 „, ф s )] X z 0
X[2 T'(6 i, ф i,6 2s = П )], где a 01 = 2A ' (6i , Фi , 62s , фs ) exP(iai ^ ;
~ 02 =[ 2 A '( 6 1i = 0, 6 2s , Ф s) ] x
x["A '(6,, ф,, 6„ = 0)]x x E 0i exp(-ia i z0)exP*ikz0>;
P'— аффинор, ассоциированный с матрицей (I - Tp)-1. Матрица T= ассоциирована с аффино- ром p , равным:
p =
Г e^ ) 2
I z o J
[2 A '( 6 2i = 0, 6 2s , ф s )] X
X [ 1 A '( 6 ii = П , 6 s = 0)].
Выражения для случая, когда плоская первичная волна набегает на рассеиватели сверху вниз получаются из (58–67), если во всех экспонентах, содержащих a i , поменять знак на противоположный. Отметим, что при условии R s >> z 0 ( R s, 6 s, ф s— сферические координаты точки наблюдения относительно центра системы рассеивателей), в выражениях (58–67) можно положить 6 1s = 6 2s = 6 s .
Выражение для совокупной свп двух рассеивателей в геометрооптическом случае может быть получено из выражения (43), выписанного в сферической системе координат по отношению к началу координат (т.е. к средине отрезка, соединяющего центры рассеивателей):
T '( 6 i , ф i , 6 s , ф s ) =
= 1 T '( 6 i , ф i , 6 s , ф s )exp( ik cos( 6 s ) z ^) +
+ 2 T ' ( 6 i , ф i , 6 s , ф s ) exp( - ik cos( 6 s ) z ^). (68)
Аналогичный результат может быть получен стандартной процедурой перехода к предельным расстояниям. Пусть R— радиус-вектор точки на- блюдения относительно начала координат, а Rm— радиус-векторы, соединяющие центры рассеивателей с точкой наблюдения (m = 1,2); rm — радиус-векторы центров рассеивателей относительно начала координат. Тогда при R, Rm ^^ ( R = |R| , Rm = |R m| ) имеем usm- = (m T' (6i, фi, 6s, Фs) exp(-ikRrm ))x exp(ikR) (69)
X p); m = 1,2, R где u sm^ — предельное поле m -го рассеивателя;
R
R = R — единичный вектор в направлении от начала координат к точке наблюдения. Суммируя поля (69) и выделяя предэкспоненту, снова придем к выражению (68) для совокупной свп двух рассматриваемых рассеивателей.
Далее может быть применена процедура вычисления совокупной свп системы произвольного числа рассеивателей, описанная выше.
Список литературы Метод решения задачи множественного рассеяния электромагнитных волн в однородной безграничной среде
- Schwarzschild K.//Math. Ann. 1902. V. 55. P. 177.
- Hönl H. Maue A.W. and Westphal K. Handbuch der Physik/Ed. by S. Flügge. Berlin: Springer-Verlag, 1961. V. 25, N 1.
- Foldy L.L. The multiple Scattering of Waves//Phys. Rev. 1945. V. 67. P. 107-119.
- Lax M. Multiple Scattering of Waves//Rev. Modern Phys. 1951. V. 270, N 3. P. 287-310.
- Lax M. Multiple Scattering of Waves II. The effective field in dense systems//Phys. Rev. 1952. V. 85. P. 261-269.
- Tsang L., Kong J.A. Scattering of Electromagnetic Waves: Advanced Topics. New York: A Wiley-Interscience Publication, 2001. 405 p.
- Huang C.-C., Tsang L., Chan C.H. and Ding K. H. Multiple Scattering Among Vias in Planar Waveguides Using Preconditioned SMCG Method//IEEE Trans. on Microwave Theory and Techniques. 2004. V. 52, N 1. P. 20-28.
- Trinks W. Zur vielfachstreung an kleinen kugeln//Ann. Phys. Leipzig, 1935. V. 22. P. 561-590.
- Twersky V. Multiple scattering of electromagnetic waves by arbitrary configurations//J. Math. Phys. 1967. V. 8, N 3. P. 589-610.
- Bruning J.H., Lo Y.T. Multiple scattering of EM Waves by Spheres: Part I -Multipole Expansion and Ray-Optical Solutions//IEEE Trans. on Antennas and Propagation. 1971. V. Ap-19, N 3. P. 378-390.
- Narasimhan M.S., Ravishankar S. Multiple scattering of EM Waves by Dielectric Spheres Located in the Near Field of a Source of Radiation//IEEE Trans. on Antennas and Propagation. 1987. V. Ap-35, N 4. P. 399-405.
- Sebak A., Antar Y. Multiple scattering by parallel conducting elliptic cylinders: TE case//IEE Proc.-Microw. Antennas Propag. 1995. V. 142, N 2. P. 178-182.
- Waterman P.C. Symmetry, Unitarity and Geometry in Electromagnetic Scattering//Phys. Rev. 1971. D. V. 3, N 4. P. 825-839.
- Peterson B., Strom S. T Matrix formulation for Electromagnetic Scattering from an Arbitrary Number of Scatterers and Representation of E(3)//Phys. Rev. D. 1973. V. 8, N 10. P. 3661-3678.
- Varadan V.K. Multiple Scattering of Acoustic, Electromagnetic and Elastic Waves//Acoustic, Electromagnetic and Elastic Wave Scattering -Focus on the T-Matrics Approach. Pergamon Press, 1980. P. 103-134.
- Шарфарец Б.П. Решение задач множественного рассеяния акустической волны на произвольном количестве рассеивателей в однородном трехмерном пространстве//Научное приборостроение. 2004. Т. 14, № 3. C. 65-75.
- Зацерковный А.В., Сергеев В.С., Шарфарец Б.П. Использование амплитуды рассеяния для решения задач дифракции волн в полупространстве//Акуст. журн. 2001. Т. 47, № 5. С. 650-656.
- Шарфарец Б.П. Поле сферического излучателя звука в идеальном волноводе//Акуст. журн. 2002. Т. 48, № 4. С. 547-551.
- Шарфарец Б.П. Метод расчета поля излучателя и поля рассеяния неоднородного включения в плоскослоистых волноводах//Акуст. журн. 2004. Т. 50, № 1. С. 107-111.
- Colton D., Kress R. Inverse Acoustic and Electromagnetic Scattering Theory. New York: Springer, 1998. 346 p.
- Colton D., Kress R. Integral Equation Methods in Scattering Theory. New York: Wiley, 1983. 271 p.
- Шарфарец Б.П. Поле направленного излучателя в слоисто-неоднородном волноводе//Акуст. журн. 1985. Т. 31, № 1. С. 119-125.
- Hansen T.B., Yaghjian A.D. Plane-Wave Theory of Time-Domain Fields: Near-Field Scanning Applications. New York: IEEE Press, 1999. 378 p.
- Morse P.M., Feshbach H. Methods of Theoretical Physics, Part I. New York: McGraw-Hill, 1953. 997 p.
- Cho S.K. Electromagnetic Scattering. New York: Springer-Verlag, 1990. 389 p.
- Tsang L., Kong J.A., Ding K.H. Scattering of Electromagnetic Waves, Vol. 1: Theory and Applications. Wiley Interscience, 2000. 426 p.
- Бреховских Л.М. Волны в слоистых средах. М.: Наука, 1973. 344 c.