WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

На правах рукописи

Широкий Александр Александрович

АППРОКСИМАЦИОННЫЕ СВОЙСТВА ТРИАНГУЛЯЦИЙ ПОВЕРХНОСТЕЙ

01.01.01 вещественный, комплексный и функциональный анализ А В Т О Р Е Ф Е Р А Т диссертации на соискание ученой степени кандидата физико-математических наук

Казань 2012

Работа выполнена на кафедре математического анализа и теории функций ФГБОУ ВПО “Волгоградский государственный университет”

Научный консультант: доктор физико-математических наук, доцент Клячин Владимир Александрович

Официальные оппоненты: доктор физико-математических наук, профессор Шабалин Павел Леонидович кандидат физико-математических наук, доцент Кесельман Владимир Михайлович

Ведущая организация: Институт математики им. С.Л. Соболева Сибирского отделения РАН

Защита состоится 17 мая 2012 года в 14 часов 30 минут на заседании диссертационного совета Д 212.081.10 при Казанском (Приволжском) федеральном университете по адресу: 420008, г. Казань, ул. Профессора Нужина, 1 / 37, ауд. 337.

С диссертацией можно ознакомиться в Научной библиотеке имени Н.И. Лобачевского Казанского (Приволжского) федерального университета по адресу:

420008, г. Казань, ул. Кремлёвская, 18.

Автореферат разослан 05 апреля 2012 года и размещён на официальном сайте Казанского (Приволжского) федерального университета: www.ksu.ru.

Учёный секретарь диссертационного совета Д 212.081.к. ф.-м. н., доцент Липачёв Е.К.

Работа выполнена при финансовой поддержке грантов РФФИ (проект 1101-97021-р_поволжье_а) и факультета математики и информационных технологий ВолГУ.

Общая характеристика работы

По своей тематике данная диссертационная работа выполнена на стыке нескольких разделов анализа. Главным вопросом исследования является вопрос об оценках степени аппроксимации производных гладких функций, заданных на поверхности, производными кусочно-линейных функций. Данная тематика тесно связана как с теорией расчётных сеток, так и с вопросами анализа, берущими своё начало от классического примера Карла Шварца (1890).

Основным объектом исследования являются триангуляции поверхностей.

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

В этом случае встаёт вопрос о точности построенной оценки, а также об устойчивости используемой разностной схемы. Например, в работах С.К. Годунова и Г.П. Прокопова1, 2 нерегулярные сетки используются для аппроксимации первых производных в составе эллиптических дифференциальных уравнений Лапласа, а также формулируются условия, которым должна удовлетворять сетка для обеспечения устойчивости результата численного моделирования. В более поздней работе С.К. Годунова, В.Т. Жукова, О.Б. Феодоритовой3 нерегулярные сетки применяются для расчёта инвариантных подпространств для симметрических гиперболических систем.

Однако первым заметным результатом, характеризующим зависимость качества приближения от способа построения разбиения стоит признать классический пример Карла Шварца4 (1890), известный также как “сапог Шварца”.

Он представляет собой семейство приближений кругового цилиндра с помощью полиэдральных поверхностей. Предельная площадь этих приближений может быть сделана произвольно большой. Эта конструкция позволила увиГодунов С.К. О расчетах конформных отображений и построений разностных сеток / С.К. Годунов, Г.П. Прокопов // Журнал вычислительной математики и математической физики. 1967. № 7 (5). С. 1031 – 1059.

Годунов С.К. О решении дифференциальных уравнений с использованием криволинейных разностных сеток / С.К. Годунов, Г.П. Прокопов // Журнал вычислительной математики и математической физики. 1968. № 8 (1). С. 28 – 46.

Годунов С.К. Метод расчета инвариантных подпространств для симметрических гиперболических уравнений / С.К. Годунов, В.Т. Жуков., О.Б. Феодоритова // Журнал вычислительной математики и математической физики. 1968. № 46 (6).

С. 1019 – 1031.

Контрпримеры в анализе / Б. Гелбаум, Дж. Олмстед. Волгоград : Платон, 1997. С. 191.

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

На самом деле геометрические параметры нерегулярной сетки оказывают непосредственное влияние на качество аппроксимации, вплоть до наличия/отсутствия сходимости. Геометрические свойства элементов нерегулярных сеток рассмаривались, например, в работах таких математиков, как S. Korotov1, V.A. Garanzha2, V.T. Rajan3 в работе последнего сформулирован ряд результатов, касающихся триангуляции Делоне в n-мерном евклидовом пространстве. В работе авторов H. Pottmann, R. Krasauskas, B. Hamann, K. Joy, W. Seibold4 исследуются кусочно-аффинные аппроксимации квадратичных функций в Rn и в том числе рассматривается вопрос об оптимальном приближении (с максимальными линейными участками, но при этом с приемлемой точностью) в R2 и отчасти в R3. J.R. Shewchuk5 поднимает вопрос о хорошем конечном линейном элементе в R2 и в R3 с точки зрения аппроксимации градиента функции. В работе Е.А. Пабат и В.А. Клячина6 исследуются свойства кусочно-линейных интерполяций поверхностей уровня функций, заданных на нерегулярных сетках. Были найдены достаточные условия на триангуляции -сетей, при которых имеет место C1-аппроксимация поверхностей уровня C2-гладких функций.

При решении ряда практических задач возникает необходимость в построении триангуляций различных поверхностей. В частности, В.М. Миклюковым7 исследовалась проблема триангуляции локально липшицевой поверхности. Учитывая “популярность” триангуляции Делоне, естественным, казалось бы, образом возникает желание построить соответствующее обобщение на случай поверхностей. Однако, несмотря на то, что алгоритм построения триангуляций Делоне в евклидовых пространствах произвольной конечной размерноKorotov S. Some geometric results for tetrahedral finite elements. // Proceedings of the International Conference NUMGRID-2010.

M. : Фолиум, 2010. С. 41 – 46.

Garanzha V.A. Discrete Extrinsic Curvatures and Approximation of Surfaces by Polar Polyhedra // Computational Mathematics and Mathematical Physic. 2010. Vol. 50. № 1. С. 65 – 92.

Rajan V.T. Optimality of the Delaunay triangulation in Rd. Discrete Computational Geometry. 1994. Vol. 12. С. 189 – 202.

Pottmann H., Krasauskas R., Hamann B., Joy K., Seibold W. On Piecewise Linear Approximation of Quadratic Functions. // Journal for Geometry and Graphics. 2000. Vol. 4. № 1. С. 31 – 53.

Shewchuk J.R. What is a Good Linear Element? Interpolation, Conditioning, and Quality Measures / J.R. Shewchuk // In Proceedings, 11th International Meshing Roundtable (September 2002). Ithaca, New York, USA : Sandia National Laboratories, 2002. С. 115 – 126.

Клячин В.А. C1-аппрроксимация поверхностей уровня функций, заданных на нерегулярных сетках / В.А. Клячин, Е.А. Пабат // Сиб. журн. индустр. матем. 2010. № 13 (2). С. 69 – 78.

Миклюков В.М. Введение в негладкий анализ. Изд. 3, перераб. / В.М. Миклюков. Волгоград : изд-во ВолГУ. 2011. С. 77.

сти был предложен ещё в 1934 году самим Б.Н. Делоне1, известные алгоритмы построения триангуляции поверхности предполагают построение обычной плоской триангуляции с последующим проецированием её на поверхность (см., например, работы А.В. Скворцова и Н.С. Мирзы2). Ясно, что проекция плоской триангуляции Делоне на поверхность в общем случае не будет проявлять никаких замечательных свойств, характерных для триангуляции Делоне. Таким образом, возникает задача построения триангуляции Делоне на поверхности с сохранением свойств, присущих её классическому варианту.

Объектом исследования данной диссертационной работы являются триангуляции поверхностей.

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

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

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

Научная новизна исследования. Все результаты, полученые в работе, являются новыми.

Результаты, выносимые на защиту

. На защиту выносятся следующие основные результаты:

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

2. Оценки погрешности аппроксимации градиентов в евклидовой метрике для триангуляции Делоне плоской области и остроугольной триангуляции в пространстве.

3. Оценки погрешности аппроксимации градиентов для остроугольной триангуляции в метрике поверхности.

Delaunay B.N. Sur la sphre vide. A la mmoire de Georges Voronoї / B.N. Delaunay // Известия АН СССР. 1934. № 6. С. 793 – 800 / Перевод с фр. А. Ю. Игумнов в сб. Записки семинара “Сверхмедленные процессы”. Вып. 1. Волгоград : Изд-во ВолГУ, 2008. C. 147 – 153.

Скворцов А.В. Триангуляция Делоне и её применение / А.В. Скворцов. Томск : Изд-во Томского ун-та, 2002. С. 96.

4. Теоремы о сходимости градиентов кусочно-аффинной аппроксимации, построенной над триангуляциями -сетей в Rn, к градиентам C2-гладкой функции.

5. Обобщение условия пустого шара на случай n-мерных гиперповерхностей в (n+1)-мерном пространстве. Оценка погрешности вычисления градиентов для триангуляций Делоне двумерных поверхностей.

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

Апробация работы. Основные результаты работы докладывались на международных и российских конференциях: девятой международной школеконференции “Теория функций, её приложения и смежные вопросы” (Казань, 1 – 7 июля 2009 г.), 15-й Саратовской зимней школе “Современные проблемы теории функций и их приложения” (Саратов, 27 января – 3 февраля 20г.), семинаре механико-математического факультета Томского государственного университета под руководством проф. А. В. Скворцова, семинаре-совещании “Сети в анизотропных пространствах” (Волгоград, 21 – 23 апреля 2011 г.), десятой международной школе-конференции “Теория функций, её приложения и смежные вопросы” (Казань, 30 июня – 6 июля 2011 г.), 16-й Саратовской зимней школе “Современные проблемы теории функций и их приложения” (Саратов, 27 января – 3 февраля 2012 г.), а также на конференциях и семинарах Волгоградского государственного университета.

Публикации. Основные результаты диссертации опубликованы в работах [1] [10] общим объёмом 2,9 пл. При этом статьи [1], [2] опубликованы в журналах, рекомендованных Высшей аттестационной комиссией (ВАК) Министерства образования и науки Российской Федерации для публикации результатов научных исследований.

Структура и объём диссертации. Диссертация изложена на 92 страницах и состоит из введения, трёх глав и списка литературы. Библиография диссертации содержит 44 наименования, включая работы автора.

Краткое содержание работы Все утверждения сохраняют принятую в основном тексте нумерацию.

Во введении даётся обзор литературы по вопросам, связанным с темой диссертации, а также кратко излагаются основные результаты диссертации.

Глава 1. “Априорные оценки погрешности аппроксимации градиента” §1.1 носит вводный характер. В нём сформулированы основные понятия и определения, необходимые для изложения результатов.

В §1.2 формулируется основной результат первой главы и следствия из него.

Пусть D Rn некоторая область в Rn. Для функций класса C1(D) через f(x) обозначим вектор градиента в точке x D, а производную по направлению произвольного вектора через скалярное произведение векторов f(x) и :

f = f(x), .

Если в Rn имеется набор точек p0, p1,..., pk, где 1 k n, то обозначим через S(p0,..., pk) k-мерный симплекс с вершинами в этих точках. Будем предполагать, что векторы p1 - p0, p2 - p0,..., pk - p0 линейно независимы.

По отношению к плоскости симплекса S произвольный вектор v Rn может быть разложен в сумму v = vT + vN, где vT ортогональная проекция вектора v на плоскость симплекса S, а vN ортогональное дополнение проекции v.

Зафиксируем произвольный симплекс S D. Для функции f C1(D) построим аффинную функцию Lf(x) вида Lf(x) = a, x + b, a Rn, b R такую, что f(pi) = Lf(pi), при i = 0,..., k и aN = 0.

Пусть вектор-функция F : D Rm непрерывна. Тогда функцию (F, t) = sup{|F (x1) - F (x2)| : x1, x2 D, и |x1 - x2| < t} назовём модулем непрерывности функции F.

Величину f(p0) = f(p0) - Lf(p0) назовём погрешностью вычисления f(p0) посредством функции Lf(x).

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

Обозначим I() = (f, t)dt + (f, ).

Теорема 1.2.1. Пусть S D некоторый треугольник с заданной в нём аффинной функцией Lf(x), построенной для некоторой непрерывно дифференцируемой в области D функции f(x). Тогда, если S и d+ обозначают S значение угла при вершине p0 и длину максимальной стороны треугольника S соответственно, то справедлива оценка:

3I(d+) - (f, d+) T S S f(p0) - Lf(p0) .

sin S Теорема 1.2.2. Пусть k = 2 и вершины треугольника S лежат на двумерной гладкой поверхности F D Rn. Предположим, что угол между плоскостью этого треугольника и касательной плоскостью к поверхности в точке p0 равен , 0 . Если функция f C1(D), то выполнено:

3I(d+) - (f, d+) S S |f(p0) - Lf(p0)| + sin sup f(x).

sin S xD Здесь обозначает градиент дифференцируемой функции в метрике поверхности F.

Теорема 1.2.3. Пусть S D некоторый тетраэдр с заданной в нём аффинной функцией Lf(x), построенной для некоторой непрерывно дифференцируемой в области D функции f(x). Обозначим через S величину угла между вектором p2 - p0 и двумерной гранью, лежащей напротив вершины p2, через S величину угла при вершине p1 треугольника p0p1p3 и через d+ S длину максимального ребра тетраэдра S. Тогда имеет место оценка:

8I(d+) T S f(p0) - Lf(p0) sin S sin S Теорема 1.2.4. Пусть вершины тетраэдра S лежат на трехмерной гладкой поверхности F D Rn. Предположим, что угол между плоскостью тетраэдра и касательной плоскостью к поверхности в точке p0 равен , 0 . Если функция f C1(D), то выполнено:

8I(d+) S |f(p0) - Lf(p0)| + sin sup f(x).

sin S sin S xD Глава 2. “Аппроксимационные свойства триангуляций” В данной главе приводятся оценки аппроксимации градиентов для различных типов триангуляций, а также формулируются теоремы о сходимости градиентов кусочно-аффинной аппроксимации к исследуемой с её помощью функции.

Пусть в области D Rn задана последовательность конечных наборов точек. Для каждого такого набора рассмотрим его триангуляцию.

Определение 2.0.1. Триангуляцией конечного набора точек назовём набор не имеющих общих внутренних точек симплексов с вершинами из этого набора, образующих покрытие участка пространства, ограниченного выпуклой оболочкой рассматриваемого набора точек.

Определение 2.0.2. Триангуляцию будем называть остроугольной, если для каждого симплекса углы между двумя любыми смежными его гранями острые.

Определение 2.0.3. Триангуляция набора точек называется триангуляцией Делоне, если описанная сфера всякого симплекса триангуляции не содержит внутри себя каких-либо точек из этого набора.

Для всякого симплекса S длину максимальной его стороны обозначим d+.

S Для некоторой триангуляции T положим d+(T ) = max d+.

S ST Будем рассматривать такие наборы точек Pm и их триангуляций Tm, для которых выполнены следующие условия:

d+(Tm) 0 при m и (2.1) > 0 m0 N : m > m0 и x D a Pm : |a - x| < . (2.2) Условие (2.2) означает, что Pm является -сетью при всех достаточно больших m.

Рассмотрим некоторую функцию f(x) класса C1(D) и для каждого m построим кусочно-аффинную функцию fm(x) такую, что fm(a) = f(a) для любой точки a Pm.

Заметим, что при выполнении условий (2.1) и (2.2) последовательность fm(x) равномерно сходится к функции f(x) на каждом компактном подмножестве K D. Но эти условия не обеспечивают сходимости производных функций fm(x) к производным функции f(x) даже почти всюду или по норме 1, в пространствах Соболева Wlocp(D). Это подтверждается хорошо известным примером Шварца.

Далее будем предполагать, что функция f : D R является C2-гладкой в области D Rn, причём будем считать, что все вторые производные ограничены в D. Тогда найдётся постоянная M > 0 такая, что модуль непрерывности её градиента (f, t) M · t.

В §2.1 рассматриваются аппроксимационные свойства триангуляций Делоне плоских областей.

Оценка погрешности аппроксимации градиента для триангуляции Делоне плоской -сети может быть получена как следствие теоремы 1.2.1.

Теорема 2.1.1. Пусть задана триангуляция Делоне T конечного набора точек, являющегося -сетью в плоской области D R2. Пусть точка p является вершиной триангуляции T такой, что для некоторого треугольника из T, имеющего вершиной эту точку, описанный круг полностью лежит в области D. Тогда справедливо следующее неравенство:

f(p) - Lf(p) 14M.

Пусть Gm множество точек, составляющее вершины и рёбра симплексов триангуляции Tm.

Следствие 2.1.1. Пусть f(x) C2-гладкая в области D R2 и над последовательностью Tm триангуляций Делоне области D, удовлетворяющих условиям (2.1) и (2.2) построена её кусочно-аффинная аппроксимация fm(x).

Тогда справедливо равенство:

lim fm(x) = f(x), x D \ Gm. (2.3) m m В §2.2 рассматриваются аппроксимационные свойства остроугольных триангуляций в пространстве.

В трёхмерном евклидовом пространстве выполнение условия Делоне не гарантирует сходимости градиентов аппроксимации к градиентам исследуемой функции. Однако существует класс триангуляций (остроугольные триангуляции), для которого можно получить конечную оценку разности градиентов как следствие теоремы 1.2.3. Если мы будем рассматривать остроугольные триангуляции -сетей с узлами в некоторой области D R3, то при 0 оценка также устремится к нулю.

Теорема 2.2.1. Пусть задана остроугольная триангуляция T конечного набора точек, являющегося -сетью в области D R3. Пусть точка p является вершиной триангуляции T. Тогда 48M f(p) - Lf(p) .

2 - Следствие 2.2.2. Пусть f(x) C2-гладкая в области D R3 и над последовательностью Tm остроугольных триангуляций области D, удовлетворяющих условиям (2.1) и (2.2) построена её кусочно-аффинная аппроксимация fm(x). Тогда справедливо равенство:

lim fm(x) = f(x), x D \ Gm. (2.4) m m В §2.3 рассматриваются аппроксимационные свойства триангуляций в Rn.

Пусть D некоторая область в Rn. Зафиксируем некоторую маленькую величину > 0. Будем говорить, что триангуляция T является триангуляцией с углами, отделёнными от нуля, если для всякого невырожденного симплекса S T ни один из его двугранных углов не будет меньше . Рассмотрим последовательность Tm таких триангуляций области D. Положим также d-(Tm) = min d-, STm S где d- длина самого короткого ребра симплекса S.

S Справедлива следующая теорема.

Теорема 2.3.1. Предположим, что для области D Rn, последовательности наборов точек Pm и их триангуляций Tm с углами, отделёнными от нуля, выполнены условия (2.1) и (2.2). Тогда для любой функции f C2(D) и любой компактно вложенной подобласти U D имеет место неравенство n max max f(x) - fm(x) d+(Tm) 1 +.

sTm, SU xS sinn+1 Следствие 2.3.1. Пусть f(x) C2-гладкая в области D Rn и над последовательностью Tm триангуляций с углами, отделёнными от нуля, области D, удовлетворяющих условиям (2.1) и (2.2) построена её кусочно-аффинная аппроксимация fm(x). Тогда справедливо равенство:

lim fm(x) = f(x), x D \ Gm. (2.9) m m В §2.4 оценка разности градиентов уточняется для функций, являющихся решениями эллиптических уравнений вида Lu = aijDiju + biDiu + cu = f, где aij = aij(x), bi = bi, c = c(x) функции со значениями из R, f C(), а любое открытое подмножество в Rn.

Зафиксируем некоторый n-мерный симплекс S Rn и свяжем с ним ортонормированный базис {eS} как результат процесса ортогонализации Граммаi Шмидта векторов {pi - p0}, i = 1,..., n. Пусть k обозначает угол между вектором pk - p0 и плоскостью, натянутой на векторы e1,..., ek-1, k = 2,..., n; 1 =, а k обозначает угол между гранью Sk и плоскостью k.

Пусть um(x) кусочно-аффинная аппроксимация u такая, что в узлах триангуляции её значения совпадают со значениями u. Тогда справедлива следующая теорема.

Теорема 2.4.2. Предположим, что для области Rn, последовательности наборов точек Pm и их триангуляций Tm выполнены условия (2.1) и (2.2). Тогда для любой ограниченной функции u C2, (), являющейся решением эллиптического уравнения Lu = f в c правой частью f C(), и любой компактно вложенной подобласти U имеет место неравенство max max |u(x) - um(x)| STm, SU xS n d+(Tm) 1 + max, (2.11) STm, SU |sin i| |cos i| i=где |f(x) - f(y)| = C1 sup |u| + sup d2 |f(x)| + sup d2+, C1 = const.

x x, y |x - y| x x x, y, x =y В §2.5 рассматриваются аппроксимационные свойства остроугольных триангуляций трёхмерных поверхностей.

Рассмотрим некоторую n-мерную гиперповерхность F в пространстве Rn+1, заданную графиком C2-гладкой функции xn+1 = u(x1,..., xn), x D Rn.

Определим отображение проекции : F D, (x1,..., xn+1) = (x1,..., xn).

Пусть P некоторый конечный набор точек в Rn+1, лежащих на поверхности F и таких, что любой симплекс с вершинами из P является невырожденным. Обозначим через p = (p) проекцию точки p в область D. Будем рассматривать триангуляцию поверхности F n-мерными симплексами, построенными на множестве P. Триангуляцией поверхности назовём такой набор T n-мерных симплексов, что 1. Каждая точка p заданного набора является вершиной одного из симплексов S T.

2. Каждая вершина любого симплекса S T является одной из точек набора P.

3. Внутренность пересечения любых двух симплексов пуста.

4. Проекция системы симплексов Sj T, j = 1,..., m является обычной триангуляцией набора точек P в Rn. Здесь P = {p | p = (p); p P }.

Будем рассматривать последовательность таких наборов точек Pm, m = 1, 2,..., расположенных на поверхности F, и соответствующих им триангуляций Tm, удовлетворяющих условиям (2.1) и (2.2). Для каждого m построим кусочно-аффинную функцию fm(x) такую, что fm(a) = f(a) для любой точки a Pm.

Теорема 2.5.1. Пусть D область, расположенная на некоторой трёхмерной поверхности F. Рассмотрим остроугольную триангуляцию T конечного набора точек P, являющегося -сетью в области D. Пусть точка p является вершиной триангуляции T, а угол между плоскостью тетраэдра S и касательной плоскостью к поверхности F в точке p. Тогда 48M |f(p) - Lf(p)| + sin sup f(x).

xD 2 - Следствие 2.5.1. Для построенной выше последовательности функций fm(x), x D над остроугольными триангуляциями Tm трёхмерной поверхности F, удовлетворяющими условиям (2.1) и (2.2), выполнено равенство:

lim fm(x) = f(x), x D \ Gm. (2.12) m m Замечание. Для остроугольных триангуляций двумерных поверхностей также можно получить подобную оценку. Однако в двумерном случае возможно построение оценки для более общего класса триангуляций триангуляций Делоне. Мы вернёмся в этому вопросу в §3.3.

Результаты данной главы были опубликованы в работах [1], [3] [6].

Глава 3. “Обобщения условия Делоне” В данной главе формулируются некоторые обобщения условия пустого шара. С точки зрения задачи аппроксимации производных такие обобщения могут оказаться полезными в зависимости от конкретной вычислительной задачи.

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

Рассмотрим в Rn семейство выпуклых множеств F (x, r), x Rn, r R и множество точек P, расположенных в некоторой области D Rn.

Пусть S произвольный невырожденный симплекс в Rn. Определим описанное множество (если оно существует) F (S) из семейства F (x, r) как множество, чья граница содержит все вершины симплекса (а значит, F (S) содержит весь симплекс в силу выпуклости F (x, r)).

Рассмотрим произвольную триангуляцию множества точек P.

Определение 3.1.1. Будем говорить, что триангуляция равномерна относительно семейства F (x, r), если для любого симплекса S этой триангуляции внутренность множества F (S) не содержит вершин других симплексов.

Определение 3.1.2. Будем говорить, что триангуляция локально равномерна относительно семейства F (x, r), если для любого симплекса S этой триангуляции внутренность множества F (S) не содержит вершин других симплексов, имеющих общую (n - 1)-мерную грань.

Теорема 3.1.1.1 Для того, чтобы для семейства выпуклых множеств F (x, r) из локальной равномерности следовала глобальная равномерность, достаточно, чтобы указанное семейство обладало свойством: для любого симплекса S существует и единственно описанное множество F (S).

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

Теорема 3.1.2. Пусть (x) строго выпуклая вниз функция, x Rn. Тогда семейство выпуклых множеств F (x, r) = {y : (y) = (x), y - x + r} удовлетворяет условию теоремы 3.1.1.

Замечание. Обычная триангуляция Делоне получается в случае (x) = |x|2.

Эта теорема была доказана в следующей работе: Клячин, В.А. Об одном обобщении условия Делоне / В.А. Клячин // Вестник Томского государственного университета, Математика и механика, 2008, № 1 (2). С. 48 – 50.

В §3.2 рассматривается вариант обобщения понятия триангуляции Делоне на случай пространства с метрикой.

Рассмотрим произвольный конечный набор точек P в некоторой области Rn с метрикой и обозначим квадрат элемента длины ds2 = gijdxidxj.

i, j Здесь gij C(), и в каждой точке x образуют положительно определённую матрицу.

Определение 3.2.1. Для произвольного симплекса S = (p1,..., pk) описанным шаром B(S) в метрике назовём такой шар B(x0, r), где x0 центр шара, а r его радиус, что pi B(x0, r) и x0 лежит в плоскости симплекса S(p0,..., pk).

Далее будем рассматривать только такие метрики, которые определены элементом длины ds и шары в метрике являются выпуклыми множествами. Также предположим, что для каждого симплекса S его описанный шар существует и единственен.

Определение 3.2.2. Если внутрь шара в метрике, описанного вокруг любого построенного симплекса триангуляции над P, не попадает ни одна из точек этого множества, то будем говорить, что триангуляция удовлетворяет условию Делоне в метрике и называется триангуляцией Делоне в метрике.

Определение 3.2.3. Будем говорить, что для двух соседних симплексов с вершинами из P выполнено условие пустоты шара, если описанный шар одного симплекса не содержит внутри себя вершин другого симплекса.

Справедлива следующая теорема.

Теорема 3.2.1. Если для любой пары соседних симплексов триангуляции T выполнено условие пустоты шара, то T триангуляция Делоне в метрике.

Для триангуляций -сетей в метрике можно построить оценку разности градиентов. Будем использовать обозначения, принятые во второй главе. Сформулируем условия, аналогичные (2.1) и (2.2):

При m d+(Tm) = max max d(pk, pl) 0, где pk, pl S, k = l; (3.1) STm k, l > 0 m0 N : m > m0 и x D a Pm : d(a, x) < . (3.2) Далее будем рассматривать триангуляции в метрике, удовлетворяющие данным условиям. Справедлива следующая теорема.

Теорема 3.2.4. Предположим, что для области D Rn с некоторой метрикой для последовательностей наборов точек Pm и соответствующих им триангуляций Tm выполнены условия (3.1) и (3.2). Пусть U D некоторая подобласть такая, что для некоторого натурального kU S D при всех m > k0.

STm Тогда при m > k0 выполнено следующее неравенство:

max max |f(x) - fm(x)|g STm,SU xS j n 1 d+(Tm) + n Mg max, STm, SU |sin i| |cos i| j=1 i=где Mg максимальное собственное значение матрицы gij.

Для триангуляций Делоне в метрике также справедливо утверждение, касающееся сходимости градиентов.

Теорема 3.2.5. Для построенной выше последовательности функций fm(x), x D над двумерными триангуляциями Делоне в Tm в метрике, удовлетворяющими условиям (3.1) и (3.2), выполнено равенство:

lim fm(x) = f(x), x D \ Gm. (3.7) m m В §3.3 строится обобщение условия Делоне на случай n-мерной гиперповерхности в (n + 1)-мерном пространстве. Для триангуляций Делоне двумерных поверхностей сформулированы и доказаны теоремы об оценке разности градиентов, а также о сходимости градиентов.

Пусть P некоторый набор точек в Rn+1, лежащих на гиперповерхности F и таких, что любой симплекс в вершинах из P является невырожденным.

Пусть S некоторый n-мерный симплекс в Rn+1.

Определение 3.3.1. Описанным шаром B(S) для S назовём (n + 1)-мерный шар, содержащий на своей границе все вершины симплекса и имеющий наименьшее возможное значение радиуса.

Определение 3.3.2. Будем говорить, что триангуляция поверхности F удовлетворяет условию Делоне, если для любого симплекса триангуляции описанный шар в смысле определения 3.3.1 не содержит внутри себя ни одной вершины триангуляции.

Определение 3.3.3. Будем говорить, что для двух n-мерных симплексов, пересекающихся по общей (n - 1)-мерной грани выполнено условие пустоты шара, если (n + 1)-мерный описанный шар одного симплекса не содержит внутри себя вершин другого симплекса.

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

Условие ). Рассмотрим произвольный n-мерный симплекс S триангуляции и построим пересечение US всех полупространств, содержащих S, и которые определяются гиперплоскостями (n + 1)-мерных сфер, являющихся пересечениями описанной сферы симплекса S и его соседнего. Таких полупространств будет не более n + 1 для каждого симплекса S. Мы предполагаем, что никакая вершина триангуляции не принадлежит внутренности US для всякого S.

Заметим, что в плоском случае, когда вершины триангуляции лежат в одной гиперплоскости, условие очевидно выполняется.

Теорема 3.3.1. Пусть n-мерная гиперповерхность F задана над выпуклой областью D Rn. Триангуляция гиперповерхности F, для которой любая пара симплексов, пересекающихся по общей (n - 1)-мерной грани, удовлетворяет условию пустоты шара, является триангуляцией Делоне.

Ранее в §2.5 была получена оценка погрешности аппроксимации градиента для остроугольных триангуляций трёхмерных поверхностей. Теперь можно из теоремы 1.2.3 получить соответствующую оценку для триангуляций Делоне двумерных поверхностей. Рассмотрим последовательность наборов точек Pm, m = 1, 2,..., расположенных на поверхности F, и соответствующих им триангуляций Tm, удовлетворяющих условиям (2.1) и (2.2). Для каждого m построим кусочно-аффинную функцию fm(x) такую, что fm(a) = f(a) для любой точки a Pm.

Теорема 3.3.4. Пусть задана триангуляция Делоне T конечного набора точек, являющегося -сетью области D, расположенной на двумерной поверхности F R3. Пусть точка p является вершиной триангуляции T такой, что для некоторого треугольника S из T, имеющего вершиной эту точку, пересечение его описанной окружности с поверхностью F полностью лежит в области D. Пусть угол между плоскостью треугольника S и касательной плоскостью к гиперповерхности F в точке p. Тогда справедливо следующее неравенство:

|f(p) - Lf(p)| 14M + sin sup f(x).

xD Следствием этой теоремы будет утверждение о сходимости градиентов кусочно-аффинной функции fm(x) к градиентам исследуемой с помощью неё функции f(x).

Следствие 3.3.1. Для построенной выше последовательности функций fm(x), x D над триангуляциями Делоне двумерной поверхности F, удовлетворяющими условиям (2.1) и (2.2), выполнено равенство:

lim fm(x) = f(x), x D \ Gm. (3.8) m m В §3.4 рассматриваются триангуляции Делоне с коэффициентом сжатия сфер.

Пусть P произвольный конечный набор точек некоторой области D Rn.

Зафиксируем число 0 < 1. Пусть S произвольный симплекс в Rn.

Определение 3.4.1. Шар BS с центром, совпадающим с центром описанного около S шара радиуса RS и радиусом RS будем называть -шаром сим плекса S. В случае R2 будем называть BS -кругом.

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

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

Для триангуляции с коэффициентом сжатия описанной окружности можно построить оценку аппроксимации градиента функции f(x), x D, C2-гладкой в области D. Здесь и далее мы будем использовать обозначения, принятые во второй главе.

Теорема 3.4.1. Пусть задана -триангуляция Делоне T конечного набора точек, являющегося -сетью в плоской области D R2. Пусть точка p является вершиной триангуляции T такой, что для некоторого треугольника из T, имеющего вершиной эту точку, его -круг полностью лежит в области D. Тогда f(p) - Lf(p) M.

Следствие 3.4.1. Пусть f(x) C2-гладкая в области D R2 и над последовательностью Tm -триангуляций Делоне области D, удовлетворяющих условиям (2.1) и (2.2) построена её кусочно-аффинная аппроксимация fm(x).

Тогда справедливо равенство:

lim fm(x) = f(x), x D \ Gm.

m m Результаты данной главы были опубликованы в работах [2], [4] [9].

Пользуясь случаем, автор выражает глубокую благодарность за полезные обсуждения и замечания своему научному руководителю, д. ф.-м. н. В.А. Клячину, а также к. ф.-м. н., доценту В.В. Попову, к. ф.-м. н. Н.М. Полубояровой и всем участникам семинара “Сверхмедленные процессы” под руководством д. ф.-м. н., профессора В.М. Миклюкова. Хотелось бы отдельно выразить свою искреннюю признательность декану факультета математики и информационных технологий ВолГУ, д. ф.-м. н., профессору А.Г. Лосеву за внимательное отношение и энергичную поддержку.

Публикации по теме диссертации Статьи в ведущих рецензируемых научных журналах, включённых в список ВАК:

[1] Широкий, А.А. Триангуляция Делоне многомерных поверхностей [Текст] / В.А. Клячин, А.А. Широкий // Вестник Самарского государственного университета. Естественно-научная серия. 2010. № 2010 / 4 (78).

С. 51 – 55. – 0,3 п. л.

[2] Широкий, А.А. Триангуляция Делоне многомерных поверхностей и её аппроксимационные свойства [Текст] / В.А. Клячин, А.А. Широкий // Известия вузов. Математика. 2010. № 1. С. 31 – 39. – 0,5 п. л.

Публикации в других изданиях:

[3] Широкий, А.А. Аппроксимационные свойства остроугольных триангуляций [Текст] / В.А. Клячин, А.А. Широкий // Материалы научной сессии, г. Волгоград, 26 – 30 апр. 2010 г. Вып. 6. Математика и информационные технологии. Волгоград : Изд-во ВолГУ, 2010. С. 40 – 44. – 0,3 п. л.

[4] Широкий, А.А. Аппроксимационные свойства триангуляции Делоне [Текст] / В.А. Клячин, А.А. Широкий // Записки семинара “Сверхмедленные процессы”. Вып. 5. Волгоград : Изд-во ВолГУ, 2010. С. 8 – 14.

– 0,5 п. л.

[5] Широкий, А.А. Триангуляция Делоне многомерных поверхностей [Текст] / В.А. Клячин, А.А. Широкий // Материалы 15-й Саратовской зимней школы “Современные проблемы теории функций и их приложения”. Саратов : Изд-во СГУ, 2010. С. 89. – 0,1 п. л.

[6] Широкий, А.А. Триангуляция Делоне многомерных поверхностей и её аппроксимационные свойства [Текст] / В.А. Клячин, А.А. Широкий // Вестник Волгоградского государственного университета. Серия 1: Математика, физика. 2009. № 12. С. 39 – 44. – 0,4 п. л.

[7] Широкий, А.А. Триангуляция Делоне многомерных поверхностей и её аппроксимационные свойства [Текст] / В.А. Клячин, А.А. Широкий // Материалы Девятой международной научной школы-конференции “Теория функций, её приложения и смежные вопросы”. Казань : Казан. мат.

об-во, 2009. С. 155 – 157. – 0,2 п. л.

[8] Широкий, А.А. Триангуляция Делоне с коэффициентом сжатия сфер [Текст] / В.А. Клячин, А.А. Широкий // XV региональная конференция молодых исследователей Волгоградской области, 9 – 12 нояб. 2010. Вып.

4. Физика и математика. Волгоград : Изд-во ВолГУ, 2011. С. 71 – 74.

– 0,3 п. л.

[9] Широкий, А.А. Условие почти пустого шара. [Текст] / А.А. Широкий // Материалы Десятой международной научной школы-конференции “Теория функций, её приложения и смежные вопросы”. Казань : Казан. мат.

об-во, 2011. С. 387 – 389. – 0,2 п. л.

[10] Широкий, А.А. Применение триангуляции Делоне двумерной поверхности к примеру Шварца [Текст] / А.А. Широкий // Материалы 16-й Саратовской зимней школы “Современные проблемы теории функций и их приложения”. Саратов : Изд-во СГУ, 2012. С. 200. – 0,1 п. л.




© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.