Сравнительная характеристика различных систем базисных функций полной области применительно к решению интегральных уравнений Фредгольма второго рода Текст научной статьи по специальности «Математика»

УДК 519.62

СРАВНИТЕЛЬНАЯ ХАРАКТЕРИСТИКА РАЗЛИЧНЫХ СИСТЕМ БАЗИСНЫХ ФУНКЦИЙ ПОЛНОЙ ОБЛАСТИ ПРИМЕНИТЕЛЬНО К РЕШЕНИЮ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ ФРЕДГОЛЬМА ВТОРОГО РОДА

© 2012 Е.А. Алашеева,1 М.Ю. Маслов2

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

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

1. Общие замечания

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

1 Алашеева Елена Александровна (allena-81@mail.ru), кафедра высшей математики Поволжского государственного университета телекоммуникаций и информатики, 443010, Российская Федерация, г. Самара, ул. Льва Толстого, 23.

2Маслов Михаил Юрьевич (mike@psati.ru), кафедра электродинамики и антенн Поволжского государственного университета телекоммуникаций и информатики, 443010, Российская Федерация, г. Самара, ул. Льва Толстого, 23.

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

2. Исходные уравнения

Пусть элементарный излучатель находится в точке с координатами г = = {х',у',х'} (см. рисунок) вблизи проводящего тела. Амплитуда тока, возбуждающего излучатель I, длина излучателя Ь Здесь и далее штрихованные координаты обозначают точки источников первичного поля. Геометрические параметры задачи приведены на рисунке.

Рис. К выводу интегральных уравнений

Здесь г = {х,у,х} — радиус-вектор точки наблюдения, г = {х',у',х'} — радиус-вектор точки источника, г = {х'',у'',х''} — радиус-вектор точки на поверхности экрана.

Пусть интегральное уравнение имеет смысл граничного условия для магнитного вектора [1] на поверхности проводящего тела

НТо - ИТз = 0, (2.1)

где НТо и ИТз — тангенциальные компоненты вектора напряженности первичного магнитного поля и поля рассеяния на поверхности

Для тангенциальной компоненты магнитного поля рассеяния справедливо выражение [2]:

НтЗ(г) = [ 38(?') х [дгааа(г,Г')]¿5\ (2.2)

Js

где 38 (г'' ) — вектор плотности тока на поверхности п — вектор нормали в точке с координатами г", 100(г,г") = ехр(—1к\т'— г''\)/(4п\г — г"\) — функция Грина для свободного пространства, к = 2п/\ — волновое число. Штрих у оператора дтав,

в (2.2) указывает на то, что дифференцирование производится по координатам Г'. Выражение (2.2) с учетом граничного условия (2.1) можно представить в виде:

n(r0) х H0(r0) = n(r0) x lim I jS(Г') х [gradG(r,r")]'dS, (2.3)

r^roJ S

где ro — текущая координата на поверхности S.

Вычисление предела в (2.3) с помощью теории потенциала приводит к следующему интегральному уравнению [3]:

П(го) х Йо(то) = 1j(T) - j П(?о) х {js(?') х [gradG(r, r'')]'}dS. (2.4)

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

ri(ro) х T2(ro) = П(го). (2.5)

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

ri(fo) х n(ro) = -T2(ro),T2(ro) х n(ro) = ri(ro). (2.6)

Тогда (2.4) распадается на два уравнения относительно двух ортогональных компонент вектора плотности тока:

T2(ro) ■ H(ro) = -1 Ti(ro)js(r) - I T2(ro) ■ {js(r') х [gradG(r,r')]'}dS. (2.7)

ti( ro) ■ H(ro) = -2?2(?o)js(r) + j Ti(ro) ■ {js(?') х [gradG(r, r")]'}dS. (2.8)

Далее в зависимости от вида элементарного излучателя левые части уравнений будут иметь различный вид [2].

3. Применение различных базисов к аппроксимации токовых функций на проводящих поверхностях

Систему базисных функций следует выбрать так, чтобы элементы матрицы импедансов [1] вычислялись по возможности быстро, а сама матрица была удобна для дальнейшего ее обращения. Удачный выбор системы базисных функций позволяет добиться более быстрой сходимости решения, что приводит к меньшему объему расчетов при заданной точности.

Как правило, чем точнее базисные функции ¡п передают характер действительного распределения тока на излучателе, тем лучше сходимость и в некоторых случаях тем более обусловлена матрица импедансов (т. е. более устойчиво решение). Традиционно базисы подразделяются на два типа [4]. К первому типу относятся базисные функции, определенные и не равные нулю во всей области определения оператора Ь, за исключением, может быть, множества точек меры нуль. Систему базисных функций этого типа мы называем базисом полной области. В отличие от них базисные функции второго типа установлены во всей области определения оператора Ь, но обращаются в нуль в части этой области. Такая система базисных функций называется базисом подобластей.

Системы базисных функций полной области

Как отмечалось раньше, это такие базисные функции /п,которые определены и отличны от нуля во всей области существования A, причем A = n In/n. Теоретически таких систем функций можно выбрать бесконечно много. Однако выбор системы базисных функций определяется многими конкретными обстоятельствами, связанными с рассматриваемой задачей. Приведем следующие примеры базисных функций (табл. 1).

Таблица 1

Примеры различных базисов полной области

Разложение функции A в ряд Вид базисных функций

Разложение Фурье: A(y) « h cos ПУ + h cos П + cos П + ... fn = COS (2n-lhy

Разложение Чебышева: A(y) « IT (y) + I2T2(y) + I3T4(y) + ... fi = To(y) = l, f2 = T2(y) = 2y2 - 1

Разложение Маклорена: A(y) « Ii + I2y2 + 1зу4 + ... fn = y2n-2

Для двумерного случая, когда решение А(х) является функцией двух переменных х _ (х1,х2), в качестве базисных функций удобно выбирать двумерные аналоги приведенных базисов, например, двумерный ряд Фурье.

Базис полной области чаще всего используется, когда искомое решение А(х) предполагается достаточно гладким и область определения интегрального оператора правильной формы.

Системы базисных функций подобластей

В одномерном случае, когда решение А(у) зависит от одной переменной, область определения оператора Ь разбивается на сегменты Ау^, и задается набор функций, равных нулю во всех сегментах, кроме какого-то одного. Приведем примеры таких базисов.

При решении двумерной задачи можно использовать аналогичный подход. Теперь функция А(х) зависит от двух переменных х _ (х1,х2), а область определения оператора Ь представляет собой некоторую плоскую фигуру Q. Функции /п(х) могут быть определены в различных подобластях области определения оператора Ь, так что каждая функция /п(х) существует лишь в некоторой подобласти Qn. В случае интегрального уравнения, когда Ь является интегральным оператором, применение таких базисных функций эквивалентно замене интеграла на конечную риманову сумму. Это можно записать следующим образом [4]:

А( х) _ < ап/п(х),х ^ Qn ,

1 0, в противном случае

4. Численное решение

Рассмотрим ИУ Фредгольма второго рода, в общем виде записываемое как: /• Ь /• ь

/ / К(х,у,п,у)Л(п,у)в,пвм + J(х,у) _ /(х,у) (4.9)

.) а ./а

Таблица 2

Примеры различных базисов подобластей

Наименование функций

Вид базисных функций

1. Базисные сплайны

сплайн нулевой степени:

- линейный сплайн:

0,У < Уг, к ,Уг < У < Уг+1, 0,У > Уг+1 0,У < Уг,

1 + У-Уi+l к +

к2 У-Уь+1 к2

Уi < У < Уг+1, Уг+1 < У < Уг+2,

-кубический сплайн:

/г = <

0,У > Уг+2 0,У < Уг,

Ш(У - Уг)3,Уг < У < Уг+1, 6к + 2к*(У -Уг) + иЬ(У - Уг)2--2к4 (У - Уг)3,Уг+1 < У < Уг+2, 6к + (Уг+з - У) + (Уг+з - У? -- 2к4 (Уг+з - У)3, Уг+2 < У < Уг+3, вк4(Уг+4 - У)3, Уг+3 < У < Уг+4, 0,У > Уг+4

2.Кусочно-синусоидальные функции:

/г =

1БШк{У*+1-У!+&""к(У-Уг), У е [Уг, Уг+к], 0,У е [- ъ,уг}() У е [Уг+1, + ъ]

3. Базис из сплайновых вейвлет

-обыкновенная вейвлета:

Фг = <

0,У е [- ъ,у2г]

{У-к12г) ,У е [У2г, У2г+1], (У2^+2-7к+бУ2^+1) ,У е [У2г+1 Уг+2], -е(У2^ + к6У-10У2^ ,У е [у2г+2,У2г+з], 6(У2^-16У+10У2^ ,У е [У2г+з,У2г+4], -(У2^ + к-бУ2^ ,У е [У2г+4,У2г+5], ,У е [У2г+5, У2г+б],

0,У е [У2г+1б, + ъ]

4. Тригонометрическая интерполяция:

Л + Б Бт к(у - Уг)+ ¡г = { +С сов к (у - Уг),У е [уг,Уг+1], 0,У е [- ъ,уг]Ц У е [Уг+1, + ъ]

с заданной функцией /(х,у) и неизвестной функцией . (х,у). При решении подобных уравнений методом Галеркина, в котором весовые и базисные функции совпадают, имеется известная проблема выбора базиса, адекватного поставленной задаче. Нужно выбрать такую систему базисных функций, чтобы матрица импе-дансов [1] была как можно меньшего размера, время ее заполнения тоже минимальным и полученное решение соответствовало физическим представлениям о нем (табл. 2).

Можно, конечно, выбрать в качестве базиса кусочно-линейный базис, или, как его еще называют, базис из сплайнов первой степени. При таком выборе мы могли бы решить задачу даже в случае довольно сложной области определения у

интегрального оператора (просто нужно было бы выбрать более густую сетку разбиения для построения базисных функций). Однако решение задачи осложняется тем, что при дискретизации интегрального уравнения получается система линейных алгебраических уравнений с плотной матрицей. А такая система решается прямыми методами, например Гаусса, довольно долго (даже учитывая возможности современных ЭВМ, порядка нескольких часов) [4].

Также можно выбрать в качестве базисных функций вейвлет-базис. В таком случае матрица импедансов будет псевдоразреженной [5; 6], что позволит экономить машинное время при решении системы, а также память ЭВМ при заполнении самой матрицы, используя разреженные технологии. Но в данной ситуации сама матрица системы линейных уравнений заполняется очень долго, т. к. довольно трудно предсказать положение "почти нулей" в матрице и, как следствие, трудно создать эффективно работающий алгоритм. Справедливости ради стоит заметить, что при небольшом усложнении формы поверхности, на которую наведено излучающее устройство, применение данного базиса становится наиболее актуальным [6]. Однако ограничимся рамками данного исследования.

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

2pi . 2pi 2pi 2pi 2pi J (x, y) = ai + cos + a3 sln + a4 cos ^^y + a5 cos cos ~Л~У + '' ''

где Л длина волны.

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

Таблица 3

Время нахождения численного решения модельной задачи с использованием различных базисов

Базисные функции Время нахождения численного решения модельной задачи

сплайны 20 ч

вейвлеты 5 ч

двумерный ряд Фурье 1 мин

Вообще исследования, проведенные авторами [7; 8], показывают, что тригонометрические базисы полной области имеют существенные преимущества в смысле ресурсоемкости реализованных алгоритмов, по сравнению с иными системами базисных функций, рассмотренными выше, при решении существенно двумерных электродинамических задач с гармоническим возбуждением, геометрия которых относительно проста (плоские рефлекторы, зеркала, образованные поверхностями второго порядка). Данное обстоятельство вполне соответствуют априорным физическим оценкам характера ожидаемого решения. Однако в случае сложной геометрии (зеркала, имеющие изломы) негармоническое (в особенности, импульсное) возбуждение, применение вейвлет-базисов, обуславливающих разреженные матрицы, представляются весьма перспективными (в случае импульсных воздействий

все входящие в уравнения (6) величины будут иметь смысл соответствующих Фурье-образов).

Выводы

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

Литература

[1] Вычислительные методы в электродинамике / под ред. Р. Митры; пер. с англ. Э.Л. Бурштейна. М.: Мир, 1977. 487 с.

[2] Никольский В.В. Антенны. М.: Связь, 1966.

[3] Владимиров В.С. Уравнения математической физики. М.: Наука, 1981. 250 с.

[4] Марчук Г.И., Агошков В.И. Введение в проекционно-сеточные методы. М.: Наука, 1981.

[5] Блатов И.А. Об алгебрах операторов с псевдоразреженными матрицами и их приложениях // Сибирский математический журнал. 1996. Т. 37. № 1. C. 36-59.

[6] Блатов И.А. О методах неполной факторизации для систем с разреженными матрицами // Журнал вычислительной математики и математической физики. 1993. Т. 33. № 7. C. 819-836.

[7] Алашеева Е.А., Блатов И.А., Маслов М.Ю. Решение задачи рассеяния электромагнитного поля элемента электрического тока проводящим экраном конечных размеров // Инфокоммуникационные технологии. 2007. № 2. C. 8-14.

[8] Алашеева Е.А., Блатов И.А., Маслов М.Ю. Решение задачи об излучении линейной структуры, расположенной вблизи идеально проводящего экрана, сводимой к двумерной системе уравнений Фредгольма второго рода // Вестник СамГУ. Естественнонаучная серия. 2010. № 2(76). C. 13-22.

Поступила в редакцию 22/IX/2012; в окончательном варианте — 22/IX/2012.

COMPARATIVE CHARACTERISTICS OF VARIOUS SYSTEMS OF BASIS FUNCTIONS OF THE FULL FIELD, WITH REGARD TO THE SOLUTION OF INTEGRAL FREDHOLM EQUATIONS OF THE

SECOND KIND

© 2012 E.A. Alasheeva? M.Y. Maslov4

The article examines some aspects of numerical solution of some problems of antenna electrodynamics, reducible to the integral equations of the second kind with the use of various bases of decomposition of the sought function.

Key words: integral equations, numerical methods, projective methods, basic functions.

Paper received 22/IX/2012. Paper accepted 22/IX/2012.

3Alasheeva Elena Alexandrovna (alasheeva@mail.ru), the Dept. of Higher Mathematics, Povolsh-skiy State University of Telecommunications and Informatics, Samara, 443010, Russian Federation.

4Maslov Mihail Yurievich (mike@psati.ru), the Dept. of Electrodynamics and Antenna, Povolsh-skiy State University of Telecommunications and Informatics, Samara, 443010, Russian Federation.