Проводящий тонкий экран в волноводной структуре: задача дифракции и задача трансмиссии

УДК 517.958

И. Л. Александрова, Н. Б. Плещинский

ПРОВОДЯЩИЙ ТОНКИЙ ЭКРАН В ВОЛНОВОДНОЙ СТРУКТУРЕ: ЗАДАЧА ДИФРАКЦИИ И ЗАДАЧА ТРАНСМИССИИ1

Аннотация. Задача дифракции и задача трансмиссии электромагнитной волны на идеально проводящем тонком экране в плоском волноводе сведены к парному сумматорному функциональному уравнению относительно коэффициентов разложения искомой волны в ряд Фурье. Парное уравнение равносильно бесконечной системе линейных алгебраических уравнений, которая решается методом усечения. Предложена вычислительная схема для решения задачи трансмиссии в случае, когда экран находится в открытом пространстве.

Ключевые слова: электромагнитные волны, задача дифракции, задача трансмиссии, плоский волновод, волны в открытом пространстве, интегральные уравнения.

I. L. Aleksandrova, N. B. Pleshchinskiy

THIN CONDUCTING SCREEN IN A WAVEGUIDE STRUCTURE: THE PROBLEM OF DIFFRACTION AND THE PROBLEM OF TRANSMISSION

Abstract. The diffraction problem and the transmission problem for an electromagnetic wave on a thin conducting screen in a plane waveguide are reduced to a paired summatory functional equation with respect to coefficients of the Fourier expansion. The paired equation is equivalent to an infinite set of linear equations. This set is solved by the truncation method. The authors obtain a computational scheme for solving the transmission problem on the screen in the open space.

Key words: electromagnetic waves, diffraction problem, transmission problem, plane waveguide, waves in open space, integral equations.

1. Постановка задач

Как правило, любая электромагнитная волна в волноводной структуре может быть представлена в виде суммы двух волн противоположной ориентации. Будем говорить, что в случае среды без потерь волна ориентирована положительно (прямая волна), если она переносит энергию или затухает в направлении оси волновода. Волна отрицательной ориентации (обратная волна) переносит энергию или затухает в противоположном направлении. Пусть в волноводе имеется неоднородность конечного размера, например, диэлектрическая вставка или проводящее тело. Если на неоднородность набегают приходящие с бесконечности электромагнитные волны, причем в общем случае с двух сторон, то появляется дифрагированное поле - волны, уходящие от неоднородности на бесконечность. Таким образом, в общем случае в процессе дифракции участвуют четыре волны: две волны, положительно и отрицательно ориентированные, с одной стороны от неоднородности и две волны, также разной ориентации, с другой стороны.

1 Работа выполнена при поддержке РФФИ, проект 12-01-97012-р_поволжье_а.

Задача дифракции электромагнитных волн на неоднородности в волноводной структуре ставится так: по набегающим на неоднородность волнам нужно найти уходящие от нее волны. Назовем задачей трансмиссии следующую задачу: по электромагнитному полю с одной стороны от неоднородности нужно найти поле с другой ее стороны (лат. 1хап8т188Ю - передача, переход).

Можно рассматривать две равноправные задачи трансмиссии: или по полю с одной стороны от неоднородности ищется поле с другой стороны, или наоборот. У задачи дифракции, вообще говоря, также есть пара - задача, когда по уходящему от неоднородности полю нужно восстановить набегающие на нее волны. Отметим, что и задача дифракции, и задача трансмиссии ставятся как задачи сопряжения для уравнений электродинамики, но имеется существенное отличие: в первом случае нужно найти их решения в двух областях, удовлетворяющие некоторым условиям на общей границе, а во втором - в одной и той же области, но решения разной ориентации.

Задачи дифракции электромагнитных волн на идеально проводящих тонких экранах относятся к классическим задачам электродинамики (см., например, [1], гл. IV, § 2). Близкими к ним по постановке и по методам исследования являются задачи дифракции на отверстиях в экранах, а также так называемые задачи «о связи объемов через отверстие». Первое (и фактически единственное) явное решение задачи дифракции на идеально проводящей полуплоскости получил А. Зоммерфельд. Морс и Рубинштейн [2] построили решение задачи дифракции на ленте в виде ряда по функциям Матье как предельный случай решения задачи дифракции на проводящем эллиптическом цилиндре. Приближенное решение задачи дифракции на ленте и щели, основанное на функциональном уравнении Винера - Хопфа, дано в книге Б. Нобла [3]. Монография А. С. Ильинского и Ю. Г. Смирнова [4] посвящена аналитическому исследованию векторных электродинамических задач на незамкнутых проводящих тонких поверхностях. Теория разрешимости соответствующих краевых задач построена методами теории псевдодифференциальных операторов.

В данной работе исследуется случай, когда неоднородность представляет собой идеально проводящую бесконечно тонкую пластину (экран), расположенную на поперечном сечении волноводной структуры. Будем обозначать через M экран, на котором Wo (x) = -м0 (x), и через N - оставшуюся часть сечения (дополнение для M), где Wl (x) = 0 . Как обычно, предельные значения (следы) компонент электромагнитного поля: Wo(x) = -щ(x) на M; Wl(x) = 0 на N, должны удовлетворять условиям: на M касательные составляющие электрического вектора обращаются в нуль, а на N касательные составляющие электрического и магнитного вектора непрерывны.

Рассматриваются две простые волноводные структуры: плоский волновод и открытое пространство (рис. 1). Предполагается, что составляющие гармонически зависящего от времени электромагнитного поля не зависят от поперечной координаты у декартовой системы координат. В этом случае, как известно, компоненты ТЕ- или ТМ-поляризованных волн выражаются через потенциальную функцию - решение двумерного уравнения Гельмгольца. Следовательно, следы компонент поля на сечении волноводов выражаются через следы потенциальной функции и ее нормальной (по отношению к сечению) производной. В дальнейшем эти следы будем помечать нижними

индексами 0 и 1, будем называть их нулевым и первым следами волны соответственно.

a \ x

u JV M . V +

0

Рис. 1. Проводящий тонкий экран в плоском волноводе и в открытом пространстве

Ограничимся случаем, когда электромагнитное поле ТЕ-поляризовано. Также предположим, что волноводные структуры слева и справа от экрана заполнены одной и той же однородной изотропной средой без потерь.

2. Электромагнитные волны в плоском волноводе

В случае ТЕ-поляризации поля потенциальная функция и(х, г) -решение уравнения Гельмгольца

Э2и Э2и ,2 / ч ^ /14

—2 + —2 + к и

Эх Эг

должна обращаться в нуль при х = 0 и при х = а. Так как функции

. . 2 . ппх 1 _

^ (х) = л-81П-, п = 1,2,...,

V а а

образуют полную ортонормированную систему функций на отрезке [0, а], то любое решение однородной краевой задачи для уравнения (1) в полосе имеет вид

I - \

и

п=1

где йп и йп - некоторые постоянные; уп - такие числа, что

у n +

=k2.

V a J

Условимся, что Reyn >0 или Imyn <0. Тогда при зависимости от времени вида eimt слагаемые с коэффициентами ип соответствуют волнам положительной ориентации, а с коэффициентами ип - волнам отрицательной ориентации. Ситуация, когда уп =0 (волны нулевой ориентации) здесь не рассматривается.

Пусть Ид(х) = и (х,0), ^(х) =—(х,0) - следы волны на сечении г = 0

Эг

волновода. Будем обозначать идп и щ п коэффициенты Фурье этих следов относительно системы функций ,?п (х). Легко получить следующее утверждение.

Лемма 1. Электромагнитная волна в плоском волноводе положительно ориентирована тогда и только тогда, когда

и\,п +1У п ид,п =0, п = 1,2,..,

или когда

а а

щ

0 0

где

1

К^, х) = _ ^ У^п

п=1 п=1 У п

Заметим, что равенства (4) представляют собой нелокальные граничные условия, связывающие друг с другом нулевой и первый следы волны в плоском волноводе. Их можно рассматривать как пару обращающих друг друга интегральных уравнений.

Ядро К_1^, х) имеет логарифмическую особенность на диагонали области определения, а ядро К^, х) - гиперсингулярное. Сходимость рядов в формулах (2) и (5) в общем случае, так же как и некоторые равенства, следует понимать в обобщенном смысле.

Для отрицательно ориентированных волн в формулах (3)-(5) нужно поставить другой знак перед постоянными уп .

Лемма 2. Ориентированная волна в плоском волноводе однозначно определяется по одному своему следу на сечении волновода. Неориентированная волна однозначно определяется по двум своим следам.

Действительно, если на сечении г = 0 однородного плоского волновода заданы функции (в общем случае - распределения) ^ (х) и щ (х), то при любом п коэффициенты представления (2) находятся из системы уравнений ип + ип = и0п, _ iуnun + iупип = щ п . Для ориентированной волны для

каждого п остается только одно уравнение.

Утверждение леммы остается в силе для любого сечения плоского волновода, а также для полубесконечного волновода, поскольку потенциальная функция поля в полуполосе естественным образом продолжается на всю полосу.

3. Плоский волновод. Сумматорные уравнения

Пусть на экран в плоском волноводе набегает слева электромагнитная волна с потенциальной функцией и (х, г). При ее дифракции появляются волны, отраженные влево и прошедшие вправо; их потенциальные функции обозначим и (х, г) и V (х, г).

Будем искать V(х, z) = и(х, г) + w(х, г), где w(х, г) - возмущение от экрана потенциальной функции. Тогда граничные условия и условия сопряжения на сечении волновода имеют вид

йО (x) + йО (x) = О, йО (x) + wo(x) = О на M ; йО (x) = Wo (x), й1 (x) = W1 (x) на N.

Теорема 1. Задача дифракции электромагнитной волны на экране в плоском волноводе эквивалентна бесконечной системе линейных алгебраических уравнений

wk ^ 2 ' уnwn 2 ' In,m Jm,k = 2', unIn,k, k = 1,2,^,

где

n=1

m=1 У m

n=1

Isn (t)sm (t) dt, Jm,k = Ism (t)sk (t) dt.

M

N

Доказательство. Из равенств

2 М>^п

п=1 п=1 п=1

Преобразуем ПСФУ к БСЛАУ с помощью интегрально-сумматорного тождества

а Ґ Л ^ 1 Л

2 Vп (х) = {I 2 Уп^п (0 2-----*т (((х) Л,

п=1 0 V п=1 /V т=1 Ут /

Из тождества (10) следует, что на N

V

1

2 wnsn

п=1 М V п=1 /V т=1 У т

xЄ [О,а].

M V n=1

1

2', уnwn 2', sm (x)In

n=1 m=1 Уm

Спроектируем новое парное уравнение (первое равенство из (9) и равенство (11)) на функции Sk (х) и получим БСЛАУ (7). #

При решении задачи дифракции по ип из БСЛАУ нужно найти wn , тогда Уп = ип + wn и ип = wn, п = 1,2,... Отметим, что интегралы (8) вычисляются аналитически, причем все рассуждения легко перенести на ситуацию, когда М - не один экран на сечении волновода, а несколько таких экранов.

А

Аналогичным образом рассматривается задача дифракции в случае, когда волна от внешнего источника набегает на экран справа или когда заданы две волны с потенциальными функциями и (х, г) и V(х, г). При этом в БСЛАУ (7) изменятся правые части уравнений, но не матрица коэффициентов. Поэтому разделять общую задачу дифракции на две подзадачи не имеет особого смысла.

Задача трансмиссии в общей постановке также может быть разделена на две подзадачи. Действительно, если заданы потенциальные функции

V (х, г) и V(х, г), то по первой из них можно найти решение задачи дифракции - волну и'(х, г), прошедшую влево, и волну х, г), отраженную вправо. Затем по разности V(х, г) - у'(х, г) нужно восстановить приходящую справа волну и (х, г) и второе слагаемое и"(х, г) в потенциальной функции волны, уходящей влево. Такое упрощение задачи трансмиссии целесообразно.

Итак, при решении упрощенной задачи трансмиссии по Vn нужно найти ип, но для этого можно использовать ту же самую БСЛАУ (7). Достаточно заменить в правой части коэффициенты ип на Vn — wn . Отметим, что числа Vn не могут быть заданы произвольно. Должно быть выполнено условие 2 Vnsn (х) = 0 на М .

п=1

Теорема 2. Задача трансмиссии электромагнитных волн на экране в плоском волноводе эквивалентна БСЛАУ

^ 21Уп^п 21 ^.т^т,] 21 In,kwn = 21 ^п-п,], ] = 1,2,---

п=1 т=1Ут п=1 п=1

4. Плоский волновод. Вычислительный эксперимент

Приближенное решение БСЛАУ (7) может быть получено методом усечения. В конечномерной СЛАУ

К Ь 1 к

—^ 2 1уп ^ 2', -п,т ^т,] = 2 2-п.] un, ] = 1.. К,

п=1 т=1Ут п=1

присутствуют два параметра усечения К и Ь .

Задача дифракции. Численное решение СЛАУ (13) легко находится, например, методом Гаусса.

Можно искать одновременно и вектор w = ^,..., wк), и вектор

V = (1,..,^), если удвоить число неизвестных и число уравнений, т.е. добавить к соотношениям (12) уравнения

—^ + Ч = Ч. ] = 1К.

Задача трансмиссии. Если заданы значения Vn и нужно определить wn и ып , то проще всего переписать уравнения (13), (14) следующим образом:

К Ь 1 к

2', у п^п 21 ^.т^т,] 21 ^п.^п =0, и] = .

п=1 т =1Ут п=1

Если в качестве исходных данных в задаче трансмиссии использовать приближенное решение задачи дифракции, то вектор и может быть восстановлен с достаточной точностью при относительно малых размерах экрана. При увеличении экрана результат становится хуже, но можно компенсировать потерю точности за счет увеличения параметра Ь . Аналогичный эффект наблюдается, если вместо (15) рассматривать СЛАУ вдвое меньшей размерности

К Ь 1 к к

21Уп ^ 21 —п,т ^т,] 21 -п,] ^ = 21 -п,] уп , ] = 1.. К,

п=1 т=1Ут п=1 п=1

которая является усечением БСЛАУ (12) или выводится непосредственно из (13).

Вычислительный эксперимент показал, что решение задачи трансмиссии устойчиво по отношению к достаточно малым возмущениям исходных данных.

Приведем простые рассуждения, уточняющие связь между задачей дифракции и задачей трансмиссии. Если задача дифракции (прямая задача) сводится к операторному уравнению вида Л^! = Ви , то при решении задачи трансмиссии (обратной задачи) нужно искать вектор w так, чтобы получилось w + и = V. Тогда новое уравнение для определения w имеет вид (Л + В)w = Ву . Следовательно, задача трансмиссии может быть решена численно примерно с той же точностью, что и задача дифракции, если при возмущении оператора Л оператором В сохраняются его свойства.

Стабилизировать счет при решении задачи трансмиссии можно и с помощью следующего приема. Если известно, что у искомого вектора и последние компоненты должны быть равны нулю: ып =0 при п > К1 (здесь К1 < К), то эти равенства можно использовать при модификации алгоритма. Нужно взять К1 первых уравнений из СЛАУ (16) и добавить к ним уравнения Wk = V], ] = К1 + 1.. К .

Решение задачи трансмиссии можно получить также методом, который часто используется при решении обратных задач. Выберем множество линейно независимых векторов и3, у = 1.. К и найдем соответствующие им решения V3 задачи дифракции. Если векторы V-3 тоже линейно независимы, то любой заданный вектор V можно разложить по этому базису. Если

1 К 1 К

V = + ... (ХкУ , то решение задачи трансмиссии и = +... а ки . Вы-

числительный эксперимент показал, что в этом случае результат тем лучше, чем больше экран.

5. Волны в открытом пространстве

Открытое пространство (плоскость в двумерном случае) также является волноводной структурой. Потенциальными функциями собственных волн являются решения уравнения Гельмгольца (1) вида и(х,г) = е^^^, если

^2 + С2 = ]2 (в общем случае £, и С - комплексные числа).

Сузим класс допустимых решений уравнения Гельмгольца. Пусть при каждом фиксированном значении г дважды дифференцируемая функция и(х,г) как функция аргумента х имеет медленный рост на бесконечности (не

может возрастать быстрее, чем полином). Тогда в выражении и( х, г) = е*&+*& параметр Ъ может принимать только вещественные значения. Кроме того, у неограниченной во всех направлениях волноводной структуры как бы появляются стенки, ось г становится ее осью, а прямые г = с - поперечными сечениями.

Будем рассматривать и(х,г) при фиксированном г как распределение медленного роста на бесконечности. Применим интегральное преобразование Фурье по переменной х и получим из уравнения Гельмгольца диффе-

д2и

ренциальное уравнение —— + (]2 — Ъ2)и (Ъ, г) = 0 (обобщенная производная

дг 2

совпадает с классической). При общее решение этого уравнения

и(Ъ, г) = С1©е—(Ъ)г + С2(^)егуг, где с^), С2(Ъ) - распределения медлен/-------------------------- д^и

ного роста и у(Ъ) = V ]2 — Ъ2 Ф 0 . При Ъ = ±] имеем —— = 0 и, следователь-

дг

но, и (Ъ, г) = С1(^) + С2(^)г .

Условимся, что или Яеу(Ъ) > 0 , или 1ту(Ъ) < 0 при ЪФ±] . Тогда элементарная волна е ^ е (Ъ)2 имеет положительную ориентацию по отноше-

—гЪх /у(£)г т>

нию к оси г, а волна е ^ е^’ - отрицательную ориентацию. В случае

Ъ = ±] будем говорить, что элементарные волны имеют нулевую ориентацию.

Таким образом,

и

-^

здесь распределение С1(Ъ) определяет волны положительной ориентации, а распределение С2(Ъ) - волны отрицательной ориентации (по отношению к оси г или к поперечному сечению волноводной структуры г = 0 ). «Лишнее» решение (при Ъ = ±] ) не зависит от г, и(х, г) = С1е—1] + С2е1]. Такое решение нулевой ориентации не переносит энергию вдоль оси г и не затухает.

Утверждения о следах потенциальных функций на прямой г = 0 такие же, как и в случае плоского волновода. Отличие в том, что они формулируются на языке образов Фурье следов волны.

Лемма 3. Электромагнитная волна в открытой плоскости положительно ориентирована тогда и только тогда, когда

«1

или когда

—^

где

—^

—^

Лемма 4. Ориентированная волна в открытой плоскости однозначно определяется по одному своему следу на прямой г = 0. Неориентированная волна однозначно определяется по двум своим следам.

Действительно, положительно ориентированная волна восстанавливается по образу Фурье нулевого следа так:

Для волн отрицательной ориентации нужно изменить во всех формулах знак у функции у © .

Если же заданы оба следа ио (х), щ (х), то имеем СЛАУ для образов Фурье следов слагаемых положительной и отрицательной ориентации

Открытая плоскость отличается от плоского волновода тем, что каждую волну в такой волноводной структуре однозначно определяет некоторое распределение (образ Фурье следа потенциальной функции), а не дискретный набор коэффициентов разложения по полной системе собственных волн (коэффициенты Фурье).

Задача дифракции электромагнитной волны на проводящем тонком экране в открытом пространстве может быть сведена к интегральным уравнениям различного вида. Удобно использовать связь между следами искомого поля на поперечном сечении волноводной структуры. Наиболее просто получить интегральное уравнение 1-го рода с логарифмической особенностью в ядре.

Пусть, как и в случае плоского волновода, на экран набегает слева электромагнитная волна с потенциальной функцией и (х, г). Будем искать потенциальные функции и (х, г) и V (х, г) волны, отраженной влево, и волны, прошедшей вправо, причем так, что V(х, г) = и(х, г) + w(х, г), где w(х, г) -

возмущение от экрана. Как и в случае плоского волновода, должны быть выполнены условия

—^

uo © + uo © = uo (£Х ui © + ui © = ui (£X Ui (£) + іy(£)«o (£) = o, Ui © — iy(£)Mo (£) = o.

6. Экран в открытом пространстве. Интегральные уравнения 1-го рода

й0 (x) + Uo (x) = 0, Uo (x) + Wo(x) = 0 на M ; Uo (x) = Wo (x), u (x) = Wi (x) на N.

Теорема 3. Задача дифракции электромагнитной волны на экране в открытой плоскости эквивалентна интегральному уравнению относительно первого следа возмущения от экрана

Доказательство. Из условий (20) следует, что м0(х) = ^0(х) и на N, и на М . Тогда из равенства (18) и его аналога для волн отрицательной ориентации следует, что -М1 (х) = ^1 (х). Поэтому ^1 (х) = Й1(х) = 0 на N, и второе из условий (19) принимает вид

М

Таким образом, для определения вспомогательной функции w(х, г) в полуплоскости г >0 имеем смешанные граничные условия

Эта задача, легко видеть, сводится к интегральному уравнению (21). Заметим, что это уравнение - парное, его следует рассматривать как уравнение 3-го рода. #

При численном решении интегрального уравнения (21) удобно использовать метод Галеркина с разложением искомой функции по полиномам Чебышева 1-го рода с весом, если ширина ленты не более, чем длина волны. Этот метод используется достаточно часто (см., например, [4, 5]), поэтому на деталях мы не останавливаемся. Отметим только, что если М = [-1,1] (этого всегда можно добиться заменой переменных), то

Вместо интегрального уравнения (21) можно рассматривать интегральное уравнение относительно образа Фурье Wl (£) первого следа Wl (х). Действительно, непосредственно из равенства (18) следует, что

Его решение должно быть таким, чтобы было выполнено и второе условие из (23).

M

Wo (x) = —Uo (x) на M; Wi (x) = 0 на N.

—^

Тогда первое условие из (23) сводится к уравнению

—^

доопределенные нулем вне отрезка [-1,1], переходят в функции Бесселя, а точнее, в функции

Поэтому приближенное решение интегрального уравнения (24) целесообразно искать в виде

Как показал вычислительный эксперимент, при использовании метода Галеркина в этом случае получаются те же самые значения коэффициентов а^ .

Построим вычислительную схему для решения задачи трансмиссии в случае, когда экран расположен в открытом пространстве. Предположим, что образы Фурье следов волны, набегающей слева на экран, имеют ограниченный носитель, расположенный в пределах отрезка [-А,А]. Тогда

В такой волне отсутствует часть элементарных гармоник, затухающих в направлении оси г (возможно даже, что затухающие слагаемые вообще отсутствуют).

Выберем на [-А, А] конечную систему линейно независимых функций Ф7 (^), 7 = 1.. Ь . Для каждой из них найдем

и построим приближенное решение интегрального уравнения (21), т.е. вектор

ствия между полиномами Чебышева и функциями Бесселя при преобразова-

k=o

а7 =(а0,..,аК), составленный из коэффициентов разложения функции w^j (^) по полиномам Чебышева с весом. В силу указанного выше соответ-

'0 ’•••’ UK

нии Фурье легко выписать коэффициенты разложения образа Фурье w^j (^) по функциям ф7 (^), 7 = 1.. Ь . Так как w1 (£) - /у(£)Wj0) (£) = 0, то коэффициенты разложения по тем же функциям образа Фурье нулевого следа Wj0) (^) равны

Ъп = (/)* А -1/*©Ф"©^.

',2 к=1 -А у

Будем рассуждать на языке коэффициентов разложения по базисным функциям образов Фурье нулевых следов волн, участвующих в волновом процессе. Векторы из коэффициентов разложения прошедших вправо волн

у0 (^) получаются как суммы единичных векторов е7 , соответствующих

следам М(0(^) = ф7(^) набегающих слева волн, и векторов Ъ7 = (Ъ^,..,ЪК).

При этом векторы е7 + Ъ будут линейно независимыми.

Если задан нулевой след Уэ(^) волны справа от экрана, то соответствующий ему вектор коэффициентов разложения по функциям ф7 (^) запишем как линейную комбинацию векторов е-1 + Ъ с некоторыми коэффициентами а7 . Тогда нулевой след искомого решения задачи трансмиссии имеет

вид

Ь

2 а 7 Ф7 (£).

7=1

Если в качестве исходных данных в задаче трансмиссии использовать приближенное решение задачи дифракции, то коэффициенты разложения функции ^1(х) в ряд по полиномам Чебышева могут быть восстановлены

с достаточной точностью. Возмущение параметра А не влияет на результат. Данный метод устойчив к небольшим возмущениям исходных данных.

7. Интегральные уравнения 2-го рода

Задача дифракции электромагнитной волны на экране в открытом пространстве может быть сведена также и к интегральным уравнениям 2-го рода.

Рассмотрим связи вида (19) между нулевым и первым следами потенциальной функции w( х, г) с учетом граничных условий (23) - равенства (22) и

Wl

N М

Чтобы перейти к интегральным уравнениям 2-го рода, нужно исключить из формул (22) и (25) одну из искомых функций (см. [6], а также [7], гл. 6).

В случае плоского волновода можно получить интегральное уравнение точно того же вида, если рассуждать не на языке коэффициентов Фурье следов потенциальных функций, а на языке самих этих функций. Действительно, ПСФУ (9) равносильно смешанным граничным условиям (23). Интегрально-сумматорное тождество (10) превращается во второе равенство (4), причем интегрирование в этом равенстве проводится только по М в силу условия Wl (х) = 0 на N . Если выражение для Wl (х) из первого равенства (4) подставим в равенство (22), то получим

wo (x) = I I wq (т)Ki (т, t) dт

m К o

K— i(t, x) dt

при х є N . С другой стороны, Wo (х) = -м0 (х) при х є М . Если это парное уравнение спроектировать на функции Sfc (х), то придем к БСЛАУ (7).

В случае плоского волновода регулярные БСЛАУ, эквивалентные задачам дифракции и трансмиссии, были равносильны интегральным уравнениям 2-го рода, заданным на всем сечении волноводной структуры. Покажем, как для открытой плоскости можно свести задачу дифракции электромагнитной волны на экране к регулярному интегральному уравнению на оси. Для этого нужно принять во внимание, что в периодическом случае коэффициенты Фурье функции заменяют ее образ Фурье, а сумма ряда Фурье - аналог обратного преобразования Фурье.

Как было установлено выше, возмущение w

Wo (х) = -м0 (х) на М; Wl (х) = 0 на N,

при этом для положительно ориентированной волны (см. вторую формулу (19)) с учетом второго граничного условия

Wo( х) = | Wl(^) К-і(ґ, х) & = | Wl(t) К-і(ґ, х) &.

—тс М

Функции под знаком интеграла заменим на обратные преобразования Фурье, и тогда на N

M

\

w0(x) = I jb I^(S )w0 Й1)<Г'¥ dSi I ^ (t—x) d^2

2_ 1 ^"'o^ І I

V2_

2_ J Y

—S2 xI (S2 —Si) d S2

d Si,

здесь и дальше

I© = | ег^ сИ, 3(£) = |Ж.

М N

Теперь перейдем к образу Фурье искомой функции. Имеем

1 1 1

^о(^)= ,---- | Wo(х)еХ Сх = ---- |Wo(х)еХ Сх + ---- |Wo(х)ех

i

л/2_ 0V~'W V2_M'r' V2_N"dX —2_

X

XIUq (x)eixS dx + —2 I Y(Si )"0 (Si) I “TS^I(S2 — Si) J(S — S2) dS2

M — ^ _ —^

d Si.

Известия высших учебных заведений. Поволжский регион

Выразим функцию м0 (х) через ее образ Фурье м0 (£) и получим

ё^ = /

W

Y

Se (-~, +~),

/(S) = 2- j«o(Si)I(S-Si)dSi.

2л

—^

Так как J(Е, — £2) = 2л5(^2 — £) — ^(^ — ^2), то это интегральное уравнение можно переписать в виде

—^0© + \^1)К(^аё^1= /©, ^ (—+-),

к

4л2

2л-

i г i

^S)1 (S-Si)- f -(FT1 (S2 -Si) 1 (S-S2) dS2

S) J y

y

Можно перейти от интегрального уравнения к БСЛАУ, если подобрать систему функций (распределений) на оси, по которым можно разлагать в ряды следы потенциальных функций всех участвующих в процессе дифракции волн. Как известно, любое распределение умеренного роста на бесконечности и его образ Фурье разлагаются в ряды по функциям Эрмита hj (х), причем

между коэффициентами разложений имеется простая связь.

Если

uo (S) = 2 ujhj (S), w0 (S) = 2 wjhj (S),

j=0 j=0

то

/ (S) = 2 Uj 2Л j hj (Si) I (S-Si) d Si.

1=0

Умножим обе части интегрального уравнения на hk (^), проинтегрируем по оси и получим БСЛАУ

—Рл + 2 ™]Кк] = 2и/, к=0д- • -1=0 1=0

где

Kkj = j jhk(S)hj(Si)K(Si,S)dSidS,

/kj =2л j jhk(S)hj(Si)I(S-Si)dSidS.

—^^

Список литературы

1. Хенл, Х. Теория дифракции / Х. Хенл, А. Мауэ, К. Вестпфаль. - М. : Мир, i964. - 428 с.

2. Morse, P. M. The diffraction of waves by ribbons and by slits / P. M. Morse, P. J. Rubinstein // Phys. Rev. - i938. - V. 54, № ii. - P. 895.

3. Нобл, Б. Метод Винера - Хопфа / Б. Нобл. - М. : ИЛ, i962. - 280 с.

4. Ильинский, А. С. Дифракция электромагнитных волн на проводящих тонких экранах (Псевдодифференциальные операторы в задачах дифракции) / А. С. Ильинский, Ю. Г. Смирнов. - М. : ИПРЖР, i996. - i76 с.

5. Смирнов, Ю. Г. Применение многочленов Чебышева к решению одномерных интегральных уравнений типа потенциала / Ю. Г. Смирнов. - Пенза : Изд-во Пенз. гос. тех. ун-та, i994.

6. Плещинский, Н. Б. Уравнение Гельмгольца в полуплоскости и скалярные задачи дифракции электромагнитных волн на плоских металлических экранах / Н. Б. Плещинский. - Препринт ПМФ-03-02. Казань : Казанск. матем. об-во, 2003. -30 с.

7. Плещинский, Н. Б. Модели и методы волноводной электродинамики : учеб. пособие / Н. Б. Плещинский. - Казань : Казан. гос. ун-т, 2008. - i04 с.

References

1. Khenl, KH. Teoriya difraktsii / KH. Khenl, A. Maue, K. Vestpfal'. - M. : Mir, i964. - 428 s.

2. Morse, P. M. The diffraction of waves by ribbons and by slits / P. M. Morse,

P. J. Rubinstein // Phys. Rev. - i938. - V. 54, № ii. - P. 895.

3. Nobl, B. Metod Vinera - Khopfa / B. Nobl. - M. : IL, i962. - 280 s.

4. Il’inskiy, A. S. Difraktsiya elektromagnitnykh voln na provodyashchikh tonkikh

ekranakh (Psevdodifferentsial'nyye operatory v zadachakh difraktsii) / A. S. Il'-inskiy, YU. G. Smirnov. - M. : IPRZHR, i996. - i76 s.

5. Smirnov, YU. G. Primeneniye mnogochlenov Chebysheva k resheniyu odnomernykh integral'nykh uravneniy tipa potentsiala / YU. G. Smirnov. - Penza : Izd-vo Penz. gos. tekh. un-ta, i994.

6. Pleshchinskiy, N. B. Uravneniye Gel'mgol'tsa v poluploskosti i skalyarnyye zadachi difraktsii elektromagnitnykh voln na ploskikh metallicheskikh ekranakh / N. B. Pleshchinskiy. - Preprint PMF-03-02. Kazan' : Kazansk. matem. ob-vo, 2003. -30 s.

7. Pleshchinskiy, N. B. Modeli i metody volnovodnoy elektrodinamiki : ucheb. posobiye / N. B. Pleshchinskiy. - Kazan' : Kazan. gos. un-t, 2008. - i04 s.

Александрова Ирина Леонидовна

ассистент, кафедра прикладной математики, Казанский федеральный университет

Aleksandrova Irina Leonidovna Assistant, sub-department of applied mathematics, Kazan Federal University

E-mail: iralexand@ksu.ru

Плещинский Николай Борисович

доктор физико-математических наук профессор, заведующий кафедрой прикладной математики, Казанский федеральный университет

E-mail: pnb@ksu.ru

УДК 517.958 Александрова, И. Л.

Проводящий тонкий экран в волноводной структуре: задача дифракции и задача трансмиссии / И. Л. Александрова, Н. Б. Плещинский // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2013. - № 1 (25). - С. 45-60.

Pleshchinskiy Nikolay Borisovich Doctor of physical and mathematical sciences, professor, head of sub-department of applied mathematics, Kazan Federal Unversity