- •Введение
- •Основные уравнения
- •Принцип перестановочной двойственности
- •Энергетические характеристики.
- •Векторные и скалярные потенциалы. Волновые уравнения.
- •Дельта –функция
- •Функция Грина
- •Некоторые сведения из функционального анализа. Функционалы и операторы в гильбертовом пространстве
- •Проекционные методы Ортогональные ряды и проекционная интерпретация.
- •Вариационные принципы. Процесс Ритца.
- •Применение метода Ритца к анализу свойств резонаторов
- •Наиболее употребительные проекционные схемы в электродинамике.
- •Метод частичных областей. (процесс Трефтца)
- •Интегральные уравнения электродинамики.
- •Вывод интегральных уравнений
- •Интегральное уравнение Фредгольма.
- •Решение интегрального уравнения
- •Решение интегрального уравнения для диафрагмы в плоском волноводе.
- •Проблема устойчивости решения.
- •Мчо с учётом особенности на ребре.
- •Метод Моментов
- •Итерационные методы
- •Дискретизационные методы Метод коллокации.
- •Сшивание в дискретных точках.
- •Разностные схемы. Сеточные методы
- •Литература
Дельта –функция
Дельта-функция определяется следующим образом
причем так, что интеграл
.
Очевидно, что пределы могут быть любые, охватывающие начало координат
Из определения -функции следует, что если – любая непрерывная функция, то
, (22)
действительно
где – сколь угодно малая величина – .
-функция, четная функция, т.е. .
Действительно в соответствии с (22)
,
заменим в этой формуле на
Для -функции
действительно
.
-функция может быть получена как предел непрерывной функции
,
так как . а .
Трехмерная -функция равна нулю везде кроме начала системы координат, и интеграл по всему пространству равен 1.
Функция Грина
Функцией Грина называется некоторая вспомогательная функция, используемая при решении дифференциальных уравнений. Пусть у нас имеется уравнение Гельмгольца
(23)
Функция Грина удовлетворяет уравнению с той же правой частью и -функцией в качестве возбуждающего тока:
(24)
Это уравнение не определяет ещё полностью функцию G, нужно поставить ещё условия на границе области или условие излучения. Это можно сделать различным образом, то есть вводить для одного и того же уравнения различные функции Грина, мы вернёмся к этому ниже.
Функция Грина зависит от координат точки , по которой производится дифференцирование, и от координат точки , в которой располагается элементарный диполь, точнее –источник. В (24) играет роль параметра. При любых граничных условиях функция Грина симметрична:
что отвечает теореме взаимности: поле в точке , возбуждаемое источником, располагаемым в точке (левая часть), равно полю в точке , возбуждаемому источником, расположенным в точке (правая часть).
Использование функции Грина для решения основано на преобразовании объёмного интеграла в поверхностный по второй формуле Грина. Умножая (24) на и ( 23 ) на и вычтем из первого выражения второе
-
--------------------------------------------------------
Проинтегрируем по произвольной области V, содержащей точку
где см. ( 22 ).
В соответствии со второй формулой (теоремой) Грина
Область V ограничена поверхностью S. Внутри S U и U непрерывны и содержится точка . В конкретных задачах S должна содержать поверхность металлических тел и внешнюю поверхность, которую мы в некоторых случаях будем устремлять в .
(25)
Эта формула показывает, что для того, чтобы получить значения функции , то есть решение уравнения (23) в какой-либо точке , надо взять функцию Грина, порождаемую источником, расположенным именно в точке . Формула (25) даёт непосредственно решение (23), если G подчинить тем же условиям, на границе области, которым должно удовлетворять U. Например, если на поверхности S, окружающей объём и содержащей все источники , то для , тогда (25) примет вид:
(26)
Если на S и ,
то В (25) тоже равен 0.
Строго можно показать, что если U и G удовлетворяют одинаковым граничным условиям, то
.
Таким образом, применяя формулу (26), мы можем считать, что задача решена если она решена для –источника, то есть известна функция Грина данной задачи.
Функция Грина в вакууме при наложении условия излучения имеет вид
(27)
Физический смысл (26) очень прост. Функция Грина (27) – это поле точечного источника. Поле произвольного распределения источников находится как суперпозиция полей точечных источников.
Векторная функция Грина
Если линейное дифференциальное уравнение векторное, то при отыскании функции Грина в правую часть, кроме -функции, необходимо ввести ещё единичный вектор. Решение уравнения представляет собой вектор с тремя составляющими. Чтобы охватить все случаи, необходимо решить уравнение при трёх ориентациях единичного вектора в правой части. Таким образом, функция Грина в векторном случае определяется девятью скалярными величинами , которые образуют тензор второго ранга. По форме запись решения через тензорную функцию Грина остаётся такой же, как и в скалярном случае, но расшифровывается эта запись подобно тому, как расшифровывается тензорная связь между векторами и (4),(5).
Исключение из этого правила составляет функция Грина уравнения для векторного потенциала или вектора Герца
Нетрудно показать, что вектор и в этом случае всюду параллельны. Поэтому решение его может быть записано через скалярную, а не векторную функцию Грина.
где G определяется формулой (27).
Простота этой формулы обусловила широкое применение векторных потенциалов и вектора Герца . В неоднородных средах функция Грина векторного потенциала представляет собой тензор. И в этом случае применение векторного потенциала не упрощает решение задачи.
Функция Грина уравнения Максвелла.
Функция Грина уравнения Максвелла можно вводить двумя способами в зависимости от того, в какое из двух уравнений вводится -источник. Запишем два первых из уравнений Максвелла
1. |
2. |
В первом варианте вводим функцию Грина элементарного электрического диполя, обозначим, создаваемое при этом поле обозначим ( и ). Заменим в первом уравнении на -функцию, умноженную на единичный вектор .
3. |
4. |
Во втором варианте функцию Грина элементарного магнитного диполя, создаваемое при этом поле обозначим (,). Заменим во втором уравнении на -функцию, умноженную на единичный вектор .
Здесь принято вносить -источник с противоположным относительно первого случая знаком.
Применим лемму Лоренца в интегральной форме, то есть векторный вариант формулы Грина.
Для этого:
Умножим скалярно 4 на , 1 на , вычтем из 4–1
(28)
Умножим скалярно 2 на , 3 на , вычтем из 2–3
(29)
Вычтем из выражения ( 28 ) выражение ( 29 ).
При этом используем векторное тождество
.
В этом случае
;
и окончательно
.
Проинтегрируем по объёму, включающему точку
Аналогично можно получить:
Если интегрирование ведётся по бесконечному объёму, то интегралы
так как в бесконечности поле равно нулю. То же самое будет при интегрировании по бесконечному объёму, включающему металлические тела.