CC BY
16ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА 2022 • Математика. Механика. Информатика • Вып. 4(59)
МАТЕМАТИКА
Научная статья УДК 519.6
DOI: 10.17072/1993-0550-2022-4-9-17
Решение интегральных уравнений Фредгольма методом замены интеграла квадратурой с двенадцатым порядком погрешности в матричном виде
Н. К. Волосова1, К. А. Волосов2, А. К. Волосова2, М. И. Карлов3, Д. Ф. Пастухов4' 5, Ю. Ф. Пастухов4' 5
Московский государственный технический университет МГТУ им. Н.Э. Баумана, Москва, Россия 2Московский Университет Транспорта (МИИТ), Москва, Россия 3Московский физико-технический университет МФТИ, Москва, Россия 4Полоцкий государственный университет, Новополоцк, Белоруссия ^Московский государственный университет им. М.В. Ломоносова, Москва, Россия
Автор, ответственный за переписку: Дмитрий Феликсович Пастухов, dmitrij.pastuhov@mail.ru
Аннотация. Предложен алгоритм численного решения уравнения Фредгольма второго рода с непрерывным ядром методом замены интеграла и матричным решением СЛАУ с квадратурной формулой двенадцатого порядка погрешности с числом интервалов интегрирования, кратным десяти. Новая формула, по сравнению с формулой Симпсона, дает 15 значащих цифр для узловых значений функции решения даже при небольшом числе интервалов 10,20 на отрезке за конечное число элементарных операций. Полученный алгоритм имеет двойную точность и минимальное время вычислений. В то время как формула Симпсона совместно с матричным методом решения СЛАУ дает только 6 значащих цифр с числом интервалов интегрирования равным двадцати. Более того, для формулы Симп-сона двойная точность недоступна (15 нулей в бесконечной норме невязки решения), так как язык FORTRAN допускает максимальные массивы матриц 200*200. Получены оценки верхней границы допустимого параметра |Х| для матрицы уравнения Фредгольма со строгим диагональным преобладанием или с небольшой нормой интегрального ядра. Ключевые слова: уравнение Фредгольма; численные методы; уравнения математической физики; матрица; интегральные уравнения
Для цитирования: Волосова Н. К., Волосов К. А., Волосова А. К., Карлов М. И., Пастухов Д. Ф., Пастухов Ю. Ф. Решение интегральный уравнений Фредгольма методом замены интеграла квадратурой с двенадцатым порядком погрешности в матричном виде // Вестник Пермского университета. Математика. Механика. Информатика. 2022. Вып. 4(59). С. 9-17. DOI: 10.17072/1993-0550-2022-4-9-17.
Статья поступила в редакцию 26.08.2022; одобрена после рецензирования 20.10.2022; принята к публикации 31.10.2022.
MATHEMATICS
Research article
Solution of the Fredholm Integral Equations Method of Replacing the Integral by a Quadrature with the Twelveth Order of Error
in Matrix Form
N.K. Volosova1' K.A. Volosov2, A. K. Volosova2, M. I. Karlov3, D. F. Pastuhov4'5' Yu. F. Pastuhov4'5
University of the Moscow State Technical University named after N.E. Bauman, Moscow, Russia
Это работа О 2022 Волосова Н. К., Волосов К. А., Волосова А. К, Карлов М. И., Пастухов Д. Ф.. Пастухов Ю. Ф. лицензируется под CC BY 4.0. Чтобы просмотреть копию этой лицензии, посетите http://creativecommons.org/licenses/by/4.0/
2Russian University of Transport (MIIT), Moscow, Russia
3Moscow Physics and Technology University MIPT, Moscow, Russia
4Polotsk State University, Novopolotsk, Belrussia
5Moscow State University named after M.V. Lomonosov, Moscow, Russia
Corresponding author: Dmitriy F. Pastukhov, dmitrij.pastuhov@mail.ru
Abstract. An algorithm for the numerical solution of the Fredholm equation of the second kind with a continuous kernel by the method of integral replacement and the matrix solution of SLAE with a quadrature formula of the twelfth order of error with a number of integration intervals divisible by ten is proposed. The new formula, compared to Simpson's formula, gives 15 significant digits for the nodal values of the solution function, even with a small number of intervals of 10.20 on a segment in a finite number of elementary operations. The resulting algorithm has double precision and minimal computation time. While Simpson's formula, together with the matrix method for solving SLAE, gives only 6 significant digits with twenty integration intervals. Moreover, double precision is not available for Simpson's formula (15 zeros in the infinite norm of the solution residual), since the FORTRAN language allows maximum matrix arrays of 200*200. Estimates are obtained for the upper bound of the admissible parameter |X| for the matrix of the Fredholm equation with strict diagonal dominance or with a small norm of the integral kernel. 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 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-4-9-17.
The article was submitted 26.08.2022; approved after reviewing 20.10.2022; accepted for publication 31.10.2022
Введение
В работе [1] описан метод замены интеграла для численного решения уравнения Фредгольма второго рода. Для этого нужно составить систему алгебраических линейных уравнений (СЛАУ), в которой неизвестными являются узловые значения функции. Примеры из задачника [1] используют квадратурную формулу Симпсона в методе замены ядра.
В данной работе мы предлагаем матричный алгоритм численного решения СЛАУ, используя квадратурную интегральную формулу с 12-м порядком погрешности [2], [3].
В работе [2] использован метод последовательных итераций, метод прост, но ограничивает допустимую область по параметру X, то есть применим при малых по модулю значениях X. Поэтому в данной работе матричный алгоритм решения СЛАУ существенно расширит область допустимого параметра X.
Постановка задачи
Рассмотрим уравнение Фредгольма 2-го рода (1), которое нужно решить методом замены интеграла, т. е. таблично на дискретной сетке х,г = 0, п найти узловые значения неизвестной функции у, г = 0, п, .
Метод описан в работе [1]:
y( x) K (x, s)y(s)ds = f (x)
(1)
где: К(х,з) е С([а,Ь]х[а,Ь]) - непрерывная функция (кусочно-непрерывная на квадрате (х, я) е [а,Ь] х [а,Ь] - ядро интегрального уравнения Фредгольма(1)), /(х) е С[а,Ь] - заданная непрерывная правая часть уравнения (1), X - известный параметр. Заменяя интеграл квадратурной интегральной формулой [1], запишем дискретный аналог (1):
п , > п
у1 С К (хг, у = / (х )Уг = Л ,
j=o
j=o
xi = a + ih, s. = a + jh; i, j = 0, n .
(2)
В формуле (2) - алгебраической системе линейных уравнений относительно у весовые
коэффициенты ,у = 0, п заданы, с. определяют точность решения системы (2) и неявно
содержат шаг интегрирования h =
b - a
n -
число интервалов интегрирования на отрезке №
В работе [2] мы использовали аналогичную формулу (3), которая содержит эффективный шаг интегрирования # в явном виде:
г п ^
y, -Wp
V j=
X yjKuCj = fi, y, = у(х, I
a
n
Ки = к(х,з}} = 0,п,х = а + И• г,И = —а,у: = /(х) , (3)
п
рИ , Ь - а
где: Н =,И =- - эквивалентный шаг
р 2 п
интегрирования, индекс кратности р связан с целыми числами р, к п = рк, р, к е N , число интервалов п нацело делится на р.
В работе [2] нами построена и использована в работе [3] квадратурная интегральная формула с двенадцатым порядком погрешности (р=10). Весовые коэффициенты в системе (3) определяются алгоритмом (4):
ph 10 h
H„==w =^- = — = 5Кh =
b — a
2 2
C =
r 16067 ■ n ■
C =-, / = 0 v / = n;
299376
Л/ГП г
C =-, /(mod 10) = 1 v /(mod 10) = 9;
74844
C3 = —
16175
99792
, /(mod 10) = 2 v /(mod 10) = 8;
C =
5675
6237
/(mod 10) = 3 v /(mod 10) = 7;
C = —
4825 5544
, /(mod 10) = 4 v /(mod 10) = 6;
У, — 5hk
X.У/K,./C/ + 0(h12)= f,i,/ = 0,n, (5)
/=1 )
К1,1 = К(х;,^)-Х = а + И• ^^ = а + И• И = —,Л = /(хг) •
п
В формуле (5) весовые коэффициенты С.
записаны в виде рациональных дробей согласно формуле (4) с двойной точностью. Определим квадратную диагональную матрицу:
[C,, / = к D . = 8 ,C, = ] /
к,/ /
8 , =
0, / * к' к J [0, / * к'
1, / = к
X D ,/ = X8,/C/ = Cj, D / = D
к=0 к=0
J,к ■
(6)
У — 5hA
( n
X /y,
V/-0
( n (
\
(
nn
о y — 5hl
X X K,jD_
--fi о y, — 5M
Л ^ Ук
X Ku X
V/'=° к=0
,к
Vк=0V/
= f, i = 0, n . (7)
(4)
Формула (7) c учетом порядка аппроксимации интеграла o(h12) (5) и весовых коэффициентов (4), эквивалентна матричной записи (8): (I — 5hX • KD)y + o(h12) = f о y = (I — 5hA • KD)— f + o(h12) . (8)
Возьмем от обеих частей уравнения (8) функцию нормы:
y = (I — 5h2-KD)—1 f ^
IHL=|(I — 5h^KDffl <||(I — 5h^KD)—^||f|L . (9)
Определение. Назовем матрицу A = I—5hl • KD матрицей уравнения Фредгольма 2-го рода с весами (4).
В формуле (9) выбраны бесконечные матричная и векторные нормы, так как они являются согласованными [1]. Матричная и векторная норма определяются формулами
n
=max XI / y\\„=y
» = () П . 1 1 / = 0,n1
= 1Z807, /(mod 10) = 5;
6 12474
= J606^,/(mod 10) = 0л0< /<n;
7 149688
10
X C/ = 2
/=0
Окончательно, с учетом (4) [3, стр. 197], п = 10к, к е N формула (3) переходит в формулу (5):
( п \
/=0
(10)
Из формул (10) следует, что для непрерывных функций, которые дискретно проектируются на двухмерную сетку (квадрат) и одномерную сетку (отрезок), бесконечные матричная и векторная нормы всегда существуют и ограничены. Из неравенства (9) следует, что бесконечная норма вектора решения у ограничена в случае
||/||м <+<» и ||(/ - 5НХ • КО)-1|| <+да .
Поэтому вектор решения непрерывен и ограничен, если ограничена норма
I - 5И1 • К^)-1 <+< . Обозначим
A = I — 5h2-KD, r = min
* i=0,n
I— X
I
/=0,n, /*i )
Теорема 1 (оценка верхней границы нормы обратной матрицы) [4].
Для матрицы А с диагональным преобладанием справедлива оценка
llA-1 < I =
1
min
i=0,n
I— X
л •
(11)
/=0,n, /*i )
n
0
n
a
x
a
Теорема 2 (оценка верхней границы модуля параметра |X|).
Пусть матрица A = I — 5hl- KD имеет строгое диагональное преобладание и пусть
(b - a) min I max K(x, s)| 1 —10h min IK(x = s)|min |C| > 0,
xe[a,i]V se[a,4]' ') xe[a,4]' 1 i=0,n' ''
тогда для допустимого параметра X справедлива оценка
0 <ll<7-7-1-
I min I max\K(x,s) 1—10h min \K(x = s)|min|C.|
V*e[a,b]V sE[a,bJ ') xe[a,bj 1 ¡=0,и' ''
(12)
Доказательство. По Теореме 1 норма обратной матрицы lA"1]! = ||(/ — 5hl- KD)-^
ограничена сверху, тогда верно
f „ Л
> 0.
r = min
— S к
V j=0,n, j*i
kl — S kl-Г = min
j=0,n, j*i
\
> 0Vi = 0, n о
Ю
kl— S kl
j=0,n, (
n
|1 — 5hlKiJC\ — 5h|l SIKi,j||Cj| > 0.
j=0,j *i
1 — 5hlKi f\ + 5h|l— 5h|ISjKi,j||Cj| > 0,
n n _
|1—5hlK,C,| — 5h|l SKIICj = Kl — S kl - Г > 0Vi = 0,n
j=0,j*i j=0,n,
n
1 +10 h\l Ki j ||Ci I — 5h| lS| Ki, j ||Cj I -
j=0
n
|1 — 5hAKiiC\ + 5h|l\K,\C\ — 5h|lS|Ki,j||Cj| > 0
j=0
5h| lSKlkl < 5 ьa) и maxi K, iSki=
j=0
j=0
5 || maxiK ,|2k = i0^ (й - a)|Ы maxK I = (й - a)|| maxK I.
И j=0,J ' 1 10к У=0,л' ' ' Мл' ' '
£сли 1 +10кЫКи\\С} -(й -а)|Ы maxiK I > 0 ^
И _
1+10k|||K .||С|. - 5к||^Ки||С| > 0у/ = 0, и (13)
' j=0
Очевидно, что формула (13) будет выполнена при любых значениях х е [а, Ъ], i = 0, и , если она будет выполнена при условии
1 + 10к\Ыmin K,.|min\С\ -(й -а)|Ыmax K . > 0 ^
i=0,n 1 i=0,n
j=0, n
1+10h| l min K ( x, s = x) min CI—(b — a) l max IK (x, s) > 0
1 se[a,b] 1 i=0,n И 1 se[a,b] 1
Последнее неравенство написано также для непрерывного ядра K(x, s).
1
0 <l<
(b — a)max|K(x,s) — 10h min |K(x = s)|min|C.|
se[a,b] 1 se[a,b] 1 i=0,n' i1
min C I = « 0,053668 .
Таким образом, для каждой строки матрицы А = I - 5И1- КО справедливо неравенство (учтем также (4) и п = 10к, к е N)
, (14)
i=0,»' ^ 299376 Далее, потребуем выполнимость (14) при любых значениях переменной х из отрезка [a,b], тогда в задаче всегда выполнено неравенство (13). Потребуем, чтобы знаменатель (14) был положительным при всех параметрах х, имеем, если
(b-a) min I max|K(x,s)| |-10h min |K(x = s)|min|C.| > 0 ^
xe[a,b\ se[a,b] 'J xe[a,b]' 1 i=0,n'
Vx e [a, b]: (b - a) max K(x, s)| -10h min \K(x = s)| min \C, | > 0
se[a,b] 1 xe[a,b] 1 i=0,n
0<W<, 1
min I maxK(x,s)| 1-10h min |K(x = s)|min|C| I
xe[a,b]^ se[a,b] 'J xe[a,bf 1 i=0,n' 'J
(15)
Теорема 2 доказана, так как формулы (15) и (12) совпадают.
Теорема 3. Пусть выполнено условие 1 n
5(b-aillmax — Kj|max\C\ < 1 , тогда обрат-
i=0,n n - , '=° "
j=0
ная матрица к матрице А = I - 5И1- КО существует, алгоритм (4), (5), (8) корректен и верно неравенство
1
5(b — a)max|C, I max1S K
v ' i=07n^ Л i=07n n S '
■ о
1 b
5(b—a jb—a) ilK ( x,s) ds
. (16)
Доказательство.
Известно, что если |< 1, то
1 ад
(I-5)-1 = I + В + В2 +... = I+ УВк .
У ' I-В Ы
В качестве матричной нормы рассмотрим бесконечную норму
В = 5АЯ-Щ Щ 5к\Щ К Л О| „< 1«
5hl max SKj I maxi C , 1 о 5 11 max ¡TK Imaxl C < 1
i=0,n , ' i=0,n у i=l» , 'i=l»
!>—a) i
n 1 i=0,n
<1о
j=0
5(b — a)lmax-¿K,, Jmax|C\ < 1 о
i=0,n n j-Q i=0,n
1 1
5(b — a)max| Cj max1 ¿Ik,. J 5(b — a)max| C| max -,—г il K (x, s)\ds
i=0,n i=0,n n j=0 1 i=0n' x=[a,b] (b — a)a
i=0,n
n
n
1
В формуле (16) интеграл
1 ь
-Л К (х, 5)1 ds представляет собой среднее
(Ь - а) а
значение модуля интегрального ядра по переменной 5. Теорема 3 доказана.
Теорема 4. Пусть выполнено условие в котором все нормы матриц - нормы Фробе-ниуса < 5к\Ц|ЦоЦ^ < 1, тогда обратная
матрица к матрице А = I - 5И1 КО существует, алгоритм (4), (5), (8) корректен и верно неравенство
11<-, * . (17)
Однако шесть примеров, приведенных ниже, решенных численно алгоритмом (4), (5), (8), не имеют диагонального преобладания и не являются монотонными или М-матрицами, а Теорема 1 применима только для одного примера.
Рассмотрим примеры из замечательного задачника по численным методам [1].
1. Решить интегральное уравнение (№ 29.3) методом замены ядра [1, стр. 137]
11
y(x) — - J xesy(s)ds = e
3,714493834 (b — a)JJ K 2(x, s)dsdx
I a a
Доказательство. Рассмотрим матрич- нием y(x) = x + ex.
с K(x,s) = xes,f(x) = ex= точным реше-
ную норму Фробениуса || A^ = X
a,.
i=0J=0
n
B = 5hb KD, (KD)=X K.fj
/=0
|(KD)J <Ц" <Щ^/ JX/KDl = /Xl(KD)il IBIIf = 5h| 411 Щ\г = {a- ]hii щ\г <
5(b — a)H 1XC41XXXx|K',/12 = 5(b — a)|¿1 к XC
|| n /=0 H n i=0 /=0
Так как n = 10 к, к = n/10
n /=0
n I 1 n
X и xk
0 V n /=0
к 10 1 10
XC/2 =J^ XC/4BIL = 5h|411 KDf <
/=0 К 10 /=0
25.
10 Pa /
(b - a)|4 JX C
и<
n 1 n
XI XIK,
i=0 V n /=0
Л
(b—a)fFJ X [ $ )
1
< 1 о
■ о
(b—a)
i5 X C.2 10 ii 1
b—a
IK !(x, s)dsdx
a 11 ьь
—(b—a)XC! J J K !(x, s)dsdx
Программа написана на языке Compaq Visual Fortran 6.6 по алгоритму (4), (5), (8) (код программы приведен в конце работы) с условием примера (№ 29.3 [4]), результаты программы представлены в таблице. Обратная матрица вызывается библиотекой msimsl.
Все переменные и функции в программе заданы с двойной точностью.
Значение узлов равномерной сетки, численные и точные узловые значения функции, модуль их разности с параметром n=20, для примера 1- № 29.3 [1]
3,714493834 (Ь - а)| |К 2(х,
I а а
Теорема 4 доказана. В работе [4] рассматриваются монотонные матрицы А (матрицы, у которых элементы обратной матрицы по-ложительныА-1>0). Там же рассматриваются квадратные М-матрицы (матрицы у которых элементы главной диагонали положительные, а остальные элементы отрицательные).
xi num yi exact yi \ = num exact \yi — yi
0.10000 00 1.0048374180 3596 1.0048374180 3596 0.0000000E +00
0.20000 00 1.0187307530 7798 1.0187307530 7798 2.2204460E -16
0.30000 00 1.0408182206 8172 1.0408182206 8172 0.0000000E +00
0.40000 00 1.0703200460 3564 1.0703200460 3564 2.2204460E -16
0.50000 00 1.1065306597 1263 1.1065306597 1263 2.2204460E -16
0.60000 00 1.1488116360 9403 1.1488116360 9403 8.8817842E -16
0.70000 00 1.1965853037 9141 1.1965853037 9141 2.2204460E -16
0.75000 00 1.2223665527 4102 1.2223665527 4101 1.1102230E -15
0.80000 00 1.2493289641 1722 1.2493289641 1722 4.4408921E -16
0.90000 00 1.3065696597 4060 1.3065696597 4060 6.6613381E -16
1.00000 0 1.3678794411 7144 1.3678794411 7144 8.8817842E -16
Norma C 1.110223024625157E-015, time resolving 4.687500000000000E-002 ms
x
0
2
0
2
1
о
bb
1
<
В примере 1 программа по алгоритму (4), (5), (8) дает значение невязки по норме Чебышева (Norma C):
1 _
<||A~lf =|(/ - 5hl- KD^I < 1 1
y
-y
= 1.11022302 4625157 -10"
min
i=0,n
- X h
1/6
= 6 < +Ю .
Если в примере 1 использовать весовые коэффициенты из формулы Симпсона, то решение матричным методом (4), (5), (8) дает норму невязки
j=0,n, j*i )
Исследуем эту же задачу с помощью Теоремы 2.
У
-У
= 2.73255722 0497256E -10"
min \C\ = 16067/2993 76 и 0,053668
i=0,n
(b-a) min I maxiK(x,s)| 1 — 10h min K(x = s)|mjn|Cj > 0 о
xe[a,b]yie[a,b] ) xe[a,b] i=0,n
что в 108 раз больше чем при использовании формулы (4). Теорема 1 не применима, так как
матриЦа A = I— 5hl•KD не имеет диагонального (1 — 0) min Сmax|1 + 2xs|] —10^ mn J1 + 2x2\ milLCI > 0
преобладания. С учетом Теоремы 4 обоснуем xe[0-1]'v se[0-1] применимость алгоритма (4), (5), (8).
20 xe[0,1]l
1/2 = ll<-
3,714493834 (b - a) JJk2(x,s)dsdx
1
1
S 0,639 •
y( x) — J (1 + 2xs)y(s)ds = — x/6 — 1/2
K (x, s) = 1 + 2xs, f (x) = — x /6 — 1 / 2, l = 1 с точным решением y(x) = x +1/2 .
Программа дает норму невязки
I, ,num , ,exac. y — y
= 2.44249065 4175344 -10"
x +
1/2.
3/^9
||— x/6 — 1/2^ 2/3 4
(1 — 0) min (1 + 2x) — 10 —H-1606^ > 0 о
xe[0,1] ' 20 299376
1
1--0,053668 > 0
2
условие выполнено, тогда верно
0 <|l = 1 < 1
I , ,, 3,714493834 л/(в2 -1)/6
3,714493834 (1 - 0)Цx1elsdsdx ^ '
V 0 0
Замечание 1. Программа показывает, что диагональное преобладание в каждой строке матрицы А = I - 5И2 КО не зависит от вида весовых коэффициентов в формуле (3), так как в каждой строке совпадают первые 6 цифр диагональных преобладаний как при вычислении с весовыми коэффициентами по формуле (4), так и при использовании формулы Симпсона с эффективным шагом:
рИ 2И Ь- а
р=2 = т = Т = = п '
2. Решить интегральное уравнение (№29.12) [1, стр. 139], применим алгоритм (4),
(5), (8).
min I max K(x, s) 1 —10h min K(x = s)min |C.|
xe[a,b]( se[a,b] ') ie[a,bf 'i=0,^ ''
-о
(1 — 0) min max|1 + 2xs| —101—0 min 1 + 2x21 min \C,11
xe[0,1] se[0,1] 1 20 *e[0,1]l li=0,n=20' i')
■ о
о
1 —10,053668
2
i 1.0276.
в задаче 2 = 1, поэтому Теорема 2 здесь также применима.
3. Решить интегральное уравнение (№30.1) [1, стр. 141], применим алгоритм (4),
(5), (8).
y(x) + J x(exs — 1)y(s)ds = es — 1
с числом интервалов интегрирования п=20. Матрица А = I - 5И2 КО имеет диагональное преобладание, по Теореме 1 получим
К(х, 5) = х(вх5 -1), /(х) = вх -1,2 = -1 с точным решением у( х) = 1.
Программа дает значение нормы невязки \ушт -уехас,\с = 4.44089209 8500626 • 10-16 с
числом интервалов интегрирования п=20. Теорема 1 не применима, так как матрица А = I - 5А2- КО не имеет диагонального преобладания.
4. Решить интегральное уравнение (№ 30.1) [1, стр. 141], применим алгоритм (4), (5), (8):
1
у(х) - | = х ^ К(х, 5) = х5, /(х) = х, 2 = 1
-1
с точным решением у(х) = 3х .
C
C
1
a a
1
1
0
C
"Л
Программа дает значение нормы невяз-
ки
У
-y
= 8.88178419 7001252 -10"
числом интервалов интегрирования п=20. Теорема 1 не применима, так как матрица А = I - 5Л2- КО не имеет диагонального преобладания.
5. Решить интегральное уравнение (№30.2) [1, стр. 141], применим алгоритм (4), (5), (8).
y( x) -1 (xs + x2 )y(s)ds = 1
K(x, s) = xs + x2, f (x) = 1,2 = 1 с точным решением y( x) = 1 + 6x2.
Программа дает значение нормы невяз-
ки
У
-y
= 3.10862446 8950438 -10"
числом интервалов интегрирования п=20. Теорема 1 не применима, так как матрица А = I - 5Л2- КО не имеет диагонального преобладания.
6. Решить интегральное уравнение (№30.8) [1, стр. 143], применим алгоритм (4), (5), (8).
y(x) - j(xs2 - x)y(s)ds = 1 + (4/3)x
K(x, s) = xs2 - x, f (x) = 1 + (4 / 3)x, 2 = 1 с точным решением y( x) = 1.
Программа дает значения нормы невяз-
ки y
,num „ exac.
-y
= 6.66133814 7750939 -10"
числом интервалов интегрирования п=20. Теорема 1 не применима, так как матрица А = I - 5Л2- КО не имеет диагонального преобладания.
Предположение: Если ослабить условия Теоремы 1, то есть попробовать доказать ее с условием
1
f Л
n
min s ai,j
i=0,n V j=0,n У
f Л
n
min S ai, j
i=0,n V j=0,n y
где
S ai,j
j=0,n
> 0, atJ > 0Vi = 0,n atJ < 0Vi = 0,n
j=0,n ' '
V j=°.n
- условное диагональное преобладание элементов >й строки матрицы А = I -5Н1-КО, т.е.
п
сумма а1,] сохраняет знак для всех строк
J=0П '
матрицы A.
Именно условное диагональное преобладание показывает программа в численном эксперименте во всех 6 приведенных примерах матрицы Фредгольма А = I - 5Л2- КО и устойчивость численного алгоритма (4), (5), (8).
В работе получены основные результаты:
1). Предложен алгоритм решения интегрального уравнения Фредгольма второго рода матричным методом с двойной точностью (4), (5), (8) для небольшого числа интервалов интегрирования n=10,20,... Весовые коэффициенты (4) алгоритма имеют существенное преимущество по точности по сравнению с применением весовых коэффициентов из формулы Симпсона.
2). Численно показано, что матричный метод (4), (5), (8) решения СЛАУ для уравнения Фредгольма второго рода имеет большую область допустимого параметра X, по сравнению с решением уравнения методом последовательных приближений.
3). Алгоритм (4), (5), (8) таблично дает ответы ко всем примерам для уравнений Фредгольма второго рода из задачника [4]. Что указывает на его практическую применимость в широкой области параметра 2 е (0,2)с R и класса непрерывных интегральных ядер K (x, s) е C([a, b]x[a, й]) . В теоремах 2-4 получены оценки верхней границы допустимого параметра |X| для численного алгоритма (4), (5), (8).
Данная задача с алгоритмом решения (4), (5), (8), как и задача [5], решены численно с двойной точностью. Отметим также работу исследователей А.Д. Чернышева, В.В. Горяй-нова, С.В. Кузнецова, О.Ю. Никифоровой [7].
Список источников
1. Бахвалов Н.С., Лапин А.В., Чижонков Е.В. Численные методы в задачах и упражнениях: учеб. пособие. М.: БИНОМ. Лаборатория знаний, 2010. 240 с.
2. Волосова Н.К., Волосов К.А., Волосова А.К. [и др.]. Решение интегральных уравнений Фредгольма с невырожденными ядрами последовательными приближениями квадратурой с десятым порядком погрешности // Тенденции развития науки и образования. 2022. № 85-2. С. 21-25. DOI 10/18411/tmio-05-2022-55/-EDN:CKXBNI.
3. Численные методы. Лекции. Численный практикум: учеб. пособие к лекционным и практическим занятиям для студентов специальности 1-400101 "Программное обеспечение информационных технологий", 1980101 "Компьютерная безопасность" / Пастухов Д.Ф., Пастухов Ю.Ф., Волосова Н.К., Волосова К.А., Волосова А.К. Ново-полоцк. М., 2021.
с
с
-1
c
с
с
с
r
n
n
n
4. Волков Ю.С., Мирошниченко В.Л. Оценки норм матриц, обратных к матрицам монотонного вида и вполне отрицательным матрицам // Сибирский математический журнал. 2009. Т. 50, № 6. С. 1248-1254.
5. Пастухов Д.Ф., Волосова Н.К., Волосова А.К. Некоторые методы передачи QR-кода в стеганографии // Мир транспорта. 2019. Т. 17, № 3(82). С. 16-39.
6. Чернышев А.Д., Горяйнов В.В., Кузнецов С.В., Никифорова О.Ю. Применение быстрых разложений для построения точных решений задачи о прогибе прямоугольной мембраны под действием переменной нагрузки // Вестник Томского государственного университета. Математика и механика. 2021. № 70. С. 127-142. DOI 10.17223/19988621/70/11.
References
1. Bakhvalov N.S., Lapin A.V., Chizhonkov E.V. Numerical methods in tasks and exercises: a tutorial. M.: BINOM. Knowledge Laboratory; 2010. 240 p. (In Russ.).
2. N. K. Volosova, K. A. Volosov, A. K. Voloso-va, eds. Solution of Fredholm integral equations with non-degenerate kernels by successive quadrature approximations with the tenth
Информация об авторах:
order of error // Trends in the development of science and education. 2022;(85-2):21-25. DOI 10/18411/trnio-05-2022-55/-EDN:CKXBNI. (In Russ).
3. Pastukhov D. F., Pastukhov Yu. F., Volosova N. K., Volosova K. A., Volosova A. K. Numerical methods. Lectures. Numerical practice. Textbook for lectures and practical exercises for students of the specialty 1-400101 "Information technology software" 1-980101 "Computer security". Novopolotsk, M., 2021. (In Russ.).
4. Volkov Yu.S., Miroshnichenko V.L. Estimates of the norms of matrices inverse to monotone matrices and completely negative matrices. Siberian Mathematical Journal. 2009;(50:6): 1248-1254. (In Russ.).
5. Pastukhov D.F., Volosova N.K., Volosova A.K. Some methods of transmitting a QR code in steganography. World of transport. 2019;(17: 3(82)):16-398. (In Russ.).
6. Chernyshev A.D., Goryainov V.V., Kuznetsov S.V., Nikiforova O.Yu. Application of fast expansions for constructing exact solutions to the problem of the deflection of a rectangular membrane under the action of a variable load. Bulletin of the Tomsk State University. Mathematics and mechanics. 2021;(70):127-142. DOI 10.17223/19988621/70/11. (In Russ.).
Наталья Константиновна Волосова - аспирант МГТУ им. Н. Э. Баумана (105005, Россия, Москва, 2-я Бауманская ул., д. 5, стр. 1), navalosova@yandex.ru, https://orcid.org/0000-0538-2445;
Константин Александрович Волосов - доктор физико-математических наук, профессор кафедры прикладной математики Университета Транспорта России (127994, ГСП-4, Россия, Москва, ул. Образцова, д. 9, стр. 9), konstantinvolosov@yandex.ru, https://orcid.org/0000-0002-7955-0587, Аи^ГО 128228;
Александра Константиновна Волосова - кандидат физико-математических наук, начальник аналитического отдела ООО "Трамплин" Университета Транспорта России (127994, ГСП-4, Россия, Москва, ул. Образцова, д. 9, стр. 9), alya01@yandex.ru, https://orcid.org/0000-0002-0538-2445, Аи^ГО 607500;
Михаил Иванович Карлов - кандидат физико-математических наук, доцент кафедры высшей математики МФТИ (141701, Россия, Московская область, г. Долгопрудный, Институтский пер., 9.), karlov@shade.msu.ru, А^^гГО 14680;
Дмитрий Феликсович Пастухов - кандидат физико-математических наук, доцент кафедры технологий программирования Полоцкого государственного университета (211440, Республика Беларусь, Витебская обл., г. Новополоцк, ул. Блохина, 29), доцент Московского государственного университета им. М.В. Ломоносова (119991, Россия, Москва, Ленинские горы, д. 1), dmitrij.pastuhov@mail.ru, https://orcid.org/0000-0003-1398-6238, А^^гГО 405101;
Юрий Феликсович Пастухов - кандидат физико-математических наук, доцент кафедры технологий программирования Полоцкого государственного университета (211440, Республика Беларусь, Витебская обл., г. Новополоцк, ул. Блохина, 29), доцент Московского государственного университета им. М.В. Ломоносова (119991, Россия, Москва, Ленинские горы, д. 1), pulsar1900@mail.ru, https://orcid.org/0000-0001-8548-6959, А^^гГО 405109.
Information about the authors:
Natalya K. Volosova - Post-graduate Student of Moscow State Technical University. N. E. Bau-man (2nd Baumanskaya street, 5-1, Moscow, Russia, 105005), navalosova@yandex.ru, 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 (Obraztsova street, 9-9, Moscow, GSP-4, Russia, 127994), konstantinvolosov@yandex.ru, https://orcid.org/0000-0002-7955-0587, AuthorID 128228;
Aleksandra K. Volosova - Candidate of Physical and Mathematical Sciences, Chief Analytical Department "Tramplin" LLC Russian University of Transport (Obraztsova street, 9-9, Moscow, GSP-4, Russia, 127994), alya01@yandex.ru, https://orcid.org/0000-0002-0538-2445, AuthorID 607500;
Mikhail I. Karlov - Candidate of Physical and Mathematical Sciences, Associate Professor of the Department of Higher Mathematics, Moscow Institute of Physics and Technology (Institutskiy per., 9, Dolgoprudny, Moscow region, Russia, 141701), karlov@shade.msu.ru, AuthorID 14680;
Dmitriy F. Pastukhov - Candidate of Physical and Mathematical Sciences, Associate Professor of Polotsk State University (Blokhin Street, 29, Novopolotsk, Vitebsk Region, Republic of Belarus, 211440), Moscow State University named after M.V. Lomonosov (Leninskiye Gory, 1, Moscow, Russia, 119991), dmitrij.pastuhov@mail.ru, https://orcid.org/0000-0003-1398-6238; AuthorID 405101;
Yuriy F. Pastukhov - Candidate of Physical and Mathematical Sciences, Associate Professor of Polotsk State University (Blokhin Street, 29, Novopolotsk, Vitebsk Region, Republic of Belarus, 211440), Moscow State University named after M.V. Lomonosov (Leninskiye Gory, 1, Moscow, Russia, 119991), pulsar1900@mail.ru, https://orcid.org/0000-0001-8548-6959, AuthorID 405109.
