И.К. Лифанов - Особые интегральные уравнения и методы их численного решения (1109002), страница 5
Текст из файла (страница 5)
где ‑
периодическая обобщенная функция.
Уравнение (28) имеет единственное решение для любой обобщенной функции , а уравнения (29) и (30) имеют решение с точностью до константы и условием их разрешимости является равенство
где под значением интеграла от по [0,2
] понимается величина
. Причем формулы обращения для этих уравнений в пространстве обобщенных функций те же, что и в пространстве
и имеют следующий вид (см. п.1). Для уравнения (28) имеем
для уравнения (29) имеем
для уравнения (30) имеем
где в (33) и (34) для С выполняется равенство
так как для разрешимости этих уравнений должно выполняться равенство (31).
Однако для построения вычислительных методов решения уравнений
(28)-(30) надо уметь оценивать близость решений этих уравнений по близости их правых частей. Для этого заметим, что если взять произвольный ряд из экспонент (обобщенную функцию) и применить к нему один из операторов L, S или H, то, получим аналогичный ряд у которого коэффициенты отличаются от коэффициентов исходного ряда множителем порядка ,
=
. Поэтому естественно возникает идея использования пространств
типа пространств Соболева обобщенных функций. Пусть
-
периодическая обобщенная функция на
, тогда
где под понимается значение обобщенной функции (линейного функционала)
на функции
=
. Обозначим через
множество таких обобщенных функций
, что
т.е. функция принадлежит пространству
.
Введем теперь для функций из скалярное произведение по формуле
Так введенное скалярное произведение делает гильбертовым пространством причем
- пространство функций с интегрируемым квадратом модуля на
Теперь рассмотрим характеристические особые интеральные
уравнения на отрезке. Для этого опять вначале напомним спектральные соотношения, полученные в п.1
где ,
,
,
. Теперь характеристические особые интегральные уравнения на отрезке записываются в виде (см. п.1)
Рассмотрим Гильбертовы пространства построенные в п.2 и возьмем пространства обобщенных функций
. Пусть теперь
является обобщенной функцией на
, т.е.
, тогда
где под понимается значение обобщенной функции (линейного функционала)
на функции
. Обозначим через
множество таких обобщенных функций
на
, что
т.е. функция принадлежит пространству
. Введем теперь для функций из
скалярное произведение по формуле
Так введенное скалярное произведение делает гильбертовым пространством причем
- пространство функций с интегрируемым квадратом модуля на [-1,1]. Теперь для любой обобщенной функции
,
, k=1-4, будем полагать, что значение
,
,
,
получается по правилу перемены местами знаков интеграла и суммы в этих операторах. Теперь соотношения (39)-(45) показывают следующее. Оператор
взаимнооднозначно отображает пространство обобщенных функций
на себя и пространство
на пространство
; оператор
отображает
на
и ядром отображения является множество констант и с таким же свойством
на
; оператор
отображает
в
и для функций
из образа выполняется соотношение
и с таким же свойством в
; оператор
отображает взаимнооднозначно
на
и
на
; оператор
отображает взаимнооднозначно
на
и
на
; и , наконец, оператор
отображает взаимнооднозначно
на
и
на
.
Важно, что теперь можно говорить о близости правых частей и соответсвующих решений, являющихся обобщенными функциями.
Теперь сделаем следующее важное замечание. Для уравнений (46)-(51) на базисных элементах соответсвующих пространств справедливы формулы обращения, полученные в п.1. Поэтому, в силу определения значений операторов ,
,
,
на обобщенных функциях соответствующих пространств, для уравнений (46)-(51) справедливы следующие формулы обращения.
Для уравнения (46) имеем
для уравнения (47) имеем
для уравнения (48) имеем
при выполнении условия (55),
для уравнения (49) имеем
для уравнения (50) имеем
для уравнения (51) имеем
Используя проведенные выше построения в пространствах обобщенных функций получим теперь некоторые точные решения характеристических особых интегральных уравнений первого рода в этих пространствах.
Например, задача бесциркуляционного обтекания кругового цилиндра потоком идеальной несжимаемой жидкости при наличии на нем устройства отсоса внешнего потока сводится к решению уравнения [3]
В равенстве (63) под функцией понимаем
периодическую дельта функцию, определяемую равенствами:
=0, при
, и
=+
, при
, и для любой функции
,
, выполняется равенство
. Из этого вытекает, что если
обобщенная функция, то так же имеем равенство
. Теперь из (33) получаем, что решением уравнения (63) является фунция
Таким образом, получаем формулу[3]
Этот же результат получается [3] если обобщенные функции и
представить рядами Фурье по функциям
и подставить эти ряды в уравнение (63). Действительно, имеем
Теперь из (1.4) следует
Используя опять (1.4) получим
С другой стороны, по определению функци полусаем
Сравнивая (67) и (68), видим справедливость равенства (65).
Из представлений рядами Фурье (66) и (68) следует, что обе функции и
принадлежат любому пространству
для
. Заметим так же следующее. Из соотношений (22), (23) следует, что равенство (65) можно почленно интегрировать и поэтому верно соотношение
где ,
;
,
;
,
, и
, а для q=0 имеем
,
;
,
;
,
. Так как свободный член ряда Фурье для функции справа в (69) должен быть равен нулю, то
. Отметим, что функция
является первообразной для функции
, обращающаяся в нуль при
.
Решение задачи о расчете входного сопротивления тонкой проволочной антенны при запитке антенны источником тока, повлекло за собой рассмотрение гиперсингулярного интегрального уравнения на отрезке, в правой части которого стоит функция с особенностью типа 1/x внутри области поиска решения [5]. Аналогом этого уравнения в периодическом случае будет уравнение
В силу формулы обращения (34) для рассматриваемого уравнения и формулы (69) получаем, что решнием уравнения (70) будет функция
для которой выполняется равенство (35).
Для приближенного метода решения уравнения (63) методом дискретных вихрей важно, что функция является пределом в смысле определения 4 последовательности
периодических функций
=
,
;
=0,
,
, т.е.
.
Теперь рассмотрим уравнение (47) на отрезке. Возьмем в его правой части функцию ,
, тогда в силу формулы (57) получим, что решением этого уравнения будет функция (С=0)
Подставляя теперь последнюю функцию в уравнение (47) получаем интересное равенство
Наконец рассмотрим еще частный случай уравнения (51)
который встречается в теории антенн. Оператор, стоящий слева в (74) отображает пространство на себя. Поэтому решение этого уравнения и его правую часть рассматриваем как элементы этого пространства. Как и при обращении уравненя (2.45) сведем уравнение (74) к кравнению
которое эквивалентно уравнению (74) при выполнении условия
где . Из формулы обращения (2.48) для уравнения (2.41) получаем, что общее решение уравнения (75) дается формулой
В силу равенства (76) получаем, С=1. Поэтому, так как , имеем
где =0, при
,
=
, при x=q,
=1, при
. Функция
является первообразной для функции
на отрезке [-1,1],
. Таким образом, решение уравнения (74) дается формулой
П.4. Метод дискретных вихрей численного решения характеристических особых интегральных уравнений.
Вначале рассмотрим численное решение характеристических особых интегральных уравнений на отрезке. Потом будем рассматривать эти уравнения в периодическом случае. Поскольку в периодическом случае при рассмотрении метода дискретных вихрей удобно будет ввести в рассмотрение точки комплексной плоскости, которые обычно принято обознать буквой [2], то и при рассмотрении метода дискретных вихрей для особых интегральных уравнений на отрезке будем использовать вместо переменной
переменную
. В классической теории сингулярных интегральных уравнений уравнения (3.47)-(3.50) записывают единным образом в виде