Материал: Sb98337

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Правая часть уравнения является решением однородного уравнения, и потому частное решение надо искать по формуле

y* x Acos 2x Bsin 2x .

Вычислим производные y* :

y* Acos 2x Bsin 2x x 2Asin 2x 2B cos 2x( A 2Bx)cos 2x (B 2Ax)sin 2x,

y* 2B cos 2x 2Asin 2x 2 A 2Bx sin 2x 2 B 2Ax cos 2x.

Подставим эти выражения в уравнение

2B cos 2x 2Asin 2x 2 A 2Bx sin 2x 2 B 2 Ax cos 2x4x Acos 2x B sin 2x sin 2x.

Проведем тождественные преобразования

4B cos 2x 4 Asin 2x sin 2x A 14 , B 0.

Общее решение неоднородного уравнения

yy0 y* C1 cos 2x C2 sin 2x 14 x cos 2x.

3.4.Системы однородных линейных дифференциальных уравнений

спостоянными коэффициентами

Решение таких систем связано с несколькими темами, которые в реше-

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

мы снова сводим решения задачи о дифференциальных уравнениях к поиску

собственных чисел и собственных векторов матрицы:

поиск собственных чисел и собственных векторов;

решение систем линейных уравнений;

решение дифференциальных уравнений;

базис и линейные комбинации.

Рассмотрим систему X t AX t , где X t – набор

x1 t , x2 t , , xn t .

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

комбинации экспонент. Точнее говоря, находим собственные числа матрицы A. Сначала предположим, что они вещественные и различные – 1, 2, , n .

Для каждого из них найдем собственный вектор Yk. Тогда базис пространства

21

решений выглядит так:

e 1xY ,

e 2xY , ,

e n xY . Это вытекает из выполне-

 

 

 

 

1

2

n

 

 

 

 

 

ния двух условий:

 

 

 

 

 

 

 

 

X k k X (из формулы для производной экспоненты);

 

 

 

 

AX k k X k (поскольку k – собственное число); поэтому X k AX k .

 

Общее

решение

записывается в

виде:

C e 1x X

1

C e 2x X

2

 

 

 

 

 

 

 

 

1

2

 

C e n x X

n

.

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

Пример. Найти частное решение однородной системы

x1 4x1 x2x2 3x1 2x2

x3 2x1 3x2 4x3,

которая удовлетворяет условию x1 0 6, x2 0 6, x3 0 24.

Найдем собственные числа матрицы: 1, 4, 5. Собственные векторы:

3

0

1

 

 

 

 

 

 

 

9

, 0

, 1 .

 

 

 

 

 

 

 

7 1 5

 

 

 

 

 

 

 

 

 

 

3

 

0

 

 

1

 

Тогда базис множества решений имеет вид: 9 et ,

0

e4t ,

1

e5t .

 

 

7

1

5

Можно проверить подстановкой, что эти векторы подходят.

 

 

3

 

 

 

0

 

 

 

1

 

Общее решение имеет вид: X t C

 

9

et C

2

 

0

e4t C

 

1

e5t .

1

 

 

 

 

 

 

3

 

 

 

 

 

7

 

 

1

 

 

5

Для нахождения частного решения находим константы из системы:

6

 

3

 

 

0

 

1

X 0

6

 

C

9

 

C

2

 

0

C

 

1 .

 

 

 

 

1

 

 

 

 

 

 

3

 

 

 

24

 

 

7

 

 

 

1

 

 

5

Отсюда C1 1,C2 2,C3

3 .

 

 

 

 

 

 

 

 

 

 

3

 

0

 

 

 

3

 

 

 

 

 

 

 

 

Ответ: X t 9 et

0

e4t

3

e5t .

 

 

 

 

 

7

2

 

15

 

 

 

 

 

 

 

22

Рассмотрим теперь случай с кратным корнем для системы линейных дифференциальных уравнений. Если собственное число имеет кратность r ≥ 2, то решение системы ищут в виде вектора

a

a

t a

tr 1

 

 

11

12

1r

 

 

 

 

 

 

 

e t ,

 

 

a

t a

tr 1

 

a

 

 

11

12

1r

 

 

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

3.5. Метод изоклин

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

y y2, y 1 2, y 1 1.

Метод состоит в том, что сначала отмечают точку с координатами (1, 2),

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

небольшое расстояние (например: 1/2), затем для этой точки находят значение производной, угловой коэффициент новой касательной, и несколько раз

повторяют процесс.

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

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

надо составить характеристическое уравнение. В зависимости от вида и зна-

ков корней, можно определить вид точки покоя.

Пусть система дифференциальных уравнений имеет вид

x a11x a12 y

 

 

 

 

 

 

 

 

 

 

 

y a21x a22 y

 

 

 

 

 

a11

a

a12

 

0.

Составим характеристическое уравнение det

a

 

 

21

 

22

 

 

23

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

Для комплексных корней с отрицательной вещественной частью полу-

чим устойчивый полюс, с положительной вещественной частью – неустойчивый полюс, с нулевой вещественной частью – центр.

24

Список рекомендуемой литературы

Ефимов А. В., Демидович Б. П. Сборник задач по математике для ВТУЗов.

Ч. 2. М.: Наука, 1986.

Макаров Б. М., Подкорытов А. Н. Лекции по вещественному анализу. СПб.: БХВ-Петербург, 2011.

Камке Э. Справочник по обыкновенным дифференциальным уравнениям. СПб.: Лань, 2005.

25