- 64-

3. Приближение функций.

Пусть на отрезке определена некоторая функцияоднако полная информация о ней недоступна. Известны лишь ее значения в конечном числе точекэтого отрезка, которые мы будем считать занумерованными в порядке возрастания:

.

Требуется по известным значениям

,

«восстановить», хотя бы приближенно, исходную функцию то есть построить на отрезкефункцию, достаточно близкую к. Функцию принято называть интерполирующей функцией, точки - узлами интерполяции.

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

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

---

1. Интерполирование.

   1. Классическая постановка задачи интерполирования.

Выберем некоторую систему функций , заданных на отрезке, и будем строитькак их линейную комбинацию:

,

где числовые коэффициенты подлежат определению, согласно условиям:

, .

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

,

или в развернутом виде:

.

Для того, чтобы коэффициенты можно было определить и притом единственным образом, необходимо и достаточно, чтобы определитель полученной системы линейных уравнений был отличен от нуля:

.

Определение.

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

Очевидно, что для однозначной разрешимости задачи интерполирования в классической постановке необходимо и достаточно, чтобы система функций была Чебышевской. Только такие системы функций мы и будем использовать в этой главе. Необходимым условием принадлежности системы функцийк Чебышевской является их линейная независимость. Однако это условие не является достаточным. Например, для системы из двух линейно независимых функций, с узлами интерполяцииопределитель

---

и данная система функций при выбранных значениях ине является Чебышевской.

2. Интерполирование полиномами.

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

.

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

с неизвестными коэффициентами .

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

, .

Определителем этой системы является определитель Ван-дер- Монда:

.

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

Задача 1.

Построить линейный полином

по заданным узлам интерполяции и соответствующим им значениям функции

и.

---

Линейная система уравнений для определения ив данном случае имеет вид:

,

.

Определитель этой системы равен . Решив систему, получим:

.

Следовательно,

.

Перепишем этот полином в несколько другой форме, выделяя ив качестве множителей

.

Геометрический образ интерполирующей функции - прямая, проходящая на плоскостичерез точки с координатамии. Уравнение этой прямой, наряду с и , можно переписать в виде:

.

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

---

3. Построение интерполяционного полинома в форме Лагранжа.

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

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

,

где полиномы степени, «ориентированные» на точкив том смысле, что

Такие полиномы легко построить:

или в развернутом виде:

Иногда нам будет удобно записывать в виде:

.

Из выражения и формул очевидно, что построенный полином действительно является интерполяционным полиномом для функциина сетке с узлами. Его принято называть интерполяционным полиномом в форме Лагранжа. Этим подчеркивается, что возможны и другие эквивалентные представления интерполяционного полинома. С одним из них мы познакомимся в следующем разделе.

---

Смотрите также файлы

• BOOK_P~2.doc

• BOOK_P~1.doc

• BOA462~1.doc

• 8.doc

• 7.doc