2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки


Правила форума


В этом разделе нельзя создавать новые темы.



Начать новую тему Ответить на тему
 
 Аппроксимация сплайнами с заданными условиями
Сообщение06.07.2011, 15:27 
Аватара пользователя


25/08/09
19
Chemnitz, Germany
Добрый день, уважаемые участники форума!

Существует следующее множество точек:
Изображение
Мне необходимо аппроксимировать данное множество c помощью прямых линий и сплайнов таким образом, чтобы в точках сшивания кусков сохранялась непрерывность первой и второй производных.

По каждой из координат сплайны задаются параметрическим полиномом 5-й степени:
$x=C_0+C_1\cdot t+C_2\cdot t^2+C_3\cdot t^3+C_4\cdot t^4 +C_5\cdot t^5$
$y=D_0+D_1\cdot t+D_2\cdot t^2+D_3\cdot t^3+D_4\cdot t^4 +D_5\cdot t^5,$ где $t\in[0;1]$

Очевидно, что первые две точки соединяются горизонтальной прямой $y=0$

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

(Подробнее)

В случае зависимости $y=f(x)$ я решаю систему линейных уравнений таким образом:

$A=\begin{pmatrix} 1 & \sum\limits_{k=0}^K x_k & \sum\limits_{k=0}^K x^2_k & \sum\limits_{k=0}^K x^3_k & \sum\limits_{k=0}^K x^4_k & \sum\limits_{k=0}^K x^5_k \\ \sum\limits_{k=0}^K x_k & \sum\limits_{k=0}^K x^2_k & \sum\limits_{k=0}^K x_k^3 & \sum\limits_{k=0}^K x^4_k & \sum\limits_{k=0}^K x^5_k & \sum\limits_{k=0}^K x^6_k \\ 1 & x_{start} & x^2_{start} & x^3_{start} & x^4_{start} & x^5_{start} \\ 0 & 1 &  2x_{start} & 3x^2_{start} & 4x^3_{start} & 5x^4_{start} \\ 0 & 0 & 2 & 6x_{start} & 12x^2_{start} & 20x^3_{start} \\ 1 & x_{end} & x^2_{end} & x^3_{end} & x^4_{end} & x^5_{end}\end{pmatrix} B=\begin{pmatrix} \sum\limits_{k=0}^K y_k \\ \sum\limits_{k=0}^K y_kx_k \\ y_{start} \\ y'_{start} \\ y''_{start} \\ y_{end} \end{pmatrix}$

Итоговая матрица коэффициентов: $C=A^{-1}\cdot B$

Здесь 1-я и 2-я строки матриц A и B определяют аппроксимацию (полиномиальную регрессию) методом наименьших квадратов, 3-я строка определяет начальную точку полинома, 4-я $-$ значение первой производной в начальной точке, 5-я $-$ значение 2-й производной в начальной точке, 6-я $-$ конечную точку полинома.

Как будет выглядеть система в случае если $x=f_1(t),$ а $y=f_2(t)$ при $t\in[0;1]?$ Причем замечу, что мне важна непрерывность до второй производной не конкретной координаты $x$ или $y,$ а именно сплайнов, которые задаются через эти полиномы?

 Профиль  
                  
 
 Re: Аппроксимация сплайнами с заданными условиями
Сообщение07.07.2011, 20:22 
Заслуженный участник
Аватара пользователя


30/01/09
7068
У Вас трудность в том, что не понятно, какому значению параметра $t$ соответствует каждая аппроксимируемая точка. Если бы такое соответствие было, то задача раскладывалась бы в две одномерные. Выходов два. Первый - самому придумать какой-нибудь простой алгоритм параметризации конечной системы точек, исходя из специфики задачи. Второй - поискать, что сделано на эту тему в специальной литературе.

 Профиль  
                  
 
 Re: Аппроксимация сплайнами с заданными условиями
Сообщение08.07.2011, 07:14 
Заслуженный участник
Аватара пользователя


23/08/07
5495
Нов-ск
Destin в сообщении #465766 писал(а):
По каждой из координат сплайны задаются параметрическим полиномом 5-й степени
Задавайте полиномом третьей степени. По каждой из координат стройте сплайн независимо. В качестве параметра возьмите длину кусочно линейного пути от точки к точке.

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 3 ] 

Модераторы: Модераторы Математики, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group