Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Лосев-КлячинМиклюков МА в КИ

.pdf
Скачиваний:
863
Добавлен:
27.05.2015
Размер:
3.66 Mб
Скачать

§9. Приближенное вычисление корней уравнений

131

ТЕОРЕМА 9.1. Предположим, что функция F дифференцируема всюду на (−∞, +) и удовлетворяет там условию (1). Тогда

i)существует и единственен корень ξ уравнения (2);

ii)итерационная последовательность {xn}, описываемая равенством (3), сходится к ξ при любом выборе начальной точки x0;

iii)для любого n > 1 выполнено

x

n

ξ

| ≤

qn

|x1 − x0

|

.

(6)

|

 

1

q

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство. Прежде всего заметим, что итерационная последовательность {xn} фундаментальна. Действительно, фиксируем произвольно ε > 0. Так как q < 1, то существует номер N(ε) такой, что при всех n > N(ε) выполнено

qn |x1 − x0| < ε.

1 − q

Таким образом, для любых m > n > N(ε), в силу (5), имеем

|xm − xn| ≤ qn |x1 − x0| < ε,

1 − q

что и требуется.

Поскольку последовательность {xn} фундаментальна, то на основании критерия Коши она сходится. Обозначим ее предел через ξ. Так как F дифференцируема, то она непрерывна, и потому из (3) вытекает

ξ = lim xn+1 = lim F (xn) = F (ξ),

n→∞ n→∞

т.е. число ξ является корнем уравнения (2).

Докажем, что этот корень является единственным. Предположим противное, т.е. для некоторого η = ξ также выполнено F (η) = η. Не умаляя общности, можем считать, что ξ < η. Согласно формуле конечных приращений Лагранжа,

имеем

0 < η − ξ = F (η) − F (ξ) = F (c)(η − ξ),

где c (ξ, η) – некоторая точка. Отсюда, в силу условия (1), находим

|η − ξ| ≤ |F (c)||η − ξ| ≤ q|η − ξ|.

132

Глава 5. Основные теоремы дифференциального исчисления

Противоречие с предположением 0 < q < 1.

Имеются модификации метода, пригодные для вычисления корня уравнения (2) в случае, когда F задана на конечном отрезке [a, b]. См., например, В.А. Ильин и Э.Г. Позняк, Основы математического анализа, часть I, стр. 392.

9.2.Метод касательных

Предположим, что уравнение F (x) = 0 имеет корень ξ на отрезке [a, b]. Воспользуемся следующим алгоритмом для приближенного нахождения корня.

Пусть: x0 – произвольная начальная точка;

x1 – точка пересечения касательной к графику F в точке x0 с осью ox;

x2 – точка пересечения касательной к графику F в точке x1 с осью ox;

...............................................................................

xn – точка пересечения касательной к графику F в точке xn−1 с осью ox.

§9. Приближенное вычисление корней уравнений

133

Найдем итерационную формулу. Уравнение касательной, проходящей через точку xn−1, имеет вид

y = F (xn−1)(x − xn−1) + F (xn−1).

Тем самым, точка пересечения этой касательной с осью y = 0, т.е. точка xn определяется по формуле

F (xn−1) xn = xn−1 F (xn−1).

9.3.Метод хорд

Предположим, что уравнение F (x) = 0 имеет корень ξ на отрезке [a, b].

Описание метода: x0 [a, b] – произвольная начальная точ-

ка; x1 – точка пересечения хорды A0B с осью y = 0; x2 – точка

пересечения хорды A1B с осью y = 0 и т.д.

УПРАЖНЕНИЕ 1. Обоснования методов касательных и хорд посмотреть самостоятельно (например, в книге В.А. Ильина и Э.Г. Позняка, Основы математического анализа, часть I, стр. 393-398).

Глава 6

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

§1. Производные многочлена и его разложение по степеням

ТЕОРЕМА 1.1. Пусть P (x) = n ckxk – произволь-

k=0

ный многочлен степени n. Пусть a R – произвольная точка. Тогда

n

P (k)(a)

 

 

k

 

(x − a)k.

(1)

k!

P (x) =

=0

 

 

 

ЗАМЕЧАНИЕ. Здесь P (k)(a) – значение производной k-го порядка в точке a при k = 1, 2, ..., n. При k = 0 – значение самой функции.

Доказательство. Произведя замену x = (x − a) + a в многочлене P (x), получаем

n

P (x) = ck[(x − a) + a]k,

k=0

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

P (x) = d0 + d1(x − a) + d2(x − a)2 + ... + dn(x − a)n. (2)

Найдем коэффициенты d0, d1, ..., dn. Подставляя в обе части

(2) x = a, имеем

P (a) = d0.

Дифференцируя равенство (2) по x, приходим к соотношению

P (x) = d1 + 2d2(x−a) + ...+ kdk(x−a)k−1 + ...+ ndn(x−a)n−1

и при x = a получаем

P (a) = d1.

§2. Формула Тейлора с остаточным членом в форме Пеано

135

Далее

P (x) = 2d2 + 3 · 2d3(x − a) + ... + k(k − 1)dk(x − a)k−2 + ...+ +n(n − 1)dn(x − a)n−2

и

P (a) = 2d2.

Продолжая далее этот процесс, находим

P (k)(a) = k!dk

или

P (k)(a)

 

 

 

 

 

dk =

 

 

 

 

 

 

 

,

(k = 0, 1, 2, . . . , n).

 

 

 

 

 

 

 

k!

 

 

 

 

 

Следовательно, из (2) будем иметь

 

 

 

 

P (a)

P (a)

P (n)(a)

 

P (x) = P (a)+

 

 

(x−a)+

 

(x−a)2+. . .+

 

 

(x−a)n.

1!

 

2!

n!

§2. Формула Тейлора с остаточным членом в форме Пеано

Рассмотрим задачу о существовании многочлена

Pn(x) = c0 + c1x + c2x2 + ... + cnxn

такого, чтобы для функции f в окрестности точки x = a выполнялось соотношение

f(x) − Pn(x) = o((x − a)n) x → a.

Было дано решение задачи при n = 1 (см., например, параграф "Дифференциал" (формула (3)), или, формула "конечных приращений" Лагранжа). Приведем ее решение в общем виде.

Согласно теореме 1.1 предыдущего параграфа, многочлен Pn(x) можно записать в виде

n

(k)

 

Pn(x) =

Pn (a)

(x − a)k.

(1)

k!

k=0

136

Глава 6. Формула Тейлора

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

n

f(k)(a)

 

 

k

 

(x − a)k.

(2)

k!

Tn(x) =

=0

 

 

 

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

ТЕОРЕМА 2.1. Если f имеет в точке x = a производную порядка n, то

f(x) − Tn(x) = o((x − a)n) при x → a.

(3)

Доказательство. Положим

ϕ(x) = f(x) − Tn(x).

Функция ϕ дифференцируема n раз, причем

ϕ(a) = ϕ (a) = ... = ϕ(n)(a) = 0.

Покажем, что отсюда вытекает свойство

ϕ(x) = o((x − a)n) при x → a

или, что то же самое,

ϕ(x) = ε(x)(x − a)n,

где ε(x) 0 при x → a.

Воспользуемся методом математической индукции. Утверждение справедливо при n = 1, поскольку из равенства

ϕ(a) = ϕ (a) = 0

и формулы

ϕ(x) − ϕ(a) = ϕ (a)(x − a) + o(x − a) при x → a

следует, что ϕ(x) = o(x − a) (x → a).

Предположим, что данное высказывание верно для n = k, и покажем, что оно верно при n = k + 1. Итак, пусть

ϕ(a) = ϕ (a) = ... = ϕ(k+1)(a) = 0.

Тогда производная ϕ дифференцируема k раз в точке x = a и, по предположению индукции, для нее выполнено

ϕ (x) = ε(x)(x − a)k = o((x − a)k).

(4)

1Тейлор Брук (18.8.1685-29.12.1731) – математик и философ. Род. в Эдмонтоне (Англия). Член Лондонского королевского общества.

ξ − a < 1. x − a

§2. Формула Тейлора с остаточным членом в форме Пеано

137

Однако, по формуле Лагранжа,

ϕ(x) = ϕ(x) − ϕ(a) = ϕ (ξ)(x − a),

где ξ точка, заключенная между a и x, и ξ → a при x → a. Поэтому из (4) получаем

ϕ(x) = ϕ (ξ)(x − a) = ε(ξ)(ξ − a)k(x − a) =

= ε(ξ)

ξ

a

k

x

(x − a)k+1 = ε1(x)(x − a)k+1,

a

 

 

 

где ε1(x) 0 при x → 0. Последнее равенство является следствием того, что точка ξ находится между точками a и x, а также неравенства

ЗАМЕЧАНИЕ. Формула (3) может быть записана в виде

f(x) = Tn(x) + o((x − a)n) при x → a.

(5)

Величина o((x−a)n) имеет смысл остаточного члена. Формула (5) называется формулой Тейлора с остаточным членом в форме Пеано2. В частном случае для n = 1 мы имеем

f(x) = f(a) + f (a)(x − a) + o(x − a) при x → a,

т.е. уже известную нам формулу.

ТЕОРЕМА 2.2. Если f имеет производную порядка n в

точке x = a и существует многочлен Pn(x) = c0 + c1x + ... + cnxn, для которого

f(x) = Pn(x) + o((x − a)n) x → a,

(6)

то Pn = Tn, то есть является многочленом Тейлора с центром в точке a.

Доказательство. Согласно теореме 2.1 имеем

f(x) = Tn(x) + o((x − a)n)

2Пеано Джузеппе (27.8.1858, Кунео, – 20.4.1932, Турин), итальянский математик. С 1890 профессор Туринского университета (Италия). Занимался изучением основных понятий и фактов анализа (вопрос о возможно более широких условиях существования решения дифференциальных уравнений, определение и содержание понятия кривой) и формально-логическим обоснованием математики. Во всеобщее употребление вошла его аксиоматика натурального ряда чисел. Первым построил пример непрерывной кривой, целиком заполняющей некоторый квадрат (1890 г.).

138

Глава 6. Формула Тейлора

и поэтому из (6) получаем

Tn(x) − Pn(x) = o((x − a)n) − o((x − a)n) = o((x − a)n).

Учитывая, что многочлен Pn можно записать в виде (1), находим

n

Pn(k)(a)

 

n

f(k)(a)

 

k

 

 

(x − a)k

 

 

(x − a)k = o((x − a)n).

k!

k!

=0

k=0

 

 

 

 

 

Подставляя в это равенство x = a, получаем

f(a) = f(0)(a) = Pn(a).

Следовательно, нулевые члены в суммах в (6) взаимно уничтожаются и мы получим

 

n Pn(k)(a)

 

 

n

f(k)(a)

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

(x − a)k = o((x − a)n)

 

k!

 

 

k!

 

 

k=1

 

(x − a)k

 

 

 

 

 

 

 

=1

 

 

 

 

 

 

 

 

 

или, поделив обе части равенства на (x − a),

 

n Pn(k)(a) x a k−1

 

 

n f(k)(a) x a k−1

o x a n−1 .

k

 

 

 

)

 

 

 

 

 

 

( )

= (( ) )

k!

(

 

 

k!

 

=1

 

 

 

 

 

 

 

 

 

 

k=1

 

 

 

 

 

 

 

 

 

Полагая здесь x = a, получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (a)

=

Pn(a)

.

 

 

 

 

 

 

 

 

 

1!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1!

 

 

 

 

Продолжая этот процесс, заключаем, что

 

 

 

 

 

 

 

 

f(k)(a)

 

 

Pn(k)(a)

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

k!

 

k!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Pn(x) = Tn(x).

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

Величина Rn(x) = f(x) − Tn(x) называется остаточным членом в формуле Тейлора. Ранее мы показали, что если f имеет в точке x = a производную порядка n, то Rn(x) = o((x − a)n) при x → a. Однако зачастую данной информации об остаточном члене бывает слишком мало, поскольку никаких границ для величины o((x − a)n) при этом

(n + 1)!

§3. Остаточный член в форме Лагранжа и в форме Коши

139

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

Предположим, что f дифференцируема (n + 1) раз на отрезке [a, a + h]. Зафиксируем x (a, a + h) и рассмотрим вспомогательную функцию

ψ(y) = f(x) − f(y) f (y)(x − y) − ... − f(n)(y)(x − y)n,

1! n!

где a ≤ y ≤ x. Ясно, что ψ(a) = Rn(x) и ψ(x) = 0.

Всюду на отрезке [a, a+h] существует производная ψ , при-

чем

 

 

 

 

 

 

 

 

 

 

 

 

ψ (y) = −f (y)

f (y)

(x − y) − f (y)

f (y)

 

 

 

 

 

 

(x − y)2

1!

2!

 

 

−f (y)(x−y) −· · ·−

f(n+1) y

)

 

 

f(n) y)

 

 

=

 

(

(x − y)n

 

 

(

 

(x − y)n−1

 

n!

 

 

(n

1)!

 

 

 

 

 

 

 

 

 

 

= f(n+1)(y)(x − y)n. n!

Выберем произвольно функцию g дифференцируемую на [a, x], для которой g = 0. К паре функций ψ и g применим обобщенную формулу конечных приращений Коши. Имеем

ψ(x) − ψ(a)

=

ψ (ξ)

,

g(x) − g(a)

g (ξ)

 

 

где ξ (a, x) – некоторая точка. Эту точку ξ можно записать в виде ξ = a + θ(x − a), где θ (0, 1) – некоторое число. Так

как ψ(x) = 0, ψ(a) = Rn(x) и ψ (ξ) =

f

(n+1)

(ξ)

(x − ξ)n, то

 

 

n!

 

 

мы получаем

 

 

 

 

 

 

 

 

 

 

 

 

R

(x) =

f(n+1)(ξ)

(x

ξ)n

g(x) − g(a)

.

(1)

n!

 

 

n

 

 

 

g (ξ)

 

 

Выбирая в (1) различные функции g, будем получать разные формы остаточного члена в формуле Тейлора.

Положим g(y) = x − y. Получим

Rn(x) = f(n+1)(ξ)(1 − θ)n(x − a)n+1 n!

остаточный член в форме Коши. Положим g(y) = (x − y)n+1. Тогда

Rn(x) = f(n+1)(ξ)(x − a)n+1

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

140

Глава 6. Формула Тейлора

§4. Примеры разложения функций по формуле Маклорена. Использование в приближенных вычислениях

Разложение вида

 

 

 

 

 

 

 

 

 

f (0)

 

f (0)

2

+ ... +

f(n)(0)

n

+ Rn(x),

f(x) = f(0) +

 

x +

 

x

 

 

x

 

 

2!

 

 

 

1!

 

 

 

 

n!

 

 

т.е. разложение функции по формуле Тейлора в окрестности точки a = 0, называется разложением (формулой) Маклорена3.

ПРИМЕР 1. Представим функцию ex по формуле Маклорена. Прежде всего заметим, что эта функция дифференци-

руема бесконечное число раз, причем (ex)(n) = ex и поэтому

ex = n xk + Rn(x). k!

k=0

ПРИМЕР 2. Пусть y = sin x. Тогда y(n)(x) = sin(x + nπ2 ) и y(n)(0) = sin(nπ2 ). Следовательно,

sin x = x −

x3

x5

n−1 xn

 

+

 

− · · · + (1)

2

 

+ Rn(x),

3!

5!

n!

где n – нечетно.

ПРИМЕР 3. Пусть y = cos x. Как и выше имеем

 

x2

x4

n xn

cos x = 1

 

+

 

− · · · + (1)

2

 

+ Rn(x),

2!

4!

n!

где n – четно.

ПРИМЕР 4. Мы определили число e как предел

e = lim n 1 .

n→∞ k! k=0

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

3Маклорен Колин (1698-14.6.1746). Род. в Килмодан (Шотландия). Член Лондонского королевского общества (1719). Первым опубликовал работу о разложении функций в степенные ряды.