Формулы приближенного дифференцирования, основанные на первой интерполяционной формуле Ньютона

Пусть функция y(x), задана таблично в равноотстоящих узлах x_i на отрезке [a,b]:

x	x₀	x₁	x₂	…	x_n
f(x)	y₀	y₁	y₂	…	y_n

Для нахождения на [a,b] производных , и т.д., функцию y приближенно заменим интерполяционным полиномом Ньютона, построенным для системы узлов :

где и .

Производя перемножение биномов, получим:

Будем дифференцировать данный многочлен как сложную функцию: . Получим:

(1)

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

(2)

Таким же способом можно вычислить производные любого порядка.

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

Иногда требуется находить производные функции в основных табличных точках x_i. В этом случае формулы численного дифференцирования упрощаются. Так как каждое табличное значение можно считать за начальное, то положим x=x₀, следовательно, t=0.

Получим:

		(3)
		(4)

Погрешность приближенного вычисления первой производной в нулевой точке будет равна

Пример 1. Приняв шаг h = 0,05, построить на отрезке [3,5; 3,6] интерполяционный полином Ньютона для функции у = е^х, заданной таблицей

х	3,50	3,55	3,60	3,65	3,70
у	33,115	34,813	36,598	38,475	40,447

Составляем таблицу разностей (Таблица 3).

Заметим, что в столбцах разностей, следуя обычной практике, мы не указываем десятичные разряды (которые ясны из столбца значений функции). Так как разности третьего порядка практически постоянны, то в формуле (15) полагаем n = 3. Приняв x ₀ = 3,50, y ₀= 33,115, будем иметь:

P ₃(x) = 33, 115+ 1,698 t + 0,087 + 0,005

или

P ₃(x) = 33, 115+ 1,698 t + 0,0435 + 0,00083 ,

где

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

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

Пример 2. Построить эмпирическую формулу для функции у, заданной таблично