Электронная библиотека Веда
Цели библиотеки
Скачать бесплатно
Доставка литературы
Доставка диссертаций
Размещение литературы
Контактные данные
Я ищу:
Библиотечный каталог российских и украинских диссертаций

Вы находитесь:
Диссертационные работы России
Физико-математические науки
Вычислительная математика

Диссертационная работа:

Роженко Александр Иосифович. Теория и алгоритмы вариационной сплайн-аппроксимации : Дис. ... д-ра физ.-мат. наук : 01.01.07 : Новосибирск, 2003 231 c. РГБ ОД, 71:05-1/136

смотреть содержание
смотреть введение
Содержание к работе:

Введение 4

Глава 1. Элементы функционального и выпуклого анализа 29

  1. Определения и обозначения 30

  2. Условия нормальной разрешимости линейного ограниченного оператора 36

  3. Эквивалентные нормы в банаховых пространствах 39

  4. Линейная независимость операторов 42

  5. Наилучшее приближение в выпуклом множестве 44

  6. Сходимость наилучших приближений 45

1.7. Пространство С\ос(0) на локально-компактном множестве . . 46
Глава 2. Теория абстрактных сплайнов 49

  1. Однозначная разрешимость смешанной задачи 53

  2. Характеризация и операторное представление сплайнов ... 58

  3. Разложение сглаживающего квазисплайна 66

  4. Выбор параметра сглаживания 70

  5. Представление невязки сглаживающего сплайна в виде суммы ряда 78

  6. Сходимость абстрактных интерполяционных сплайнов .... 80

  7. Сплайны на подпространствах: вариационная постановка и сходимость 88

Глава 3. Алгоритмы построения сплайнов 97

  1. Воспроизводящее отображение полугильбертова пространства 100

  2. Структура пространства сплайнов 106

  3. Алгоритм построения смешанного сплайна 109

  4. Построение полиномиального сплайна нечетной степени ... 114

  5. Построение полиномиального сплайна четной степени .... 121

  6. Использование разложения по базису ?-сплайнов 123

  7. Построение jB-сплайнового разложения полиномиального сплайна 127

  8. Построение гетерогенного полиномиального сплайна 130

  9. Построение вариационного L-сплайна 133

Оглавление З

Глава 4. Сплайн-аппроксимация функций многих перемен
ных
142

  1. 1>т-сплайн 145

  2. Сходимость )т-сплайнов 147

  3. Алгоритм построения 1)то-сплайна на подпространстве .... 155

  4. Сплайны Дюшона 158

  5. Алгоритм построения аналитического сплайна 163

  6. Построение аналитического сплайна на вырожденной сетке . 169

  7. Метод наименьших квадратов в сплайн-аппроксимации . . . 177

  8. Интервальная сплайн-интерполяция 180

Глава 5. Сплайн-аппроксимация в тензорном произведении

гильбертовых пространств 194

  1. Тензорные произведения пространств 195

  2. Нормально разрешимые операторы в тензорном произведении пространств 201

  3. Сплайны в тензорном произведении пространств 204

  4. Тензорные произведения функциональных пространств .... 210

Заключение 215

Литература 217

Приложение. Библиотека сплайн-аппроксимации JSpline+ . 225

Описание ядра библиотеки 226

Описание раздела линейной алгебры 227

Описание раздела сплайн-аппроксимации 230

Введение к работе:

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

Принято считать, что теория сплайнов берет свое начало с работы Шёнберга [115], в которой было введено понятие сплайна как функции одной переменной, "склеенной" из кубических многочленов. На основе этого приема были предложены различные варианты аппроксимаций: полиномиальные сплайны высоких степеней, тригонометрические сплайны, L-сплайны. Рассматривались также сплайны с разнородными условиями интерполяции (гетерогенные сплайны) и различными типами краевых условий [12, 27, 69, 70, 117]. Усилиями зарубежных и советских авторов (Алберг, Нилсон, Уолш, B.C. Рябенький, Ю.С. Завьялов, СБ. Стечкин, Ю.Н. Субботин и др.) были получены оценки сходимости таких сплайнов к функциям различной гладкости с привлечением различных типов краевых условий и способов сгущения сетки. Задача одномерной полиномиальной сплайн-интерполяции была обобщена на многомерный случай интерполяции на прямоугольных сетках [1, 82, 27, 40].

Принципиально новый виток развития теории сплайнов начался с работы Холлидея [91], в которой был найден вариационный принцип для кубического сплайна. Аттья в [74, 75] обобщил понятие сплайна, рассматривая его как решение задачи аппроксимации в гильбертовом пространстве, минимизирующее некоторый выпуклый функционал типа полунормы. Новый подход оказался весьма продуктивным. Были получены в общей форме условия существования, единственности и характеризации сплайнов [41]; найдены вариационные принципы для биполиномиальных интерполяционных и сглаживающих сплайнов [24-26, 30]; разработаны алгоритмы построения сплайнов методом воспроизводящих ядер [76, 96]; доказаны общие теоремы сходимости сплайнов [13, 92]; показана связь между сплайнами и оптимальными в смысле Сарда аппроксимациями линейных функционалов [112, 116]; выявлена тесная связь сплайнов с теорії-

Введение

ей поперечников [36] и теорией регуляризации некорректно поставленных
* задач [50].

Следующий важный шаг в развитии теории сплайнов был сделан Дю-шоном [86-89]. Он исследовал задачу построения І)ш-сплайна в Rn с произвольно расположенными узлами интерполяции, а также получил алгоритм и оценки сходимости такой интерполяции. С этого момента берет свое начало теория сплайнов многих переменных на хаотических сетках. Отметим в этой связи работы А.Ю. Бежаева и В.А. Василенко [79-81], О.В. Матвеева [43-46], М.И. Игнатова и А.Б. Певного [29], Пауэла [102], Мэдича и Нельсона [98-100], Шабака и By [114], Лайта и Вейна [95].

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

Для изложения результатов нам потребуются некоторые обозначения
используемые в работе (более полно все необходимые обозначения вводят-
ся в гл. 1, а также частично в других главах). Через X, Y, Z обознача
ем гильбертовы (иногда банаховы) пространства. Используем обычные
обозначения для скалярного произведения и нормы в X: {-,-)х и || * ||х-
Нижний индекс часто опускаем, если понятно из контекста о каком ска
лярном произведении или норме идет речь. Через А Є L(X, Y) обозначаем
линейный ограниченный оператор А, действующий из X bY. Его ядро и
образ обозначаем через N(A) и R(A) соответственно. Важную роль в те
ории сплайнов играют нормально разрешимые операторы, т.е. такие, что
R(A) замкнутое подпространство. Полунорму ||Лж||у, порожденную опе-
,^ ратором А, обозначаем также через ||ж||л- Аналогичные обозначения ис-

пользуем для скалярного произведения. Прямая сумма XY пространств X и Y есть пространство пар вида х фу с покомпонентным сложением и умножением на скаляр. В X У естественным образом вводится р-норма:

ІІ*ву||х.у = (||*|& + ||»||НІ/'.

Обычно используем 2-норму. Прямая сумма операторов А : X —> У" и
В : X -> Z есть оператор AqB'.X-^YQZ, действующий по правилу:
ф В)х = Ах Вх. 2-норма, порождаемая прямой суммой операторов,
(* имеет вид

Введение

Млев = (IMIi + МІГ2 = (ЯМІ + \\Bx\\%f\

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

Гл. 1 вводная. В ней приводится математический аппарат, а также доказываются некоторые факты из смежных областей математики (функционального и выпуклого анализа), используемые в работе. Некоторые из доказанных здесь утверждений будут процитированы ниже.

Теоретические исследования автора содержатся в гл. 2, 4, 5. Практические аспекты иследований автора сконцентрированы в гл. 3, 4.

Начнем с теоретических результатов.

1. Задача смешанной сплайн-аппроксимации изучается в гл. 2. Смешанный сплайн есть решение следующей задачи:

сга = arg min а||Тж||2 + ІІЛ^ж — ^2ІІ2, (0.1)

где Г Є L(X, Y), Ai Є L(X, Zi), A^ Є L(X, Z2), все пространства гильбертовы, z\ Є R(Ai), zi Є ^2) ск > 0 - параметр сглаживания. Данная задача объединяет в себе задачи интерполяции (А\ = А, А2 = 0) и сглаживания (Лі = 0, Л.2 = А). Постановка (0.1) была предложена А.Ю. Бежаевым и В.А. Василенко [80]. В этой же работе приведены условия однозначной разрешимости задачи (0.1) в практической ситуации конечномерности ядра N(T) оператора Т, получено условие характеризации и система операторных уравнений для смешанного сплайна. Автор продолжил эти исследования и доказал теорему о необходимых и достаточных условиях однозначной разрешимости задачи (0.1):

Теорема 2.1.3. Пусть оператор Т нормально разрешим (R(T) замкнуто б Y). Тогда следующие утверждения эквиваленты:

(а) существует оператор В Є L(X, W), действующий в некоторое
гильбертово пространство W, такой, что N(Ai)
С N(B) и норма
Нж1|гвл2 = (||^ж1|2+ Н^ж1|2 + 1И2Сс||2) ' эквивалентна норме ЦягЦх/

(б) смешанная задача (0.1) однозначно разрешима при любых z\ z^ Є
Д(Лі) Z2;

Введение

(в) задача и = arg min а\\Ти —f\\2-h\\A2U — g\\2 однозначно разрешима

uGN(Ai)

при любых f ф g Є Y ф ^2/

(г) на подпространстве N(A\) нормы || ||тел2 w || * ||х эквивалентны.

Данная теорема обобщает результаты П.-Ж. Лорана [41, теоремы 4.4.1-4.4.4] на смешанный случай.

Доказательство этой теоремы, а также других утверждений о смешанных сплайнах базируется на технике сведения задачи (0.1) к эквивалентной задаче наилучшего приближения в выпуклом множестве. В этом заключается отличие предлагаемого метода доказательства от методов, используемых другими авторами (данная методика близка к методике, используемой в работах В.А. Морозова [48, 50]).

В основе этой техники лежит следующая теорема об эквивалентных нормах в банаховых пространствах, доказанная автором:

Теорема 1.3.3. Справедливы утверждения:

(а) если нормы ]| \\ав и || ||х эквивалентны, то N(A) П N(B) — {0}
и N(A) + N(B) замкнуто;

(б) если на подпространстве N(A) нормы \\ \\в и \\ \\х эквивалентны
и R(A) замкнуто, то нормы
|| \\а@в и \\ ||х эквивалентны;

(в) если R(A), R(B) и N(A) + N(B) замкнуты и N(A) П N(B) = {0},
то нормы
|| ||лев и || \\х эквивалентны;

(г) если R(A) замкнуто, N(A) конечномерно и N{A)C\N(B) = {0}, то
нормы
|| \\ав и || ||лг эквивалентны.

2. Ясно, что любой сглаживающий сплайн является интерполяционным для некоторого другого вектора измерений z. В монографии П.-Ж. Лорана [41, теорема 4.6.4] рассматривается и обратная задача: пусть имеется интерполяционный сплайн и задано а > 0; существует ли вектор 5, для которого этот сплайн является сглаживающим. Обозначая через А+ и А+ операторы, сопоставляющие вектору z Є R(A) интерполяционный и сглаживающий сплайны соответственно, данную задачу можно сформулировать так: совпадают ли образы операторов А+ и А+? Ответ на этот вопрос

Введение

положительный. Множества интерполяционных и сглаживающих сплай-
'fc нов при фиксированном А совпадают и образуют пространство сплайнов

Spl(.A) С X, характеризующееся условиями

VueN(A) (Тх,Ти) = 0.

Автор ставит аналогичный вопрос об операторе А, сопоставляющем
вектору z = ziZ2 Є R(Ai) ф Z2 сплайн аа - решение задачи (0.1). Ответ

на этот вопрос дает следующая теорема, предлагающая также способ продолжения оператора A + на все пространство Z\ ф Z2\

Теорема 2.2.10. Пусть R(T) и R(Ai) замкнуты и задача (0.1) однозначно разрешима при любых допустимых исходных данных. Обозначим А = Аі ф А2. Тогда:

(а) A* : R(A\) ф Z2> Spl(A) - линейный ограниченный оператор;

(б) оператор А+ удовлетворяет тождеству

to-

A+z = A+Vz, (0.2)

где V - ортопроектор на R(Ai) Ф R(A2), z = z\ Ф z2 Є і?(^4і) Ф Z2;

(в) оператор А непрерывно продолжается по формуле (0.2) на все пространство Z = Z\ 0 ^2, nmo соответствует замене задачи (0.1) wa задачу

да = arg min «ll^dl2 + \\А2х — z2\\2',

W (г) если (Т, Л) - сплайн-пара и операторы А\ и А2 линейно независимы

(R(Ai 0 А2) = R{Ai) ф R(A2)), то сужение оператора А+ на R(A) есть непрерывный изоморфизм между R(A) и Бр\(А).

Пару (Г, А) здесь и далее называем сплайн-парой, если R(T), R(A) и N{T) + N(A) замкнуты, N(T) П N(A) = {0}.

Пункт (г) этой теоремы можно выразить словами так: если (Г, А) -
сплайн-пара и операторы А\ и А2 линейно независимы, то любой элемент
пространства сплайнов Spl(-A) является смепіанньїм сплайном при любом
,* а > 0 и подходящем (однозначном) выборе z Є R{A).

Введение

«є

3. Задача смешанной сплайн-аппроксимации изучалась автором также в контексте сплайнов на подпространстве. В.А. Василенко назвал сплайн, построенный методом конечных элементов, сплайном на подпространстве [14].

Пусть {^г}г>о - семейство замкнутых подпространств в X. Смешанным сплайном на подпространстве ЕТ называется решение задачи

oj = arg min а\\Тх\\2 + \\А2х - z2\\2. (0.3)

Здесь предполагается, что Аї1{х{) П Ет ф 0, т.е. условия интерполяции А\х = z\ не противоречивы на Ет.

Введем эквивалентную норму || ||* = || \\у/бтва2 п0 теореме 2.1.3, возьмем ортопроектор Рт на Ет в норме || ||* и введем оператор сплайн-интерполяции S:

Su = arg min а||Тж||2 + \\А2х\\2.

Axx—A\u

Автор доказывает следующую теорему сходимости сплайнов <т„ к сга:

Теорема 2.7.9. Пусть подпространства ЕТ предельно плотны в X при т > 0 и найдется такое tq > 0, что при т < 7 R(AiPT) = R(Ai) и выполнено одно из условий:

(а) dim R(A\) < со;

(б) ||(7-Рт)5||.<С7<1.

Тогда ||5 — аа\\х при т -> 0.

Доказательство этой теоремы базируется на доказанной автором теореме о сходимости наилучших приближений в замкнутых выпуклых множествах. Пусть М С X - замкнутое выпуклое множество и Рм ' X —» X -оператор, сопоставляющий каждому элементу / Є X наилучшее приближение h Є X на множестве М по формуле Рм/ = Л, = argmin^M \\х — /||2.

Теорема 1.6.1. Пусть Мі С X - непустые замкнутые выпуклые множества и fi X - некоторые элементы, і Є IN. Предположим, что fi—tf при г —> оо и множества Мі удовлетворяют одному из условий:

Введение

(а) Mi+1 С Мі и П Мг = М;

»1N

(б) Мі С М и существует последовательность Х{ Є М», сходящаяся к

Тогда последовательность hi = PmJi сильно сходится к h — Рм/, т.е. \\Ы — h\\x -) 0 при г -> оо.

4. Следующая тема исследований автора - выбор параметра сглаживания а при построении сглаживающего сплайна

аа = argmin а||Тж||2 + \\Ах - zf. (0.4)

Задача сглаживания исследовалась многими математики, начиная с Райн-ша [103, 104]. Мы ограничимся здесь только исследованиями поведения сглаживающего сплайна при изменении а и алгоритмом, базирующемся на критерии невязки (см. монографию В.А. Морозова [50]). Работы автора в этом направлении [108, 61, 63] содержат несколько новых результатов. Предположим, что (Г, А) - сплайн-пара, и введем скалярное произведение

(и, v)„ = {щу)ТфА = (Tu,Tv)Y + {Au,Av)z,

порождающее эквивалентную норму в X. Представим пространство X в виде

X = N(A) 0 N(T) Х, Х = (N(A) + N(T))±, (0.5)

где ортогональное дополнение берется относительно скалярного произведения (, )*.

Теорема 2.3.3. Пусть (Т, А) - сплайн-пара и операторы А*, Т* сопряжены к А, Т относительно скалярного произведения (, )*. Тогда:

(а) операторы А*А и Т*Т перестановочны, причем А*А + Т*Т = їх;

(б) разложение (0.5) приводит операторы А*А и Т*Т к виду

А*А(щ + и2 + из) = ОщА+ 1щт)Щ + Ли3, Т*Т(щ + и2 + щ) = А)Щ + ОщТ)и2 + Тщ,

где щ Є N(A), и2 Є N(T), щ Є X;

Введение

(в) А,Т Є L(X) -перестановочные, самосопряженные, положительно определенные изоморфизмы, А-\-Т' = їх, \\Л\\ < 1, \\Т\\ < 1.

Здесь через и Оу обозначены тождественный и нулевой операторы bV.

Теорема 2.3.4. Рассмотрим задачу построения сглаживающего квазисплайна

аа = argmin а\\Тх — у\\2 + \\Ах z\\2.

Положим

o-qo = are min \\Ax z\\2, сто = arg min \\Tx - y\\2.

xu(A*A)-l(A*z)

Тогда: gq — (7^ Є X; oa = qa = gq + ra, причем qa и ra принадлежат X и определяются из уравнений

(I + aA~lrr)qa = qo = o-Q- a^, (I + a~lT~lA)ra = Гоо = о-,» - cr0.

Из теоремы 2.3.4 легко выводится сходимость сглаживающего квазисплайна сга к предельным элементам ctq и егдо с оценками 0(a) и 0(а~1) соответственно, причем при ctq ф (Too эти оценки не улучшаемы по порядку [63]. Факт сходимости квазисплайна оа к предельным элементам хорошо известен (см., напр., [50]). По-видимому, была ранее известна и оценка сходимости сплайна о~а к <то со скоростью 0(a), хотя в [17, 50] получена только оценка 0(а112).

Один из методов выбора параметра сглаживания в задаче (0.4) заключается в решении уравнения <р(а) = \\Ааа — z\\ = є относительно а. Это уравнение называют критерием невязки. В работе В.И. Гордоновой и В.А. Морозова [23] предлагается выбирать параметр сглаживания, решая эквивалентное уравнение

ГІР) = є', (0.6)

где /3 = а-1, ф(/3) = <р(а). Если ^(0) ф VK00)? то ф3(Р) при s > 0 -строго монотонно убывающая выпуклая вниз функция, при — 1 < s < 0 -

Введение

строго монотонно возрастающая выпуклая вверх функция, и скорость сходимости метода Ньютона максимальна при s = — 1.

Автору удалось получить формулу шага метода Ньютона для решения уравнения (0.6) при s = — 1 в терминах оператора невязки Raz — z Ао~а:

Теорема 2.4.4. Итерация метода Ньютона для решения уравнения (0.6) при s = — 1 имеет вид

1-и>(аА) , . {Raz,R2az) (Raz,R2az)

(р(оск)/є-ш(ак) (Raz,Raz) (pl{a)

Данная теорема дает возможность построения универсального алгоритма выбора параметра сглаживания. Этот алгоритм реализован автором в библиотеке программ JSpline+ [93].

Интересное применение имеет полученная автором формула дифференцирования оператора невязки Ra:

Лемма 2.4.3. Для оператора невязки Ra Є L(Z), задаваемого правилом Raz — z — Acra = (I — AA*)z, имеет место следующая формула дифференцирования по параметру а:

DkRa = (_i)*+i^Д(/ - Ra), к > 1.

На основании этой леммы с привлечением теоремы 2.3.3 автору удалось доказать следующую теорему о разложении невязки:

Теорема 2.5.1. Пусть «о > 0. Тогда невязка z Ааа представима в виде суммы ряда

[k=0\ а0 і абсолютно сходящегося на отрезке [0,2«о]

Ааао, (0.7)

Похожее разложение (только самого сплайна аа, а не невязки) рассматривалось ранее А.Ю. Бежаевым. Так в [80, Theorem 12.10] доказана сходимость разложения сплайна сга в ряд Тейлора, но на открытом интервале (0,2ао). Автору в теореме 2.5.1 удалось расширить интервал сходимости ряда для Ааа и явно представить разложение по степеням оператора Ra.

Введение

Автором также исследовано поведение сплайна а вблизи бесконечно-
щ сти. Полагая Qp = I Ri/p, получаем точно такие же правила дифферен-

цирования оператора Qp по j3 как в лемме 2.4.3, из которых, с привлечением теоремы 2.3.3, выводится следующая

Теорема 2.5.3. Пусть / > 0. Тогда невязка Z ~~ Aaijp представимо, в виде суммы рлда

* г - Аат = \ Е (1^] (* - М/а), (0-8)

Ц=оч A) J \

абсолютно сходящегося на отрезке [0,2/].

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

шаге метода.

Теоремы 2.5.1, 2.5.3 дают возможность построения более эффективного алгоритма выбора параметра сглаживания за один шаг: выбираем начальный коэффициент сглаживания ао; выполняем один раз построение и декомпозицию матрицы сплайновой системы; анализируем, с какой стороны от решения находимся и в зависимости от этого используем разложение (0.7) или (0.8) для поиска параметра сглаживания. При этом потребуется рассчитать начальный отрезок достаточной длины от соответствующего ряда, получить формулу вычисления невязки (квадрат

* которой есть полином, зависящий от а или /3) и приближенно найти ре
шение получившегося уравнения. Эта же идея может быть использована
в итерационных алгоритмах построения сглаживающих сплайнов.

5. Задачи сплайн-аппроксимации функций многих переменных

представляют как теоретический, так и практический интерес. Теоре
тические вопросы включают в себя проверку корректности (однознач
ной разрешимости) задачи сплайн-аппроксимации, доказательство сходи
мости интерполяционных сплайнов и получение оценок сходимости. Ав-
Л" тор ограничивается только случаем, когда аппроксимируемая функция

Введение

принадлежит пространству, в котором ищется сплайн. Более общие конструкции (продолжение оператора сплайн-интерполяции в банаховы пространства и оценки сходимости) рассматривались, например, в работах О.В. Матвеева [43, 44].

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

Построение сплайна методом конечных элементов (сплайны на подпространствах в терминологии В.А. Василенко). Здесь различаются два варианта построения конечных элементов: когда носители конечных элементов строятся на основе сетки узлов интерполяции и когда носители не связаны с сеткой. Второй вариант предпочтительнее, поскольку мы не можем выбирать сетку узлов интерполяции, а значит не можем гарантировать качественную триангуляцию области определения функции если используем сетку узлов интерполяции как основу для построения конечных элементов (в этом отличие МКЭ для сплайнов от МКЭ в задачах математической физики).

Аналитические методы, использующие представление сплайна через воспроизводящее ядро полугильбертова пространства (X, || \\т)-Пример таких сплайнов - сплайны Дюшона. Ограничения данного подхода: воспроизводящее ядро известно только в частных случаях.

Построение сплайна на регулярной сетке, являющейся прямым произведением одномерных сеток. Это сплайн-аппроксимация в тензорном произведении пространств.

6. )т-сплайны в произвольной области представляют в основном теоретический интерес. Задача >т-аппроксимации обычно рассматривается либо в области с липшицевой границей, либо в Яп. Автор в [58] распространил постановку задачи Х)ш-аппроксимации на произвольную область в Rn. При этом в качестве X надо брать L(Q) - пространство функций в ГІ, m-e частные производные которых интегрируемы с квадратом.

Доказательство сходимости )т-сплайнов в произвольной области получено .автором. Ранее В.А. Василенко было доказано [80], что если область О, ограничена и имеет липшицеву границу, то при произвольном

Введение

способе сгущения конечных є-сеток в fi 1)т-сплайны, интерполирующие функцию / Є W^fi), сходятся к / в норме пространства W^fi). Автору удалось убрать все ограничения, выделенные в предыдущем предложении курсивом, и доказать следующую теорему:

Теорема 4.2.3. Пусть последовательность множеств {ш{ С и? П S1}jgbst предельно плотна в и, т > п/2 и f Є L^) - интерполируемая функция. Тогда последовательность интерполяционных сплайнов 0{ — SWJ сходится к сплайну a = Suf в норме пространства Ь(П).

Здесь Suоператор, сопоставляющий функции / Є 1^(0.) интерполяционный )т-сплайн а с условиями интерполяции cr(t) = f(t), t Є w.

Данная теорема есть следствие другой, более общей теоремы, доказанной автором в [107, 61]:

Теорема 4.2.2. Пусть fi С Rn ~ локально-компактное множество, X непрерывно вложено в CioC(fi) при некотором 0 < а < 1. Предположим, что оператор Т Є L(X, Y) имеет конечномерное ядро и замкнутый образ; множество w С П содержит L-набор точек для N(T); последовательность множеств {ш{ С ШПГ2}іЄдо предельно плотна в и. Тогда найдется такое число г'о > 0, что при і > г'о задача

о і = argmin{||Tz||2, х Є X, x(t) = f(t), t Є 0}

однозначно разрешима и последовательность {оу} сходится в норме пространства X к

а = argmin{||T:r||2, х Є X, x(t) = /(і), t Є и}

при г —> со.

Здесь Cfoc($l) - пространство непрерывных функций в fi, локально удовлетворяющих условию Гельдера с параметром 0 < а < 1. Множество fi С Rn назовем локально компактным, если для любой точки t Є Q найдется компактная относительная окрестность fit этой точки в П. В свою очередь, множество fit С fi называют относительной окрестностью точки t в П, если существует множество Dt С Rn - окрестность точки t в Rn, такое, что fit = DtC\ fi.

Введение

Замечание. Замкнутые и открытые множества в Rn являются локально
«I компактными. Автор в 1.7 дает эквивалентные определения локально-

компактного множества и доказывает, что пространство Cioc(^) функций, непрерывных в П, является пространством Фреше.

Теорема 4.2.2, в свою очередь, выводится из теоремы о сходимости

интерполяций на сгущающемся семействе функционалов, усиливающей

теорему В.А. Василенко о сходимости на сгущающихся є-сетях (ср. [80]).

Пусть Г2 С Rn - некоторое множество, и задано параметрическое семейство функционалов Ф(П) = {(pt Є X', t Є ГХ}. Далее, пусть cv -некоторое подмножество Q, и / 6 X - произвольный элемент. Положим

M(f,u>) = {х Є X : cpt{x) = t{f), t Є и)}.

Автором доказана в [107, 61] следующая

Теорема 2.6.7. Пусть N{T) конечномерно, множество Ct локально-компактно и семейство функционалов Ф(П) непрерывно по параметру, т. е. \\s (pt\\x' -* 0 при s —* і) s,i Є П. Пусть также задача

а = arg Л»? ч игя?и2

однозначно разрешима.

Тогда если последовательность множеств {ші С о7ПП}геж предельно плотна в ш, то существует такое г'о Є IN, что при і > г'о задача

U{ = arg min ||Тя;||2

однозначно разрешима, и ||ст,- — о~\\х —> 0 npw і —> оо.

7. Оценки сходимости Г>т-сплайнов в ограниченной области с лип-шицевой границей были получены А.Ю. Бежаевым в [4] в виде

\\Dkx\\Lp(n) < Chm-k-n^n^\\Dmx\\L2{ii) (0.9)

при р>2 ж т — к — п/2 -f п/р > 0 (т — к — п/2 > 0 при р = оо). Здесь

х Є Ь(0,) - произвольная функция, имеющая h-сеть нулей в П и констан-

? та С > 0 не зависит от ж и выбора /і-сети. Автор в [58] распространил эти

Введение

оценки на ограниченные области, удовлетворяющие слабому условию конуса (в частности, на области с разрезами, выколотыми точками, острыми пиками, направленными внутрь), а в [61] - на (возможно неограниченные) области являющиеся объединением конечного числа областей типа Ь (удовлетворяющих слабому условию конуса и допускающих продолжение за пределы области с сохранением класса гладкости). В работах О.В. Матвеева [43, 44], получен более широкий спектр оценок сходимости І)то-сплайнов в различных метриках, однако рассмотрение ограничивается только областями, удовлетворяющими сильному условию конуса.

8. Сплайны Дюшона. Французский математик Дюшон предложил в [87] сплайны, являющиеся обобщением -От-сплайнов в ^(Rn). Он рассмотрел пространство D~mHr(Rn), состоящее из обобщенных функций, ттг-е производные которых принадлежат "почти Соболевскому" классу ^ГП) функций умеренного роста с ограниченной полунормой

ч1/2

Ы\нг=(!Кп\т\\Гх\*<1т) ,

которая при г < п/2 является нормой. Здесь Т - преобразование Фурье. Воспроизводящее ядро пространства D~mHr{Rv) относительно полунормы || |І)ш есть радиальная базисная функция

\s — |27ln|s — t\, 7 - целое,
\s — t\27 иначе,

где 7 = m + r n/2, 0 < 7 < ra, [7] - целая часть 7, a

1/2

G7(M) = (-1)W+1 <

(n \

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

- /)||ь,(П) < Ch-k-»'2+n'*\\h - f\\wm,

если П - ограниченная область, удовлетворяющая слабому условию конуса, s > п/2, s - к - п/2 + п/р > О, р > 2 [114, 95].

Введение

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

изведении сеток. Это хорошо известные задачи биполиномиального восполнения [1, 24-27, 30]. Такие сплайны, называемые тензорными, достаточно легко строятся: алгоритм редуцируется к сериям одномерных задач. Наряду с задачами, имеющими тензорную природу (параллелепипедаль-ные сетки узлов интерполяции), можно рассматривать и задачу в тензор-ном произведении пространств с произвольной сеткой узлов интерполяции. Для этого необходимо построить нормально разрешимый оператор, чтобы использовать его в энергетическом функционале.

Следуя В.А. Морозову [48], назовем операторы А и В взаимно дополнительными, если оператор А Ф В порождает эквивалентную норму в X. Будем также называть оператор А дополнительным к В (В дополнительным к А), если А и В взаимно дополнительны.

Рассмотрим тензорное произведение гильбертовых пространств X = Xi.. п. Автор доказал следующую теорему о нормально разрешимом операторе в X:

Теорема 5.2.3. Пусть Xi, Yi, W{ - гильбертовы пространства, Ті Є L(Xi,Yi), Ві Є L(X{,Wi) - некоторые операторы, і — 1,...,п. Если операторы Ті нормально разрешимы и операторы В{ дополнительны к Ті, то оператор Т = [(Ті фВ\) ... пфВп)] Q [В\ <8>... Вп] нормально разрешим и N(T) = N(Ti) ... N(Tn).

Здесь операция G означает, что из прямой суммы тензорных произве
дений операторов, получающейся в результате раскрытия скобок, следует
** исключить компоненту В\ ... В
п.

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

Операторы Т Є L(X,Y) и Б Є L(X,W) назовем строго взаимно дополнительными, если они взаимно дополнительны и линейно независимы. Это означает, что пространство X разлагается в прямую сумму

Введение

N(T) + N(B), а операторы T и В нормально разрешимы.

Пусть 7 _ воспроизводящее отображение пространства (X, || \\т) и тг -воспроизводящее отображение пространства (X, || ||в). Возьмем проектор Р Є 1/(Х) со свойствами: N(P) = N(T), R(P) = N(B). Построим двойственный ему оператор Р' Є L(X') по правилу (Р'(р)(х) — (р(Рх), х Є X, Є X'. Рассмотрим также проектор Q = І—Р и двойственный ему оператор Q'. Через U0 обозначаем аннулятор подпространства U С X, т.е. подпространство линейных функционалов из X', равных нулю на U.

Автор доказывает следующую теорему о регуляризованном воспроизводящем отображении:

Теорема 5.3.5. Регуляризованное воспроизводящее отображение -у = P'jP' + QttQ' удовлетворяет следующим условиям:

(а) \/tp Є N(T)\ хЄХ (Вх,В^уср) = 0;

*

(б) У ер Є N(B), хеХ {Тх, Tjip) = 0;

(в) -у ~ воспроизводящее отображение (X, || \\т)\
/
(г) 7 ~~ воспроизводящее отображение (X, || ||в);

(д) 7 ~ воспроизводящее отображение (X, || ||тв);

(е) отображение непрерывно, симметрично и действует на все про
странство
X.

Пусть теперь Xi, Yi, W{ - гильбертовы пространства, Т{ Є L(Xi,Yt)
и В{ Є L(Xi, W{) - строго взаимно дополнительные операторы, 7» Є
1»(Хг-, Хг) - регуляризованные воспроизводящие отображения относитель-
^ HO (Xi, || ||т,, || ||в.), г = 1,. . . , П.

Теорема 5.3.7. Отображение у = 71 7п является воспроизводящим отображением пространства X = Xi Хп относительно нормы, порожденной оператором Т — (Ті jBi) ... п ф 5n), гі относительно полунормы, порожденной оператором Т = Т в [#i ... J5n].

На этом собственно круг теоретических исследований автора заканчи
вается. Перейдем к практическим исследованиям, посвященным алгорит-
» мам построения сплайнов.

Введение

10. Центральное место в практических исследованиях автора занимает алгоритм построения смешанного сплайна с весами. Пусть А\ =

(fi ф ... Ф <рмj -А-2 = Ум+1 Ф Ф и функционалы у>г Є -X7, г = 1,..., N, линейно независимы. Предположим, что образ оператора Т замкнут, ядро конечномерно с базисом щ,... , гг#- и операторы Т и А = А\ ф Аг порождают эквивалентную норму в X.

Пусть гильбертово пространство X непрерывно вложено в C\oc(i), О С Rn, G(s,t) - воспроизводящее ядро пространства X относительно полунормы || Цу. Рассмотрим задачу

>і(б-а) = 2ї, г = 1,..., М,

i=M+l Х

построения смешанного сплайна с весами рг > 0. Решение этой задачи имеет вид

N К

^a(s) = 53 ^mG{s, -) + 53 /W(s).

г=1 г=1

Векторы неизвестных коэффициентов Л = (Лі,..., \n)T> № = (а*ъ ---5 А4*:)7 находятся из невырожденной системы линейных уравнений

С + аР U\ (\\ (z\ ,плп.

* о) („)-(о)- (ai0)

где С ж U - NxN- и ЛГхііГ-матрицьі с коэффициентами

cij = iptщк = <рі{щ), г, j = 1,..., АГ, к - 1,..., К

[ірі действует на G по переменной s, a диагональная матрица с коэффициентами

0 при і = 1,..., М,

Рг При І = М + 1, . . . , N.

Отметим, что коэффициенты вектора Л не зависят от выбора воспроизводящего ядра G (хотя это ядро и неединственно, но при замене его на другое воспроизводящее ядро в решении изменяются только коэффициенты вектора /і).

Введение

Систему уравнений (0.10) можно преобразовать к системе меньшего
Л размера с симметричной положительно-определенной матрицей, если по-

строить некоторый линейный оператор Н : RN —> RN~K так, чтобы выполнялось условие N(H) = R(U).

Теорема 3.3.6. Смешанный сплайн аа имеет вид

N-К К

1 і=1 г=1

где фі = T,f=i hij(fj, і = 1,..., N К.

Вектор v находится из системы уравнений

Aav = Hz (0.11)

с симметричной, положительно-определенной матрицей Аа = А 4-
аНРН*, причем матрица А = НСН* симметрична и положительно
определена. Коэффициенты матрицы А не зависят от выбора воспро-
^
изводящего ядра G и задаются формулой а^ = (фіЄ,^С)т = (wi,Wj)y,

причем функции W{ = TtJjiG Є Y также не зависят от выбора воспроизводящего ядра.

Сплайн аа удовлетворяет условиям интерполяции с вектором za = z — ctPH*v, а вектор неизвестных коэффициентов /2 можно найти из условий интерполяции Ааа = za и представления Таа = T.f^{K fi^i, либо из системы уравнений U/л = z — (С + aP)H*v.

Преобразование системы (0.10) к виду (0.11) хорошо известно. Обычно в качестве Н берется разреженная матрица конечных разностей, постро-

ф- енных на достаточном количестве L-наборов узлов. Еще используются

построения с шаблонами, число узлов в которых на единицу больше [101]. Автор предложил иной метод, основанный на чисто алгебраических преобразованиях системы (0.10) к виду (0.11). Этот метод был реализован Д.С. Кротовым [38] под научным руководством автора. Отметим также, что аналогичный подход практически в то же время был предложен Ша-баком [113].

Идея предложенного автором алгоритма заключается в построении матрицы Н в виде J2Q, где матрица Q приводит матрицу U к верхне-тре-

^ угольной матрице с помощью, например, ортогональных преобразований,

Введение

а матрица J2 понижает размерность вектора, отбрасывая его первые К
щ
компонент (см. 4.5). Отметим, что использование ортогональных пре-

образований оптимально в смысле сопсЬ, а именно, число обусловленности сопсЬА матрицы А совпадает с числом обусловленности оператора ТС\щщі-, являющегося сужением оператора VC на подпространство R(U)L, где V - ортопроектор на Я(С)Х [62].

В практических задачах трудно гарантировать, что узлы сетки распо-
^ ложены "хорошо". Может возникнуть ситуация неединственности регле-

ния из-за вырождения сетки (например, сетка узлов попала в некоторое собственное афинное подпространство в И„). Автор в [62] распространил предложенный алгоритм на случай вырождения сеток (см. п. 4.6.4). Для этого он разработал алгоритм модифицированного Q-R-разложения, позволяющий вычислять нормальное псевдорешение систем уравнений с прямоугольной матрицей. Этот алгоритм реализован под руководством автора в библиотеке JSpline+ [93].

$" 11. Классические методы построения интерполяционного, сглаживающе-

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

Интересный подход к решению больших задач использовался А.Ю. Бе-жаевым в комплекте программ LocalGreen в [97]. Исходное "облако" точек с помощью последовательной дихотомии гиперплоскостями разбивается на перекрывающиеся подмножества приемлемого размера (порядка 1000 узлов вместо исходных сотен тысяч), в которых осуществляется сплайн-аппроксимация. Затем полученные решения в подобластях "склеиваются" с помощью разбиения единицы, т.е. выполняется композиция сплайнов (оценки сходимости метода композиции сплайнов получены Д.С. Кротовым [37] под руководством автора диссертации). Отметим также эксперименты М.И. Игнатова и А.Б. Певного [29] с регулярным разбиением множества узлов на подмножества для двумерной задачи.

Отметим также другую специфическую черту больших задач: измере-

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

Введение

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

низкая, поэтому в базис часто добавляют радиальные функции (нейрон
ные сети типа РБФ [51]). Функцию G(s,t) называют радиальной, если
G(s, t) = g(\s 1\). На функцию G(s, t) дополнительно накладывается тре
бование условной положительной определенности [102]. Это означает, что
G есть воспроизводящее ядро некоторого полугильбертова пространства
*~ функций [98, 99]. Например, в качестве G можно брать воспроизводящее

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

Идея использования метода наименьших квадратов в сплайн-аппро-оксимации не нова. Она предлагалась еще де Бором [83] для приближения функции одного переменного с помощью линейной комбинации В-сплайнов. Такой способ дает приемлемые результаты, если в сетке узлов измерений нет больших пропусков. В противном случае в промежутках между узлами интерполяции будут "провалы", что приведет к плохим результатам.

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

Формально, пусть в Rn задана сетка узлов интерполяции ujs = {s^Li и некоторая сетка узлов LOt = {tj}f=i- Будем искать сплайн

a(s) = XjGis, tj) + m(s), (0.12)

удовлетворяющий условиям ортогональности

Ё XjUiitj) = 0, і = 1,...,К, (0.13)

и минимально квадратично уклоняющийся от условий интерполяции
щ) a(si) = Z{, і = 1,..., М. (0-!4)

Введение

Здесь щ - базисные функции N(T).

Данная задача есть задача наименьших квадратов с линейными ограничениями. Для ее решения поступаем следующим образом. Построим оператор Н : RN RN~K в виде J2Q, удовлетворяющий условию N(H) = R(Ut), Ut = Ш). Воспользовавшись заменой переменных Л = H*v, сводим данную задачу к задаче наименьших квадратов без ограничений в пространстве векторов v Є RN~K Ф RK. Последнюю задачу решаем, воспользовавшись Qil-разложением с накоплением (см. п. 4.7.1).

12. Другой подход к решению задач с неточными данными заключается в построении сплайна, проходящего вблизи интерполируемых значений: в узлах сетки задаются доверительные интервалы, через которые должно проходить решение. Это задача в выпуклом множестве, изучавшаяся П.-Ж. Лораном [41]. Он получил необходимые и достаточные условия характеризации решения. Основная особенность решения такой задачи заключается в разделении множества узлов на "узлы прилипания", в которых решение "прилипает" к верхнему или нижнему ограничению, и свободные узлы, в которых решение проходит строго между ограничений. Если некоторый узел сетки свободный, то соответствующий ему коэффициент в представлении (0.12) равен нулю. Таким образом, в представлении сплайна участвуют базисные функции только в узлах прилипания, и можно надеяться получить решение, в котором задействовано существенно меньше узлов сетки. Подход, основанный на поиске узлов прилипания сплайна, был реализован А.В. Ковалковым для одномерного сплайна (см. [18]). Аналогичные идеи были запрограммированы М.И. Игнатовым для сплайна многих переменных (см. [29]).

Пусть X = Х($1), Y = У (О) - гильбертовы пространства функций в ГІ С Лп, Г Є L(X, Y) - нормально-разрешимый оператор с конечномерным ядром и X непрерывно вложено в Cioc(^). Пусть в узлах конечной сетки ш = {ii,..., t^} С ft заданы доверительные интервалы [zf, zf], zj < zf, і — 1,..., N, в которые должен попадать сплайн.

Задачу

а = arg min ||Г*||2, (0.15)

агЄЛ/ц,(г-+)

Mw(z~,z+) = {x Є X : zf < x(U) < zt, і = 1,..., N},

Введение

назовем задачей интервальной сплайн-интерполяции.

# Как было отмечено выше, интервальный сплайн характеризуется уз
лами прилипания к ограничениям. П.-Ж. Лоран получил условия един
ственности решения, в которых существенно, что сетка узлов интерполя
ции удовлетворяет условию Хаара на N(T). Для полиномов многих пере
менных это условие обычно не выполняется. По этой причине, единствен
ность решения при построении интервального сплайна многих перемен-

** ных гарантировать невозможно.

Решение задачи (0.15) имеет вид (0.12), (0.13), причем справедливо следующее: если o~(ti) = zf, то Aj > 0; если а(и) — zf, то А* < 0; если zf < ст(и) < zt, то А; = 0 [41].

Узлы сетки, для которых Af ф 0, назовем активными, а узлы сетки, в которых о{и) — z~ или а(и) = л*, назовем узлами прилипания. Ясно, что сетка узлов прилипания содержит в себе сетку активных узлов.

Следующие теоремы доказаны автором:

* Теорема 4.8.2. Пусть сетка uj содержит L-набор для N(T). Тогда
множество S
u(z~,z+) решений задачи (0.15) есть замкнутый выпуклый
ограниченный многогранник в пространстве сплайнов
Spl(a;), параллель
ный N(T), т. е. если сгі,сг2
Є Su(z~,z+), то Є N(T).

Теорема 4.8.3. Пусть сетка и содерэюит L-набор для N{T). Тогда для любых 0"і,о"2 Є Stj(z~,z+) множества активных узлов сплайнов о~\ и &2 совпадают. Обозначая через ш - множество активных узлов, имеем cri(t) = сг2() для всех t Є wac.

Если к тому же множество шас содержит L-набор для N(T), то
Ч
решение задачи (0.15) единственно.

Теорема 4.8.4. Пусть сетка и) содержит L-набор для N(T). Тогда найдется такое решение а Є Sw{z~,z+), что его множество узлов прилипания о;с1(<т) содержит L-набор для N(T).

Теорема 4.8.4 дает обоснование алгоритма поиска решения путем пе
ребора узлов прилипания сплайна. Отметим первый шаг предлагаемого
в п. 4.8.5 алгоритма, отсутствующий в алгоритмах других авторов: по-
,* строение начального приближения, в результате которого либо находит-

Введение

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

сетки узлов прилипания также отличается от методов, использованных А.В. Ковалковым (см. [18]) и М.И. Игнатовым (см. [29]). Вместо того, чтобы сразу искать подсетку узлов прилипания для сплайна с заданными доверительными интервалами (и рисковать получить решение с числом узлов, сравнимым с числом узлов сетки), предлагается сначала ослабить

* ограничения, а затем постепенно их сжимать. При этом можно контроли
ровать число узлов прилипания и строить более сложные критерии оста
нова процесса приближения.

Метод решения промежуточной задачи квадратичного программирования на подсетке также отличается от методов, использованных другими авторами. Решение задачи ищется методом проекции градиента. Автор получил крайне простые формулы вычисления сплайна, являющегося проекцией градиента. Разобьем подсетку й, на которой ищем решение промежуточной задачи, на три подсетки: подсетка u)q состоит из узлов прилипания, в которых условия характеризации предыдущего приближения выполняются; подсетка а>і состоит из узлов прилипания, в которых условия характеризации предыдущего приближения нарушаются; подсетка 0)2 содержит остальные узлы подсетки ш.

Теорема 4.8.7. Сплайн

tju! k=l

коэффициенты которого Аф/t задают вектор проекции градиента, есть решение задачи сплайн-интерполяции

о = arginin{||Tz||2, x(tk) = 7*, tk Є a)},

Г 0 при tk Є шо,

[ -А& при tk Є wi U a>2-

Здесь Afc — коэффициенты предыдущего приближения. Предложенный алгоритм построения сплайна был реализован

* И.А. Молотковой [47] под руководством автора диссертации. Описание

Введение

данного алгоритма завершается рассмотрением случая существенной не
фі единственности решения, когда вся сетка узлов не содержит L-набор. При
водится модификация алгоритма для этого случая.

13. В приложении к диссертации приводится описание библиотеки сплайн-аппроксимации JSpline+ [21, 93], созданной под руководством и при непосредственном участии автора студентами НГУ А.В. Галковым,

^ О.А. Лихачевым, А.Е. Никишкиной, Н.Ф. Фурсовой, а также аспирантом

ИВМиМГ СО РАН Д.В. Петраковым.

В настоящий момент в библиотеке реализовано большое количество алгоритмов линейной алгебры, в том числе модифицированное QR-разло-жение в трех вариантах, Q-R-разложение с накоплением, решение системы линейных уравнений с блочной матрицей вида (0.10) путем сведения ее к системе (0.11), решение системы (0.10) с произвольной правой частью. В разделе сплайн-аппроксимации реализованы алгоритмы одномерной полиномиальной сплайн-аппроксимации сплайнами четной и нечетной

* степеней. Для этих сплайнов имеются режимы интерполяции, сглажива-

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

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

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

Введение

#

Математические утверждения (теоремы, леммы и т.д.) не имеют собственной нумерации. При ссылке на них используются номера пунктов, в которых эти утверждения приведены. Часто формулировки теорем состоят из нескольких утверждений, имеющих буквенную нумерацию (а), (б) и т.д. При ссылке на конкретное утверждение такой теоремы к номеру пункта добавляется соответствующая буква. Так, например, фраза "теорема 1.2.4 (в)" означает ссылку на утверждение (в) теоремы, приведенной в п. 1.2.4. Нумерация формул двойная и не зависит от разбиения текста по параграфам. Конец доказательства помечается символом , прижатым к правому краю страницы.

Подобные работы
Милосердов Владимир Владимирович
Дискретно-стохастические численные алгоритмы со сплайн-восполнениями
Медведева Мария Викторовна
О вложениях классов функций Н w в классы функций ограниченной обобщенной вариации и некоторые вопросы теории аппроксимации функций
Бирюков Лев Николаевич
Распределение нулей функций из обобщённых пространств Бергмана и некоторые применения в теории аппроксимации
Сорокин Владимир Николаевич
Аппроксимации Эрмита-Паде и обобщенные системы Никишина в теории диофантовых приближений и в теории динамических систем
Сорокин Владимир Николаевич
Аппроксимации Эрмита-Паде и обобщенные системы Никишина в теории диофантовых приближений и в теории динамических систем
Лазарев Александр Алексеевич
Методы и алгоритмы решения задач теории расписаний для одного и нескольких приборов и их применение для задач комбинаторной оптимизации
Садыков Руслан Равильевич
Алгоритмы решения задач теории расписаний для одного прибора с критериями Lmax и EwjUj
Чехович Юрий Викторович
Элементы алгебраической теории синтеза обучаемых алгоритмов выделения трендов
Зайков Геннадий Алексеевич
Численные алгоритмы расчета краевых задач теории упругости для составных систем специального вида
Измакова Ольга Анатольевна
Рекуррентные алгоритмы обучения и самообучения в теории распознавания образов

© Научная электронная библиотека «Веда», 2003-2013.
info@lib.ua-ru.net