Интерполяция и экстраполяция сигналов [39, 55].

Полиномиальная интерполяция. Одним из основных типов точечной аппроксимации является интерполирование. Сущность его заключается в том, что функции y_i = s(x_i) сопоставляется интерполяционный многочлен

f(x) = а₀ + а₁х + а₂х² + … + a_nxⁿ = a_i·xⁱ, (14.2.1)

принимающий в точках x_i те же значения y_i, что и функция f(x). Для выполнения полиномиальной интерполяции достаточно по выражению (14.2.1) составить систему линейных уравнений для n узловых точек и определить n значений коэффициентов a_i.

При N точках функции y_i максимальная степень интерполяционного многочлена n=N-1, и в этом случае говорят о глобальной интерполяции с прохождением f(x) через все значения точек y_i. Однако в этом случае при большом количестве узлов получается очень высокая степень многочлена. Кроме того, экспериментальные табличные данные могут содержать ошибки измерений, а глобальная интерполяция повторит все допущенные при измерениях ошибки. Для исключения этого фактора стараются подобрать аппроксимирующий многочлен как можно меньшей степени (как правило, n=1, 2, 3), график которого проходит близко от узловых точек. На практике мерой отклонения многочлена f(x) от заданной функции на множестве точек (x_i,y_i) является величина s среднеквадратичного приближения s² = S_i(f(x_i)-y_i)², минимальное значение которой обеспечивается подбором коэффициентов a_i.

Пример выполнения глобальной интерполяции приведен на рис.14.2.1. Равномерной дискретизации данных для интерполяции не требуется. Максимальная степень полинома на практике обычно устанавливается не более 8-10, а большие массивы данных интерполируются последовательными локальными частями.

Рис. 14.2.1. Интерполяция данных.

Линейная и квадратичная интерполяция являются самыми простыми способами обработки таблиц и выполняются по уравнениям:

f(x)_лин = а₀ + а₁х. f(x)_кв = а₀+а₁х+а₂х².

При линейной интерполяции вычисления дополнительных точек выполняются по линейной зависимости. Графически это означает простое соединение узловых точек отрезками прямых. В системе Mathcad для этого используется функция linterp(X,Y,x), где X и Y – вектора узловых точек. Функция linterp(X,Y,x) возвращает значение функции при её линейной интерполяции по заданным аргументам х. При небольшом числе узловых точек (менее 10) линейная интерполяция оказывается довольно грубой. Для целей экстраполяции функция linterp не предназначена.

Для практического использования более удобны формулы аппроксимации, не требующие предварительного определения коэффициентов аппроксимирующих полиномов. К числу таких формул относится интерполяционных многочлен по Лагранжу /40/. При аппроксимации функции у(х) многочленом n-й степени Y(x):

Y(x) = + +…

…+ . (14.2.2)

Пример интерполяции по Лагранжу приведен на рис. 14.2.2.

Рис. 14.2.2. Интерполяция по Лагранжу.

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

Метод построения кривых был предложен де Кастелье в 1959 году и основан на разбиении отрезков, соединяющих исходные точки в отношении t, а затем в рекурсивном повторении этого процесса для полученных отрезков:

(t) = (1-t) (t) + t (t), (14.2.3)

где нижний индекс - номер точки, верхний индекс - уровень разбиения. Уравнение кривой n-ого порядка задается уравнением: P(t) = (t).

Рис. 14.2.3.

Для примера построим кривую для трех опорных точек P₀(t=0), P₁ и P₂(t=1) на интервале t=[0, 1] (рис. 14.2.3).

Для каждого tÎ[0, 1] определим точку :

= (1-t)P₀ + tP₁, = (1-t)P₁ + tP₂,

= (1-t) + t =(1-t)²P₀(t)+2t(1-t)P₁(t)+t²P₂(t),

и тем самым получим кривую второго порядка.

Рис. 14.2.3.

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

=(1-t)P₀+tP₁, =(1-t)P₁+tP₂, =(1-t)P₂+tP₃,

=(1-t) +t =(1-t)²P₀(t)+2t(1-t)P₁(t)+t²P₂(t), =(1-t) +t =(1-t)²P₁(t)+2t(1-t)P₂(t)+t²P₃(t),

= (1-t) + t =(1-t)³P₀(t)+3t(1-t)²P₁(t)+ 3t²(1-t)P₂(t)+t³P₃(t).

Общая аналитическая запись для кривых Безье по N+1 опорной точке:

P^N(t) = P_i , = tⁱ (1-t)^N-i,

= N!/(i! (N-i)!) – биномиальные коэффициенты.

Кривые Безье всегда проходят через начальную P₀ и конечную P_N точки. Если рассматривать опорные точки в противоположном порядке, то форма кривой не изменяется. Если опорные точки лежат на одной прямой, то кривая Безье вырождается в отрезок, соединяющий эти точки. Степень многочлена, представляющего кривую в аналитическом виде, на 1 меньше числа точек.

Дата добавления: 2015-07-19