Ряд Тейлора

Материал из testwiki
Перейти к навигации Перейти к поиску

Ряд Те́йлора — разложение функции в бесконечную сумму степенных функций. Частный случай разложения в ряд Тейлора в нулевой точке называется рядом Маклорена.

Ряд Тейлора был известен задолго до публикаций Брука Тейлора[1] — его использовали ещё в XIV веке в Индии[2], а также в XVII веке Грегори и Ньютон.

Ряды Тейлора применяются при аппроксимации функции многочленами. В частности, линеаризация уравнений происходит путём разложения в ряд Тейлора и отсечения всех членов выше первого порядка.

Обобщением понятия ряда Тейлора в функциональном анализе является ряд Фантапье.

Определение

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

n=0kf(n)(a)n!(xa)n=f(a)+f(a)(xa)+f(2)(a)2!(xa)2++f(k)(a)k!(xa)k,

используемая в приближённых вычислениях, как обобщение следствия теоремы Лагранжа о среднем значении дифференцируемой функции:

при xa=h0 верно f(x)=f(a+h)=f(a)+f(a)h+O(h2)f(a)+f(a)h=f(a)+f(a)(xa).

При записи суммы использованы обозначение f(0)(x)=f(x) и соглашение о произведении по пустому множеству: 0!=1, (xa)0=1.

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

f(x)=n=0+f(n)(a)n!(xa)n=n=0+φn(x;a) с общим членом φn(x;a)=f(n)(a)n!(xa)n, зависящим от параметра a.

Другими словами, рядом Тейлора функции f(x) в точке a называется ряд разложения функции по положительным степеням двучлена (xa):

f(x)=f(a)+f(a)(xa)+f(a)2!(xa)2++f(n)(a)n!(xa)n+.[3]

Как указано ниже в примерах, наличия бесконечной дифференцируемости функции f(x) в окрестности точки a не достаточно, чтобы ряд Тейлора сходился к самой функции где-либо, кроме самой точки a.

3. Рядом Тейлора в точке a функции f(z) комплексной переменной z, удовлетворяющей в некоторой окрестности U точки a условиям Коши — Римана, называется степенной ряд

f(z)=n=0+f(n)(a)n!(za)n.

В отличие от вещественного случая, из условий следует, что найдётся такое значение радиуса R>0, что в DR={z:|zz0|<R}U ряд сходится к функции f(z).

4. В случае a=0 ряд

f(x)=n=0+f(n)(0)n!xn

называется рядом Маклорена.

Аналитическая функция

1. Функция f(x) вещественной переменной x называется аналитической в точке x=a, если существуют такой радиус R>0 и такие коэффициенты ck=ck(a)=ck(a;f), k=0,1,2,, что f(x) может быть представлена в виде сходящегося на интервале (aR;a+R) степенного ряда: k=0+ck(xa)k, то есть x(aR;a+R) limn+k=0nck(xa)k=f(x).

Функция называется аналитической на промежутке (на множестве), если она является аналитической в каждой точке этого промежутка (множества).

2. Степенной ряд k=0+ck(za)k на любом компактном подмножестве K области сходимости DR={z:|zz0|<R} допускает почленное дифференцирование любое количество раз.

Если в k-ю производную функции k=0+ck(za)k подставить z=a, то получится ckk!.

Таким образом, для аналитической в точке a функции f(z) для некоторого R>0 всюду в DR={z:|zz0|<R} является верным представление f(z)=k=0+f(k)(a)k!(za)k.

Следствие. Функция f(x) вещественной переменной x является аналитической в точке a тогда и только тогда, когда она равна своему ряду Тейлора с параметром a на некотором открытом интервале, содержащем точку a.

3. Вопрос: будет ли для произвольной бесконечно дифференцируемой в точке a функции f(x) вещественного переменного x её ряд Тейлора k=0+f(k)(a)k!(xa)k сходиться к f(x) всюду на каком-нибудь интервале (aR;a+R), то есть представима ли f(x) этим рядом?

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

Примеры. Функции вещественной переменной f2(x)={e1x2,x00,x=0, f+(x)={e1x,x>00,x0, fv(x)={e1|x|,x00,x=0 являются бесконечно дифференцируемыми в точке x=0, причём все эти производные равны нулю.

Следовательно, ряды Тейлора всех этих функций с параметром a=0 тождественно равны нулю. Однако, для любого R>0 в окрестности (R;+R) точки a=0 найдутся точки, в которых функции отличны от 0. Таким образом, эти функции не являются в точке a=0 аналитическими.

Шаблон:Начало скрытого блока

Доказательство проведём для функции f(x)=f2(x)={e1x2,x00,x=0, предложенной Огюстеном Луи Коши.

Функция exp(1z2), является аналитической функцией комплексной переменной для всех z{0}.

Для z0 очевидно, что ddzexp(1z2)=exp(1z2)(2z3).

Функция f(x) для x — это «исправленная» функция exp(1x2), x{0}, дополненная пределами слева limx0,x<0exp(1x2)=0 и справа limx0,x>0exp(1x2)=0 в точке x=0.

Найдём производную функции f(x) в точке x=0. По определению: f(0)=limΔx0,Δx{0}f(0+Δx)f(0)Δx=limh0,h{0}f(h)0h=00=limh0,h{0}f(h)h=limh0,h{0}2f(h)h3.

Поскольку для x(0;1) выполняется 0<e1x2<e1x, то докажем, что для произвольного α>0 верно limx0,x>0e1xxα=0.

Применение правила Лопиталя непосредственно к частям

limx0,x>0e1x=limx0,x>0xα=0 не приводит к результату.

Выполним замену переменной: 1x=t:

limx0,x>0e1xxα=limt+tαet=++=limt+αtα1et.

Пусть k=α. Применяя правило Лопиталя k раз, в числителе получим либо (при α=k) константу k!, либо (при α<k) бесконечно малую α(α1)(αk+1)tαk:

limt+tαet=++==limt+α(α1)(αk+1)tαket=0.

Таким образом,

f(0)=limh0,h{0}2f(h)h3=0.

Найдём (для x0) несколько начальных производных функции f(x):

f(x)=2f(x)x3
f(x)=(2f(x)x3)=2(f(x)1x3+f(x)(1x3))=2(2f(x)x31x3+f(x)(1x3))=2f(x)(2x63x4)
f(x)=(2f(x)(2x63x4))=4f(x)(2x93x7+6x56x7)

И так далее. Во всех случаях, очевидно, получается произведение f(x) на сумму целых отрицательных степеней x. Конечная сумма бесконечно малых является бесконечно малой. Таким образом, limx0,x{0}f(k)(x)=0.

Вычисляя последовательно по определению (как выше) производные f(x) в точке x=0, обнаруживаем, что все производные в точке x=0 равны нулю. Шаблон:Конец скрытого блока

Примером гладкой функции, не являющейся аналитической ни в одной точке своей области определения, служит функция Фабиуса.

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

Ряд Тейлора, являясь степенным рядом, имеет в качестве области сходимости круг (с центром в точке a) для случая комплексной переменной и интервал (с центром в точке a) — для случая вещественной переменной.

1. Например, функция f(x)=11x может быть разложена в ряд Тейлора следующим образом: 11x=k=0xk (это известная формула суммы бесконечной убывающей геометрической прогрессии). Однако если функция 11x определена для всех действительных чисел, кроме точки x=1, то ряд k=0xk сходится только при условии |x|<1.

2. Радиус сходимости ряда Тейлора можно определить, например, по формуле Даламбера:

R=limk|f(k)(a)k!f(k+1)(a)(k+1)!|=limk|f(k)(a)f(k+1)(a)(k+1)|.

3. Рассмотрим для примера экспоненциальную функцию ex. Поскольку любая производная экспоненциальной функции равна самой функции в любой точке, то радиус сходимости экспоненциальной функции равен R=limk|eaea(k+1)|=limk(k+1)=. Значит, ряд Тейлора экспоненциальной функции сходится на всей оси x для любого параметра a.


4. От параметра — точки разложения a ряда Тейлора — зависит область его сходимости.

Например, разложим в общем случае (для произвольного a) в ряд Тейлора функцию f(x)=11x: f(x)=11x=11ak=0(xa1a)k.

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

Действительно,

11ak=0(xa1a)k=11a11(xa1a)=11x.

Область сходимости ряда может быть задана неравенством |xa1a|<1. И теперь эта область зависит от a. Например, для a=0 ряд сходится при x(1;1). Для a=0,5 ряд сходится при x(0;1).

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

Предположим, что функция f(x) имеет все производные до n+1-го порядка включительно в некотором промежутке, содержащем точку x=a. Найдем многочлен Pn(x) степени не выше n, значение которого в точке x=a равняется значению функции f(x) в этой точке, а значения его производных до n-го порядка включительно в точке x=a равняются значениям соответствующих производных от функции f(x) в этой точке.

Достаточно легко доказать, что такой многочлен имеет вид Pn(x)=k=0nf(k)(a)k!(xa)k, то есть это n-я частичная сумма ряда Тейлора функции f(x). Разница между функцией f(x) и многочленом Pn(x) называется остаточным членом и обозначается Rn(x)=f(x)Pn(x). Формула f(x)=Pn(x)+Rn(x) называется формулой Тейлора[4]. Остаточный член дифференцируем n+1 раз в рассматриваемой окрестности точки a. Формула Тейлора используется при доказательстве большого числа теорем в дифференциальном исчислении. Говоря нестрого, формула Тейлора показывает поведение функции в окрестности некоторой точки.

Теорема: Шаблон:Рамка Если функция f(x) имеет n+1 производную на отрезке с концами a и x, то для произвольного положительного числа p найдётся точка ξ, лежащая между a и x, такая, что

f(x)=k=0nf(k)(a)k!(xa)k+(xaxξ)p(xξ)n+1n!pf(n+1)(ξ).

Шаблон:Конец рамки

Это формула Тейлора с остаточным членом в общей форме (форма Шлёмильха — Роша).

Различные формы остаточного члена

В форме Лагранжа:

Rn(x)=(xa)n+1(n+1)!f(n+1)[a+θ(xa)]p=n+1;0<θ<1

В форме Коши:

Rn(x)=(xa)n+1(1θ)nn!f(n+1)[a+θ(xa)]p=1;0<θ<1

В интегральной форме:

Rn(x)=1n!ax(xt)nf(n+1)(t)dt

Шаблон:Начало скрытого блока

Методом интегрирования по частям получим Rn(x)=1n!ax(xt)nf(n+1)(t)dt=1n!ax(xt)ndf(n)(t)=1n!((xt)nf(n)(t))|ax1n!axf(n)(t)d(xt)n==1(n1)!ax(xt)n1f(n)(t)dt(xa)nf(n)(a)n!=...=axf(t)dtk=1nf(k)(a)(xa)kk!=f(x)k=0nf(k)(a)(xa)kk!
откуда
f(x)=k=0nf(k)(a)(xa)kk!+Rn(x)

Шаблон:Конец скрытого блока Ослабим предположения:

  • Пусть функция f(x) имеет n1 производную в некоторой окрестности точки a и n-ю производную в самой точке a, тогда:
В асимптотической форме (форме Пеано, локальной форме):
Rn(x)=o[(xa)n]

Шаблон:Начало скрытого блока

Поскольку Rn(a)=Rn(a)=Rn(a)=...=Rn(a)(n)=0, то предел отношения Rn(x)(xa)n при x, стремящемся к a, может быть найден по правилу Лопиталя: limxaRn(x)(xa)n=limxaRn(x)((xa)n)=limxaRn(x)((xa)n)=...=limxaRn(x)(n)((xa)n)(n)=Rn(a)(n)n!=0
Поскольку исходный предел равен нулю, это значит, что остаточный член Rn(x) является бесконечно малой функцией более высокого порядка, чем (xa)n, при xa. А это и есть определение о-малого.

Шаблон:Конец скрытого блока

Критерий аналитичности функции

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

Во-первых, необходимым условием аналитичности функции является сходимость ряда Тейлора в некоторой непрерывной области. Действительно, если ряд Тейлора сходится всего в одной точке, то это точка x=a, потому что в ней ряд Тейлора сходится всегда. Но тогда ряд Тейлора равен функции f(x) только в этой единственной точке, а значит, данная функция не будет аналитической.

Во-вторых, по формуле Тейлора в ряд Тейлора с остаточным членом может быть разложена любая (а не только аналитическая) функция, бесконечно дифференцируемая в окрестности, содержащей точку a. Пусть ряд Тейлора с параметром a такой функции сходится в этой окрестности. Если существует предел каждой из двух последовательностей, то предел суммы этих последовательностей равен сумме их пределов. Тогда для всех x из окрестности a по формуле Тейлора можно записать limnRn(x)=limn(f(x)Pn(x))=f(x)limnPn(x), где limnPn(x) — ряд Тейлора.

Очевидно, что функция f(x) является аналитической в точке a тогда и только тогда, если в указанной окрестности точки a существует непрерывная область X такая, что для всех xX остаточный член её разложения по формуле Тейлора стремится к нулю с ростом n: limnRn(x)=0.

В качестве примера рассмотрим экспоненциальную функцию ex. Её ряд Тейлора сходится на всей оси x для любых параметров a. Докажем теперь, что эта функция является аналитической во всех точках a.

Остаточный член разложения этой функции в форме Лагранжа имеет вид Rn(x)=(xa)n+1(n+1)!eξn, где ξn — некоторое число, заключенное между x и a (не произвольное, но и не известное). Тогда, очевидно,

limnRn(x)=limn(xa)n+1(n+1)!eξnMlimn(xa)n+1(n+1)!=0

Здесь используется, что на фиксированном промежутке экспонента ограничена некоторым числом M

Причем, как видно, предел остаточного члена равен нулю для любых x и a.

Ряды Маклорена некоторых функций

  • Экспонента: ex=1+x1!+x22!+x33!+=n=0xnn!,x.
  • Натуральный логарифмряд Меркатора»): ln(1+x)=xx22+x33=n=0(1)nxn+1(n+1)=n=1(1)n1xnn для всех 1<x1.
  • Биномиальное разложение: (1+x)α=1+n=1(αn)xn, для всех |x|<1 и всех комплексных α, где (αn)=k=1nαk+1k=α(α1)(αn+1)n! — обобщённые биномиальные коэффициенты.
    • Квадратный корень[5]: 1+x=1+12x1124x2+113246x311352468x4+11357246810x5=n=0(1)n+1(2n)!22n(2n1)(n!)2xn=n=0(1)n+1(2n+1)!!(2n)!!(2n+1)(2n1)xn для всех |x|1.
    • Обратный квадратный корень[5]: 11+x=112x+1324x2135246x3+13572468x413579246810x5+=n=0(1)n(2n)!22n(n!)2xn=n=0(1)n(2n+1)!!(2n)!!(2n+1)xn для всех 1<x1.
    • Шаблон:Iw:
      • 11x=1+x+x2+x3+=n=0xn для всех |x|<1.
      • 1(1x)2=1+2x+3x2+4x3+=n=1nxn1 для всех |x|<1.
      • 1(1x)3=1+3x+6x2+10x3+=n=2(n1)n2xn2 для всех |x|<1.
      • Конечный геометрический ряд: 1xm+11x=n=0mxn для всех x=1, m0.
  • Тригонометрические функции[5][6]:
    • Синус: sinx=xx33!+x55!x77!+x99!=n=0(1)nx2n+1(2n+1)!,x.
    • Косинус: cosx=1x22!+x44!x66!+x88!=n=0(1)nx2n(2n)!,x.
    • Тангенс: tg x=x+13x3+215x5+17315x7+=n=1(1)n+122n(22n1)B2n(2n)!x2n1 для всех |x|<π2, где B2n — числа Бернулли.
    • Котангенс: ctg x=x113x145x32945x5+=n=0(1)n22nB2n(2n)!x2n1 для всех 0<|x|<π, где B2n — числа Бернулли.
    • Секанс: secx=1+12x2+524x4+61720x6+=n=0(1)nE2n(2n)!x2n для всех |x|<π2, где E2n — числа Эйлера.
    • Косеканс: cosecx=x1+16x+7360x3+3115120x5+=n=0(1)n+12(22n11)B2n(2n)!x2n1 для всех 0<|x|<π, где B2n — числа Бернулли.
  • Обратные тригонометрические функции[5][7]:
    • Арксинус: arcsinx=x+123x3+13245x5+1352467x7+ =n=0(2n+1)!!(2n)!!x2n+1(2n+1)2=n=0(2n)!22n(n!)2x2n+12n+1 для всех |x|1[8].
    • Арккосинус: arccosx=π2arcsinx=π2x123x313245x51352467x7=π2n=0(2n+1)!!(2n)!!x2n+1(2n+1)2=π2n=0(2n)!22n(n!)2x2n+12n+1 для всех |x|1.
    • Арктангенс: arctgx=xx33+x55x77+ =n=1(1)n12n1x2n1 для всех |x|1.
    • Арккотангенс: arcctgx=π2arctgx=π2x+x33x55+x77=π2n=1(1)n12n1x2n1 для всех |x|1.
  • Гиперболические функции[5][9]:
    • Гиперболический синус: shx=x+x33!+x55!+x77!+=n=0x2n+1(2n+1)!,x.
    • Гиперболический косинус: chx=1+x22!+x44!+x66!+=n=0x2n(2n)!,x.
    • Гиперболический тангенс: thx=xx33+2x51517x7315+=n=122nB2n(22n1)(2n)!x2n1 для всех |x|<π2.
    • Гиперболический котангенс: cthx=x1+13x145x3+2945x5=n=022nB2n(2n)!x2n1 для всех 0<|x|<π.
    • Гиперболический секанс: sechx=112x2+524x461720x6+=n=0E2n(2n)!x2n для всех |x|<π2.
    • Гиперболический косеканс: cosechx=x116x+7360x33115120x5+=n=02(22n11)B2n(2n)!x2n1 для всех 0<|x|<π.
  • Обратные гиперболические функции[5][10]:
    • Гиперболический арксинус: arshx=x123x3+13245x51352467x7+ =n=0(1)n(2n+1)!!(2n)!!x2n+1(2n+1)2=n=0(1)n(2n)!22n(n!)2x2n+12n+1 для всех |x|1.
    • Гиперболический арктангенс: arthx=x+13x3+15x5+17x7+=n=0x2n+12n+1 для всех |x|<1.
  • W-функция Ламберта: W0(x)=xx2+3x328x43+125x524=n=1(n)n1n!xn,|x|1/e.
  • Другие функции:
    • ex1x=1+n=1en!n!xn[11]

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

Пусть функция f(x,y) имеет непрерывные производные до (n+1)-го порядка включительно в некоторой окрестности точки (x0,y0). Введём дифференциальный оператор

T=(xx0)x+(yy0)y.

Тогда разложение (формула Тейлора) функции f(x,y) по степеням (xx0)p(yy0)q для p+qn в окрестности точки (x0,y0) будет иметь вид

f(x,y)=k=0nTkf(x0,y0)k!+Rn(x,y),

где Rn(x,y) — остаточный член в форме Лагранжа:

Rn(x,y)=T(n+1)f(ξ,ζ)(n+1)!, ξ[x0,x], ζ[y0,y]

Следует иметь в виду, что операторы x и y в Tk действуют только на функцию f(x,y), но не на (xx0) и/или (yy0).

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

В случае функции одной переменной T=(xx0)ddx.

Формула Тейлора многих переменных

Для получения формулы Тейлора функции n переменных f(x1,x2,...xn), которая в некоторой окрестности точки (a1,a2,...,an) имеет непрерывные производные до (m+1)-го порядка включительно, введём дифференциальный оператор

T=(x1a1)x1+(x2a2)x2+...+(xnan)xn.

Тогда разложение (формула Тейлора) функции по степеням (xiai)ki в окрестности точки (a1,a2,...,an) имеет вид

f(x1,x2,...xn)=k=0mTkf(a1,a2,...,an)k!+Rm(x1,x2,...xn),

где Rm(x1,x2,...xn) — остаточный член порядка (m+1).

Для функции n переменных, бесконечно дифференцируемой в некоторой окрестности точки (a1,a2,...,an), ряд Тейлора имеет вид:

f(x1,x2,...xn)=k=01k!i1=1ni2=1n...ik=1nkf(a1,a2,...,an)xi1xi2...xik(xi1a1)(xi2a2)...(xinan).

В другой форме ряд Тейлора можно записать таким образом:

f(x1,x2,...xn)=k=0k1=0k2=0...kn=0k1+k2+...+kn=k1k1!k2!...kn!kf(a1,a2,...,an)x1k1x2k2...xnkn(x1a1)k1(x2a2)k2...(xnan)kn.

Пример разложения в ряд Маклорена функции трёх переменных

Найдём выражение для разложения в ряд Тейлора функции трёх переменных x, y и z в окрестности точки (0,0,0) до второго порядка малости. Оператор T будет иметь вид

T=xx+yy+zz.

Разложение в ряд Тейлора запишется в виде

f(x,y,z)=k=02Tkf0k!+R2(x,y,z)=
=(1+T+T22)f0+R2(x,y,z);

Учитывая, что

T2=x22x2+y22y2+z22z2+2xy2xy+2xz2xz+2yz2yz,

получим

f(x,y,z)=f0+xf0x+yf0y+zf0z+x222f0x2+y222f0y2+z222f0z2+
+xy2f0xy+xz2f0xz+yz2f0yz+R2(x,y,z).

Например, при f(x,y,z)=ex+y+z,

f(x,y,z)=1+x+y+z+x22+y22+z22+xy+xz+yz+R2(x,y,z).

Примечания

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

Литература

Шаблон:Последовательности и ряды Шаблон:Дифференциальное исчисление

  1. Taylor, Brook, Methodus Incrementorum Directa et Inversa [Direct and Reverse Methods of Incrementation] (London, 1715), pages 21-23 (Proposition VII, Theorem 3, Corollary 2). Translated into English in D. J. Struik, A Source Book in Mathematics 1200—1800 (Cambridge, Massachusetts: Harvard University Press, 1969), pages 329—332.
  2. Gupta R. C. The Madhava-Gregory series, Math. Education 7 (1973), B67-B70.
  3. Запорожец Г. И. «Руководство к решению задач по математическому анализу» — С. 371
  4. Шаблон:Книга
  5. 5,0 5,1 5,2 5,3 5,4 5,5 Градштейн И. С., Рыжик И. М. Таблицы интегралов, сумм, рядов и произведений Шаблон:Wayback. — 4-е изд. — М.: Наука, 1963.
  6. Шаблон:АбрамовицСтиган
  7. Шаблон:АбрамовицСтиган
  8. При значении Шаблон:Math, близком к 1, эта расчётная формула сходится медленно, т.е. даёт большую погрешность при приближении функции суммой первых нескольких членов ряда. Поэтому можно воспользоваться формулой arcsinx=arccos1x2, где arccosx=π2arcsinx.
  9. Шаблон:АбрамовицСтиган
  10. Шаблон:АбрамовицСтиган
  11. Шаблон:Cite web