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

Математические заметки СВФУ Июль—сентябрь, 2015. Том 22, № 3

УДК 517.518.87

АСИМПТОТИЧЕСКИ ОПТИМАЛЬНЫЕ ФУНКЦИОНАЛЫ ПОГРЕШНОСТИ РЕШЕТЧАТЫХ КУБАТУРНЫХ ФОРМУЛ Ц. Ж. Юмова

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

Ключевые слова: решетчатые кубатурные формулы, периодический функционал погрешности, анизотропное пространство Соболева.

Ts. Zh. Yumova. The asymptotical optimal error functionals of lattice cubature formulas.

Abstract: With the using of functional analytical method of the theory cubature were minimized the norm of periodic error functional, the quality of the constructed formulas for the family of elements of anisotropic space had been improved.

Keywords: lattice cubature formulas, periodic error functional, anisotropic Sobolev spaces.

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

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

Основные результаты исследования задач оптимизации кубатурных формул на классах функций анизотропных пространств Инеодинаковых вдоль разных координатных направлений из-за гладкости функций, рассматривались в работах Ц. Б. Шойнжурова [2,3], М. Д. Рамазанова [4]. В частности, М. Д. Рамазанов исследовал кубатурные формулы на невесовом пространстве

© 2015 Юмова Ц. Ж.

периодических функций с единичным кубом в качестве основного периода. Однако при таком определении нормы функции возникали определенные трудности при периодическом продолжении функции на единичный куб. Ц. Б. Шой-нжуров продолжил функции из рассматриваемой области на все пространство, «избавившись» от ограничений, что позволило к периодической на всем пространстве функции применить преобразование Фурье.

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

Предварительные сведения и обозначения. Пусть к = 1,... ,п, Ьи >

0 — шаги решетки, хи — узлы формулы, Си — коэффициенты формулы, ти — гладкость функции вдоль выбранных координатных направлений, т = (тох,

т2,..., тп), т* = п/ ^ т-1^, Ь = diag(h1, ..., Ьп) — матрица периодов,

Ан = {х € Еп, 0 < хи < Ьи, к =1, 2,..., п} — фундаментальный параллелепипед с длинами ребер Ьи, А^ = Ьп = det Ь =0, А = {х € Еп, 0 < хи < 1, к = 1, 2,..., п} — фундаментальный единичный куб.

В работе исследуются кубатурные формулы на классах периодических функций анизотропного пространства с матрицей периодов Ь, у которых каждая обобщенная частная производная обладает свойством гладкости, зависящей от координатных направлений. Построенные формулы будут использованы для целого семейства подынтегральных элементов анизотропного пространства ~\¥™{ЕП) с естественной нормой

1 /р

М^Мт

+ \Отк<^(х)\р ¿х

к= 1

<.

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

N

^Си^х^)= ]ТС^(х - х(и)),<^(х) . и=1 \и=1 )

Дельта-функции имеют смысл при воздействии на непрерывные пробные функции, отсюда требование вложения основного пространства в пространство непрерывных функций ]¥™(ЕП) С С(ЕП), обеспечиваемое условием вложения

п

р — т-1 > 0 [5]. Это вложение непрерывно, т. е. функционал погрешности к=1 _ кубатурной формулы не только линеен, но и ограничен на ]¥™(ЕП). Знание

численной мажоранты для его нормы в сопряженном пространстве И7!™ (Еп)

р

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

Известно [6], что дифференциальные свойства анизотропного пространства не одинаковы по направлениям: оно полное при 1 < р < то, сепарабельно при 1 < р < то и рефлексивно и равномерно выпукло при 1 < р < то.

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

п

ных. Дифференциальный оператор Ь(П) = ^ ( — 1)ткЕ2гПк, входящий в такое

к=0

уравнение, порождается видом нормы функции в основном пространстве. Доказано (см., например, [3]), что любой финитный линейный функционал 1(х) из Ш™ представим в виде

/и

Dmku * Dmktpdx,

7--Г»

где if £ W™, и £ W™ и выполняется равенство

L(D)u=l(x), lewf. (1)

Известно [7], что фундаментальное решение £2ш(х) оператора L(D) неединственно, оно определяется с точностью до слагаемого €9,ш{х), являющегося произвольным решением однородного уравнения L(D)e= 0. Для того чтобы функция £2т(х) £ W™ (£«), удовлетворяющая уравнению

L(D)£2m(x) = 6(х), (2)

была фундаментальным решением уравнения (1), необходимо и достаточно, чтобы ее преобразование Фурье удовлетворяло уравнению

L(2mOF\e2Wi(x)} = 1,

где

ии

L(2n^) = ]T(-1)mk (2ni& )2mk + 1 = 1 + )2mk.

k=1 k=1 Фундаментальное решение оператора L(D), в дальнейшем ядро представляет собой бесконечно дифференцируемую при x = 0 суммируемую в Еи функцию и имеет вид

1

£2 т(х) = F~

E(-1)mfe |2п& | 2mk

k=0

Благодаря введению Ц. Б. Шойнжуровым [2] специальной нормы, для которой соответствующий дифференциальный оператор хорошо изучен и описан в

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

п

вложения 1 — р-1 Е т-1 > 0 были получены общие представления периодиче-к=1

ского функционала погрешности = 1 — Фо(^-1ж):

Лг

к=0

(Бтк<Ао(^-1ж) + С(0))<^(ж) ¿ж

и соответствующей ей экстремальной функции

^^ е2пгЬ 17(х-у)

Е |27Г7| 2т* СТ 7=0 1 "

(3)

С\

Лг

к=0 в=о 1 + Е (2пгЛ.-1в3)2т

3=1

(0) к Ск

1

р-1

х зяп С(0) ^

е-2пгЬ 1 ву

+ ^г С(0М ¿у

в=0 1 + £(2п^-1вз )2

3=1

где С к — решение системы уравнений

Л

Е

к=1

(—)тк е-2пгвкХк

+ С ,(0)

Е (2пгвз)2т

3=0

1

р-1

П

Е (2пгвз)2т

3=0

+ С,(0) | ¿ж = 0.

(4)

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

ХЛг

+ Кк С0)

к=0 в=0 1 + Е(2п^-1в3 )2т'

3=1

• \ 1/р' ¿ж . (5)

Асимптотическое выражение нормы оптимального функционала погрешности. В пространстве периодических функций выделим главный член нормы функционала погрешности (5), не зависящий от к = + Щ + • • • + Л.2 при стремлении последнего к нулю.

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

ж

шт* (л г)

р

Ьт = Ьт = • • • = Ьт = нт .

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

(6)

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

Ьк = Ьт /тк, N = Ыгп*/гпк, к = 1, 2,..., п.

Учитывая требования согласованности порядка сходимости с шагом решетки и гладкостью функции по координатным направлениям из условия (6), преобразуем интеграл в правой части (5):

J

Лг

Е

к=0

х ^ (-2тгОгк13к)тке-2тНк13кХк

/ У п

в=0 Е (2пiЬkвj)2т ¿=0

+ьтк Ск

(0)

¿х

= ьт*рьп /

(-2т£к )тк е-2пгвк Хк (0

+ С

Л к=0 в=0 )2т

¿=0

¿х. (7)

Равенство

(-2п^кГк е-2пгвх = Ьт* е-2пгвх имеет место при нулевом значении параметра к, и тогда интеграл (7) перепишется в виде

3--

Е

Ьт* е-2пгвкХк

Л в=0 Ь2т + £(2^)2т*

¿х

¿=1

+Е Е

(-2п^кГке-2пгвкХк , ^(0

+ С

к=1 в=0 Ь2т* + £ (2п^)2 ¿=1

¿х.

Теорема. Если в пространстве Ж^Д^) выполнены условия (2), (3),

1 1

€ 1 -р^1 Ет^1 > 1<Р<°°1 - + - = 1,

к=1 р р

Ьк — шаг решетки, норма функционала ?0(Ь-1х) определяется равенством (5), где параметры Ск0) — решения системы уравнений (4), то норма оптимального периодического функционала погрешности при Ь ^ 0 имеет асимптотическое выражение

Ь р

Е Е

(-2т£к )тк е-2пгвк Хк (0

+ С

чЛ к=1 в=0 Ь2т* + £ (2™^)2т

¿=1

1/р'

(1+ О(Ь)). (8)

р

р

р

р

Доказательство. Учитывая условия (6) и полагая в интеграле (5) -1Хк = ук, к = 1, 2,. .., п, получим

Н'о^МН^д.)

р' \ 1/р'

Е Е : " ; '

(0)

к=0 в=0 Е 1вз)

г 3=0

¿ж 1

у = ^ 1ж ¿ж = ж = ^у А^ ^ А

-1вкУк N \ 1/Р

(г ~ | иткГ<(

(/ е ( Е '-"т..+ ^) *

Л

/С:

=0 4в=0 Е

к=0

Преобразуем знаменатель:

к=0 к=0 к=0

п

Е (2^-1вк)2тк = Е |2пЛ-1&|2тк = Е |2пвк|2ткК

0

п

Е |2пвк|2тк^т*-2тк. (9) к=0

Далее,

Е |2пвк|2тк-2тк = ^2т* + |2пв1|2т1 -2т1 + |2пв2|2т2-2т2 к=0

+ ••• + |2пвп-1 |2т"-1 ^т*-2т"-1 + |2пвп|2т* > |2пвк|2т*. (10) С учетом (9) и (10) имеем

к=0 V в=0 — Е (2пгвк)2тк -2т к=0

(Г п / (—2ПУв )тк й2т*-тк р-2п®вкУк \ 1/Р

Е (Е( *»—^т^ +

Л к=Лв=0 — Е (2пгвк)2тк к=0

1 /р'

/е^Г-'^Е {-2т^тке^кУк + срАчу

\Л к=0 ^в=0 — Е (2пгвк)2тк^т*-2тк ' >

к=0

1/р'

< ^¿»»--.^( / (£ (-^Г-«"**» + с»') *)

к=0 ^ Чв=0 — Е (2пг&)2тк/ /(11)

Л

к=0

Перепишем цепочку (11):

|Ь(А-.)||' < /а» Е / Е Е + а

к=0Л к=0 в=0 1^Е(2пвк)2 к=1

р

к

(—2пгво)т° е-2пгв°х°

ЙГЛ" / /г2"1*1'' ■

А к=о в=о 1 + £ (2пвк)2тк

г к=1

+ Со

¿у

+ КЛ" / Л

(2т* — т1 )р'

Аг

к=0

^ (-2п^^1)т1 е-2пгв1Х1

-п-+ С1

в=о 1+£(2пвк)2т к=1

¿У

+ • • • + кл" / 53 Л

(2т*—т„)р'

Аг

к=0

^ (-2ТГ 1рп)тпе~2^пхп в=о 1 + £ (2пвк)2тк

¿У

к=1

А к=о в=0 1 + £ (2пвк)2тк к=1

¿у(1 + Й1Лр ).

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

||Го(Л

— 1х)1 р

< КЛ" / 53 Л(т*+")р' к=0

А

^ (-2п^вк)тке-2пгвкХк

-п--^

в=о 1 + £ (2пвк)2т к=1

¿у(1 + 0(ЛР )),

чА к=1 в=0 £(2пгв)2 ¿=1

1/р'

(1+ О(Л)).

Таким образом, норма (4) при Л ^ 0 имеет асимптотическое выражение (8). Теорема доказана.

ЛИТЕРАТУРА

1. Соболев С. Л. Введение в теорию кубатурных формул. М.: Наука, 1974.

2. Шойнжуров Ц. Б. Некоторые вопросы теории кубатурных формул в пространстве Wpm(E п) // Теория кубатурных формул и приложения функционального анализа к задачам математической физики. Новосибирск: Наука, 1980. С. 302—306.

3. Шойнжуров Ц. Б. Оценка нормы функционала погрешности кубатурных формул в различных функциональных пространствах. Улан-Удэ: Изд-во БНЦ СО РАН, 2005.

4. Рамазанов М. Д. Лекции по теории приближенного интегрирования. Уфа: Изд-во Башкир. ун-та, 1973.

5. Никольский С. М. Приближение функций многих переменных и теоремы вложения. М.:

Наука, 1977.

6. Бесов О. В., Ильин В. П., Никольский С. М. Интегральные представления функций и теоремы вложения. М.: Наука, 1975.

7. Владимиров В. С. Уравнения математической физики. М.: Наука, 1988.

Р

р

Р

Р

Р

8. Юмова Ц. Ж. Вычисление параметров функционалов погрешностей кубатурных формул. Saarbrücken: LAMBERT Acad. Publ., 2011.

Статья поступила 20 августа 2015 г. Юмова Цыренханда Жзмбэевна

Восточно-Сибирский гос. университет технологий и управления, ул. Ключевская, 40в, Улан-Удэ 670013 syum@mail .ru