Численное моделирование строгой математической дифракции на прямоугольном клине Текст научной статьи по специальности «Математика»

%

(N

тН

О

(N

УДК 535.42; 517.9 ББК 22.311

А.В. Боровский А.Л.Галкин

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ СТРОГОЙ МАТЕМАТИЧЕСКОЙ ДИФРАКЦИИ НА ПРЯМОУГОЛЬНОМ КЛИНЕ

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

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

A.V. Borovsky A.L. Galkin

NUMERICAL SIMULATION OF RIGOROUS MATHEMATICAL DIFFRACTION ON RECTANGULAR WEDGE

SSI

2sê Es =

Ép

НА

Ой"

°¡<s ^ 1

i

^ Л £

« FI 8sî ЕнО а

£м§|

МО |

аи-i

нее

К

н

н

о

н

PQ

СО

Н

The article deals with the problem of rigorous diffraction on a perfectly conducting rectangular wedge (step). The authors propose the exact solution to the Helmholtz equation for this problem in case when the incident wave has one reflection (on horizontal surface of the wedge), and calculate the intensity of the diffracting electromagnetic wave with TE-polarization. It is shown that a vertical wall increases the diffraction of the original wave in the angle of incidence close to the vertical and demonstrates the first diffraction peak of the line that separates light and shadow. Reducing the Sommerfeld angle of the incident wave is makes the solution for the step similar to the Sommerfeld’s solution for the half-plane; in this case the impact of the vertical wall on the effect of diffraction decreases.

Keywords: diffraction on rectangular wedge; numerical simulation of mathematical diffraction.

Математическая модель дифракции (рассеяния световой волны на физическом объекте) была сформулирована более 100 лет назад Г. Кирхгофом и А. Зоммерфельдом [13; 14]. Она сводилась к решению краевой задачи для уравнения Гельмгольца в области, ограниченной одной или несколькими поверхностями, и с условиями Зоммерфельда на бесконечности для рассеянной волны. В связи с математической сложностью данной задачи Кирхгоф предложил сделать ряд упрощений. Впоследствии эти предельные случаи стали называться дифракцией Фраунгофера (описывает дальнюю зону светорассеяния) и дифракцией Френеля (описывает ближнюю зону светорассеяния). Решение задачи дифракции в рамках строгой модели впервые удалось получить А. Зоммерфельду [9]. Им был

© А.В. Боровский, А.Л. Галкин, 2012

%

(N

тН

О

(N

SS1

2нё Яч >

MhSI

HS5

n3i

Ен^ i

ой-

°i<s ^ 1

5s s

^ Л £ ё§|

« ff I

8SI

"Jo

ЕнО а

£м§ц

МО f £И4

нее

«

н

н

о

н

PQ

со

Н

предложен метод контурных интегралов в комплексной плоскости для решения двухмерной краевой задачи дифракции. Впоследствии большой вклад в исследование этого метода внесли Ф. Франк и Р. Мизес [12]. Мы будем называть данный метод — методом Зоммерфельда-Франка-Мизеса (ЗФМ). В дальнейшем задачами строгой дифракции занимались М. Борн и Э. Вольф, Г. Гринберг и другие авторы [1; 8]. Постановка задачи строгой дифракции имеется в курсе теоретической физики Л. Ландау и Е. Лифшица, однако, ее решение не приводится, что говорит о незавершенности исследований [10].

В данной работе рассмотрим строгую математическую дифракцию на прямоугольном металлическом клине. Как ни странно, решение этой задачи отсутствует в трудах классиков теоретической физики, упомянутых выше. Хотя, именно эта задача представляет повышенный практический интерес, поскольку моделирует взаимодействие световых пучков с краями оправ и диафрагм в оптических устройствах. На важность этой проблемы указано в [6]. Моделирование строгой дифракции на прямоугольном клине (иногда такой клин называют ступенькой) представляет так же познавательный интерес, так как исследуется строгое решение уравнений Максвелла, описывающих взаимодействие электромагнитного поля с реальным физическим объектом.

В данной работе авторы опираются на выполненные ранее исследования дифракционного интеграла ЗФМ [2; 3] и численное моделирование строгой математической дифракции на полуплоскости для электромагнитных волн, обладающих TE- и TH-поляризациями [4; 5].

Постановка задачи. Рассмотрим клин с углом при вершине 5 (для прямоугольного клина 5 = п / 2), рис. 1, на который падает внешняя плоская электромагнитная волна, характеризуемая волновым вектором k. В полярной системе координат вектор k задается модулем k = |k| = 2п / X, где X — длина волны, и полярным углом а. Точка наблюдения задается модулем радиус-вектора r = |r| и полярным углом ф.

Выражение для падающей плоской волны в полярных координатах дается следующей формулой [1; 9; 12]

U+ (р, 0) = exp(ipcos0), р = kr, 0=ф-а.

(1)

Рис. 1. Дифракция плоской волны на проводящем клине в области I

Краевая задача дифракции ставится следующим образом: требуется найти решение двухмерного уравнения Гельмгольца

1 д ( dU) 1 d2U Т7 „

р— I +—т-^г + U = 0

р др! г др I р2 дф2

(2)

%

тН

о

<М

«К1

2нё Яч >

нйе

«31

оК

^ 1

I

^ Л £ ё§! ь Е

« Я41

8Н! й§ о

ЕнО а

£м§|

МО I £И4

нее

и

н

н

о

н

м

со

н

в области углов 0 < ф < 2п - 5 при любых значениях радиуса 0 < р < да. На поверхностях клина ставится условие равенства нулю решения

и\Ф= 0 = 0, и| Ф= 2л-8 = 0 (3)

В силу линейности уравнения Гельмгольца решение краевой задачи можно разложить на сумму решений геометрической оптики и дифракционной составляющей

и = иё + Ud. (4)

На дифракционную часть решения и(1 ставится условие Зоммерфель-да — эта часть решения должна стремиться к исчезающей на бесконечности цилиндрической волне

(5)

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

и(р, 0) = ехр^(я: + рcos0)), р = ^, 0=ф-(2п-а),

что дает

Ц(р, 0) = -exp(iрcos(ф + а)).

Таким образом, сумма падающей и отраженной волн геометрической оптики записывается в виде:

и = и+ (р, ф-а) + и-(р, ф+а) = exp(iрcos(ф-a))-exp(iрcos(ф + a)).

В силу четности косинуса, решение геометрической оптики обращается в нуль на верхней грани клина

и |ф= 0 = и+ (р, -а) + и-(р, +а) = exp(iрcos(-а))-exp(iр^(+а)).

Дифракционную составляющую и11 целесообразно искать в виде разности двух функций, описывающих искажение падающей и отраженной волн геометрической оптики

Ud (ф, а) = V(ф - а) - ^ф + а).

Построим решение краевой задачи, используя функции ^0) с периодом Т = 2(2п - 5). Такие периодические функции должны быть четными. Если решение будет найдено, то, согласно теореме единственности для уравнения Гельмгольца [7; 11], оно будет единственным решением краевой задачи.

На верхней грани клина решение обращается в нуль

Ц (0, а) = ^-а) - ^+а) = 0, в силу четности функции ^0). На нижней грани клина

иЛ (2п - 5, а) = ^(2я - 5) - а) - ^(2я - 5) + а) = 0, в силу зеркальной симметрии функции V относительно половинки периода. Перечисленными свойствами обладает дифракционный контурный интеграл ЗФМ и^ (0) с периодом 2пр, где р = Т / 2п. (Теорию контурного интеграла ЗФМ см. в [2; 3]).

Общие вопросы. Аналитическое выражение для дифракционного контурного интеграла ЗФМ, обладающего периодом 2пп и представляющего сумму решения геометрической оптики и дифракционной добавки, рассмотрено в [3]. Это выражение имеет вид:

%

(N

тН

О

(N

SS1

2sé Яа >

MhSI

HSS

n3i

Ен^ i

ой-

°¡<s ^ 1

5s i

^ Л £ ё§|

« Р I

8sí й§ о

ЕнО а

£м§ц

МО f £И4

нее

«

н

н

о

н

PQ

со

Н

i

X —

n

UC (р, 0) = exp(¿p cos0) X 1+n]T1cos U 1 - mY| exp {-i nml sin {e-^os0 Н2П (^

n

2n

U (0) = -;

n

n

(6)

lim Uc (р, 0) =

р ^ » n ’ '

|exp(ipcos0), 0| < n i 0, 01 > n

Здесь Hf/П — вторая функция Ханкеля с дробным индексом. В теории цилиндрических функций определены первая и вторая функции Ханкеля. Они отличаются комплексным сопряжением. Действительная часть второй функции Ханкеля равна функции Бесселя, а мнимая часть равна функции Неймана с отрицательным знаком. Выражение в квадратных скобках удовлетворяет «центральному тождеству теории дифракции для клина» (ЦТД), доказательство которого приведено, например, в [3].

1+Icos К1 - n))exp {-i In)sin (n ]!exp(-i4 cos0 )Hmn®^s

In, 0| <n ¡0, 01 >n

(7)

При рассмотрении задачи дифракции на прямоугольном клине исходным является дифракционный интеграл с периодом 6п. Для него целочисленный индекс п принимает значение п = 3. При этом, целые числа т равны т = 1, 2. Выражение для дифракционного интеграла ЗФМ в этом случае принимает следующий вид:

U3 (р, 0) = exp(ip cos0) X

1

X — 2

1 +

л/э Г Vs

2

Тэ г 2

V

2

V

1

— i 2

,Vs

2

Г 2

1) — cosí —

2 i ( 3

Г 1_)

cos —0

(3 J

(8)

Решение уравнения Гельмгольца для случая клина с прямым углом строится с применением теоремы о понижении периода. Период интеграла (8) понижается с 6п до 3п. Окончательное строгое решение для задачи дифракции на прямоугольном клине имеет вид:

/3е(р, ф-у) -ис^

1 (9)

Здесь (у = а - п) — угол Зоммерфельда. Функция (9) обращается в нуль на горизонтальной и вертикальной гранях ступеньки, т.е. при значениях угла наблюдения ф = 0, 3п / 2.

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

Е(г, ^ = ReU(r)cos(юt) + 1ти(г^т(ш^, (10)

и(р,ф,Y) = [UC(р, ф-у)-U(р, ф+у)]1 +

+ [U(р, 3п + (ф-у))-UC(р, 3п + (ф+у))]2.

m = 1

р

р

%

тН

о

<М

«К1

2нё Яч >

нйе

«31

оК

^ 1

I

^ Л £ ё§! ь Е

« Я" I

8Н! й§ о

ул* ЕнО а

£м§|

МО I £И4

нее

и

н

н

о

н

м

со

н

усреднив по быстрым осцилляциям вектор Пойнтинга, получим следующее выражение для интенсивности электромагнитной волны [3; 4; 5]

I =

1

„ ТТд 1ти Т ттдЯви Ив и-------1ти

Эф

Эф

+

„ Г7Э 1ти Т Г7ЭЯви Яви--------1ти

др

др

(11)

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

Ниже представим результаты численного моделирования задачи строгой дифракции для электромагнитной волны ТЕ поляризации на прямоугольном клине. Математическая модель для нее описывается выражениями (8)-(11).

Методика вычислений. Формулы модели представляют значительную сложность для вычислений. Цилиндрическая функция Ханкеля содержит в мнимой части функцию Неймана, которая имеет особенность в нуле. Поэтому формулы модели содержат несобственные в нуле интегралы от осциллирующих функций, одна из которых является цилиндрической. Раскрытие интегрируемой особенности мнимой части функции Неймана представляло главную сложность проводимых вычислений. Кроме того, выражение для интенсивности содержит производные от решения задачи (9), вычисление которых представляет отдельную проблему численных методов.

Для проведения численных расчетов применялись пакеты объектного программирования «Маткад» и «Математика». Первый пакет позволяет проводить вычисления несобственных интегралов непосредственно, при этом, время расчетов не является оптимальным. Для применения второго пакета были проделаны следующие оценки. В интервале [0, 1] справедливы асимптотики с константами Ац, Вц

1 2

1тНЦ»© - В„), М = 3.3,

мнимая экспонента в (8) изменяется слабо, поэтому относительная погрешность в связи с обрезанием особенности интеграла на уровне е равна

5 = оМ 1

11т НЦ2)(^

£ (1 + £2МАм /Вм(1 -м)/(1 + М))

1 + Ам /Вц(1 -м)/(1 + м) .

При е = 10-6 51/з = 1,3 ■ 10-4; 52/3 = 8,9 ■ 10-3, т.е. погрешность вычисления интегралов становится меньше 1%, но наибольшая величина 1тНЦ2)(е) достигает значения 1тН^е) = 6,84103. Для построения интерполяционных функций в пакете «Математика» на интервале [е, 1] далее применяется линейная интерполяция на логарифмической сетке {т10-й}, т = 1, 2, ..., 9; й = 0, 1, ..., 7. Увеличение порядка интерполяции из-за отмеченного сильного перепада функций Ханкеля при подходе к особенности приводит к появлению ложной пичковой структуры в численном решении. Для определения производных, входящих в выражение (11) для интенсивности, используется дифференцирование по ф

2

2

2

с4

тН

О

(М

«К1

2нё Ма >

М^^!3|

и«!

ёК

н- 1

*о!

ь Е

«¡г I

§н! й§ о

и^«

^м!ц

М0|

РнИ-1

К®«

и

н

н

о

н

м

со

н

численного решения при фиксированном р, производные по р определяются как центральные разности по двум разнесенным слоям по р. Разработанная методика для расчетов с использованием пакета «Математика» позволила заметно сократить время расчетов.

Тестирование вычислений. Тестирование вычислительных программ проводилось следующим образом. Во-первых, была проверена идентичность функции Ханкеля в обоих пакетах. Во-вторых, была осуществлена проверка вычисления несобственного интеграла, которая показала, что первоначальные проблемы в пакете Математика, связанные с потерей точности при учете особенности для функции Неймана в нуле, преодолены. В-третьих, была протестирована на аналитическом примере эффективность счета интегралов от двух осциллирующих функций с затуханием. Данная проверка дала положительный результат. В-четвертых, было проверено асимптотическое выполнение в численных расчетах центрального тождества теории дифракции (7). Результаты такой проверки для функции Ф(р, 0) = и3С(р, 0)вхр(-фcos0) представлены на рис. 2 и 3. Вычисления асимптотически совпадают с аналитическими результатами. В-пятых, было проверено совпадение решения (9) при р ^ да, ф ^ 0 с решением задачи строгой дифракции на полуплоскости [2; 5]. Параллельно была выяснена роль скобок [□]1 и [П]2 в (9). Результаты этого исследования представлены на рис. 4 и 5. Выводы следующие: основной вклад в решение вносит первая скобка; при стремлении ф ^ 0 вклад второй скобки несуществен; при стремлении ф ^ 3п / 2 (вблизи вертикальной стенки) вклад второй скобки становится существенным; вторая скобка обеспечивает обращение в нуль решения уравнения Гельмгольца на вертикальной стенке. Исследования тестового характера подтвердили правильность математической модели и численных расчетов.

Рис. 2. Асимптотика функции Ф(р, 9) при р ^ да для углов |9| > п, 9 = 3п + п / 3, (а) действительная часть, (б) мнимая часть функции Ф(р, 9)

100

б

Рис. 3. Асимптотика функции Ф(р, 9) при р ^ да для углов |9| < п, 9 = п / 3, (а) действительная часть, (б) мнимая часть функции Ф(р, 9)

б

а

а

2

с4

тН

О

(М

«К1 Ма г

н!®1

яЗ|

Ен< |

Ъ<1

и- 1 ё§| «¡л

Он®

15881

ьм!| ио г ли» К®«

и

н

н

о

н

м

со

н

Численное моделирование. Ниже представим результаты численного моделирования строгой дифракции на ступеньке (на прямоугольном клине) для случая, когда угол Зоммерфельда для падающей волны ненамного отличается от прямого угла: у = п / 2 - 0,05. При этом, следует ожидать значительного влияния вертикальной стенки. Расчеты действительной и мнимой частей решения уравнения Гельмгольца для этого случая уже были представлены на рис. 4 (в) и 5 (в). Из них видны три характерные области углов, в которых динамика волн выглядит по-разному. В диапазоне углов 0 < ф < у ~ п / 2 имеет место интерференция падающей и отраженной от верхней грани ступеньки волн. Размах амплитуды интерференции превышает единицу. В область углов у < ф < п + у попадает только падающая волна, отраженной волны здесь не возникает. Интерференции в данной области углов нет. Размах амплитуды волны равен единице. В диапазоне углов у + п < ф < 3п / 2 имеет место геометрическая тень. Эффект дифракции размывает амплитуду поля в районе перехода к геометрической тени. Наличие вертикальной стенки отжимает электромагнитное поле, и максимум амплитуды смещается влево от положения строгой границы свет-тень. Поведение амплитуды отличается плавным переходом от максимального значения к нулю.

Рис. 4. Вклад скобок 1 и 2 и их сумма. Дифракция на ступеньке. Действительная часть решения уравнения Гельмгольца для случая р = 50, у = -0,05 + п / 2. (600 точек): (а) скобка — 1, (б) скобка — 2, (в) — сумма

На рис. 6 представлено угловое распределение интенсивности для рассматриваемого случая. Падающая и отраженная волны при углах Зоммер-фельда близких к прямому гасят друг друга в результате интерференции. В области углов у < ф < 3п / 2 прошедшая волна сильно возмущена эффектом дифракции. Дифракция усиливается при приближении к вертикальной стенке. Если в плоской невозмущенной волне интенсивность равна единице, то в результате дифракции интенсивность в основном максимуме возрастает в 1,7 раза. Наличие стенки довольно сильно отжимает главный пик дифракции влево. Дифракция на ступеньке развивается сильнее, чем на полуплоскости, на что указывает нижний рисунок. Интенсивность главного пика для ступеньки 1,7, а для полуплоскости 1,4.

в

2

с4

тН

О

(М

часть решения уравнения Гельмгольца для случая р = 50, у = -0,05 + п / 2. (600 точек): (а) скобка — 1, (б) скобка — 2, (в) — сумма

«К1

Ма > МнВ|

яЗ|

вК

Н- 1

ь Е

«¡г I

§н! й§ о

и^«

^м||

М0|

РнИ-1

К®«

И

н

н

о

н

м

со

н

Рис. 6. Дифракция на ступеньке (а) и полуплоскости (б).

Угловое распределение интенсивности для случая: р = 50, у = п / 2 - 0,05

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

Рис. 7. Дифракция на ступеньке. Интенсивность для случая р = 50, у = п / 2 - 0,4 (а), у = п / 2 - 0,8 (б)

%

(N

тН

О

(N

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

SS1

Д S =

КнИ

нД

НА

Ой"

^ 1

5s ü

«Л*

h Е

8 ff I

8SI

"So ЕнО а

£м§|

МО I

нее

К

н

н

о

н

PQ

со

Н

Список использованной литературы

1. Борн M. Основы оптики / M. Борн, Э. Вольф. — М.: Наука, 1973. — 720 с.

2. Боровский А.В. Дифракционный интеграл Франка-Зоммерфельда / А.В. Боровский, А.Н. Бородин. — Иркутск: Изд-во БГУЭП, 2012. — 64 с.

3. Боровский А.В. Математическая теория дифракции плоской волны на идеально-проводящем клине / А.В. Боровский, А.Н. Бородин, А.Л. Галкин. — Иркутск: Изд-во БГУЭП, 2013. — 104 с.

4. Боровский А.В. Численное моделирование дифракции TH волны на идеально проводящей полуплоскости / А.В. Боровский, А.Н. Бородин, А.Л. Галкин // Известия Иркутской государственной экономической академии. — 2012. — № 5. — С. 124-129.

5. Боровский А.В. Численное моделирование эффектов дифракции / А.В. Боровский, А.Н. Бородин, А.Л. Галкин // Известия Иркутской государственной экономической академии (электронный журнал). — 2012. — № 4. — URL: http://eizvestia.isea.ru/reader/article.aspx?id=13859.

6. Бородин А.Н. Измерение светорассеяния в зависимости от радиуса кривизны края апертурной диафрагмы / А.Н. Бородин, А.Н. Малов, С.А. Чупра-ков // Компьютерная оптика. — 2005. — Вып. 28. — С. 87-88.

7. Владимиров В.С. Уравнения математической физики / В.С. Владимиров. — М.: Наука, 1982. — 512 с.

8. Гринберг Г. Избранные вопросы математической теории электрических и магнитных явлений / Г. Гринберг. — М.; Л.: Изд-во АН СССР, 1948. —727 с.

9. Зоммерфельд А. Оптика / А. Зоммерфельд. — М.: Иностр. лит-ра, 1953. — 486 с.

10. Ландау Л.Д. Теоретическая физика: в 10 т. / Л.Д. Ландау, Е.М. Лиф-шиц. — Т. 8: Электродинамика сплошных сред. — М.: Наука, 1982.— 621 с.

11. Тихонов А.Н. Уравнения математической физики / А.Н. Тихонов, А.А. Самарский. — М.: Наука, 1972. — 735 с.

12. Франк Ф. Дифференциальные и интегральные уравнения математической физики / Ф. Франк, Р. Мизес. — М.; Л.: ОНТИ, 1937. — 998 с.

13. Kirchhoff G. // Ann.d.physik. — 1883. — № 18. — P. 663.

14. Sommerfeld A. Mathematische Theorie der Diffraction / A. Sommerfeld // Mathematische Annalen. — 1896. — Vol. 47. — Pp. 317-374.

References

1. Born M. Osnovy optiki / M. Born, E. Vol’f. — M.: Nauka, 1973. — 720 s.

2. Borovskii A.V. Difraktsionnyi integral Franka-Zommerfel’da / A.V. Boro-vskii, A.N. Borodin. — Irkutsk: Izd-vo BGUEP, 2012. — 64 s.

%

(N

тН

О

(N

3. Borovskii A.V. Matematicheskaya teoriya difraktsii ploskoi volny na ideal’no-provodyashchem kline / A.V. Borovskii, A.N. Borodin, A.L. Galkin. — Irkutsk: Izd-vo BGUEP, 2013. — 104 s.

4. Borovskii A.V. Chislennoe modelirovanie difraktsii TH volny na ideal’no provodyashchei poluploskosti / A.V. Borovskii, A.N. Borodin, A.L. Galkin // Izvestiya Irkutskoi gosudarstvennoi ekonomicheskoi akademii. — 2012. — № 5. — S.124-129.

5. Borovskii A.V. Chislennoe modelirovanie effektov difraktsii / A.V. Borovskii, A.N. Borodin, A.L. Galkin // Izvestiya Irkutskoi gosudarstvennoi ekonomicheskoi akademii (elektronnyi zhurnal). — 2012. — № 4. — URL: http://eizvestia. isea.ru/reader/article.aspx?id=13859.

6. Borodin A.N. Izmerenie svetorasseyaniya v zavisimosti ot radiusa kriviz-ny kraya aperturnoi diafragmy / A.N. Borodin, A.N. Malov, S.A. Chuprakov // Komp’yuternaya optika. — 2005. — Vyp. 28. — S. 87-88.

7. Vladimirov V.S. Uravneniya matematicheskoi fiziki / V.S. Vladimirov. — M.: Nauka, 1982. — 512 s.

8. Grinberg G. Izbrannye voprosy matematicheskoi teorii elektricheskikh i mag-nitnykh yavlenii / G. Grinberg. — M.; L.: Izd-vo AN SSSR, 1948. -727 s.

9. Zommerfel’d A. Optika / A. Zommerfel’d. — M.: Inostr. lit-ra, 1953. — 486 s.

10. Landau L.D. Teoreticheskaya fizika: v 10 t. / L.D. Landau, E.M. Lifshits. — T. 8: Elektrodinamika sploshnykh sred. — M.: Nauka, 1982. — 621 s.

11. Tikhonov A.N. Uravneniya matematicheskoi fiziki / A.N. Tikhonov, A.A. Samarskii. — M.: Nauka, 1972. — 735 s.

12. Frank F. Differentsial’nye i integral’nye uravneniya matematicheskoi fiziki / F. Frank, R. Mizes. — M.; L.: ONTI, 1937. — 998 s.

13. Kirchhoff G. // Ann.d.physik. — 1883. — № 18. — P. 663.

14. Sommerfeld A. Mathematische Theorie der Diffraction / A. Sommerfeld // Mathematische Annalen. — 1896. — Vol. 47. — Pp. 317-374.

SSI

2sl

KS г

MhSI

«31

I

^ 1

Ps ü

S I

8s! ИН 25 о

СнО!

£м§||

МО f

нее

н

н

о

н

PQ

со

н

Информация об авторах

Боровский Андрей Викторович — доктор физико-математических наук, профессор, кафедра информатики и кибернетики, Байкальский государственный университет экономики и права, г. Иркутск, e-mail: [email protected].

Галкин Андрей Леонидович — доктор физико-математических наук, заведующий сектором, Институт общей физики РАН им. А.М. Прохорова, г. Москва, e-mail: [email protected].

Authors

Borovsky Andrey Victorovich — Doctor of Physical and Mathematical Sciences, Professor, Chair of Computer Science and Cybernetics, Baikal State University of Economics and Law, Irkutsk, e-mail: [email protected].

Galkin Andrey Leonidovich — Doctor of Physical and Mathematical Sciences, Head of Sector, General Physics Institute of RAS named after A.M. Prokhorov, Moscow, e-mail: [email protected].