Приближение функции многочленами Тейлора

Пусть функция y = f(x) определена в окрестности точки a и имеет в этой окрестности n + 1 производную. Тогда в этой окрестности справедлива формула Тейлора:

f(x) = c0 + c1(x - a) + c2(x - a)2 + … + cn(x - a )n + Rn(x) = Tn(x) + Rn(x),

где

ck =

Tn(x) - многочлен Тейлора:

Tn(x)= c0 + c1(x - a) + c2(x - a)2 + … + cn(x - a )n, (4.1)

Rn(x) - остаточный член формулы Тейлора. Его можно записать различными способами, например, в форме Лагранжа:

Rn(x)= , a x.

Многочлен Тейлора (4.1) обладает свойством, что в точке x = a все его производные до порядка n включительно совпадают с соответствующими производными функции f, т. е.

T(a)= f(k)(a), k = 0, 1, …, n.

В этом легко убедиться, дифференцируя Tn(x). Благодаря этому свойству многочлен Тейлора хорошо приближает функцию f в окрестности точки a. Погрешность приближения составляет

|f(x) - Tn(x)| = |Rn(x)|,

т. е. задавая некоторую точность > 0, можно определить окрестность точки a и значение n из условия:

|Rn(x)| = < . (4.2)

Пример 4.1.

Найдем приближение функции y = sinx многочленом Тейлора в окрестности точки a = 0. Воспользуемся известным выражением для k-ой производной функции sinx:

(sinx)(k) = sin x + k (4.3)

Применяя последовательно формулу (4.3), получим:

f(0) = sin0 = 0;

f (0) = cos(0) = 1;

f"(0) = -sin0 = 0;

f(2k-1)(0) = sin (2k - 1) = (-1)k - 1 ;

f(2k)(0) = 0;

f(2k+1)() = (-1)kcos.

Следовательно, многочлен Тейлора для функции y = sinx для n = 2k имеет вид:

sinx = x - + … + (-1)k - 1 + R2k(x),

R2k(x) = (-1)k .

Зададим = 10 -4 и отрезок [-,]. Определим n =2k из неравенства:

|R2k(x)| = < < < = 10-4.

Таким образом, на отрезке -, функция y = sinx с точностью до = 10-4 равна многочлену 5-ой степени:

sinx = x - + = x - 0.1667x3 + 0.0083x5.

Пример 4.2.

Найдем приближение функции y = ex многочленом Тейлора на отрезке [0, 1] с точностью = 10 -5.

Выберем a = ?, т. е в середине отрезка. При этом величина погрешности в левой части (4.2) принимает минимальное значение. Из математического анализа известно, что для k-ой производной от ex справедливо равенство:

(ex)(k) = ex.

Поэтому

(ea)(k) = ea = e1/2,

Следовательно, многочлен Тейлора для функции y = ex имеет вид:

ex = e1/2 + e1/2(x - ?) + (x - ?)2 + … + (x - ?)n+ Rn(x),

При этом, учитывая, что x [0, 1], получим оценку погрешности:

|Rn(x)| < . (4.4)

Составим таблицу погрешностей, вычисленных по формуле (4.4):

 
n  
Rn 0.057 0.0071 0.00071 0.000059 0.0000043  
             

Таким образом, следует взять n = 6.