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

Методическое пособие 388

.pdf
Скачиваний:
4
Добавлен:
30.04.2022
Размер:
1.11 Mб
Скачать

Следовательно, интервалом сходимости является интервал

2 x 2 2 или 4 x 0 .

Теперь надо исследовать поведение ряда на концах интервала.

При x 4 имеем ряд

 

2

n

2

 

 

n 1

, который сходится по при-

 

 

 

( 1)

 

n

знаку Лейбница.

n 1 n2n 1

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2n

 

 

 

 

При x 0 получаем расходящийся ряд

 

 

2

 

1

− это гармониче-

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

n 1

n2

 

n 1n

 

ский ряд.

Следовательно, областью сходимости является полуинтервал [ 4,0) .

2.2. Свойства степенных рядов

В дальнейшем мы будем пользоваться основными свойствами степенных рядов, которые приведём без доказательства.

1. Сумма S(x) степенного ряда (9.4) является непрерывной функцией в интервале сходимости ( R , R ).

2. Степенные ряды

 

 

xn è

 

 

a

n

b xn , имеющие радиусы сходимости со-

 

n 1

 

n 1

n

 

 

 

 

ответственно R1 и R 2 , можно почленно складывать, вычитать и умножать. Ра-

диус сходимости суммы, разности и произведения этих рядов не меньше, чем меньшее из чисел R1 и R 2 .

3. Степенной ряд внутри интервала сходимости можно почленно дифференцировать; при этом для ряда

S(x) a

0

a x a

2

x2

... a

n

xn ...

(9.8)

 

1

 

 

 

 

при R x R выполняется равенство

 

2a2 x 3a3 x

2

... n an x

n 1

... .

(9.9)

S (x) a1

 

 

4. Степенной ряд можно почленно интегрировать на каждом отрезке, расположенном внутри интервала сходимости; при этом для ряда (9.4) приR x R выполняется равенство

x

x

 

x

x

 

t n dt ... .

(9.10)

S(x)dx a

0

dx a tdt ... a

n

a

a

1

a

 

 

 

a

 

 

 

Ряды (9.9) и (9.10) имеют тот же радиус сходимости, что и исходный степенной ряд.

Эти свойства широко используются в теоретических исследованиях и в приближенных вычислениях.

131

§ 3. Разложение функций в степенные ряды 3.1. Разложение в степенные ряды элементарных функций

Рассмотрим ряд xn 1 x x2 x3 ... xn . , который представляет

n 1

собой геометрическую прогрессию со знаменателем q x . Известно, что он

сходится при

 

x

 

1, а его сумма равна S

 

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

x

 

 

 

 

 

Таким образом, можно записать

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1 x x2 ... xn

... .

(9.11)

 

 

 

 

1 x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это равенство означает, что функция f (x)

 

 

1

представлена в виде

1

x

 

 

 

 

 

 

 

 

 

 

 

 

 

степенного ряда (говорят, что функция разложена в степенной ряд), радиус схо-

димости которого R 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставив в равенстве (9.11) вместо x t , получим разложение в сте-

пенной ряд функции f (t)

 

 

 

1

 

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1 t t 2 t3

... ( 1)n t n ...,

(9.12)

 

 

 

 

 

1 t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для которого также R 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1, получаем

 

 

Интегрируя почленно ряд (9.12) по отрезку [0, x], где

 

x

 

 

 

 

 

 

x

dt

x [1 t t2 t3 ... ( 1)n tn ...]dt или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 1 t

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln(1 x) x

x2

 

 

x3

 

x4

... ( 1)

n xn 1

... при

 

x

 

1.

(9.13)

 

 

 

 

 

 

 

2

 

 

3

 

4

 

(n 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При x 1 ряд также сходится, а при x 1 ряд расходится. Следовательно, областью сходимости ряда (9.13) является полуинтервал ( 1,1].

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

Для приложений важно уметь всякую функцию f (x) разлагать в степен-

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

Предположим, что в интервале (x0 R , x0 R ) функция f (x) разлагается в степенной ряд по степеням (x x0 ) , т. е.

f (x) a0 a1(x x0 ) a2 (x x0 )2 a3 (x x0 )3 ... an (x x0 )n .... (9.14)

132

Согласно свойствам степенных рядов, его можно почленнодифференцировать:

 

 

a1 2a2 (x

x0 ) 3a3 (x x0 )

2

 

nan (x x0 )

n 1

,

 

 

 

 

 

f (x)

 

 

 

 

 

 

 

 

 

f

 

1 2 a2 1 2 3a3 (x x0 )

 

 

n(n 1)(x x0 )

n 2

 

,

 

 

 

 

 

(x)

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

n(n 1)(n 2)an (x x0 )

n 3

 

,

(x) 1 2 3a3 1 2 3 4a4 (x x0 )

 

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

 

 

 

 

 

 

f (n) (x) 1 2 3 ...

(n 1)na

1 2

3...

 

n(n 1)(x x ) ...

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

Полагая в этих равенствах и в равенстве (9.14) x x0 , получаем

 

 

 

 

1

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

(n)

 

 

 

 

a0 f

(x0 ) , a1 1!

2! f

 

 

 

3! f

 

 

 

n!

f

 

(x0 ).

 

 

f (x0 ), a2

(x0 ),a3

(x0 ),an

 

 

 

Подставляя эти равенства в (9.14), получаем ряд

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x) f (x

)

f (x0 )

(x x )

f (x0 )

(x x )2

...

 

 

f (n) (x0 )

(x x

 

)n

...

,

 

(9.15)

 

 

 

 

 

 

 

0

1!

 

0

 

2!

 

 

 

0

 

 

 

 

n!

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

который называется рядом Тейлора функции f (x) . Его коэффициенты назы-

ваются коэффициентами Тейлора.

Если в ряде Тейлора положить x0 0 , то получим частный случай ряда Тейлора, который называется рядом Маклорена:

 

 

f

 

 

 

 

2

 

f

(n)

(0)

 

n

 

f (x)

f (0)

(0)

x

f (0)

x

...

 

x

... .

 

1!

2!

 

 

n!

 

 

 

 

 

 

 

 

 

 

 

 

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

Но отсюда еще не следует, что он будет сходиться к данной функции; он может сходиться, но не к функции f (x). Так, например, функция

 

 

1

 

 

 

 

x

2

,

при

x 0,

f (x) e

 

 

 

 

,

 

при

x 0

0

 

имеет в точке x 0 производные всех порядков, причем f (n) (0) 0 при всех n . Ряд Маклорена имеет вид

0 10! x 20! x2 ... n0! xn ... .

Он сходится, но его сумма S(x) в любой точке равна нулю, а не f (x). Пусть функция f (x) в интервале (x0 R, x0 R) имеет производные любо-

гопорядка.

133

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

f (x) f (x

0

)

f (x0 )

(x x

0

) ...

f (n) (x0 )

(x x

0

)n r (x) ,

 

 

 

 

1!

 

 

n!

n

 

 

 

 

 

 

 

 

где

rn (x)

f (n 1) [x0 (x x0 )]

(x x0 )n 1 , который называется остаточным чле-

 

 

 

(n 1)!

 

ном ряда в форме Тейлора.

 

 

Из равенства f (x) Sn (x) rn (x) следует, что ряд Тейлора сходится то-

гда и только тогда, когда lim

r (x) 0.

 

 

n

n

 

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

ся по существу к определению значений x , при которых rn (x) 0.

Если сделать это трудно, то следует каким-нибудь способом убедиться, что написанный ряд Тейлора сходится к данной функции.

Теорема. Если в окрестности точки x0 все производные функции f (x) ограничены одним и тем же числом M 0 , т. е. f (n) (x) M , то для всех x из

этой окрестности ряд Тейлора для этой функции сходится, и его сумма равна f (x). (без доказательства).

4.1. Разложение некоторых элементарных функций в ряд Тейлора (Маклорена)

Для разложения функции f (x) в ряд Маклорена (9.15) нужно:

1)найти производные f (x) , f (x) , ..., f (n) (x),...;

2)вычислить значения производных в точке x0 0;

3)написать ряд (9.16) для данной функции и найти интервал сходимости;

4)найти интервал ( R, R ) , в котором остаточный член ряда Маклорена

rn (x) 0 при n .

Если такой интервал существует, то в нем функция f (x) и сумма ряда

совпадают.

Получим разложения в ряд Маклорена основных элементарных функций. 1. Разложение в ряд Маклорена функции f (x) e x .

Получаем

1)f (x) ex , f (x) ex ,..., f (n) (x) ex ,... ;

2)f (0) 1, f (0) 1,... f (n) (0) 1,...;

3)

ex ~1

 

x

 

x2

...

xn

...;

(9.16)

 

 

 

2!

 

 

1!

 

 

n!

 

 

 

 

an

 

 

 

(n 1)!

 

lim n 1 , т. е. ( , );

4)

R lim

 

 

 

lim

 

an 1

 

n!

 

 

n

 

n

 

 

 

n

 

 

 

 

 

 

 

 

 

134

5) имеем для всех x ( R, R) , f (n) (x) ex eR M , то есть все про-

изводные в этом интервале ограничены одним и тем же числом M eR . Следовательно, по рассмотренной выше теореме nlim rn 0.Таким образом, в разло-

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

 

 

 

 

 

 

 

 

 

 

е

x

1

x

 

 

x2

...

xn

... .

 

 

 

 

 

 

 

 

 

 

 

 

1!

 

2!

n!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Разложение в ряд функции f (x) sin x.

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1)

(x) cos x sin(x

2 ) ,

f (x) sin x sin(x 2

2 ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x) sin(x

2 ),...

 

 

 

 

 

 

 

 

 

 

 

……………………………………..

 

 

 

 

 

 

 

……………………………………..

 

 

 

 

 

 

 

f n (x) sin(x n )

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

………………………………………

 

 

 

 

 

 

 

 

0

 

,

n 0,2,4,...

 

 

 

 

 

 

 

 

 

 

2)

 

 

 

 

 

 

 

n 3,7,11,,...

 

 

 

 

 

 

 

f (n) (0) 1 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1,5,9,...

 

 

 

 

 

 

 

 

 

 

 

1 ,

 

 

 

 

 

 

 

 

 

 

3)

sin x ~ x

 

x3

 

 

 

x5

... ( 1)

n

x

2n 1

 

... .

 

 

3!

 

5!

 

(2n 1)!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Нетрудно проверить, что R и ряд сходится на всей числовой оси, т. е. при всех x ( , ).

4) любая производная функции f (x) sin x по модулю не превосходит

 

 

f (n) (x)

 

 

)

 

1 .

 

 

 

единицы, т. е.

 

 

 

sin(x n

 

Следовательно, согласно теореме,

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

имеет место разложение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

2n 1

 

 

 

 

 

 

sin x ( 1)n

 

 

.

 

 

 

 

 

(2n 1)!

 

 

 

 

 

 

n 1

 

 

3. Разложение в ряд функции f (x) cos x .

Формулу разложения этой функции в ряд Маклорена можно получить аналогично. Однако проще получить разложение функции cos x , воспользовавшись свойством 3 степенных рядов.

135

Продифференцировав предыдущий ряд почленно ряд, получим

cos x 1

x

2

 

x

4

... ( 1)n

x

2n

...

 

( 1)n

x

 

.

 

 

 

 

 

 

 

 

2n!

2n!

2!

4!

 

 

n 1

 

 

Этот ряд также сходится к данной функции на интервале ( , ) .

1. Разложениеврядфункции f (x) (1 x)a , где a любоевещественноечисло. Имеем

1)f (x) a(1 x)a 1 , f (x) a(a 1)(1 x)a 1,... f (n) (x)

a(a 1)...(a n 1)(1 x)a n ,.n N.

2)

f (0) 1, f

 

a,... f

(n)

(0)

a(a

1)...(a n 1),...

 

 

 

 

(0)

 

 

 

 

3)

(1 x)

a

1

ax

a(a 1)

x

2

...

a(a 1)(a 2)...(a n 1)

x

n

... .

 

2!

 

 

n!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Можно доказать, что этот ряд сходится к данной функции при x 1.

Этот ряд называется биномиальным. Если a n N , то все члены ряда с (n 1) го номераранынулю, таккаксодержатмножитель (a n) n n 0.

ВэтомслучаеполучаемизвестнуюформулубиномаНьютона:

(1 x)

n

1

 

n

x

n(n 1)

x

2

...

n(n 1) (n k 1)

x

n k

... x

n

. (9.17)

 

1!

2!

 

k!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим некоторые частные случаи биномиального ряда.

5.При a 1:

1( 1)n xn . 1 x n 1

 

 

1

 

 

 

 

 

 

 

1

 

 

n 1 1 3 5 ... (2n 3)

n

6.

При a

2

:

 

1 x 1

 

2

x ( 1)

 

 

n

 

x .

 

 

 

 

 

 

 

 

 

n 1

 

 

2 n!

 

 

 

 

 

1

 

1

 

 

 

 

 

 

n

 

При a

 

 

 

 

 

n 1 3 5 ...(2n 1)

 

7.

2

:

 

 

 

 

1 ( 1)

 

 

n

 

x .

1

 

x

 

2

n!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

8.

Разложение в ряд функции f (x) ln(1 x) .

 

 

Разложение этой функции получим путем интегрирования ряда (9.17) в промежутке от 0 до x при x 1:

ln(1 x)

 

( 1)n

x

n 1

, ( 1 x 1).

 

 

 

n 1

 

n 1

 

 

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

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

136

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

Пусть требуется вычислить значение функции f (x) при x x1 с заданной

точностью 0.

Если функцию f (x) в интервале ( R, R) можно разложить в степенной

ряд f (x) a

0

a

x a

2

x2 ...a

n

xn ... , то точное значение

f

(x ) при x (

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

равно сумме этого ряда при x x , т. е. f (x

1

) a

0

ax

1

a

2

x2

...a

n

xn

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

а приближенное − частичной сумме Sn (x1) , т. е.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x ) S

n

(x ) a

0

a x a

2

x21 ... a

n

xn1.

 

 

 

 

 

 

 

 

1

 

1

 

1 1

 

 

 

 

 

 

 

 

 

 

 

R, R)

...,

Точность этого равенства увеличивается с увеличением числа членов. Абсолютная погрешность этого приближенного равенства равна модулю

остатка ряда, т. е.

 

f (x ) S

n

(x

)

 

 

 

r (x

)

 

, где r (x ) a

xn 1

an 2

... .

 

 

 

 

 

1

 

 

1

 

 

 

 

n 1

 

 

n

1

n 1 1

n 2

 

Итак, ошибку

 

f (x1) Sn (x1 )

 

можно найти, оценив остаток rn (x1) ряда.

 

 

Для знакочередующихся рядов, согласно признаку Лейбница, имеем rn (x1 ) un 1 (x1 ) un 2 (x1 ) ... un 1 (x1 ) .

Пример 1. Вычислить 3 124 с точностью до 0,0001.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

Решение: преобразуем корень 3 124 5(1

 

 

 

 

 

)

 

3 .

 

 

 

 

 

 

124

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Применяя формулу (9.17) при n

1 , получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 124 5[1 (

1)

1

 

 

( 3)(

3

 

1)

(

1

 

)2

...]

 

 

1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

3 124

 

 

 

 

 

 

 

 

124

 

 

 

 

 

 

 

5(1

 

1

 

 

2

...) 5

 

 

5

 

 

 

 

 

5

 

 

....

 

3

124

2 2

372

9

62 124

 

 

 

3 124

 

 

 

 

 

 

 

 

 

Поскольку это знакочередующийся ряд и

 

 

 

 

 

 

5

 

 

 

 

 

 

 

5

0,0001,

9 124 124

69192

 

 

 

 

 

 

 

 

 

 

 

 

 

 

то 3 124 5 3725 5 0,0134 4,9866.

Пример 2. Вычислить число e с точностью до 0,001. Решение. Подставляя x 1 в формулу (9.16) получим

e 1 1!1 2!1 ... n1! ...

137

Возьмём n слагаемых и оценим ошибку rn (x) :

 

 

 

 

 

 

 

 

r (x)

 

 

 

1

 

 

 

 

 

 

1

 

 

...

1

 

1

 

1

 

 

 

 

1

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

(n

1)!

 

(n 2)!

 

 

 

(n

 

 

n

2 (n 2)(n 3)

 

 

 

 

 

 

 

 

1)!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

1

,

те. .

rn (x)

1

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(n 1)!

1

 

 

1

 

 

n! n

n!n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Остается подобрать наименьшее натуральное число n , чтобы выполнялось неравенство n1!n 0.001. Нетрудно проверить, что это неравенство выпол-

няется при n 6. Поэтому имеем

e 1 2!1 3!1 4!1 5!1 6!1 2,718.

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

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

b

Пусть, например, требуется вычислить f (x)dx с точностью до 0.

a

Если подынтегральную функцию f (x) можно разложить в ряд Тейлора

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

1

Пример. Вычислить определенный интеграл 3 e x2 dx с точностью до 0,001.

0

Решение. Применить для вычисления этого интеграла формулу НьютонаЛейбница мы не можем, так как первообразная для e x2 не выражается в элементарных функциях ( e x2 dx неберущийся интеграл).

Поэтому разложим подынтегральную функцию в ряд Маклорена, заменяя

в формуле (9.16) x на ( x2 ) :

 

 

 

 

 

 

e x2

1

x2

 

x4

 

x6

...,

x ( , ).

 

 

 

 

1!

2!

3!

 

 

138

Интегрируем этот ряд на отрезке [0,

1] :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

1

 

 

1

 

x 2

 

x 4

 

 

 

 

1

1

 

 

1

 

 

 

3 e x2 dx

3 (1

 

 

...) dx

 

 

... .

 

 

1!

2!

 

 

3 1! 33

5

2! 35

 

 

0

 

 

0

 

 

 

 

 

 

3

 

 

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

так как

 

1

 

0,001 ,а

1

 

 

 

0,001, то с точностью до 0,001 имеем

5

 

5

 

 

 

 

3

 

 

2! 3

 

 

 

3 1! 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 e x 2 dx

1

 

1

0, 3210.

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

3

 

3 3

 

 

 

 

 

 

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

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

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

5.3.1. Способ последовательного дифференцирования

Пусть требуется решить уравнение

y

 

f (x, y, y ),

(9.18)

 

 

 

удовлетворяющее начальным условиям

 

 

 

y

 

x x0

y0 ,

y

 

x x0

y0

 

.

 

 

(9.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение y y(x) уравнения (9.18) ищем в виде ряда Тейлора:

 

 

y y(x

)

y (x0 )

(x x )

y (x0 )

(x x )2

...

 

 

y(n) (x )

(x x

) ...

(9.20)

 

 

 

 

0

 

1!

0

2!

0

 

 

 

n!

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Первые два коэффициента находим из начальных условий (9.19). Подставив в уравнение (9.18) значения x x0 , y y0 , y y0 , находим

третий коэффициент: y (x0 ) f (x0 ; y0 ; y0 ) . Значения y (x0 ), y(n) (x0 ),... нахо-

дим путем последовательного дифференцирования уравнения (9.18) по x и вычисления производных при x x0 . Найденные значения производных (коэффи-

циентов) подставляем в равенство (9.20).

Пример. Решить уравнение y

 

x

2

y

2

, y( 1) 2 ,

 

 

1

; найти

 

 

 

y ( 1)

2

пять членов (отличных от нуля) разложения в ряд Тейлора.

139

Решение. Будем искать решение уравнения в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y y( 1)

y ( 1)

(x 1)

y

( 1)

(x 1)2

 

y ( 1)

(x 1)3

...

 

 

2!

 

1!

 

 

1

 

3!

 

 

 

 

 

 

 

Здесь y( 1) 2 ,

 

. Находим

 

( 1)

2

 

2

2

5. Для нахож-

y ( 1)

2

y ( 1)

 

 

дения последующих коэффициентов дифференцируем данное дифференциальное уравнение:

y 2x 2 yy ,

y(4) 2 2( y )2 2yy ,

y(5) 4y y 2y y 2yy 6y y 2yy ,....

При x 1 имеем

y ( 1) 2 2 2 12 0,

y(4) ( 1) 2 2 14 2 2 5 22,5,

y(5) ( 1) 6

1

5 2 2 0 15 , ...

 

2

 

Подставляя найденные коэффициенты в искомый ряд, получим y 2 12 (x 1) 52 (x 1)2 1615 (x 1)4 18 (x 1)5 ....

ПРИЛОЖЕНИЕ

§ 6. Способ неопределенных коэффициентов решения дифференциальных уравнений

Пусть, например, требуется решить задачу Коши:

y p (x) y q(x) y f (x) ; y(x0 ) y0 , y (x0 ) y0.

Предполагая, что коэффициенты p(x), q(x) и свободный член f (x) разлагаются в ряды по степеням (x x0 ) , сходящимися в некотором интервале (x0 R, x0 R) , искомое решение y y(x) ищем в виде степенного ряда с неопределенными коэффициентами:

y c

c

(x x

) c

(x x

)2 ... c

(x x

)n ...

(9.22)

0

1

0

2

0

n

0

 

 

Коэффициенты c0 y0, c1 y0 определяютсяприпомощиначальныхусловий.

Для нахождения последующих коэффициентов дифференцируем ряд (9.22) два раза (каков порядок ДУ) и подставляем выражения для функции y и

ее производных в уравнение (9.21), заменив в нем p(x), q(x), f (x) их разложе-

140