ОБОБЩЕНИЕ МЕТОДА ПЕТРОВА-ГАЛЕРКИНА ДЛЯ РЕШЕНИЯ СИСТЕМЫ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ ФРЕДГОЛЬМА ВТОРОГО РОДА Текст научной статьи по специальности «Математика»

_ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА_

2023 • Математика. Механика. Информатика • Вып. 1(60)

«Математика»

Научная статья УДК 519.6

DOI: 10.17072/1993-0550-2023-1-5-14

Обобщение метода Петрова-Галеркина для решения системы интегральных уравнений Фредгольма второго рода

Н. К. Волосова1, К. А. Волосов2, А. К. Волосова2, М. И. Карлов3, Д. Ф. Пастухов4, Ю. Ф. Пастухов4

Московский государственный технический университет МГТУ им. Н.Э. Баумана, Москва, Россия

2 Российский Университет Транспорта, Москва, Россия

3Московский физико-технический университет, Москва, Россия

4Полоцкий государственный университет, Новополоцк, Республика Беларусь

Автор, ответственный за переписку: Дмитрий Феликсович Пастухов, [email protected]

Аннотация. Рассматривается задача численного решения линейной системы m интегральных уравнений Фредгольма второго рода. Впервые предложено обобщение проекционного метода Петрова-Галеркина для решения данной задачи. Количество координатных функций n двух линейно независимых систем может быть равно, больше или меньше m - количества интегральных уравнений. Преимущество данного алгоритма заключается в том, что он не чувствителен к малости параметров X в системе интегральных уравнений. Алгоритм требует правильного выбора двух линейно независимых систем координатных функций и их числа. Определена диагональная и антидиагональная задача. Для антидиагональной задачи из двух уравнений Фредгольма алгоритм решения сведен к матричному решению. Решены два примера для антидиагональной задачи из двух уравнений Фредгольма, в которых численные решения задачи совпадают с точными решениями системы. Доказаны две теоремы для достаточных условий корректности предложенных численных алгоритмов в двух случаях. В первом случае рассматривается антидиагональная задача с двумя уравнениями Фредгольма второго рода. Вторая теорема рассматривает условия корректности для диагональной задачи общего вида. Несомненно, предложенный алгоритм будет полезен в задачах механики и вычислительной математики.

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

Для цитирования: Волосова Н. К., Волосов К. А., Волосова А. К., Карлов М. И., Пастухов Д. Ф., Пастухов Ю. Ф. Обобщение метода Петрова-Галеркина для решения системы интегральных уравнений Фредгольма второго рода // Вестник Пермского университета. Математика. Механика. Информатика. 2023. Вып. 1(60). С. 5-14. DOI: 10.17072/1993-0550-2023-1-5-14.

Статья поступила в редакцию 05.01.2023; одобрена после рецензирования 27.01.2023; принята к публикации 15.03.2023.

«Mathematics»

Research article

The Petrov-Galerkin Method Generalization for Solving a Second Kind Fredholm Integral Equations System

N.K. Volosova1, K.A. Volosov2, A.K. Volosova2, M.I. Karlov3, D.F. Pastuhov4, Yu.F. Pastuhov4

University of the Moscow State Technical University named after N.E. Bauman, Moscow, Russia

Эта работа О 2023 Волосова Н.К., Волосов К.А., Волосова А.К., Карлов М.И., Пастухов Д.Ф., Пастухов Ю.Ф. лицензируется под СС BY 4.0. Чтобы просмотреть копию этой лицензии, посетите http://creativecommons.org/licenses/by/4.0/.

2 Russian University of Transport, Moscow, Russia

3 Moscow Physics and Technology University, Moscow, Russia

4 Polotsk State University, Novopolotsk, Belrussia

Corresponding author: Dmitriy F. Pastukhov, [email protected]

Abstract. The numerical solution problem of a second kind Fredholm m integral equations linear system is considered. The Petrov-Galerkin projection method generalization for solving this problem is proposed the first time. n is a coordinate functions number of two linearly independent systems. It can be equal, more or less then the m - the integral equations number. This algorithm advantage is that it is not sensitive to the parameters X smallness in the integral equations system. The algorithm requires the correct choice of two linearly independent coordinate functions systems and their number. The solution algorithm is reduced to a matrix solution for an antidiagonal problem. Two examples are solved for the antidiagonal problem, and numerical solutions of the problem coincide with the exact solutions. Two theorems are proved for sufficient conditions for the proposed numerical algorithms correctness in two cases. In the first case, an antidiagonal problem is considered with two Fredholm equations of the second kind. The second theorem considers well-posedness conditions for a general diagonal problem. Undoubtedly, the proposed algorithm will be useful in mechanics and computational mathematics problems.

Keywords: Fredholm equation; numerical methods; equations of mathematical physics; matrix; integral equations

For citation: Volosova N. K., Volosov K. A., Volosova A. K., Karlov M. I., Pastuhov D. F., Pastuhov Yu. F. The Petrov-Galerkin Method Generalization for Solving a Second Kind Fredholm Integral Equations System. Bulletin of Perm University. Mathematics. Mechanics. Computer Science. 2023;1(60):5-14. (In Russ.). DOI: 10.17072/1993-0550-2023-1-5-14.

The article was submitted 05.01.2023; approved after reviewing 27.01.2023; accepted for publication 15.03.2023.

Введение

Постановка задачи

В работах [1], [2], [3] описана постановка задачи для интегрального уравнения Фред-гольма второго рода и различные ее методы численного решения. Одним из важных методов является проекционный метод Петрова-Галеркина, применяемый для одного интегрального уравнения Фредгольма в конечномерных пространствах линейно-независимых координатных функций. Уравнение Фред-гольма второго рода подробно изучено в работе [4], в частности, для малых значений X и для больших значений параметра X.

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

В данной работе мы обобщили метод Петрова-Галеркина на систему интегральных уравнений Фредгольма. В последнее время уравнения Фредгольма второго рода активно применяются в конформных отображениях [5], в механике [6]. В работе [7] для численного решения уравнения Фредгольма с двойной точностью при небольших параметрах X применен метод замены интеграла в матричной форме.

Рассмотрим систему линейных уравнений Фредгольма 2-го рода, которую будем решать проекционным методом Петрова-Галеркина. Метод решения одного интегрального уравнения Фредгольма (1) описан в работе [1], [2], [3].

y(x) - Ц K(x, s)y(s)ds = f (x).

(1)

где: К (х, 5) е С ([а, Ь] х[а, Ь]) - непрерывная функция (кусочно-непрерывная на квадрате (х, 5) е [а,Ь] х [а,Ь] - ядро интегрального уравнения Фредгольма (1)), /(х) е С[а, Ь] - заданная непрерывная правая часть уравнения (1), X - известный параметр. Обобщим уравнение (1) на случай m интегральных уравнений и m неизвестных квадратично-интегрируемых функций у (х) е ¿2 [а, Ь], I = 1, т [1, 2]. Правые части уравнений (х), I = 1, т и параметры X,], ] = 1,т в задаче (2) заданы. Число интегральных ядер К ^ (х, 5) е С([а, Ь] х [а, Ь]), I, ] = 1, т в системе интегральных уравнений (2) равно

2

m .

m b _

yt (x) - £ Цj J Ki,j(x s)yj(s)ds = fi(xX i, j = lm . (2)

j=i

a

a

Определение 1. Задачу (2), у которой интегральные ядра имеют только равные индексы i = j, назовем диагональной. Очевидно, что каждая неизвестная функция y (x) в диагональной задаче входит только в свое i-е уравнение, а система интегральных уравнений разбивается на m независимых уравнений с разделенными функциями y (x), i = 1, m .

Например, явление квазиупругости в механике сводится к диагональной задаче из двух интегральных уравнений [8], в которой функции механического напряжения и относительной деформации связаны дополнительно линейным законом Гука.

Определение 2. Задачу (2), у которой интегральные ядра содержат произвольные индексы i < j(i > j), назовем смешанной. Если при этом отсутствуют ядра Кt(x,s) и параметры ХИ с равными индексами i = j , то задачу (2) назовем антидиагональной.

Сначала решим задачу (2) методом Петрова-Галеркина в общем виде. Затем в качестве частного случая рассмотрим антидиагональную систему двух интегральных уравнений.

Как и в работе [1], выберем две системы линейно независимых функций

{P(x)}"i=i,{wi(x)}h .

По первой системе координатных функций {р (x)}"i=1 запишем разложение неизвестных функций решения y. (x), j = 1, m , для удобства добавляя каждый раз функцию правой части системы (2) f (x), i = 1, m .

Решение системы уравнений (2) ищем в виде (3):

п ___

У,- (x) = f (x) + XCjPj (x), i = 1 m, j = 1,n . (3)

j=1

Коэффициенты разложения С\, i = 1, m, j = 1, n в формуле (3) неизвестны ( n■ m неизвестных). Подставим искомые функции y (x) в систему интегральных уравнений (2):

n m b ( n^ \

f (x) + X CjPj (x) J К,, (x, s)l fj (s) + X Ср (s) \ds =

j=i

j=i

= f (x) « R(x) = Xj(x) -

j=1

m b f n Л _

-Y^jJKjs)lfj(s) +Xj(s)\ds, k = 1,n . (4)

j=i

В формуле (4) функция (х), г = 1, т представляет собой невязку /-го уравнения в системе уравнений Фредгольма (2). Потребуем ортогональности невязки (х), г = 1, т всем функциям линейно-независимой второй системы:

Ж, (х)}"^ < ,у>= 0, г = 1, т, I = 1, п . (5)

Запишем формулу (5) - п ■ т условий более подробно с учетом формулы (4), в первой сумме заменим индекс ] на к.

п Ъ

XС'к \фк (х)У, (х)^х-

к=1 а

т ( п \

-ХЯ(х)Ки(х,з)1 /;(з) +Хфк(з) \dsdx = 0 о

]=1 а а V к=1 )

п Ъ _т 1г Ъ Ъ

X Ск\?к (хЖ (х^х! X Ск\\У1 (х)К1г j (х, з)Рк (з^х=

= (xKj(x,s)fj(s)dsdx, i, j = l,m,k,l = 1,n . (6)

j=1 a a

Введем в формуле (6) обозначения:

b b

ak ,i = JPk (xWi (x)dx j J (x, spk (s)dsdx

a a a b b _ _

fj'1 =JJ^i (x)Kt j (x,s)fj (s)dsdx,i, j = 1,m,k,l = 1,n .

a a

Тогда формула (6) примет вид

n j=m,k=n m

~ ~ i,l _ V 1

1СЛ1 - Ycj , i, j = 1,m; k,l = 1,n . (7)

j=i

В формуле (7) всего т ■ п условий для определения т ■ п элементов неизвестной матрицы коэффициентов разложения. Коэффициенты а ¡'к,, = 1,п , ; j,i = 1,т; к,I = 1,п , найдены, а матрица параметров Я,^j = 1,т задана в задаче (2). Система условий (7) содержит

г г,

-— . -— трехиндексные ] ■ и четырехиндексные коэффициенты. По нижним индексам j, к(^ проводится суммирование.

Как частный пример задачи (2) рассмотрим антидиагональную систему двух интегральных уравнений Фредгольма с двумя неизвестными функциями и(х), у(х):

b

u(x) — \ 2JК г(x,s)v(s)ds = f (x),

1,2 ¡^ 1,2 a b

v( x) — ^iJ К j( x, s)u(s)ds = f2 (x)

a

a

a a

k =1

l,k=1

a

a

Решение (8) ищем в виде (9)

п

1( х) = Л( х) + Х 1 (х),

3=1

п

(9)

v(х) = /2(х) + Х С]?3 (х)

=1

Невязки системы (8) имеют вид

п Ь ( п

Я =ХС\р1(х)-Л,21К12(х,5) /2(5) + ХСр(5)

а

3=1

3=1

п и п

К2 =£ С)р} ( х)-Х2,1 \ К 2,1 (х, 5) /1(5) +Х С1^ (5)

а

=1

Л •

=1

Условия ортогональности невязок имеют вид

<Я!,щ>= 0,к = 1п;<Я2,^>= 0,к = 1,п . (10) Запишем 2« условий формулы (10)

XС1 ]Ук(х)РзШх — Х,21Кц(х,5)Щк(х) /2(5) С)ф1 (5) й5йх = 0

1=1 о о V 1=1 >

п Ь Ь ( п ^

X С2 (х)р (х^х —Х211К2,1 (х, 5)^к (х) /1(5) + Х С'р (5) ё5ёх = 0

V я1 ;

Ь Ь

'к- .!.)- 2,1

а а

ЬЬ

С2

а а а а Ь Ь _

. =Ц К !(х, 5~)Рj (5)*Ук (х)^5^х, к,/ = 1, п .

а а

X а/;—х1,2 ХЬк,А2=Х1,2 /к2

1=1 1=1

п п

X ак,1'С] — Х2,1 X ^,1'С) = Х21/к =1 =1

(11)

Система двух уравнений (11) допускает матричную запись и ее решение:

(ас1 —а вс 2 = а р2 [ас 2—4 бс 1 = 4 р1 .

Решая (12), получим такой результат

(12)

С1 = А"1 (4 Р2 + А ВС2), АС2 — 4 БС1 =Хд Р АС2 — X б(А— (4 Р2 + А ВС2 ))=4 Р1«

с2 = (а — ХАл °А чв) ' (х,1 р1 + 4,4,2 ^а _1р2). (13)

Аналогично, получим, используя систему (12):

(ас1 —а вс 2 = а р 2 (ас 2—4 бс1 = x р 1, с 2 = а-1 x р 1 +4,1 бс 1) ас1—а вс 2 = а р 2« ас1 —а в(а—(4 р 1 + 4! бс 1 ))=а р 2« С1 = (а — ХгАВА-1Б)1 (аР2 +4гАВА^Р1) . (14)

В формулах (13), (14) элементы матриц А, В, Б и векторов Б1, Б2 определяются формулами (11).

Рассмотрим численный пример

1

и( х) — | х5У(5)^5 = х,

(15)

Для упрощения последних формул введем обозначения:

Ь Ь Ь

ак 1 =\^к(х)Р1(х)^х /к = 11К 2,1(х, 5)^к(х)Л (5)^х

у( х) — | (х + 5)м(5)^5 = 2 х — 2/3

0

А = А = 1, К г(х, 5) = х5, К21 (х, 5) = х + 5 , / (х) = х, /2 (х) = 2 х — 2/3, а = 0, Ь = 1.

Пример (15) имеет тот же алгоритм решения, что и задача (8). Вид функций / (х) = х, /2 (х) = 2х — 2/3 показывает, что линейно-независимые системы состоят из двух функций:

р (х) = 1, р (х) = х}, {^ (х) = 1,1^2 (х) = х}, т = п = 2 .

По формулам (11) находим матричные элементы для СЛАУ:

ак,1 = \Ук(хр(хУх к,1' =1,2,

0

1 1

а 1 = (хр (х^х = | ^х = х 0 = 1 ,

0 0

1 1 ^2 1

а12 (х)р2 (х)^х = | хйх = — =

0 0 2 0

1 1 х3 1

а22 = (хр2 (х^х = | х 2^х =—

0 0 3 0

1 1 _

Ь^; = Ц х5р (5)^ (х)dsdх, к,1 = 1,2 ,

0 0

11 11 ^ = Ц х5р (5)^ (х)dsdх =| | хsdsdх= -

0 0 0 0 2

Ь12 = 11 хяр (5)^ (х)dsdх = 11 хs2dsdх =— 0 0 0 0 3

1 1 1 1 2 Ь21 = Ц хр (5)^ (х)dsdх =| | х2-d-dх = —

а

2,1

0

00

00

0

Ь22 = Ц хяр2 (5)^2 (х)с5Сх = Ц х2я2 СяСх = — 0 0 0 0 3

1 1 _

= Ц(х + 5)р (5)щк (х)d5dх, к, 1 = 1,2 ,

0 0 1 1

¿и = [[(х + 5)d5dх = [| х5 + — | ¿х = | ^ + х | = 1

1 { 2 3 5 5

0 0 1 1

А = _[_[(х + 5)5d5dх = [| х— + — | ¿Х =| — + х| = —.

0 0 1 1

^ =[[(х + 5)хd5dх = х25 + х— | ¿Х =|~ + =

4 3

12

12

¿„ - = Г 17х + 5)5хd5dх = |"| х2--I--I ¿х =1--I--

2,2 [[ [I 23 1 I 6 6 ,0

0 0 1 1

Г1 = [ [ (х + 5~Щк (х)/ (5^х = [ [ (х + 5~Щк (х)5d5dх, к = 1,2 , 0 0 0 0

0 0 1 1

1 ( 5 2 5 3 I1

/I = [[(х + s)sdsdх = [| х— + — ¿х =

00 0 V 23 У0

х I1

V 4 + 3 ,0

7

12'

1 1

„2 „„Л1 Л.3 „л1

, 5 х5 '-+-

2 3

¿х =

хх

Т ~в

/I = [[(х+s)хsdsdх = х

00 0 V

11 _

/к = [ [ хщк (х) /2 (5)С5Сх = [ [ хщк (х)(25 — 2 / 3)С5Сх к = 1,2,

11

/2 = [ [ х5(25 — 2 / 3)Ых=

1/251 х 2хs^1

3 6

Сх=

,0

/,2 =[[ х25(25 — 2/3)с15ёх =[

°' 253 х2 2х2 5 2

2х х' ~ "б

' 2х3 х3

= 6 , 1

3

6

9 9 I 9

Используя найденные коэффициенты, составим СЛАУ (11):

п п _

X ак,Р) — Х,2 X Ьк,1С' = Х1,2/к , k, п = 1,2

1=1 п

=1 п

О

X ак, 1С 2 — А X Ск,1С1 = Х2,1/к , k, п = 1,2 =1

а00C0 + a02C0 -(Ь00C0 + Ь12С2 )= /12 1С1 + а22С1 — (Ь20C0 + Ь22С2 ) = ./2

с| + а2 (¿1+ с1)—/

О

С+° с° —|1 С2 +1С221 =

2 2 V 4 1 6

° С + ° С° —|1С2 +1С221= 1

2 1 3 2 V 6 1 9

9

-- С11=1

12 2 I 12

С2 + -С22 —| С1 + — С1 1 = —

О

(16)

1 с2+1С 22 — С+1 С° 1 =1

2 1 3 2 V12 1 3 21 3

2

С1 = 0, с1 = 1, С2 = 2, с\ = 1

подстановкой убеждаемся, что найденные коэффициенты

удовлетворяют

2

С1 = 0, с1 = 1, С2 = -, с\ = 1

3'

СЛАУ (16).

Согласно формуле (3), решение задачи (15) имеет вид

и(х) = / (х) + С р (х) + С\р2 (х) = х + х = 2х,

(17)

2 2 2 2 г(х) = /2 (х)+С р (х) + С2р (х) = 2х — + х = 3х

Подставляя численно найденные решения (17) в условие примера (15), получим, что

Г ,|1

2х — [ х535С5 = 2х — х5 | = х = / (х)

3х — [(х + s)2sds = 3х — | х52 + 2 53 | = 2х — — = /2 (х)

3

численное решение (17) является точным в примере (15).

Замечание 1. В некоторых частных случаях, используя дополнительную информацию о правых частях задачи вида (2) или (8), точное решение задачи вида (2) или (8) можно получить с меньшим набором координатных функций. Найдем численное решение примера (18), используя по одной координатной функции (поскольку правые части (18) содержат всего одну координатную функцию степенного вида р (х) = х, щ (х) = х . В этом случае матрицы А, В, Б и векторы Б1, Б2 являются числами.

и (х) — [ х5У(5)С5 = х,

0

1 , у( х) — [ хи(5)С5 = 2 х

0

А = А = 1, К г(х, 5) = х5, К21 (х, 5) = х , / (х) = х, /2 (х) = 2х, а = 0, Ь = 1.

(18)

Оц = [щ (хрр (х)Сх = [ х2 ¿х = —

0 0 3 0 1 1

Ь 1 = [ [ К 2 (х, 5)р (5)щ (х)С5Сх

0 0

1 1 53 = [[ х5 • 5 • хС5Сх=-

0 0 3

11 11 С1 =11К1 (х, 5)р (5)щ (х)С^Сх = [ [ х • 5 • хС5Сх =

1

6 ,

0

0

0 0

1 1

0

0

2

0

0

0

11

0

0 0

11

х

5

3

2

= Ц КхЛх ¿)У1(х)/2(!!)^х =

а а

1 1

= Ц (xs)x2sdsdx = 2 -

= Ц к2л( ^ ¡¡Уу&Ж-)^ =

Л х ■ х ■ ^^=||

, 2

о V 2 Ло

dx=

V 6 л,

Запишем СЛАУ, используя формулу (11),

Я2 =¿¡,1 = 1 .

¡апс; -ЪС = Х [аиС1 - dllCl = А

О <

1с1 -1С2 = -

3 1 9 1 9

1 с - -1 с ; = 1 3 1 6 1 6

о

с = 1, с- = 1.

По формуле (9) получим

|и( х) = / (х)+С\р (х) = х+1- х = 2х, [у(х) = /2 (х) + С\р (х) = 2х+1 ■ х = 3х

Сделаем проверку:

Г ,|1

2х-I х-3= 2х - х- = х = / (х)

о

3х -1 х2^- = 3х - (х-21 = 2х = /2 (х)

о

В примере (18) найденное численное решение и(х) = 2х,у(х) = 3х является точным. По сравнению с примером (15), в котором вычислено 16 интегральных коэффициентов, в примере (18), мы вычислили всего 5 интегральных коэффициентов, то есть в три раза меньше.

Замечание 2. Предположим, что точное решение задачи (2) имеет разложение по системе координатных функций в виде бесконечного ряда или с конечным числом слагаемых Ы>п, превышающим число координатных функций {р (х)}п=1; п < N .

Тогда точное решение имеет вид [1],

[2], [3]:

N _ _

у (х)ехс = А (х) + x(х), г = 1,т, j = 1,N (18)

j=1

или

j=1

У1 (х)пит = А (х) + X Cрj (х) г = 1 т j = 1 П

С учетом формулы (3) получим невязку для каждой функции решения задачи (2) соответственно в случаях (18) и (19):

N _

Г (х) = Уг (х)ехс - (х)пит = XС)р (х), г = 1 т , (20)

j=п+1

да _

Г (х) = Уг (х)ехс - У (х)пит =X^ г = 1 т . (21)

j=n^-1

Если системы линейно-независимых функций {р(х)}п=1 ,{у1 (х)}п=1 являются полиномами, то, как показано в работах [1], [2], [3] они выбираются, начиная с самых грубых функций (дающих максимальный вклад погрешности без их использования), например, в виде степенных мономов (п+1 функция в системе)

ро (х) = У о (х) = 1, р (х) = У1(х) = х,...,Рп (х) = у (х) = хп}.

Тогда невязки решения в формулах (20), (21) дают одинаковую асимптотику решения:

ехас( ' ишт\\

--О

Xср (х)

У^1

--О(хп+1 = 1,т. (22)

У(хХхас, = /(х) + Xср(х\ г = 1,т,] = 1,2,3,... (19)

j=l

Рассмотрим достаточные условия корректности алгоритма (13), (14) численного решения задачи (8).

Теорема 1. Достаточные условия корректности численного решения задачи (8). Пусть матрица A обратима (существуют матрица А-1обратная к А) и выполнено условие

А"ТIIБ11И - Ч < 1, тогда алгоритм (13), (14) корректен и верна оценка

||(а-ЛцЛЦОА-^Ч < 1Ц ч = \Я2\|Я2Д||А-1|2\\Б\Р\\ и н 1 - 4

||(А -ЛцЛцБА-'Б"^] <Ы . (23)

II II 1 - ч

Доказательство

Согласно формулам (13), (14) задачи (8)

С1 = (а-\2Л2Л БА-1вУ Я Р2 + АА1 БА~1Р1),

с2 = (а-я2ДА,2 °а _1б)-1(Яд р 1 +ЯЛ.2 °а 2)

векторы С1 С2 коэффициентов разложения функций решения и(х), у(х) по системе координатных функций {р (х)}'п=1 существуют тогда и только тогда, если существуют обратные матрицы А 1 и (а -Я12Я21БА1о) ,

(а-яА,2°а 1бУ .

о о

х

оо

и

Матрица А обратима по условию теоремы 1, что обеспечивает существование векторов

ЦF2 +ЦЦBA lF1) , (x2lF1 +ЦЦDA~lF2) в

формулах (13), (14) с конечной нормой

||Ц F2 +Ц2Ц BA^F '|| < Цр1+Ц ДЦ 2DA 4F 2|| < да.

Отметим, что условие

Ц2Щ||A"H IBUDI - q < 1 является общим для

матриц (a -Ц2ЦBA_1d)-1 , (a -ЦЦDA_1в) 1 ,

поэтому достаточно доказать обратимость одной из матриц.

Так как матрица Z представима в виде произведения двух матриц

Z - A - ЦЦBAlD - a([ - ЦЦ A"1BA"1D) , то если матрица Z обратима, получим Z 1 =(l-Ц2ЦA lBA lDjlA- . Здесь I - единичная квадратная матрица размерности (nx n). По условию Теоремы 1 выполнено условие

Ц 2Щд| |A-f|| B||||d|| - q < 1, следовательно,

(i - Ц2Ц1 a~1ba~1d )-1 = jj Цц) (A-1 BA-d)

k=0

Z 1 ) (A-1BA-1D)k 1 A-1.

матрицу Z :

IIZ "1 -

jj (4,24,1 )k (a-'ba-'d)

ЁЦЛ |4д| )k i|| A'1

Ail -

\a-1\\ -

qk A1 - 1

A-1

j-1

4

R(x) -j j (x)Ц K,,.(x, s)f (s)+ j j (s)

j-1 a V j-1

A

Требуем ортогональности невязок всем п линейно-независимым функциям второй системы {щ (х)}=1 < Я,. ,щ >= 0, г = 1, т, I = 1, п

п Ь Ь Ь ( п Л

X С [ р (х)щ (х)Сх—4 [[К / (х, (х) / (5)+X Ср (5) ¿С = 0,

1=1 а а а V 1=1

г = 1, т, 1,1 = 1, п

С учетом последней формулы и условия |ц||ц|ЦаЦ^!B||||d|| - q < 1 оценим по норме

ук=0 , ■■ 1—q Теорема 1 доказана.

Сравним условия корректности в диагональной задаче (24) из системы уравнений Фредгольма второго рода с условиями корректности численного алгоритма для одного интегрального уравнения (частный случай системы уравнений (24) при г=т=1):

Ь _

у (х) — Хг [ К . (х, 5)У1 (5)С5 = / (х), ' = 1, т . (24)

а

Система (24) имеет решение и невязку:

п ___

У1 (х) = / (х) + X Ср1 (х), ' = 1, т, 1 = 1, п ,

ds, i -1, m.

п Ь п ""

X С. [ р] (х)щ (х)Сх — 4 X С] [ [ К ; (х, 5)щ (х)р. (5)С5Сх =

1=1 а 1=1 а а

4 4 _ _

= 4[[ К ,(х, 5)щ (х)/ (5)С5Сх, г = 1, т, 1 = 1, п (25)

а а

Ь Ь Ь

а] ^ = [ р (х)щ (х)Сх — А [ [ К, (х, 5)щ (х)р (5)с5Сх

а а а

' = 1, т; 1,1 = 1, п

Ь 4 _ _

Л = ХЦК, (х, 5)щ (х)/ (5)С5Сх, г = 1, т, / = 1, п .

а а

Тогда формула (25) примет простой вид:

А'С' = Р' О С' = (а' Р',г = 1т . (26)

В формуле (26) Аг (п х п), г = 1, т матрица коэффициентов с номером г, С (п х1), г = 1, т -вектор-столбец размерности п с номером г. Рг, г = 1, т - вектор-столбец правой части в г-й строке системы интегральных уравнений (высоты п). Все указанные коэффициенты определяются формулой (25). Как показано в работах [1], [2], [3], для одного интегрального уравнения Фредгольма в системе уравнений (26), решение существует и единственно тогда и только тогда, если det(A) ф 0, соответственно

в (26) det(A' )ф 0, ' = 1, т , что дополнительно наклады вает ограничения. Ограничения на элементы матриц А', г = 1, т, в том числе на величину параметров ХИ и норму интегральных

ядер ЦК,,!, г = 1, т .

Разобьем матрицу а' (п х п), ' = 1, т на разность двух матриц и введем обозначения:

Ь ЬЬ

Ь\] = [р(х)щ(х)Сх = Ь ¿¡^ = [[Ки (х,5)щ(х)р(5)С5Сх,

а а а

а]. 1 = Ьи — Х,'<, О а' = В — Х,'Д,' = 1т . (27) Докажем Теорему 2.

Теорема 2 (достаточные условия корректности алгоритма (25), (26)). Пусть матрица В в (27) обратима и выполнены условия А В 11 Бг I < qi < 1,' = 1,т, mааxqi=q<1 ,

k-0

k -0

k

2

-

СО

тогда алгоритм (26) вычисления коэффициентов Сг, г = 1, т корректен. Верны оценки:

[ä)-'|| ,i = 1,m .

-q

(28)

Доказательство

По условию Теоремы 2 матрица B об-

ратима, тогда

A' = B -A'D = b(I -\ß lDt),i = 1

= 1, m .

Тогда формула (26) С' = (а')1 ,г = 1, т корректна, если и только если обратимы матрицы

а' = б(i -я,.,, б), г = 1т,

а' = г, = Б( - ЛцБ^Б' ) 2г 1 =(/-Я'В'О, )-1 Б1.

По условию Теоремы 2 |я.|||БЦДг]|<ч,-<1, имеем:

I1 = (I )-1 б-1 = ^б-'д ) \Бч, г = 1т.

Последнюю матрицу Zi 1 оценим по норме IK "1 ^Е^МИ DJ f B "I <ijr q,k Ii B II <

<IS* B-I=1

вi

, i = 1, m .

q

Теорема 2 доказана.

По сравнению с одним интегральным уравнением Фредгольма второго рода в диагональной задаче с системой интегральных урав-нений(24) (Теорема 2) требуется не только обратимость матрицы Б с одним условием ||б^Ц- ч < 1, но и т неравенств (условий-

ограничений U. J 5 1 Юг < qt,i = 1,m, max q, = q < 1).

Ill llll II i=1,m

В работе получены основные результаты:

1) Поставлена задача для системы интегральных уравнений Фредгольма второго рода с непрерывными ядрами на квадрате Кj(x,s) e Ci[a,b]x [a, b]), i , j = 1, m . Введено

определение диагональной задачи и антидиагональной задачи.

2) Впервые для задачи общего вида (2) предложен алгоритм ее численного решения обобщенным проекционным методом Петрова-Галеркина (3)-(7). Идея состоит в требовании ортогональности m

невязок интегральных уравнений системы к n функциям второй линейно-независимой системы (х)}.=1 . Предложен алгоритм решения антидиагональной системы двух уравнений Фредгольма второго рода (8)-(14). Формулы (13), (14) матричного вида для коэффициентов разложения C1,C2 е R" решения u(x),v(x) по первой системе {р (х)}".=1 функций. Указана формула для невязки неизвестных функций между точным и численным решением.

3) Используя алгоритм (8)-(14) численно решены два примера (15) и (17), в которых численные решения примеров совпали с точными решениями.

4) Доказаны Теоремы 1, 2 (достаточные условия корректности численных алгоритмов (13), (14) для антидиагональной задачи с двумя уравнениями Фредгольма и (25), (26) для диагональной задачи (24) общего вида).

Отметим также, что решение краевых задач в частных производных иногда сводится к решению интегральных уравнений [5], [6], [9]. Последние 10 лет активно изучаются уравнения с производными дробного порядка, сводимые к интегральным уравнениям [10].

Список источников

1. Бахвалов Н.С., Лапин А.В., Чижонков Е.В. Численные методы в задачах и упражнениях: учеб. пособие. М.: БИНОМ. Лаборатория знаний, 2010. 240 с.

2. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. 7-е изд. М.: БИНОМ. Лаборатория знаний, 2011. 636 с.

3. Бахвалов Н.С. Численные методы в задачах и упражнениях / Н.С. Бахвалов, А.В. Лапин, Е.В. Чижонков. М.: БИНОМ. Лаборатория знаний, 2013. 240 с. ISBN 978-59963-2266-4.

4. Васильева А.Б. Интегральные уравнения: учебник / А.Б. Васильева, Н.А. Тихонов; А.Б. Васильева, Н.А. Тихонов. Изд. 3-е, стер. СПб. [и др.]: Лань, 2009. 159 с. ISBN 978-5-8114-0911-2.

5. Полянский И.С., Логинов К.О. Приближенный метод решения задачи конформного отображения произвольного многоугольника на единичный круг // Вестник Удмуртского университета. Математика. Ме-

k=0

k=0

k=0

GO

ханика. Компьютерные науки. 2022. Т. 32, № 1. С. 107-129. DOI 10.35634/vm220108.

6. Юденков А.В., Володченков А.М. Устойчивость математических моделей основных задач анизотропной теории упругости // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. 2020. Т. 30, № 1. С. 112-124. DOI 10.35634/vm200108.

7. Волосова Н.К., Волосов К.А., Волосова А.К. [и др.] Решение интегральных уравнений Фредгольма методом замены интеграла квадратурой с двенадцатым порядком погрешности в матричном виде // Вестник Пермского университета. Математика. Механика. Информатика. 2022. Вып. 4(59). С. 9-17. DOI 10.17072/1993-0550-2022-4-9-17.

8. Ильюшин А.А. Труды. Теория термовязко-упругости. Т. 3. М.: ФИЗМАТЛИТ, 2007. 288 с.

9. Волосов К.А. Конструкция решений квазилинейных уравнений с частными производными // Сибирский журнал индустриальной математики. 2008. Т. 11, № 2(34). С. 29-39.

10. Кумыкова С.К., Эржибова Ф.А., Гучаева З.Х. Задача типа задачи Бицадзе-Самар-ского для уравнения смешанного типа // Современные наукоемкие технологии. 2016. № 9-1. С. 73-78.

References

1. Bahvalov N.S., Lapin A.V., CHizhonkov E.V. Chislennye metody v zadachah i uprazhneni-yah: ucheb. posobie. M.: BINOM. Labora-toriya znanij; 2010. 240 s. (In Russ.).

2. Bahvalov N.S., Zhidkov N.P., Kobel'kov G.M. Chislennye metody. 7-e izd. M.: BINOM. La-boratoriya znanij; 2011. 636 s. (In Russ.).

3. Bahvalov N.S. Chislennye metody v zadachah i uprazhneniyah / N.S. Bahvalov, A.V. Lapin, E.V. Chizhonkov. M.: BINOM. Laboratoriya

znanij, 2013; 240 s. ISBN 978-5-9963-22664. (In Russ.).

4. Vasil'eva A.B. Integral'nye uravneniya: uchebnik / A.B. Vasil'eva, N.A. Tihonov; A.B. Vasil'eva, N.A. Tihonov. Izd. 3-e, ster. SPb. [i dr.]: Lan'; 2009. 159 s. ISBN 978-58114-0911-2. (In Russ.).

5. Polyanskij I.S., Loginov K.O. Priblizhennyj metod resheniya zadachi konformnogo oto-brazheniya proizvol'nogo mnogougol'nika na edinichnyj krug. Vestnik Udmurtskogo uni-versiteta. Matematika. Mekhanika. Komp'yuternye nauki. 2022;32(1):107-129. (In Russ.). DOI 10.35634/vm220108.

6. Yudenkov A.V., Volodchenkov A.M. Ustojchivost' matematicheskih modelej os-novnyh zadach anizotropnoj teorii uprugosti. Vestnik Udmurtsko-go universiteta. Matemat-ika. Mekhanika. Komp'-yuternye nauki. 2020;30(1): 112-124. (In Russ.). DOI 10.35634/vm200108.

7. Volosova N.K., Volosov K.A., Volosova A.K., Karlov M.I., Pastuhov Yu.F. Solution of the Fredholm Integral Equations Method of Replacing the Integral by a Quadrature With the Twelveth Order of Error in Matrix Form. Bulletin of Perm University. Mathematics. Mechanics. Computer Science. 2022;4(59):9-17. (In Russ.). DOI: 10.17072/1993-0550-2022-49-17.

8. Il'yushin A.A. Trudy. Teoriya termovyazkouprugosti. T. 3. M.: FIZMAT-LIT; 2007. 288 s. (In Russ.).

9. Volosov K.A. Konstrukciya reshenij kvazi-linejnyh uravnenij s chastnymi proizvodnymi. Sibirskij zhurnal industrial'noj matematiki. 2008; 11(2(34)): 29-39. (In Russ.).

10. Kumykova S.K., Erzhibova F.A., Guchaeva Z.H. Zadacha tipa zadachi Bicadze - Samar-skogo dlya uravneniya smeshannogo tipa. Sovremennye naukoemkie tekhnologii. 2016;(9-1):73-78. (In Russ.).

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

Наталья Константиновна Волосова - аспирант МГТУ им. Н. Э. Баумана (105005, Россия, Москва, 2-я Бауманская ул., д. 5, стр. 1), [email protected], https://orcid.org/0000-0538-2445;

Константин Александрович Волосов - доктор физико-математических наук, профессор кафедры прикладной математики Российского Университета Транспорта (127994, ГСП-4, Россия, Москва, ул. Образцова, д. 9, стр. 9), [email protected], https://orcid.org/0000-0002-7955-0587, AuthorID 128228;

Александра Константиновна Волосова - кандидат физико-математических наук, начальник аналитического отдела ООО "Трамплин" Российского Университета Транспорта (127994, ГСП-4, Россия, Москва, ул. Образцова, д. 9, стр. 9), [email protected], https://orcid.org/0000-0002-0538-2445, AuthorID 607500;

Михаил Иванович Карлов - кандидат физико-математических наук, доцент кафедры высшей математики МФТИ (141701, Россия, Московская область, г. Долгопрудный, Институтский пер., 9.), [email protected], AuthorlD 14680;

Дмитрий Феликсович Пастухов - кандидат физико-математических наук, доцент кафедры технологий программирования Полоцкого государственного университета (211440, Республика Беларусь, Витебская обл., г. Новополоцк, ул. Блохина, 29), [email protected], https://orcid.org/0000-0003-1398-6238, AuthorlD 405101;

Юрий Феликсович Пастухов - кандидат физико-математических наук, доцент кафедры технологий программирования Полоцкого государственного университета (211440, Республика Беларусь, Витебская обл., г. Новополоцк, ул. Блохина, 29), [email protected], https://orcid.org/0000-0001-8548-6959, AuthorlD 405109.

Information about the authors:

Natalya K. Volosova - Post-graduate Student of Moscow State Technical University. N. E. Bauman (5-1, 2nd Baumanskaya Street, Moscow, Russia, 105005), [email protected], https://orcid.org/0000-0538-2445;

Konstantin A. Volosov - Doctor of Physical and Mathematical Sciences, Professor of the Department of Applied Mathematics of the Russian University of Transport (9-9, Obraztsova Street, Moscow, GSP-4, Russia, 127994), [email protected], https://orcid.org/0000-0002-7955-0587, AuthorlD 128228;

Aleksandra K. Volosova - Candidate of Physical and Mathematical Sciences, Chief Analytical Department "Tramplin" LLC Russian University of Transport (9-9, Obraztsova Street, Moscow, GSP-4, Russia, 127994), [email protected], https://orcid.org/0000-0002-0538-2445, AuthorlD 607500;

Mikhail I. Karlov - Candidate of Physical and Mathematical Sciences, Associate Professor of the Department of Higher Mathematics, Moscow Institute of Physics and Technology (9, Institutskiy per., Dolgoprudny, Moscow region, Russia, 141701), [email protected], AuthorlD 14680;

Dmitriy F. Pastukhov - Candidate of Physical and Mathematical Sciences, Associate Professor of Polotsk State University (29, Blokhin Street, Novopolotsk, Vitebsk Region, Republic of Belarus, 211440), [email protected], https://orcid.org/0000-0003-1398-6238; AuthorlD 405101;

Yuriy F. Pastukhov - Candidate of Physical and Mathematical Sciences, Associate Professor of Polotsk State University (29, Blokhin Street, Novopolotsk, Vitebsk Region, Republic of Belarus, 211440), [email protected], https://orcid.org/0000-0001-8548-6959, AuthorlD 405109.