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

2085

.pdf
Скачиваний:
3
Добавлен:
15.11.2022
Размер:
1.15 Mб
Скачать

Примеры решения задач Пример 1. Исходя из определения устойчивости по Ляпунову,

выяснить устойчивость решения дифференциального уравнения

dx a x с начальным условием x 1 0. dt t

Решение. Разделяем переменные dx a dt и получаем общее x t

решение x Cta . Частное решение, удовлетворяющее данному начальному условию, есть xч 0. Нетрудно установить, что любое другое частное решение, удовлетворяющее начальному условию x 1 x0 0, имеет вид xч x0ta . Разность произ-

вольного частного решения и частного решения при данном на-

чальном условии

равна

 

 

 

 

 

 

 

a

0 x0t

a

. Рассмотрим

 

xч xч x0t

 

 

 

различные случаи постоянной a:

 

 

 

 

 

 

 

 

 

1) Если a 0, то

 

 

 

 

 

 

 

x0t

a

 

 

 

x0

 

, а модуль разности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xч xч

 

 

 

 

 

 

 

 

 

 

 

x0 0

 

 

 

 

x0

 

 

 

 

 

 

 

начальных условий

 

 

 

 

 

 

. Следовательно, при

 

 

 

 

 

по определению решение устойчиво, но не асимптотически устойчиво.

2) Если a 0, то

 

 

 

 

x0t

a

t

a

 

x0

 

0приt . Следо-

 

 

 

 

 

xч xч

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вательно, решение асимптотически устойчиво.

 

 

3) Если a 0, то

 

 

 

 

 

x0t

a

 

 

 

t

a

 

x0

 

приt .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xч xч

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, решение неустойчиво.

Пример 2. Исследовать на устойчивость все положения равно-

 

dx

2x y x3

 

 

 

 

 

 

 

весия системы уравнений

 

 

.

dt

 

 

dy

x 2y 3x

5

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

91

 

Решение. Положения равновесия данной системы определяют-

 

 

2x y x

3

0

 

ся из системы уравнений

 

 

Отсюда находим

 

 

 

 

.

 

x 2y 3x5

0

 

 

 

 

 

 

 

 

три положения равновесия O 0,0 ,

A 1,1 и

B 1, 1 . Иссле-

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

производные

 

2x y x3

2 3x2 ,

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

2x y x3 1,

 

 

 

x 2y 3x5 1 15x4,

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

x

 

 

 

 

 

x 2y 3x5 2.

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

O 0,0 получаем

 

 

 

 

 

 

1)

Для точки

a 2,

b 1,

c 1,

d 2.

Поэтому соответствующая система первого приближения имеет

 

 

dx

2x y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вид

dt

 

 

 

.

Ее характеристическое

уравнение

 

 

 

 

 

 

 

 

dy

x 2y

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

2

 

1

 

2 4 5 0.

Корни этого

уравнения

 

 

 

 

1

2

 

 

 

 

 

 

 

 

1,2

2 i . Значит,

положение

равновесия O 0,0

является

устойчивым фокусом.

 

 

 

2) Для точки A 1,1

получаем a 1, b 1, c 14, d 2.

Поэтому соответствующая система первого приближения име-

 

 

d x 1

 

 

x 1

 

y 1

 

 

 

 

 

 

ет вид

 

dt

 

 

 

 

 

. Ее характеристическое уравне-

 

 

 

 

 

 

 

 

 

 

 

 

d y 1

 

x 1 2 y 1

 

 

 

 

 

14

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

92

ние

 

1

 

1

2 16 0.

Корни

этого уравнения

 

 

 

14

2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

0. Значит, положение равновесия

65

0

и

2

65

1

 

2

 

 

 

 

2

 

 

 

 

 

A 1,1

 

 

 

 

 

 

 

 

 

 

неустойчиво и является седлом.

 

3) Для точки

B 1, 1 получаем

a 1,

b 1, c 14, d 2.

Поэтому соответствующая система первого приближения имеет

 

 

 

 

 

d x 1

 

 

 

x 1

 

y 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вид

 

 

 

dt

 

 

 

 

 

 

. Ее характеристическое уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d y 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14 x 1 2 y 1

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

1

 

1

 

 

 

2 16 0. Корни этого уравнения

65

0

 

 

 

 

14

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

1

 

 

 

 

B 1, 1

 

и 2

65

 

0. Значит, положение равновесия

не-

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

устойчиво и является седлом.

Задачи и упражнения для самостоятельного решения

В задачах №№ 1, 2 исходя из определения устойчивости по Ляпунову, выяснить устойчивость решения дифференциальных

уравнений dx f t,x с начальным условием x 0 0. dt

1)

dx

1 t x.

2)

dx

sin2 x.

 

 

 

dt

 

dt

В задачах №№ 3–10 исследовать тип положения равновесия систем и построить фазовые траектории.

3)dx 3x 2y dt

5)

dx

2y,

dy

 

dt

 

dt

,

dy

x 4y.

4)

dx

x,

dy

x 2y.

 

 

 

 

 

 

dt

 

 

 

dt

dt

2x 3y .

6)

dx

x y ,

 

dy

x y .

 

 

 

 

 

 

 

dt

 

dt

93

7)

dx

4x 2y,

dy

 

2x y.

8)

dx

0,

 

dy

x y.

 

dt

 

 

dt

 

dt

 

 

 

 

dt

 

 

 

 

 

 

9)

dx

x y ,

 

dy

x y.

10)

dx

 

x,

 

dy

y.

 

 

 

dt

 

 

dt

 

 

 

dt

 

 

 

dt

В задачах №№ 11, 12 исследовать на устойчивость все положения равновесия систем.

 

dx

 

 

 

dx

 

11)

 

 

 

x x y 2

 

, 12)

 

 

 

sin y

 

 

dt

dt

 

 

.

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

y 1 x

 

 

 

 

 

 

sin x

 

 

 

dt

 

 

dt

 

Форма отчетности: конспект, устный опрос.

ЗАНЯТИЕ № 26

ПРИБЛИЖЕННОЕ РЕШЕНИЕ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

Литература: [33], c. 14-16.

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

Контрольные вопросы и задания

1. В чем состоит метод Эйлера приближенного решения дифференциального уравнения вида y f x, y ? Опишите ал-

горитм метода?

2.Что такое ломаная Эйлера? Как она соотносится с интегральной кривой?

3.Каков алгоритм метода Адамса приближенного решения дифференциальных уравнений?

4.Как выводится формула Адамса? Как при этом используется формула Тейлора?

94

5.Какой метод более точный: метод Эйлера или метод Адамса?

6.В чем заключается метод Рунге-Кутта приближенного решения дифференциальных уравнений?

7.Какие Вы еще знаете методы численного решения дифференциальных уравнений?

8.Как ищется приближенное решение систем дифференциальных уравнений первого порядка?

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

Форма отчета: программа и результаты счета.

ЗАНЯТИЕ № 27

ПРИМЕНЕНИЕ СТЕПЕННЫХ РЯДОВ

КВЫЧИСЛЕНИЮ ОПРЕДЕЛЕННЫХ ИНТЕГРАЛОВ

ИИНТЕГРИРОВАНИЮ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

Литература: [12], c. 290-294

.

Контрольные вопросы и задания

1.При каких условиях степенной ряд можно почленноинтегрировать и интеграл от суммы равен сумме интегралов?

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

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

4.Как найти приближенное решение дифференциального уравнения в виде степенного ряда?

5.Если решение представимо в виде ряда Тейлора, как находится приближенное решение?

6.В чем состоит метод решения дифференциального уравнения с помощью степенных рядов?

95

Примеры решения задач Пример 1. С помощью разложения в степенной ряд вычислить

1 2

 

dt

 

 

интеграл I

 

 

с точностью до 0,0001.

 

 

 

1 t4

0

 

 

Решение. Разложим подынтегральную функцию в биномиаль-

ный

 

ряд

1 x m 1

m

 

x

m m 1

x2

,

полагая в

нем

 

1!

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

2!

1 3

 

 

 

 

 

 

 

 

 

 

x t

4

,

m

:

 

 

 

 

 

 

 

 

 

1

 

t

4

 

 

t

8

. Этот ряд схо-

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2 2 2

 

 

 

 

 

 

1 t4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дится при

t

1. Интегрируя его, найдем

 

 

 

 

 

 

 

 

1 2

 

 

 

 

 

 

 

 

 

 

1 2

 

 

 

 

 

1

 

4

 

 

 

 

3

 

8

 

 

 

 

 

 

 

 

 

 

 

 

t5

t9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

1

 

t

 

 

 

 

 

 

t

 

 

 

 

dt t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

10

 

24

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

0,50000 0,00313 0,00008

 

 

 

 

 

 

 

 

 

 

24 29

 

 

 

 

2

 

10 25

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

Полученный результат представляет знакочередующийся сходящийся числовой ряд. Взяв сумму первых двух его членов, получим приближенное значение интеграла с заданной точностью: I 0,50000 0,00313 0,49687 0,4969, так как абсолютное значение третьего члена меньше 0,0001. Заметим здесь, что промежуточные вычисления проводятся с одним лишним знаком после запятой.

Пример 2. Найти решение дифференциального уравнения

y y y 0 с начальными условиями y 0 1, y 0 0. x

Решение. Полагая, что искомое решение представляет сходящийся степенной ряд

96

y a0 a1x a2x2 akxk akxk ,

k 0

найдем ряды для y и y его почленным дифференцированием

y a1 2a2x 3a3x2 k 1 ak 1xk k 2 ak 2xk 1

,

y 2a2 6a3x 12a4x2 k 2 k 1 ak 2xk .

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

первых коэффициентов:

y 0 a

1,

y 0 a

0. Подстав-

ляя ряды для y , y и

 

0

 

1

 

y в исходное уравнение и сделав приве-

дение подобных слагаемых, получим

1 4a2 9a3x a2 16a4 x2 a3 25a5 x3

ak k2 4k 4 ak 2 xk 0.

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

x0

1 22a

0

 

 

 

 

 

 

2

 

 

 

x1

32a3 0

 

 

x2

a

2

 

42a

 

0

 

 

 

 

 

4

 

 

x3

a

 

52a

0

 

 

 

 

3

 

5

 

 

 

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

 

 

 

 

 

 

 

xk

a

k

k 2 2 a

0

 

 

 

 

 

k 2

 

 

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

 

 

 

 

 

 

 

из которой определяем значения остальных коэффициентов:

a a a a

 

0,

a

1

,

a

1

,

 

22

2242

3 5 7

2m 1

 

2

 

4

 

 

1

 

 

 

 

1 m

 

 

1 m

 

a

, …,

a

 

 

 

 

,…

 

 

 

 

 

6

 

 

2m

 

224262 2m 2

 

4m m! 2

 

224262

 

 

 

 

 

 

 

 

97

 

 

 

 

 

 

Таким образом, искомое частное решение дифференциального уравнения имеет вид

 

x

2

 

x

4

 

1

m

x

2m

 

1

m

x

2m

y 1

 

 

 

 

 

 

 

 

 

4 1! 2

42 2! 2

4m m! 2

 

 

 

 

 

 

 

m 0 4m m! 2

.

Пример 3. Найти первые шесть членов разложения в ряд решения уравнения y xsin y , удовлетворяющего начальным ус-

ловиям y 1 0, y 1 .

2

Решение. Будем искать решение уравнения в виде ряда Тейло-

 

 

 

 

y

k

1

x 1 k . Подставляя в исходное уравнение

ра:

y x

 

 

 

 

 

 

 

 

 

 

k 0

 

k!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

1 и

y 1

,

находим

y 1 1 sin

1. Далее,

последо-

 

 

0

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вательно дифференцируя уравнение, имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y sin y xy cos y ,

 

 

y 1 1,

 

 

 

 

 

 

y

IV

y

 

 

 

y

 

 

 

 

 

 

 

 

sin y

 

 

 

 

 

 

 

cosy

 

cosy

xy

cosy xy y

 

 

 

 

2y cos y xy cos y x y 2 sin y ,

yIV 1 1,

yV 2y cos y 2 y 2 sin y y cos y xyIV cosy

xy y sin y y 2 sin y 2xy y sin y x y 3 cosy

3y cos y 3 y 2 sin y 3xy y sin y xyIV cosy

x y 3 cos y , yV 1 3 3 6.

Так как первое слагаемое ряда y 1 0, то вычислим еще

yVI 3yIV cos y 3y y sin y 6y y sin y 3 y 3 cosy3y y sin y 3x y 2 sin y 3xy yIV sin y 3x y 2 y cos y

98

yIV cos y xyV cos y xyIV y sin y y 3 cosy

3x y 2 y cos y x y 4 sin y ,

yVI 1 3 6 3 3 3 1 1 10.

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

y x x 1 1 x 1 2 1 x 1 3 1 x 1 4

2

2

6

24

1 x 1 5 1 x 1 6 .

20 72

Задачи и упражнения для самостоятельного решения

Решить задачи №№ 2930, 2931, 2935, 2936 [13]; 12.325, 12.326, 12.327, 12.329, 12.332 [19].

Форма отчетности: устный опрос, контрольная работа.

ЗАНЯТИЕ № 28

РАЗЛОЖЕНИЕ В РЯД ФУРЬЕ ФУНКЦИЙ, ЗАДАННЫХ НА ИНТЕРВАЛЕ (0, l)

Литература: [12], 331-333; [20], 257-259; [36], c. 31-35.

Контрольные вопросы и задания

1.Какой ряд называется тригонометрическим?

2.Сформулируйте достаточный признак разложимости функций в ряд Фурье.

3.Как разложить в ряд Фурье периодическую функцию с пе-

риодом 2 ?

4. Выведите формулы для вычисления коэффициентов ряда Фурье, если функция имеет период 2l .

99

5. Как разложить в ряд Фурье функцию, заданную на интервале

0, l ?

6.Что означает: "продолжить функцию четным образом", "нечетным образом"?

7.Как будет выглядеть график функции, продолженной четным или нечетным образом?

8.Как при таких продолжениях определяются коэффициенты разложения в ряд Фурье?

Примеры решения задач Пример 1. Разложить в ряд Фурье функцию, заданную на ин-

тервале 0, 2 :

1,

0

x 1

f x

 

.

 

x, 1

x 2

Решение. Продолжаем периодически данную функцию на всю числовую ось (рис. 14).

y

2

1

-4

-2

0

1

2

4

x

Рис. 14 Вычисляем коэффициенты ряда Фурье по формулам , учитывая, что T l 2:

 

 

2

2

1

 

2

 

10

 

x

2

 

 

2

5

 

 

 

 

 

 

 

 

a0

 

f x dx 1dx xdx x

 

 

 

 

 

,

 

 

 

 

 

 

2

 

2

 

 

2

 

 

0

0

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

1

2

 

 

 

 

 

 

 

ak f x cosk xdx

cosk xdx xcosk xdx

0

 

 

 

 

0

1

 

 

 

 

 

 

 

100

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]