Институт вычислительной математики и математической геофизики СОРАН



Всероссийская конференция по вычислительной математике КВМ-2011


Тезисы докладов


Аппроксимация функций и кубатурные формулы

Ненасыщаемый алгоритм решетчатых кубатурных формул для наилучших решеток

Рахматуллин Д.Я.

Учреждение Российской академии наук Институт математики с вычислительным центром Уфимского научного центра РАН (Уфа)

В конце 60-х годов прошлого века С.Л. Соболев ([1], [2]) предложил алгоритм формул высокой точности, объединяющий подходы функционального анализа и алгебраический. Именно, гладкость интегрантов задавалась принадлежностью их конкретному банаховому функциональному пространству B(Omega), качество формулы определялось нормой функционала погрешности с оптимизацией по коэффициентам. А сам алгоритм строился как сумма локальных формул для интегрирования по элементарным ячейкам решетки с помощью алгебраических формул с теми же узлами, точно интегрирующих многочлены до некоторой степени, связанной с гладкостью интегрантов. С.Л. Соболев установил, что при N, устремленном к бесконечности, последовательность функционалов погрешностей таких формул обладает свойством асимптотической оптимальности на пространствах интегрантов B(Omega)= L_2^(m)(Omega).

Свойство асимптотической оптимальности было установлено и для многих других банаховых пространств, обычно употребляемых в вычислительной математике. Важнейшей особенностью соболевского алгоритма является свойство пограничного слоя: узлам, удаленным от границы области на расстояние больше O(N^{-1/n}) соответствуют одинаковые коэффициенты. Это на порядок уменьшает объем вычислительной работы и позволяет установить для таких формул фактическую эквивалентность асимптотической и порядковой оптимальности на соболевских пространствах W_p^m.

Заложенная в алгоритме алгебраическая точность локальных формул на всех многочленах до некоторой степени M проявилась в условной "асимптотической ненасыщаемости" алгоритма. Именно, построенная последовательность кубатурных формул остается асимптотически (при N, стремящемся к бесконечности) оптимальной на каждом пространстве $W_{p}^{m}$ с m из (n/p, M), то есть в ослабленной форме удовлетворяет введенному К. И. Бабенко определению ненасыщаемости. (У К. И. Бабенко ([3], [4]) универсальна порядковая оптимальность без ограничения сверху на гладкость интегрантов). Ненасыщаемость - важное свойство вычислительных алгоритмов, позволяющее созданным по ним программам автоматически настраиваться на проявляющиеся характеристики гладкостей параметров задач, обеспечивая наилучшие скорости сходимостей аппроксимаций.

М.Д. Рамазановым был предложен ([5]) новый алгоритм решетчатых кубатурных формул, обладающих свойством ограниченного пограничного слоя и являющихся асимптотически ненасыщаемыми на всех пространствах $W_{2}^{m} (mathbb{R}^{n})$ с m>n/2. Этот алгоритм ограничен кубической решеткой узлов.

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

Литература.

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

2. Sobolev, S. L. The Theory of Cubature Formulas (Mathematics and Its Applications) / S. L. Sobolev, V. L. Vaskevich. - 1st edition. - Springer, 1997. - URL: http://www.amazon.com/Theory-Cubature-Formulas-Mathematics-Applications/dp/0792346319/ref=sr_1_1?ie=UTF8&s=books&qid=1274770504&sr=8-1. - ISBN 0792346319.

3. Бабенко, К. И. Основы численного анализа / К. И. Бабенко. - 2-е изд. - М. ; Ижевск: РХД, 2002. - 848 с. - ISBN 5-93972-162-1.

4. Belykh, V. N. On the best approximation properties of {$C^infty$}-smooth functions on an interval of the real axis (to the phenomenon of unsaturated numerical methods) / V. N. Belykh// Siberian Mathematical Journal. - 2005. - Vol. 46, No. 3. - P. 373-385. - URL: http://dx.doi.org/10.1007/s11202-005-0040-z.

5. Рамазанов, М. Д. Новый алгоритм асимптотически оптимальных решетчатых кубатурных формул / М. Д. Рамазанов// Уфимский математический журнал. - 2010. - С. 61-80.

Примечание. Тезисы докладов публикуются в авторской редакции



Ваши комментарии
Обратная связь
[ICT SBRAS]
[Головная страница]
[Конференции]

© 1996-2000, Институт вычислительных технологий СО РАН, Новосибирск
© 1996-2000, Сибирское отделение Российской академии наук, Новосибирск
    Дата последней модификации: 06-Jul-2012 (11:49:22)