Реферат: Некоторые дополнительные вычислительные методы

2

0,5

0,503

0,503

-0,688

-0,686

------

0,503

0,504

------

------

0,0015

0,0005

3. Интерполирование и экстраполирование

Задача интерполирования состоит в том, чтобы по значениям функции f(x) в нескольких точках отрезка восстановить ее значения в остальных точках данного отрезка. Разумеется, такая постановка задачи допускает сколь угодно много решений. Задача интерполирования возникает, например, в том случае, когда известны результаты измерений y_k = f(x_k ) некоторой физической величины f(x) в точках x_k , k = 0, 1,…, n и требуется определить ее значение в других точках. Интерполирование используется также при необходимости сгущения таблиц, когда вычисление значений f(x) по точным формулам трудоемко. Иногда возникает необходимость приближенной замены (аппроксимации) данной функции (обычно заданной таблицей) другими функциями, которые легче вычислить. При обработке эмпирических (экспериментальных) зависимостей, результаты обычно представлены в табличном или графическом виде. Задача заключается в аналитическом представлении искомой функциональной зависимости, то есть в подборе формулы, корректно описывающей экспериментальные данные.

Интерполирование с помощью многочленов

Пусть функциональная зависимость задана таблицей y_{0 =} f(x₀ ); …, y₁ = f(x₁ ); …, y_{n =} f(x_n ). Обычно задача интерполирования формулируется так: найти многочлен P(x) = P_n (x) степени не выше n, значения которого в точках x_i (i = 0, 1 2,…, n) совпадают со значениями данной функции, то есть P(x_i ) = y_i . Геометрически это означает, что нужно найти алгебраическую кривую вида проходящую через заданную систему точек М_i (x_i , y_i ) (см. рис. 4). Многочлен Р(х) называется интерполяционным многочленом. Точки x_i (i = 0, 1, 2,…, n) называются узлами интерполяции.

Для любой непрерывной функции f(x) сформулированная задача имеет единственное решение. Действительно, для отыскания коэффициентов а₀ , а₁ , а₂ ,…, а_n получаем систему линейных уравнений определитель которой отличен от нуля, если среди точек x_i (i = 0, 1, 2,…, n) нет совпадающих. Решение системы можно записать различным образом. Однако наиболее употребительна запись интерполяционного многочлена в форме Лагранжа или в форме Ньютона.

Инт ерполяционный многочлен Лагранжа

Пусть на отрезке [a,b ] некоторая функция f(x ) задана лишь в некоторых точках , т.е. известны ее значения , которые, собирают в таблицу:

x	x0	x1	...	xn
f(x)	y0	y1	...	yn

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

Этот многочлен называется многочленом Лагранжа .

Его основные свойства:

1) это - многочлен степени ;

2) , т.е. многочлен Лагранжа имеет в точках те же значения, что и функция ;

3) если фиксировать любое число то окажется выполненным неравенство

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

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

Пример. Построить интерполяционный многочлен Лагранжа для функции заданной таблицей

x	1	2	3	5
y	1	5	14	81

И найти значение функции при x=4.

Решение. Используя формулу Лагранжа найдем:

После некоторых преобразований получим Тогда f(4)≈L₃ (4)=36,5.

Интерполяционные многочлены Стирлинга и Бесселя