Научная статья на тему 'Поперечники некоторых множеств дифференцируемых функций'

Поперечники некоторых множеств дифференцируемых функций Текст научной статьи по специальности «Математика»

CC BY
160
31
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
АППРОКСИМАЦИЯ / ЛОКАЛЬНЫЕ СПЛАЙНЫ / ПОПЕРЕЧНИКИ БАБЕНКО И КОЛМОГОРОВА / APPROXIMATION / LOCAL SPLINES / BABENKO AND KOLMOGOROV WIDTHS

Аннотация научной статьи по математике, автор научной работы — Бойков Илья Владимирович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Бойков Илья Владимирович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Поперечники некоторых множеств дифференцируемых функций»

УДК 518.5

И. В. Бойков

ПОПЕРЕЧНИКИ НЕКОТОРЫХ МНОЖЕСТВ ДИФФЕРЕНЦИРУЕМЫХ ФУНКЦИЙ

Аннотация. Вычислены поперечники Бабенко и Колмогорова функциональных множеств Quy (Q,M) и Q"y (AM), где Q = [-1,1]1, l = 1,2,..., r и и -натуральные числа, у - действительное неотрицательное число. Построены локальные сплайны, которые являются оптимальным методом аппроксимации функциональных множеств Q-у (q, m ) и Q-у (Q, M).

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

Abstract. Evaluated Babenko and Kolmogorov widths of functional sets QUy(Q,M) and Q.у(Q,M), where Q = [-1,1]1, l = 1,2,..., r and и are natural numbers, у is a real nonnegative number. Constructed local splines which are the optimal in order methods for approximation of functional classes Q.y (Q, M) and QU,y (Q, M).

Keywords: approximation, local splines, Babenko and Kolmogorov widths.

Введение

Класс функций Qr (Q,M) был введен К. И. Бабенко, и им же была сформулирована задача вычисления поперечников Бабенко и Колмогорова на этом классе функций [1]. Эта задача была решена автором [2]. Им же были введены классы функций Qr у (Q,M), Qryp (Q,M), Br у (Q,M) и вычислены

поперечники Бабенко и Колмогорова на этих классах функций [3, 4]. В работах [2-4] были также построены локальные сплайны вычисления функций из классов Qr, Qrу, ВГу. Показано, что эти сплайны являются

наилучшим по порядку по точности методом приближения функций из классов Qr, Qr^, Br ,у.

Интерес к наилучшей аппроксимации классов функций Qr, Qr у, Br у

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

В данной работе результаты [3, 4] распространяются на более широкие

классы функций - классы функций Q-у(Q,M), Q-у(Q,M), Q = [-1,1]1, l =1, 2,... Для этих классов функций вычислены поперечники Бабенко и Колмогорова и построены локальные сплайны.

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

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

Определение 1 [1]. Пусть L - множество n -мерных линейных подпространств пространства В. Выражение

dn (X, В) = шГ sup inf llx - -I,

Ln xeX ueLn

где последний inf берется по всем подпространствам Ln размерности n, определяет n -поперечник Колмогорова.

Определение 2 [1]. Пусть % - множество всех n -мерных линейных подпространств пространства В, Map(X, %) - совокупность всех

непрерывных отображений вида П : X ^ X, где X е%. Выражение

dn'(X,В) = inf sup \\x -П(x)||,

(Ln ,П) xeX

где inf берется по всевозможным парам (L, П), состоящим из n -мерного

линейного пространства L с В и непрерывного отображения П : X ^ L, определяет линейный n -поперечник Колмогорова.

Определение 3 [1]. Пусть %e Rn. Выражение

5n (X)= inf sup diam П-1П(x),

(П-.X ^ Rn ) xeX

где inf берется по всем непрерывным отображениям П : X ^ Rn, определяет n -поперечник Бабенко.

Приведем определения классов Qr, Q^, QUу, QUу.

В работе К. И. Бабенко [1] введен класс функций Qr (Q, М).

Определение 4. Пусть Q = [-1,1]1,l = 1,2,... Функция ф(Х1,..., xi) принадлежит классу Qr (Q, M), если выполнены условия

max

xsQ

< M

Э^ф(Х1,..., xi)ГдХ.1... dxJ

xsQ при 0 <| v |< r;

Э|v|ф(x1,..., xi )ldxV... dx^1 < M/(d(x, T))|v|-r, x eQ \ Г,

при r <| v |< 2r +1, где x = (xb..., xi), v = (v1,..., vi), | v |= V1 + ... + vi, d(x,Г) -расстояние от точки x до границы Г области Q, вычисляемое по формуле

d(x,Г) = min min(| -1 - xt |, |1 - x,■ |).

1<i<i

Приводимые ниже классы функций являются обобщениями Qr (Q, M).

Определение 5 [3]. Пусть Q = [-1,1/, i = 1,2,... Функция ф (x1,..., xi) принадлежит классу Qr у (Q,M), если выполнены условия

тах

хеО

Э^ф( х)/Эх^... Эхр

< М

при 0 <| V |< г;

Э^ф(х)/Эх^...Эх^ <М/((х,Г)) г С , хеО\]

при Г <| V |< 5, где 5 = г + [у] + 1, У = [у] + Ц, 0< Ц <1, С = 1 “Ц при у нецелом, 5 = г + у при у целом.

Определение 6. Пусть й = [-1,1]1, I = 1,2,..., у и и -

неотрицательные целые числа. Множество (й,М) состоит из функций

/(^, ^2,..., ^), удовлетворяющих условиям

тах

геО

< М, 0 <| V |< г -1;

Э^ф(г )/Эг11... Эг^

Э^ф(г )/Эг[1... Эц1 < М (1+ |іп иб (г, Г)|), г е О \ Г,| V |= г; Эмф(г )/Эг[1... Эц1 < М (1+1 іп и-1б (г, Г)|)/ ((г, Г) )|-г, г е О \ Г, г <| V |< 5,

где 5 = г + у.

Наряду с классом функций Qrу (О, М) введем класс функций

б“у (О, м ).

Определение 7. Пусть й = [-1,1]1, I = 1,2,..., и - натуральное число, у - нецелое число. Класс ОТ-у (й, М) состоит из функций, удовлетворяющих следующим условиям:

тах

геО

Э^ф(г )/Эг^... Эг;

< М, 0 <| V |< г,

Э^ф(г )/Эг^... Эг;'

<

М

~г -С

(1+|1п иб (г, Г) |), г <| V |< 5,

где 5 = г + [у] + ц, у = [у] + ц, 0 < ц < 1, С = 1 -Ц.

г / (с) *

Через Тг (/,[а,Ь], с) = ^ ——— (г - с) обозначим отрезок ряда Тейлора

к=0

функции /(г), определенный в области [а,Ь] по степеням (г - с), с е [а,Ь].

1 Поперечники функциональных множеств QUу, QUу, О = [-1,1]

Теорема 1. Пусть О = [-1,1/ , г, и, у - положительные целые числа, 5 = г + у. Справедлива оценка

8И ^,у (О,М)) ^ бп ^“у (О,М), С) ^ п-5.

Доказательство. Вначале оценим снизу величину 8п^^у(О,М)). С этой целью заметим, что класс функций Qrу (О,М) вложен в класс

функций QUу(О,М). Для класса функций Qrу(О,М) известна [3, 4] оценка снизу поперечника Бабенко:

8п ^г ,у (О, М)) > Ап-5.

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

8п (ёи,у (О, М)) > 8п (<2Г ,у (О, М)) > Ап-5. (1)

Оценка снизу получена.

Построим локальные сплайны, аппроксимирующие функции из класса

QUУ (О, М) с точностью 0(п-5).

В отдельности опишем построение в случае и = 1 и и > 2.

Пусть и = 1. Разделим сегмент [-1,1] на 2N частей точками

г* = -1 + (k/N)V, і* = 1 - (*/Ю”, к = 0,1,..., N, V = 5/(5 - у). Введем обозначения: а* =[гк,гк+1] А2 =[ік+1,тк] к = 0,1,...,Ж-1.

1 2

Сегменты А0, А0 разделим на более мелкие сегменты А0,у = [г0,у,г0,у+1] А0,у =[т0,у+1,т0,у], І = 0,1,---,М0 -1 где

г0, І = г0 + (г1 - г0) у/М0, т0, І = т0 - (і0 -11) у/М0, у = 0,1,..., М0, М0=[1пЖ ].

Построим интерполяционные полиномы таким образом, чтобы в число узлов интерполяции входили концы сегментов А0 і , І = 1,2,

і = 0,1,..., М0 -1, а* , І = 1,2, к = 0,1,..., N -1.

Полином Р5 (/,[а,Ь]), интерполирующий функцию /(г) е QUу (О,М) на сегменте [а, Ь], строится следующим образом.

Обозначим через С к, * = 1,2,..., 5 узлы полинома Чебышева первого рода степени 5. Отобразим сегмент [С1, С5 ] с [-1,1] на сегмент [а,Ь] таким образом, чтобы точки С1 и С5 перешли в точки а и Ь соответственно. Образы точек С1, . ., С5 на сегменте [а, Ь] обозначим через С1,..., С'5.

Интерполяционный полином, построенный по узлам С1, . . ., С5, обозначим через Р5 (/,[а, Ь]).

На сегментах Аду, І = 1,2, у = 0,1,...,М0 -1, А*, І = 1,2, к = 0,1,..., N -1, функция / (г) е ОТ-у (О, М) аппроксимируется интерполяционными полиномами Р5 (/, А0 у), Р5 (/, А*), І = 1,2, у = 0,1,..., М0 -1, к = 0,1,..., N -1. Сплайн, составленный из этих полиномов, обозначим через

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Ь (г).

Можно показать, что

II/(') - Л V)| 1с([-Ц], < А»-5.

Отсюда следует, что

^ (О., у (й, М), С) < Ап-5. (2)

Используя неравенство [5] 8п (X, С) < 2dn (X, С) и оценки (1) и (2), завершаем доказательство теоремы при и = 1.

Разделим сегмент [-1,1] на 2N частей точками tk = -І +

N

т— =1-1 — I , k = 0, І,..., N, v = s/( s -у). Обозначим через Д— и д)2

k N

i2

сегменты Д— =[t—, t—+і] и Д — =[т—+1, т— ], к = 0,1,..., N -1. Введем числа

M0=[lnu/rN] и M— = [ln(u-i)/s(N/k)], к = 1,2,... N. В случае, если

ln(u-i)/s( N/k) й 1, полагаем M— =1. Каждый из сегментов Д— и Д^ разделим на M— равных частей, к = 0,1,..., N -1. Сегменты, полученные в результате

деления, обозначим через Д— j, і = 1,2, к = 0,1,..., N -1, і = 0,1,..., к -1.

В каждом сегменте Д— ^ функцию f (t) будем приближать интерполяционным полиномом Ps (f, Д— і), і = 1,2, к = 0,1,..., N - І, і = 0,1,...,M— -1. Сплайн, составленный из полиномов Ps(f,Д— ^), і = 1, 2, к = 0,1,..., N -1, і = 0,1,..., M— -1, обозначим через f^ (t).

Оценим точность аппроксимации функции f (t) сплайном f^ (t).

Пусть t є Д|0, і = 1,2. Очевидно,

\\f (t) - fN (t c (д0 0) й AEs-1 (f, Д0,0 ) s,

где Д0 0 =[t0 0, t0 i], Es (f, [a,b]) - наилучшее приближение функции f(t) на сегменте [a, b] полиномом s -го порядка, Аs - константа Лебега.

Используя отрезок ряда Тейлора Tr-i(f, Д0 0,-1), покажем, что

Es-i(f, Д0,0) й f (t) -Tr-i(f, Д0,0,-1)

C (Д0,0)

й

------— max

(r - i)! гєДІ

0,0

j f(r )(T)(t -T)r 1 d т

-1

M

-----г— max

(r - i)! /єлІ

0,0

j (i+|lnu(i + т) |)(t -т)г-i dт

-1

й

< B(hi0 | lnuh00 I) < B

NvM0

NvM0

<

< В----1---(1п иЫ + 1пи 1п N ) = Б—,

N 1п т N

где % = Vмo, к = г1 -го-

Здесь и всюду ниже через В обозначаются константы, независящие от N и от функции /(г).

Аналогичным образом оцениваются значения Ег^(/, д) у),

у = 1,2,--.,М) -1 и Ег — (/, Д^-), у = 1,2,--.,М) -1.

Так как константа Лебега А5 < В 1п 5, то окончательно имеем

|/(0 -Л(О||с<д0) < BN-s, <= 1,2-

Оценим ||/(г) - /А> (г)|С(д^ У): / = 1,2, к = 1,2,-,,, N -1, у = 0,1,,--, к -1-Очевидно,

II/(г) - /М (г )|| с (дк])

<

B(tk,1 - tk,0)S Г N^VY

2s-is!

Г 1 + ln u -іГ N \vY < B Г hk 1

І І k, ) I kM

s'N^

N

k I 11+ln^71 <

< B

N ) І N

i

ln N

k

(и -1)/s

N

,vY

N

1 + lnu-ij I <

(k + 0)(v-i) s-vY B

< B-----------------<-----•

Ns Ns

Аналогичным образом оцениваются нормы

Ilf (t) - fN (t) С (Ai ), k = 1, 2, •••, N -1, j = 1, 2, •••, k -1,

v k ,j’

||f (t) - /n (t)||с(A2 ), k = 1, 2, •••, N -1, j = 0,1, •••, k -І Собирая полученные выше оценки, имеем

Ilf (t) - /n (t)|С([-1,1]) <BN

Оценим число узлов, используемых при построении локального сплайна^

Для этого оценим число сегментов Дк у, 1 = 1,2, к = 0,1,---, N -1,

у = 0,1,---, к -1-

Очевидно, в случае, когда и -1 < 5,

N-1

т = 2^Мк <2

к=0

Г и N-1 и- — Л 1пг N + V 1п 5 — к=1 к

г

< 2

N-1

N

\

1пг N + V 1п— к=1 к

= 2

N-1

1п г N + (N -1) 1п N - V 1п к

к=2

< 2

N -1

1п

N + (N - 1)1п N - | 1п хёх

< 2

1п г N + (N -1) 1п N - (N -1) 1п(N -1) + N

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

= 2

1пг N + 1п

N N -1

N-1

Л

+ N

= BN -

Рассмотрим общий случай- Пусть q = [(и -1)/5] + 1 - Тогда

N -1

т = 2 V Мк < 2

к=0

' и N-1 и-1 Л Г N-1 Л

1пг N + V 1п

N

к=1

< В

= В

Г N-1

—+ V

V к=2

1п q

N

Л Г N

< В

N

/

Л Г

ёх

= В

N + V 1п q

V к=1 1

N

Л

< BN -

Следовательно, общее число узлов, используемых при построении локального сплайна /— (г), равно п = 5т = BN -Отсюда следует, что

||/(г) - А (г)|С([.ц], < ^-

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

ёп (еи,у (й, М), С) < Ап-5 -

(3)

Используя неравенство 5п (X) < 2ёп (X, С) и оценки (1), (3), завершаем доказательство теоремы в общем случае-

Аналогичным образом доказывается следующее утверждение -Теорема 2- Пусть й = [-1,1], г, и - положительные целые числа, у -положительное рациональное число, 5 = г + [у ] + 1 - Справедлива оценка

5п (0“ (й,М)) ^ ёп (ви,у (й,М), С) ^ п-5 -

У

2 Поперечники множества функций 0й у (й,М), й = [-1,1] , I = 2, 3,...

Теорема 3. Пусть й = [-1,1]1, I = 2,3,..., и = 1,2,..., V = 5/(5-у). Справедливы оценки

5n (QU,y (й,M)) П dn (QU,y (й,M), С)

п

n- s/l, v < l/(l - 1),

l Л s/l

I , v = l/(l -1), r > ul, s > (u - i)L n )

Доказательство. Вначале найдем оценку снизу поперечника Бабенко 5И (0йу(й,М)). Для этого заметим, что множество функций Qrу(й,М)

вложено в Qu y (й, M )• Известно [3, 4], что

5n (Qr ,y (й, M)) >

n-s/l, v < l/(l -1), s/l

, v = l/(l - !)•

ln n

n

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

5n (QU,y (й, M)) > 5n (Qr ,y (й, M)) > A

n-s/l, v < l/(l -1),

ln n

s/l

, v = l/(l - !)•

Оценка снизу получена-

Приступим к построению непрерывных локальных сплайнов, аппроксимирующих функции /(г) е ОУ.у (й,М), г = (1,..., г/ ).

Вначале рассмотрим случай, когда V < //(/ -1).

Обозначим через Дк множество точек г = (г1,..., г/) таких, что расстояние от г до границы Г области й удовлетворяет неравенству

—■j < d (t, Г) < І I , k = 0,1, •••, N -1, v = s/(s -y )•

Здесь d(t,Г) = min min(| ti -1|,| tt +1|)

1<i<l

Каждую область Ak покроем кубами и параллелепипедами Aik с гранями, параллельными координатным плоскостям, и с ребрами, длины

ii,•••, il

которых равны hk =

k +1 !v Г k^v

N

N

k = 0, i,_, N -І При этом в число

І •••, il

^-i «1,-

k = 0,i,_ N - L

вершин кубов (или параллелепипедов) Д. входят вершины кубов (или

1---, ‘I

параллелепипедов) Дк

То обстоятельство, что среди кубов Дк і при каждом значении к

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

М0=1 +

(1п N )и/г + = М

Nї(и -1)/5

1п—

к

. Каждый из квадратов Дк ■

11,..., 11

покроем Мк квадратами, которые обозначим через Дк ■ .

1..., ■.

I; ..., ]/

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

Обозначим через Р* * (/, [а^ 61;...; а/, Ь/ ]) интерполяционный

полином, который определяется формулой

Р*, ..., 5 (/ ,[аЬ Ь1;-..; а/, Ьг ]) = Р^Ч-Л р!1 [ / ,[а/, Ь ]],[а/-ь Ьг -1]],.--,[аь 61]], где полином Р* (/, [а,Ь]) определен в предыдущем параграфе, а верхний

индекс в выражении Рг (/, [а.,Ь. ]) определяет переменную, по которой проводится интерполяция-

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Построение сплайна начнем с куба Д— 1. В этом кубе функцию /(г) приближаем интерполяционным полиномом Р* * (/, Д— 1). Построив этот

полином, переходим к приближению функции / (г) в кубах и

N - 2

параллелепипедах Д. ■ ■ ■ . В каждом из кубов (и параллелепипедов)

‘1, ---, 7, Л, ---, ^/

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

Г ГД*"2. . . 1

-,511, ^..., l2, Л— .і J,

полиномом Р5 51 /, Д. ■ ■ ■ I, который строится следующим

образом. В узлах интерполяционного полинома Р5 5I /, Д. 2 ■ ■

’■■■’ ^ г1,...,7,д,...^і

не расположенных на гранях куба Д. 1, берутся значения функции /(ґ),

а в узлах, расположенных на соответствующих гранях куба Д. 1, берутся

значения полинома Р5 5 (/, Д. 1). Сначала строятся интерполяционные

полиномы Р5 51 /, Д. 2 ■ -їв области Д. 2, аналогичным

’■■■’ ^ гl,..., Іі, Jl,..., .і )

образом строятся интерполяционные полиномы в области Д .-3 и

последовательно во всех областях до Д0 включительно.

Сплайн, составленный из полиномов

Р*,~* С Г •<..., і;.,.... ] •к = 0-1-^.-1-

обозначим через /. (ґ1, ..., ґі).

После проведения громоздких вычислений можно показать, что

II f (t)-fN(t)iic(Q)< AN~s

и что число узлов n, используемых при построении сплайна fN (t), равно

n П Nl.

Из непрерывности локального сплайна fN (t) и двух последних оценок следует, что

dn (Q,Y (Q,M), c) < An-s/l.

Из этого неравенства, неравенства Sn (X, С) < 2dn (X, С) и оценки снизу поперечника Бабенко следует справедливость теоремы при v < l/(l -1).

Аналогичным образом исследуется случай, когда v = l/(l -1).

Теорема доказана.

Рассмотрим случай, когда и = 1, v = s/(s - у), v > l/(l -1).

В этом случае справедливо утверждение, которое приведем без доказательства.

Теорема 4. Пусть Q = [-1,1jl, l = 2,3,..., и = 1, v = s/(s- у),

v > l/(l -1). Справедлива оценка

Sn (Q.,y (Q,M)) n dn (Q1,y (Q,M), С) n (s-nnl-1}.

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

1. Бабенко, К. И. О некоторых задачах теории приближений и численного анализа / К. И. Бабенко // Успехи математических наук. - 1985. - Т. 40. - Вып. 1. -С. 3-28.

2. Бойков, И. В. Оптимальные по точности алгоритмы вычисления интегралов / И. В. Бойков // Оптимальные методы вычислений и их применение : межвузовский сборник научных трудов. - Вып. 8. - Пенза : Изд-во Пенз. политехн. ин-та, 1987. - С. 4-22.

3. Бойков, И. В. Аппроксимация некоторых классов функций локальными сплайнами / И. В. Бойков // Журнал вычислительной математики и математической физики. - 1998. - Т. 38. - № 1. - C. 25-33.

4. Бойков, И. В. Оптимальные методы приближения функций и вычисления интегралов / И. В. Бойков. - Пенза : Изд-во ПензГУ, 2007. - 236 с.

5. Теоретические основы и конструирование численных алгоритмов задач математической физики / под ред. К. И. Бабенко. - М. : Наука, 1979. - 196 с.

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

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

Boykov Ilya Vladimirovich Doctor of Science (in Mathematics), professor, head of sub-department of highest and applied mathematics, Penza State University

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

Поперечники некоторых множеств дифференцируемых функций /

И- В- Бойков // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2009. - № 1 (9). - С. 44-54.

i Надоели баннеры? Вы всегда можете отключить рекламу.