CC BY
39
Краткие сообщения
УДК 681.327.22:621.376.57
П.П. Кравченко, И.С. Павлова
МОДИФИЦИРОВАННЫЕ АЛГОРИТМЫ ИНТЕРПОЛЯЦИИ СПЛАЙНАМИ НА ОСНОВЕ ОПТИМИЗИРОВАННЫХ ДЕЛЬТА-ПРЕОБРАЗОВАНИЙ ВТОРОГО ПОРЯДКА
Рассматривается вопрос построения на основе дельта-преобразований второго порядка алгоритмов локальных сплайнов, характеризующихся низкой вычислительной трудоемкостью, и использование которых позволяет решать вопросы формирования требуемого характера изменения интерполируемой кривой между узлами интерполяции.
Оптимизированные дельта - преобразования второго порядка; сплайн; интерполяция.
P.P. Kravchenko, I.S. Pavlova
THE MODIFIED ALGORITHMS OF INTERPOLATION BY SPLINES ON THE BASIS OF OPTIMIZED DELTA-TRANSFORMATIONS OF THE SECOND
ORDER
In the article the question of construction algorithms of the local splines based on the delta -transformations of the second order, characterized by low computational complexity, and the use of which enables us to solve problems of forming the required character of the variation of the interpolated curve between the points of interpolation is considered.
Optimized delta - transformations of the second order; a spline; an interpolation.
При построении графиков, двух - и трехмерных объектов, рельефов поверхностей в реальном масштабе времени используются интерполяционные и аппроксимирующие сплайны. На практике особый интерес представляет использование интерполяционных сплайнов, характеризующихся низкой вычислительной трудоемкостью и возможностями адаптации к особенностям интерполируемой траектории.
Среди известных сплайнов с низкой вычислительной трудоемкостью можно выделить кубические [1] и на основе оптимизированных дельта-преобразований второго порядка [2]. В частности, использование алгоритмов на основе дельтапреобразований второго порядка позволяет обеспечить предельно низкую вычислительную трудоёмкость (3 базовых операции сложения) на этапе построения собственно сплайна.
В данной статье рассматривается модифицированный алгоритм построения сплайнов с использованием дельта-преобразований второго порядка, предоставляющий возможности изменения характера интерполируемой кривой, а также уменьшения исходного количества сплайнов.
Сущность алгоритма заключается в формировании сплайна по трем отсчетам, один из которых, расположенный между двумя крайними узлами, вводится как дополнительный. Итоговая интерполяционная кривая проходит через заданные узлы, а так же через вспомогательный отсчет или вблизи от данного отсчета.
Краткие сообщения
Рассматриваемый алгоритм включает в себя два основных этапа: вычисление параметров сплайна (подготовительный этап) и собственно построение сплайно-вой кривой.
Введем обозначение пар координат узловых точек в виде ((t„ , Уп , yn X (tc„ , ycn ), On+1, Уп+1 )) , где n - номер узла иHтеPп0ляЦии, tn - независимая переменная, yn - отсчет в n-ом узле, ycn - дополнительный отсчет, y„ -производная в n-ом узле.
Параметры сплайна (квант преобразования Pn и расстояние от узла (n + 1)
до точки переключения знака кванта преобразования Dn ) на подготовительном этапе вычисляются с использованием приводимых ниже формул:
-L - sign (L )JL2 - 4G V
p __ „ о у „ J у „ „ „ •
0,5T„2
D_
( Уп+1 - ycn ) / hn + hA I2 - yn + T„P„
2Pn
где
Ln = Уп+1 - Уп - 0, 5Tn (Уп + (Уп+1 - Уеп )/ hn ) - 0-25hn ((Уп+1 - Усп )/ hn - Уп У;
Gn = (Tn/2 - hn /4)2 - 2(hn /4)2; Fn = -0.25(( Уп+1 - ycn)/hn - Уп )2;
hn = {n+1 - ten ; 0 < hn < 0 8Tn ; Tn = tn+1 - tn .
Для (n+1) - го сплайна производная yn+1 определяется в виде:
yn+1 = (Уп+1 - Усп)/ К + КР„ /2.
Пусть vt - шаг дискретизации независимой переменной, NA - количество дискретных шагов в диапазоне [t ; t„+1 - Dn ], NB - количество дискретных шагов
в диапазоне [t„; t„+1 - Dn ].
Пусть на подготовительном этапе помимо отмеченных выше значений Pn и Dn определены также c„ = P„Vt2, VSN = y„+1V t, VS0 = y„Vt [2]. Подготовительный этап может быть в вычислительном процессе совмещен во времени с собственно построением сплайновой кривой или предшествовать построению сплайновой кривой с сохранением определяемых параметром.
Из узла n+1 строится кривая при движении влево с шагом V t на расстояние D„ до точки переключения кванта преобразования (NA шагов). Расчет выполняется по формулам:
YNa = yn+1 ;
VSN , = VSN - c*;
Na-i- 1 Na-i „ ’
где i = 0,... Na 1, NЛ _
Yna -i-1 _ Yn> -i -VSna -i-1 - 0,5c* D„
Vt
Далее из узла п строится кривая при движении вправо с шагом V; на расстояние (Тп - Вп) до точки переключения кванта преобразования (МБ шагов). Расчет выполняется по формулам:
^0 = Уп ;
ТО,+1 = ТО, + сП;
F _ Y + vs.+, - 0,5c*.
і+1 і і+1 ? „ ■
где 1 = ^Nb -1, nb _
Tn - Dn
nn
Vt
При использовании алгоритмов с целью формирования определенной направленности сплайна, дополнительный отсчет вводится на основе реальных данных об интерполируемом участке (например, известного аналитического представления интерполируемой функции или нечетких требований к формированию участка кривой и т. п.).
Для обеспечения достаточного качества сплайновой интерполяции целесообразно наряду с рассмотренными алгоритмами использовать также алгоритмы без дополнительного отсчета [2].
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Квасов Б.И. Методы изогеометрической аппроксимации сплайнами. - М.: Ижевск: НИЦ «Регулярная и хаотическая динамика», Институт компьютерных исследований, 2006. - 416 с.
2. Кравченко П.П. Основы теории оптимизированных дельта - преобразований второго порядка. Цифровое управление, сжатие, параллельная обработка информации: Монография. - Таганрог: Изд-во ТТИ ЮФУ, 2008. - 192 с.
Кравченко Павел Павлович
Технологический институт федерального государственного автономного образовательного учреждения высшего профессионального образования «Южный федеральный университет» в г. Таганроге.
Е-mail: kravch@tsure.ru.
347928, .Таганрог, пер. Некрасовский, 44.
Тел.: 88634314945.
Павлова Ирина Сергеевна
Е-mail: irina_p.84@mail.ru.
Kravchenko Pavel Pavlovich
Taganrog Institute of Technology - Federal State-Owned Educational Establishment of Higher Vocational Education “Southern Federal University”.
Е-mail: kravch@tsure.ru.
44, Nekrasovski, Taganrog, 347928, Russia.
Phone: +7(8634314945.
Pavlova Irina Sergeevna Е-mail: irina_p.84@mail.ru.
