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

ISSN 2074-1871 Уфимский математический журнал. Том 12. № 1 (2020). С. 43-55.

УДК 519.64.7

ОБ АСИМПТОТИЧЕСКОЙ СХОДИМОСТИ ПОЛИНОМИАЛЬНОГО МЕТОДА КОЛЛОКАЦИЙ

ДЛЯ ОДНОГО КЛАССА СИНГУЛЯРНЫХ ИНТЕГРО-ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

А.И. ФЕДОТОВ

Аннотация. Среди приближенных методов решения операторных уравнений наиболее употребительными являются методы коллокаций и Галеркина. Каждый из них имеет свои преимущества и недостатки. Так, методы Галеркина применяются для уравнений в гильбертовых пространствах, и погрешности решений, полученных этими методами, имеют порядки наилучших приближений точных решений. Однако методы Галеркина не всегда конструктивны, так как для их реализации требуется вычислять интегралы, что не всегда можно сделать явно. Методы коллокаций применяются для уравнений, заданных в пространствах непрерывных функций, и поэтому всегда конструктивны. Однако погрешности решений, полученные методами коллокаций, по порядку обычно меньше, чем порядки наилучшего приближения точного решения.

В данной работе для одного класса сингулярных интегро-дифференциальных уравнений на отрезке обоснован полиномиальный метод коллокации. Для обоснования впервые для таких уравнений была применена методика сведения обоснования метода коллокаций к обоснованию метода Галеркина. Для периодического случая такая методика была впервые использована автором для обоснования метода коллокаций для сингулярных интегро-дифференциальных и псевдодифференциальных уравнений. Для уравнений, заданных на разомкнутом контуре, эта методика использована впервые. Кроме того, впервые доказана ограниченность нормы интерполяционного оператора Лагранжа в пространствах Соболева s > 1/2, с весом Чебышева второго рода. Именно этот результат позволил показать, что и для уравнений в непериодическом случае полиномиальный метод коллокаций обеспечивает такую же скорость сходимости, что и метод Галеркина.

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

Mathematics Subject Classification: 65R20

1. Введение

Арнольд и Венланд предложили в [1] оригинальную методику обоснования метода еплайн-коллокаций для периодических псевдодифференциальных уравнений в пространствах Соболева. Обоснование основано на эквивалентности метода сплайн-коллокаций и модифицированного метода Галеркина-Петрова и обосновании последнего путем сведения его к стандартному методу Галеркина. В работах |2| |0| этот подход был использован для обоснования метода сплайн-коллокацни для различных классов сингулярных интегральных и псевдодифференциальных уравнений. Было показано, что сильная эллиптичность

A.I. Fedotov, On asymptotic convergence of polynomial collocation method for one class

of singular integro-differential equations.

©Федотов А.И. 2020.

Поступила 24 июня 2019 г.

является достаточным, а в некоторых случаях (см, [7],[8]) и необходимым условием сходимости метода сплайн-коллокаций,

В работе [9] автором был обоснован полиномиальный метод коллокаций для широкого класса сингулярных интегро-дифференциальных, периодических псевдодифференциальных и систем псевдодифференциальных уравнений в пространствах Соболева, Результаты этой работы показывают, что полиномиальный метод коллокаций сходится для более широкого класса сингулярных уравнений, чем метод сплайн-коллокаций, А именно, показано, что полиномиальный метод коллокаций сходится для всех эллиптических, а не только сильно эллиптических уравнений. Более того, скорость сходимости полиномиального метода коллокаций возрастает с ростом гладкости точного решения неограниченно, в то время как рост скорости сходимости метода сплайн-коллокаций ограничен порядком используемых сплайнов,

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

2. Постановка задачи Рассмотрим сингулярное интегро-дифференциальное уравнение

*(*) + - /1 гг^т = у®, | г |< 1, (1)

■к 1 - т (т - г)

с условием

Г1 х(т)йт

т

2

0. (2)

Здесь х - искомая, а у - известные функции, заданные на интервале отрезке [-1,1] и интервале (-1,1) соответственно, А - данное действительное число, а сингулярный интеграл понимается в смысле главного значения по Коши-Лебегу.

Вспомогательные результаты

В этом разделе приведены три леммы, необходимые для дальнейшего изложения. Доказательство первой леммы имеется, например, в [10], второй - в [11]. Результаты третьей леммы являются новыми, поэтому она приведена с доказательством.

Лемма 1. Обозначим И и V линейные операторы, действующие из банахова, пространства, X в банахово пространство У. Предположим, что оператор И обратим, и выполнено условие \\У Цх ^у ||И-1||у < 1. Тогда опера тор И + V : X ^ У также

обратим, и справедлива оценка

||(И + V )-1||

у

\И-1Цу

1

НИ-!

У

Вновь обозначим X и У банаховы пространства, и пусть Хп С Х,Уп С У,п их подпространства. Рассмотрим уравнения

Кх = у, К : X ^ У,

^п^п Уп, -Кп : Хп ^ ^п , ^ 1 , 2, где К и Кп, п = 1, 2,..., линейные ограниченные операторы.

1, 2,...,

(3)

(4)

Лемма 2. Предположим, что оператор К : X ^ Y обратим, и операторы, Кп, п = 1, 2,..., равномерно сходятся к нем,у

\\к - Kn\\Xn^Y ^ О при п ^ ж.

Если dimХп = dim Yn, п = 1, 2,..., то для всех п, удовлетворяющих условию

Яп = \\К-1\\у\\К - Kn\\Xn^Y < 1,

приближенные уравнения (4) имеют единственные решения х*п Е Хп для, любых правых частей уп Е Yn, и верна, оценка

\\х* - <\\х < W^l^ (\\у - Уп\\г + qn\\y\\y),

1 Qn

где х* = К-1у - 'точное решение уравнения (3).

В дальнейшем будем, как обычно, обозначать N множество натуральных чисел, No множество натуральных чисел, дополненных нулем, a R множество действительных чисел.

Обозначим Ti(t) = cos(/ arceos t), l Е N0, t Е [-1,1], систему полиномов Чебышева первого рода ортогональных на [-1,1] с весом p(t) = (1 - í2)-1'2, t Е [-1,1]. Обозначим

, . sin((/ + 1) arceos t) , „т г

Ui(t)= " ^-—^, I Е No, t Е [-1,1],

sin(arccos t)

систему полиномов Чебышева второго рода ортогональных на [-1,1] с весом q(t) = (1 - t2)1'2, t Е [-1, 1].

Обозначим Hp+1 пространство Соболева порядка s + 1 Е R с тесом р, т.е. замыкание множества всех гладких действительнозначных функций, определенных на отрезке [-1,1], относительно нормы

1'2

\NIя-1 ЧЕ6s+1)x2(D\ ,1= \ ЕN, (5)

.я- = {Еi2(s+1)-2(o) ,1Ч1 I = 0

UGNq ) ^ ' '

а ^

ж(0) = 1 I р(т)х(т)с1т, х(1) = 2 I р(т)х(т)Т (т)с1т, I Е Н, К К

коэффициенты Фурье функции х по системе полиномов {Т ■ В пространстве Н+1 определим скалярное произведение

(f, 9)н+ = Е ?(3+1)WW1), f, 9 Е Щ+1.

le N0

С введенным скалярным произведением поетранетво становится гильбертовым пространством, причем норма (5) выражается через скалярное произведение

\NIЩ+1 = \¡(х,х)щ+1, ж Е нр+1.

Обозначим Н^ пространство Соболева порядка в Е К с весом д, т.е. замыкание множества всех гладких действительнозначных функций, определенных на интервале (-1,1) относительно нормы

\\у\\ч = {Е^то} , (6)

иеНо )

где

2 Г1

Ш = Ч <1(г)у(г)и1(г)(1т, I Е N0,

коэффициенты Фурье функции у по системе полиномов { В пространстве Щ также

определим скалярное произведение

<Лд)щ = Т.(1 + 1)23?Ш)> еЩ.

le.No

С введенным скалярным произведением пространство Н^ становится гильбертовым пространством, а норма (6) выражается через это скалярное произведение

\н- = л <У,У)н-, У е щ.

Всюду в дальнейшем будем полагать выполненным условие 5 > 1/2.; при котором (см., напр., [13]) пространство Щ вкладывается в пространство непрерывных функций, а пространство Щ+1 вкладывается в пространство функций, первая производная которых непрерывна.

Зафиксируем п е N0 и обозначим

п

(Рпу)0) = ^ уЦк) & (г), ге [-1,1], (7)

к=0

интерполяционный полином Лагранжа функции у е по узлам

-(к + 1)

и = со8—--, к = 0,1,..., п. (8)

п + 2

Здесь

а= —г, к = 0,1,...,п, ге [-1,1],

- фундаментальные полиномы, соответствующие узлам (8), В [14] автором доказана ограниченность нормы оператора Лагранжа в паре пространств Соболева ( з > 1/2. Следующая лемма устанавливает ограниченность нормы оператора Лагранжа Рп в паре пространств Соболева (Щ,Щ), з > 1/2.

Лемма 3. Для любых п е ее К, 5 > 1/2, справедлива оценка

<у/1 + ((2 з),

где ^(Ь) = 3-г - дзета-функция Римана, ограниченная и убывающая, при Ь> 1. Доказательство. Возьмем произвольную функцию у е Нд, з > 1/2. Используя равенство

Iип+1(1) - (п + 2)Тп+2(1)

1 - г 2

и известное (см., напр., [12]) соотношение

1 Г1 <Ат)и,,+1(т)Лт = п е N0, ( е [-1,1], (10)

КЗ-1 Г -г

вычислим коэффициенты Фурье полинома (7), При 0 ^ I ^ п получим

( Р5)(1) = - [1 <1(г)(РпУ)(т)иг(т)<1т = - /1 Ч(т) £ у(гк)(7п+1^Т](!\ ъ]-1 -к}-1 (т- и )к+1(1к)

= Шп+1(Ч - (п + 2)1п+2(Ч, п е N0, * е [-1,1], (9)

у уи)(2 I'1 д(г)ип+1(г)(и1 (т) -и( fc ))<!т + 2Т ) Г1 д(т)ип+1(т)<1т \

+-т к Ч-1 ^-Ш+Ж))

Так как полином ( Ui(r) — Ui(tk))/(т — tk), т G [—1,1], является полиномом степени I — 1 < п +1, то в силу ортогональности

i1 т) — Ui(t k= 0, k = 0,1, ...л,

7-i (— tkwi+itk) , , , ,

Теперь, используя соотношения (9) и (10), найдем

. 2 к(к+1)

1 f1 Я(тЫУг =sin п + 2 ,, = 0,1,...,п.

к J-1 (т — tk)U'n+i(tk) п + 2 ' ' ' '

окончательно получим

( РпУ )(/) = Е y(tk )sln2 ^П+^Ui (tk), 0 (11)

п + 2 п + 2

k=0

п <

k, = 0, 1, ... , п

Un+1, имеем

( p-)( П 2 Г1 ( )( p )( )U ( w 2 " r1 q(r)Un+i(T)Ui(T)dT

( Pny)(l) = - q(T)(Pny)(T)Ui(T)dr = - J>y(tkW —;--—T—

К J-1 К J-1 (r — tk)Un+1(tk)

2 A (. ) [1 g(r)(Un+1(r) — Un+1(tk))Ui(r)dr < . = ~l^y(t k П -+-T~\tTi—ТЛ\-,п<1.

К k=0 J-1 (T — fk)Un+1(tk)

А так как ( Un+1(r) — Un+1(tk))/(r — tk), к = 0,1,..., п, полиномы степени меньше I, то в силу ортогональности

г1 q(t)( Un+1(r) — Un+1(t k ))Ui(r)dr =0к = 01

J-1 (r — tk)Un +1(tk) , , ,..., ,

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

( рПу)(0 = 0, п<1. (12)

По определению нормы в пространстве Н с учетом равенств (12)

1 I Pny I I Н. = Е (¿ + 1)2s( Pny № (13)

Теперь вычислим коэффициенты ( Pny)(l) только для 1,0 ^ I ^ п. В (11) заменим 0 ^ I ^ п, найдем

( Pny)(0 = ^ £ Е y(^)sln2 ^^k^U^^k)Ui(tk)

k=0 meNo

2 v^ n(m +1)(k + 1) + 1)(k + 1) (m)--> sin —----sin

E ^(m)

п + 2 п + 2 п + 2

meNo k=0

1 ^ ^ ( K(m — l)(k + 1) K(m + / + 2)(k + 1)\

— > y(m) у cos--cos-

+ п + 2 п + 2

к—П N '

1 E ^(m—1)

п + 2

meN

( . (2п + 3)^(m — (/ + 1)) sin

_2(п + 2)

. K(m — (/ + 1))

sin-:-:-

2(п + 2)

. (2п + 3)п(т +(1 + 1)) \ sm . .

2(п + 2)

. п(т +(1 + 1)) sm ■

0 <1<п.

2(п + 2) /

Представляя числители выражений в последних скобках в виде

. (2п±3^тт-1±±)1 • ( dMiw ж(т—(l + i))

sm-———^-= sm п(т — (1 + 1)) cos

2(п + 2)

, ,, . п(т — (I + 1)) . . .. — cos п(т — (I + 1)) sin-—-—-= sin п(т — (/ + 1)) cos

2(п + 2)

п(т — (1 + 1))

2(п + 2)

+ (—1Г- sin П(т - (1 + ')) 1 ; 2(п + 2)

(2п + 3)п(т + (1 + 1)) sin —

2(п + 2)

т Е No, 0 ^п,

2(п + 2)

sinn^ + (I + 1)) cos

п(т + (1 + 1))

/ /7 • п(т + (I + 1)) , ,, — cosn^ + (I + 1)) sin-—-—-= sin^(т + (I + 1)) cos

2(п + 2)

п(т + (l + 1))

получим

2(п + 2)

+ (—1)т- sin П(т +(1\1)).

1 ; 2(п + 2)

(

2(п + 2)

те N0, 0 ^п,

(РпУ )(1) =

2(п + 2)

Y1 У(т — 1)

meN

sinn^ — (1 + 1)) п(т — (1 + 1)) cos

sin

п(т — (1 + 1)) 2(п + 2)

2(п + 2) \

sinn^ +(1 + 1)) п(т +(1 + 1)) cos

. п(т +(1 + 1))

sin -;-г-

2(п + 2)

2(п + 2)

, 0

Так как

sinn^ — (1 + 1)) п(т — (1 + 1)) cos

sin

!

п(т — (I + 1)) 2(п + 2)

2(п + 2)

(—1)^2(п + 2), т — (1 + 1) = 2(п + 2)fi, 0, т — (l + 1) = 2(п + 2)ц,

sinn^ +(1 + 1)) п(т +(1 + 1)) cos

V Е No,

{

. п(т +(1 + 1))

sin-;-г-

2(п + 2)

2(п + 2)

(—1Г2(п + 2), т +(1 + 1) = 2(п + 2)ц, N 0, т +(1 + 1) = 2(п + 2)ii, V Е iN

то

( ?Пут = Е (—1УУ(2(п + 2) V + l) + Y.—T-1 У(2(п + 2)V — I — 2)), 0

^eNo

Возвращаясь к (13), найдем \\РпУ\\% =

^eN

Е >2- (Е

1<г<п+1 VeNo

П У(—1Ту(2(п + 2)^ + 1 — 1)

1

+ Y,(—1)"-1У(2(п + 2)а — I — 1)\ ^ 2 £ 12s( £ (—1)"у(2(п + 2)а + 1 — Л

"eN J 1<г<п+1 V"eNo /

+2 Е Р fa—1)"-1 У(2(п + 2)V — I — 1))

1<г<п+1 V^eN )

2

_ 13(2(п + 2)а + 1 — 1)s V

2 Е (Е1>" j^ + ^+i — iy тп+ +l — 1>)

ЫЫп+1 \"eNo —!—'-i—-'- /

+2 т. [п—1т-1 **п1))

ЫКп+1 \"eN \ \ ' jh" ! у

У ( у -—-27 У (2(п + 2) v + 1 — 1)2stf (2(п + 2)а + 1 — 1))

k+iKko (2(п + 2)V + 1 — 1)2S ¿No ----1 7

E fe (2(п + 2)1 - / — 1)2- ^(2(п + 2)V — г — 1)2 *^2(2(п + 2)V — * — 1)) iMu+1 \"eN K K ' "eN /

2

^ \Ы\?т2 max > -

\ЫМп+1"^0 (2(п + 2)V + 1 —

+ max

1 <Kn +1 —/

"eNo

l2s

(2(п + 2)а — I — 1)2°

Оценим no отдельности максимумы сумм в последнем выражении. Для первой суммы имеем

12 s

max > -(14)

"5o (2(п + 2)а + I — 1)2'

(Ш+Г)" max V (2а + -Ы)'

\п + 2 1<г<п+1 ^ V п + 2 v 7 "eNo v-!—7

«1+а +^ Г«^

"eN "eN

Для второй суммы найдем

12 s

(2(п + 2)V — I — 1Г (15)

2

п + 1 2 + 1 -2

^ - max >[ 2а--^ > (2а — 1)

\п + 2 п + 2) J

"eN "eN

Подставляя (14) и (15) в (13), получим

\\РпУ\\2щ ^ (1 + С(2S))\М\^, s> 2. Лемма 3 доказана, □

Обозначим Еп(у)^ наилучшее приближение по норме пространства Н^ функции у е Н^ алгебраическими полиномами. Известно, что наилучшее приближение функции в гильбертовом пространстве доставляет отрезок ее ряда Фурье

Еп(.у)1 = 11 у-Япу\\щ, шт= Е ти(*), ге (-1,1). Следствие 1. Для, любой функции у е в > 1/2, и любого п е н0 верна оценка

- РпУ\\т < (1 + у/1 + <;(2з))Еп(у)

Доказательство. Возьмем произвольные функцию у е в > 1/2, и число п е н0. Цепочка неравенств

\\У - РпУ\\т < \\У - ЯпУ\\т + \\ЯпУ - РпУ\\н° < Еп(у)1

+ \\ Рп\\н^т \\У-ЯпУ\\нц < (1 + у/1 + <(2*))Еп(у)

1 х( ) х

-1 /1 - г2

X г1

У-1 /1 -

доказывает требуемую оценку, □

4. Анализ разрешимости Задачу (1), (2) запишем в виде операторного уравнения

Кх = Их + Ух = у, К : X ^ У,

X = {х еЩ+ | /Щ = 0}, У = Щ.

^ , *е (-1,1).

- т2(т - *)' 1 ' ;

Теорема 1. Для, всех IX 1< 1, задача, (1), (2) однозначно разрешима при любой правой части, у еУ, и верна, оценка

\\К-1\\у^х ^ (1- I X I)-1.

Доказательство. Покажем вначале, что оператор И : X ^ У обратим и шн х^-у = \ \ И \ \ -у ^х = 1- Действительно, возьмем произвольные функции х е X и у еУ и запишем их в виде их рядов Фурье в соответствующих пространствах

х(*) = Ех()т tе [-1,1], у(1) = Еу(1)Ш, 1 е (-1,1). ген гено

Уравнение

Бх = у, И : X ^ У, (16)

в этом случае будет иметь вид бесконечной системы уравнений

1х(1)иг-1(г) = у(/ - ВД-1(*), 1е н, 1е (-1,1),

а его решением будет функция

х(1) = Е 1-1у(1 - 1)Тг(г), ге [-1,1]. ген

Из произвольности выбора элемента у е У следует обратимость оператора И : X ^ У. Вычислим норму операторов И : X ^ У и И-1 : У ^ X. Для произвольного элемента х е X имеем

\ \ вх\\у = ^ +1)25((/ + 1)х(1 +1))2 = Е г2(*+1)х2(0 = \ \ х|| X,

гено ген

для произвольного элемента у е У найдем

\\о-1у \\х = Е12{з+1)(1 -1у(1 -1))2 = ВI + 1)2зу2(1) = М2.

Это и означает, что \\Б\\х= \\Б-1\\у^х = 1-

Найдем теперь норму оператора V : X ^ У. Вновь возьмем произвольный элемент х е X

х(г) = ^Хс(1)т1(г), г е [-1,1],

1&Ч

и применим к нему оператор V. Так как (ем, напр., [12])

1 [1 = и-1а), I е N [-1,1],

то

( Vх)(t) = -/1 = А^Х(1)и-1(1), г е (-1,1),

и норма функции Vх в пространстве У оценивается неравенством

\\VxW2 = -2 Е(I + 1)2*хс2(1 + 1) ^ А^12^+1)хс2(1) = \2\\х\\2Х,

1&Чо 1&1

поэтому \\У \\ - По лемме 1 оператор К = Б + V : X ^ У обратим для всех А,

| - | < 1

\\К-1\^х ^ (1-1 А |)-1, | А |< 1. Теорема 1 доказана, □

5. Метод Галеркина

Зафиксируем п е N0. Приближенное решение задачи (1), (2) будем искать в виде отрезка ряда Фурье

хп+1^)= Е хсп+1(1)Тг(1), 1е [-1,1], (17)

неизвестные коэффициенты хп+1(1), 1 = 1, 2,..., п+1.; которого найдем по методу Галеркина из системы уравнений

1хсп+1(1) + Ахп(1) = у(1 - 1), 1 ^ I ^ п + 1, (18)

где

2 Г1

У(1) = ~ Я(г)у(т)и1(т)(1т, 0 ^ I ^ п, - -1

коэффициенты Фурье функции у по системе полиномов {и{}1 ем0.

Теорема 2. Для любых фиксированных А е К | А |< 1 и п е N0 система, уравнений (18) метода Галеркина решения задачи (1), (2) имеет единственное решение

Гп+1(1) = (1 + А)-1у(1 - 1), 1 ^ I ^ п +1,

и приближенные решения

х*п+1^)= £ х*п+1(1)т), ге [-1,1],

сходятся, к точному решению х* задачи, (1), (2) со скоростью

\\х* -х*п\\х ^ (1- | А \)-1Еп(уХ.

Доказательство. Обозначим

Хга = spanjTiYn = span{Ui}?=0

подпространства пространств X и YeooTBeTCTBenHO, Систему уравнений (18) запишем в виде операторного уравнения

KnXn+1 = Qn(DXn+1 + Vxn+1) = QnV, Kn :Xn ^ Yn, (19)

По теореме 1 оператор K в условиях теоремы 2 обратим, К роме того, dim Xn = dim Yn и K — Kn = 0 na Xn, n E No- Поэтому по лемме 2 операторное уравнений (19) имеет единственное решение

x*n+1W = £ (1 + Л)-1^(г — 1)TiW' te [—!, !],

1<i<n+1

при любой правой части Qny E Ynа, и погрешность приближенного решения оценивается неравенством

I I X* — xn+11 |x ^ |\к-1|\Y^X\\y — Qnyh ^ (1— | Л l)-1En(y)sq. Теорема 2 доказана, □

При всех преимуществах метода Галеркина он имеет один существенный недостаток -неконструктивность. Действительно, для вычисления коэффициентов Фурье правой части уравнения (1) необходимо брать интегралы, которые не от всех функций вычисляются в явном виде. Этого недостатка лишен метод коллокаций, который, однако, в ряде пространств, например, гельдеровских или непрерывных функций, имеет худшую скорость сходимости чем метод Галеркина,

В следующем разделе показано, что в пространствах Соболева скорость сходимости метода коллокаций не уступает скорости сходимости метода Галеркина,

6. Метод коллокаций

Вновь зафиксируем n E N. Приближенное решение задачи (1), (2) как и по методу Галеркина будем искать в виде отрезка ряда Фурье (17), но его неизвестные коэффициенты {Xn+ 1(/)}t+11 найдем теперь по методу коллокаций из системы уравнений

(Dxn+1)(tk) + ( Vxn+1)(4) = y(tк),k = 0,1,...,n, (20)

по узлам (8)

Обозначая w — Kxn+1 — у, можно метод Галеркина (18) записать в виде системы уравнений

2 Í'1

- q(r)w(T)Ui(T)dr = 0, 1 = 0,1,...,n, (21)

я J-1

а метод коллокаций (20) - в виде системы уравнений

w(tk) = 0, k = 0,1,..., n. (22)

Аппроксимируем интегралы (21) интерполяционными квадратурными суммами

2 Г q(r)( Pnw)(r)Ui(r)dr = —¿w(t k )Ui (t k )sin2 *(k + 1), l = 0,1,...,n, я J-1 n + 2 k=0 n + 2

и обозначим

2 [1 2 n я(к + 1)

П = - q(T)w(T)Ui(r)dT---Vw(t k )Ui(t k )sin2 ( + J , 1 = 0,1,...,n,

я J-1 n + 2 k=0 n + 2

остаточные члены этих квадратурных сумм. Из чисел {г}^ образуем полином

п

(Кпт)Ц) = ЕпигЦ), 1& [-1,1].

1=о

Теперь запишем метод Галеркина (18) для искомой функции т — Япт

2 Г1

к

q(r)(w - Rnw)(r)Ul(r)dr = 0, 1 = 0,1,..., п, (23)

'-1

Систему уравнений (23) назовем модифицированным методом Галеркина-Петрова для задачи (1), (2),

Лемма 4. Метод коллокаций (20) и модифицированный метод Галеркина-Петрова (23) эквивалентны в том, смысле, что равенства, (22) выполняются тогда, и только тогда, когда, выполняются, равенства, (23).

Доказательство. Равенства (23) представим в виде

2 Г1 2 Г1

- д(т)(т — Кпт)(т)и1(т)(1т = - д(т)т(т)и1 (т)(т — п

2 ^ , ч • 2^(к + 1)

Twitк)Ut(tк) sin2 ( +1 =0, 1 = 0,1, ...,п.

f J г) -U /

п + 2^ п + 2

к=0

Теперь равенства (23) следуют из равенств (22) тривиально.

Пусть выполнены равенства (23), Матрица ( Ui(tk))]l^o невырождена, поэтому однородная система уравнений

Y^w(tкmtк) sin2 =0, 1 = 0,1, ...,п.

к=0

имеет только нулевое решение

2^(к+1)

Так как

w( tk) sin -= 0, к = 0,1,...,п.

Vk) п + 2 ' ' ' '

sin2 п(к + 1 =0,к = 0,1,...,п, то w(tk) = 0, к = 0,1,...,п. п + 2

Лемма 4 доказана, □

Лемма 5. Для, любой функции т & Щ и любого п & N верна, оценка

\\Ппт\\щ ^ у/1 + фз)Еп(т)а9.

Доказательство. Зафиксируем функцию т & Н^ и число п & М0. Коэффициенты Г[, 1 = 0,1,.., п, это первые п + 1 коэффициенты Фурье функции т — Рпт, поэтому по лемме 3 имеем

\\Rnw\\Hs = \\Qn(w - Pnw)\\H'q = \\Pn(Qnw - w)||H| w 1 + C(2s)En(w)¡. Лемма 5 доказана, □

Теорема 3. Для любых фиксированных А Е Е, | А |< 1 и п Е N0 система уравнений полиномиального метода коллокаций (20) имеет единственное решение и

приближенные решения

<+1(*)= Е ®Пц(0т 1Е [-1, 1],

сходятся, к точному решению х* задачи, (1), (2) со скоростью

I I х* — хПЦх ^ (1- | А \)-1Еп(у)1.

Доказательство. Систему уравнений (20) полиномиального метода коллокаций запишем, следуя леммме 4, в виде системы уравнений (23) модифицированного метода Галеркина-Петрова, В операторной форме система уравнений (23) будет иметь вид Qnw = QRnw. Делая обратную замену w = Кхп+1 — у, получим уравнение

QKхn+l = Qn(y + Rnw) (24)

метода Галеркина для уравнения

Бх + Ух = у + Кп IV.

По лемме 2 оператор Кп = QnK обратим в паре пространств (Хп, Уп), и погрешность приближенного решения хП+1 уравнения (24) по методу Галеркина оценивается неравенством

I I х* — хП+1 I\х ^ (1— | А |)-1\\у—КаШ — QnУ + QnRnW|\у

< (1— I А |)-1(\\у — Qny||Y + \\КпШ — QnRnW|\у). Так как Кп = Qn — QnPn, то КПш — QnRnw = 0, поэтому

I I х* — хП+1\\х ^ (1— | а |)-1Еп(^.

Теорема 3 доказана, □

7. Замечание

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

СПИСОК ЛИТЕРАТУРЫ

1. D. N. Arnold, W. L. Wendland On the asymptotic convergence of collocation methods // Math, of Сотр. 41, № 164. 1983. P. 349-381.

2. D. N. Arnold, W. L. Wendland The convergence of spline collocation method for strongly elliptic equations on curves // Number. Math. 46. 1985. P. 317-341.

3. J. Elshner On spline collocation for singular integral equations on an interval // Semin. Anal. Oper. Equat. and Numer. Anal. 1985/86. 1986. P. 31-54.

4. J. Elshner On spline approximation for a class of integral equations III // Semin. Anal. Oper. Equat. and Numer. Anal. 1986/87. 1987. P. 25-40.

5. J. Elshner On spline approximation for singular integral equations on an interval // Preprint Akad. Wiss. DDR, Karl Weierstrass Inst. Math. 4. 1987.

6. J. Elshner On spline approximation for singular integral equationson an interval // Math. Nachr. 139 1988. P. 309-319.

7. S. B. Prossdorf, G. Schmidt,- A finite element collocation method for singular integral equations 11 Math. Nachr. 100 1981. P. 33-60.

8. S. B. Prossdorf Recent results in numerical analysis for singular integral equations // Proc. 9th Conf. Probl. and Meth. Math. Phvs. (TPM) Karl-Marx-Stadt, June 27 - July 1. 1988. P. 224-234.

9. A. I. Fedotov On the asymptotic convergence of the polynomial collocation method for singular integral equations and periodic pseudodifferential equations// Archivum Mathematicum. Tomus 38, № 1. 2002. P. 1-13.

10. Красносельский M. А., Вайникко Г. \!.. Забрейко П. П., Рутицкий Я. Б., Стеценко В. Я. Приближенное решение операторных уравнений, М.: Наука. 1969. 418 с.

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

12. Градштейн И. С., Рыжик И. М. Таблицы интегралов, сумм, рядов и произведений, М.: Наука. 1971. 1108 с.

13. Тейлор М. Псевдодифференциалъные операторы, М.: Мир. 1984. 472 с.

14. Федотов А. И. Оценка нормы интерполционного оператора, Лагранжа в многомерном пространстве Соболева с весом // Математические заметки. 99, № 5. 2016. С. 752-763.

Александр Иванович Федотов,

Казанский филиал московского социально-гуманитарного института, ул. Столярова, 3, 420030, г. Казань, Россия E-mail: [email protected]