Теорема Тейлора

Материал из testwiki
Версия от 12:09, 3 марта 2022; imported>JamTypeX (Аналитичность)
(разн.) ← Предыдущая версия | Текущая версия (разн.) | Следующая версия → (разн.)
Перейти к навигации Перейти к поиску
Экспоненциальная функция y = ex (сплошная красная линия) и соответствующий многочлен Тейлора четвёртого порядка (штрих-пунктирная зелёная линия) вблизи начала координат
Эта статья о многочленах Тейлора дифференцируемых функций. О рядах Тейлора аналитических функций см. соответствующую статью.

Теорема Тейлора даёт приближение к функции, дифференцируемой k раз, вблизи данной точки с помощью многочлена Тейлора k-го порядка. Для аналитических функций многочлен Тейлора в данной точке является частичной суммой их ряда Тейлора, который, в свою очередь, полностью определяет функцию в некоторой окрестности точки. Точное содержание теоремы Тейлора до настоящего времени не согласовано. Конечно, существует несколько версий теоремы, применимых в различных ситуациях, и некоторые из этих версий содержат оценки ошибки, возникающей при приближении функции с помощью многочлена Тейлора.

Эта теорема названа в честь математика Брука Тейлора, который сформулировал одну из её версий в 1712 году. Явное выражение для ошибки приближения было дано намного позже Жозефом Лагранжем. Ранее, в 1671 году, Джеймсом Грегори уже было упомянуто следствие из теоремы.

Теорема Тейлора позволяет овладеть приёмами вычислений начального уровня, и она является одним из центральных элементарных инструментов в математическом анализе. При изучении математики она является начальной точкой для изучения асимптотического анализа. Теорема также используется в математической физике. Она также обобщается на функции нескольких переменных и векторные функции f:nm для любых размерностей n и m. Это обобщение теоремы Тейлора является базовым для определения так называемых струй, которые появляются в дифференциальной геометрии и в теории дифференциальных уравнений с частными производными.

Предпосылки для введения теоремы

График Шаблон:Nowrap (голубого цвета) с его линейным приближением Шаблон:Nowrap (красным цветом) в точке a = 0.

Если вещественно-значимая функция f(х) является дифференцируемой в точке a, то она имеет линейное приближение в точке a. Это означает, что существует функция h1 такая, что

f(x)=f(a)+f(a)(xa)+h1(x)(xa),limxah1(x)=0.

Здесь

P1(x)=f(a)+f(a)(xa) 

это линейное приближение функции f в точке a. График функции Шаблон:Nowrap является касательной к графику функции f в точке Шаблон:Nowrap. Ошибка приближения такова

R1(x)=f(x)P1(x)=h1(x)(xa). 

Заметим, что ошибка приближается к нулю немного быстрее, чем разница Шаблон:Nowrap приближается к нулю по мере того, как x стремится к a.

График Шаблон:Nowrap (голубого цвета) с квадратичным приближением Шаблон:Nowrap (красного цвета) в точке a = 0. Заметны значительные улучшения приближения.

Если мы ищем лучшее приближение f, мы можем использовать многочлен второй степени вместо линейной функции. Вместо нахождения производной от f в точке a, мы можем найти две производные, получив таким образом многочлен, который так же как и f возрастает (или убывает), и так же как и f имеет выпуклость (или вогнутость) в точке a. Многочлен второй степени (квадратный многочлен) в этом случае будет выглядеть следующим образом:

P2(x)=f(a)+f(a)(xa)+f(a)2(xa)2.

Теорема Тейлора позволяет убедиться, что квадратичное приближение является, в достаточно малой окрестности точки a, лучшим приближением, чем линейное. В частности,

f(x)=P2(x)+h2(x)(xa)2,limxah2(x)=0.

Здесь ошибка приближения такова

R2(x)=f(x)P2(x)=h2(x)(xa)2 

которая, при ограниченном характере h2, приближается к нулю быстрее, чем приближается к нулю Шаблон:Nowrap по мере того, как x стремится к a.

Приближение функции f(x) = 1/(1 + x2) с помощью многочленов Pk порядка k = 1, …, 16 относительно точки x = 0 (красный) и точки x = 1 (салатовый цвет). Приближение вообще не улучшается за пределами (-1,1) и (1-√2,1+√2), соответственно.

Таким образом, мы будем продолжать получать более хорошие приближения к f, если будем использовать многочлены всё более высокой степени. В общем, ошибка в приближении функции с помощью полиномов порядка k будет приближаться к нулю немного быстрее, чем приближается к нулю Шаблон:Nowrap по мере того как x стремится к a.

Это следствие имеет асимптотическую природу: оно лишь говорит нам, что ошибка Rk приближения с помощью многочленов Тейлора k-го порядка Pk приближается к нулю быстрее, чем ненулевой многочлен k-го порядка по мере того как xa. Оно не говорит нам, насколько велика ошибка в любой окрестности центра приближения, но для этого существует формула для остатка (приведена ниже).

Наиболее полные версии теоремы Тейлора как правило приводят к равномерным оценкам ошибки приближения в малой окрестности центра приближения, но эти оценки не являются адекватными для окрестностей, которые слишком велики, даже если функция f является аналитической. В этой ситуации следует выбирать несколько многочленов Тейлора с разными центрами приближения, чтобы иметь надёжное Тейлорово приближение к исходной функции (см. Анимированный рисунок выше). Возможна также ситуация, когда возрастание порядка многочлена не увеличивает качество приближения вообще, даже если функция f дифференцируется бесконечное число раз. Такой пример приведён ниже.

Теорема Тейлора для функций от одной вещественной переменной

Формулировка теоремы

Точная формулировка большинства базовых версий теоремы такова.

Шаблон:Quotation

Многочлен, возникающий в теореме Тейлора, является многочленом Тейлора k-го порядка

Pk(x)=f(a)+f(a)(xa)+f(a)2!(xa)2++f(k)(a)k!(xa)k

функции f в точке a.

Теорема Тейлора описывает асимптотическое поведение остаточного члена

 Rk(x)=f(x)Pk(x),

который является ошибкой при нахождении приближения функции f с помощью многочленов Тейлора. Используя «O» большое и «o» малое теорему Тейлора можно сформулировать так

Rk(x)=o(|xa|k),xa.

Формулы для остатка

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

Шаблон:Quotation

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

Можно так же найти и другие выражения для остатка. Например, если G(t) является непрерывной на закрытом интервале и дифференцируемой с нестремящейся к нулю производной на открытом интервале между a и x, то

Rk(x)=f(k+1)(ξ)k!(xξ)kG(x)G(a)G(ξ)

для некоторого числа ξ между a и x. Эта версия охватывает формы Лагранжа и Коши как частные случаи, и выводится с помощью теоремы Коши о среднем значении (расширенной версии теоремы Лагранжа о среднем значении).

Запись формулы для остатка в интегральной форме является более общей, чем предыдущие формулы, и требует понимания интегральной теории Лебега. Однако она сохраняется также для интеграла Римана при условии, что производная порядка (k+1) от f является непрерывной на закрытом интервале [a,x].

Шаблон:Quotation

Вследствие абсолютной непрерывности f(k) на закрытом интервале между a и x, её производная f(k+1) существует как L1-функция, и это следствие может быть получено с помощью формальных вычислений с использованием теоремы Ньютона — Лейбница и интегрирования по частям.

Оценки остатка

На практике часто бывает полезно численно оценить величину остаточного члена приближения Тейлора.

Будем считать, что f является (k+1)-раз непрерывно дифференцируемой на интервале I, содержащем a. Будем считать, что существуют действительные постоянные числа q и Q такие, что

qf(k+1)(x)Q

на всём протяжении I. Тогда остаточный член удовлетворяет неравенству[1]

q(xa)k+1(k+1)!Rk(x)Q(xa)k+1(k+1)!,

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

|f(k+1)(x)|M

на интервале Шаблон:Nowrap с некоторым r>0, то

|Rk(x)|M|xa|k+1(k+1)!Mrk+1(k+1)!

для всех Шаблон:Nowrap Второе неравенство называется равномерной оценкой, потому что она сохраняет равномерность для всех x на интервале Шаблон:Nowrap

Пример

Приближение ex (голубой) с помощью многочленов Тейлора Pk порядка k=1,…,7 с центром в точке x=0 (красный).

Допустим, мы хотим найти приближение функции Шаблон:Nowrap на интервале Шаблон:Nowrap и убедиться, что ошибка не превышает значения 10−5. В этом примере считаем, что нам известны следующие свойства экспоненциальной функции:

(*)e0=1,ddxex=ex,ex>0,x.

Из этих свойств следует, что Шаблон:Nowrap для всех k, и в частности, Шаблон:Nowrap. Отсюда следует, что многочлен Тейлора k-го порядка функции f в точке 0 и его остаточного члена в форме Лагранжа даётся формулой

Pk(x)=1+x+x22!++xkk!,Rk(x)=eξ(k+1)!xk+1,

где ξ — это некоторое число между 0 и x. Поскольку ex возрастает согласно (*), мы можем использовать ex ≤ 1 для x ∈ [−1, 0], чтобы оценить остаток на подынтервале [−1, 0]. Для нахождения верхней границы значения остатка на интервале [0,1], можем использовать свойство Шаблон:Nowrap для 0<ξ<x, чтобы оценить

ex=1+x+eξ2x2<1+x+ex2x2,0<x1

используя многочлен Тейлора второго порядка. Выражая из этого неравенства ex, приходим к выводу, что

ex1+x1x22=21+x2x24,0x1

приняв, что числитель принимает максимальное из всех своих возможных значений, а знаменатель принимает минимальное из всех своих возможных значений. Используя эти оценки значений ex, мы видим, что

|Rk(x)|4|x|k+1(k+1)!4(k+1)!,1x1,

и требуемая точность определённо достигается в том случае, когда

4(k+1)!<1054105<(k+1)!k7.

(где факториал 7!=5 040 и 8!=40 320.) В конечном счёте, теорема Тейлора приводит к приближению

ex=1+x+x22!++x77!+R7(x),|R7(x)|<105,1x1.

Отметим, что это приближение позволяет вычислить значение e≈2.71828 с точностью до пятого знака после запятой.

Аналитичность

Разложение Тейлора для вещественных аналитических функций

Пусть I является открытым интервалом. По определению, функция f:I является вещественной аналитической, если она на данном участке определена сходимостью степенного ряда. Это означает, что для каждого aI существует некоторое r > 0 и последовательность коэффициентов ckR такая, что Шаблон:Nowrap и

f(x)=k=0ck(xa)k=c0+c1(xa)+c2(xa)2+,|xa|<r.

В общем, радиус сходимости степенного ряда может быть вычислен по Шаблон:Не переведено

1R=lim supk|ck|1k.

Этот результат основан на сравнении с бесконечно убывающей геометрической прогрессией, и тот же самый метод показывает, что если степенной ряд, разложенный по a, сходится для некоторого bR, он должен сходиться равномерно на закрытом интервале Шаблон:Nowrap, где rb = |ba|. Здесь мы только рассмотрели сходимость степенного ряда, и не исключено, что область Шаблон:Nowrap расширяется за пределы области определения I функции f.

Многочлен Тейлора от вещественной аналитической функции f в точке a

Pk(x)=j=0kcj(xa)j,cj=f(j)(a)j!

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

Rk(x)=j=k+1cj(xa)j=(xa)khk(x),|xa|<r.

Здесь функция

hk:(ar,a+r);hk(x)=(xa)j=0ck+1+j(xa)j

также является аналитической, поскольку её степенной ряд имеет тот же радиус сходимости, что и исходный ряд. При условии, что Шаблон:NowrapI и r < R, все эти ряды сходятся равномерно на интервале Шаблон:Nowrap. Конечно, в случае аналитических функций можно оценить остаточный член Rk(x) путём «обрезания» последовательности производных f′(a) в центре приближения, но при использовании комплексного анализа появляются и другие возможности, которые описаны ниже.

Теорема Тейлора и сходимость ряда Тейлора

Существует разногласие между многочленами Тейлора дифференцируемых функций и рядами Тейлора аналитических функций. Можно рассматривать (справедливо) ряд Тейлора

f(x)k=0ck(xa)k=c0+c1(xa)+c2(xa)2+

бесконечное число раз дифференцируемой функции f:RR как её «многочлен Тейлора бесконечно большого порядка» в точке a. Теперь оценка остатка многочлена Тейлора подразумевает, что для любого порядка k и для любого r>0 существует постоянная Шаблон:Nowrap такая, что

(*)|Rk(x)|Mk,r|xa|k+1(k+1)!

для каждого x∈(a-r, a+r). Иногда эти постоянные могут быть выбраны таким образом, что Шаблон:Nowrap, когда Шаблон:Nowrap и Шаблон:Nowrap остаётся неизменной. Тогда ряд Тейлора функции f сходится равномерно к некоторой аналитической функции

Tf:(ar,a+r);Tf(x)=k=0f(k)(a)k!(xa)k.

Тут важно упомянуть тонкий момент. Возможна ситуация, когда бесконечное число раз дифференцируемая функция f имеет ряд Тейлора в точке a, который сходится в некоторой открытой окрестности точки a, но предельная функция Tf отличается от f. Важным примером этого феномена является такой

f:;f(x)={e1x2,x>0,0,x0.

Используя цепное правило можно показать индуктивно, что для любого порядка k,

f(k)(x)={pk(x)x3ke1x2,x>00,x0

для некоторого многочлена pk. Функция e1x2 стремится к нулю быстрее, чем любой полином, по мере того как Шаблон:Nowrap, тогда f является бесконечное число раз дифференцируемой и Шаблон:Nowrap для каждого положительного целого k. Теперь оценки для остатка многочлена Тейлора функции f показывают, что ряд Тейлора сходится равномерно к нулевой функции на всей действительной числовой оси. Не будет ошибки в следующих утверждениях:

  • Ряд Тейлора функции f сходится равномерно к нулевой функции Tf(x)=0.
  • Нулевая функция является аналитической, и каждый коэффициент её ряда Тейлора равен нулю.
  • Функция f является бесконечное число раз дифференцируемой, но не аналитической.
  • Для любого kN и r>0 существует Mk, r>0 такое, что остаточный член многочлена Тейлора k-го порядка функции f удовлетворяет условию (*).

Теорема Тейлора в комплексном анализе

Теорема Тейлора обобщает функции f:, которые являются комплексно дифференцируемыми на открытом подмножестве UC комплексной плоскости. Однако её полезность снижена другими теоремами комплексного анализа, а именно: более полные версии подобных результатов могут быть выведены для комплексно дифференцируемых функций f : UC с использованием интегральной формулы Коши как показано ниже.

Пусть r > 0 такое, что замкнутый круг B(z, r) ∪ S(z, r) содержится в U. Тогда интегральная формула Коши с положительной параметризацией Шаблон:Nowrap окружности S(z, r) с Шаблон:Nowrap даёт

f(z)=12πiγf(w)wzdw,f(z)=12πiγf(w)(wz)2dw,,f(k)(z)=k!2πiγf(w)(wz)k+1dw.

Здесь все подынтегральные выражения являются непрерывными на окружности S(z, r), что обосновывает Шаблон:Не переведено. В частности, если f является один раз комплексно дифференцируемой на открытом множестве U, то она фактически бесконечное число раз комплексно дифференцируема на U. Имеем оценку Коши[2]

|f(k)(z)|k!2πγMr|wz|k+1dw=k!Mrrk,Mr=max|wc|=r|f(w)|

для любого zU и r > 0 такой, что B(z, r) ∪ S(c, r) ⊂ U. Эти оценки подразумевают, что комплексный ряд Тейлора

f(z)k=0f(k)(c)k!(zc)k

функции f сходится равномерно в любом круге B(c, r) ⊂ U с S(c, r) ⊂ U в некоторой функции Tf. Кроме того, используя формулу интегрирования по контуру для производных f(k)(c),

Tf(z)= k=0(zc)k2πiγf(w)(wc)k+1dw=12πiγf(w)wck=0(zcwc)kdw= 12πiγf(w)wc(11zcwc)dw=12πiγf(w)wzdw=f(z),

таким образом, любая комплексно дифференцируемая функция f на открытом множестве UC является комплексно аналитической. Всё то, что было написано выше для вещественных аналитических функций справедливо также и для комплексных аналитических функций, где открытый интервал I заменён на открытое подмножество UC и a-центрированные интервалы (ar, a + r) заменена на c-центрированные круги B(c, r). В частности, разложение Тейлора сохраняется в виде

f(z)=Pk(z)+Rk(z),Pk(z)=j=0kf(k)(c)k!(zc)k,

где остаточный член Rk является комплексно аналитическим. При рассмотрении рядов Тейлора методы комплексного анализа позволяют получить несколько более мощные результаты. Например, используя интегральную формулу для любого положительно ориентированную жорданову кривую γ которая параметризирует границу ∂WU области WU, можно получить выражение для производных Шаблон:Nowrap как показано выше, и слегка изменив расчёты для Шаблон:Nowrap, прийти к точной формуле

Rk(z)=j=k+1(zc)j2πiγf(w)(wc)j+1dw=(zc)k+12πiγf(w)dw(wc)k+1(wz),zW.

Важная особенность здесь состоит в том, что качество приближения с помощью многочлена Тейлора в области WU является мажорируемым значениями функции f на границе ∂WU. Так же, применяя оценки Коши к выражению остатка Ряда, получаем равномерные оценки

|Rk(z)|j=k+1Mr|zc|jrj=Mrrk+1|zc|k+11|zc|rMrβk+11β,|zc|rβ<1.

Пример

График комплексной функции f(z) = 1/(1 + z2). Модуль показан высотой подъёма и аргумент показан цветом: циан=0, синий=π/3, фиолетовый=2π/3, красный=π, жёлтый=4π/3, зелёный=5π/3.

Функция f:RR, определяемая уравнением

f(x)=11+x2

является вещественной аналитической, то есть, в данной области определяется её рядом Тейлора. Один из рисунков, приведённых выше, показывает, что некоторые очень просто задаваемые функции не могут быть выражены с помощью приближения Тейлора в окрестности центра приближения, если эта окрестность слишком велика. Это свойство легко понять в рамках комплексного анализа. Более конкретно, функция f расширяется до мероморфной функции

f:{}{};f(z)=11+z2

на компактифицированной комплексной плоскости. Она имеет простые оси в точках z=i и z=−i, и она всюду аналитическая. Её ряд Тейлора, имеющий центром z0, сходится на любом круге B(z0,r) с r<|z-z0|, где тот же ряд Тейлора сходится при zC. Вследствие этого ряд Тейлора функции f, имеющий центром точку 0, сходится на B(0,1) и он не сходится для любого zC с |z|>1 вследствие имеющихся осей в точках i и −i. По тем же причинам ряд Тейлора функции f, имеющий центром точку 1, сходится на B(1,√2) и не сходится для любого zC с |z-1|>√2.

Обобщения теоремы Тейлора

Высшие порядки дифференцируемости

Функция f:RnR является дифференцируемой в точке aRn тогда и только тогда, когда существует линейная форма L : RnR и функция h : RnR такая, что

f(𝒙)=f(𝒂)+L(𝒙𝒂)+h(𝒙)(𝒙𝒂),lim𝒙𝒂h(𝒙)=0.

Если этот случай имеет место, то L = df(a) является дифференциалом функции f в точке a. Кроме того, когда частные производные функции f существуют в точке a, то дифференциал f в точке a даётся формулой

df(𝒂)(𝒗)=fx1(𝒂)v1++fxn(𝒂)vn.

Вводя мультииндекс, запишем

|α|=α1++αn,α!=α1!αn!,𝒙α=x1α1xnαn

для αNn и xRn. Если все частные производные k-го порядка функции Шаблон:Nowrap являются непрерывными в Шаблон:Nowrap, то, по теореме Клеро, можно изменить порядок смешанных производных в точке a, тогда запись

Dαf=|α|fx1α1xnαn,|α|k

для частных производных высших порядков является правомерной в этой ситуации. То же самое является верным, если все частные производные (k − 1)-го порядка функции f существуют в некоторой окрестности точки a и являются дифференцируемыми в точке a. Тогда можно сказать, что функция f является k раз дифференцируемой в точке a .

Теорема Тейлора для функций многих переменных

Шаблон:Quotation

Если функция Шаблон:Nowrap является k+1 раз непрерывно дифференцируемой в замкнутом шаре B, то можно получить точную формулу для остатка разложения Тейлора до частных производных Шаблон:Nowrap порядка от f в этой окрестности. А именно

f(𝒙)=|α|=0kDαf(𝒂)α!(𝒙𝒂)α+|β|=k+1Rβ(𝒙)(𝒙𝒂)β,Rβ(𝒙)=|β|β!01(1t)|β|1Dβf(𝒂+t(𝒙𝒂))dt.

В этом случае, вследствие непрерывности частных производных (k+1)-го порядка на компактном множестве B, непосредственно получаем

|Rβ(𝒙)||β|β!max|α|=|β|max𝒚B|Dαf(𝒚)|,𝒙B.

Доказательства

Доказательство теоремы Тейлора для одной вещественной переменной

Пусть[3]

hk(x)={f(x)P(x)(xa)kx=a0x=a

где, как указано в формулировке теоремы Тейлора,

P(x)=f(a)+f(a)(xa)+f(a)2!(xa)2++f(k)(a)k!(xa)k.

Достаточно показать, что

limxahk(x)=0.

Доказательство основано на повторяющемся применении правила Лопиталя. Заметим, что каждое Шаблон:Nowrap, f(j)(a)=P(j)(a). Отсюда каждая следующая производная числителя функции hk(x) стремится к нулю в точке x=a, и то же самое справедливо для знаменателя. Тогда

limxaf(x)P(x)(xa)k=limxaddx(f(x)P(x))ddx(xa)k==limxadk1dxk1(f(x)P(x))dk1dxk1(xa)k=1k!limxaf(k1)(x)P(k1)(x)xa=1k!(f(k)(a)P(k)(a))=0

где переход от предпоследнего выражения к последнему следует из определения производной в точке x = a.

Примечания

Шаблон:Примечания

Источники

Ссылки