Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
2ой семестр / Математика / Теория / Магазинников Л.И. - Функции комплексного переменного.pdf
Скачиваний:
92
Добавлен:
18.07.2023
Размер:
1.21 Mб
Скачать
dxn−2

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

b

R

Если подынтегральная функция в интеграле f(x)dx разлагает-

a

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

1

Пример 3.29. Вычислить I = R cos x2dx с точностью α = 0,001.

0

Решение. Разлагая функцию в ряд Тейлора по степеням x, по-

лучаем ряд cos x2 = 1 −

x4

x8

x12

 

+

 

 

+ · · ·, сходящийся на всей

2!

4!

6!

оси. Интегрируя почленно этот ряд в пределах от 0 до 1, находим

1

 

 

1

 

 

1

 

 

1

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

5

·

2!

+ 9

·

4!

13

·

6!

 

 

cos x2dx = 1 −

 

 

 

+ · · ·. Мы получили знакоче-

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

редующийся ряд, четвёртое слагаемое его

 

< 0,001. Поэтому

13 · 720

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R cos x2dx = 1 − 0,1 + 0,00463 ≈ 0,905.

0

3.5.4. Интегрирование дифференциальных уравнений

Требуется найти решение y = y(x) дифференциального уравне-

ния

y′′ = F (x, y, y),

(3.59)

 

удовлетворяющее начальным условиям y(x0) = a0, y(x0) = a1. Бу-

дем считать, что в окрестности начальных данных уравнение (3.59) удовлетворяет условиям теоремы существования и единственности, а также, что решение y = y(x) представимо в виде суммы ряда Тей-

лора

y(x) = y(x0) + y(x0)(x − x0) + · · · + y(n)(x0)(x − x0)n + · · · . (3.60) n!

Чтобы найти ряд (3.60) нужно вычислить его коэффициенты y(x0), y(x0), y′′(x0), . . . , y(n)(x0), . . . , причём первые два определены начальными условиями y(x0) = a0, y(x0) = a1. Укажем два метода

отыскания этих коэффициентов. Значения всех производных функции y(x) в точке x0 можно найти последовательным дифференцированием по x функции F (x, y(x), y(x)):

y′′′(x0) = dF (x0, a0, a1) , · · · , y(n)(x0) = dn−2F (x0, a0, a1) , · · ·. dx

82

Пример 3.30. Найти пять членов разложением в ряд Тейлора решения y(x) дифференциального уравнения y′′ = xyy, удовлетворя-

ющего условию y(0) = 1, y(0) = 1.

Решение. Имеем y′′ = xyy, y′′(0) = 0; y′′′ = yy+ x(y)2 + xyy′′, y′′′(0) = 1; y(IV ) = 2(y)2 + 2yy′′ + 3xyy′′ + xyy′′′, y(IV )(0) = 2; y(V ) = 9yy′′ + 3yy′′′ + 3x(y′′)2 + 4xyy′′′ + xyy(IV ), y(V )(0) = 3. Сле-

довательно, y(x) = 1 + x +

x3

2x4

 

3x5

 

+

 

 

+

 

 

+ · · ·.

3!

4!

5!

Метод неопределённых коэффициентов заключается в следующем. Ищем решение y(x) в виде ряда

 

y(x) = X an(x − x0)n

(3.61)

n=0

с неопределёнными коэффициентами an, которые требуется найти.

Дифференцируя дважды ряд (3.61) почленно, получаем

X

X

y= nan(x − x0)n−1, y′′ =

n(n − 1)an(x − x0)n−2. (3.62)

n=1

n=2

Затем разлагаем в ряд Тейлора функцию F (x, y, y), используя при

этом ряды (3.62). Пусть

 

 

 

 

 

X

 

F (x, y, y) =

bn(x − x0)n.

(3.63)

 

 

n=0

 

Вносим (3.62) и (3.63)

в

исходное уравнение.

Получим

 

 

X

 

X

 

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

=

bn(x − x0)n. Приравнивая ко-

n=2

 

n=0

 

эффициенты при одинаковых степенях (x − x0), мы и получим

необходимые соотношения для определения неизвестных коэффициентов an, причём коэффициенты a0 и a1 находим из начальных

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

Пример 3.31. Найти общее решение уравнения y′′ = y методом

неопределённых коэффициентов.

X

Решение. Ищем решение в виде ряда y(x) =

anxn, тогда

 

n=0

 

 

 

X

X

 

y′′ =

n(n − 1)anxn−2, следовательно,

n(n − 1)anxn−2 =

n=2

n=2

 

 

 

 

X

=anxn. Приравнивая коэффициенты при одинаковых степе-

n=0

нях x, получаем n(n − 1)an = an−2, n = 2, 3, . . .. Отсюда находим

83

 

a0

 

 

 

a1

 

 

 

 

a0

 

 

 

 

 

 

a1

a2 =

 

, a3 =

 

 

,

a4 =

 

 

 

, a5 =

 

 

, · · ·,

1 · 2

1 · 2 · 3

1 · 2 · 3 · 4

1 · 2 · 3 · 4 · 5

a2n =

a0

 

, a2n−1

=

 

a1

 

. Таким образом, решение можно за-

 

 

 

 

 

 

 

(2n)!

 

 

(2n − 1)!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2n

 

x2n−1

 

 

 

 

 

 

 

 

 

X

 

 

 

 

X

 

 

 

 

 

 

 

писать в виде y(x) = a0

(2n)!

 

+ a1

(2n

 

1)!

 

= a0ch x + a1sh x,

 

 

 

 

 

 

n=0

 

 

 

 

n=1

 

 

 

 

 

 

 

 

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

3.5.5. Применение рядов Тейлора к отысканию пределов и производных

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x cos x −

 

sin

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 3.32. Найти

lim

x

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x→+0

1 − cos x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Разложим в ряд Тейлора числитель и знаменатель дан-

ной дроби:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x3

 

x5

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

x3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x cos x − x sin x =

x −

 

 

 

+

 

 

− · · · − x −

 

 

+

 

 

 

 

− · · · ;

2!

 

4!

3!

5!

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

x4

 

 

 

 

 

 

 

 

 

x cos x −

 

sin

 

 

=

1

cos x =

 

 

+

· · ·

. Поэтому

 

lim

x

x

 

2!

 

 

4!

 

 

 

 

 

 

 

cos x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

+0

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

x2

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

 

 

 

 

 

+

 

 

 

x3 + · · ·

=

 

 

 

 

 

+

 

 

 

x + · · ·

=

1

.

=

3!

2!

5!

 

lim

 

6

2

5!

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

x→+0

 

x4

 

 

 

 

 

 

 

x→+0

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

+ · · ·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ · · ·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2!

4!

 

 

 

 

 

 

 

 

 

 

 

 

2

 

4!

 

 

 

 

 

 

 

 

 

 

Пример 3.33. Найти f(V )(0), если f(x) = e2x sin 3x.

Решение. Ряд Тейлора для f(x) можно получить как произведение рядов для функций e2x и sin 3x:

f(x) = 1 + 2x +

22x2

23x3

 

 

 

24x4

 

 

 

 

 

3

3x3

35x5

 

 

 

 

 

+

 

 

 

+

 

 

 

+ · · · 3x −

 

 

 

+

 

 

− · · · .

2!

3!

 

4!

 

 

3!

 

5!

 

Коэффициент a

 

 

при x5

равен

35

 

22 ·

33

+

24 · 3

=

81

 

 

9+2 =

 

199

.

 

 

 

 

 

 

 

 

 

 

5

 

 

 

2!

 

40

 

 

 

 

 

 

 

 

 

5!

3!

·

4!

 

 

 

 

 

 

40

 

 

f(V )(0)

 

 

 

 

199

 

 

 

 

 

 

−199 · 120

 

 

 

 

 

 

Следовательно,

 

=

 

, т.е. f(V )(0) =

=

597.

 

 

5!

 

 

 

 

 

 

 

 

 

 

 

40

 

 

 

 

 

 

 

 

 

 

40

 

 

 

 

 

 

84