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

УДК 517.392

И. В. Бойков, А. Н. Тында, П. С. Краснощеков

ОПТИМАЛЬНЫЕ ПО ТОЧНОСТИ МЕТОДЫ РЕШЕНИЯ НЕКОТОРЫХ КЛАССОВ СЛАБОСИНГУЛЯРНЫХ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ ВОЛЬТЕРРА

Аннотация. Целью работы является построение оптимальных по порядку по точности методов решения интегральных уравнений Вольтерра со слабосингулярными ядрами из различных классов функций. Так как вопросы построения оптимальных по точности алгоритмов тесно связаны с оптимальной аппроксимацией функций из классов решений, то в работе применяется аппарат поперечников Бабенко и Колмогорова компактов из классов решений как в одномерном, так и в многомерном случаях. Вычислены порядки поперечников Бабенко и Колмогорова компактов из классов решений одномерных и многомерных уравнений с рассматриваемыми слабосингулярными ядрами. Построены оптимальные по порядку точности численные методы сплайн-коллокационного типа. Полученные теоретические оценки подтверждаются приведенными в заключении численными примерами решения двумерных интегральных уравнений Вольтерра.

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

I. V. Boykov, A. N. Tynda, P. S. Krasnoshchyokov

NUMERICAL METHODS OF OPTIMAL ACCURACY FOR WEAKLY SINGULAR VOLTERRA INTEGRAL EQUATIONS

Abstract. Objective: the main aim of this paper is the construction of the optimal with respect to accuracy order methods for weakly singular Volterra integral equations of different types. Methods: since the question of construction of the accuracy-optimal numerical methods is closely related with the optimal approximation problem, the authors apply the technique of the Babenko and Kolmogorov n-widths of compact sets from appropriate classes of functions. Results: the orders of the Babenko and Kolmogorov n-widths of compact sets from some classes of functions for one-dimensional and multidimensional cases are evaluated. The special local splines realizing the optimal estimates are also constructed. The optimal (with respect to accuracy order) spline-collocation methods are suggested. Conclusions: the obtained theoretical estimates are verified by the numerical examples for 2-D Volterra integral equations adduced in the paper.

Key words: Volterra integral equations, optimal algorithms, Babenko and Kolmogorov n-widths, weakly singular kernels, collocation method.

Введение

Интегральные уравнения Вольтерра имеют огромное количество приложений в экономике, медицине, экологии [1, 2]. Численным методам решения интегральных уравнений Вольтерра, включая уравнения Абеля - Вольтерра, посвящено большое количество работ, в частности [1-7].

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

Статья организована следующим образом. В главе 1 мы вводим используемые классы функций и доказываем некоторые утверждения относительно гладкости точных решений ИУВ. Глава 2 посвящена вычислению поперечников Бабенко и Колмогорова компактов из введенных классов функций. Здесь же строятся специальные локальные сплайны, реализующие оптимальные оценки. В главе 3 описан проекционный метод для многомерных ИУВ, основанный на аппроксимации точных решений такими сплайнами. Численный пример решения двумерного ИУВ приведен в главе 4.

1. Определения и вспомогательные утверждения

1.1. Классы функций

Определение 1.1. Через Qry(ДМ), О = [0,г/, I = 1,2,. обозначать класс функций / (^,.., ^), определенных удовлетворяющих следующим условиям:

на

будем Q и

Эм/

dtV

at;1 -at;1

r <|v |< s, t е{П \ Го},

M

(P(t, Го))1

где М - некоторая константа, 0<М , t = (tl,,.,Ц); V = (11,..,^),

| V |= VI +--+ VI; 5 = г + у, £, = 0, если у - целое; 5 = г + [у] +1, у = [у] + ц,

0<ц<1, £, = 1-Ц, если у - нецелое; Го - пересечение границы Г области

О со множеством координатных плоскостей; р^,Г0) = тт | ti |.

1

Замечание 1.1. В одномерном случае (I = 1) Г0 определяется как точка

t = 0.

ГГ I

Определение 1.2. Через Qг^(О,М), О = [0,Г] , I = 1,2,., обозначим класс функций /(^,...,tl), определенных на О и удовлетворяющих следующим условиям:

ЭМ/ tl)

dtV1 -dtVi

dt} -dt/ r <| v |< s, t ф 0,

M

(P(t,0)y

где М - некоторая постоянная, t = (1,..,1); V = (11,. .,11), | V |= 11 +--------------+ VI;

5 = г + у, £, = 0, если у - целое; 5 = г + [у] +1, у = [у] + ц, 0<ц<1, С = 1-Ц,

если у - нецелое; р^,0) = тт | % |.

к=1,...,1

Замечание 1.2. Очевидно, что в одномерном случае (I = 1) классы

Г ГГ

функций Qг у (О, М) и Qг у (О, М) эквивалентны.

Определение 1.3. Пусть О = [0,т], I = 1,2,., г = 1,2,.., 0< у< 1.

Г

Функция /(^,..,tl) принадлежит классу Бгу(О,А), если выполняются следующие неравенства:

jf (tl,., tl )j< a,

djvjf ti)

dlvf ^l,. ti)

dtV1 -.а^

<

dtV1 ... dtVi

Avjjvj

jvj—r-1+Y

< Ajvj j vjjvj, о<| vj<r, tєО,

, tє {О\ Го}, r <j v j< ~,

где А - константа, не зависящая от | V |.

Определение 1.4. Пусть О = [0,т], I = 1,2,., г = 1,2,., 0<у< 1.

ГГ

Функция /(^,..,tl) принадлежит классу Бгу(О,А), если выполняются следующие неравенства:

jf (tl,., ti )j< a,

djvjf ti)

djvjf ti)

at;1... dt}

atp dtvll

Ajvj j v j

vj—r-1+Y

< Ajvj j v jjvj, о<| v j< r,

, t Ф о, r <j v j< ra,

где А - константа, не зависящая от | V |.

Замечание 1.3. В одномерном случае ^ = 1) классы функций Г ГГ

Бгу(О, А) и Бгу(О, А) также эквивалентны.

1.2. Гладкость решений

Для простоты изложения рассмотрим сначала одномерное интегральное уравнение:

x(t) - jH (t, т)x(T)dт = f (t), о < t < T.

Ниже доказана лемма, определяющая свойства гладкости точного решения уравнения (1).

Г

Лемма 1.1. 1. Пусть Н(^т) е Qr^([0,Т],1) по каждой переменной

Г

в отдельности, /

Г

уравнения (1) принадлежит классу Qr у([0,Т],М).

2. Пусть H(t,т) е Бгу([0,T], A1) по каждой переменной в отдельности,

*

/

*

уравнения (1) такое, что x(t) е Br у ([0,T], A).

Доказательство. Хорошо известно, что точное решение (1) является

te[0,T] te[0,T]

непрерывной функцией. Пусть С0 = max | x(t)|, Fk = max | /(k)(t)|

Hk =

max

д kH

dtk1dTk2

0 < к1, к2 < к , к1 + к2= к, 1 < к < г . Через Мк

(г,т)е[0,Т

будем обозначать константы, зависящие только от порядка к.

Для доказательства утверждений данной леммы формально

t

продифференцируем выражение х(t) = / ^) + |н (t, т) х(т)ё т:

x'(t) = /'(t) + jH't (t, т) x(x)d т + H (t, t) x(t).

Так как правая часть этой формулы есть непрерывная функция, то существует непрерывная производная х'^) точного решения х(^), которая

может быть оценена следующим образом: | х'^) |< ^ + Н0С0 + Н1С0Т = М1. Дифференцируя (2) еще раз, имеем

х''(t) = /'V) + (Н + НТ) х^) + Н (t, t) х' (t) +

+jH't't (t, т) x(т)d т + H't (t, t) x(t).

Принимая во внимание, что

Hk, Fk =

0(1), ^ 1

0

k—r—Z

k > r,

и используя предыдущую оценку для | х'^) |, получаем

| х''(t) |< ^2 + 2Н1С0 + Н0(^1 + Н0С0 + Н1С0Т) + Н2С0Т + Н1С0 =

= ^2 + Н2С0Т + 3Н1С0 + Н 0( ^1 + Н 0С0 + Н1С0Т) = М 2.

Далее:

| х'"(0 |< ^3 + 7Н2С0 + 5Н11 х' | +Н0 | х'' | +Н3С0Т = М3.

Таким образом, дифференцируя (1) к раз (к < 5), можем заключить, что существует х(к )^) и справедлива оценка

x(k )(t) j< \

Mk, k<r

Mk

r < k < s.

Лк-г Г

Следовательно, х(t) е Qry([0,Т],М) с М = тах{Мк ,1}, т.е. первое

к=1,5

утверждение Леммы 1.1 доказано.

г

Рассмотрим теперь класс Бгу([0,Т]). Нетрудно заметить, что | х ;(г)|

ограничена. Оценка для производной порядка V, г < V < ^, имеет следующий вид:

| х^^) |< С V +С2У С+т +4 С32—^"1—,

1 '' ^-г-1+у ^-г-1+7 ^-г-2+7 ^-г-1+7 ' 7

где константы С1, С2, С3 зависят только от А1, А2 .

Г

Следовательно, точное решение х(г) е Бг у([0,Т], А).

Утверждения относительно гладкости точных решений слабосингулярных ИУВ распространяются на случай многомерных интегральных уравнений следующего вида:

tl г1

0 0

где t ), т = (т1,..,т), 0< г^...,tl <Т ; слабосингулярные ядра

g (г - т) могут иметь вид

g

или

g Ц ) =

Применяя аналогичную технику доказательства с гораздо более

объемными выкладками, заключаем:

1. Если g(?[,...,^ ) имеет вид (5), к(г,т) имеет непрерывные частные

Г

производные до порядка 5 = г + [у] +1 и /(л) е Qry (О,1), то точное решение

Г

х(г) уравнения (4) принадлежит классу Qry(О,М) с у = 5 - г -а . Если же

ГГ к(г, т) - аналитическая функция, /(г) е Бг у (О,1), то х(г) е Бг у (О, А).

2. Если g(?[,...,^) имеет вид (6), к(г,т) имеет непрерывные частные

ГГ

производные до порядка 5 = г + [у] + 1 и /(г) е Qry (О,1), то точное решение

ГГ

х(г) уравнения (4) принадлежит классу Qry(О,М) с у = 5 - г -а . Если же

ГГ ГГ

к(г, т) - аналитическая функция, /(г) е Бг у (О,1), то х(г) е Бг у (О, А).

Заметим, что свойства гладкости точных решений многомерных слабосингулярных интегральных уравнений Фредгольма были исследованы Г. М. Вайникко в работе [8].

2. Оптимальное восстановление функций из классов Qr y

Для построения оптимальных по точности методов численного решения ИУВ необходимо сначала разработать оптимальные методы аппроксимации функций из соответствующих классов. Отметим, что в работе [9] получены оценки поперечников Бабенко и Колмогорова для достаточно широкого набора классов функций, которым принадлежат решения интегральных уравнений Вольтерра со степенными и логарифмическими особенностями ядер.

В данной главе вычисляются n -поперечники Бабенко и Колмогорова

r rr r

компактов из классов функций Qry (ДM), Qry (ДM), Бгу (Д A) и

rr

Бг у (Д A), а также строятся локальные сплайны, реализующие оптимальные

(по порядкам погрешности) аппроксимации. Полученные в данной главе результаты являются распространением результатов работ [9-13].

Напомним определения n -поперечников Бабенко и Колмогорова.

Пусть Б - банахово пространство, X с Б - компактное множество, П : X ^ X - отображение X с Б на конечномерное пространство X.

Определение 2.1 [14]. Пусть L - множество n -мерных линейных подпространств пространства Б ; n -поперечник Колмогорова dn (X, Б) определяется выражением

dn

L XEXuel!1

где внешний инфимум берется по всем n -мерным подпространствам L.

Определение 2.2 [14]. n -поперечник Бабенко 5n (X) определяется выражением

5n (X) = inf supdiam П-1П(х),

n:X ^RnxEX

где инфимум вычисляется по всем непрерывным отображениям П : X ^ Rn .

Если инфимум в (7) достигается для некоторого L, то такое подпространство называется экстремальным подпространством.

Поперечники играют важную роль в численном анализе и теории аппроксимации, поскольку тесно связаны с такими проблемами, как е -сложность интегрирования и аппроксимации, оптимальное дифференцирование, оптимальная аппроксимация решений операторных уравнений. Подробное исследование этих проблем в свете общей теории оптимальных алгоритмов проведено в [15].

В данной работе через A и Ak, k = 1,2,., будем обозначать некоторые положительные константы, не зависящие от N .

Теорема 2.1. Пусть О = [0,Т]1, І = 1,2,.. Справедливы оценки (б*,т (О,М), С (О)) ~ 5И (е;,т (О,М), С (О)) ~ гп,

где єп ~ n s , если l = 1, и

Г(s-у)/(і-1) v >.

l

n-s/l (ln n)s/l, v =

- s/l

v <

l -1 l ; l -1;

l

l -1,

если І Ф1. Здесь V = 5/(5 -у) (здесь и ниже запись вида є~ X означает асимптотическую эквивалентность, т.е. є = О(Х)).

Доказательство. Для того чтобы оценить точную нижнюю грань поперечников Бабенко, разделим область О на части А к, к = 0,1,., N — 1. Здесь Ак означает множество точек из О, удовлетворяющих условиям

N1 T<p(t,Го)<fkN1 I T, k = 0,1,.,N-1, t = (tl,.,ti).

'к +1У ( к

Пусть Ик = І N І Т — І N І Т, к = 0,1,., N — 1. Каждую из областей

Ак покроем кубами и параллелепипедами Ак і с ребрами. превосходящими величин Ик и параллельными координатным осям.

В каждом кубе Ак . =[ак ,ак ,, ;...;ак ,ак м ], к = 0,1, — 1, построим

і1 ,.,іі і1 і1+1 іі іі+1

функцию следующего вида:

(('1 — <*<+1 -'1)-('і — а1!. Ха* — ',))5

А________1 1_______________1 1_________ ' рАк •

И5(2і—1)((к +1)/ N )п , і1,.,іі ’

не

. (t) =

l1,.,ll

0, t gA'

k.

ll,...,ll.

к г

Константа А выбирается из условия ¥. . е Qг у (О,М).

Оценим теперь максимум функции ¥к . (г). Очевидно, что

l1,—,ll

max j ^k , (t)j> A1 . k +1

;єО l1,.,ll I k +1

si N Уу = A (k + е)(v-1)s 1

0< Є <1.

Значение V выбирается так, чтобы тах | ¥. . (г) | не зависел от

г 11,-,Ч

номера к . Для этого достаточно выполнения условия V = 5/(5 - у).

є ~

n

Поэтому

l1,",ll м5

Определим функцию (г) формулой

N -1

¥х(г)= I I Ьм ,¥к л«• , ^1

к=011,''' ,1/

Применяя теорему Борсука [14], получаем

5n (Qr% (Q, M), С (Q)) ,

r Ns

где п - количество кубов Дк I , покрывающих О . Нетрудно видеть, что для л = 2,3,..

п/-1

N-1 k=0

гI г

N-1

= 2' I

k=0

/-1

(k + 1)v - kv

N-1

~ Nv(l-1) + 2l I k=1

Nv - k'

l-1

v-1

Nv(l-1), l v > ; l-1

Nl ln N, l v l -1;

Nl, v < l . l-1

Следовательно,

S„ (Q,.T (Q, M), С (Q)) >

Я-(»-7)/(1-1), v

l - 1 n-s/l (ln n)s/l, v = —; l -1 l

-s/l

v <•

l-1

Построим непрерывный локальный сплайн, реализующий оценку (11). Сначала более подробно опишем случай I = 1.

Пусть N и п - целые числа, связанные соотношением

к

п = (М -1)(5 -1) + 5 . Введем разбиение сегмента [0,Т] точками v^ = Т(^)9 ,

к = 0,1,...,М, где q = 5/(5-у), если у - целое, или q = 5/(5 - [у] -1), если у -нецелое.

n

№ 2 (26), 2013 Физико-математические науки. Математика

Обозначим через Ак сегменты Ак =^к,Vk+l], к = 0,1,.,N — 1. Пусть

Ьк = ^+1 + ^ ^+1 — ^ у ^ 2 + 2 Уу,

і = 1,2,..,5 — 2; ^0= Vk, Ьк—1 = Vk+l; к = 1,2,.,N — 1,

Ь0 = ^ ^ + V1 — V0 ^ у = і 2 г 2 Ь0 = „ Ь0 = V Ьі 2 + 2 , і 1,2,.,г 2; Ь0 V0, ?г—1 ^,

где у і и - нули полиномов Лежандра степеней (5 — 2) и (г — 2) соответственно.

Через Р5 (/, Ак) обозначим оператор, заменяющий функцию /('),

г ёД. , интерполяционным многочленом степени (5 -1) для к = 0,м -1,

‘г к

построенным по узлам с,у .

Через /м (г) обозначим локальный сплайн, определенный на отрезке [0, Т] и составленный из полиномов Р5 (/, Д.), к = 0,1,., N -1.

Нетрудно заметить, что

Ilf (t) - fN (t)||C[Ak ] < AN-s, k = 1, N -1. Для сегмента A0 имеем оценку

Ilf (t) - fN(t)|C[A0] < ^ v0) = A2

r f T \r

T

Nq

_ A3 _

Nqr

Таким образом, ||/(г) - /м (г)|С[0Т] - АМ 5 . Поскольку общее число

функционалов п, используемых при построении сплайна /м (г), оценивается как п ~ N , получаем

5п ^ (О))~ п-5.

Случай I = 2,3,..

Введем оператор

Р[/,[а.,61;.;а,,Ъг]] = Р1 •••РЛ ,

где РЛ - интерполяционный полином, действующий по переменной г., I = 1,2,.,I, и построенный, как и ранее, в одномерном случае I = 1. Построение непрерывного локального сплайна начинаестя с области Д N-1.

В этой области функция /(г.,.,г,) заменяется интерполяционным полиномом Р5 [/, ДN-.]. Для построения локального сплайна в ДN-2

покроем эту область кубами и параллелепипедами Д^-2. . Заметим, что их

11,..,11

ребра не превосходят hN-2 . Здесь вершины ДN-1, расположенные на границе

с ДN-2, также являются вершинами соответствующих кубов из множества Д^-2. .

В Д^-2. функция /(г.,.,г,) аппроксимируется интерполяционными

11,—,11

полиномами Р., [ /, ДN-2. ].

5 11--,11

Заметим, что мы интерполируем функцию Р5 [ /, Д N -1] вместо / (г.,., г1) в г ё {Д / Р|Д N -1} (это необходимо для соблюдения условия

непрерывности сплайна).

В областях Д., к = 0,1,.,N - 3, сплайн строится затем аналогичным образом.

Г

Обозначим через / (г.,.,) сплайн, составленный из полиномов Р5 [/, Д. , ]. Очевидно, что

f

< AN-s .

Из этой оценки и неравенства (11) следует, что выполняется правая часть соотношения (8). Используя неравенство 5п - 2dn [14], завершаем доказательство.

Теорема 2.2. Пусть О = [0,Т]. Справедливы оценки

5п

Доказательство. Построим сначала непрерывный локальный сплайн, реализующий оценку (13). Это позволит получить верхнюю оценку

Г

п-поперечника Колмогорова dn (Бгу (О, А), С (О)).

Отрезок [0,Т] разделим на N +1 частей Д. = \у.^.+.], к = 0,.,N,

узлами V) = 0, Vk =2к-1-Т, к = 1,.,N +1.

Пусть

Рк = ^+1 + ^ ^+1 - ^ у ^ 2 + 2 Уу,

у = 1,2,...,шк -2; Й = Vk, £.-1 = Vk+l; к = 0,1,.,",

где у у - нули многочленов Чебышева первого рода степени ш. - 2 ; щ = г, ш. =[1° к (г + 1-у)АТ ] +1, к = 1,2,.,N .

Через р. (/, Дк) обозначим, как и ранее, оператор, интерполирующий функцию /(г), г ё д. , с помощью полиномов степени (шк -1), построенных по узлам , к = 0^. Обозначим также hk = Vk+l - V., к = 0,1,., N.

Через ^ (г) обозначим локальный сплайн, определенный на отрезке [0,Т] и составленный из полиномов Р. (/, Д.), к = 0,1,.,N.

Нетрудно видеть, что оценка погрешности аппроксимации на сегменте г ёД(0 имеет вид

/ Т у+1-У а

Поскольку степени (шк -1) интерполяционных полиномов растут пропорционально номеру к сегмента, необходимо оценить константу Лебега Хт. для системы узлов (14). Эта константа необходима для определения

погрешности аппроксимации на участках Д., к = 1, N .

Хорошо известно, что константа Лебега не зависит от длины отрезка, а только от распределения узлов на нем. Поэтому для простоты обозначений рассмотрим отрезок [-а, а] и узлы

гу = ауу, у = 1,2,.,т. - 2; г0 = -а, щ- = а,

где у у - нули многочленов Чебышева первого рода степени (т. - 2).

Тогда

mk -1

^mk = max I

?g[—a,a] y=o

Vmk -1,i

mk -1

где у,..-ц (г ) = П (- —--------фундаментальные многочлены.

у=0,у*-

Пусть I Ф 0, I Ф т. -1, тогда

<

I Vmk -1,i(t) I

Ф mk -1,i

Ф mk -1,i

max 1 9m, -1,i(t)1, fe[-a,a]

где

Ф m> -1,i

"k 1>,v' (ti -t1)^(ti -ti-1)(ti -ti+1)^(ti -tmk-2)5

9m,, -1,i (t)=- k

k ’ (ti -t0)(ti -tmk-1)'

Очевидно, что

max I9mk-1,i(t)I=—4---------772 2 a2

fe[-a,a] k |(ti + a)(ti - a)| a2 -1 a2 -a2co„2

2

1

1

1 - cos2-

2

Принимая во внимание, что

sin2

1

mk, Vl'.

2

max j fM, -і,о(t)j~ max j Уm, -1,m, -1(t) j~ mk, te[—a,a] ;є[—a,a]

имеем

mk -2

V. < A4m2 I

l=1

Ф m. -1,i

+ j Vm. -1,0 (t)j + j Vm. —1,mk -1(t)j O(mk ln mk ).

Следовательно, погрешность аппроксимации ||f(t) - fN (t)|C[a ] =

k = 1,N — 1 может быть оценена следующим образом (см., например, [16]):

vk+1 - vk

где q

\\f (t) - fN (t)||C[A. ] < M 5(r +1 -у)к

2mq

k

а-г-;+у ’

+1

максимальный порядок производных,

используемых для оценки погрешности.

Продолжая предыдущее неравенство, получаем

A6mk: ln mk

f T2k-1-N

Ilf (t) - fN (t )| C[ Ak ] <

2mk

Aqqq

q-r —1+у

4

A7 k2 ln kTqAqqq

mq 2(N+2-k)q2(k-1-N)(q-r-1+у) ASk2ln k TqAqqq

ASk2ln k qq

Ask ln k

2

Agk ln k

<

r+1—у 2N(r+1-у) 2 9

Следовательно, для достаточно больших значений N справедливы оценки

II/(г) - Л(г)С[дк ] < 2NЩ-Y).

Таким образом, для всего отрезка [0,Т ] получаем

Общее число п функционалов, используемых при построении сплайна, оценивается как

N N 10

п = 2т. = 2 "9 к(г+1 -у)АТ ~ н2.

к=0 к=0

Поэтому

бп ^ (о, а), с (о)) >

Неравенство (16) позволяет нам определить верхнюю границу п-поперечника Колмогорова

dn (бг т (О, А), С (О)) < ГА1\ .

2^ п

Принимая во внимание неравенство 5п < 2dn, завершаем доказательство.

Заметим, что вместо системы (14) также можно использовать другую систему узлов:

где у у - нули многочленов Чебышева первого рода степени (тк — 1); то = г, тк = [к(г +1 — у)АТ] +1, к = 1,2,..," .

Это позволит исключить дополнительный множитель т2 в оценке

константы Лебега Хт . Однако замкнутая система узлов (14) более подходит

к

из практических соображений при построении численного решения ИУВ с помощью проекционного метода, описанного в разделе 3.

Теорема 2.3. Пусть О = [0,т], I = 2,3,. Тогда справедливы

следующие оценки:

5И {<2Г^

Доказательство. Для того чтобы оценить точную нижнюю грань оп ^гу (О, М), С (О)), покроем область О кубами следующим образом. Куб

А1 = Д1 1 является пересечением областей

0 < ;і <1 N I T

n• •• n

0 < ;і <1 NIT

V = s/(s - у), если у - целое, V = s/(s - [у] -1), если у - нецелое.

2 2 Область Д определяется затем как Д = Д2' \ Д.", где

Ак' = ,;і): 0<tl,,;і <| n I Tі,

Ak" = ,;і): 0<tl,—,;і <| n I Tі.

Эта область также покрывается кубами и параллелепипедами А

;,—,li

ребра которых параллельны осям координат и не превышают

( 2 У ( 1 У

\ = І -^ \ Т — І ^ I Т . Дальнейшее построение проводится по аналогии.

Каждая область Ак = А к' \ А(к—1)'', к = 3,., N — 1, покрывается кубами

величины

и параллелепипедами Д. . , ребра которых не превосходят

н--1=[ N 1Т - (^ ГТ-

В области Д. определяется функция ¥к . , затем в области О г1 ,-',г/ г1,-,г/

вводится функция (г.,-, ^) (по аналогии с доказательством теоремы 2.1).

А1

Затем получаем оценку | (г.,...,^) |>--.

N

Определим число п параллелепипедов Д. . . Очевидно, что

г1 ,"•,*/

N-1

n ~ I

к=1

к +1

N

l-1

к +1V fk^v

N

N

N-1

~ I

к=1

l-1

v-1

N-1

~ I kl-1 ~ N. к=i

ГГ -е/1

Следовательно, 5п (<2гу (О, М), С (О)) > Ап . Построение локального

ГГ

сплайна /N (г.,.,^) и дальнейшие рассуждения проводятся по аналогии с доказательством теоремы 2.1. Теорема доказана.

Теорема 2.4. Пусть О = [0,Т]1, I = 2,3,-, 0<у< 1. Тогда справедливы следующие оценки:

5n

1

п\^г ,y v

(r+1-Y)/(l-1) •

Доказательство. Обозначим через Ад множество точек ґеО таких,

что 0 < р(ґ, Гд) <2—"т, а через Ак, к = 1,2,.,", множество точек ґ еО , удовлетворяющих неравенствам

2 к—1 2 к

—Т <р(ґ, Г0) < ^Т.

Покроем каждую из областей Ак, к = 0,1,.,N, кубами Ак . . Ребра

\,"-,11

этих кубов параллельны ребрам О, не меньше Нк и не больше, чем 2Нк, где

k-1

N

-T, к = 0,1,..,N-1.

Оценим общее число элементов А. . , покрывающих область О.

г1 ,"',11

Очевидно, что

N

n ~ Z к=1

1-

k-1

N

l-1

2к 2k—1

2N 2N

N

~ Z(2

к=1

N-к+1

l-1

N 0

~ Z2

k=1

кк

2l -1 -1

,(N+1)(l-1) — 2l-1

Таким образом, п ~2*(1 1).

Повторяя рассуждения статьи [12], получаем оценку 5п

и заключаем, что

1

5n (В*,у (П, А), С (П)) >

(r+1-Y)/(l-1) •

Построение непрерывного локального сплайна, реализующего оценку (19), производится подобно построению, приведенному в теореме 2.1.

10 N (г +1 -у) АТ

Параметр 5 здесь равен 5 =

+1.

Принимая во внимание хорошо известное неравенство 5п < 2ёп, связывающее п -поперечники Бабенко и Колмогорова, завершаем

доказательство.

Теорема 2.5. Пусть О = [0,Т/, I = 2,3,. Справедливы оценки 5п (<; (О, А), С (О)) ~ ёп (Б~ (О, А), С (О)) ~

1

21+Уп

Доказательство проводится аналогично доказательству теоремы 2.4.

3. Приближенное решение многомерных ИУВ

Рассмотрим многомерные интегральные уравнения Вольтерра вида

11 \

о о

где t = (1,..,/), т = (Т1,...,Т/), 0 < ^,..,tl < Т ; слабосингулярные ядра

g^ -т) которых имеют следующую форму:

/, ,ч ,г+а ,г+а

£ ('ь-. Ч ) = ґ1 '/

3.1. Численный алгоритм

Приближенное решение уравнения (21) будем искать в виде

сплайна хN ('1,., ) с неизвестными значениями 4 (^к^ ^к X

г1 і

(£к,.,£к)е Ак , к = 0,1,.," — 1, в узлах сетки.

'1 7 Ч ’•••’ ч

*

Построение сплайна Х" ('1,., ґ/) описано в главе 2 и зависит от рассматриваемого класса функций.

Значения хN (£к ,•••, ) в каждом кубе Ак . , к = 0,1,., N — 1,

г1 Ч г1,--,7

определяются последовательно с помощью метода сплайн-коллокации из

систем линейных уравнений

(I - K) Pn

x(t), А к : - PjV x(t), А k

(; Ч I; ,..., 4

-Pi

N

J- J PN [h(t, т)]g (t -т) Pn

x(т), А.

d т, Аk

l"- ’7

Ч,--,7

= Pn [f^^k . (t), Ak . ].

^-V'V .1,.,./

Здесь ^ - оператор проектирования на множество локальных

сплайнов

вида ХN ('1,.,'/); /. . ('1,., ґ/) - новая правая часть уравнения

г1,--,7

(21), включающая интегралы по областям А^ . , і = 0,1,., к,

г1,---,7

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

Все интегралы в (23) вычисляются приближенно с помощью составных кубатурных формул типа Гаусса, построенных на неравномерных сетках.

Возможны несколько способов обхода подобластей Дк . ,

к = 0,1,...," — 1. Один из таких способов, допускающий распараллеливание вычислительного процесса, предложен в работе [17].

3.2. Обоснование сходимости

Перепишем уравнение (21) и проекионный метод (23) в операторном

виде:

x - Kx = f, K: X ^ X, X с С

где X - одно из множеств Qry(П,М) или Бгу(П,А); - множества

соответствующих локальных сплайнов.

Поскольку однородное интегральное уравнение Вольтерра х — Кх = 0 имеет только тривиальное решение, оператор I — К - инъективный. Следовательно, оператор I — К имеет ограниченный обратный оператор

(I — К) 1 : X ^ X. Для достаточно больших значений N имеем оценки

-1

-1

С

-1

<

С

-1

С

< 2

-1

С

= A (const),

если

llK - Р'4

С

<

1

-1

С

Покажем, что последняя оценка имеет место для достаточно больших N. Так как у(^) = (Кх)^) е X и X - плотное множество в С(П) (что

Г Г

справедливо для классов QrY(П,М) и Бгу(П,А) ), то получаем

где є N ^ 0 при N ^га. Поэтому ||К —

1

-1

начиная

-1

с достаточно больших значений N .

Таким образом, операторы (I — PNK) 1 существуют и равномерно ограничены, а уравнение (25) имеет единственное решение для всех достаточно больших N . Принимая во внимание, что PNX ^ х при N

для всех х е X, применим оператор проектирования PN к левой и правой частям уравнения (24):

х — PNKx = PNf + х — PNX.

Вычитая это уравнение из (25), получим

(I — PNK)(XN — х) = PNX — х,

(XN — х) = (I — PNK)—1( PNX — х).

Следовательно,

— х\\С < А^х — х\\С <£N

Таким образом, точность приближенного решения, полученного с помощью проекционного метода (25), определяется точностью £ N (X) аппроксимации функций из X локальными сплайнами.

С другой стороны, из теорем главы 2 следует, что для функций из X порядок оценки (26) не может быть улучшен (теоремы 2.1 и 2.2). Следовательно, алгоритм (23) является оптимальным по порядку точности на

Г Г

классах функций Qry (П,М) и Вгу (П, А).

В работе [5] доказано, что такие численные методы для ИУВ также являются оптимальными по порядку по сложности.

Численное решение многомерных ИУВ с оптимальной точностью требует выполнения огромного числа арифметических операций. В работе [17] на примере двумерных интегральных уравнений рассмотрена задача ускорения вычислительного процесса за счет использования многопроцессорных систем.

Замечание 3.1. Предложенный алгоритм также может быть применен к численному решению многомерных ИУВ вида

Н ч

х

0 0

t = (^,...,Ц), т = (Т1,.,Т/), tе [0,Г]1, с коэффициентами из классов Qry (П, М) или Бгу (П, А).

4. Численная иллюстрация

В качестве примера рассмотрим интегральное уравнение на классе

Г

функций В2 0 5(П,1) (входные функции в данном случае также являются

Г

функциями класса Qry (П,1) с г = 2 и некоторым значением у):

^2^1 — — — 25п2/16^6

х

10485/6

0 0

2

где (^і,(2)еО = [0,1] • Точное решение уравнения (28) имеет вид х((Ъ(2) = ((1(2)2 5 •

* *

Поскольку Бгу (О, А) с 2гу (О,М), то к уравнению (28) были

применены два различных алгоритма аппроксимации для каждого из этих классов (их описания приведены в главе 2). Результаты приведены в табл. 1, 2.

Таблица 1

*

Погрешность аппроксимации на классе Q2 2 5

N 1 2 3 5 10 15 20

Є1 7,17 e-3 1,12 e-5 2,17 e-7 4,89 e-8 6,13 e-9 6,32 e-11 5,89 e-13

Є2 0,011 6,15 e-4 9,05 e-6 6,73 e-7 2,39 e-8 6,84 e-10 2,87 e-11

* Погрешность аппроксимации на классе B205 (Д1) Таблица 2

N 1 2 3 5 10 15 20

Є1 6,25 e-4 1,04 e-7 3,10 e-8 5,99 e-9 5,23 e-10 7,31 e-13 6,19 e-15

Є2 0,002 7,45 e-6 9,73 e-7 7,62 e-8 1,48 e-9 6,35 e-12 3,17 e-13

В табл. 1, 2 N - число подобластей в основном разбиении Q;

£1 = max | x(ti, tj) - Xn (tt, tj) | - погрешность в узлах сетки;

ij

£2 =||x(ti, t2) - xn(ti, t2^|c(Q) - погрешность в Q .

Список литературы

1. Baker, C. T. H. A perspective on the numerical treatment of Volterra equations / C. T. H. Baker // J. Comp. Appl. Math. - 2000. - Vol. 125. - Р. 217-249.

2. Brunner, H. Collocation methods for Volterra integral and related functional differential equations / H. Brunner. - Cambridge : Cambridge University Press, 2004.

3. Brunner, H. The piecewise polynomial collocation method for nonlinear weakly singular Volterra equation / H. Brunner, A. Pedas, G. Vainikko // Math. Comp. - 1999. -Vol. 68, № 227. - Р. 1079-1095.

4. Diogo, T. Collocation methods for second-kind Volterra integral equations with weakly singular kernels / T. Diogo, S. McKee, T. Tang // Proc. Roy. Soc. Edin. - 1994. -Vol. 124A. - Р. 199-210.

5. Tynda, A. N. Numerical algorithms of optimal complexity for weakly singular Volterra integral equations / A. N. Tynda // Comp. Meth. Appl. Math. - 2006. - Vol. 6, № 4. - P. 436-442.

6. Tynda, A. N. Spline-collocation technique for 2D weakly singular Volterra integral equations / A. N. Tynda // Bulletin of Middle-Volga Math. Society. - 2008. - Vol. 10, № 2. - Р. 68-78.

7. Verlan, A. F. Integral equations: methods, algorithms, programms / A. F. Verlan, V. S. Sizikov. - Kiev : Naukova Dumka, 1986.

8. Vainikko, G. M. On the smoothness of solution of multidimensional weakly singular integral equations / G. M. Vainikko // Mat. Sbornik. - 1989. - Vol. 180, № 12. -Р. 1709-1723.

9. Бойков, И. В. Поперечники соболевских классов функций с особенностями на границе / И. В. Бойков, А. Н. Тында // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2013. - № 1. - С. 61-81.

10. Boikov, I. V. The optimal methods of approximation of the functions and computing the integrals / I. V. Boikov. - Penza : Penza State University Publishing House, 2007. - 236 p.

11. Boikov, I. V. Approximation of Some Classes of Functions by Local Splines / I. V. Boikov // Computational Mathematics and Mathematical Physics. - 1998. - Vol. 38, № 1. - Р. 21-30.

12. Boikov, I. V. The optimal algorithms of recovery of the functions and computing of the integrals on a class of infinitely differentiable functions / I. V. Boikov // Izvestia Vuzov. Matematika. - 1998. - № 9. - Р. 14-20.

13. Boikov, I. V. Accuracy-optimal approximate methods for solving Volterra integral equations / I. V. Boikov and A. N. Tynda // Differential Equations. - 2002. - Vol. 38, № 9. - Р. 1305-1313.

14. Babenko, K. I. Theoretical Foundations and Construction of Numerical Algorithms for Problems in Mathematical Physics / K. I. Babenko. - Moscow : Nauka, 1979.

15. Traub, J. F. A General Theory of Optimal Algorithms / J. F. Traub and

H. Wozniakowski. - New York : Academic Press, 1980.

16. Dziadyk, V. K. Introduction in Theory of Uniform Approximation of the Functions by Polynomials / V. K. Dziadyk. - Moscow : Nauka, 1977. - 512p.

17. Boikov, I. V. Methods of optimal accuracy for approximate solution of second-kind weakly singular Volterra integral equations for multiprocessor computers / I. V. Boikov and A. N. Tynda // Proc. of the ICCM (Novosibirsk, 2002). - Novosibirsk, 2002. -Р. 381-388.

References

1. Baker C. T. H. J. Comp. Appl. Math. 2000, vol. 125, pp. 217-249.

2. Brunner H. Collocation methods for Volterra integral and related functional differential equations. Cambridge: Cambridge University Press, 2004.

3. Brunner, H., Pedas A., Vainikko G. Math. Comp. 1999, vol. 68, no. 227, pp. 10791095.

4. Diogo T., McKee S., Tang T. Proc. Roy. Soc. Edin. 1994, vol. 124A, pp. 199-210.

5. Tynda A. Comp. Meth. Appl. Math. 2006, vol. 6, no. 4, pp. 436-442.

6. Tynda A. N. Bulletin of Middle-Volga Math. Society. 2008, vol. 10, no. 2, pp. 68-78.

7. Verlan A. F., Sizikov V. S. Integral equations: methods, algorithms, programms. Kiev: Naukova Dumka, 1986.

8. Vainikko G. M. Mat. Sbornik. 1989, vol. 180, no. 12, pp. 1709-1723.

9. Boykov I. V., Tynda A. N. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Fiziko-matematicheskie nauki [University proceedings. Volga region. Physics and mathematics sciences]. 2013, no. 1, pp. 61-81.

10. Boikov I. V. The optimal methods of approximation of the functions and computing the integrals. Penza: Penza State University Publishing House, 2007, 236 p.

11. Boikov I. V. Computational Mathematics and Mathematical Physics. 1998, vol. 38, no.

I, pp. 21-30.

12. Boikov I. V. Izvestia Vuzov. Matematika. 1998, no. 9, pp. 14-20.

13. Boikov I. V., Tynda A. N. Differential Equations. 2002, vol. 38, no. 9, pp. 1305-1313.

14. Babenko K. I. Theoretical Foundations and Construction of Numerical Algorithms for Problems in Mathematical Physics. Moscow: Nauka, 1979.

15. Traub J. F., Wozniakowski H. A General Theory of Optimal Algorithms. New York: Academic Press, 1980.

16. Dziadyk V. K. Introduction in Theory of Uniform Approximation of the Functions by Polynomials. Moscow: Nauka, 1977, 512 p.

17. Boikov I. V., Tynda A. N. Proc. of the ICCM (Novosibirsk, 2002). Novosibirsk, 2002, pp. 381-388.

Бойков Илья Владимирович

доктор физико-математических наук, профессор, заведующий кафедрой высшей и прикладной математики, Пензенский государственный университет

Б-шаіІ: math@pnzgu.ru

Тында Александр Николаевич

кандидат физико-математических наук, доцент, кафедра высшей и прикладной математики, Пензенский государственный университет

Б-шаіІ: tynda@pnzgu.ru

Краснощеков Павел Сергеевич доктор физико-математических наук, профессор, действительный член Российской академии наук, главный научный сотрудник, Вычислительный центр имени А. А. Дородницына Российской академии наук

Б-шаіІ: math@pnzgu.ru

Boykov Il'ya Vladimirovich Doctor of physical and mathematical sciences, professor, head of sub-department of higher and applied mathematics,

Penza State University

Tynda Aleksandr Nikolaevich Candidat of physical and mathematical sciences, associate professor, sub-department of higher and applied mathematics, Penza State University

Krasnoshchekov Pavel Sergeevich Doctor of physical and mathematical sciences, professor, full member of the Russian Academy of Sciences, principal scientist, Computing center named after A.A. Dorodnitsyn of the Russian Academy of Sciences

УДК 517.619 Бойков, И. В.

Оптимальные по точности методы решения некоторых классов слабосингулярных интегральных уравнений Вольтерра / И. В. Бойков, А. Н. Тында, П. С. Краснощеков // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2013. - № 2 (26). -С.87-107.