Купить Matlab  |  Mathematica  |  Mathcad  |  Maple  |  Statistica  |  Другие пакеты Поиск по сайту
Internet-класс  |  Примеры  |  Методики  |  Форум  |  Download
https://hub.exponenta.ru/


Курс введения в вычислительную математику.
Готовые занятия

Список курсов ВМ

 
Занятие 8
Теоретический материал Теоретическая справка Примеры Задачи для самостоятельного решения Контрольные вопросы

Постановка задачи приближения функции по методу интерполяции ~ Интерполяционный многочлен Лагранжа ~ Разделенные разности~ Интерполяционный многочлен Ньютона~ Оценка погрешности интерполяции

 

Постановка задачи интерполяции функций.

Пусть функция y = f(x) задана таблицей своих значений:

image002.gif(1118 bytes), i=0,1,...n. Требуется найти многочлен степени n, такой, что значения функции и многочлена в точках таблицы совпадают:

image004.gif(328 bytes), i=0,1,... n.

Справедлива теорема о существовании и единственности интерполяционного многочлена.

 

Одна из форм записи интерполяционного многочлена - многочлен Лагранжа:

image006.gif(437 bytes), где

image008.gif(1296 bytes)

Многочлен image010.gif(229 bytes)представляет собой многочлен степени n , удовлетворяющий условию

image012.gif(447 bytes) .

Таким образом, степень многочлена image014.gif(189 bytes) равна n и при image016.gif(206 bytes) в сумме обращаются в нуль все слагаемые, кроме слагаемого с номером image018.gif(200 bytes), равного image020.gif(181 bytes). Поэтому многочлен Лагранжа является интерполяционным.

В начало страницы

ПРИМЕР 1.Построение многочлена Лагранжа.

 

Другая форма записи интерполяционного многочлена - интерполяционный многочлен Ньютона с разделенными разностями. Пусть функция f задана с произвольным шагом и

точки таблицы занумерованы в произвольном порядке. Величины

image022.gif(534 bytes) называют разделенными разностями первого порядка. Разделенные разности второго порядка определяются формулой:

image024.gif(663 bytes).

Определение разделенной разности порядка image026.gif(210 bytes)таково:

image028.gif(720 bytes).

Используя разделенные разности, интерполяционный многочлен Ньютона можно записать в следующем виде:

0301.gif (1360 bytes)

0302.gif (1180 bytes)

В начало страницы

ПРИМЕР 2. Построение интерполяционного многочлена Ньютона с разделенными разностями.

Величину image032.gif(397 bytes)называют погрешностью интерполяции или остаточным членом интерполяции.

 

Оценка погрешности интерполяции.

Если функция n+1 раз на отрезке [a,b] , содержащем узлы интерполяции image034.gif(181 bytes), i=0,1,...n, то для погрешности интерполяции справедлива оценка:

image036.gif(722 bytes). Здесь image038.gif(520 bytes), image040.gif(372 bytes).

Эта оценка показывает, что для достаточно гладкой функции при фиксированной степени интерполяционного многочлена погрешность интерполяции стремится к нулю не медленнее, чем величина, пропорциональная image042.gif(225 bytes). Этот факт формулируют так: интерполяционный многочлен степени n аппроксимирует функцию с (n+1) порядком точности относительно image044.gif(211 bytes).

В начало страницы

ПРИМЕР 3. Использование остаточного члена интерполяции.

В практическом плане формула Ньютона обладает преимуществами перед формулой Лагранжа. Предположим, что в необходимо увеличить степень многочлена на единицу, добавив в таблицу еще один узел image046.gif(197 bytes). При использовании формулы Лагранжа это приводит к необходимости вычислять каждое слагаемое заново. Для вычисления image048.gif(239 bytes)достаточно добавить к image050.gif(225 bytes)лишь очередное слагаемое: image052.gif(496 bytes). Если функция f достаточно гладкая, то справедливо приближенное равенство image054.gif(433 bytes). Это равенство можно использовать для практической оценки погрешности интерполяции :image056.gif(376 bytes).

В начало страницы

Примеры Задачи для самостоятельного решения Контрольные вопросы
| На первую страницу | Поиск | Купить Matlab

Исправляем ошибки: Нашли опечатку? Выделите ее мышкой и нажмите Ctrl+Enter


Copyright © 1993-2024. Компания Softline. Все права защищены.

Дата последнего обновления информации на сайте: 04.03.17
Сайт начал работу 01.09.00