Меню
Главная
Авторизация/Регистрация
 
Главная arrow Информатика arrow ИНФОРМАТИКА.
Посмотреть оригинал

Интерполяционные многочлены Ньютона

Рассмотрим понятие конечных разностей.

Пусть задана функция у=f{x) на отрезке [х0, х„], который разбит на п одинаковых отрезков (случай равноотстоящих значений аргумента): Ax=h = const. Для каждого узла х0, х, 0 + /г, ..., х„ () + п h определены значения функции в виде

Введем понятие конечных разностей.

Конечные разности первого порядка

Конечные разности второго порядка Аналогично определяются конечные разности высших порядков:

Конечные разности функций удобно располагать в таблицах, которые могут быть диагональными (табл. 5.1) или горизонтальными (табл. 5.2).

Диагональная таблица

Таблица 5.1

X,

У,

Лу/

Д2у,-

А3у/

<1

А"У/

Хо

Уо

Ауо

х.

У

А2у0

Ау,

Д3Уо

*2

У2

А‘у,

А4у0

Ду2

Д-yi

А5у0

*3

Уз

<1

А4у,

Ауз

А3у2

х4

У4

А2у3

Ау4

*5

_В_

Горизонтальная таблица

Таблица 5.2

х,-

Я

АУ/

А2у,-

А3у,

<1

а5у,

Хо

Уо

Ау0

Д2уо

А3уо

Д4у<>

А5Уо

Х|

У

Ау,

А2у,

А3у,

а4у.

х2

У2

Ау2

А2у2

Д3У2

Хз

Уз

Ау3

А2у3

Х4

У4

Ау4

Х5

_й_

Первая интерполяционная формула Ньютона

Пусть для функции у=/(х) заданы значения у, =/(х,) для равностоящих значений независимых переменных:

где h - шаг интерполяции.

Необходимо найти полином Р„{х) степени нс выше п, принимающий в точках (узлах) х, значения:

Интерполирующий полином ищется в виде:

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

Полагаем в (5.13) х=х0, т. к. второе, третье и другие слагаемые равны 0, то

Найдем коэффициент а{.

Приэс=Х1 получим:

Для определения а2 составим конечную разность второго порядка. При х=х2 получим:

Аналогично можно найти другие коэффициенты. Общая формула имеет вид:

Подставляя эти выражения в формулу (5.13), получаем:

где х„ ух - узлы интерполяции; х - текущая переменная; h - разность между двумя узлами интерполяции; h - величина постоянная, т. е. узлы интерполяции равно отстоят друг от друга.

Этот многочлен называют интерполяционным полиномом Ньютона для интерполяции в начале таблицы (интерполирование «вперед»), или первым полиномом Ньютона.

Для практического использования этот полином записывают в преобразованном виде, вводя обозначение t=(х - x0)/h, тогда

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

Блок-схема алгоритма метода Ньютона для интерполирования «вперед» приведена на рис. 5.3, программа - в приложении.

Пример 5.3. Дана таблица значений теплоемкости вещества в зависимости от температуры Cp=f{T) (табл. 5.3).

Таблица 5.3

х (7)

300

400

500

600

Y (Ср)

52,88

65,61

78,07

99,24

Построить интерполяционный многочлен Ньютона для заданных значений функции:

Составим таблицу конечных разностей функции (табл. 5.4).

Таблица 5.4

X

У

4)’

АТ

aV

300

52,88

12,73

-0,27

8,98

400

65,61

12,46

8,71

500

78,07

21,17

600

99,24

Воспользуемся формулой (5.16):

Блок-схема алгоритма метода Ньютона для интерполирования «вперед»

Рис. 5.3. Блок-схема алгоритма метода Ньютона для интерполирования «вперед»

После выполнения преобразований получим интерполяционный многочлен вида:

Полином имеет третью степень и дает возможность вычисления при помощи найденной формулы значения у для неизвестного х.

Пример 5.4. В табл. 5.3.1 приведены значения теплоемкости в зависимости от температуры. Определить значение теплоемкости в точке Г=450 К.

Воспользуемся первой интерполяционной формулой Ньютона. Конечные разности рассчитаны в предыдущем примере (табл. 5.3.2), запишем интерполяционный многочлен при х=450 К:

Таким образом, теплоемкость при температуре 450 К будет

Значение теплоемкости при Г=450 К получили такое же, что и рассчитанное по формуле Лагранжа.

Вторая интерполяционная формула Ньютона

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

Коэффициенты а0, аь ..., а„ определяем из условия:

Полагаем в (5.18) х=х„, тогда

Полагаем х=х„_|, тогда следовательно,

Если x=xn-2i то

Аналогично можно найти другие коэффициенты многочлена (5.18):

Подставляя эти выражения в формулу (5.18), получим вторую интерполяционную формулу Ньютона, или многочлен Ньютона для интерполирования «назад»:

Введем обозначения:

Произведя замену в (5.19), получим:

Это вторая формула Ньютона для интерполирования «назад».

Пример 5.5. Вычислить теплоемкость (см. табл. 5.3) для температуры Г=550 К.

Воспользуемся второй формулой Ньютона (5.19) и соответствующими конечными разностями (см. табл. 5.4):

Следовательно, значение теплоемкости при температуре 550 К равно

 
Посмотреть оригинал
Если Вы заметили ошибку в тексте выделите слово и нажмите Shift + Enter
< Предыдущая   СОДЕРЖАНИЕ   Следующая >
 

Популярные страницы