О численном решении интегральных уравнений с частными интегралами Текст научной статьи по специальности «Математика»

MS С 65R20

О ЧИСЛЕННОМ РЕШЕНИИ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ С ЧАСТНЫМИ ИНТЕГРАЛАМИ

В.А. Калитвин

Липецкий государственный педагогический университет, ул. Ленина, 42, Липецк, 398020, Россия, e-mail: kalitvin@gmail.com

Аннотация. Построены а;п'оритмы численного решения .линейных и нелинейных интегральных уравнений с частными интегралами, содержащими постоянные и неременные пределы интегрирования. Алгоритмы основаны на применении метода механических квадратур. Приведены теоремы о сходимости этого метода для рассматриваемых классов уравнений.

Ключевые слова: интегральные уравнения, уравнения с частными интегралами, интегро-дифференциальные уравнения Барбашина, метод механических квадратур.

1. Введение. К линейным интегральным уравнениям с частными интегралами

х(Ь,в) = / 1(Ь,в,т)х(т,в)йт + / т(Ь, в,а)х(Ь,а)йа+

о а ос / -• \

р Ь пд, ^

+ и{Ь, 8,т,о)х{г,о)йтйо + f (1,в)

ас

и их частным случаям приводятся задачи интетро-дифференциальных уравнений Бар-башина, механики сплошных сред и ряда других прикладных задач 11—5|.

К нелинейным интегральным уравнениям с частными интегралами вида

рЬ рЬ

х(Ь,в) = / с(т, в)х(т, в)йт + / / к(т, в,а,х(т,а))йтйа + f ^,в) (2)

о а о а о с

приводится задача Коши для нелинейных интетро-дифференциальных уравнений Бар-башина |2|,

В общем случае найти явное решение интегральных уравнений (1) и (2) не представляется возможным. Поэтому важное значение имеет разработка приближенных и численных методов решения этих уравнений.

Применение хорошо известных численных методов решения линейных интегральных уравнений второго рода к интегральным уравнениям (1) требует осторожности, так как известные обоснования этих методов часто связаны с предположением о компактности интегральных операторов, содержащихся в таких уравнениях, которой не обладают частично интегральные операторы, определяемые первым и вторым слагаемыми в правой части уравнения (1), даже в случае ненулевых непрерывных ядер 1(1, в, т)

Работа поддержана Минобрнауки России (проект №1.4407.2011)

и т(Ь,в,а). В частности, при обосновании метода механических квадратур для решения линейных интегральных уравнений второго рода |6-8| используется компактность линейных интегральных операторов, содержащихся в этих уравнениях.

Аналогичные проблемы возникают при обосновании численных методов решения нелинейного интегрального уравнения (2).

В связи с этими обстоятельствами в работе строятся алгоритмы численного решения линейных и нелинейных интегральных уравнений с частными интегралами, основанные па применении метода механических квадратур, и приводятся теоремы о сходимости этого метода дня уравнений рассматриваемых классов.

2. Численное решение линейного интегрального уравнения с частными интегралами. Будем рассматривать интегральное уравнение

с частными интегралами, где £ € [а,Ь],в Е [с, , заданные функции с(т,в),к(т,в,а), f (Ь,в) и К(г,в) непрерывны по совокупности переменных, а интегралы понимаются в смысле Лебега. Уравнение (3) есть частный случай линейного уравнения (1) с частными интегралами.

Пусть Б = {(т,в) : а < т < Ь < Ь,в Е [с, /!]}. Через С (Б) обозначим множество непрерывных на треугольнике Б функций с супремум нормой, а через X — множество функций из С (Б), имеющих непрерывную частную производную по Ь. С (Б) и X — банаховы пространства относительно норм ||х||с(д) = 1х(Ь, в)| и ||х||х =

тахр(1х(Ь, в)| + |х£(£,в)|) соответственно.

Под решением уравнения (3) будем понимать непрерывную функцию х(Ь, в), подстановка которой в уравнение (3) обращает это уравнение в тождество.

В силу [4,5] уравнение (3) имеет единственное решение х Е С (Б). Тогда справедливо тождество

Так как подынтегральные функции в правой части тождества (4) непрерывны, то она дифференцируема по Дифференцируя по £ тождество (4), получим тождество

Следовательно, решение интегрального уравнения (3) является решением иптетро-дифференциального уравнения Барбашина (ИДУБ)

х(Ь,в) = с(т, в)х(т, в)йт + I / к(т, в,а)х(т,а)йтйа + f (Ь,в)

(4)

(6)

с начальным условием

х(а, в) = f (а, в).

(7)

Как и выше, под решением ИДУБ (6) здесь и далее понимается функция х Е X, подстановка которой в уравнение (6) обращает его в тождество. х

(5) и равенство (7). Интегрируя обе части тождества (5) по отрезку [а,{\ и учитывая условие (7), получим тождество (4), которое показывает, что решение задачи Коши

(6)/(7) является решением интегрального уравнения (3). Таким образом, интегральное уравнение (3) и задача Коши (6)/(7) эквивалентны.

Задача Коши (6)/(7), очевидно, эквивалентна следующему двумерному интегральному уравнению:

x(t, s)

nt t-d

r(t, s, т, a)x(т, a)dтda + g(t, s) = (Rx)(t, s) + g(t, s)

(8)

где

тit's'т'a)= к(т,^°)ехр{. [ci‘,s

g(t,s)= й(т,^ехр{ / + f (a,s)exp{ / c(C,s

Уравнение (8) является линейным интегральным уравнением с непрерывным ядром r(t, s, т, а) и непрерывной функцией g(t, s). Оно имеет единственное решение в C(D), так как спектральный радиус компактного в C(D) интегрального оператора R равен нулю |1,4,5|, и это решение можно найти методом последовательных приближений.

Дня численного решения уравнения (8) могут быть использованы многочисленные методы решения линейных интегральных уравнений, в частности, метод механических квадратур.

При применении метода механических квадратур к уравнению (8) отрезки [а,Ь] и [c, d] разбиваются на части точками

tp = а + ph(p = 0,1,..., Р,а + Ph < Ь < (P + 1)h) ,

sq = c + qg(q = 0,1,...,Q,c + Qg < d< (Q + 1)g),

в уравнении (8) заменим t и s на tp и sq соответственно, а интеграл вычислим по формуле

p Q

r(tp, sq, т, a)x(т, a)dтda = hg ЕЕ Ipqijrpqijx(ti, sj) + r:

i=0 j=0

pq

(9)

где = т(ір,вд), а Трд — остаток в формуле (9). В результате, интегральное уравнение (8) заменяется системой уравнений относительно неизвестных

х{и, ) (г = 0,; і = 0,1,...,0).

t

Отбрасывая в этой системе уравнений остатки, получим систему уравнений для приближенных значений Хрд функции х в точках (ір, вд)

гДе 9ря = 9(1р, 8д), & дрд — погрешности вычислений для уравнений системы (9) с хрд.

Аналогично в |9| доказывается

Теорема 1. Пусть в формула (9) остатки стремятся к нулю равномерно относительно р, д при к,д ^ 0, |I — А < ж н погрешности вычислений стремятся к нулю равномерно относительно р, д при к,д ^ 0. Тоща при всех достаточно малых к и д приближенное решение хРд (р = 0,1,... ,Р; д = 0,1,... ,О) может быть найдено из системы (10), причем для любого заданного £ > 0 существуют такие ко и д0, что при к < к0 и д < до 1хрд - х^р,8д)! <£ (р = 0,1,...,Р; д = 0,1,...,О).

Таким образом, при численном решении линейного интегрального уравнения (3) с частными интегралами целесообразно перейти к эквивалентному линейному двумерному интегральному уравнению (8) и численно решать это уравнение с применением методов численного решения линейных интегральных уравнений, например, с применением метода механических квадратур.

Отметим, что при сделанных выше предположениях описанный процесс численного решения линейного интегрального уравнения (3) непосредственно применяется и к численному решению задачи Коши дня линейного иптегро-дифференциального уравнения Барбашипа.

3. Численное решение нелинейного интегрального уравнения с частными интегралами. Результаты раздела 1 естественным образом распространяются па нелинейные интегральные уравнения с частными интегралами (2), где Ь Е [а,Ь\,в Е [о,&\. и Е (-ж, +ж), заданные функции с(т, в),к(т, в,а,и), f (Ь, в) и функция непре-

рывны по совокупности переменных, функция к(т, в, а, и) удовлетворяет условию Липшица

а интегралы понимаются в смысле Лебега. Линейное интегральное уравнение (3) с частными интегралами есть частный случай уравнения (2). Оно получается из (2) при к(т, в, а, и) = к(т, в, а)и.

Под решением уравнения (2) будем понимать непрерывную функцию х(Ь, в), подстановка которой в уравнение (2) обращает это уравнение в тождество.

В силу [4] уравнение (2) имеет единственное решение х Е С (Б). Также как в разделе

1 показывается, что уравнение (2) эквивалентно интегро-дифференциальпому уравнению Барбашипа

р Я

(р = 0,1,...,р; д = 0,1,-..,Я), (10)

і=0 ]=0

\к(т, в, а, и) — к(т, в, <7,ь)\ < N\и — V

(и)

с начальным условном (7), где иод решенном И ДУБ (11) здесь н далее понимается функция х Е X, подстановка которой в уравнение (11) обращает это уравнение в тождество.

Задача Коши (11)/(7), очевидно, эквивалентна следующему двумерному интегральному уравнению:

х(Ь,в) = / / г(Ь, в,т,а,х(т,а))йтйа + д(Ь,в) = (Кх)(Ь,в) + д(Ь,в) , (12)

о а о с

I

га,5,г,а,и) = ехР { / <«,.)*,

Т

I I

дЦ. в) = ! !^(т. в) ехР { J с(С, в)^}^т + /(а, в) ехр { J с(£,

аТ

Уравнение (12) является интегральным уравнением Урысона с непрерывным ядром г(Ь, в, т, а, и) и непрерывной функцией д(Ь, в), оно имеет единственное решение в С(V) и это решение можно найти методом последовательных приближений |4|,

Дня численного решения интегрального уравнения Урысона (12) могут быть использованы различные методы решения интегральных уравнений, в частности, метод механических квадратур |3-7|,

При применении метода механических квадратур к уравнению (12) отрезки [а, Ь] и [с, разбиваются та части точками Ьр и вд (см. раздел 2), в уравнепни (12) £ и в заменяются на Ьр и вд соответственно, а интеграл вычисляется по «квадратурной» формуле

1>гр га р ' Я'

I I г(1р,вд,т,а,х(т,а))йтйа = кд ЕЕ 1рдг] Грдг] (х (£г ,в] )) + грд , (13)

^а ^с г=0 j=0

где трд^(х(и, в^)) = г(Ьр, вд,и, в],х(и, в^)), а грд — остаток в формуле (13). В результате, интегральное уравнение (12) заменяется системой уравнений относительно неизвестных

х(и, вз) (1 = 0,; 3 = 0,■

Отбрасывая в этой системе уравнений остатки, получим систему нелинейных уравнений для приближенных значений хрд функции х в точках (Ьр, вд)

р Я

хрд кд "У Л^У ^ Трдг]грдц (х1] ) + дрд + $рд (р 0, 1, ■ ■ ■ , Р; д 0, 1, ■ ■ ■ , Q) , (1^)

г=0 j=0

гДе дрд = д(ЬР, вд), а 6рд — погрешности вычислений для уравнений системы (14) с неизвестными хрд,

Если теперь в «квадратурной» формуле (13) остатки стремятся к пуню равномерно относительно р, д при к,д ^ 0, 0 < трд^ < А < той погрешности вычислений стремятся

к нулю равномерно относительно р, д при к,д ^ 0, то при всех достаточно малых к и д приближенное решение хрд (р = 0, 1,■■■,P; д = 0,1,■■■,Q) может быть найдено из системы (14), причем для любого заданного е > 0 существуют такие к0 и д0, что при к < ко и д < до

Таким образом, при численном решении нелинейного интегрального уравнения (2) с частными интегралами целесообразно перейти к эквивалентному нелинейному двумерному интегральному уравнению (12) и численно решать это уравнение с применением методов численного решения нелинейных интегральных уравнений, например, с применением метода механических квадратур. В частности, в силу теоремы 8 (для многомерного случая) из |10|, примененной к уравнению (12), справедлива Теорема 2. Пусть выполнены условия:

а) в «квадратурной» формуле (13) остатки стремятся к нулю равномерно относительно р, д при к,д ^ 0 и 0 < < А < то;

б) погрешности вычислений стремятся к нулю равномерно относительно р, д при

к,д ^ 0;

в) существует непрерывная частная производная дк(т, в, а, и)/ди.

Тогда при всех достаточно малых к и д приближенное решение хрд (р = 0,1, ■ ■ ■, Р; д =

0,1, ■ ■ ■, Q) может быть найдено из системы (14), причем для любого заданного е > 0

ко до, к < ко д < до

следующая оценка скорости сходимости:

(р = 0,1, ■ ■ ■, Р, д = 0,1, ■ ■ ■ ^), а х*(т,а) — единственное решение уравнения (12).

В заключение отметим, что, при сделанных выше предположениях, описанный процесс численного решения нелинейного интегрального уравнения (2) непосредственно применяется и к численному решению задачи Коши дня нелинейного И ДУБ (11).

1. Калитвин А.С. Линейные операторы с частными интегралами / Воронеж: ЦЧКИ, 2000.

252 с.

2. Appell .J.M., Ivalitvin A.S., Zabrejko P.P. Partial Integral Operators and Integra-Differential

Equations / New York-Basel: Marcel Dekker Inc., 2000. 560 p.

3. Appell .J., Ivalitvin A.S., Nashed M.Z. On some partial integral equations arising in the

mechanics of solids /7 Zeitschr. Ang. Math. Mech. 1999. 79; 10. S.703 713.

4. Калитвин А.С., Калитвин В.А. Интегральные уравнения Вольтерра и Вольтерра-Фредх'ольма

с частными интегралами / Липецк: ЛГПУ, 2006. 178 с.

хрд

x{tp,sg)! <£ ip = 0,1,...,Р; q = 0,1,...,Q).

(15)

CiSpq <

max \xpg xitp, sg)| < C2^PQ

0<p<P, 0<g<Q 14 ‘ 4 ^

где C\ и c2 - некоторые постоянные

£pq

max 0<p<P, 0<g<Q

rpg(zpgzpg^) r(tp,sg,T,*J,x &))

Литература

5. Калитвин А.С., Фролова Е.В. Линейные операторы с частными интегралами. C-теория /

Липецк: ЛГПУ, 2004. 196 с.

6. Положительные решения операторных уравнений / М.А. Красносельский, Г.М. Вай-

никко, П.П. Забрейко, Я.Б. Рутицкий, В.Я. Стеценко / М.: Наука. Главная редакция физико-математической литературы, 1969. 456 с.

7. Михлин С.Г. Некоторые вопросы теории погрешностей / Л.: ЛГУ, 1988. 336 с.

8. Даугавет И.К. Теория приближенных методов. Линейные уравнения: 2-е изд., перераб.

и доп. / СПб.: БХВ-Потсрбург, 2006. 288 с.

9. Ivalitvin V.A. On the numerical solution of Barbashin’s intcgro-diffcrential equations with

Python application /7 .Journal of Mathematical Sciences. 2013. 188; 3. P.250-255.

10. Вайникко Г.М. Возмущенный метод Галеркина и общая теория приближенных методов для нелинейных уравнений /7 Журнал вычислительной математики и математической физики. 1967. " 7; 4. C.723-75L

ON NUMERICAL SOLUTION OF INTEGRAL EQUATIONS WITH PARTIAL INTEGRALS V.A. Kalitvin

Lipetsk State Pedagogical University,

Lenina St., 42, Lipetsk, 398020, Russia, e-mail: kalitvin@gmail.com

Abstract. The solution schemes of linear and nonlinear integral equations with partial integrals containing fixed and variable bounds of integration are constructed. These schemes based on application of the mechanical quadrature method. The theorems about convergence of this method for classes of equations under consideration are reduced.

Key words: integral equations, equations with partial integrals, Barbashin’s integro-differential equations, mechanical quadratures method.