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

____________УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА

Том 156, кн. 2 Физико-математические науки

2014

УДК 519.642

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

Ю.Р. Агачев, А.И. Леонов, И.П. Семенов

Аннотация

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

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

В работе методом квадратур решаются интегральные уравнения вида

+1 /»

Kx = x(t)+ n(s)h(t, s)x(s) ds = y(t), -1 +1 /» -i < t < i, (1)

Kx = x(t) + ^(t) h(t, s)x(s) ds = y(t), -1 -i < t < i, (2)

где у, € Li( —1,1), h(t, s), y(t) - данные, x(t) - искомая функции, а также многомерные аналоги уравнений (1), (2).

Обоснование метода механических квадратур для интегральных уравнений Фредгольма второго рода (в одномерном и многомерном случаях) проведено в ряде работ (см., например, [1, 2] и приведенную там библиографию). Как известно, этот метод является одним из наиболее простых прямых методов решения интегральных уравнений. При этом, как правило, вычислительная схема метода строится в предположении непрерывности ядра интегрального оператора и правой части уравнения.

Наиболее общие результаты по методу механических квадратур получены в работах [1, 2]. В них обоснование метода проведено по-существу в случае, когда коэффициенты уравнения R-интегрируемы в своих областях определения. Однако если коэффициенты имеют разрывы второго рода, то результаты из [1, 2] перестают быть справедливыми. Этого в некоторых случаях можно избежать, если в основу метода квадратур положить так называемую интервальную квадратурную формулу (см., например, работы [3, 4] и приведенную там библиографию). Здесь мы воспользуемся наиболее простой интервальной формулой, основанной на сплайн-функциях нулевой степени (кусочно-постоянных функциях).

5

6

Ю.Р. АГАЧЕВ И ДР.

1. Некоторые вспомогательные результаты

1°. Пусть p(t) - фиксированная весовая функция на отрезке [— 1,1], L2,p = = L2,p(-1, 1) - пространство квадратично-суммируемых на [—1,1] с весом p(t) функций. В случае функции двух переменных, заданной в квадрате [-1,1]2 и квадратично суммируемой на нем с весом q(t, s), соответствующее пространство будем обозначать L2,q((-1,1)2). Известно, что L2,p с нормой

г V 11/2

\\x\\2,p = S p(t)\x(t)\2dt\ , x е L2,p,

-1

является банаховым пространством.

В пространстве L2,p введем в рассмотрение интегральный модуль непрерывно-

сти:1

U(x; S)2,p = sup ||x(- + n) - x()\\b2 J-1,1), x е L2,p.

\n\<S

Отметим здесь, что все основные свойства модуля непрерывности выполняются; в частности, w(x; S)2,p ^ 0 при S ^ 0 для любой функции x е L2,p.

2° . На отрезке [-1,1] введем в рассмотрение сетку равноотстоящих узлов

• 1 — t0 < t1 < • • • < tn — + 1, tk

2k

n

- 1,

n - 1 е N.

(3)

Пусть Sn есть оператор, который любой суммируемой на [-1,1] функции z(t) ставит в соответствие «интерполяционный в среднем» сплайн (Snz)(t) нулевой

степени на сетке (3), то есть кусочно-постоянную функцию вида

(S<0z)(t) =^2Фк (z)^k (t), Фк (Z) k=1

tk

n J z(t) dt,

tk-1

(4)

где

Имеет место

Фк (t)

1 , t е ( tk- 1 , tk ] , o, t е (tk-1,tk].

Лемма 1. Пусть весовая функция p(t) такова, что 1/p(t) также является весовой. Тогда для любой функции z е L2,p(-1, +1) имеет место неравенство

- Snzh,p ^(z; -) .

V n) 2,p

Доказательство. Обозначим через S погрешность аппроксимации, то есть S =\\z - Slz\\2,р . Тогда

n tk 2 n ^

s2 = ^ P(t)\z(t) - фк(z)\2 dt = p(t) [z(t) - z(T)] dT

k=11L k=11— t— 1

dt.

1 Считаем, что функция продолжена вне квадрата [—1, 1]2 нулем. В случае периодических функций модуль непрерывности можно ввести по формуле

ш(х; 6)2,р = sup ||x(- + n) — x(-)\\b2 „(—l,l—n) •

z

2

МЕТОД МЕХАНИЧЕСКИХ КВАДРАТУР

7

Отсюда, применяя неравенство Коши-Буняковского, получим

tk

tk

S2 < Пу - 2 ^

P(t)

k=1

tk-1

tk-1

\z(t) — z(t)|2 dr dt.

После замены r на t + r, изменения порядка интегрирования и последующей замены в первом интеграле r на —т приходим к следующей оценке:

2/п , tk

n ^ / tk

s < ny f dA ( p(t)\z(t—r)—z(t)\2dt+

k=1 0 Uk-i+r

tk T \

+ p(t)\z(t + r) — z(t)\2 dO <

tk-i 2

1

< 2 sup / p(t)\z(t + r) — z(t)\2 dt = 2 u(z;2/n)2 p,

ТИС2/n J

откуда и следует утверждение леммы. □

Следствие. В условиях леммы 1 операторы бП : ^2,Р ^ ^2Р, задаваемые формулой (4), ограничены по норме в совокупности.

Справедливость следствия вытекает из того, что в условиях леммы 1 операторы Бп сильно сходятся к единичному оператору в пространстве L2,p, и поэтому по теореме Банаха-Штейнхауса (см., например, [5, с. 129]) ограничены по норме в совокупности в указанном пространстве.

Замечание. При p(t) = 1 утверждение леммы 1 в периодическом случае имеется в работе [6].

3° . Обозначим через H оператор, задаваемый интегралом в левой части уравнения (1), то есть

(Hx)(t) = (H (hx))(t) 1

+ 1

/ p(s)h(t, s)x(s) ds.

1

Лемма 2. Пусть функция р G Li( — 1,1), p(t) = \p(t)\, —1 < t < 1, а функция h(t, s) удовлетворяет условию

h G L2,q(( — 1,1)2), q(t, s) = p(t)p(s).

Тогда оператор H действует из пространства L2,p в L2,p и является вполне непрерывным.

Отметим, что этот результат с учетом критерия компактности множества в пространстве L2,p вытекает из неравенств для нормы и интегрального модуля непрерывности в соответствующем пространстве:

\\Hx\\2,p < \\h\\2,q ■ \\x\\2,p, x G L2,p,

u(Hx; s)2,p < Ut(h; S)2,q ■ \\x\\2,p, x G L2,p, где ut(h; S)2,q - интегральный модуль непрерывности по переменной t функции

h(t, s) G L2,q (( —1 1)2) .

8

Ю.Р. АГАЧЕВ И ДР.

Замечание. Лемма 2 фактически означает, что вполне непрерывный в пространстве L2,p интегральный оператор H может иметь ядро, содержащее не только фиксированные особенности по внутренней переменной, но и подвижные особенности определенного интегрируемого порядка. Такого свойства у ядра не может быть, если рассматривать интегральный оператор (а следовательно, и уравнение) в пространстве непрерывных функций, как это делается в большинстве работ, посвященных методу квадратур.

2. Метод механических квадратур

1°. Сначала займемся построением вычислительной схемы метода квадратур и ее обоснованием для уравнения (1).

Для интеграла вида

+1

/ p(s)z(s) ds

-1

рассмотрим интервальную квадратурную формулу, полученную заменой функции z(s) на сплайн (Snz)(t):

tk

^(s)z(s) ds

k = 1

Akфк (z)? Ak

^(s) ds,

tk-i

1

где функционалы Фк определены в (4).

Приближения к решению уравнения (1) будем искать в виде сплайна

(5)

Xn(t) = £ Ск фк (t), (6)

k=1

неизвестные коэффициенты которого найдем из системы

Cj + 53 Ak Ф3 (^k (h))ck = Фз (V), j = 1, 2,...,n, k=1

(7)

где фk - функционал Фk относительно параметра s. Для вычислительной схемы (6), (7) имеет место

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

1) Р: 1/И € L1(-1, +1);

2) У € L2,p(-1, +1) ; h € L2,q (( — 1, +1)2) ; P(t) = \P(t)\,q(t, s) = P(t)P(s) ;

3) уравнение (1) имеет единственное решение при любой правой части из L2,p . Тогда система (7) однозначно разрешима при всех n, начиная с некоторого натурального no . Приближенные решения, построенные по формуле (6), сходятся к точному решению уравнения (1) в метрике L2,p со скоростью

|Х - Хп\\2,р

0<ш[ у;

1

n

+ Щ

2,Р

h; —

n

+ шs

2,q

Следствие. Если в условиях теоремы 1 существуют производные у' € L2 р, К, hS € L2,q, то скорость сходимости приближенных решений к точному характеризуется порядковой оценкой

\Х - Хп\\2,р

= O 1 I.

n

МЕТОД МЕХАНИЧЕСКИХ КВАДРАТУР

9

Доказательство. Будем рассматривать уравнение (1) в пространстве X = = L2,p (предположения 1) и 2) обеспечивают действие оператора K в пространстве X). Согласно лемме 2 в пространстве X уравнение (1) относится к уравнению второго рода с вполне непрерывным оператором. Поэтому с учетом теории Фредгольма для таких уравнений из условия 3) теоремы следует существование ограниченного двустороннего обратного оператора K-1 : X ^ X.

Далее, в пространстве X выберем при произвольно фиксированном натуральном n подпространство Xn, состоящее из элементов вида (6). Нетрудно доказать, что система линейных алгебраических уравнений (7) эквивалентна операторному уравнению

Knxn = xn + PnHPn (hXn) = РпУ (xn G Xn)i (8)

где оператор проектирования Pn = бП .

Поэтому для доказательства разрешимости уравнения (8) и сходимости приближенных решений (9) к точному решению уравнения (1) мы можем воспользоваться общей теорией приближенных методов функционального анализа (см., например, теорему 7 гл. I монографии [7]). Согласно этой теории нам достаточно показать для указанных уравнений близость правых частей и операторов на подпространстве Xn .

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

Sn = \\y — Pny\\x <V2w(y\2/n)2p ^ 0, n ^ ж.

Для доказательства близости соответствующих операторов уравнений (1) и (8) возьмем произвольный элемент xn G Xn и рассмотрим vn = ||Kxn — Knxn||x. С учетом очевидного свойства

Pn(uv) = PnuPnV, u,v G X,

и леммы 1 имеем

Vn = \\HXn — PnHPtthXn )\\x <

< \\HXn — PnHXn\\x + \\PnHXn — PnHPn(h)xn\\x <

< V2.w(Hxn; 2/n)2,p + \\Pn(H(h — Pn(h))xn\\x <

< V2Wt(h;2/n)2,q ■ \\xn\\x + \\Pn\\L2,P^L2,P ■ \\H(h — Pnh)xn\\x <

< V2 Ut(h; 2/n)2,q ■ \\xn\\x + HPn||L2,p^L2,„ ■ \\h — Pnh\\2,q ■ \\Xn\\x.

n\\L2,p^L2,p ||^ V'* J п'Ч'Ь'П ',q ' llJ'n|x + II Pn II L2,p yL2,p ■ \\h — Pn hl\2 ,q ' ll^n

Отсюда и из леммы 1 и ее следствия непосредственно вытекает оценка

vn

Oy^t(h; l/n)2,q + Us(h; 1/n)2,q} ■ \\xn\\x■

Тем самым нами доказана близость операторов уравнений (1) и (8) на подпространстве Xn :

£n = \\K — Kn\\xn—x = o\ut(h] 1/n)2,q + Us(h; 1/n)2,q} ^ 0, n ^ ж.

Таким образом, все предположения теоремы 7 гл. I монографии [7] выполнены. Остальное очевидно. □

2° . Теперь займемся построением вычислительной схемы метода квадратур для уравнения (2).

10

Ю.Р. АГАЧЕВ И ДР.

С этой целью введем в рассмотрение новые функции z(t) и g(t), связанные с искомой функцией и правой частью уравнения (2) соотношениями x(t) = p(t)z(t) и y(t) = p(t)g(t). Тогда уравнение (2) преобразуется к уравнению

+1

Kz = z(t) + j p(s)h(t, s)z(s) ds = g(t), — 1 < t < 1, (9)

-1

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

П

Zn (t) =^2 Ck фк О, (10)

к=1

неизвестные коэффициенты которого найдем с учетом интервальной квадратурной формулы (5) из системы

n

cj + 53 Акфз (®k(h))ck =фз (g), j = 1, 2,...,П. (11)

к = 1

Если система (11) имеет решение, то приближенное решение уравнения (2) восстанавливаем с учетом связи между решениями уравнений (2) и (9):

xn(t) = g(t)zn(t). (12)

Имеет место следующая

Теорема 2. Пусть выполнены предположения:

1) Р> 1/И € L1( — 1, +^;

2) y € L2,p( — 1, +^ ; h € L2,q (( —1, +1)2) ; P(t) = 1P\P(t)I q(tl s) = 1/(P(t)P(s)) ;

3) однородное уравнение, соответствующее (2), имеет лишь тривиальное решение.

Тогда для достаточно больших n система (11) однозначно разрешима. Приближенные решения, построенные по формулам (12), (10, сходятся к точному решению уравнения (2) в метрике L2,p со скоростью

\\x — Xn II 2,р

у/и;

1

n

+ Ш

2,1/р

h ; —

n

+

2,q

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

llz — zn\\2,1/p = o\u(g;1 ) + Шг(h;1 ) + us(h;1 ) V (13)

l V nJ 2,1/p \ nJ 2,q \ nJ 2,q)

Учитывая связь между решениями уравнений (2) и (9), находим

||x — Xn\\2,p = \z — zn\\2,1/p.

Отсюда и из (13) непосредственно вытекает утверждение теоремы 2.

□

МЕТОД МЕХАНИЧЕСКИХ КВАДРАТУР

11

3° . Следует отметить, что вычислительную схему метода механических квадратур можно построить непосредственно и для исходного уравнения (2). Для этого введем оператор Hi, задаваемый интегральным слагаемым в левой части этого уравнения, то есть

(Hi x)(t)

+ i

(H1(hx))(t) = p(t) h(t,s)x(s) ds.

-1

Оператор Hi обладает свойствами, аналогичными H.

Лемма 3. Пусть функции р, 1/р G Li(-1, 1), p(t) = 1/\p(t)\, — 1 < t < 1, а функция h(t, s) удовлетворяет условию

h G L2,q([—1,1]2), q(t, s) = \p(t)p(s)\.

Тогда оператор Hi действует из пространства L2,p в L2,p и является вполне непрерывным.

Заметим, что ядро интегрального оператора Hi : L2,p ^ L2,p также может иметь подвижные особенности определенного интегрируемого порядка в силу вы-

+i

бора специальной “““# функц"". Кроме того ™п>ала / z(s) ds »“еем

-i

равенство

\ 2 n,

z(s) ds = - (z).

_i k=i

Приближения к решению уравнения (2) будем искать в виде сплайна (6), неизвестные коэффициенты которого найдем из системы

2

cj + n Фз (^У2Фз (Фк (h))ck = Фз (у\ 3 = 1, 2,.../>

k = i

(14)

Имеет место следующая

Теорема 3. Пусть выполнены предположения:

1) ^, 1/И G Li( — 1, ;

2) У G L2,p( 1, +1) ; h G L2,q (( —1, +1)2) ; p(t) = 1/\^(t)1, q(t, s) = \^(t)^(s)\ ;

3) уравнение (2) имеет единственное решение при любой правой части из L2,p . Тогда система (14) однозначно разрешима при всех n, хотя бы достаточно

больших. Приближенные решения, построенные по формуле (6), сходятся к точному решению уравнения (2) в метрике L2,p со скоростью

\\x — Xnh,p

о<М у;

1

n

+ Шг 2,Р

h; —

n

+ Ш 2,q

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

2. Если функции h(t,s) и y(t) непрерывны на отрезке [—1,1], то из вычислительных схем (6)—(7) и (12), (10), (11) можно получить с обоснованием известные

12

Ю.Р. АГАЧЕВ И ДР.

схемы метода механических квадратур для уравнений вида (1) и (2) (см., например, [1, 2, 8]). Для этого достаточно функционалы Фд, аппроксимировать хорошо известными интерполяционными квадратурными формулами, в частности малыми формулами прямоугольников и трапеций.

3. В случае уравнений (1) и (2) с гладкими исходными данными построенные вычислительные схемы неэффективны. В этом случае необходимо применять прямые методы, которые не обладают свойством насыщаемости или построены на основе сплайнов высоких степеней. Однако тогда простотой вычислительной схемы метода в определенной степени придется пренебречь.

3. О методе кубатур для многомерных уравнений

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

1 1

x(t, s) + j j ^i(r)^(a)h(t, s,r,a)x(r,a) dr da = y(t,s), t,s G [—1,1],

-1-1 1 1

x(t,s) + M1(t)M2 (s) J j h(t, s,r,a)x(r,a) dr da = y(t,s), t,s G [— 1, 1]. 11

(15)

(16)

Предполагается, как и в одномерном случае, что функции ^ на промежутке [-1,1] имеют особенности интегрируемого порядка.

Для интеграла вида

1 1

/ / ji(r,a)z(r,a) dr da

-1-1

рассмотрим интервальную кубатурную формулу, полученную заменой функции z(r,a) на двумерный сплайн нулевой степени (Snmz)(r,a) = (SnooS^mz)(r,a),

ё0 a0

где Ьпж и m означают одномерные операторы, построенные в предыдущем

пункте, по первому и второму аргументам по равномерным разбиениям сегмента [— 1,1] с шагом 2/n и 2/m соответственно2:

1 1

tk sj

11

П 110 л л

ц(т, a)z(r, a) dr da « ЕЕ AkjФк] (z), Akj = / /x(r,a) dr da. (17)

k=1j=1

tk-1 sj- 1

Здесь tk = 2k/n-1, sj = 2j/m-1, а функционалы Фк] определены по аналогии с (4):

t-k sj

Фкj(z) = Пт J j z(t,s) dtds.

tk-1 s j-1

Приближения к решению уравнения (15) будем искать в виде двумерного сплайна

xnm (t,s) = J2J2 ckj Фкj (, s) (18)

k=1j=1

где фкj есть характеристическая функция прямоугольника [tk-1,tk] х [sj-1, sj].

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

МЕТОД МЕХАНИЧЕСКИХ КВАДРАТУР

13

Неизвестные коэффициенты {ckj} найдем из системы

n m

Cil + Akj *«jh))ckj = &it(y), i = 1, 2,...,n, l = 1, 2,... ,m, (19)

k=ij=i

где Ф^’^ - функционал Фу по параметрам т, а.

Теорема 4. Пусть выполнены предположения:

1) Mii ^l^ii №■> 1/М2 € Li( —1 +1);

2) y € L2,p(( —1, +1)2) = L2,p ; h € L2,q(( —1 +1)4) ; p(t,s)= lMi(t)M2(s)| ;

q(t,s,T,a) = p(t,s)p(T,a);

3) уравнение (15) имеет единственное решение при любой правой части из

L2,p ■

Тогда система (19) однозначно разрешима при всех n, m, хотя бы достаточно больших■ Приближенные решения, построенные по формуле (18), сходятся к точному решению уравнения (15) в метрике L2P со скоростью

\x-XnmW 2,р = O\u[y; h; + <xT

h N

2,q

где N = min(n, m), а uts(h; S)2,q есть модуль непрерывности функции h по переменным (t,s), вычисленный в точке S ■

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

Лемма 4. Пусть весовая функция p(t,s) такова, что 1/p(t,s) также является весовой■ Тогда для любой функции z € L2,p(( — 1, +1)2) имеет место неравенство

llz — $ПтА2,р < 2^z; N0 .

Вычислительную схему метода кубатур для уравнения (16) можно построить точно так же, как это было сделано в одномерном случае. Введем в рассмотрение новые функции z(t,s) и g(t,s), связанные с искомой функцией и правой частью уравнения (16) соотношениями x(t,s) = ni(t)g,2(s)z(t, s) и y(t, s) = g,i(t)^2 (s)g(t, s). Тогда уравнение (16) преобразуется к уравнению

+i+i

Kz = z(t,s) + J j ni(T)g2(v)h(t, s,T,a)z(T,a) dT da = g(t,s), —1 < t,s < 1, (20) -i -i

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

Приближения к решению уравнения (20) будем искать в виде двумерного сплайна

nm

znm (t, s) = ЕЕ Ckj Фkj (t,s)7 (21)

k=ij=i

неизвестные коэффициенты которого найдем с учетом интервальной кубатурной формулы (17) из системы

nm

Cil + ЕЕ A kj ФИ Akj’A (h))ckj = фи (g), i = 1,2,...,n, l = 1, 2,...,m. (22)

k=i j=i

14

Ю.Р. АГАЧЕВ И ДР.

Если система (22) имеет решение, то приближенное решение уравнения (16) восстанавливаем с учетом связи между решениями уравнений (16) и (20):

xnm(t: s) = s) • (23)

Имеет место следующая

Теорема 5. Пусть выполнены предположения:

1) Мь 1/Мъ М2, 1/М2 е Li( — 1, +1);

2) У € L2,p((-1, +1)2), h е L2,q((-1, +1)4), P(t,s) = 1/\Mi(t)p2(a)\, q(t,s,T,o) = = 1/(P(t,s)p(T,a));

3) однородное уравнение, соответствующее (16), имеет лишь нулевое решение.

Тогда система (22) однозначно разрешима при всех n, m достаточно больших. Приближенные решения, построенные по формулам (23), (21), сходятся к точному решению уравнения (16) в метрике L2,p со скоростью

IIх хпт^2,р

0< и[у/м;

1

N

+ иts

2,1/р

h; N

+ ита

2,q

Замечание. Вычислительные схемы в общем m-мерном случае для уравнений вида (1) и (2), заданных в параллепипедальной области, строятся на основе m-мерных сплайнов нулевой степени, представляющих собой суперпозицию m одномерных сплайнов нулевой степени. При этом обоснование методов проводится с использованием m-мерного аналога леммы 4 с постоянной 2т при модуле непрерывности.

Summary

Yu.R. Agachev, A.I. Leonov, I.P. Semenov. Solution of a Class of Integral Equations with Fixed Singularities in the Kernel by the Mechanical Quadrature Method.

Computational schemes for the mechanical quadrature method were constructed on the basis of a simplest interval quadrature formula for Fredholm integral equations of the second kind with fixed integrable singularities in the kernel on the internal and external variables. The theoretical and functional justification of these schemes was given. In particular, the convergence of the method in the space of functions square-integrable in the interval of integration with a weight depending on the characteristics of the integral operator's kernel was proved. The results obtained for the one-dimensional case were also extended to the multidimensional case.

Keywords: integral equation, fixed singularities in the kernel, weighted Lebesgue space, interval quadrature formula, method of mechanical quadratures, convergence of a method.

Литература

1. Вайникко Г.М. О сходимости метода механических квадратур для интегральных уравнений с разрывными ядрами // Сиб. матем. журн. - 1971. - Т. XII, № 1. -С. 40-53.

2. Габдулхаев Б.Г. К численному решению интегральных уравнений методом механических квадратур // Изв. вузов. Матем. - 1972. - № 12. - С. 23-39.

3. Шарипов Р.Н. Наилучшие интервальные квадратурные формулы для классов Липшица // Конструктивная теория функций и функциональный анализ: Сб. ст. - Казань: Изд-во Казан. ун-та, 1983. - Вып. 4. - С. 124-132.

МЕТОД МЕХАНИЧЕСКИХ КВАДРАТУР

15

4. Milovanovic G.V., Cvetkovic A.S. Nonstandard Gaussian quadrature formulae based on operator values // Adv. Comput. Math. - 2010. - V. 32, No 4. - P. 431-486. - doi: 10.1007/s10444-009-9114-y.

5. Треногин В.А. Функциональный анализ. - М.: Наука, 1980. - 496 с.

6. Габдулхаев Б.Г. Сплайн-методы решения одного класса сингулярных интегро-дифференциальных уравнений // Изв. вузов. Матем. - 1975. - № 6. - С 14-24.

7. Габдулхаев Б.Г. Оптимальные аппроксимации решений линейных задач. - Казань: Изд-во Казан. ун-та, 1980. - 232 с.

8. Канторович Л.В., Акилов Г.П. Функциональный анализ. - М.: Наука, 1984. - 752 с.

Поступила в редакцию

28.03.14

Агачев Юрий Романович - кандидат физико-математических наук, доцент кафедры теории функций и приближений, Казанский (Приволжский) федеральный университет, г. Казань, Россия.

E-mail: [email protected]

Леонов Александр Иванович - доцент кафедры высшей математики, Казанский государственный архитектурно-строительный университет, г. Казань, Россия.

Семенов Иван Петрович - кандидат физико-математических наук, профессор кафедры высшей математики, Казанский государственный архитектурно-строительный университет, г. Казань, Россия.