Последовательности типа Грегори с неотрицательными коэффициентами

Вычислительные технологии

Том 10, Специальный выпуск, 2005

ПОСЛЕДОВАТЕЛЬНОСТИ ТИПА ГРЕГОРИ С НЕОТРИЦАТЕЛЬНЫМИ КОЭФФИЦИЕНТАМИ*

В. И. Половинкин, Л. В. ПоловинкинА Красноярский государственный технический университет, Россия e-mail: polovinkin@fivt.krasn.ru

Sequences of the Gregory type quadrature formulae with various degree of accuracy and nonnegative coefficients are proposed.

Данная работа по своей тематике связана с теорией приближенного интегрирования С.Л. Соболева [1, 2]. Он ввел формулы с регулярным пограничным слоем. Одним из авторов этой статьи были определены формулы с пограничным слоем [3-5]. Класс этих формул шире класса формул с регулярным пограничным слоем. Было установлено, что существуют кубатурные формулы с пограничным слоем, образующие асимптотически оптимальные последовательности в классах Lp^Q), pm > n, где О — ограниченная область n-мерного пространства, [6, 7], и асимптотически оптимальные в классах Lpm)(a, b) последовательности квадратурных формул с пограничным слоем при p G (1, то] [4, 5].

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

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

В.Л. Васкевич показал, что при достаточно больших m у квадратурных формул Грегори не все коэффициенты неотрицательны. Этот результат вытекает, например, из следствия 10.1 в [2, с. 342].

Важными характеристиками последовательностей типа Грегори являются введенные в [4, 5] им сопутствующие числа. Асимптотическая оптимальность любой последовательности типа Грегори в некотором классе l!^ (a, b) будет иметь место лишь при определенном значении ее сопутствующего числа.

В работе показано, что для любых чисел œ и параметра m, характеризующего точность рассматриваемых квадратурных формул, существует последовательность типа Грегори с неотрицательными коэффициентами и сопутствующим числом œ. Отсюда вытекает

* Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (грант № 03-01-00703).

© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2005.

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

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

С целью удобства изложения в настоящей статье используются термины, связанные не непосредственно с квадратурными формулами, а с их функционалами ошибок. Например, вместо последовательностей квадратурных формул типа Грегори будут рассматриваться последовательности их функционалов ошибок — последовательности типа Грегори. Далее, а, Ь, т, п, к будут такие числа, что а < Ь,т,п — натуральные, к = (Ь — а)п-1.

Разностные операторы Дг порядков г = 1, т вводятся так:

(Д(х), /(х)) = (А1/)(х) = /(х +1) — /(х),..., (Дт/)(х) =

= (Д(Д(т-1)/ ))(х) = (Д(т-1)/ )(х + 1) — (Д(т-1) / )(х).

Если т, п, 5 — числа, 5 — натуральное число, т = 0, то оператор Д*(хт-1 — п) определим так:

(Д*(хт-1 — п), /(х)) = т(Д*(х), /(тх + тп)).

Когда оператор ^ = т или ^ = т—1, Км будет означать наибольшее из числа сочетаний

С = м.

Определение 1. Последовательность функционалов {¡н} называется последовательностью типа Грегори, если 1Н имеют вид

Ь ( í п-г-1 í

а I г=0 г^+1 г=0 J

где т, п — натуральные числа; 2£ < п; аг, вг,г = 0, постоянные и выполняются условия

Определение 2. Если 1Н — функционалы вида (1), то числа к, каг,квг,г = 0,£, называются коэффициентами 1Н.

Определение 3. Последовательность типа Грегори {¡н} называется последовательностью с неотрицательными коэффициентами, если у всех функционалов Iн из этой последовательности каждый коэффициент неотрицателен.

Теорема 1. При любом т существуют последовательности типа Грегори с неот-рицательныи коэффициентами.

Лемма 1. Пусть заданы последовательность типа Грегори вида (1) и число г €{0,£}. Тогда существуют натуральные числа й > п1 > п, £ + й и последовательность типа Грегори {рн}

й п-1

г=0 г=й+1 г=0

Ь

п = го-

При этом г = 0, й, — коэффициенты такие, что:

а) при г € {0, £}\г а = а^;

б) а > 0, когда г € {£ +1, й};

в) аТ = аг + (Кт)-1.

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

Пусть натуральное число 5 таково, что коэффициенты функционалов 1Н, соответствующие точкам а + Кг + Квг, г = 1, т, равны К. Например, возьмем 5 = £ +1 — г и при этом будем считать п > г + те + Положим

ж=0

Опираясь на известную формулу [9, с. 1026]

т

к=0

можно проверить непосредственно, что {рвида [6] удовлетворяет утверждению леммы с й = г + те, п1 = г + т- + £ +1.

Следовательно, лемма верна.

Применяя в случае необходимости лемму 1 несколько раз, можно убедиться в справ-деливости такого утверждения.

Лемма 2. Пусть заданы последовательность типа Грегори вида (1) и число г € {0, £}. Тогда существует последовательность типа Грегори {рн} вида (3) такая, что аТ > 0, а также выполняются условия (4) и (5).

Прямым ее следствием является такой ее результат.

Лемма 3. Пусть задана последовательность типа Грегори {Iвида (1). Тогда существует последовательность типа Грегори {рн} вида (3) такая, что все а,г = 0,й, неотрицательны.

Аналогично лемме 3 может быть установлена справедливость следующего утверждения.

Лемма 4. Пусть задана последовательность типа Грегори вида (1). Тогда существуют натуральные числа й,п1,й < п1, и последовательность типа Грегори {рн} с рн вида

ь ( t и-4-1 а

(Л/)=/ /(х) й(х) — к\^аг/(а + Кг)+ ^ /(а + Кг) + ^ Д(6 — Кг) I,п = пТТГО

а I ¿=0 г=<+1 ¿=0 )

такая, что все в > 0 при г = 0, й.

Доказательства лемм 3 и 4 различаются главным образом тем, что вместо функционалов вида (6) надо рассматривать функционалы

") , f

ж=0

при необходимым образом выбранных s.

Пусть теперь {1h} — произвольная последовательность типа Грегори, удовлетворяющая (1) и (2), например последовательность функционалов ошибок квадратурных формул Грегори, точных на многочленах степени ниже m. Применяя к ней последовательно леммы 3 и 4, получим последовательность типа Грегори, удовлетворяющую утверждению теоремы. Следовательно, теорема верна.

Определение 4. Сопутствующим числом последовательности типа Грегори вида (1) {1h} называется величина

ж =

k^x

В работах [4, 5] показано, что у любой последовательности типа Грегори существует сопутствующее число.

Следующий результат обобщает теорему (1) и верен при любом m.

Теорема 2. Каково бы ни было число ж, существует последовательность типа Грегори с неотрицательными коэффициентами, у которой ж будет сопутствующим числом.

Доказательство. Пусть задано число ж. Выберем некоторую последовательность типа Грегори {1h} вида (1) c неотрицательными коэффициентами. Согласно теореме (1) такая последовательность существует.

Обозначим через ж(1) сопутствующее число {1h}. Рассмотрим вначале случай, когда

Пусть t — число, соответствующее 1h в формуле (1). Будем считать, что n > 2(t + m). Это условие не повлияет на справедливость теоремы.

Определим функционалы ih, ph из равенств

ж=0

h ; h , / |\-1,

ph = 1h +

Непосредственно проверяется равенство нулю функционалов ih на многочленах степени ниже m. Отсюда, а также из формул (10) и (11) следует,что {ph} будет последовательностью типа Грегори.

Из формулы (7), где m заменено на m - 1, и из неравенства (9) вытекает, что {ph} — последовательность функционалов с неотрицательными коэффициентами.

Формулы (11) и (8) показывают, что сопутствующее число ж(р) последовательности {ph} таково:

ж

h^0

Из равенства (10) имеем

n— t—m— 1 / / \ s

(¿h(x), xm) = (Am—11 (^ - j - lj - Am—1 (^ - j ) ,

n—t—m— 1

I Ix — a \

m

^A^^a - ^ , xm^ = (n - 2t - m - 1)hm+1m!,

откуда следует, что

lim{h—m(£h(x), xm)} = lim{((b - a)h—1 - 2t - m - 1)hm!} = (b - a)m!. (13) h^ü h^ü

Сравнивая формулы (12) и (13), получаем, что ж(р) = ж, а следовательно, при ж, удовлетворяющих неравенству (9), теорема верна.  Если ж не удовлетворяет неравенству (9), то надо выбрать конечный набор чисел ж1 ... ж5 так, чтобы

—1

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

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

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

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

[2] Соболев С.Л., Васкевич В.Л. Кубатурные формулы. Новосибирск: Изд-во ИМ СО РАН, 1996.

[3] Половинкин В.И. Последовательности функционалов с пограничным слоем // Сиб. мат. журн. 1974. Т. 15, № 2. С. 413-429.

[4] Половинкин В.И. Последовательности квадратурных формул с пограничным слоем // Теория кубатурных формул и приложения функционального анализа к задачам математической физики: Тр. семинара им. С.Л. Соболева: Сб. науч. тр. / АН СССР. Сиб. отд-ние. Ин-т математики. 1977. № 1. С. 149-158.

[5] Половинкин В.И. Последовательности квадратурных формул с пограничным слоем и последовательности типа Грегори // Квадратурные и кубатурные формулы. Решение функциональных уравнений. (Методы вычислений): Сб. науч. тр. / Ленингр. ун-т. 1981. Вып. 12. С. 7-25.

[6] Половинкин В.И. Асимптотическая оптимальность последовательностей формул с регулярным пограничным слоем при нечетных т // Сиб. мат. журн. 1975. Т. 16, № 2. С. 328-335.

[7] Половинкин В.И. Асимптотически наилучшие последовательности кубатурных формул // Сиб. мат. журн. 1975. Т. 16, № 6. С. 1255-1262.

[8] Половинкин В.И. Решетчатые кубатурные формулы с положительными коэффициентами // Теория кубатурных формул и приложения функционального анализа к задачам математической физики: Тр. семинара им. С.Л. Соболева: Сб. науч. тр. / АН СССР. Сиб. отд-ние. Ин-т математики. 1978. № 1. С. 79-87.

[9] Математическая энциклопедия. Т. 2. М.: Сов. энциклопедия. 1979.

Поступила в редакцию 2 ноября 2005 г.