Компьютерная графика / МАШ_ГРАФИКА

1.1 Построение сплайнов по однократным узлам

Допустим, плоская кривая у(х) задана набором значений у(x_i) = у_i в узлах x_i (i = 0,1 ,...,n). Необходимо построить интерполяционный сплайн наименьшей степени, прохо-дящий через узлы x_i : S(x_i) = у_i, (i = 0,1 ,...,n). Если кривая является пространственной, то наряду с функцией у(х) должна быть задана зависимость z(x). Задача интерпо-лирования по оси z решается точно также , как и по оси у.

Рассмотрим отрезок [x_i,x_i₊₁]. В его граничных точках заданы по одному геометрическому условию: у(x_i) = у_i , у(x_i₊₁) = у_i₊₁ . Поскольку узлы являются однократными, то наименьшая степень полинома S_i₊₁(x) (i = 0,1 ,...,n-1) равна единице. Как показано в Главе 2 (2.11), интерполяцион-ный многочлен Лагранжа для S_i₊₁(x) имеет вид:

S_i+1(x)= у_i (x_i+1–x)/(x_i+1 -x_i)+ у_i+1 (x-x_i) /(x_i+1 -x_i). (3.1 а)

В векторной форме:

S_i+1(x) = (Ф(x),Y);

Ф(x) = ( (x_i+1 –x) / (x_i+1 -x_i), (x-x_i) / (x_i+1 -x_i));Y= ( у_i , у_i+1 ).

Переходя к безразмерной переменной t=(x-x_i) / (x_i₊₁ -x_i), у которой t(x_i)=0; t(x_i₊₁)=1, для функций Лагранжа и полинома S_i₊₁(t) получим следующие выражения:

Одномерный сплайн рассмотренного вида является лома-ной, приближённо заменяющей исходную кривую у(х).

Интерполирование двухмерными сплайнами поверх-ностей по однократным узлам производится аналогично. Допустим, поверхность z(x,у), которая задана набором зна-чений z(x_i,у_j) = z_i_j в узлах двухмерной сетки (x_i,у_j) (i = 0,1,...,n; j=0,1,...,m). Необходимо построить интерполяци-онный сплайн S(x,y) наименьшей степени, проходящий че-рез узлы (x_i,у_j) : S(x_i,у_j) = z_i_j (i = 0,1 ,..., n; j = 0,1 ,..., m).

Рассмотрим прямоугольный участок сетки [x_i,x_i₊₁] [у_j, y_j₊₁]. В его граничных точках заданы по одному геомет-рическому условию: z(x_i ,у_j )= z_ij , z(x_i₊₁,у_j ) =z₍_i₊₁₎_j , z(x_i , у_j₊₁)= z_i₍_j₊₁₎, z(x_i₊₁,у_j₊₁) =z₍_i₊₁₎₍_j₊₁₎. Минимальные степени квадратичной формы S₍_i₊₁₎,₍_j₊₁₎(x,y) по переменным х,у рав-ны единице. Как показано в Главе 2 (2.22), явное решение для S₍_i₊₁₎,₍_j₊₁₎(x,y) при помощи функций Лагранжа имеет вид:

S₍_I_+1),(_j₊₁₎ (x,y) =Ф(x) ZФ(у); (3.2 а)

Ф(x)=((x_i₊₁ –x) / (x_i₊₁ -x_i), (x-x_i) / (x_i₊₁ -x_i));

Ф(у)=((y_j₊₁ –y) / (y_j₊₁ -y_j), (y-y_j) / (y_j₊₁ -y_j));

Переходя, как и в одномерном случае, к безразмерным переменным t_х= (x - x_i)/(x_i₊₁ - x_i) и t_у = (у - y_j)/(y_j₊₁- y_j) и

используя матрицу М_Л , решение для S₍_i₊₁₎,₍_j₊₁₎(x,y) можно представить в виде:

Построенный двухмерный сплайн S(x,y) = { S₍_i₊₁₎,₍_j₊₁₎(x,y), (i = 0,1 ,..., n-1; j = 0,1 ,..., m-1) } задаёт множество про-странственных четырёхугольников (в общем случае –не-плоских), заменяющих исходную поверхность z(x,у) таким образом, что они совпадают с заданными значениями в уз-лах интерполирования (x_i ,у_j).

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

Содержание