Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
108
Добавлен:
26.05.2014
Размер:
361.3 Кб
Скачать
an, мы можем утверждать,

ВАВТ; 2 курс; м. ан; весна 2003; shurik-f.narod.ru

Л.10. Произвольные числовые ряды.

Ряды с произвольные членами. Абсолютная и условная сходимость.

Ряд, в котором есть как положительные, так и отрицательные члены, называется знакопеременным рядом (не путать со знакочередующимся рядом, речь о котором ниже). Сложность ситуации в этом случае (в отличие от положительных рядов) состоит в том, что последовательность частичных сумм теперь не обязана возрастать; если положительные ряды могли только или сходиться, или расходиться к бесконечности, то произвольные ряды могут и просто расходиться

(пример: ряд 1 − 1 + 1 − 1 + . . .). Более того, если сходимость положительного ряда является следствием того, что общий

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

Чтобы разобраться во всем этом, нам нужно ввести два понятия.

+∞

Ðÿä P

n=1

+∞

Ðÿä P

n=1

 

+∞

 

nP

an называется сходящимся абсолютно, если сходится ряд

=1 |an|.

an называется сходящимся условно, если сам он сходится, а ряд из абсолютных величин расходится.

Если рассуждать чисто логически, надо рассмотреть следующие четыре случая.

1.

Ðÿäû

 

 

сходятся;

 

 

 

 

 

 

an è

an

 

 

 

2.

Ðÿäû

Pan è

P||an||

 

 

 

 

 

 

 

3.

Ðÿä

 

расходятся;

 

 

 

 

 

a сходится, а ряд

|

a

n|

расходится;

4.

Ðÿä

 

Pn

P

 

 

 

 

 

 

 

Pan расходится, а

ðÿä

 

 

|an|

сходится.

 

 

P

 

 

 

 

 

P

 

 

 

P

 

 

Проще всего с третьей ситуацией: она в точности описывается термином условная сходимость . Кстати, сам термин намекает на наше скептическое отношение к такого рода сходимости: вроде бы ряд сходится, но сходится натужно,условно ; можно сказать, что такой ряд сходится из-за того, что члены ряда ведут себя примерно как рак и щука ; а вот когда они объединяют свои усилия и начинают тянуть в одну сторону, сумма оказывается бесконечностью.

Чтобы разобраться с остальными случаями, докажем теорему.

Теорема о сходимости абсолютно сходящегося ряда. Если сходится ряд из модулей, то и сам ряд сходится.

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

Доказательство. Пусть M и N суммы положительных слагаемых ряда и модулей отрицательных слагаемых соответственно. Тогда сумма ряда из модулей всех членов ряда равна M + N; поскольку по предположению ряд сходится абсолютно,

P

M +N конечное число, следовательно, M и N тоже конечные числа. Возвращаясь к ряду что сумма этого ряда конечна и равна M − N.

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

Возвращаемся к условно сходящимся рядам. Поскольку ряд из модулей в этом случае расходится, M +N = ∞. Какими

могут быть тогда M è N? Если они конечны, то это означает абсолютную сходимость. Если M = ∞, à N является

конечным числом (или наоборот), то это означает что в одну сторону суммарно мы сдвигаемся бесконечно далеко, а в другую на конечный отрезок, следовательно, ряд будет расходиться к бесконечности. И только если M = N = ∞,

возможна условная сходимость (а возможна и расходимость).

 

 

 

 

 

 

 

Полезно держать в голове примеры абсолютно сходящихся рядов и условно сходящихся рядов. Пример абсолютно

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

 

P

(−1)

n

P

1

 

минусы, например, через один. Так, абсолютно сходящимся будет ряд

 

 

 

 

 

забыли его? Пример условно сходящегося ряда: P (−1)n

 

 

 

 

 

2n

 

, поскольку ряд

2n сходится. Не

n . То, что ряд из модулей расходится, очевидно ведь получается гармонический ряд. То, что сам ряд сходится (и тем самым сходится условно), Вы чуть позже сможете доказать самостоятельно с помощью признака Лейбница.

Свойства сходящихся и абсолютно сходящихся рядов.

1. Åñëè ðÿäû

 

P

P λA +P

an è

bn сходятся (и их суммы равны A è B соответственно), то и ряд

(λan + µbn) сходится, причем

его сумма равна

µB.

 

Это свойство очень легко доказывается (примерно так же, как и аналогичное свойство для определенных интегралов): надо перейти к пределу в равенстве, связывающем частичные суммы рядов.

21

{cn} удовлетворяет следующим условиям:

2.(Это свойство я уже приводил в конце 8-й лекции; сейчас я его записываю для полноты коллекции). Если у ряда изменить, отбросить или добавить конечное число членов, то измененный ряд будет сходиться или расходиться в зависимости от того, сходится или расходится первоначальный ряд.

Далее свойства абсолютно сходящихся рядов; на условно сходящиеся рядов они не переносятся!

3.Если произвольным образом переставить местами члены абсолютно сходящегося ряда, то он не перестанет быть абсолютно сходящимся рядом с той же самой суммой.

4. Åñëè ðÿäû +∞ an è

+∞ bn абсолютно сходятся (к своим суммам A è B), то ряд, составленный из всевозможных

P

k P

 

 

 

n=1

n=1

 

 

 

попарных произведений a bm, расположенных в произвольном порядке, также абсолютно сходится (к AB).

Такой ряд называется

рядов

P

P

Свойства 3 и 4 доказывать на экзамене не

 

 

произведением

 

an è

bn.

 

 

обязательно.

Знакочередующиеся ряды. Признак Лейбница. Оценка остатка знакочередующегося ряда.

Знакочередующимся называется ряд, у которого соседние члены имеют разные знаки. Считая, что первый член ряда больше нуля, можем записать ряд в виде

c1 − c2 + c3 − . . . ,

ãäå c1, c2, . . . положительные числа.

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

Признак Лейбница. Предположим, что последовательность 1. она монотонно убывает (то есть c1 > c2 > c3 > . . .);

2. lim cn = 0.

n→+∞

Тогда знакочередующийся ряд c1 − c2 + c3 − . . . + (−1)n+1cn + . . . сходится.

Короче этот признак можно сформулировать так: знакочередующийся ряд, общий член которого по абсолютной величине монотонно стремится к нулю, сходится. Конечно, в этом признаке не утверждается, что исследуемый ряд будет сходиться абсолютно, хотя не утверждается и обратное; образно говоря, признак гарантирует, что исследуемый знакочередующийся ряд, удовлетворяющий условиям 1 2, сходится как минимум условно.

Доказательство. Изображаем частичные суммы ряда на числовой прямой. Так как мы предположили, что c1 > 0, первая частичная сумма S1 = c1 находится справа от нуля. Так как c1 > c2 > 0, вторая частичная сумма S2 = c1 − c2 будет находиться слева от S1, но справа от нуля. Далее, при переходе от S2 к S3 мы снова делаем шаг вправо, но меньший, чем перед этим шаг налево. Поэтому S3 будет находиться между S2 и S1. Аналогично доказываем, что S4 находится

между S2 и S3. Таким образом, частичные суммы с четными номерами образуют возрастающую последовательность, а частичные суммы с нечетными номерами убывающую:

S2 < S4 < . . . < . . . < S3 < S1

Образное сравнение: навстречу друг другу несутся два поезда, причем расстояние между ними стремится к нулю: S2n−1 − S2n = c2n → 0. Вывод: где-то эти поезда встретятся, а у наших частичных сумм, что четных, что нечетных, один и тот

же предел. Этот предел будет пределом всей последовательности частичных сумм, что и доказывает, что ряд сходится. Кстати, иногда ряд, удовлетворяющий условиям признака Лейбница, называется рядом Лейбница. Тогда признак

Лейбница может быть сформулирован так:

Ряд Лейбница сходится.

+∞

Заканчивая лекцию, дадим оценку остатка ряда Лейбница. Напоминание. Остатком ðÿäà P an называется ряд,

n=1

полученный из первоначального отбрасыванием первых n членов (то есть тех, которые входят в Sn):

+∞

X

ak,

k=n+1

а также сумма этого ряда. Ясно, что если ряд сходится (к числу S), а Rn обозначение его n-го остатка, то Sn + Rn = S,

поэтому остаток сходящегося ряда стремится к нулю при n, стремящемся к бесконечности. В случае ряда Лейбница про остаток ряда можно утверждать также следующее:

|Rn| < cn+1

Иными словами, остаток Rn ряда Лейбница по модулю меньше модуля первого отброшенного члена. Задача. Доказать последнее утверждение самостоятельно.

22

областью сходимости

ВАВТ; 2 курс; м. ан; весна 2003; shurik-f.narod.ru

Л.11. Функциональные ряды.

Функциональный ряд. Область сходимости.

Ðÿä

u1(x) + u2(x) + . . . + un(x) + . . . ,

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

ýòîò x во все функции, получаем уже числовой ряд. Поэтому к функциональному ряду нужно относиться как к множеству

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

тем самым от x, сумма функционального ряда является функцией (обозначим ее через S(x)), определенной на области сходимости этого ряда. Естественно частичные суммы функционального ряда обозначить через Sn(x).

Правильная сходимость функционального ряда.

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

Говорят, что функциональный ряд u1(x) + u2(x) + . . . + un(x) + . . . сходится правильно на множестве D (естественно, лежащем в области сходимости ряда), если существует положительный сходящийся числовой ряд c1 + c2 + . . . + cn + . . .

такой, что для каждого x D è n N выполнено

|un(x)| ≤ cn.

Смысл этого понятия состоит в том, что с помощью ряда P

cn мы контролируем скорость убывания членов функционального ряда сразу для всех x D. Какие полезные выводы можно сделать, если есть правильная сходимость, мы узнаем совсем

скоро, но сначала пример, когда правильной сходимости нет, и это приводит к прискорбным последствиям.

Пример. D = [0; 1]; u1(x) = x; u2(x) = x2 −x; u3(x) = x3 −x2; . . .; un(x) = xn −xn−1; . . .. Тогда Sn(x) = u1(x) + u2(x) +

. . . un(x) = x + x2 − x + x3 − x2 + . . . xn − xn−1 = xn. Следовательно,

n→+∞

n

n→+∞

1,

S(x) = lim

S (x) =

lim xn =

0,

åñëè

åñëè

x[0; 1) x = 1

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

Непрерывность суммы правильно сходящегося ряда. Почленное интегрирование и дифференцирование функциональных рядов.

1. Сумма правильно сходящегося на множестве D функционального ряда, состоящего из непрерывных функций, является непрерывной на D функцией.

2. Если функции un(x) непрерывны на отрезке [a, b], причем ряд из этих функций сходится правильно на [a, b], òî ðÿä

Zax u1(x) dx + Zax u2(x) dx + . . . + Zax un(x) dx + . . .

 

 

также сходится правильно на [a, b], причем если S(x) сумма исходного ряда, то

ax S(x) dx сумма ряда из интегралов.

 

un(x) имеют непрерывные

R

 

 

un(x) сходится, а ряд,

составленный из производных, сходится правильно в области D, òî ðÿä

 

un(x) можно

 

P

3. Более сложное свойство. Если функции

 

 

производные, ряд

 

что означает, что суммой ряда из производных является производная

P

 

почленно дифференцировать,

 

 

 

 

 

суммы первоначального ряда:

S0(x) = u01(x) + u02(x) + . . . u0n(x) + . . .

23

ВАВТ; 2 курс; м. ан; весна 2003; shurik-f.narod.ru

Л.12. Степенные ряды.

 

 

Степенной ряд. Теорема Абеля.

 

 

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

 

Степенной ряд ýòî ðÿä âèäà

 

 

c0 + c1(x − a) + c2(x − a)2 + .

. . + cn(x − a)n + . . . ,

(3)

ãäå cn è a постоянные числа. Частный случай такого ряда ряд

 

c0 + c1x + c2x2 + . . .

+ cnxn + . . .

(4)

Изучение произвольного степенного ряда можно всегда свести к изучению этого более простого ряда: для этого нужно сделать замену x − a = t. Давайте так и поступим.

Первое простое наблюдение: такой ряд всегда сходится как минимум в одной точке: при x = 0, получается ряд

c0 + 0 + 0 + . . .. Как оказывается, иногда область сходимости на этом исчерпывается, иногда ряд сходится на всей

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

Теорема Абеля. Если ряд (4) сходится в точке x0, то он сходится (причем абсолютно) во всех точках, находящихся на более близком расстоянии от начала координат, чем точка x0.

Иными словами, если |x| < |x0|, то в точке x ряд сходится. Задача. Сформулировать теорему Абеля для ряда (3).

Для простоты докажем эту теорему в предположении, что в точке x0 ряд сходится абсолютно. Собственно говоря, в этом случае и доказывать нечего: абсолютная сходимость в точке x мгновенно следует из признака сравнения: |cnxn| =

|cn| · |x|n ≤ |cn| · |x0|n = |cnxn0 |. Желающие могут найти доказательство в общем случае в умных книжках или же попробовать доказать самостоятельно. Между прочим, это не так уж и сложно.

Простое следствие из теоремы Абеля: если ряд (4) расходится в точке x0, то он расходится во всех точках, находящихся на большем расстоянии от начала координат, чем x0. Продумайте самостоятельно это утверждение.

Вывод. Область сходимости степенного ряда (4) может быть только следующей: {0}, (−∞; +∞); [−R; R], (−R; R),

(−R; R] èëè [−R; R). При подготовке к экзамену попробуйте придумать примеры таких рядов. Считая точку отрезком

нулевой длины, а числовую прямую интервалом бесконечной длины, видим, что область сходимости всегда задается числом R (условно будем считать, что бесконечность также является числом) таким, что интервал (−R; R) входит в нее,

интервалы (−∞; −R) è (R; +∞) не входят, а концевые точки могут как входить, так и не входить.

Интервал и радиус сходимости степенного ряда. Нахождение радиуса сходимости. Область сходимости степенного ряда.

Радиус сходимости степенного ряда это то самое число R, о котором шла речь в предыдущем абзаце. Немного повторяясь, можно сказать, что радиус сходимости это такое число R, ÷òî åñëè |x| < R, то ряд сходится, а если |x| > R, то расходится. Интервал (−R; R) называется интервалом сходимости. Таким образом, область сходимости и интервал сходимости это почти одно и то же, да не совсем: концевые точки −R è R могут входить в область сходимости (одна

или обе), а могут не входить.

 

 

 

 

 

 

 

Иногда для нахождения области сходимости сначала ищут радиус (а вместе с ним интервал) сходимости. Предположим,

что мы исследуем ряд с помощью признака Даламбера. Имеем: D = lim

˛

cn+1xn+1

˛

= |x| lim

˛

cn+1

˛. Обратите внимание,

cnxn

cn

 

˛

 

˛

 

˛

 

˛

 

˛

 

˛

 

˛

 

˛

 

˛

 

˛

 

˛

 

˛

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

сказать, сходится или не сходится ряд в этой точке: если D < 1, òî åñòü |x| lim ˛

cn

˛

< 1, иными словами, |x| < lim

˛cn+1

˛

(обратили внимание на то, дробь перевернулась?), то ряд сходится, если же

˛

cn+1

 

 

 

cn

î˛

D >˛

1, то ряд расходится. Это говорит˛

 

˛

 

˛

 

˛

 

˛

том, что радиус сходимости можно попробовать поискать по формуле

˛

 

˛

 

˛

 

˛

R = lim

˛cn+1

˛

 

˛

cn

˛

 

˛

 

˛

 

˛

 

˛

Рассуждая аналогично (используя признак Коши), получаем еще один способ нахождения радиуса сходимости:

1 R = lim p

n |cn|

Правильная сходимость степенного ряда на любом отрезке, лежащем внутри интервала сходимости.

Теорема. Степенной ряд сходится правильно на любом отрезке, целиком лежащем в интервале сходимости.

24

[−a; a] (−R; R). Выберем точку x0,

Внимание! Здесь говорится об отрезке, лежащем в интервале сходимости, а не в области сходимости. Этим подчеркивается, что при приближении к концевым точкам ряд может вести себя все хуже и хуже (сходиться все натужнее). То, что мы не имеем права приближаться к концевым точкам сколь угодно близко, в доказательстве используется самым непосредственным образом.

Доказательство. Пусть требуется доказать правильную сходимость на отрезке

лежащую строго между a è R. Благодаря теореме Абеля мы можем утверждать, что не только ряд

cnx0n

сходится,

 

|cnx0

| сходится. Если точка x [−a; a]

 

|x| < |x0|

 

|cnx

n

| < |cnx0 |. Âñå

 

P

 

íî è ðÿä

n

 

 

 

 

n

 

 

 

P

 

 

, òî

 

, а тогда

 

 

 

вышесказанное и означает

 

 

 

 

 

 

 

правильную сходимость на отрезке [−a; a].

Возникает естественный вопрос, так ли уж важно, чтобы между отрезком [−a; a] и концевыми точками оставался

зазор, или это требуется только при нашем способе доказательства? Иными словами, обязана ли быть правильная сходимость на всей области сходимости или хотя бы на интервале сходимости? Ответом на это может послужить следующий пример. Рассмотрим (надеюсь всем известный) ряд Pxn/n, радиус сходимости которого равен единице,

а область сходимости является полуинтервалом [−1; 1). Если взять отрезок [−1; a] при любом a (−1; 1), то на нем максимальное значение модуля n-го члена равно 1/n, а ряд из этих чисел расходится (ведь это же знаменитый гармонический

ряд!), и поэтому наш ряд не может сходиться правильно; нет правильной сходимости и на интервале (−1; 1). Следующий вопрос, который может возникнуть: может быть, правильной сходимости на всей области сходимости нет

ни у какого степенного ряда? Ответ дает ряд

P

xn/n2, областью сходимости которого служит отрезок [ 1; 1], причем

сходится он всюду правильно, поскольку на

 

 

˛xn/n2˛

≤ 1/n2, àряд с такими

 

области сходимости справедливо неравенство

˛

˛

 

 

 

 

 

 

членами сходится.

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

Подчеркнем, теорема говорит о непрерывности даже в концевых точках (конечно, если там ряд сходится). Впрочем, доказательство мы проведем только для точек интервала сходимости.

Итак, пусть точка x0 лежит в интервале сходимости. Рассмотрим отрезок [−a; a] такой, что точка x0 лежит внутри

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

Между прочим, смотрите, как интересно: на всем интервале сходимости правильной сходимости может не быть, однако непрерывной сумма степенного ряда на интервале сходимости будет всегда!

Почленное интегрирование и дифференцирование степенных рядов.

 

 

 

 

Теорема. Ряд, составленный из производных членов степенного ряда, имеет тот же интервал сходимости.

Следует отметить, что при этом область сходимости меняться может.

 

Px

n−1/n сходится на полуинтервале

x

 

P

n 2 сходится на отрезке

[

 

1; ; ряд из производных

 

 

 

 

 

 

 

Пример. Ряд

 

x /n

 

1]

 

 

 

x

 

[ 1; 1); ðÿä

P

n/n сходится на полуинтервале [

 

 

1; 1); ряд из производных

P

 

n−1

(−1; 1).

 

 

 

 

 

 

 

 

 

 

 

Проведем доказательство в предположении, что существует lim ˛

cn

˛ (аналогичное доказательство можно провести,

 

 

 

 

 

 

 

 

˛

cn+1

˛

 

 

 

 

 

 

 

 

 

 

 

˛

 

˛

 

 

 

p

если существует предел lim n |cn|; как Вы уже конечно сообразили, в первом случае при доказательстве используется признак Даламбера, а во втором признак Коши); в общем случае доказательство сложнее.

Обозначим радиусы сходимости исходного ряда и ряда из производных через R1 è R2 соответственно. Если исходный

ðÿä èìåë âèä Pcnxn, то продифференцированный Pncnxn−1, тогда

R1

= lim

˛cn+1

˛

; R2 = lim ˛

(n + 1)cn+1

˛

= lim n + 1 lim

˛cn+1

˛

= lim

˛cn+1

˛ = R1

 

 

˛

cn

˛

˛

ncn

˛

 

n

˛

cn

˛

 

˛

cn

˛

 

 

˛

 

˛

˛

 

˛

 

 

 

˛

 

˛

 

˛

 

˛

 

 

˛

 

˛

˛

 

˛

 

 

 

˛

 

˛

 

˛

 

˛

Доказательство завершено.

Из теорем о почленном интегрировании и почленном дифференцировании правильно сходящегося функционального ряда следуют аналогичные теоремы для степенных рядов.

Теорема. В интервале сходимости степенной ряд можно почленно интегрировать.

Наиболее красивая формулировка звучит так: для любого x из интервала сходимости ряда Pcnxn с суммой S(x) ðÿä

x

n

 

P(cn/(n + 1))x

n+1 сходится, и его сумма равна

x

(x) dx.

 

PR0 cnx

 

dx =

+∞

n сходится на интервале

 

R0

S, и его сумма (как сумма бесконечной геометрической

Пример. Как известно, ряд

n=0 x

 

x + x

 

/2 + x

 

 

(−1; 1)

( 1; 1), и его сумма равна

прогрессии) равна 1/(1 x);

 

P

 

2

3

/3 + . . .

 

 

 

 

 

поэтому ряд

 

 

 

 

 

 

также сходится на интервале

 

 

 

 

 

 

 

 

 

 

 

− ln(1 − x).

Теорема. В интервале сходимости степенной ряд можно почленно дифференцировать.

Пример. Дифференцируя ряд из предыдущего примера, получаем, что сумма ряда 1 + 2x + 3x2 + 4x3 + . . . равна

1/(1 − x)2.

Замечание. Поскольку любой степенной ряд можно почленно продифференцировать без изменения интервала сходимости, никто не может нам помешать повторить эту приятную процедуру дюбое число раз. Таким образом, сумма степенного ряда является бесконечно дифференцируемой функцией в интервале сходимости ряда.

25

ВАВТ; 2 курс; м. ан; весна 2003; shurik-f.narod.ru

Л.13. Ряд Тейлора.

Формула Тейлора с остаточным членом в форме Лагранжа.

Первоначальные наблюдения. Рассмотрим многочлен P (x) = 2x3

 

5x2 + 7x + 9, записав его в обратном порядке: P (x) =

9 + 7x − 5x2 + 2x3. Найдем значения этого многочлена и его

производных в нуле:

P (0) = 9 (это его свободный член);

 

 

 

 

 

 

 

 

P 0(x) = 6x2

10x+7; P 0(0) = 7 (это коэффициент при x); P 00(x) = 12x

10; P 00(0) =

 

10 (это удвоенный коэффициент

 

 

 

 

 

 

 

 

 

 

 

 

 

очевидно, пойдут одни нули. Поэтому

ïðè x2

);

 

 

 

 

 

 

 

 

 

 

); далее,

 

 

 

P 000(x) = 12x; P 000(0) = 12 (это ушестеренный коэффициент при x3

 

 

 

 

многочлен можно записать в виде P (x) = P (0) + P 0(0)x +

P 00 (0)

x2 +

P 000 (0)

x3.

 

 

 

 

 

 

 

 

 

2

 

 

 

 

6

 

 

 

 

 

 

Аналогичный вывод можно сделать относительно любого многочлена: если P (x) = c0 + c1x + c2x2 + . . . + cnxn, òî

P (0) = c0; P 0(0) = c1; P 00(0) = 2c2; . . .; P (n)(0) = n!cn; таким образом, P (x) = P (0) + P 0(0)x +

P 00(0)

x2

+ . . . +

P (n)

xn.

2

 

 

 

 

n!

Для тех, кто забыл, напоминаю, что n! это не эн восклицательный знак , а n-факториал , то есть произведение натуральных чисел от 1 äî n:

n! = 1 · 2 · . . . · n

Возникает естественный вопрос: что же, производные в нуле важнее, чем в других точках? Можно ли представить многочлен в таком виде, чтобы использовались производные не в нуле, а в какой-то другой точке x0? Оказывается,

 

, òî

x − x0. Проведите сами такое

 

P (n)(x0)

можно, только тогда многочлен нужно разложить по степеням

 

рассуждение, доказав, что

åñëè P (x) = b0 + b1(x − x0) + b2(x − x0)2 + . . . + bn(x − x0)n

 

b0 = P (x0); b1 = P 0(x0); . . .; bn =

 

 

 

n! .

А теперь переходим к общему рассуждению.

 

 

 

 

 

Пусть нам дана функция f(x), которую мы хотим в окрестности точки x0 хотя бы приближенно представить в виде

многочлена. Если про функцию известно лишь, что она непрерывна, то таким приближением служит значение функции в этой точке:

f(x) ≈ f(x0);

геометрически это означает, что график функции f(x) в окрестности точки x0 мы приближаем горизонтальной прямой y = f(x0). Если непрерывна не только сама функция, но и ее первая производная, то можно приблизить функцию точнее,

если проводить не горизонтальную прямую, а касательную к графику в точке x0, которая, как известно, имеет уравнение y = f(x0) + f0(x0)(x − x0). Таким образом, получаем более точное приближение нашей функции:

f(x) ≈ f(x0) + f0(x0)(x − x0).

Кстати, можно получить не приближенное, а абсолютно точное выражение функции f(x) с помощью производной; для этого надо применить теорему Лагранжа (помните, во втором семестре у Вас были Ролль-Лагранж-Коши?):

f(x) = f(x0) + f0(ξ)(x − x0),

ãäå ξ некоторая точка, лежащая между x0 è x. Конечно, если мы меняем точку x, точка ξ также меняется.

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

f(x) = c0 + c1(x − x0) + c2(x − x0)2 + c3(x − x0)3 + . . . + cn(x − x0)n + . . .

Чему равны коэффициенты cn? Рассуждение повторяет то, которое было в начале лекции. Подставляем вместо x число

x0, в результате все слагаемые, за исключением первого, обратятся в ноль. Следовательно, f(x0) = c0. Дифференцируя ряд (помните, на прошлой лекции мы говорили о возможности дифференцирования степенного ряда?), получаем

f0(x) = c1 + 2c2(x − x0) + 3c3(x − x0)2 + . . . + ncn(x − x0)n−1 + (n + 1)cn+1(x − x0)n + . . .

Поэтому f0(x0) = c1. Аналогично, сосчитав n-ю производную и подставив в нее x0, получаем f(n)(x0) = n! · cn.

Итак, если функцию можно разложить в степенной ряд по степеням (x − x0), то это разложение обязательно следующее:

f(x) = f(x0) + f0

f00

(x0)

(x − x0)2 +

f000(x0)

(x − x0)3 + . . . +

fn(x0)

(x − x0)n + . . .

(x0)(x − x0) +

 

 

 

 

 

 

 

 

2

 

6

 

n!

Но пусть ряд пока подождет; разберемся сначала с конечными суммами. Предположим, что функция f(x) в точке x0

имеет производные до n-й включительно (конечно, отсюда следует, что производные с первой по

(n − 1)-ую существуют

даже в некоторой окрестности этой точки (чуть позже мы потребуем, чтобы существовала даже

n + 1-я производная, и

не только в точке x0, но и в ее окрестности, но пока это не требуется). Составим сумму

 

 

 

f00

(x0)

 

f000(x0)

 

fn(x0)

 

Tn(x) = f(x0) + f0(x0)(x − x0) +

 

 

 

(x − x0)2 +

 

 

(x − x0)3 + . . . +

 

 

(x − x0)n

 

2

 

6

 

n!

26

(n + 1)!
достаточные условия

Естественно, хочется надеяться, что эта сумма мало отличается от функции f(x), хотя бы в очень маленькой окрестности x0, но в общем случае это утверждать нельзя. Независимо от этого, если разность между f(x) è Tn(x) обозначить через Rn(x), получается формула

f(x) = Tn + Rn

Эта формула называется формулой Тейлора для функции f(x), Tn(x) многочленом Тейлора, а Rn(x) остаточным

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

Теорема. Если в некоторой окрестности точки x0 функция f(x) имеет (n + 1)-ю производную, то остаточный член формулы Тейлора для точек этого интервала может быть записан в виде

f(n+1)(ξ)

Rn(x) = (n + 1)! (x − x0)n+1,

ãäå ξ некоторая точка, лежащая между x0 è x.

Остаточный член, записанный в таком виде, называется остаточным членом в форме Лагранжа. Естественно, для каждого x точка ξ, вообще говоря, своя. Обратите внимание, что остаточный член в форме Лагранжа имеет форму,

аналогичную всем членам многочлена Тейлора, только производная считается не в точке x0, а в некоторой другой точке.

Ряд Тейлора и Маклорена.

Пусть дана функция f(x), имеющая в точке x0 производные âñåõ порядков. Тогда мы можем составить степенной ряд

f(x0) + f0

f00

(x0)

(x − x0)2 +

f000(x0)

(x − x0)3 + . . . +

f(n)(x0)

(x − x0)n + . . . ,

(x0)(x − x0) +

 

 

 

 

 

 

 

2

 

6

 

n!

который называется рядом Тейлора функции f(x) по степеням (x − x0). В частном случае, когда x0 = 0, ряд Тейлора

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

функции расходится во всех точках, кроме x0. Все эти вопросы выходят за пределы нашего курса.

Условия разложимости функции в ряд Тейлора.

Сначала очевидное утверждение.

Утверждение, гордо именуемое критерием разложимости функции в ряд Тейлора. Ряд Тейлора функции f(x)

сходится к f(x) тогда и только тогда, когда остаточный член формулы Тейлора стремится к нулю. Более содержательными являются такой разложимости.

Теорема. Если существуют окрестность точки x0 и число C, такие, что все производные функции f(x) существуют в этой окрестности и в этой окрестности по абсолютной величине не превосходят C, то ряд Тейлора сходится к f(x) â

этой окрестности.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство. Запишем остаточный член формулы Тейлора в форме Лагранжа (впрочем, никаких других форм у

нас в курсе не было): Rn(x) =

 

f(n+1)(ξ)

(x

 

x0)n+1. Имеем:

 

Rn(x)

 

=

f(n+1)(ξ)

(x

 

x0)n+1

 

=

|f(n+1)(ξ)|

 

x

 

x0

n+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(n + 1)!

 

|

|

| (n + 1)!

|

|

 

C |x − x0|n+1

 

 

 

 

 

 

 

 

 

 

 

 

 

(n + 1)!

 

|

 

0 ïðè n

+

. Если последнее утверждение кажется не вполне очевидным, можно объяснить его так:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

обозначим для простоты |x − x0| через a. Åñëè a < 1, то утверждение просто очевидно, поскольку тогда числитель стремится к нулю (а знаменатель стремится к бесконечности), и поэтому дробь и подавно стремится к нулю. Если

a = 1, то дробь стремится к нулю из-за стремления знаменателя к бесконечности. Если же a > 1, то ситуация сложнее:

к бесконечности стремятся и числитель, и знаменатель, поэтому предел зависит от того, кто стремится к бесконечностибыстрее (если Вы помните, в такой ситуации говорят о наличии неопределенности, в данном случае типа ). Äëÿ

того, чтобы понять, что в нашем случае предел равен нулю, будем брать сразу достаточно большие n, например, n > 2a;

тогда при переходе от n ê n + 1 числитель дроби умножается на a, а знаменатель на число, большее 2a, тем самым дробь уменьшается более, чем в 2 раза. Надеюсь, что теперь никто не сомневается в утверждении.

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

Разложим в ряд Маклорена некоторые из элементарных функций. Начнем с ex.

 

 

 

1. f(x) = ex; x0 = 0; f(0) = e0 = 1; f0(x) = ex; f0(0) = e0 = 1, и т. д. Следовательно, получаем ряд 1 + x +

x2

+

x3

+

 

 

. . . + xn + . . . = +∞ xn

 

2

6

 

 

0! понимают единицу. Докажем, что этот ряд сходится к ex на всей числовой прямой.

n!

nP

, ãäå ïîä

 

 

 

 

 

n!

 

 

 

 

 

 

=0

 

 

 

 

 

 

27

Зафиксируем для этого произвольное R > 0 и рассмотрим отрезок [

 

R; R]. Поскольку все производные нашей функции

равны ex, все они, будучи во-первых положительными (поэтому

модуль от них можно не брать), не превосходят на

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

отрезке [−R; R] числа e

R. Следовательно, в силу только что доказанной теоремы ряд Тейлора сходится к

x на отрезке

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x сходитсяeê

e

x íà âñåé

[−R; R]. Поскольку R произвольное число, мы можем утверждать, что ряд Тейлора функции e

 

 

 

 

 

 

 

 

числовой прямой.

 

 

 

 

 

 

f0(x)

 

cos x; f0(0) = 1; f00(x)

 

 

sin x; f00(0) = 0; f000(x) =

 

 

 

 

 

 

2. f(x) = sin x; x

 

= 0

;

 

;

=

=

cos x;

 

000

 

 

 

 

0

 

f(0) = sin 0 = 0

 

 

 

 

,

 

 

 

 

 

 

 

x

3

 

x

5

 

 

 

 

; далее снова повторяются числа

 

,

 

,

 

 

, и так далее, до бесконечности. Получается ряд

 

 

 

 

 

 

 

f

(x) = −1

0

1

0

−1

x −

6

 

+

120

 

2n+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. . . + (−1)n

x

+ . . .. Он сходится к sin x на всей числовой прямой, так как производные синуса, являясь плюс-минус

(2n+1)!

синусами и косинусами, не превосходят по модулю единицы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3. f(x) = cos x. Здесь можно рассуждать двумя способами: или повторить тот же путь, который мы прошли, исследуя

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

 

2

 

 

4

6

+ . . . + (−1)n

2n

случае получится ряд 1 − x2

+ x24

x

x

+ . . .

720

(2n)!

 

1

 

 

 

 

 

 

 

 

 

 

 

 

4.

f(x) =

 

 

. Здесь ситуация совсем простая, если помнить сумму бесконечной геометрической прогрессии: 1 +

1 − x

x + x2

+ x3 + . . . + xn + . . . =

 

 

1

 

 

 

 

 

|x| < 1; при прочих x ряд расходится.

1

x при условии, что

 

1

 

 

 

 

 

 

 

 

 

 

 

5.

f(x) =

 

. Ясно, что эта функция является суммой ряда 1 − x + x2 − x3 + . . . + (−1)nxn + . . .

1 + x

6.

f(x) = (1 + x)α. Здесь α любое действительное (не обязательно целое) число. Если α = n натуральное число,

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

(1 + x)1 = 1 + x; (1 + x)2 = 1 + 2x + x2; (1 + x)3 = 1 + 3x + 3x2 + x3.

Для произвольного натурального n имеем:

(1 + x)n = 1 + nx +

n(n − 1)

x2

+

n(n − 1)(n − 2)

x3

+ . . . +

n(n − 1)(n − 2) . . . (n − k + 1)

xk

+ . . . + xn

2

 

6

 

 

k!

 

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

В самом деле, для f(x) = (1 + x)α имеем: f(0) = 1; f0(;x) = α(1 + x)α−1; f0(0) = α; f00(x) = α(α − 1)xα−2; f00(0) = α(α − 1);

;

(n)

 

= α(α

− 1)(α − 2) . . . (α − n + 1)(1 + x)

α−n

 

(n)(0) = α(α

1)(α

2) . . . (α

n + 1); . . .

 

. . . f

 

(x)

 

 

f α имеет вид

 

 

 

 

 

Следовательно, ряд Тейлора функции f(x) = (1 + x)

 

 

 

 

 

 

 

 

 

 

 

1 + αx +

α(α − 1)

x2 +

α(α − 1)(α − 2)

x3 + . . . +

α(α − 1)(α − 2) . . . (α − n + 1)

xn

+ . . .

 

 

 

 

2

 

6

 

 

 

 

 

 

n!

 

 

 

 

Называется этот ряд биномиальным рядом. Доказывать, что он сходится к (1 + x)α, мы не будем, принимая это на

веру, а вот интервал сходимости найти можно. Исключаем случай натурального n, когда биномиальный ряд состоит из

конечного числа членов и поэтому сходится к (1 + x)α всюду. Для всех прочих n радиус сходимости равен 1. Пожалуй,

доказательство этого я оставлю Вам в качестве задачи. Разложения, которые нужно знать:

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

x2

 

x3

 

 

 

 

 

xn

 

 

 

 

 

 

 

(x (−∞; +∞))

 

 

 

 

 

 

e

 

 

= 1 + x +

 

 

 

 

+

 

 

 

 

 

+ . . . +

 

 

+ . . .

 

 

 

 

 

 

 

 

 

 

 

2

 

 

6

 

n!

 

 

 

 

 

 

 

 

 

 

 

 

x3

 

 

 

 

x5

 

 

 

 

 

 

 

 

 

 

 

 

n x2n+1

 

 

 

 

 

 

 

 

 

 

 

 

 

sin x = x −

 

 

 

+

 

 

 

 

 

− . . . + (−1)

 

 

 

 

+ . . .

(x (−∞; +∞))

 

 

 

 

 

6

120

 

(2n + 1)!

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

x4

 

 

 

 

 

x6

 

 

 

 

 

 

 

 

 

 

 

n x2n

 

 

 

 

 

 

 

 

 

 

 

cos x = 1 −

 

 

 

 

+

 

 

 

 

 

 

+ . . . + (−1)

 

+ . . .

(x (−∞; +∞))

 

 

 

 

 

2

24

720

(2n)!

 

 

 

 

 

 

 

 

 

 

1

 

 

= 1 + x + x

2

+ x

3

+ . . . + x

n

+ . . . (x (−1; 1))

 

 

 

 

 

 

 

 

 

1 − x

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

3

 

 

 

 

 

 

 

 

n

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 1

− x + x − x + . . . + (−1) x

 

 

+ . . .

(x (−1; 1))

 

 

 

 

 

 

1 + x

 

 

 

 

 

 

(1 + x)α = 1 + αx +

α(α − 1)

x2

+

α(α − 1)(α − 2)

x3 + . . . +

α(α − 1)(α − 2) . . . (α − n + 1)

xn + . . . (x

 

(

1; 1))

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

x3

 

 

 

 

 

x4

 

 

 

 

 

 

 

 

 

 

 

 

xn

 

 

 

 

 

 

 

ln(1 + x) = x −

 

 

+

 

 

 

 

 

 

+ . . . + (−1)n−1

 

+ . . .

(x (−1; 1))

 

 

 

 

 

2

3

4

 

n

 

 

 

 

Мы здесь приводим интервалы сходимости; попытайтесь разобраться сами с поведением этих рядов на концах.

 

Задача. Пользуясь тем, что arctg x = R0x

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + x2 , доказать, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x3

 

x5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arctg x = x −

 

+

 

 

− . . .

 

 

 

 

(x [−1; 1])

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

5

 

 

 

 

 

 

 

 

 

28

ВАВТ; 2 курс; м. ан; весна 2003; shurik-f.narod.ru

Л.14. Приложения степенных рядов.

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

Оценка погрешностей.

Приближенное вычисление значений функций.

Пример 1. Вычислим число e с точностью до 10−5. Для этого пишем ряд Тейлора для функции ex:

ex = 1 + x +

x2

+

x3

+ . . . +

xn

+ . . . ,

 

 

n!

2

6

 

 

подставляем в него вместо x единицу:

e = 1 + 1 + 12 + 16 + . . . + n1! + . . . ,

после чего начинаем думать, сколько слагаемых нужно оставить, чтобы добиться нужной точности. Имеем: |Rn(1)| =

1

 

 

1

 

1

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

1

 

 

 

1

 

 

 

 

 

+

 

 

+ . . . =

 

 

 

 

1 +

 

 

 

+

 

 

 

 

+ . . .« <

 

 

 

1 +

 

+

 

+

 

 

 

+ . . .«.

 

(n + 1)!

(n + 2)!

(n + 1)!

n + 2

(n + 2)(n + 3)

(n + 1)!

n + 2

(n + 2)2

(n + 2)3

В скобке у нас получилась геометрическая прогрессия со знаменателем

q =

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. Для простоты заменим его на большее

n + 2

число на

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 , после чего прогрессию просуммируем. Напомним, что сумма бесконечной убывающей геометрической

 

 

 

 

прогрессии равна первому члену, деленному на 1 минус знаменатель , поэтому в нашем случае получается

1

 

= 2.

 

1 − 1/2

Поэтому |Rn(1)| <

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10−5

, подберем

n так, чтобы

2

 

 

< 10−5

, òî

(n + 1)! . Поскольку ошибка должна быть меньше

 

 

 

 

 

 

 

(n + 1)!

 

åñòü (n + 1)! > 200000. Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1! = 1; 2! = 2; 3! = 2! · 3 = 6; 4! = 3! · 4 = 24; 5! = 4! · 5 = 120;

6! = 5! · 6 = 720; 7! = 6! · 7 = 5040;

8! = 7! · 8 = 40320;

9! = 8! · 9 = 362880 > 200000.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, n + 1 = 9; n = 8, то есть с точностью 10−5 число e равно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e = 1 + 1 +

1

 

+

1

+

 

1

+

 

1

+

 

1

+

1

 

+

1

 

= 2 +

 

28961

= 2.7182787 . . .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5040

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

6

 

24

120

720

 

40320

 

 

40320

 

 

 

 

 

 

 

 

 

 

 

 

 

Ошибка меньше 9!1 = 0.0000027 . . ., поэтому e (2.718276; 2.71828144). Тем самым мы можем гарантировать четыре правильные цифры после запятой, пятая же может быть как семеркой, так и восьмеркой. Для справки:

e = 2.718281828459045 . . .

.

Как видите, сложнее всего было сделать правильную оценку остатка ряда.

Пример 2. Вычислить число π с точностью до 10−5.

Воспользуемся выписанным на прошлой лекции разложением функции arctg x:

 

x3

x5

arctg x = x −

 

+

 

− . . . (x [−1; 1])

3

5

Подставим, например, x = 1. Òàê êàê arctg 1 =

π

 

 

 

 

 

 

 

 

4 , получаем

 

 

 

 

 

 

 

 

 

 

 

 

π = 4 −

4

+

4

4

+ . . .

 

 

 

 

 

 

3

5

7

Поскольку этот ряд является рядом Лейбница (знакочередующийся и общий член стремится к нулю), |Rn| меньше модуля

первого отброшенного члена. Поэтому n надо выбирать из условия

4

<

1

 

2n + 1 > 400000. Таким

2n + 1

100000 , òî åñòü

образом, можно отбросить все слагаемые, начиная с

4

. Просуммировав оставшиеся, получим число π с нужной

400001

точностью.

 

 

 

 

 

 

 

Задача. Довести пример до конца (шутка).

29

Замечание. Существуют ряды, сходящиеся к числу π намного быстрее.

Для справки:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π = 3.1415926536 . . .

 

 

 

 

 

Приближенное вычисление определенных интегралов.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 1.

Вычислить интеграл

1

ex − 1

dx с точностью до 10−5.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

x2

 

x3

xn

R0

x

ex

1

 

 

 

x x2

 

 

xn−1

 

 

e

 

= 1 + x +

 

+

 

+ . . . +

 

+ . . .

 

 

 

 

= 1 +

 

+

 

 

+ . . . +

 

 

 

+ . . .

 

 

 

 

 

 

 

 

 

 

n!

 

 

 

2

6

 

n!

 

 

 

 

 

x

 

 

 

2

 

6

 

 

 

 

 

 

Почленно интегрируя этот ряд, получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

ex − 1

dx = 1 +

 

1

 

+

 

 

1

 

+ . . . +

1

 

+ . . .

 

 

 

 

 

 

 

 

 

 

 

· 2

 

 

·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z0

x

 

 

 

 

2

3

6

 

n · n!

Оценим остаток ряда:

 

 

 

 

(n + 1)!

 

n + 1 + (n + 2)(n + 2) +

 

(n + 3)(n + 2)(n + 3)

|Rn| < (n + 1)(n + 1)! + (n + 2)(n + 2)! + . . . =

 

 

 

 

1

 

 

 

1

 

 

 

1

1

 

 

1

 

 

1

 

(n + 1)!

n + 1

+ (n + 1)2

+ (n + 1)3 + . . .«

=

(n + 1)! 1 − 1/(n + 1) =

(n + 1)! n

= n · (n + 1)! .

1

 

1

1

 

1

 

 

 

1

 

1/(n + 1)

1 1

 

1

 

«

+ . . . <

Имеем:

1 · 2! = 2; 2 · 3! = 12; 3 · 4! = 72;

4 · 5! = 600; 5 · 6! = 3600;

6 · 7! = 30240;

7 · 8! = 282240 > 105.

 

 

 

 

 

 

 

 

 

 

 

 

 

Поэтому

1

ex − 1

dx =

1

+

1

+

 

1

+

1

 

+

1

 

+

 

1

 

 

 

 

=

45331

= 0.55993386 . . . Ошибка меньше

 

1

 

 

= 0.0000035,

поэтому 1 eRx

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 x

 

2

 

72 600

 

. Итак, подтверждаются все пять цифр после запятой:

 

 

 

 

 

 

.

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3600

 

 

30240

 

75600

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

282240

 

 

 

 

 

R0

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

(0.5599303; 0.5599373)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.55993

 

 

 

 

 

 

Приближенное решение дифференциальных уравнений.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример. Найти первые четыре члена разложения в ряд Маклорена решения задачи Коши

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y00 = yy0e−x;

y(0) = 1;

 

 

y0(0) = 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1-й способ. Подставляя в уравнение x = 0, находим y00(0):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y00(0) = y(0)y0(0)e0 = 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Продифференцируем уравнение: y000 = (y0)2e−x + yy00e−x − yy0e−x и подставим x = 0:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y000(0) = 1 + 1 − 1 = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Еще раз продифференцируем уравнение:

 

 

(4)

 

 

 

 

0

 

00

 

−x

0

2

 

 

−x

0

 

00

 

−x

 

 

 

 

000

 

−x

−yy

00

 

−x

−(y

0

2

 

−x

−yy

00

 

−x

 

yy0e−x и подставим x = 0:

 

 

 

 

 

 

 

 

 

y

 

 

= 2y

y

 

 

e

 

 

(y

)

e

 

 

 

 

+ y

y

 

e

 

 

+ yy

 

 

e

 

 

 

 

 

e

 

 

)

e

 

 

 

e

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(4) = 2 − 1 + 1 + 1 − 1 − 1 − 1 + 1 = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y0(0)

 

 

 

 

y00

(0)

 

2

 

y000(0)

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку разложение в ряд Маклорена имеет вид y(x) = y(0) +

 

 

 

x +

 

 

 

 

x

 

+

 

 

 

 

 

x

 

+ . . ., получаем:

 

 

 

 

 

 

1!

 

2!

 

 

 

 

 

3!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x) = 1 + x +

x2

+

x3

 

+ . . .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2-й способ. Пусть

 

 

 

 

 

 

 

2

 

 

 

3

 

 

 

 

; тогда

y

0

= c1

+ 2c2x + 3c3x

2

 

 

 

 

;

 

00

= 2c2

+ 3

 

 

 

 

 

 

. Напишем

 

 

 

 

y = c0 + c1x + c2x + c3x +2. . .

x

3

 

 

 

 

 

 

 

+ . . . y

 

 

· 2x + . . .

 

 

 

 

также разложение функции e−x

:

e−x = 1 − x +

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ . . .. Начальные условия дают нам первые два коэффициента

 

 

2

6

 

разложения: c0 = 1; c1 = 1. Далее, выписанные разложения подставляем в уравнение и приравниваем коэффициенты

при одинаковых степенях x, начиная с младших. Константа в левой части равна 2c2; в правой части равна 1 · 1 · 1 = 1; 2c2 = 1; c2 = 1/2. Коэффициенты при первой степени x в левой и правой частях уравнения равны 6c3 è 1 + 2c2 − 1 = 1

соответственно, откуда c3 = 1/6. Получаются те же четыре члена разложения, которые уже были получены первым способом.

Замечание. Можно обратить внимание на то, что полученные члены разложения функции y(x) совпадают с членами

разложения функции ex. Непосредственная подстановка в уравнение подтверждает, что ex действительно является решением задачи Коши.

30

Соседние файлы в папке Математический анализ. Лекции.