Processing math: 100%

12.7 Формула Тейлора

В одномерном случае формула Тейлора с остатком в форме Лагранжа содержится в следующей теореме.

Теорема. Пусть функция γ на отрезке [α,β] имеет непрерывные производные до порядка q включительно, а на интервале (α,β) существует производная порядка q+1. Тогда справедливо равенство γ(β)γ(α)=
γ(α)1!(βα)+γ(α)2!(βα)2++γ(q)(α)q!(βα)q+γ(q+1)(ξ)q+1!(βα)q+1,
где ξ — некоторая точка из интервала (α;β).

Аналог этой теоремы в многомерном случае может иметь следующий вид.

Теорема. Пусть действительная функция f класса Cq+1 на открытом множестве ERn и пусть отрезок [a,a+h]E. Тогда справедливо равенство: f(a+h)f(a)=

=ni=1fxi(a)hi+12!ni,j=12fxixj(a)hihj++

+1q!ni1,,iq=1qfxi1xiq(a)hi1hiq+Rq,

где Rq=1q+1!ni1,,iq+1=1q+1fxi1xiq+1(a+θh)hi1hiq+1,

а θ — некоторое число из отрезка [0,1].

Положим γ(t)=f(a+th)(0t1). Ранее была доказана лемма(12.4 стр.283), согласно которой функция γ дифференцируема и её производная

γ(t)=ni=1fxi(a+th)hi,(0t1).

Снова применяя эту лемму получим

γ(t)=ni,j=12fxixj(a+th)hihj.

По индукции получаем

γp(t)=ni,,ip=1pfxi1xip(a+th)hi1hip,(0pq+1)

Применяя теперь формулу Тейлора для функции γ, находим
γ(1)γ(0)=γ(0)+12!γ(0)++1q!γ(q)(0)+1(q+1)!γ(q+1)(θ),
где 0θ1 Если воспользуемся найденными выражениями для производных функции γ и учтем, что γ(1)γ(0)=f(a+h)f(a), то получим требуемое равенство.

Примеры решения задач

  1. Найти разложение функции по формуле Тейлора второго порядка в окрестности точки M0(2,1). f(x,y)=x22x3y, используя найденное разложение, найти приближенное значение функции в точке M(2,05;0,98).

    Решение

    Вычислим все необходимое для решения:
    f(2;1)=22223=2;
    fx=2x2x3y+x22x3yln2;fy=3x22x3yln2
    fx(2;1)=2223(2+2ln2)=2(1+ln2); fy(2;1)=34223ln2=61ln2
    2fx2=2x3y(2+4xln2+x2ln22);2fxy=3x2x3y(2+xln2)ln2;2fy2=9x22x3yln22;
    2fx2(2;1)=1+4ln2+2ln22;2fxy(2;1)=6(1+ln2)ln2;2fy2(2;1)=18ln22;

    Запишем формулу Тейлора второго порядка для заданной функции:f(x;y)2+(2(1+ln2)(x2)6ln2(y1))+ 12[(1+4ln2+2ln22)(x2)212(1+ln2)ln2(x2)(y1)+18ln22(y1)2]

    Найдём приближенное значение f(2,05;0,98)2,087.

Формула Тейлора для действительных функций

Пройдите этот тест, чтобы проверить свои знания по прочитанной теме.

Список использованной литературы

О приближенном вычислении с помощью формулы Тейлора

Если остаток в формуле Тейлора |r_{n}(x_{0},x)|< \alpha _{0} &s=1, то формулу Тейлора для многочлена можно записать так:  f(x)\approx f(x_{0})+\frac{f'(x_{0})}{1!}(x-x_{0})+\frac{f»(x_{0})}{2!}(x-x_{0})^{2}+…+\frac{f^{(n)}(x_{0})}{n!}(x-x_{0})^{n} &s=1.

Важна форма записи остаточного члена:

r_{n}(x_{0},x)=\frac{f^{(n+1)}(\xi )}{(n+1)!}(x-x_{0})^{n+1} &s=1.

r_{n}(x_{0},x) &s=1 — определяет погрешность формулы. Если же f(x) &s=1 вычисляется по формуле при конкретном числовом значении x &s=1, то может оказаться, что слагаемые в этой формуле сами вычисляются приближённо. Тогда погрешность результата будет состоять из погрешности слагаемых и погрешности формулы. Если вычислять все слагаемые с одинаковой точностью \alpha _{0} &s=1 (погрешностью формулы), то общая погрешность результата равна (n+2)\alpha _{0} &s=1.

Пусть \alpha &s=1 — заранее известная точность результата. Тогда следует преобразовать \alpha _{0} &s=1 так, чтобы обеспечить выполнение неравенства   (n+2)\alpha _{0}\leq\alpha &s=1, то есть \alpha_{0}\leq\frac{\alpha}{n+2} &s=1. При достаточно малых n &s=1, например, n\leq8 &s=1: \alpha_{0}=\frac{\alpha}{10}\leq\frac{\alpha}{n+2} &s=1.

Обычно точность вычислений \alpha &s=1 задается в виде: \alpha=10^{-m} \Rightarrow \alpha_{0}=10^{-(m+1)} &s=1. Это значит, что вычисления нужно проводить с одним запасным знаком. Мы установили, что один запасной знак обеспечит требуемую точность при n\leq8 &s=1.

Пример

Вычислить e^{0,1} &s=1 с точностью до \alpha=0,001=10^{-3} &s=1.

Решение

Оценкой определим, в какой точке удобнее раскладывать исходную функцию (найдём ближайшую к необходимой точку, где известно точное значение функции):

0\leq0,1\leq0,5 \Rightarrow x\in[0;0,5] &s=1

Выпишем формулу Тейлора:

e^{x}=1+x+\frac{1}{2}x^{2}+…+\frac{x^{n}}{n!}+\frac{e^{\xi }}{(n+1)!}x^{n+1} &s=1;

Выполним вычисление по формуле Тейлора, разложив функцию в точке x_{0}=0 &s=1

Выполним оценку погрешности:

r_{n}(0,x)=\left | \frac{e^{\xi }x^{n+1}}{(n+1)!} \right |=\frac{e^{\xi} \left | x \right |^{n+1}}{(n+1)!}\leq\frac{\sqrt{e}x^{n+1}}{(n+1)!} \leq \frac{2x^{n+1}}{(n+1)!} &s=2

Оценим сверху:

\frac{2x^{n+1}}{(n+1)!}\leq \frac{1}{10} &s=2

Перенесём 2 в правую часть и выполним обозначение:

\frac{x^{n+1}}{(n+1)!}\leq 0.5*10^{-1}\alpha=\frac{\alpha}{20} &s=2.

Эта запись удобна тем, что вычисляя последовательность слагаемых U_{k}=\frac{x^{k}}{k!} &s=1 мы имеем возможность одновременно видеть достигнута ли требуемая точность.

По условию:

\alpha=10^{-3} &s=1

Подставим в оценку, сделанную ранее:

\frac{x^{n+1}}{(n+1)!}\leq 0,00005 &s=2

Для U_{k}=\frac{x^{k}}{k!} &s=1 полагаем k=0,1,2,… &s=1

x=0,1U0=1;U1=0,1;U2=0,005;

U_{3}=0,0002; U_{4}=0,00005 &s=1 — выбранное значение k подходит.

e^{0,1}\approx 1+0,1+0,005+0,0002+0,00005=1,105 &s=1

e^{0,1}\approx1,105 &s=1

Неравенство \frac{x^{n+1}}{(n+1)!}\leq0,00005 &s=2 оказалось выполненным при k=n+1=4 &s=1, n=3 &s=1.

Источники: