МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ
им. М.В. ЛОМОНОСОВА
Факультет вычислительной математики и кибернетики
И.К. Лифанов
ОСОБЫЕ ИНТЕГРАЛЬНЫЕ УРАВНЕНИЯ
И МЕТОДЫ
ИХ ЧИСЛЕННОГО РЕШЕНИЯ
Учебное пособие
Издано при поддержке образовательной программы
«Формирование системы инновационного образования в МГУ»
МАКС
ПРЕССА
МОСКВА-2006
УДК 517.977
ББК 22.193
Л64
Печатается по решению редакционно-издательского совета
факультета вычислительной математики и кибернетики ШУ им. М.В. Ломоносова
Лифанов И.К.
Л64 Особые интегральные уравнения и методы их численного
решения: Учебное пособие. - М.: Издательский отдел факуль-
тета ВМиК МГУ им. М.В. Ломоносова (лицензия ИД N 05899
от 24.09.2001 г.); МАКС Пресс, 2006. - 70 с.
ISBN 5-89407-291-3
ISBN-13: 978-5-317-01843-6
ISBN-10: 5-317-01843-9
В спецкурсе излагается теория одномерных особых интегральных
уравнений первого рода в периодическом случае и на отрезке с помощью
элементарных средств без использования достаточно сложного аппарата
краевых задач теории аналитических функций. Излагается также метод
дискретных вихрей для численного решения этих уравнений.
УДК 517.977
ББК 22.193
ISBN 5-89407-291-3
ISBN-13: 978-5-317-01843-6
ISBN-10: 5-317-01843-9
© Факультет вычислительной математики
и кибернетики МГУ им. М.В. Ломоносова, 2006
© Лифанов И.К., 2006
Кратко об авторе.
Иван Кузьмич Лифанов - доктор физико-математических наук,
профессор, заслуженный деятель науки Российской Федерации, почетный
профессор Военно-воздушной инженерной академии им. Н.Е. Жуковского,
почетный доктор Одесской национальной академии связи им. А.С. Попова.
Работает на кафедре вычислительных технологий и моделирования
факультета вычислительной математики и кибернетики Московского
государственного университета им. М.В. Ломоносова, заведует кафедрой
высшей математики в Военно-воздушной инженерной академии им. Н.Е.
Жуковского и является главным научным сотрудником Института
вычислительной математики РАН.
Его научные интересы: теоретические исследования и численные
методы для особых интегральных уравнений, математическое моделирование
в аэродинамике, электродинамике, дифракции волн.
Имеет более 290 научных публикаций, в том числе 11 монографий из
которых 3 на английском языке и одна на украинском.
Более 20 лет совместно с профессором Е.В. Захаровым руководит
научно-исследовательским семинаром «Интегральные уравнения и их
приложения» на факультете ВМиК в МГУ.
Один из организаторов Международных симпозиумов «Методы
дискретных особенностей в задачах математической физики», которые
регулярно проводятся уже 23 года.
Председатель докторского диссертационного совета в ВВИА им. Н.Е.
Жуковского и член такого совета в ИВМ РАН.
Заместитель председателя научно-методического совета по математике
Минобрнауки РФ.
Член редколлегий ряда научных журналов.
Член Национального комитета по теоретической и прикладной
механики РАН.
В этом спецкурсе мы хотим изложить теорию одномерных особых
интегральных уравнений первого рода в периодическом случае и на отрезке с
помощью элементарных средств без использования достаточно сложного
аппарата краевых задач теории аналитических функций. Потом изложить для
этих уравнений метод дискретных вихрей их численного решения.
Текст разбит на четыре пункта. В каждом пункте нумерация сквозная.
Если ссылка идет на формулу из другого пункта, то указывается номер
пункта и формулы, а если из того же пункта, то указывается только номер
формулы.
Содержание дано в виде вопросов для подготовки к сдаче спецкурса. К
каждому вопросу даны страницы в пособии, где излагается этот вопрос, а в
круглых скобках указана литература, где можно посмотреть этот вопрос
более подробно. Сами вопросы выделены в тексте жирным шрифтом.
п.1. Особые интегралы и уравнения с ними в периодическом случае.
Рассмотрим вначале интегральный оператор с ядром Гильберта
ЭДЖ) = ~- \ctg^-^g(<p)dg>	(1.1)
хЛ' 0	2
который по определению полагаем равным
X |c/g g(<p)d<p = lim( j ctg g{(p)d(p + f ctg g(<p)d<p).	(1.2)
2л '	2	In*	2	{	2
Так как функция ctg^ нечетна на отрезке [-я, я], то справедливо равенство
jctg	dq, =: о, (pQ е [0,2я].	(1.3)
о
Можно доказать, что оператор S(g(^),p0) определен для любой 2я
периодической функции g(^), удовлетворяющей условию Гельдера степени
а,0<а<1, на отрезке [0,2л], т.е. для любых р, и (р2 из отрезка [0,2 л]
выполняется соотношение |g(^2) - g(^ )| £	- Ф\ Г , где А>0 некоторая
константа.
Утверждение Ы.Справедливы следующие спектральные
соотношения для интегрального оператора с ядром Гильберта
J__ jCfgein<pd(p - -i• sign(ri)eln9^,n = 0,+l,±2,...
(1-4)
где sign (и) =
0,и = 0,
1 2x ф ф
— J<^g	~ cos ntpdcp = sin n<p0, n = 0,1,2,...
1 2n	~
— fc/g——sin =-cos«^0,и = 1,2,... .
2л 0J	2
Действительно, при целом положительном п имеем
ctg^^ =	-(?*)"] =
= i(e'f + е’*	+...++	) =
= i(e‘"v +	+... + 2е',’е'{"~')л + е”й’).
(1-5)
(1-6)
-е
(1-7)
Переходя в равенстве (7) к комплексно сопряженным величинам, получаем
(е-»’ -е-”"')= =-,(e-"F +2е-'("-1”’е-|л +... + 2е"”е’'("Ч)л +e"”ft).	(1.8)
Таким образом, для отрицательных целых чисел п имеем
~е,п^) =-i(e,n,p +2e,(n~')^pe,9’<, +,., + 2е19>е‘{п^ +е,п^°).
(1.9)
4
Интегрируя тождества (7) и (8) и меняя местами аргументы под котангенсом
получаем равенство (4). Для получения равенства (5) надо взять равенство (4)
при п=1 и п==-1, сложить их и разделить на 2. Аналогичным образом
получается равенство (6).
Так как оператор Гильберта обладает свойством линейности, т.е.
S{af (<р) + bg(<p), ) = aS(f (ср), ?Q) + bS (g(q>\ (pQ),	(1.10)
то знак конечной суммы для функций = 0,±1,±2,..., и знак оператора
Гильберта можно менять местами. Поэтому оператор Гильберта переводит
тригонометрические полиномы в тригонометрические полиномы. В силу
этого этот оператор можно распространить на пространство
L2 комплекснозначных 2я периодических функций на [0, 2я] с
интегрируемым квадратом модуля. В этом пространстве вводится норма
|/|2 = $f , f)> где под знаком квадратного корня стоит скалярный квадрат
функции £ Скалярное произведение функций f и g в этом пространстве
2я
вводится следующим образом: (/,g)- jf(<p)g(p№p- Указанное пространство
о
L2 является Гильбертовым и поэтому в нем имеется базис, состоящий из
функций y/n=-^=einv, и = 0,±1,±2,.... Пусть теперь /(<?>)- 2лг периодическая
функция из L2, тогда
/W = £./('>>„, /(”)=(/(?’).¥'»)=-Т==
лег	л/2лг 0
(1.11)
причем выполняется соотношение
I4=(Z|/(«)|V<«>.	(1.12)
neZl I
Для этих функций будем по определению полагать
ЯЛЛРо) =	\^g~^f(4>)d4> =
2тг 0	2
X	lim £ f{n)\ctg^-~^-4/n(q>)d<p =
о	лег	H--N Q
1““ S/(«)(~“'«"(»)>„(Po) = Hm	^/.(«>„(f»o)-	(1-13)
n--N	П^-N	п-~а>
Последнее равенство справедливо в силу того, что последний ряд сходится в
L2. Действительно, имеем
Zl/.ool = £ |л”)| <о°-	О-14)
neZl I heZ,h*oI I
Итак, по определению имеем
_L	=	(1.15)
2 «ez	П6г
Таким образом, если функция g(<p) принадлежит пространству Z2, то
функция /(^0), определяемая равенством
5
/(₽o)-^- fag~-^g(^)dp	(1.16)
2zr 0	2
так же принадлежит пространству L2.
Утверждение 1.2. При п * 0, п е Z, справедливо равенство
е‘Ю =	(1-17)
2тг 0	2
Действительно
-у-	-isign(n)~ fag'fa~,-e'"n‘dp0 =
2я- J 2	2л- 0J 2
isign(n)(~isign(n))em9 = emtp,
Теперь заметим, что в равенствах (4) и (17) вместо функции ew”, «^0,«gZ,
можно поставить функцию к О).
Обратимся теперь к характеристическому интегральному
уравнению первого рода с ядром Гильберта. Так называется уравнение
вида
у- fag	g{<p)d<p = f(<p0), <»„ е [0,2 л].	(1.18)
2лг 0	2
Утверждение 1.3. Пусть f(<pQ) принадлежит пространству L2. Тогда,
для того чтобы уравнение (18) имело решение из L2 необходимо выполнение
равенства
]/(^г0.	(1.19)
О
Действительно, пусть функция g(p), принадлежащая пространству 12,
является решением уравнения (18). В этом случае должно выполняться
равенство
2л1 0	2	2л J 2 neZ
fag—^e^dp = Yt(,-sign(n))g(n)e“'n = /(<»„).	(1.20)
neZ 2ЛГ q 2	ие2
Из последнего равенства следует равенство (19).
Пусть теперь функция f(<p0) в (18) принадлежит пространству L2 и для
нее выполняется равенство (19). Тогда в силу равенства (17) и определения
интеграла с ядром Гильберта от функции из L2 получаем, что функция g(<p),
определяемая равенством
,	1 2? , Ч> -Ч>0 ,, . ,	1 гг t Ч> ~Ч>« v " .	,
2л 0J 2	2л 0J 2
^l{-sign(ifj}f(n)e""r	(1-21)
6
является решением уравнения (18). Ясно, что в силу равенство (3) функция
g(p)+C для произвольной константы С так же является решением уравнения
(18).
Приведенные выше рассуждения показывают справедливость
следующей теоремы.
Теорема 1.1. Пусть в уравнении (18) функция /(^0) принадлежит
пространству Z2 и для нее выполняется равенство (19). Тогда общее решение
этого уравнения, так же принадлежащее L2, дается формулой
8W=~~1- р8	— f(<Po )d<pQ + С.	(1.22)
2л- о 2
Теперь рассмотрим интегральный оператор с логарифмической
особенностью. Так будем называть оператор
Ц£(<Р),<Р0) = - flnlsin^™^(fp)d<p, <pQ е [0,2лг].	(1.23)
о I 2 Г
Если g(cp) является непрерывной функцией на [0,2л-], то интеграл в (23) имеет
логарифмическую особенность при ср = и поэтому определен. Обозначим
через G(<p) некоторую первообразную функции Тогда интегрируя по
частям получим равенство
~ fln|sin^-^|g(^^	1о*+У"
л" J I * \ лг I 2 I	Zn ’	2
Так как первое слагаемое справа в последней формуле равно нулю, в силу
периодичности функции G(<p) (предполагается, что g(<p) не является
константой), то получаем равенство
” jln|sin	fctg G(<p)d<p, <р0 е [0,2л-].	(1.24)
Л о I 2 I 2Л о 2
Интеграл справа в (24) определен, так как функция G(<p) удовлетворяет
условию Гельдера при « = 1, т.е. условию Липшица. Поэтому равенство (24)
справедливо. Возьмем g(^)=ew<9, тогда можно взять G(tp)=(iriy~}eintp. Теперь в
силу равенства (4) получаем соотношение
-	п = ±1,±2„...	(1.25)
» i I 2 | п
Наконец, осталось вычислить значение 1(1, р0). Это значение имеется в
справочниках (см. [1] стр. 107), но мы здесь дадим вычисление этого
значения, приведенное в учебном пособии Гандель Ю.В. Введение в методы
вычисления сингулярных и гиперсингулярных интегралов, Харьков-2001, 92
с. В силу периодичности функции Injsin с периодом 2л- получаем
1 I	/9|	1 я	i	zj|	1*1	/?|
- flnsinH<70 = — flnsin-Lz0 + - fin cos-dO =
о I	2I	/r J	I	2Г I	2|
7
— J(ln|sin 0| - In 2)dO - — Jlnjsin -1/0 - In 2.	(1.26)
яг о	2яг - |	2|
Таким образом, получаем
- flnjsin^^L = -21n2.	(1.27)
* 0 I 2 I
В силу формулы (24) и формул (5), (6) получаем справедливость
соотношений
I	ф /Л|	J
— IIn sin ——— cos ntpdtp = — cos nq>, П= 1,2,...
«; I 2 |	n
(1.28)
— fIn sin ——- sin ncpd(p = - - sin nq>, n= 1,2,...	(1.29)
* о I 2 I	n
Наконец, отметим, что из соотношений (4) и (25), (27) следует
справедливость следующей теоремы.
Теорема 1.2. Пусть функция /(^0) принадлежит пространству L29 т.е.
Л^о) е й^1, на [0,2яг], и функция g(p) удовлетворяет равенству
- finjsin^-^lg(p)<ty,= /(p0), % е [0,2л-].	(1.30)
«J I 2 |
Тогда §((/>) е 12 на [0,2л-] и справедливо следующее равенство
4~~ jtalsin—^|g(?>)af<? = - М-Inlsin<^-^g(jp)d(p = /'(p»),^ e [0,2л-].	(1.31)
‘/«’о * о I 2 Г Я }d<p0 I 2 Г
Теперь рассмотрим гнперсингулярный оператор. Так будем называть
оператор вида
H(g(^),<»0) = 7-е[0,2л].	(1.32)
2
Вначале дадим определение 4я Я(1,^о), что будем понимать в следующем
смысле:
2f d0
»Jsin’^£
2
= lim (
- dO л < s
[0,2л-]\О(^л) sjn 2 У0---	Z
2
= КтГ2ctg - 2ctg + 2ctg в,,~П + 2ctg - 4ctg =
z z	z	z zj
= 2ctg^ - я) - 2«g^y = 0 •
Теперь для 2л- H(g(^),^0) по определению будем полагать
(1-33)
f gWdf> = Hm г g(4>)dq>
о sin2	е‘>0 [О,2»)\0(«о.«) sin2
2 L	2
4g((?0)^g|
(1.34)
Справедливы следующие утверждения.
8
Утверждение 1.4. Оператор	определен для любой функции
g(<p), принадлежащей классу Нх(а), т.е. g’(p) удовлетворяет условию
Гельдера степени а,0 < а < 1, на отрезке [0,2яг].
Действительно, имеем
2Г	= 2r(gW-g(^o)-g'(^o)(<g-PoW +	Л d(p +
oJsin2^ 0J	sin2^-z£	°
2	2	2
g(yo)J^--^,Poe[Q,2g].	(1.35)
»sm2«
2
Первое слагаемое справа в (35) существует как абсолютно сходящийся
интеграл, так как из того, что g’(p) е Нх(а), имеем
jg(^)-g(^o)“g (^o)(^-^o)|	•	0-36)
Второе слагаемое равно нулю в силу (33), третье слагаемое существует как
сингулярный интеграл.
Утверждение 1.5. Пусть функция g(^) - 2л- -периодическая и
принадлежит классу Нх(а) на [0,2тг]. Тогда для интеграла (34) справедлива
следующая формула интегрирования по частям:
= ^yctg^dg,.	(1.37)
2
Действительно, первообразной для функции 1 /sin2 —- — является
функция 2crg^-^. Поэтому, интегрируя по частям, получаем
«^гЧ>о-<1>	2	J	2
МП ----
2
Первое слагаемое справа равно нулю в силу периодичности функции g(#>).
Итак, для оператора Я(я(<р),р0) имеем
#(g(P),<M = T~ j =--L	g>Q е[0,2тг].	(1.38)
4л~ »sin2<?0~^	2*o	2
2
Теперь из формул (4)-(6) и (38) получаем следующие спектральные
соотношения
J_ Г—£—- -nsign(n)e,f,*°, п = 0,±1,±2,..., (р. е [0,2яг],	(1.39)
4<in2£oz£
2
9
1 2? sinwcx/® . л~ глл ,
— J-----= -Hsmn<pQin = 0,1,2,..., <pQ e[0,2лг],	(1.40)
sm __
I -C-°~"	= -ncos«^0,n - 0,1,2,..., <p0 g [0,2лг].	0-41)
Sm 2
Отметим, что из равенства (38) следует
Утверждение 1.6. Пусть функция /(^0) принадлежит пространству Z2
и функция g(<p) удовлетворяет равенству
г- /	=	е [0’2?г] ’	(1 -42)
2
тогда функция g(^) принадлежит пространству W2.
А из равенств (4) и (39) следует
Теорема 1.3. Пусть функция f \<рй) принадлежит пространству L2, т.е.
f(<p0) g W2i на [0,2я], и функция g(<p) удовлетворяет равенству
~~ fagg(<p)d<p = f{<p0[0,2я].	(1.43)
2# 0	2
Тогда справедливо равенство
fag~r^g(^)d9>=~ \~-ctg^~-g(<p)d(i>^f '((pn), р0 е[0,2я-] (1.44)
«0оi 2	2* о «00	2
Действительно, из равенств (4) и (39) имеем
“V- fag-,‘,—e‘",,d<p== -^-(-isign(n)emn‘)=nsign{n)e^ =
d<?<> 2jt 0J 2	d<p„
_ 1	1 2Н 0^0е^	[02jr] и = 0±1±2 .	(1.45)
2
Из последнего равенства следует утверждение теоремы, так как если
/(^0)е^2, то она и ее производная представляются сходящимися рядами
Фурье.
Обратимся теперь к рассмотрению интегральных уравнений. Вначале
рассмотрим характеристическое интегральное уравнение первого рода с
логарифмической особенностью, т.е. уравнение вида
- flnlsin^~^g{(p)d(p = f{<pQ) t <p0 G [0,2л-] .	(1.46)
Из спектральных соотношений (25), (27) следует, что любая функция
g(^)eZ2 является решением уравнения (46) для некоторой
функции /(^0) g W2. Действительно, пусть g(^) g L2, т.е.
=	=	(1.47)
neZ	yZft 0
10
где у/п = ~^=ein,p, причем выполняется соотношение
л/2яг
kL=(Z|g(”)| У <°°-
neZl I
Тогда
- flnlsin^—-|g(p)«ty = - jlnlsin^—^|^g(w>rt(^>/^ =
я i I 2 Г я J I 2 |S
(1.48)
£g(n)~ flnsin^—^L(^)^= ^g(n)(--5ign(«>B(^0)) + g(0)(-21n2>0(<?>0) =
S о I 2 I	n
f(<Po),<Po e[0,2^].	(1.49)
Из последнего равенства следует, что /(^0) g W}. Из него же и соотношений
(25) и (39) следует справедливость следующей теоремы.
Теорема 1.4. Для любой функции /(^0) уравнение (46) имеет
единственное решение, которое задается равенством
g(p) = J- J—, <Р е [0,2яг],	(1.50)
4*- Agin^-^O 2In20
2
и принадлежит пространству L2.
Теперь рассмотрим характеристическое гиперсингулярное
интегральное уравнение первого рода, т.е. уравнение вида
1 I g(yp	ft,e[Q)2jr].	(1.51)
4/г * • г (Pq~(P
0 sin ——~
2
Из спектральных соотношений (39) следует, что любая функция g(^)elT21
является решением уравнения (46) для некоторой функции/(^0) е£2.
Действительно, пусть g(^) g W2, т.е. g’(^) е L2. Последнее означает, что
g' (ф)=(X w)=X w=Xin w	0 -52)
neZ	nt=Z	neZ
и последний ряд является сходящимся в L2. Отсюда имеем
1 f gWd<P = 1 Г ________________= у (п. 1 Г	=
osin2^? 4*0J Sin2^»z£	4^osin2^>Z?
2	2	2
^(-nsignin^gin^^y^f^).	(1.53)
neZ
Из формулы (52) следует, что функция /(^0)в формуле (53) принадлежит
пространству L2. Последние рассуждения показывают справедливость
теоремы.
11
Теорема 1.5. Пусть в уравнении (51) функция /(^0)gZ2. Тогда
условием существования решения этого уравнения является равенство (19),
т.е.
о
его решение определено с точностью до константы и задается формулой
g(<?) = - flnlsin^——l/(p0)^0 + ~ •	(1 -54)
Ч I 2 Г 2яг
Действительно, если /(^0)е£2, то из формулы (25) следует, что при
п ф 0 из равенства (54) получаем
g(”) = -	/(и) .	(1.55)
nsjgn(n)
Подставляя значения g(w) из формулы (55) в формулу (53) получаем
справедливость теоремы.
П.2. Особые интегралы и уравнения с ними на отрезке.
Как и для особых интегралов в периодическом случае рассмотрение
начнем с оператора Коши на отрезке [-1,1],
S(g(x),x0)=-	e(-1,1),	(2.1)
Я' Д х0 - X
который по определению полагаем равным
1 [gQfrfc = l|im[ g(*)^ + [
J xo-x X£xo-xj
При g(x)=l получаем
1 r dx
7rjxo -X
=	Xo G (—1,1) .
я 1- x0
(2.2)
(23)
Из (2) и (3) следует, что оператор £(g(x),x0) определен для любой функции
g(x), которая удовлетворяет условию Гельдера степени а,0 < а < 1, на отрезке
[-1,1], т.е. g(x)f=Ha. Более того, он определен для любой функции g(x),
которая может быть представлена в виде произведения функций 2(х) g На на
[-1,1 ] и функции р(х) = (1 + х)’и (1 - х) р.
Для того, чтобы рассмотреть спектральные соотношения для оператора
5(g(x),x0) вначале обратимся опять к интегралу Гильберта. Справедливо
следующее утверждение.
Утверждение 2.1. Если функция g(<p)eHa на [-я, я] и является
нечетной, то выполняется равенство
<2-4>
-1	2	0Jcosp0-cosp
a ecnng(^) является четной, то
pg^-^g(?)<ty== -2sin^0 J-	, <p0 е[-я,я]. (2.5)
4	2	'cos<pn-cos®
12
Доказательство проведем для равенства (4), так как для равенства (5)
оно аналогично. Имеем
pg g((p)d<p = pg - g(q>)d<p + fc/gg{(p)d(p. (2.6)
-я 2	-к 2	о 2
Сделаем в первом интеграле справа в (6) замену переменной <р = ~0,
поменяем местами пределы интегрирования и обозначим опять О через (р.
Получим
|c/g g(g>)dq>- J(cl£	- - ctg ^~)g(<p)dg> -
-n 1	0	2	2
f___-______г________2_____2_g(^=_2 f
JCOS^o -cos<p
(2.7)
sin^sin^
О
2	2
Полагая теперь в (5) g(^>)=cos«^ и воспользовавшись равенством (1.5),
получим
1 7 cos npdp _ sinH^0
лJcos^0-cos^9 sin«p0 ’	’ ’
Взяв теперь в (4) g(^)=sin«^ из (1.6) получим
(2.8)
lj™”^^=cos п=1д...	,
л ; cosp0 -cos0>
Сделаем теперь в равенствах (8), (9) замену переменной cos^=x, cos^0=x{
^^е[0л], х,х0 е[-1,1]. Тогда получим
1 Vcos(«arccosx)dr = sin(«arccosx0)
71-х2(х0-х) sin(arccosx0) ’
(2.9)
(2.10)
1 г sin(arccosx) z ч . ~
— ]........................=cos(warccosx0), П~ 1,2,...
(2.И)
Известно (см. Бейтмен Г., Эрдейн А. Высшие транцендентные функции, т. II,
М., Наука, 1966,296 с.), что функция
Гп (х) = cos(« arccos х), п=0,1,2,...	(2.12)
является многочленом степени п и называется многочленом Чебышева
первого рода, а функция
^iW=^(«arccosx)5n=b25
sin(arccosx)
Является многочленом степени п-1 и называется многочленом Чебышева
второго рода. Доказательство того, что функции Тй(х) и Un_x(x) являются
многочленами, получается следующим образом [7]. Обозначим arccosx = 0 и
напомним, что
(2.13)
13
cos n0 + isin n0 - (cos 0 + isin 0)n -
- wsn 0 + /С’ cos**-10sin 0 - C2 cos**-2 0sin2 0 -iC3 cos"-3 0 sin3 0 +...
Разделяя вещественную и мнимую части, получаем тождества
cos п0 ~ cos'* 0 - С2 cos'*-2 0(\ - cos2 0) +...	(2.15)
= С\ cos”10 - С3 cos**-3 0(1 - cos2 0) +....	(2.16)
sin0
Теперь соотношения (10) и (11) можно записать так
-	=-17^(х0), х0 е[-1,1], п=0,1,2,...	(2.17)
^-Wl-x^-x)
1	х0 6[-1,1], п=1,2,...	(2.18)
Равенства (17) и (18) при х0 равным -1 или 1 надо понимать в том смысле,
что надо сделать в этих равенствах замену cos^>=x, cos^0=x0 и
воспользоваться равенствами (8) и (9). Равенства (17) и (18) и будем называть
спектральными соотношениями для интеграла с ядром Коши на отрезке в
классах функций, обращающихся в нуль или в бесконечность на концах
отрезка.
Теперь, опираясь на соотношения (17) и (18), покажем справедливость
еще следующих спектральных соотношений для интеграла с ядром
Коши на отрезке
1	=_g (х ); х0 е[-1,1], п=0,1,2,...	(2.19)
Л 4 V1 + х х0 - X
1 |Е*й^ = Р(х)> х0 е[-1,1], п=0,1,2,...	(2.20)
ЯГ 4 Vi “X х0-х
где Ри(х) = [Тя+,(х)-Ти(х)]/(1-х), e„(x) = (7K(x)-t/„_](x), хе[-1,1].
Соотношение (19) легко следует из соотношения (17).
Соотношение (20) получается из следующей цепочки равенств
Qn(x)d* = 1 f fl + x ^fl(x)-Un 1(x)dx =
x0-x	xQ-x
sin((n + l)arccosx) sin(warccosx)
1 p/l-x2 sin(arccosx) sin(arccosx) 1 *r(sin(n + l)ff - sin Kff)sin p =
л 4 1 - x	x0 - x	л J (1 - cos0>)(cos$!>o - cos g>)
/	1 x • Ф ~ - Ф Ф	~ z L Ф
. я 2cos(w + -)©sm — 2sm -cos— . я 2cos(w + -)g?cos~
1 f----------f	2d(p=
я * cos^0 -COS
2 sin2 у (cos^0 - cos ф)
14
1 *	/	14	• z 14	-	2sin(w + — )fi>0 cos-
1 fcos(w + l)« + cosw^ , _ _sin(w + l)^0 sinn^0. v 2	2 _
— 1-----------------rfy — _ (-------------------j —---------------—---- —
cosp0 -cosp	sinp0 smp0	2 sin	cos—°
2	2
_ 2 sin(n + )^0 sin _ COs(w +1)^0 - cos n<pQ _ cos((w +1) arccos xQ) - cos(n arccos x()) _
2sjn2?i	l-cosp0	l-x0
2
x° G[_U]jп=о,1Д...	(2.21)
1 x0
Эта цепочка равенств получается следующим образом. Вначале пользуемся
определением функции (7и(х0), затем делаем в интеграле замену cos^=x,
cos^0=x0h потом пользуемся известными преобразованиями
тригонометрических выражений.
Рассмотрим теперь интегральный оператор с логарифмической
особенностью на отрезке, т.е. оператор вида
£(g(*),*o)=“ Jlnlx0 -5 х0 е [-U] •	(2.22)
Будем вначале предполагать, что функция g(x) представляется в виде
произведения функции Л(х)еНа на [-1,1] и функции р(х) = (1 + х)-и(1-х)Л
где0<у,//<1. В этом случае подынтегральная функция в (22) будет
абсолютно интегрируемой на [-1,1] для любого хое[-1,1]. При
сформулированном условии на функцию g(x) справедливо следующее
утверждение.
Утверждение 2.2. Для интеграла в (22) справедлива формула
интегрирования по частям, т.е. справедлива формула
(ln|x0 -x|g(x)ak = ln|x0 — х|Сг(х)/\4- Г- —	, х0 & (-1,1),	(2.23)
1	-1Х0-Х
где <7(х) некоторая первообразная для функции g(x) на [-1,1]
Используя формулу (23) можно доказать следующее соотношение
- f-Дт 1п1х«	6 (-U).» = * 1>2>-	(2-24)
я-Дл/1-х2	п
Действительно, возьмем в равенстве (23) g(x)=cos(«arccosх)/д/1-хт,
тогда можно взять G(x)=- isin(n arccosх). При этом получаем G(-1)=G(1)=O.
п
Поэтому получаем следующее равенство
sin(narccosx)^
flnlx -xl dx= -- Vsin(w33X005_ 1 Г_______________sin(arccosx) _
J ° ТГ^Х7 «_1 х0~х nJ х0-х
15
_1	х0 е[-1,1], П=1,2,...
п д хс - х	п
которое и доказывает равенство (24).
Отдельно докажем равенство
1 1	1
— |-^====г ln|xQ ~ х|с& = - In 2, х0 е (-1,1), и = 0.
Действительно, имеем
(2.25)
(2.26)
- Г _I-- -lnjx0 -x|dr=— [ln|cos^0	[inh sin —- sin + -L/ff =*
Ячл/1-х2	^0	*0 I 2	2 I
ln2 + — flnlsin^\d(p + — [ln|sin-0-^tfy?=ln2 + — flnlsin———Ь^=-1п2, x0 e (-1,1).
^0 I 2 I I 2	*-11	2 I
В последней цепочке равенств вначале сделали замену переменной cos^=x,
cosp0=x0J затем воспользовались формулой разности косинусов, потом
формулой (1.27).
Таким образом, равенства (24) и (26) дают спектральные соотношения
для оператора с логарифмической особенностью.
Теперь рассмотрим гиперсингулярный оператор. Так будем
называть оператор вида
tf(g(x),x0)=lj-£^*	ое(-1,1).	(2.27)
Вначале дадим определение яН(1,х0), что будем понимать в следующем
смысле:
f—=Iim[ f —	- + Г —^-5- - -]=---— + -—, х„ е (-1,1). (2.28)
-i(x0-x) J (х0-х)	J.(xo-x)2 е х0+1 х„-1
Теперь для п H(g(x),x0) по определению будем полагать
= lim(7J _ Wj e (_U). (2,29)
.Jcxo-x)2 '-*0 д (x0-x)2 Jt(x0-x)2	5
Докажем теперь следующие два утверждения.
Утверждение 2.3. Интеграл (29) существует для любой функции
g(x) е Н1а на отрезке [а,6], т.е. g’(x) е На классу Гёльдера степени а на [а, />].
Действительно, преобразуем формулу (29) следующим образом:
V g(x)^ = lim г	gW ~ g(^o) - g'(x0 )(* - х0 )dx
J fv —	s—>0	*
-11ло л/	L(-l.lIM>(xo,e)
f	f .*--«32.
[~1,1]\0(ло,г) (X Xo)	[-l.lJXOfxo.e) X0 X £
g(x) - g(*0) - g\x0 )(x - x0 )dx
(х-Х0):
— lim	[
₽->о	J
|Д~1ДО(х0,г)
(х-х0):
16
+ g(x0)lim
I
J-lJMw)
dx
(x-x0)2
= A(*o) + g(*o)A(*o) -g'(*o)M*o) •
(2.30)
В формуле (30) интеграл /Дхо) существует как несобственный
интеграл, так как g(x) е Н1а , и поэтому
|g (*) - g(x0) - g'(x0)(x - х0)| < А |х - х0|1+*.	(2.31)
Интеграл /2(х0) существует в силу формулы (28), а интеграл 73(х0)
существует в смысле главного значения по Коши, см. (3). Утверждение 3
доказано.
Утверждение 2.4. Пусть функция g(x) е Н{а на отрезке [a,Z>]. Тогда для
интеграла (29) справедлива следующая формула интегрирования по частям
f g(x)<& = g(l) g(-l) Vg'Qfrfr	32)
-1 (Х ““ Х0 )	Х0 “ 1 Х0 + 1 -1	— Х
Действительно, выполнив в формуле (29) интегрирование по частям,
получим:
f _ lim g(*) ° £ , 8(x) 1 _ Г g\x)dx 2g(x0)
_l(x0 ’X)	C->0^X0- X_1	X0~~X Xo+e [-l,i]\O(x0,£) X0~X s
_ g(l) g(-l) | Пп1Гgfa -g) | g(*o +g) 2g(x0)>|
xo -1 xo +1	£	£	£ J
-lim[ J	=	(2.33)
e	) Xo X ) X0 1 X0 +1
Интеграл /3(х0) существует в смысле главного значения по Коши. Так
как g(x) с #[(«), то
|g(x0+s) + g(x0-₽)-2g(x0)| =
= |[g(*o + «)-g(x0)]-[g(x0)-g(xl, -₽)]*2АЕ,+а
и поэтому предел 74(х0) равен нулю для любого хое (а,Ь). Утверждение 4
доказано.
Замечание 2.1. Если g(a) = g(z>) = o и g(x)eHla на [a,Z>], то формула
(32) получает вид:
г g(x)dx
д(*-*о)2
хое(-1,1).
xQ-x
(2.35)
Замечание 2.2. Формула (35) справедлива и в том случае, если g(-l)
=g(l)=0, a g(x)eH'a на [a,Z>], т.е. производная g'(x)eHua на [-1,1] что
означает:
g’M=_____________
' (х-аУ(Ь-хУ’
где у/(х) е На на [-1,1], a v,// < 1.
(2.36)
17
Если функция g(x) удовлетворяет требованиям, указанным в
Замечании 2.2, то будем говорить, что g(x)G^w),‘ на [-1,1].
Покажем теперь, что для интеграла (29) выполняется следующее
спектральное соотношение.
Утверждение 2.5. Справедливо следующее спектральное соотношение:
£,^^ = ^(»+l)t/.(x0), хб(-1>1);	(2.37)
v = sin((n + l)arecosx)
п sin(arccosx) ’
где U„ (х) - полином Чебышёва второго рода.
Действительно, используя представление ип (х) и равенство
sin(arccosx) = Vl-x2, можно написать:
f^l-x2{7B(x)tZr _ уsin((n+l)arccosx)(fr
J (*-х0)2 J (х~х0)2	’	(	}
Теперь в силу формулы интегрирования по частям (32), так как
sin ((и + l)arccos(l)) - sin ((и + l)arccos(-l)) = 0, получаем:
k (*-*о)2	” + 1	Vl-x2(x„-x)
=(”+1)£-г4г^=-ф+1)^^)’ ^(-’>0
vl-x2 (х0-х)
(2.39)
Последнее равенство справедливо в силу соответствующего спектрального
соотношения (17) для сингулярного интеграла на отрезке. Утверждение 5
доказано.
Теперь обратимся к формулам обращения для характеристических
уравнений первого рода с особыми интегралами. Для этого, в силу
спектральных соотношений (17>(20), (25), (26), (37), эти характеристические
уравнения запишем в следующем виде
Lp (g(x), М = - (~/===7 lnlxo " »|g(x)* = /<хо ), *o e (- M) >	(2.40)
5A(g(x),X„) = l f Д g(X)<fa=/(x0), Xo € (-1,1), Я-Д71-Х x0-x	(2.41)
^(g(x),x„) =	•* = /(x0), x0 e( 1,1), л-'	x„-x	(2.42)
Sp, (g(x),x0) = 1 f EH № = f(x ), x 6 (_U), n V1+x x0 - x	(2.43)
c z z ч 4 If 11 + X g(x)dx ^(g(x),x0)= IJ sv =/(x0), Xo e( I,I), ЯГ VI — x x0 — X	(2.44)
18
HP1 (g(x),xa) = - [Л3*2	’ х» e (-14> ’	(2-45)
Я’Д	(*0“*)
где p, (x) = -jj—, p2 (x) = p-' (x), P3 (X) = Jpi, Pt (x) = />;' (x).
71-X2	Vl + x
Для дальнеших построений напомним сначала понятие пространства L2 р
функций на отрезке [-1,1] оси ОХ, квадрат модуля которых интегрируем на
этом отрезке с весом р = р(х), т.е. функций /(х), хе [-1,1], таких, что
[p(x)\f(xf dx, где функция р = р(х)>0 почти всюду на [-1,1]. Естетсвенным
-I
образом вводится скалярное произведение функций в этом пространстве
(/(*)>#(*)) = Jp(x)/(x)g(x>j5c. После этого пространство Ь2р становится
-1
гильбертовым и поэтому в нем можно взять полную ортонормированную
систему функций ^ир, п=0,1,2,..., являющуюся базисом в нем, т.е.
= где <У,„=0 при п*т и <\„=1 при n=m. Любая
-1
функция /(x)eZ2 p представляется рядом Фурье по системе этих функций,
сходящимся по норме L2p:
/(*) = S f№Vn.p W > Я") =	(2.46)
n=0	-1
kit = Рх*)|/(*)Га=е|/(и)| <+°°-
В пространстве L2p при р -рДх) полной ортонормированной системой
является система функций т;(х) = л/2/яТи(х), п=1,2,..., и ГДх) = л/17яТ0(х), при
д = р2(х) - система функций ВД = Т2/яВД, п=0,1,2,..., при р = р3(х) -
система функций P;(x) = V17^P„(x), при р = р4(х) - система функций
(х) =	• Теперь в силу (46) видно, что если g(x) е 12 р, то
gW = Jim S sMWk'P W •	(2-47)
Jk=O
Поэтому значение ZA(g(x),x0), SA(g(x),x0), к-1,2,3,4, #A(g(x),x0) будем
определять по правилу перемены местами знаков интеграла и суммы в этих
операторах. Отсюда в силу спектральных соотношений (17)-(20), (25), (26),
(37), получаем, что оператор ZA(g(x),x0) взаимооднозначно отображает
пространство L2 p в Z2 pi, оператор SA(g(x),x0) отображает пространство
на Ь2 Рг, оператор SA(g(x),x0) отображает пространство L2 pi в Л2 д, оператор
SA(g(x),x0) отображает пространство Z2pj на £2л, оператор SA(g(x),x0)
отображает пространство ъ2р* на £2рз, оператор HA(g(x),x0) отображает
19
пространство Z2p2 на z2 ft.
Теперь можно сформулировать результаты об обращении сингулярных
интегральных уравнений (40)-(45).
Теорема 2.1. Пусть функция f(x)eL2pi. Тогда уравнение (41) имеет
решение с точность до константы принадлежащее L2pi, которое получается
по формуле
g(x)=-l	+с, х е (-1,1),	(2.48)
х-х0
и удовлетворяет условию
ljg^=C	(249)
Действительно, пусть
Я*) = SЯ«>л,Л W > f(”) = jX (x)f(x)y/nP2(x)dx	(2.50)
л«0	-1
uU = Ь(*шГ*=еМ <+®.
,Р -1	л=01 I
Тогда по (18) получаем
1 'rJl^Sf(x }dx	1
gW=_±	+c=-l --------=5----------+c=
Я’-’ X-XQ	Я’Д	X~XQ
V*	(*0)^0
л=0 Я -\ X — Xq
+c — S /(«X+1,A (*) +C.
n«0
(2.51)
Используя (51), получаем
1 f 1 g(*)^ Удл)1 [ 1
x0-x S л.’л/1-х2	x<,~x
=	<xo) _/(*o) -
n=0
так как (см. (17) при n=0)
If с _*L = 0.
tf-iVl-x2
Аналогично доказываются следующие теоремы.
Теорема 2.2. Пусть функция /(х) е Z2 А. Тогда при выполнении условия
1'[Я^ = 0	(2.52)
tf_Wl-x2
уравнение (42) имеет единственное решение принадлежащее пространству
L1Pi, задаваемое формулой
g(x)=-1	, X е (-1,1).	(2.53)
-1 л/1 ^0 (^ Х§ )
Условие (52) следует из того, при g(x)eZ2p2 функция SP1(g(x),x0) не
содержит свободного члена.
Теорема 2.3. Пусть функция /(х) g Z2 л. Тогда уравнение (43) имеет
20
единственное решение принадлежащее пространству L1(h, задаваемое
формулой
&(*)=--, хе(-1Д).	(2.54)
(х-*о)
Теорема 2.4. Пусть функция /(х) е Z2 А. Тогда уравнение (44) имеет
единственное решение принадлежащее пространству l2Pt, задаваемое
формулой
g(x)=--	хе(-1,1).
7Г.	(Х~ХО)
Теперь рассмотрим характеристическое интегральное уравнение
первого рода с логарифмической особенностью (39). Вначале заметим, что
в силу равенств (24) и (26) получаем справедливость равенств
/“с==Т1Фо -^T„(x)<fc=- J-=L==~ln|x0 -x|T„(x)dr =
««оя-471-х2	я-Дл/1-х2 dx<>
1 f /,	x„ е[-Ц], П=1,2,„.
Л-г/1-х2(х0-х) dx0 n
~T~~ f /— - lnlxo -x|a5r= — ' [ ' ln|x0
1	4-|ni*» -xiTo(x)<* =- j—
*_J.71-x2 dx„	л J,71-x2(x„-
(2-55)
(2.56)
-х|Т0(х)&-
(2.57)
Пусть производная /(x) функции f(x) удовлетворяет условию Ha, 0<а<1.
Пусть функция g(x) является решением уравнения (40). Из равенств (56) и
(57) следует, что это уравнение можно почленно продифференцировать.
Получим равенство
- s(x)dx=f'M, ХО 6(-1,1).	(2.58)
*Л71-х2 х0-х
Из уравнения (58) функция g(x) определяется с точностью до константы (см.
теорему 1). Поэтому надо задать значение интеграла стоящего слева в (49).
Для этого умножим почленно уравнение (40) на функцию 1/*71-*о и
возьмем интеграл от обеих частей по отрезку [-1,1], получим
-x|g(x)<fe--	J 1 ln|x0 -x|dx0 -
»_JJ1-X2 ’T-iJl-X2
-In2-	= 1
Д'.171-х2	*-171-х2 ’
Из последнего равенства получаем
1 'tg(.x)dx = _ 1 1 ‘r/(x)dc
л-Д 71-х2	ln2*J^i_x2 ’
Так как /(х) е L1M, то по теореме 1 в силу формулы (48) получаем
(2.59)
(2.60)
21
g(x)=_li ijJW
я-Д x-xQ ln2 ffjVl-x2
В первом интеграле справа в (61) произведем интегрирование по частям,
получим
li~Xn
x~xQ
(2.62)
J1-X2
-------) /(x0)dr0 =
х-х0
-----------f(xQ)dxQ, хе(-Ц).
х0 (х х0)
Для оправдания равенства (62) напомним, что по определению полагаем
к ./fa)<fe0=iim( j r2;xt°. --/(x0)A0-^HZzw); (2.63)
iy\-xl (х-хпу	7/п?х₽чл 1-xJ (х-Хл)2	E
где хе(-1,1), Z/O(x,s) ==[-1,х-£]и[х + £,1].
Таким образом, из формул (61), (62) следует справедливость теоремы.
Теорема 2.5. Пусть функция /(х) е Я1 а, 0 < а < 1. Тогда уравнение (40)
имеет единственное решение, задаваемое формулой
«(*) ="-	"Л1	хе		^2-64>
^4л/1-х02(х-х0)2	In 2 * J, 71-х2
Наконец рассмотрим характеристическое гиперсингулярное
интегральное уравнение первого рода на отрезке (45). Обозначая
g*(x) = Vl-x2g(x) запишем это уравнение в виде
1	= f(x ), х0 6 (-1,1).	(2.65)
я-Д(х0-х)2
Используя замечание 1, запишем его в виде
_ 1 jCg'W)^ = /(х ), х е(_Ц),	(2.66)
^4 Хо-Х
где функцию (g*(x)) можно записать в виде (ge(x))’=A(x)M-x2. Пусть
теперь в уравнении (45) g(x) е Н1а, тогда Я(х) е На. Если теперь в уравнении
(66) функция /(х) е L2 Рг, то из теоремы 1 следует, что
(g.(x)) =a(x)/Vr^=-4=T р^-*0**0,	(2-67)
я-л/1—ж2 -1	х~хо
так как
lj^*=lj(g.(x))A=O.	(2.68)
*4 71-х2 я’-,
В силу того, что g’(-l) = 0 и с учетом равенства (68) имеем
-1	j *	. (2.69)
-I ^Л/1-Т _] т хо	Л _1	-1\1 —т (т-х0)
Если теперь во внутреннем интеграле справа в (69) сделать замену
22
ф
переменного т = cos <р, а затем tg = Л, то получим
(2.70)
,f(x6)dx„. хе (-1,1).
g’W— fin
я-д
Предыдущие рассуждения в силу формул (67) и (70) показывают
справедливость следующей теоремы.
Теорема 2.6. Пусть функция /(х) е 1г>Рг на [-1,1]. Тогда уравнение (45)
имеет единственное решение, задаваемое формулой
g(x)=T7r=^jln
Х~Х(
f(x0)dx0, хе (-1,1). (2.71)
Заметим, что в силу формулы (67) и спектральных соотношений (18) и
(37) получаем, что функция g(x) определяемая формулой (71) удовлетворяет
условию (71-x2g(x)) 71-х2 е L2tPi , а сама функция g(x) е Z2 Pz.
Замечание 2.3. Стандартный гиперсингулярный оператор НР2 (g(x)9X0)
переводит многочлен Чебышева второго рода С\(х)в многочлен -(п+1)?7л(х).
Теперь исходя из соотношений (24) и (26), а так же из формулы обращения
(64) для уравнения с логарифмической особенностью получаем
гиперсингулярный оператор который переводит многочлен Чебышева
первого рода Тп(х) в многочлен -пТл(х), т.е. верно соотношение
(2.73)
(2.74)
(2-75)
1 fl	7.(х0)&,=->»Г.(х), X е [-1,1], 11=0,1,2,...	(2.72)
^l-x2(x-x0)2
Равенство (72) еще можно получить заметив, что справедливы равенства
1~**о	= 71~'*о_______х0
71 -х2(х-х0)2 (х-х0)2 (х-х0)71-х2 ’
W = |(£/я+,«-^-,(х)), п=0,1,2,...,
хГ„(х) = 1(Г„1(х) + Гл.1(х)),п=0,1,2,...
Аналогично можно заметить, стандартный оператор с логарифмической
особенностью ZA(g(x),x0) переводит многочлен Чебышева первого рода
Тл(х)в многочлен --Ти(х), п=1,2,... . Исходя теперь из соотношения (37), а
п
так же из формулы обращения (71) для уравнения с гиперсингулярной
особенностью получаем оператор с логарифмической особенностью который
переводит многочлен Чебышева второго рода Un(x) в многочлен - —.ид*),
п + \
т.е. верно соотношение
Jin-------I ги„(х,)ф, = -™1/л(х),хе[-1,1],п=0,1,2,... (2.76)
fl-71 ~х2 -1 1 —хх0 + л/1-х 2 71-х2	л + 1
23
П. 3. Обобщенные функции на Гильбертовых пространствах и особые
интегральные уравнения.
В предыдущих двух пунктах были рассмотрены особые интегральные
уравнения в наиболее широких, используемых до настоящего времени,
пространствах. Пространствах £2 периодических функций с интегрируемым
квадратом модуля и пространствах Llp функций с интегрируемым с весом
р(х) квадратом модуля. Однако в последнее время некоторые прикладные
задачи привели к необходимости рассмотреть особые интегральные
уравнения в классах обобщенных функций на отрезке.
В качестве примера рассмотрим задачу обтекания профиля с отсосом
внешнего потока. Пусть кривая L является контуром непроницаемого
профиля крыла, безотрывно обтекаемого потенциальным потоком
несжимаемой жидкости со скоростью Uo(M), т.е.
ио(М)=\/Ф^(М)=(дФй(М)/дх) 1+(дФ0(М)/ду) j, /, j-орты осей Ох, Оу, М-точка
полскости Оху, а функция Фо (Л/) гармонична на всей плоскости.
Непроницаемость профиля и безотрывность его обтекания в данной точке
означают, что частица жидкости в этой точке имеет только касательную к
профилю скорость. Предполагаем, что кривая L задана параметрически
х = %(/), у = y(t)9 fe[0,/], с направлением обхода по часовой стрелке, если
кривая L замкнутая, а параметр t является длиной дуги кривой L. Функции
х ’(О,У (О удовлетворяют условию Гельдера степени а на [0,7], т.е.
х ’(0,У (0 ^На на [0,/], и так как параметр t является длиной дуги на Z, то
х’2(0 + у’2(0 = 1, / е[0,7]. Если кривая L замкнутая, то функции х = х(г),у = у(г) и
их производные необходимых порядков периодические с периодом t, причем
для простоты будем полагать / = 2я. Предположим еще, что в точке Mq(xq,yq)
профиля L (xq=x(tq),yq=y(tq), tqe(0,l)) происходит всасывание внешнего
потока внутрь оболочки профиля. Требуется найти поле скоростей в
жидкости, возмущенное профилем.
Будем моделировать профиль вихревым слоем [1] интенсивности
у(/) = в точке M(t) = Л7(х(г),у(О) кривой L. Тогда профиль возмущает в
любой точке M0=M(xQ9y0) плоскости скорость Kz(Af0), определяемую
формулой
К,(М.)=2-	, M„eL,	(3.1)
Гмм.
xl(M,M0)=x0yx(M,M0) = y(I-y(t),
rM	|=| i Xi (М,мо) + j у, (M, Ma
В каждой точке M кривой L определен орт касательной
тм~х\t) i + y (f) j. Выберем направление орта нормали к кривой L в этой
24
точке по формуле пм=-у (t)i+x (О j и будем обозначать ту сторону кривой
Z, куда направлен вектор пм, знаком “+” или (Z+), а противоположную
_»±
знаком или (Z‘). Тогда известны [1] формулы для скорости ГХ(МО) от
вихревого слоя в точках MQ кривой L при подходе к ним с соответствующей
стороны
^(Л/о)=— f Z±^r(t0),Ma&L. (3.2)
I	2
Устройство всасывания внешнего потока внутрь оболочки профиля
будем моделировать стоком интенсивности Q в точке Mq(xq9yq} согласно
экспериментам (см. [314] в [1] ). Поле скоростей И7(Л/0) , порожденное
стоком, имеет особенность в точке Mq(xq9yq) и задается формулой
(3.3)
Всасывание потока внутрь оболочки профиля означает, что полное
поле скоростей U(M0)= UQ(M0)+Vy(M0)+Vq(M0), возникающее в этой задаче,
имеет особенность типа (3) в окрестности точки Mq(xq9yq) с внешней
стороны профиля (z+) и не имеет особенности (является гладким ) с
противоположной стороны профиля (z~). Среди трех слагаемых в U(M0) два
известны и неизвестным является только Vr(M0). Для его нахождения
достаточно найти /(0 =	, которое находят из условия непроницаемости
профиля и безотрывности обтекания. Если это условие выполнять на (Z+), то
придется пропустить точку Mq(xq,yq)9 так как в ней не выполняется
непроницаемость профиля. Поэтому будем выполнять это условие на (L ).
Получим равенство
с7-(Л/о)и^0=0, MogZ,
или
¥-г(М0)пи^-и-0(М0)пи,-У;(М0) пи„ Мо<=1.	(3.4)
Равенство (4) можно записать в виде
f):V'(<’'°>+	t„ e[O,Z],	(3.5)
где	-и;(М0) nUa-V~(M0) Hn„	х^М,М0),	ух(М,Мй),
Уравнение (5) справедливо для замкнутой и разомкнутой кривой L. Если
кривая L разомкнутая, то уравнение (5) может быть записано в виде
25
- у-	ч e (о./),	(3.6)
2л-ого-^ i
а если контур L замкнутый, то в виде
_ J_ Lgt2z(0«*+2k2(r0,0Z(r)<*=/(г0), г0 е[0,2я).	(3.7)
™ о 2	о
Если х” (t)t у" (t) е На на [0,1] ([0,2?г] для замкнутой кривой), то можно показать
[1], что в уравнениях (6), (7) ядра K^t) и Х2(/0,г) так же принадлежат
классу На на соответствующих множествах.
Рассмотрим более подробно правую часть в уравнении (5). Здесь
функция /(г0) состоит из двух слагаемых. Функция U~(MQ)nMo является
«хорошей» в окрестности точки Mq(xq,yq), так как она лежит на участке
гладкости кривой L. Изучим теперь функцию ^‘(A/o) лм0. Из физики
движения жидкости ясно, что эта функция описывает распределенную
плотность количества жидкости, протекающей через точки кривой L.
Отсюда видно, что для точек М0 g L и Л/о * Mq имеем У~(МЪ) пМо^Vq(Mb) пМ()
= (2xY'Q<oAM^Mq) , ^{Мй,Мч)^г~2^ ^МаМя,пм0). В традиционных курсах по
математической физики доказано, что если x\t\y"(t)^Ha на [0,1] , то
функция а>х(Мй,Мч) еНа равномерно по обеим переменным при
соответствующем доопределении ее в точке MQ= Mq . Если кривая L имеет в
точке Mq касательную, то из сущности источника следует, что через точку
Mq со стороны L~ протекает QI2 жидкости (направление движения точек
жидкости в точку Mq со стороны Г имеет острый угол с вектором ~пМч ).
Поэтому для точек кривой L можно написать V~(M0) пМо = (2n)~xQa)}(MQ,Mq)
+ (0/2)£(Л/о-Л/9), где MQ - произвольная точка кривой L и функция
#(Mt-Mq) определяется равенствами: #(Mo-Mq)=O при MoeL, MQ*Mq, и
Л(Л/0 -AQ^+оопри Мо= Mq el, а также
р(м0-м,)<й0=1
Теперь ясно, что поскольку в правой части уравнения (5) стоит обобщенная
функция, то решение тоже надо искать в классе обобщенных функций. Как
понимать эти обобщенные функции и как понимать особый интеграл от
таких функций покажем ниже.
Один вариант обобщенных функций на Гильбертовых
пространствах. Рассмотрим бесконечномерное вещественное Гильбертово
пространство Н (для комплексного пространства ниже следующие
построения остаются справедливыми с соответствующими изменениями).
26
Операцию скалярного произведения векторов Н запишем как (/,g).
Пусть система векторов	является ортонормированным базисом
в пространстве Н, т.е. удовлетворяет условиям
, 8}J = 1, при i=j и 8t J = 0, при i'Ф j,	(3.8)
тогда любой вектор f из Я представляется в виде
/=Ха^*’	<3-9>
*=1
где для коэффициентов ак выполняется соотношение
СО
2Х<«>-
*=1
При этом норма вектора f определяется по формуле
И=ЖТ)=^БГ	(з-ю)
Определим теперь в пространстве И основное множество S векторов.
Определение 3.1. Вектор f пространства Н принадлежит основному
множеству S, если его разложение по базису имеет только конечное число
коэффициентов отличных от нуля, т. е. для разложения (9) существует такое
натуральное число К, что
ak=0,k>K.	(3.11)
Отметим, что множество S=S(H) является линейным и всюду плотным в
пространстве Н, но не является замкнутым в Н. Действительно, любой вектор
f из Н, имеющий в представлении (9) бесконечно много отличных от нуля
коэффициентов, является пределом в метрике пространства Н
последовательности векторов
(3-12)
ы
Поэтому введем во множестве S другое понятие сходимости векторов.
Определение 3.2. Будем говорить, что последовательность векторов
/(1),/(2),...,/(и),... из множества S сходится к вектору f конечномерно, если
существует такой общий номер К, что	
а{”} =0,к>К,п = 1,2>... ,	(3.13)
и для любого п=1,2,... выполняется условие Нт|дг£и> - а. I = 0, к=1,2,...	(3.14)
Кратко это будем записывать так	(3.15)
Итак, другими словами, последовательность векторов /(1),/(2),...,/(я),... из S
сходится к вектору f из Н, если все они лежат в общем конечномерном
подпространстве из Н и выполняется соотношение (14). Таким образом,
27
вектор f лежит в том же конечномерном подпространстве и, следовательно,
принадлежит множеству S. Это означает, что конечномерная сходимость в S
обеспечивает замкнутость S. Теперь множество S с введенной
конечномерной сходимостью будем называть пространством S основных
векторов. Можно еще сказать, что в S введена конечномерная
покоординатная сходимость.
Введем теперь понятие обобщенной функции.
Определение 3.3. Обобщенной функцией называется любой линейный
непрерывный функционал F на пространстве основных векторов S. Значение
функционала F на основном векторе f будем записывать F(f) или (F,f). При
этом непрерывность функционала понимаем следующим образом.
Функционал F называется непрерывным на S, если из конечномерной
сходимости последовательности /(,),/(2),...,/(я),... векторов пространства S к
вектору f из того же пространства следует, что выполняется соотношение
\im,(F,fw) = (F,f) .	(3.16)
Множество всех обобщенных функций обозначим через 5’ = S' (Я).
Множество S' является линейным, если линейную комбинацию AF+yG
обобщенных функций F и G определить как функционал, действующий по
формуле
(AF + /X7,/) = 2(F,/) + XG,/), feS.	(3.17)
Следуя методике в [8], можно показать, что функционал AF + pG линейный и
непрерывный на S, т.е. принадлежит S'.
Интересно отметить следующее. Покоординатная сходимость в S такова,
что справедлива следующая теорема.
Теорема 3.1. Любой линейный функционал на S является непрерывным на
S в смысле определения 3.
Доказательство. Действительно, пусть последовательность
/(1)векторов пространства S конечномерно сходится к вектору f
из S. Это означает, что для коэффициентов разложения {«"'’} векторов /(я)
существует некоторое такое К, что а(кп) =0,k >Л> = 1,2,.... Поэтому, если F
линейный на S функционал, то имеем
= lim£ )а‘”> =	= (F.f erf'.) =	.
Таким образом, теорема 1 справедлива.
Отметим, что любой ряд G = '£bky/k) а, следовательно, и любой элемент
пространства Н, определяет на S линейный функционал по следующему
к
правилу. Пусть f = £ fky/k е S, тогда по определению полагаем
*=!
28
G<J) = (G,f) = ^bkft.	(3.18)
*=1
Функционал G действительно является линейным, так как
G(V+A?) = ZMV, +wk) = &>tft	=^(Л+ /<?(«)•
4=1	4=1	4=1
Из теоремы 1 следует, что функционал G является непрерывным на S, а,
следовательно, является обобщенной функцией на Н.
Отметим теперь, что значение обобщенной функции F на элементе f из S
определяется как
F(f) = (F,f) = ^ft(F^t).	(3.19)
4=1
Покажем теперь справедливость следующей теоремы.
Теорема 3.2. Пусть задана некоторая обобщенная функция F на Н. Тогда
ряд
G =	(3.20)
4=1
является обобщенной функцией из S', совпадающей с F.
Доказательство. Пусть f произвольный элемент из S. Тогда по
определению имеем (см. (18) и (19))
с(л=Е(^^)/*=^(/)-
4-1
Итак, получаем важный вывод, что множество произвольных рядов вида
^Ькук и множество обобщенных функций на Н, при нашем построении,
*=i
совпадают.
Определим теперь сходимость в S как слабую сходимость
последовательности функционалов.
Определение 3.4. Последовательность обобщенных функций F15F2,...,Fn,...
из S сходится к обобщенной функции F из Sесли для любого вектора f из S
имеем (Fn,f) (F,f\n ->оо. В этом случае мы будем писать	в
S.
Линейное множество S с введенной в нем сходимостью назовем
пространством обобщенных функций S .
Теперь можно доказать теорему.
Теорема 3.3. Пусть последовательность F1SF2,...,F„,... из S такова, что для
каждого вектора f из S числовая последовательность (F„,/) сходится при
п -* оо. Тогда функционал F на S, определенный равенством
(F,f) = lim(F„f), feS,	(3.21)
29
также является линейным и непрерывным на S, т. е. FeS .
Доказательство. В силу теоремы 1 достаточно доказать линейность
функционала F. Действительно, имеем
(F,af + ffg) = \im(Fn,af + Pg) = lim(a(F„ /) + P{Fn, g)) =
П-ЮО	n-^X}
alim(F„/) + ^lim(F„g) = a(F,f) + 0(F,g)
п-юа	и->оо
где а и p произвольные числа, a f и g элементы пространства S.
Теорема 3 показывает, что пространство обобщенных функций S является
полным относительно введенной в нем сходимости.
Обратимся теперь к рассмотрению характеристических особых
интегральных уравнений в пространствах обобщенных функций. Причем
начнем с особых интегральных уравнений на классе периодических
функций. Для этого вначале напомним спектральные соотношения,
полученные в п.1
— | In Isin ——-le^dp -
я о I I
где sign(w) =
1,п> О,
0,и = 0,
sign(«)ew°, и = ±1,±2,...
п
-21п2,и = 0,
(3-22)
_L [ctg~—~e'n,pd(p = -i • sign(n)e,n<Pa ,n = 0,±l,±2,...
2zr i 2
(3.23)
_L_ f—e —= ~n sign( rc)ew°, n - 0, ± 1, ±2,...	(3.24)
4*o\in
Рассмотрим теперь Гильбертово пространство L2 комплекснозначных
2яг периодических функций с интегрируемым квадратом модуля,
построенное в п.1. В этом пространстве ортонормированным базисом

является система функций у/„ =
и = 0,±1,±2,....
Ясно, что для этих
функций также справедливы спектральные соотношения (22)-(24). Пусть
теперь s’(L2) пространство обобщенных функций построенное на
пространстве L2. Теперь для любой обобщенной функции g(^) = ^g(«)^M из
S'(L2)будем полагать
1 I	I	-|	д
Jlnsin—\g(p)d<p—~	-sign(n)g(n)y/„(<i>li)-2\n2g(O)y/<)(.g>0'),Q^)
о I 2 I sX ”
Z0=+l,±2,...
30
S(g(<p),<p0)-~ [ctg^~g(<p)d<p=-Y‘ ’«”(”)g(">.(₽0).	(3-26)
2^ о	nez
H(g(<P),<P<>)=~ f-±—~tg(p)d9>=-Yn'siS"^^n')4/^9>o)-	(3-27)
2
Из соотношений (25)-(27) можно сделать следующие выводы для уравнений
Дг(Лл) = - ftalsin^—= /(%),% е[0,2л],	(3.28)
п 5 I 2 I
S(g(₽),9>а) = -±- (crg^~g(<f>)dp = /(<?„),?<, е [0,2?т],	(3.29)
хЛ' 0	2
g(g(y),g>o)-j-J	= /(%)> W> е [0,2^],	(3.30)
4jr»sin2^
где f((p) - 2л периодическая обобщенная функция.
Уравнение (28) имеет единственное решение для любой обобщенной
функции /(<?), а уравнения (29) и (30) имеют решение с точностью до
константы и условием их разрешимости является равенство
J/W^ = 0,	(3.31)
о
где под значением интеграла от f(q>) по [0,2 л-] понимается величина
72л /(0). Причем формулы обращения для этих уравнений в пространстве
обобщенных функций те же, что и в пространстве L2 и имеют следующий
вид (см. п.1). Для уравнения (28) имеем
g(<р) = ~~ j	,	(3.32)
°sm2-^0	211,2 о
2
для уравнения (29) имеем
gW = -у-	,	(3.33)
2л ;	2	2л
для уравнения (30) имеем
g(<P)- ~ jlnlsin^——If (&№<?,, +^-,	(3.34)
я-0 I 2 Г 2л
где в (33) и (34) для С выполняется равенство
jgGW-C,	(3.35)
о
так как для разрешимости этих уравнений должно выполняться равенство
(31).
Однако для построения вычислительных методов решения уравнений
(28)-(30) надо уметь оценивать близость решений этих уравнений по
близости их правых частей. Для этого заметим, что если взять произвольный
ряд из экспонент (обобщенную функцию) и применить к нему один из
31
операторов L, S или H, то, получим аналогичный ряд у которого
коэффициенты отличаются от коэффициентов исходного ряда множителем
порядка	ия=тах{1, |и|}. Поэтому естественно возникает идея использования
пространств ЯАтипа пространств Соболева обобщенных функций. Пусть
2п периодическая обобщенная функция на £2, тогда
/(?>) =	/W = </(P).^.)=4=	(3.36)
neZ	о
где под (/(<?), y/J понимается значение обобщенной функции (линейного
функционала) f(<p) на функции (j/n=~^==e ing>. Обозначим через Нл
yl2n
множество таких обобщенных функций что
Ил=(ЕИя")|1 2)*<«>>	(з.з7)
«eZ I I
т.е. функция /Чр) = £Х/(я>л(^) принадлежит пространству L2.
neZ
Введем теперь для функций из Ял скалярное произведение по формуле
=	f(n)g(n).	(3.38)
neZ
Так введенное скалярное произведение делает Нл гильбертовым
пространством причем Я0^ - пространство функций с интегрируемым
квадратом модуля на [0,2я-].
Теперь рассмотрим характеристические особые интеральные
уравнения на отрезке. Для этого опять вначале напомним спектральные
соотношения, полученные в п.1
1 J- L=lnjx0 -х|Т„(х)А = -Я„(х0),х„ е (-1,1),« = 1,2,...,	(3.39)
=	*о 6(41), « = 0,1,2,...,	(3.41)
л-4-Vl-x2 х„-х
1	=т„(х„), х„ е(-1.1), « = 1,2,...,	(3.42)
XQ-X
1 jlE*	M , x e (_U), „ = 0,1,2,...,	(3.43)
1 jE* = x0 e(-l,l), « = 0,1,2,...,	(3.44)
-	= -(« + !)(/„(x0), x„ e(-l,l), « = 0,1,2,...,	(3.45)
Я’-!	(x„-x)2
где	T„(x) = cos(narccosx),	J7„(x) = sin((« + l)arccosx)/sin(arccosx),
QAx) = U„(x)-U_,(x). Теперь характеристические
32
особые интегральные уравнения на отрезке записываются в виде (см. п.1)
(g(x),x0) = - J“r="lnlxo -x|g(x)dv = /(х0),	е (-U),	(3.46)
s*(g(x),xj = -	= /(х0), х0 е (-1,1),	(3.47)
А	Я’ЛТГ-Х2 х0-х
SP1(g(x),xt,) = - Ь/1-х2	= /(х0), х0 е(-1,1),	(3.48)
п _*	х0 - X
SA(g(x),x0) = l jJp^^x^ = /(x0), хое(-1,1),	(3.49)
п Д V1 + х х0 - х
S,.(g(x),x0) = 1 Ш*= /(x0), x0 e (-1,1),	(3.50)
л _•} V1 - x x0 - x
#ft(g(x),x0) = - R/l-x2 g(*^; =/(x0),x<) e(-1,1),	(3.51)
*.i	(x0-x)2
где pj (x) = -г- - - , p2 (x) = P)'1 (x), p3 (x) =	, p4 (x) = p(’ (x).
Vl-X2	Vl + x
Рассмотрим Гильбертовы пространства построенные в п.2 и возьмем
пространства обобщенных функций S'(L2p). Пусть теперь /(х) является
обобщенной функцией на L2p, т.е. /(х) е S'p = S (L2 p), тогда
/W = 2}/(«K.,(x), /(«) = (/(х),^„.р(х)) = Jp(x)/(x>^(x)dc, (3.52)
п=0	-1
где под (f(x),y/„p(x)) понимается значение обобщенной функции (линейного
функционала) f(x) на функции у/лр(х). Обозначим через множество
таких обобщенных функций f(x) на p , что
Ил., =(2Ур(”)|2)^<00’	(3-53)
neZ I I
т.е. функция /*(х) = ^вл f(n)vn(>(x) принадлежит пространству L2p. Введем
»€Z
теперь для функций из Нр скалярное произведение по формуле
(/,«) =	f{n)g{n).	(3.54)
яе2
Так введенное скалярное произведение делает Нр гильбертовым
пространством причем HQp = L2p - пространство функций с интегрируемым
квадратом модуля на [-1,1]. Теперь для любой обобщенной функции
=	=	k=l-4, будем полагать, что значение
SPk{g{x),x^, ^ = 1,2,3,4, Hj£(x\x^ получается по правилу
перемены местами знаков интеграла и суммы в этих операторах. Теперь
соотношения (39)-(45) показывают следующее. Оператор L^g^x)^
взаимнооднозначно отображает пространство обобщенных функций S (12 Р1)
33
на себя и пространство Н\ на пространство Я^+1; оператор SA(g(x),x0)
отображает S'(L2Pi) на S’(I2pi) и ядром отображения является множество
констант и с таким же свойством на ЯД; оператор 5^(g(x),x0)
отображает 8 (L2pi) в S\L2pi) и для функций /(х) из образа выполняется
соотношение
J-^La = 0	(3.55)
и с таким же свойством ЯА в ЯА; оператор £A(g(x),x0) отображает
взаимнооднозначно S’(L2p3) на S'(L2At) и Я* на Я*; оператор SP4(g(x),x0)
отображает взаимнооднозначно 8(L2pa) на S'(L2pi) и НРа на Я*; и , наконец,
оператор Я^(х),х0) отображает взаимнооднозначно S'(L2p2) на S (L2pj) и
Я* на Я^1.
Важно, что теперь можно говорить о близости правых частей и
соответсвующих решений, являющихся обобщенными функциями.
Теперь сделаем следующее важное замечание. Для уравнений (46)-(51)
на базисных элементах соответсвующих пространств справедливы формулы
обращения, полученные в п.1. Поэтому, в силу определения значений
операторов IA(g(x),x0), SA(g(x),x0), к = 1,2,3,4, #A(g(x),x0) на обобщенных
функциях соответствующих пространств, для уравнений (46)-(51)
справедливы следующие формулы обращения.
Для уравнения (46) имеем
gW = -()-, --{-°-----/(x0)A0--L	хе(-1Д),	(3.56)
* ^^(x-xrf	1п2
для уравнения (47) имеем
g(x) = _1	+ £, х е (-1,1),	(3.57)
Я Д	X - х0	я
'[	=	(3.58)
для уравнения (48) имеем
g(x)=-- f	Хб(_1Д)>	(3.59)
Л-I^I-XqX-Xq
при выполнении условия (55),
для уравнения (49) имеем
g(x) = -- j |p^-/(x0)dx0, х е (-1,1),	(3.60)
для уравнения (50) имеем
gW =	, хе (-1,1),	(3.61)
7Г Др + Х0
для уравнения (51) имеем
34
g« =
Pn
f(xo)dxtj, xe(-l,l).
(3.62)
(3.63)
(3-65)
x-x(
1-XXO +71-J
Используя проведенные выше построения в пространствах обобщенных
функций получим теперь некоторые точные решения характеристических
особых интегральных уравнений первого рода в этих пространствах.
Например, задача бесциркуляционного обтекания кругового цилиндра
потоком идеальной несжимаемой жидкости при наличии на нем устройства
отсоса внешнего потока сводится к решению уравнения [3]
2- jc#^+^g(«>)rf(?> = <?(₽0
хтг *	2	2л
В равенстве (63) под функцией 8{(p-q) понимаем 2л- периодическую дельта
функцию, определяемую равенствами: S(<p-q)=O, при (p*q, и
при (р = q,(p,q е[Ъ,2л\, И ДЛЯ любой фуНКЦИИ у/п((р), п = 0,±1,±2,...,
2п
выполняется равенство ft(<p-q)wMdp = y/n(q). Из этого вытекает, что если
о
w((p} обобщенная функция, то так же имеем равенство - q)w((p)d(p = w(^).
о
Теперь из (33) получаем, что решением уравнения (63) является фунция
g(<p) = -~ctg~r-+£-	(3.64)
2л 2 In
Таким образом, получаем формулу[3]
-у- fag^—r^^-crg^—d? = #(9>0е[0,2^].
2л ’	2 2л 2	2л
Этот же результат получается [3] если обобщенные функции $(<p-q)n
~—ctg¥^L представить рядами Фурье по функциям ^/л =-71=е'и<₽и
2л 2	^2л
подставить эти ряды в уравнение (63). Действительно, имеем
- 2- ctg = / W=Е 7 Wv'» W >
2Я 2-	„eZ
=	ff(9>)e-m”df>=-2- fctg^ ^-e^dp.
у2л о	2л 0	2 -\/2т~
Теперь из (1.4) следует
- ~ c‘g=Е 'W( «М	(р)» g е [0>2я-] •
Zfl- 2 neZ
Используя опять (1.4) получим
2- jczg^—^(”2-rtg--^)^= 2_ jctg^±Yisign(n)^„(q)i//„^)d9>=
2я '	2	2я-	2 2п J 2
^isign(n)^„(d)/- \ctg^~^4/H(<P)d<p =
»€Z	о	2
ие7,и#0
С другой стороны, по определению <?)функци полусаем
(3.66)
(3.67)
35
(3.68)
3(<р - «) - 2 ’*'«	(^o ) - Z V'» (^ (Po )+г •
neZ	neZ ,«*O
Сравнивая (67) и (68), видим справедливость равенства (65).
Из представлений рядами Фурье (66) и (68) следует, что обе функции
и принадлежат любому пространству /Гдля Л<-~.
2п 2	2
Заметим так же следующее. Из соотношений (22), (23) следует, что равенство
(65) можно почленно интегрировать и поэтому верно соотношение
-- flnisin^^|^c/g^^^ = -^ + F(^0-^)) + C ^о^е[О,2я-], (3.69)
Я1 j I 2	11П 2	2.ТГ
где - #)) = 0, §<(p<q\ F(8((p-q)} = \^ <p-q\ F(8(<p-q)) =	q<(p<2n^ И
0 < $ < 2яг, а для q=0 имеем F(^)) = 0, ^> = 0; F(^)) = |, 0<p<2?r; F(£(p)) = l,
(p-2n. Так как свободный член ряда Фурье для функции справа в (69)
должен быть равен нулю, то С = — (q - я). Отметим, что функция F(8(<p - qj)
2я
является первообразной для функции 8(<p-q), обращающаяся в нуль при
9> = 0.
Решение задачи о расчете входного сопротивления тонкой
проволочной антенны при запитке антенны источником тока, повлекло за
собой рассмотрение гиперсингулярного интегрального уравнения на отрезке,
в правой части которого стоит функция с особенностью типа 1/х внутри
области поиска решения [5]. Аналогом этого уравнения в периодическом
случае будет уравнение
у- J-----, <p„,q е[0,2я].
4” °sin2y°-p	2
(3.70)
2
В силу формулы обращения (34) для рассматриваемого уравнения и формулы
(69) получаем, что решнием уравнения (70) будет функция
gW =	F(8(# ~ «)) + С, + £, <р 6 [0,2 я],
2я	2я
для которой выполняется равенство (35).
Для приближенного метода решения уравнения (63) методом
дискретных вихрей важно, что функция 8((p-q) является пределом в смысле
определения 4 последовательности 2я- периодических функций 8h((p-q)-\,
h
8h(<p-q)=0,	+	$б[0,2л), т.е.
^m8h{q>-q) = s(<p-q).
Теперь рассмотрим уравнение (47) на отрезке. Возьмем в его правой
части функцию /(х0) = 8(х0 -q), qe (-1,1), тогда в силу формулы (57) получим,
что решением этого уравнения будет функция (С=0)
(3.71)
36
, ч 1 V1-?2
g(x) = —.
я x-q
Подставляя теперь последнюю функцию в уравнение (47)
интересное равенство
(3.72)
получаем
J г—-у	(	)dx - 8(х0 q) 9 х0 е ( 1,1).
471-х2(х0-х) * Х'Я
Наконец рассмотрим еще частный случай уравнения (51)
1 L/Г-? g(x)dx - 1 Y in
/г	(х0 х) х0 q
который встречается в теории антенн. Оператор, стоящий слева
отображает пространство S (L2pi) на себя. Поэтому решение этого уравнения
(3.73)
(3-74)
в (74)
и его правую часть рассматриваем как элементы этого пространства. Как и
при обращении уравненя (2.45) сведем уравнение (74) к кравнению
_1	=	x0,ge(-l,D,	(3.75)
я-Д х0-х x0-q
которое эквивалентно уравнению (74) при выполнении условия
1 f(g*(x))A=0,	(3.76)
Я -1
где g’(x) = 71-jc?g(x). Из формулы обращения (2.48) для уравнения (2.41)
получаем, что общее решение уравнения (75) дается формулой
(g‘ W)’ = -я^(х - q) +	, х,q е (-1,1).	(3.77)
л/1~х
В силу равенства (76) получаем, С=1. Поэтому, так как g’(-l) = 0, имеем
g’(x) = Vr-?g(x) = arcsinx + ^-^F(^(x-^)), х,#е(-1,1),	(3.78)
где F(^(x-^))=0, при 0<х<<7, F(S(x-q))=^ при x=q, F(£(x-#))=1, при
q<x<\. Функция F(8(x-q)} является первообразной для функции 8(x-q) на
отрезке [-1,1], qe(-1,1). Таким образом, решение уравнения (74) дается
формулой
g(x) = jJ-y (arcsinх + ~ - nF(8(x - q))), x,qe (-1,1).	(3.79)
VI - x	2
П.4. Метод дискретных вихрей численного решения
характеристических особых интегральных уравнений.
Вначале рассмотрим численное решение характеристических
особых интегральных уравнений на отрезке. Потом будем рассматривать
эти уравнения в периодическом случае. Поскольку в периодическом случае
при рассмотрении метода дискретных вихрей удобно будет ввести в
рассмотрение точки комплексной плоскости, которые обычно принято
обознать буквой t [2], то и при рассмотрении метода дискретных вихрей для
особых интегральных уравнений на отрезке будем использовать вместо
переменной х переменную t. В классической теории сингулярных
37
интегральных уравнений уравнения (3.47)-(3.50) записывают единным
образом в виде
fy— =/('.Мое (41).
-1 *0	*
Теперь используя формулы обращения для уравнений
соотвтствующие решения уравнения (1) можно записать в виде
(4.1)
(3.47)43.49),
ф(/) = -Ла(') К'('о)
*	L-1
где Vi = l, vo=v,1=O,
^(O = p3(O = J^, Я,(') = «(') = -7=^=7, Я.,(г) = Я1-1(/) = р2(0,
V*+l	yjl-t2
и будем называть их соответственно решениями индекса к -1
неограниченными на обоих концах отрезка, индекса к = -1 ограниченными на
обоих концах отрезка и индекса к - 0 ограниченными в точке 1. Решение
индекса к = 0 с особенностью вида ^(t+l)l(\~t) рассматривать не будем, так
как из приведенных рассуждений будет видно, что надо изменить в этом
случае.
Пусть множества Е = {ч,к = 1,и Ео =^,у = ОД,...,л} образуют
каноническое разбиение отрезка [-1,1] с шагом й, т.е. tk = -l + kh, tQj=tj+~,
2
k,j= 0,, h =------.Справедлива следующая
п + 1
Теорема 4.1. Пусть функция f(t) принадлежит классу На на [-1,1].
Тогда между решениями систем линейных алгебраических уравнений
(43)
-vKCx ,fe(-l,l),
'о
(4.2)
Г07
л,
n-1,
(4-4)
к~1 lOj~lk
/о» + £	j = 0,1,.п,	(4.5)
*-i Ч
и решением ф(/) индекса к = 0,1,-1 соответственно в(2) уравнения
(1)выполняется неравенство
* = 1 ",	(4.6)
в котором величина 6n(tk) удовлетворяет неравенствам:
1) для всех точек tk е [-1 + 5,1 - 5], где 8 > 0 сколь угодно мало,
Д>0;	(4.7)
2) для всех точек tk е [-1,1]
Л>0.	(4.8)
i=1
38
Доказательство. Заметим вначале, что в силу результатов о
квадратурных формулах для сингулярного интеграла на отрезке [2] системы
(3)-(5) аппроксимируют уравнение (1). Для определителей систем (3)-(5)
получим
Р^=ЛИД(ГИ), к = 0,1,-1,
	1	1		1 1
	Ан А	Аи ~		Ai~A	Ан _ А
4,” =	1 Ан» “А	1 А)я ” Ai	, А’1”’ =	1 1 Ал-1 “А	А)я-1 ~^п 1 ... 1
1
А» — А
1
Чп
1
А)0 ~
1
Несложно показать, что
nnk-и nnU-'J
д(и) _ <Р^п	1<(-*.')<т<р<л-<(*-)
~	«<(*) п ,	.
П nta-'J
m-=l-<(-r) р=1
(4.9)
(4.Ю)
где £(х) = 1 при х>0 и £(х) = 0 при х<0. Покажем это для Д(я), так как для
остальных будет аналогично, с некоторыми очевидными изменениями.
Вычтем последнюю строчку в определителе д(оя) из всех предыдущих и
из каждого столбца вынесем множитель 1/(Ал-^)? Л = 1,...,л, а из каждой
строчки множитель	т = 1,..., п -1. Получим
Теперь вычтем в последнем определителе последний столбец из всех
предыдущих, вынесем множители из каждой строчки и столбца и разложим
полученный после этого определитель по последней строчке. Получим
Ц(А)л “А)жХАя ”А«)
д(;) =	---—— Д(я-1).
Метод математической индукции заканчивает доказательство формулы
С-
Из формулы (10) видно, что при любом п. Применяя правило
Крамера решения систем линейных алгебраических уравнений, получим
39
(4.Н)
ф <t )=^LI	)+v tr Cc
ФДk)	£1?	h ф C’
где vK определена в (2), a <(x) в (10); D^k и Л(^ получаются соответственно
из D™ и д(”} заменой к -го столбца на столбец из свободных членов правой
части системы, a получается из вычеркиванием j -й строчки и к -го
столбца. Для Д(^ так же как и для Л1^, получаем
nrik-'J ППк-d
ISm <pSn	1-^(г)<т</7^л-^(лг)
д(и) _ я, р*к________m,p*j
~	п-^к) п .
(4.12)
т=1-^( к-) р--1
m*j р*к
Поэтому формула (12) получит вид
Ч>М = ~^ Т	к = 1,
(4-13)

*k
где
= П (*<>« ~ **) М ” ’
Ж=1	/
/ т*к
^0,0; =П^ ~^)/r“AlmX
л>=1	/ т=1
/ т*к
r(«) _ rW 1	г(«) _ т(я) L * )
~*0,к	9 11,0; “ -*0.07 Von l0j h
lQn lk
/(«) —T^lt -t 1 7(л) _ 7(»)	1
l~\k ~ 10,k \lk *00/J •*4,0; ~20,0;	_	*
*0;	f00
Далее, имеем
|45=|л-.^,	(4.14)
p*->=nf1+TTri’ ^ = n[1+7frl> po = 1’
m=lV	*k J	m=k+\\ ‘ni lk J
(4-15)
^-1=Й| 1+г*нН = flf1+~г |> po.»=L
m=4	*0;	^0т J	w=;+iy ^0j	^0т J
Напомним, что в рассматриваемом случае
fk =-^ + khi tok~tk -
к	ок к 2	п + 1
^-h=tQm-hk=Km-k).	(4.16)
Следовательно, формулы (14)-(16) дают
p-.=n4-V”-‘=nfi+" •	<4-17>
т) til т)
В [262 из [1]]имеется следующая формула из теории гамма-функции:
40
(l + /3)(2 + /?)-(n + ffl . rS , С(
(4.18)
которую можно записать так:
Так как нас будет интересовать формула (18) с точностью до величин
порядка л, то можно написать ГТ f1 + £ 1	+ О((„ +1/-‘). Ш т) Г(1 + /7) Полагая теперь р = ± 1/2, получим Л 1 =-Ц- + о(*‘,/2)> р. Ь = <'n~k^V2 +o((n-k + V)111). Г(1/2)	’ "	Г(3/2)	V	7 ’ Напомним известные формулы Г(1/2) = -Пг, Г(а +1) = Г(а)а, Г(3/2) = -Д/2. Таким образом,	(4.18-) (4.19)
Л-1Л-»	+ °(k''l2(r! ~к+i)‘V2)+ °(*3/2 («- *+1)V2) • Аналогично имеем	(4.20)
_ _	2 1 J+1/2 J?. 1YVY . lY*2') р	Р	— 	 1 J	<	1 fJ	» J		In 	 » 1 _ I	I	
+	2J I	2) J S( iY/2/	1Y3/2>1 (y 2J <	J	2) J Отметим, что в силу формулы (16) имеем = ylh(n~k + V) ,	-а =л/Й,	(4.21)
7Щ; = 7Л(П-)+1/2), ^07-a=^(j+l/2).	(4.22) Из формул (22) следует, что формулы (20)и (21) можно записать в виде Л -d	 к	
(4.24)
? -nf, *
0’<м) J ч^-q*2 у
Аналогичные рассуждения показывают, что для к -1-1 имеем
(4-25)
41
(4.26)
Подставляя теперь формулы (25), (26) в формулу (13) и пользуясь
рассуждениями, аналогичными тем, которые были проведены при
исследовании квадратурных формул метода дискретных вихрей для
сингулярного интеграла на отрезке [2],получим
(4-27)
* а	VZ0
где величина |0г(/4)| удовлетворяет неравенствам ( 7) и (8). Теорема 1
доказана.
Замечание 4.1. Поясним смысл неизвестной /Ол в системе
. Эта система без /Од переопределена (число уравнений больше числа
неизвестных) и, как правило, несовместна. С другой стороны, эта система
должна аппроксимировать уравнение (1), имеющее единственное решение
индекса к = -1 при выполнении условия
(4.28)
-1 л/1 ~ Z0
Поэтому уровень рассогласованности системы (5) без /Ои должен
понижаться, и, следовательно, если эта система с /Ои совместна, то limy0 --0.
П->00
Но это надо доказать. Таким образом, /Ои делает определенной систему (5) и
поэтому будем называть ее регуляризирующим фактором. Для нахождения
/Ои опять воспользуемся правилом Крамера:
/0.	=	+	(4.29)
Таким образом, = 0 тогда и только тогда, когда решение индекса к = -1
И-»®
для уравнения (1) существует. Следовательно, поведение уОл при расчетах
является индикатором наличия решения индекса к = -1.
Замечание 4.2. Во многих приложениях (в аэродинамике, теории
упругости и т.д.) часто требуется вычислить не саму функцию <?(/), а
интеграл |у/(г)ф(О«й, где на [-1,1]. Из неравенств (7) и (8) следует, что
-1
для вычисления указанного интеграла можно воспользоваться формулой
прямоугольников по точкам tk, k = причем брать в этих точках не
функцию ф(0 , а значение ф„(^), т.е. выполняется неравенство
42
(430)
Действительно, имеем
1	Я	"1+Л	I
|И')ф(')** - S)'•’» (*» >А - J И0ф(0<* +
+ /И0ф(0Л-^И'*)ф('*)А + Ё^»)(ф(^)-ф,(/*)> s
-1+Л	Ы	*=1
£=1
где М = тах|<//(0|.
Теперь рассмотрим метод дискретных вихрей для численного
решения характеристического интегрального уравнения первого рода с
логарифмической особенностью на отрезке. Это уравнение теперь
запишем в виде
1 jln|z0 - t\<p{t)dt = /(/„), t, 6 (-1,1).	(4.31)
Будем предполагать, что f '(t}eHa на отрезке [-1,1]. Используя теперь
рассуждения п.2 при получении формулы обращения уравнения (2.40)
получаем, что уравнение (31) эквивалентно в смысле разыскания решения
системе
=	(432)
* -1 h -t
= C =	(4.33)
Поэтому для численного решения методом дискретных вихрей уравнения
(31) надо рассмотреть следующую систему линейных алгебраических
уравнений (СЛАУ) (см. (4))
=	=	(4.34)
Iqj
Справедлива следующая
Теорема 4.2 Пусть f '(t)&Ha на отрезке [-1,1]. Тогда между решением
^(0 уравнения (31) и решением СЛАУ (34) выполняется соотношение (6) в
котором величина вп (tk), k = 1,..., п 9 удовлетворяет соотношениям (7) и (8).
Для доказательства теоремы 2 заметим, что левая часть системы (34)
совпадает с левой частью системы (4). Поэтому имеем
43
Л ;=i л

#’(<о;)* ! 1
1 JM( 1 \f(t)dt
+/w("in2j7fZ7
Теперь, как и при получении формулы (27), имеем
1	1	1 'rf(t)dt
-1 (к - 'о In2 я-yjl- tl j 4\- t*
+^и(^),хе(-1,1),
(435)
(436)
где величина 0n{tk), к = 1,...,п, удовлетворяет соотношениям (7) и (8).
Сравнение с формулой (2.61) и доказывает теорему 2, так как из этой
формулы следует, что точное решение уравнения (31) в точке tk получается
по формуле
^*)= —*=»•-"• <4-37)
-tk -1 -Л ln2
Обратимся теперь к численному решению характеристического
гиперсингулярного интегрального уравнения первого рода на отрезке.
Запишем это уравнение в виде
-	h е (-1.1)•	(4.38)
*-iQo-O
Как видно из формулы (3.62), уравнение (38) имеет единственное решение.
Пусть множества £ = {/г,лг~ 1,и Ео = {гО7,у = 0,1,...,л} образуют такое
разбиение отрезка [—1,1] с шагом Л, что tk =-1 + (Л-1)Л, z07=r7+p
2
£ = 1,...,и;у = 1,1 , h = -—. Теперь заменим уравнение (38) следующей
п-1
СЛАУ
L ^<1 = j = l,-,n-l.	(4.39)
Я-*=1	' Vo;-О
Исходя из определения гиперсингулярного интеграла на отрезке, систему
(39) запишем в виде
--------—5—] = Ж), / =	(4.40)
toJ~'h+} hj-fk
Как следует из результатов [2] для квадратурных формул метода дискретных
вихрей для гиперсингулярного интеграла на отрезке, СЛАУ (40)
аппроксимирует уравнение (38) так же как СЛАУ (3)-(5) аппроксимируют
уравнение (1).
Оказывается справедлива следующая теорема.
Теорема 4.3. Пусть функция f(t)eHa на отрезке [-1,1]. Тогда между
решением системы (40) и решением уравнения (38) выполняется
соотношение
к = \,...,п-\,	(4.41)
а так же соотношение
44
(4.42)
I	" I
где полагается ^и('оо) = Ри(*ой) = °, a величина On(tk),k = 2,...,и-1,
удовлетворяет соотношениям (7), (8) для к = 2,..., п -1.
Доказательство. Вначале заметим, что система (40)
системе
1	h _ А. ) •
h taj-tk
y<Pn^k)-<Pn^-x) h = 0
S h
эквивалентна
(4.43)
п-\,
где, напомним, ₽„((«,) = <»„(/„„) = 0.
Система (43) совпадает с системой (4) и поэтому имеем
. 1,. 1........................................(144)
h	h j^h tk tQj
Теперь сравнение формулы (44) с формулой (2.67) и теорема 1 дают
доказательство соотношения (42). Далее из формулы (44) следует, что
*=’..."•	(4-45)
/=1 п у=1 П t, -lQj
Теперь учет формул (2.69) и (25), (26) показывают справедливость
соотношения (41).
Теперь сформулируем метод дискретных вихрей для численного
решения полного сингулярного интегрального уравнения первого рода
на отрезке, т.е. для уравнения
f*(W)4>('H = Я'о) •	(4-46)
-1 *0	* -I
Будем предполагать пока, что функции f(t) и k(tOft) принадлежат классу Яна
своих областях определения.
Разрешив уравнение (46) относительно его характеристической части,
получим, что оно эквивалентно в смысле разыскания решений индекса
г уравнению типа Фредгольма второго рода [1]
Ф(О+ /^г(г,т)ф(т)^т = /и(0-	(4.47)
где
*	-1	‘-1о
причем Г, = л, То = Т_j = 0 и при к = -1 должно выполняться условие
fx:!(t)f(t)a = jf f/?.-;(r)t(/,T)<*	.	(4.48)
45
Отметим, что ядро N^t) имеет вид
<4-49>
где а и р равны 0 или 1/2, а функция Ф(/,т) непрерывна на множестве [-
1,1]х[-1,1]. Отметим, что к интегральному уравнению второго рода с ядрами
вида (49) применима полностью теория Фредгольма построения
приближенного решения и нахождения его решении [35 из [2]]. Однако
уравнение (47) можно непосредственно свести к уравнению Фредгольма
второго рода с непрерывным ядром с помощью соответствующей замены
переменной. Это замечание позволит систему линейных алгебраических
уравнений для уравнения (46), получаемую в рассматриваемом численном
методе, эквивалентным образом преобразовать в систему линейных
алгебраических уравнений для уравнения Фредгольма второго рода,
эквивалентного уравнению (46) в данном классе решений.
Справедлива следующая
Теорема 4.4. Пусть в уравнении (46) функции f(t) и k(t^t)
принадлежат классу Н на множествах [-1,1] и [- 1,1]х[-1,1]
соответственно и это уравнение имеет единственное решение в
соответствующем данному индексу к классе функций (для к = 1 считаем
заданным значение интеграла от решения). Тогда между решениями систем
линейных алгебраических уравнений
= J =	(4.50)
t-l *07 '»	*-1
+	= /ОоД J = 1.."-1,
.°' ‘	(4.51)
^ф„(4)Л=С,
.*=1
Го.+Хт1^+£*('<17>/‘ХО*)л=/Оо;)> 7=0,1,...,и,	(4.52)
*=1 Ч Л=1
и соответствующими решениями уравнения (46) выполняется соотношение
(6), в котором величина 6(tk) удовлетворяет неравенствам (7) и (8). Здесь
множества Е =	= и Ео = {/Оу,у = 0,1,...,и} образуют каноническое
разбиение отрезка [-1,1].
Доказательство. В системах (50>(52) оставим слева слагаемые,
соответствующие характеристическому сингулярному интегральному
уравнению, а все остальное перенесем вправо. Используя результаты
теоремы 1 получим, что рассматриваемые системы эквивалентны системам
(*=0,1,-1)
=	(4-53)
т=1
где
46
n	lk~4j
7Ah)=-W1%Ш^-ткс,
n U=1«-r) lk~TOj
определение £(x) cm. (10).
Дальнейшее доказательство проведем более подробно для к - 0, так как
в остальных случаях оно аналогично. Из формулы (53) видно, что если
умножить обе части системы (50) на множитель (1-^)1/4(1+/А)3/4, потом
произведение Ф„«Э на этот множитель обозначить через $„(?*) и
рассматривать вновь полученную систему линейных алгебраических
уравнений, то она аппроксимирует интегральное уравнение Фредгольма
второго рода с ограниченным ядром
ФЖ)+	=ХЛ)>	(4.54)
о
где
Л) =	г(7, )Х1 - Г(Г, )>'*(! + /(О г,
ХЛ)=+г(оГ,
?(л)=ф(г(г1 )Xi - r<ri )УМ (i+^(^1 )У/4,
,	_ г
1 (1-г),/4(1 + г)5'4 ’ Л J(l-r)1M(l + 7)3'4 ~Tl(T) ’
'f dt
J(l-r)1/4(l + r)3/4 *
Причем, как следует из формулы(53), порядок аппроксимации будет иметь
вид
+1)-'4 +&^+—(4.55)
Из теории численных методов для интегральных уравнений
Фредгольма второго рода [49 из [2]] с непрерывным ядром следует, что
порядок аппроксимации Ф„(Л) функции Ф(/) такой же. Возвращаясь теперь к
функциям ф(0 и фя(О, получаем справедливость сформулированной
теоремы для к = 0. Для к = 1 и -1 она доказывается аналогично.
Опираясь на теорему 4 рассмотрим теперь метод дискретных вихрей
численного решения характеристического сингулярного интегрального
уравнения первого рода на системе отрезков, т.е. уравнения
= /(',),	(4-56)
L *0	*
где L является совокупностью / штук непересекающихся отрезков
[4,
Идея дальнейших рассуждений будет состоять в представлении
уравнения (56) как системы сингулярных интегральных уравнений на отрезке
[-1,1]. Поэтому, в соответствии с терминологией в теории таких систем [201
47
из [1]], будем говорить, что решение ф(г) уравнения (56) имеет индекс
=	-1,0-1,m = l,...,Z если оно: не ограничено на обоих концах; не
ограничено на одном конце; ограничено на обоих концах отрезка [Лда,Ям].
Будем это решение обозначать через ф*(0- Рассмотрим отображение gM(r)
отрезка [-1,1] на отрезок [Ат,Вт\, где
и =	(4.57)
Обозначим
<₽,. (о=ч>, (4^
Будем использовать равномерное разбиение на каждом из отрезков
= На отрезке [Ат,Вт] выберем каноническое разбиение с
шагом h множествами Ет = {tk9k=l,...,nm } и Ет0={с0„/ = 0,1,.т = 1,...,/.
tn	т V mtK7	7	7 Щ*	m4U v	77	7 mJ 7	7	7
Тогда справедлива
Теорема 4.5. Пусть функция f(t) &Н на L. Тогда между решением
системы линейных алгебраических уравнении
(4.59)
=№«)Ст> У =	»1 = 1,
*=1
и решением ф/г) уравнения (56), для которого известны значения
интегралов по тем отрезкам, составляющим L, на которых оно имеет
индекс 1, выполняется соотношение (6)/
Доказательство. С помощью отображений (57) уравнение (56)
можно рассматривать как систему I сингулярных интегральных уравнений на
[-1,1], которая при соответствующих дополнительных условиях (известно
значение интеграла от решения на тех отрезках из Z, на которых решение не
ограничено на обоих концах, т.е. на которых оно имеет индекс 1), имеет
единственное решение. Следовательно, эта система эквивалентна системе
интегральных уравнений Фредгольма второго рода, которая также имеет
единственное решение. Поэтому, повторив в дискретном виде процесс
перехода к системе интегральных уравнений Фредгольма второго рода,
получим, что система линейных алгебраических уравнений (59) эквивалентна
системе линейных алгебраических уравнений для этой системы
интегральных уравнений Фредгольма второго рода. Такой переход возможен
в силу того, что разрешимы системы (3)-(5) при любом к = 1,0-1.
Рассмотрим теперь применение метода дискретных вихрей для
численного решения характеристического сингулярного интегрального
уравнения первого рода на отрезке, когда в правой части имеется дельта
функция, т.е. к уравнению
48
= fM+QS(t.-q), t,,q € (-1,1).
-1 *0	*
(4.60)
Таким образом, уравнение (60) рассматривается в пространстве обобщенных
функций. Используя формулы обращения для уравнений (3.47)-(3.49),
соотвтствующие решения уравнения (60) можно записать в виде
ф(/) = —fc\t0)£^+R;'(q)-2—vrC*	(-1,1),(4.61)
«	t -10	t — q
где Vi = 1, v0 =	= 0,
яо(О = А(<)=Л^> Л,(О=А(О = -Д=, Л1(О = Л,-,(О = Р2(О.
п + 1	л/1-r
Напомним, что решение индекса к--1 существует при выполнении условия
1f/(Q-±g^-g)lft = O.	(4.62)
-1	-71-г2
Из равенства (62) следует, что если функция	на отрезке [-1,1], то
решение индекса *-=-1 для уравнения (60) существует при
\4^л-	(4-63)
-1Л/1-Г2
Для применения метода дискретных вихрей для численного решения
уравнения (60) возьмем на отрезке [-1,1] множества Е и е09 образующие
каноническое разбиение этого отрезка, и будем предполагать, что точка
qeE0 при некотором j=j^, т.е. q= /0 , и введем функцию ^(г0^)по
правилу: 6h(t0-q)=±, t, efo-^ + ^L и Sh(t„-<?)=0, t<,e[q-^,q + ^]. Теперь
п	2.	2	2	2
заменим уравнение (61) следующей СЛАУ
-q)t J = l,...,n,	(4.64)
Чу ~lk
S7^=/M+etf*(', -?), 7=i,...,«-1,	(4.65)
*=1 l0j ~lk
k=\
/ол+^7!^ = Ло/)+е^(^-9),7=0,1,-,п •	(4.66)
*=1 lQj ~?k
Справедлива следующая теорема.
Теорема 4.6. Пусть функция	на отрезке [-1,1]. Тогда между
решением <р„(О системы (64), (65) или (66) и решением cp(t) уравнения (60)
соответственно индекса к =0,1,-1 выполняется соотношение
) - <Рп(h	k = 1, п>	(4-67)
в котором величина 0n(tk) удовлетворяет неравенствам
1) для всех tk е [-1 + 8,q -£]и+ £,1 -, где 8 > 0 сколь угодно мало,
0„(tt)<csh^, Л>0 ,	(4.68)
49
2) для всех точек tk е [-1,1]
£X(O/^CAS Л.>0,
Jt=l
где С5,С - некоторые константы, не зависящие от п.
Доказательство, Так как матрицы систем (64)-(66) такие же
соответствующих систем (3)-(5), то получаем
S2^(tk)+vKX-I^C, k = l
п
(4.69)
как у
И,
(4.70)
где
Л Ж(-ОЛ z*“*oy	h я
Запишем теперь решение <p(t) в формуле (61) в виде
F(0=«>u(0+^(0+-J?r(rKC,
Л
Q
~tQJ4
(4.71)
где
п	t-q
Теперь рассуждения в теореме 1 показывают, что модуль
М') и SiMr
S,
Я	-1	* ~ <0
разности
(tk) удовлетворяет соотношению (6), а модуль разности р2Д0 и
(tk) удовлетворяет соотношению (67). Теорема 6 доказана.
Теперь рассмотрим применение метода дискретных вихрей для
численного решения характеристического гиперсингулярного
интегрального уравнения первого рода на отрезке, когда в правой части
имеется дельта функция, т.е. к уравнению
-1 Vo О
Таким образом, уравнение (72) рассматривается в пространстве обобщенных
функций. Используя формулу обращения (3.62) для этого уравнения,
получаем
р(0 = -Г1п
яг '
(4.72)

(4.73)
Чтобы применить метод дискретных вихрей к численному решению
уравнения (72), как и для уравнения (38), пусть множества Е =	= 1,...,п} и
£0 =	j = 0,1,...,и} образуют такое разбиение отрезка [-1,1] с шагом Л, что
tk --1 + (к- 1)й,	Го =г+ —, к =	= 1,...,и-1	, Л- ^ . и
2	и-1
предполагать, что точка qeEQ при некотором т.е. q= , и
* = 1,
будем
введем
50
функцию 3„(t0-g)no правилу: <?„(/<,-9)=|, t0 E[q~,g + %], и 3„(to-g)=O,
п	2 Z
h h
/0 ё to -	• Теперь заменим уравнение (72) следующей СЛАУ
-(*.*)[—7----------—Ц-] = /(*<>,)+6<?Х, -«), j = 1 П-1. (4.74)
П *“1	*0;	**+1	*0/ ‘к
Оказывается справедлива следующая теорема.
Теорема 4.7. Пусть функция f(t)eHa на отрезке [-1,1]. Тогда между
решением системы (74) и решением уравнения (72) выполняется
соотношение
k = \,...,n-\,k*kq	(4.75)
а так же соотношение
ф'(^)~к = 2...„л-1,	(4.76)
где полагается ри(*оо) = ^и(*ои) = 0> а величина O„(tk),k = 2,...,и1,
удовлетворяет соотношению (67) для к = 2,..., п -1.
Для доказательства теоремы 7 надо воспользоваться доказательством
теоремы 3 с учетом специфики правой части уравнения (72) и
рассуждениями при доказательстве теоремы 6.
Рассмотрим теперь применение метода дискретных вихрей для
численного решения характеристического сингулярного интегрального
уравнения первого рода на отрезке, когда в правой части имеется
функция, имеющая в некоторой точке отрезка особенность типа 1/х.
Конкретно рассмотрим уравнение
=	(4.77)
-1 *о * Iq Я
Ясно, что решением этого уравнения будет функция <p(t) = te (-1,1).
Если мы желаем рассматривать для уравнения (77) решения индекса
к = 0,1,—1, то дельта функцию надо рассматривать как произведение RK(t)y/K(t),
J^(i) = P3(r) = J^, «,(») = А(Г) = -Д=, J?.,w = я;’(О = Р2(П, где функция
V/ + 1	л/1-/2
является элементом соответствующего пространства S\l2R*) обобщенных
функций, а правую часть в (77) надо рассматривать как элемент пространства
S *(12Г,) обобщенных функций. Заметим, что для уравнения (77) существует
решение индекса -1, так как верно равенство
к=4--------= 0, ^€(-1.1).	(4.78)
а общее решение индекса 1 задается формулой
р(0 = 8(t - q)+—.	(4.79)
TTyll-t2
Численное решение уравнения (77) рассмотрим в случае индекса 1
как наиболее интересным с точки зрения приложений. Возьмем на отрезке [-
51
1,1] множества Е и Ео, образующие каноническое разбиение этого отрезка, и
будем предполагать, что точка qeE при некотором к=Л?, т.е. q= tk .
Заменим уравнение (77) следующей системой
=	7 = 1...(4.80)
‘О/ *к *0j Q
2р„(4)Л = 1 + С.
к-1
Для упрощения рассуждений возьмем 0-1. Тогда в силу теоремы 1 имеем
ч>М=	тг , L;	(4.81)
Запишем последнюю формулу в виде
‘•‘Л <482)
Перепишем формулу (82) в виде
-рй’--------------------->Lr)’ к*к<\ =ч-
h tk-q%h tt-t0) q-tOj
Теперь из формул (25) и (26) и равенства [1]
(4.83)
= /,
(4-84)
Следует, что
**м», =«>
(4.85)
где величина удовлетворяет неравенствам (68) и (69).
Теперь рассмотрим равенство (83). Используя равенство [262 из [2]]
У—Ц- = —	(4.86)
й(2* + 1)2 8
можно доказать, что
?’»(Ч)=|+О(1)-	(4-87)
Введем теперь функцию ^(0, fe(-l,l), по правилу < (0
+	(»„•(/)=	Тогда проведенные выше
рассуждения показывают, что в смысле обобщенных функций, имеем
соотношение
lim <рп (г) = р(0 = 3(t -q)-.	(4.88)
ял/1-Г2
Наконец рассмотрим применение метода дискретных вихрей
для численного решения характеристического гиперсингулярного
52
интегрального уравнения первого рода на отрезке, когда в правой части
имеется функция, имеющая в некоторой точке отрезка особенность типа
1/х. Конкретно рассмотрим уравнение
/(г^=А’'“’9еН,1)-	(4-89)
-Wo "“О
Используя решение (3.79) для уравнения (3.74), получим, что решение
уравнения (89) будет задаваться по формуле
tp(t)=~ arcsin t + - F(8(t - q))	(4.90)
я 2
где £(£(/-<7))=0, при 0</<^, F(£(jtf-#))=i, при t=q, F(^(r-^))=1, при
q < t< 1.
Чтобы применить метод дискретных вихрей к численному решению
уравнения (89), как и для уравнения (38), пусть множества Е = {^,лг = 1,...,и} и
£0 = ^О7,7 = 0,1,...,и) образуют такое разбиение отрезка [-1,1] с шагом Л, что
=-1 + (£-1)А, Го =Г+-,	£ = l,...,n;j = l,...,«-l	,	£ = и будем
2	п-1
предполагать, что точка qeE при некотором к=^, т.е. q= tk . Заменим
теперь уравнение (89), как и уравнение (38), системой
, чг	1	1 з	1
Z^(z0i)[”—-----ГТГ3 = ГТ“’
А-1	*0j	*к+1 hj »*	Я
Повторяя теперь в дискретном виде процесс получения решения уравнения
(3.74) и доказательство сходимости численного решения к точному для
уравнения (77), получим, что, в смысле сходимости обобщенных функций,
(4.91)
имеем
lim <рп (/) = - lim [р* (т)dr = <p(f) =—arcsin t + - - F(8(t - q)),	(4.92)
»->=□	л-но J	Я	2
где это та функция, которая имеется в формуле (88).
Замечание 4.3. Чтобы говорить о скорости сходимости в
пространствах обобщенных функций, надо рассматривать эти уравнения в
паре соответствующих пространств Н^. Например, если ввести функцию
g(/) = -^SL, то уравнение (89) относительно функции g(i) надо
л/1-/2
рассматривать в паре пространств Н*_ и Н^х для 2 < |, так как правая часть
этого уравнения лежит в любом пространстве я*, 2<-^. Действительно,
будем рассматривать функцию f(t) - —как элемент пространства S'(l2 ),
t-я
тогда
/(O = -1-=Z/(«)t/:(O, /(«)= ]p2W(t)u;(i)dt=
t Я n=0	-1
53
(_L	?	^)-/2л =—J2n cos((n + l)arccos<7) ,п=0,1,...	(4.93)
л-Д t-g
Из формулы (93) следует, что для любого Л < -- функция
(4.94)
л=О
принадлежит пространству ,
Обратимся теперь к рассмотрению метода дискретных вихрей для
особых интегральных уравнений в периодическом случае.
Вначале рассмотрим применение метода дискретных вихрей к
численному решению характеристического сингулярного интегрального
уравнения первого рода на окружности, т.е. к уравнению
(4-95)
L *0 *
в котором L является окружностью радиуса единица с центром в начале
координат. Пусть множества E = {tk,k = 1,...,и} и Ео = {/Оу, у = образуют
каноническое разбиение окружности, т.е. точки tk разбивают окружность на
равные части, а точка является серединой дуги между точками t} и 17+1.
Справедлива следующая
Теорема 4.8. Пусть функция на L. Тогда между решением
системы линейных алгебраических уравнений
=	7 = 1.«.	(4-96)
i-1 “Ц
где ак = tk+i-tk9 tn+l = и решением уравнения (95) ( см. [1])
ф(г) = __1	(4.97)
я 1 t~‘«
выполняется соотношение
к = \,...,п,	(4.98)
в котором величина 0(tk) удовлетворяет неравенству
fl(G)SoQr) 0<Л<1.	(4.99)
Доказательство. Система (96) совпадает	с системой (3), если в
последней заменить h на ак. Поэтому те же рассуждения дадут
ф.(г»)=-1СЁг/йт^’	(4Л00)
ак j=l “j *к l0j
где -/07+]--rOj, £ =	^ол+1 ~^01 *
Так как теперь L- окружность, то к изучению множителей 1$ и 1$}
придется подойти иначе, нежели в теореме 1. Напомним, что
tk=e'\ t.k=e,{^nln} =eie^\ к = \,...,п.
Поэтому можно написать
54
(4.101)
Так как точки гА,Л = 1,...,и, разбивают окружность L на равные части, а
tOk является серединой дуги (/А+1,/А), то 0„-&к = 2л(т-к)/п,
0От-вк=2я(т-к)/п+я/п. Вводя теперь перенумерацию, с учетом
периодичности функции ёв можно написать
п-1	п-1
jX») =	р(и) —
т-0
Для вычисления	заметим, что числа ет2х1п, т = 0,1,...,п-1, являются
корнями и-й степени из числа z = 1, т.е.
zn -1 = fj(z - ёт2*1п) = (z- 1)П (z -J*2*'")
т-0	т=1
ИЛИ
-1 = z’-' + z"-1 +...+1 = JJ (г - е"2’'”).
Z~ 1	m=l
Последнее равенство есть тождество. Поэтому, устремляя z к единице, в
пределе получим
1Ьп^*=Я=/Й>.	(4.102)
Z->1 Z-1
Для вычисления Р2(;} заметим, что числа е'(*/л+ж'2,г/л), w = 0,1,..., и-1, являются
корнями и-й степени из числа z = -1, т.е.
z” + l = fj(z-e,(jr/”+w‘2jr/n)).
да “О
Последнее равенство верно при любом г, поэтому при z -1 получаем
Р2(л)=2.	(4.103)
Формулы (101)-(103) показывают, что
(4.104)
1	( . it . яЛ .1 У1)
=  —— sin — + l cos- = Z — + <91 - .
nsinWnV n n) n \nj
Аналогично можно показать, что
/а=ч--	(4Л05)
и bj п \п)
Из формул (104), (105)следует, что
(4Л06)
J=1 *к l0j J=l \nJ 1к l0j
55
(4.107)
Учитывая теперь результаты для квадратурных формул метода дискретных
вихрей на окружности [1], видим справедливость теоремы 8. Используя
теперь рассуждения, проведенные при доказательстве теоремы 8, можно
рассмотреть применение метода дискретных вихрей к численному
решению характеристического интегрального уравнения с ядром
Гильберта первого рода
2^fClg 2	=
Выберем на отрезке [0,2тг] точки 0к, к =	, которые, интерпретируемые как
точки единичной окружности Z, разбивают ее на п равных частей; ,
к = 1,..., п, делит пополам дугу (0к, 0А+1).
Напомним, что уравнение (107) имеет решение только при условии
j/(0)d0 = o,	(4.108)
о
которое и будем считать выполненным. Для выделения единственного
решения надо задать значение решения в некоторой точке либо значение
интеграла от решения (последнее более часто встречается в приложениях).
Справедлива следующая теорема.
Теорема 4.9. Пусть функция на [0,2 я], /(0)-/(2я), и для
нее выполняется равенство (108). Тогда между решением системы линейных
алгебраических уравнений
Го. + 7“т =
2-71 )t=i	2	П
1 ХГ /л \ п
2яы п
и решением ф(0) уравнения (107), задаваемым формулой
2я- nJ 2
(4.109)
(4.110)
при условии
1	25Г
±faff)d0 = C,	(4.111)
выполняется соотношение
|ф(0*) - Фй(0*)| О{пл In п),	(4.112)
где (ОД], если п произвольно и f(0)eН(а), Л = г+а, если п нечетно и
fr\0}eH(a}.
Доказательство. Просуммировав первые п уравнений в системе
(109) и учитывая равенство(108), получим
/.»=y-t/(^)-.	(4.113)
2Л т=1
Отсюда следует, что у0„ -> 0 при п -> оо тогда и только тогда, когда уравнение
(107) имеет решение.
56
Идея дальнейших рассуждений состоит в сведении системы (109) к
системе вида (96) для уравнения на окружности с помощью равенства
itk 1 вк	1 .	/Л 1 1 ЛЧ
77-7= о	,	V	(4.114)
h lo 2	2	2
Умножим последнее равенство в системе (109) на (-/) и прибавим ко всем
первым п уравнениям. Учитывая равенство (113), получим после умножения
обеих частей на я
i ХП /п \ 2л- 1	/л \
2	*=i И 2 к=1	2	п
2 к-1	п
или
= /(?0„), т = 1,...,п,
к=1 1Ьт~!к
гдс1„=е‘в‘,11>я = е1^,
(4.115)
1 и
Ф»('*) = Фл(^),^=2^/и,Л/о„)= #(^)-|Ё/(^)—
2g п _
Система (115) совпадает с системой (96) и поэтому ее решение дается
формулой (100). В силу равенств (104), (105) получаем
фж)=-ЛГт^,	(4.И6)
* «3
где Z>„ =2яй0ш/л.
Таким образом, опять воспользовавшись равенствами
£c«g ^'"У = 0, к=1,.. ,п
т=1	2
и (114), имеем
ГЛ \	1	1 , @к	_Г/Л \ 1 V4 Г/Л \ ^Л’ . -,1 2#
ф.(0»)=-—----------------------------------------;лСН =
Я- й\2	2 2J[L 2tf «J J П
+C.	(4.117)
2Л> m-1	2	П
Сравнивая формулы (НО) и (117), и учитывая свойства квадратурных
формул метода дискретных вихрей интеграла с ядом Гильберта [1], видим
справедливость теоремы.
Обсудим еще вопрос применения метода дискретных вихрей к
численному решению характеристического интегрального уравнения
первого рода с логарифмической особенностью в периодическом случае,
т.е. к уравнению
- fln|sin*^A(tW? = f(Oo), 0, е[0,2я].	(4.118)
I 2 Г
Как было показано в п.З, уравнение (118) имеет единственное решение для
любой правой части из пространства обобщенных функций S (Z2) и
эквивалентно в смысле разыскания решения системе равенств
57
~ jctg^g(0)dff = f(e„), 0oe[O,2>r],	(4.119)
i 2я	। 2x
— jg(0)d0 = ~±_ \f(0)d9.
Zn ’	4я1п2 *
Теперь в силу теоремы 9 для численного решения уравнения (118) по методу
дискретных вихрей надо взять СЛАУ
Го. +	—‘«.(Я*)2*- = f т =
„	1, ”	(4.120)
= -7-Г2
2п^х п 4я1п2 *
Теорема 9 дает оценку близости в точках вк между решением уравнения
(118) и решением системы (120).
Теорема 9 дает возможность рассмотреть так же применение метода
дискретных вихрей к численному решению характеристического
гиперсингулярного интегрального уравнения первого рода в
периодическом случае, т.е. к уравнению
~ J-----|-дг«9)^=/(во),0О)9е[О,2^,	(4.121)
4»oJsin*3>Z*
2
где функция /(^принадлежит множеству На на отрезке [0,2я-].
Уравнение (121), как и уравнение (107), имеет решение с точностью до
константы, а условием существования решения является равенство (108).
Поэтому, по аналогии с уравнением (107), для уравнения (121) методе
дискретных вихрей применяют следующим образом. На отрезке [0,2 п] берем
два множества точек Е = {ок,к -1 п} и Ео = fy0J j =	где точки
0к> к = 1,..п, которые, интерпретируемые как точки единичной окружности £,
разбивают ее на п равных частей; 0Qk, к = 1,..., п, делит пополам дугу (0к, 0к+}).
Теперь уравнение (121) заменяем СЛАУ
Го. +	=/(^) ,j=l.-,n, (4.122)
4jr^ /in2
2
где
Интегрируя интегралы в системе (122) запишем ее в виде
58
jrt	/7—о ________a
Го. +T-Sg„(6'0J(cZg-^!—i)=/(«0 ), j=l,.,.,n, (4.123)
2я” А=1	2	2
+Zg»(^)*=c.
2^r *ж1
Система (123) эквивалентным образом может быть записана в виде
Iй Q —Q	'bn
Го. +^Zc/«-?y-4-g^(^)Y=/(^)’j=1’-A	(4.124)
^£(-^(0.))—= 0,
2* ы п 2я И
где	g’M =	, gn{0M) = g,(^„).
2Я7 И
Действительно, первые п+1 уравнений в (124) совпадают с системой (109)
относительно переменных -g*(0*). Поэтому она имеет единственное
решение относительно этих переменных. Следовательно, величины
g„(0oi),* = l.представляются в виде g„(0o*)=4* +gA)>bl,.,n, где
- известные числа, и поэтому из последнего уравнения системы (124) найдем
g„(0o„), а значит и все остальные величины gn(0OJt),i = l,.Теперь из
приведенных рассуждений и теоремы 9 получаем справедливость следующей
теоремы.
Теорема 4.10. Пусть в уравнении (121) функция Д0О) принадлежит
множеству На на отрезке [0,2я] и удовлетворяет равенству (108). Тогда
между решением уравнения (121), удовлетворяющего равенству (111), и
решением системы (123) выполняются неравенства
* = U,«,	(4.125)
|g^)-g;^)|<02(+), к = 1,...,п,	(4.126)
где 0<	<1.
Рассмотрим теперь применение метода дискретных вихрей для
особых интегральных уравнений в случае, когда в правой части имеется
дельта функция. Начнем с уравнения (3.63) с ядром Гильберта, т.е с
уравнения
+ 2jctg^g(ff)d0 = g(ffoЮМ,	(4-127)
2.Л £	2	2.П
59
которое имеет решение
g(0) = -2-^^+£.	(4.128)
Возьмем на отрезке [0,2л-] два множества точек Е = {<9k,к = 1,и
eq = Ku = 1,•••>«}, как для уравнения (121), и будем полагать, что q = 0Qj .
Теперь уравнение (127) заменим СЛАУ
Yo.	....п, (4.129)
2 л £i	2	п	2л
E±g,(0t)h = C,
к=1
где функция <?А(0Оу) определена сразу после формулы (3.71). Поскольку левая
часть системы (129) совпадает с левой частью системы (109), то ( см. (117))
*.(*.)=	+ СЛ = 1,..,л.	(4.130)
2 л 7^	2	2 л п
В силу определения функции 8h(0Qj) и равенства
>i 2 п
окончательно получаем равенство
g,(^)=-2-c®^+c.	(4.131)
2.Л Z.
Из сравнения формул (128) и (131) видим соотношение между точным и
численным решениями.
Рассмотрим еще гиперсингулярное интегральное уравнение в
периодическом случае, когда в правой части имеется дельта функция
~ f-----^—ng(0)d0=g(0o	(4.132)
2
Используя множества Е = {0к,к = и £0 = {0О/ j = как для уравнения
(127), для численного решения уравнения (132) надо взять СЛАУ
Го»	j=ls---,n, (4.133)
2л ji-i	2	2	2л
60
M 4=]
Рассмотрим еще применение метода дискретных вихрей для особых
интегральных уравнений в случае, когда в правой части имеется
функция J-c/g——-.Начнем с уравнения ядром Гильберта
2л 2
у- ]ctg^-^-g(0)M=^-ctg^—^-, 0„,q е [0,2л-].	(4.134)
Ясно, что решением уравнения (134) является функция
g(0)=<?(0-9) - А. + С,в,Ч е [0,2л-].	(4.135)
2Л
Для численного решения уравнения (134) возьмем на отрезке [0,2 л-] два
множества точек Е = {0k,k = и EQ = fy0Jj =как для уравнения (121),
и будем полагать, что q = 0к, т.е. q^E. Нам удобно будет полагать, что
q е (0,2яг). Теперь заменим уравнение (134) СЛАУ
/о.	,j=l,...,n,	(4.136)
2Л	2	П 2Л 2
^g^}h = C.
2П
В силу симметрии точек 0Qjj = 1,...,и, по отношению к любой точке
0А,£ = 1,...,и, просуммировав первые п уравнений системы (136), получим
равенство
Гол =0-
Поскольку левая часть системы (136) совпадает с левой частью системы
(109), то (см. (117))
«.(*.)=	(4.137)
2л yz 2	2л-	2 п
Для вычисления суммы Sn k используем формулу
O-0q 0Q-q в-q,	<90-а\ .
Cfg —— ctg	= - ctg	(ctg -------Ctg +1
при 0 ф q, в силу которой имеем
61
(4.138)
2л 2 2л 2	2 п 2л 2л
При к = kq, т.е. 0к -0к имеем
„ 1 1
2?г	2?г
(4Л39)
2 n 2л 2л . 2 % - &кв n
sm ----------------------
2
Iй 1
^sin2 &QJ ~&к
2л 1
---7Г’ Л=Л-
п 2л
2
Из выражений (138), (139) следует, что сумму Sn k можно представить в виде
2л
8пк^ к*к , 8пк =ДгУ / - л , к = ка.
n,k	q, n,k4	е -ек п9 «
sin2 ——--L
2
Покажем, что последовательность функций 8п,к(0) = 8пк> Oe(0Qk^fOQk),
к = \,...,п, где полагаем 0^ =0Qn, сходится к функции 8(0 ~q). Действительно,
имеем
„ _ n 1 Л 1 4лг2
----------------"Z---п---г-пО(1)-++<*) при п-^ ОО,
• 2л 4л . 7&Qj -vk n
S1n -----------------~
2
где через 0(1) обозначена величина порядка 1, т.е. ограниченная величина.
Далее, так как
при х -> 0, то, учитывая известное равенство [1, с. 334]
1	^л2
Й(2Л + 1)2" 8 ’
получаем, что
62
р».*(0)^=tfj ыще =	= i+0(l).
0	й и 4/r2~ . 2 0^-0kq n2 n
sm 2
Введем теперь функцию g;(<9) = g/^),6>e(6>0^,6>0J, Л = 1,...,и. Тогда
проведенные выше рассуждения показывают справедливость следующей
теоремы.
Теорема 4.11. Для решения g(0) уравнения (134), определяемого
формулой (135), и функции &(0), получаемой из решения g„(0k), к =1,...,м,
системы (136), выполняются соотношения
|g(0J-g;(M=lsW-S.<M=O> * = U,», ***,,	(4.140)

П вЧк	0Q1	1
^>4-1	3)4-1	М
(4.141)
Наконец рассмотрим гиперсингулярное интегральное уравнение
первого рода
0 sm ----
2
(4.142)
Для численного решения уравнения (142) надо взять множества Е и Ей, как и
для системы (136), и взять следующую СЛАУ
Го» +	g^^-ctg^=^ctg^l, j=l,...,n, (4.143)
2я-	2	2 In 2
T-L«-(0«* )a=c-
zn
Исследовать систему (143) надо так же как систему (123) с учетом
особенности правой части.
В заключение рассмотрим метод дискретных вихрей численного
решения задачи обтекания кусочно-гладкого простого контура.
Исторически метод дискретных вихрей был вначале построен для
циркуляционного и без циркуляционного обтекания идеальной несжимаемой
жидкостью гладкого разомкнутого контура. Математическое обоснование
этого метода, в этом случае, было дано в монографии Белоцерковский С.М.,
Лифанов И.К. Численные методы в сингулярных интегральных уравнениях,
М., Наука, 1985, 256 с. Однако потребности практики побудили авторов
указанной монографии со своими учениками, на основе численного
63
эксперимента и сравнении результатов этих экспериментов с результатами
физических экспериментов, распространить этот метод на широкий класс
контуров. Однако математического обоснования в этом общем случае пока
нет. Решение этой задачи было бы хорошим трамплином для молодого
математика в численных методах.
Пусть контур L обтекаемого профиля задается параметрически:
x = x(r), y=y{t\ re[0,7], и он находится в стационарном потоке идеальной
несжимаемой жидкости. Выполняя условие не протекания в точках
М0(х0,у0)=Л/0(х(г0),Я*0)) контура L, приходим к уравнению
(4Л44)
2*°
Если контур L является гладким разомкнутым, то уравнение (144) является
сингулярным интегральным уравнением первого рода на отрезке вида (46) и,
следовательно, метод дискретных вихрей в этом случае строится
следующим образом. На отрезке [0,7] параметра I берут каноническое
разбиение, состоящее из множеств:	E = {tk,k = \,...,n},
tk=kh,	h = l/(n + V), и £0 =(Г0т,/и = 0,1,...,и}, tQm=tm+h!2, m = 0,l,...,n. В
точках Mk=(xktyk),	xk=x(tk), yk=-y(tk) контура L помещаем
дискретные вихри интенсивности Гк, а точки М0т(х0т,у0т), /и = 0,1,...,и,
хОт=х(гОда), Лт=>’(7от) берем расчетными. Теперь также как для тонкого
слабоизогнутого профиля, т.е. для уравнения (1) в зависимости от
рассматриваемой задачи: циркуляционной, бесциркуляционной или
безударной (если она осуществима), уравнение (144) заменяем
соответственно системой линейных алгебраических уравнений
т = 1,(4.145)
Л-1
я	(4.146)
ЁЛ=0, т = п,
Л=1
7о.+ЕЛ®'=/.> m = 0,l,...,n,	(4.147)
Л=1
где
=------------------------;  —  ........,
2fl- /v'2 . ./2	2
Уот,1 ’
= («О» -^)2 + (л- -л)2. Л. = -и
Система (145) дает решение уравнения (144) неограниченное на кромке
профиля, соответствующей значению параметра t = 0. Признаком
существования безударного обтекания является стремление к нулю при
п ~>оо регуляризирующей переменной уОп.
В рассматриваемом случае теорема 1 дает математическое обоснование
метода дискретных вихрей для гладкого разомкнутого профиля,
64
описываемого системами (145)-(147). Эти системы получаются из
выполнения условия не протекания от системы дискретных вихрей и
набегающего потока в расчетных точках, выбор которых диктуется Б-
условием метода дискретных вихрей [1].
Замечание 4.3. После решения систем (145)-(147) можно находить
все аэродинамические характеристики обтекаемого профиля, используя в
дискретном виде соответствующие формулы гл.9 в [1]. При этом, если в этих
формулах используется функция у(х), то дискретизацию формулы удобнее
проводить по точкам Мк и полагать ук = y(tk) = rk/ESk9 где &Sk- расстояние на
контуре L между точками и MQk, а если используется функция
= то дискретизацию формулы удобнее проводить по точкам MQk и
к
полагать gw =	,так как g(0) = 0 в циркуляционной задаче и g(0) = g(Z) -0 в
Л=1
бесциркуляционной.
Замечание 4.4. Если контур L является простым разомкнутым, но
кусочно-гладким, то точки Мк расположения дискретных вихрей надо
выбирать так, чтобы угловые точки контура L входили в их число (Рис. 4.1).
Рис. 4.1 При моделировании тонкого Рис. 4.2 Система дискретных
профиля с угловой точкой дискретный вихрей и расчетных точек на
вихрь	замкнутом
помещается в эту точку	гладком профиле
Пусть теперь контур L является гладким замкнутым (Рис. 4.2), т.е.
х(0) = х(1), у(0) = у(1) и орт касательного вектора непрерывен. Тогда уравнение
(144) является сингулярным интегральным уравнением с ядром Гильберта
вида (3.7). Его решение определено с точностью до константы, которая
определяется тем, что в аэродинамике в этом случае рассматривается
бесциркуляционная задача, и поэтому
pds = O.	(4.148)
Поэтому в соответствии с результатами для уравнения с ядром Гильберта
метод дискретных вихрей в этом случае строится следующим образом. На
отрезке [0,/] теперь берем следующие два множества точек: E = {tk,k =
tk=(k-\)h, k = \,...,n9 h-Un И Eo = (ГОм,ти = 1,...,и}, tOm =tm +h/2, ю = 1,...,и.
Теперь опять в точках Mk(xk,yk) помещаем дискретные вихри интенсивности
65
Гк9 а точки MQm(xQm,yQm) берем расчетными точками. Уравнение (144)
заменяем следующей системой линейных алгебраических уравнений:
т=1,£л = о, (4.149)
*=1
где со”- такая же, как в системах (145)-(147).
66
Содержание.
Вопросы для сдачи экзамена по спецкурсу.
1.	Интегральный оператор с ядром Гильберта. Спектральные соотношения
([2,с. 43-44])................................................2-3
2.	Характеристическое интегральное уравнение первого рода с ядром
Гильберта. Формула обращения..................................3-5
3.	Интегральный	оператор	с	логарифмической	особенностью	в
периодическом случае. Спектральные соотношения ([1, с. 105-107])...5-6
4.	Интегральный	оператор	с	гиперсингулярной	особенностью	в
периодическом случае. Спектральные соотношения ([1, с. 116-117]).. .6-8
5.	Характеристическое интегральное уравнение первого рода с
логарифмической особенностью в периодическом случае. Формула
обращения......................................................8-9
6.	Характеристическое интегральное уравнение первого рода с
гиперсингулярной особенностью в периодическом случае. Формула
обращения.....................................................9-10
7.	Интегральный оператор с ядром Коши на отрезке. Спектральные
соотношения в классах функций, обращающихся в нуль или в
бесконечность на концах отрезка ([2, с.46-47])...............10-12
8.	Интегральный оператор с ядром Коши на отрезке. Спектральные
соотношения в классе функций, обращающихся в нуль на одном из
концов отрезка и в бесконечность на другом конце отрезка......12-13
9.	Интегральный оператор с логарифмической особенностью на отрезке.
Спектральные соотношения в классе функций, обращающихся в
бесконечность на концах отрезка..............................13-14
10.	Гиперсингулярный интегральный оператор на отрезке. Спектральные
соотношения в классе функций, обращающихся в нуль на концах отрезка
([2, с. 55-56, Теорема 4.2.4])...............................14-16
11.	Характеристическое интегральное уравнение первого рода с ядром
Коши на отрезке. Формула обращения............................16-19
12.	Характеристическое	интегральное	уравнение	первого	рода	с
логарифмической особенностью. Формула обращения.............19-20
13.	Характеристическое	интегральное	уравнение	первого	рода	с
гиперсингулярной особенностью на отрезке. Формула обращения...20-21
14.	3адача обтекания профиля с отсосом внешнего потока, сингулярные
интегральные уравнения и обобщенные функции [3, с. 404-405]..22-25
15.	Один вариант обобщенных функций на Гильбертовых
пространствах................................................25-28
16.	Характеристические особые (с логарифмической особенностью,
сингулярные и гиперсингулярные) интегральные уравнения первого
рода в пространствах обобщенных функций. Формулы
обращения....................................................28-32
17.	Некоторые точные решения характеристических особых интегральных
уравнений первого рода в пространствах обобщенных функций.....32-35
67
18.	Метод дискретных вихрей численного решения характеристического
сингулярного интегрального уравнения первого рода на отрезке ([1, с.
341-345]).....................................................36-41
19.	Метод дискретных вихрей численного решения характеристического
интегрального уравнения первого рода с логарифмической
особенностью на отрезке.......................................41-42
20.	Метод дискретных вихрей численного решения характеристического
гиперсингулярного интегрального уравнения первого рода на
отрезке........................................................42-43
21,Метод дискретных вихрей численного решения полного сингулярного
интегрального уравнения первого рода на отрезке.............43-45
22.	Метод дискретных вихрей численного решения характеристического
сингулярного интегрального уравнения первого рода на системе
отрезков ([1, с. 364-365])....................................46-47
23.	Метод дискретных вихрей численного решения характеристического
сингулярного интегрального уравнения первого рода на отрезке, когда в
правой части имеется дельта функция ([4]).....................47-48
24.	Метод дискретных вихрей численного решения характеристического
гиперсингулярного интегрального уравнения первого рода на отрезке,
когда в правой части имеется дельта функция ([4]).............48-48
25.	Метод дискретных вихрей численного решения характеристического
сингулярного интегрального уравнения первого рода на отрезке, когда в
правой части имеется функция, имеющая в некоторой точке отрезка
особенность типа 1/х. ([5])...................................48-51
26.	Метод дискретных вихрей численного решения характеристического
гиперсингулярного интегрального уравнения первого рода на отрезке,
когда в правой части имеется функция, имеющая в некоторой точке
отрезка особенность типа 1/х ([5])............................51-52
27.	Метод дискретных вихрей численного решения характеристического
сингулярного интегрального уравнения первого рода на окружности ([1,
с. 366-367, теорема 17.4.1.]).................................52-54
28.	Метод дискретных вихрей численного решения характеристического
интегрального уравнения с ядром Гильберта первого рода ([1, с. 371-375,
теорема 17.5.1 ]).............................................54-55
29.	Метод дискретных вихрей численного решения характеристического
гиперсингулярного интегрального уравнения первого рода в
периодическом случае..........................................56-57
ЗО.	Метод дискретных вихрей численного решения характеристического
интегрального уравнения с ядром Гильберта первого рода, когда в
правой части имеется дельта функция...........................58-58
31.	Метод дискретных вихрей численного решения характеристического
гиперсингулярного интегрального уравнения первого рода в
периодическом случае, когда в правой части имеется дельта
функция.......................................................59-59
68
32.	Метод дискретных вихрей численного решения характеристического
интегрального уравнения с ядром Гильберта первого рода в случае,
когда в правой части имеется функция —ctg— ~q ([6, п.4])........59-61
2яг 2
33.	Метод дискретных вихрей численного решения характеристического
гиперсингулярного интегрального уравнения первого рода в
периодическом случае, когда в правой части имеется функция
^‘g^1...........................................................61-61
2л 2
34.	Метод дискретных вихрей численного решения задачи обтекания
кусочно-гладкого простого контура ([1, с.453-455]).............62-64,
Литература
1.	Лифанов И.К. Метод сингулярных интегральных уравнений и
численный эксперимент, Москва, ТОО «Янус», 1995.
2.	Довгий С.А., Лифанов И.К. Методы решения интегральных уравнений.
Теория и приложения, Киев, Наукова Думка, 2002.
3.	Вайникко Г.М., Лифанов И.К., Полтавский Л.Н. Численные методы в
гиперсингулярных интегральных уравнениях и их приложения,
Москва, «Янус-К», 2001.
4.	Вайникко Г.М., Лебедева Н.В., Лифанов И.К. Численное решение
сингулярного и гиперсингулярного интегральных уравнений на отрезке
и дельта функция, Математический сборник, 2002, т. 193, № 10, с. 3-16.
5.	Лифанов И.К., Ненашев А.С. Исследование некоторых
вычислительных схем для гиперсингулярного интегрального уравнения
на отрезке, Дифференциальные уравнения, 2005, т. 41, № 9, с. 1270-
1275.
6.	Лифанов И.К. Об одном случае численного решения особых
интегральных уравнений в периодическом случае, Дифференциальные
уравнения, 2006, т. 42, № 9.
7.	Гандель Ю.В. Введение в методы вычисления сингулярных и
гиперсингулярных интегралов, Харьков-2001, 92 с.
Учебное издание
ЛИФАНОВ Иван Кузьмич
ОСОБЫЕ ИНТЕГРАЛЬНЫЕ УРАВНЕНИЯ
И МЕТОДЫ ИХ ЧИСЛЕННОГО РЕШЕНИЯ
Учебное пособие
Издательский отдел
Факультета вычислительной математики и кибернетики
МГУ им. М.В. Ломоносова
Лицензия ИД N 05899 от 24.09.01 г.
119992, ГСП-2, Москва, Ленинские горы, МГУ им. М.В. Ломоносова,
2-й учебный корпус
Напечатано с готового оригинал-макета
в издательстве ООО «МАКС Пресс»
Лицензия ИД N 00510 от 01.12.99 г.
Подписано к печати 19.12.2006 г.
Формат 60x90 1/16. Усл.печ.л. 4,5. Тираж 100 экз. Заказ 899.
119992, ГСП-2, Москва, Ленинские горы, МГУ им. М.В. Ломоносова,
2-й учебный корпус, 627 к.
Тел. 939-3890,939-3891. Тел./Факс 939-3891.