Методы интерполирования
Начальные сведения об интерполяции. Обзор некоторых методов.
Добавлена: 2008-02-29 Просмотров:373 | Рейтинг:0.00

Решение задачи методами интерполирования предполагает выполнение условия

P(x_i) = f(x_i) при  x_i in [a,b]  (1)

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

P_n(x) = sum{i=0}{n}{a_i x^i}  (2)

Задача (1) имеет решение, если степень полинома n = m-1, где m - количество точек интервала [a,b], в которых задана функция f(x). Таким образом, многочлен n-ой степени может обеспечить совпадение с приближаемой функцией f(x) в ( n+1) точке конечного интервала. Однако, поведение разности P(x) - f(x) в точках x in [a,b] и  x<>x_i при построении полинома не оговаривается. При определённых ограничениях, накладываемых на f(x), и разумном выборе координат узлов сетки xi величина {max}under{x} delim{|}{P_n (x) - f(x)}{|} {right}under{n right infty} 0. Выполнением этого условия определяется свойство равномерной сходимости интерполяционного процесса.

Можно показать, что величина уклонения delim{|}{f(x)-P_n(x)}{|}, x in [a,b]  зависит, в частности, от расположения узлов интерполяции xi на интервале [a,b]. Близкой к оптимальной, в смысле минимума функции {max}under{x} delim{|}{f(x)-P_n(x)}{|} , при x in [a,b], является, так называемая, чебышевская сетка. Координаты её узлов определяются выражением

x_k = {a + b} / 2 + { b - a } / 2  cos {(2k + 1 ) pi } / { 2 ( n + 1 ) } ( 3 )

Эта сетка обладает рядом замечательных свойств. Так, если функция f(x) обладает непрерывной первой производной, то интерполяционный процесс на чебышевской сетке всегда сходится. Эта сетка также позволяет: 

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

Однако, следует иметь в виду, что достоинства этой сетки не предопределяют
абсолютную целесообразность её использования. Это связано с тем, что: 

    • она может оказаться не всегда удобной для вычисления соответствующих значений функции f(x) и
    • сходимость интерполяционного процесса может потребовать построения другой сетки. 

В общем случае, для каждой функции существует своя система сеток, на которых интерполяционный процесс сходится.

Рассмотрим различные алгоритмы построения многочлена P(x). Важно помнить, что какой бы алгоритм не был использован и, следовательно, в какой бы форме не был представлен многочлен степени n, при фиксированной сетке решение задачи (1 ) будет единственным.

Рейтинг@Mail.ru Rambler's Top100