Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
КРАТКИЙ КУРС ЛЕКЦИЙ ПО ВЫСШЕЙ МАТЕМАТИКЕ. Ч2.pdf
Скачиваний:
177
Добавлен:
27.05.2019
Размер:
1.54 Mб
Скачать

65

Имеют место следующие утверждения:

1)Каждому k -кратному действительному корню r характеристического уравнения соответствует k частных решений вида:

erx , xerx , , xk 1erx

2) Каждой паре t -кратных комплексно сопряженных корней r1 = α + βi , r2 =α βi характеристического уравнения соответствует 2t частных

решений вида:

 

 

eαx cos βx ,

xeαx cos βx ,

xt 1eαx cos βx ,

eαx sin βx ,

xeαx sin βx ,

xt 1eαx sin βx .

Общая сумма кратностей всех корней должна равняться степени характеристического уравнения n .

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

8.1Общие определения. Нормальные системы уравнений

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

Всегда предполагается, что число уравнений равно числу неизвестных функций. Независимую переменную будем обозначать буквой t , а неизвестные функции этой переменной или через x1 (t) , x2 (t) , …, xn (t) , или, если их

не больше трех, через x(t) , y(t) , z(t) . Производные по t обозначаются точ-

ками.

Решением системы дифференциальных уравнений называется совокупность функций x1 = x1 (t) , x2 = x2 (t) , …, xn = xn (t) , которая при подстановке

в каждое из уравнений превращает его в тождество.

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

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

 

 

 

x1 = f1 (t, x1 , x2 , K, xn ),

 

 

&

 

 

 

 

&

= f2 (t, x1

, x2

, K, xn ),

 

x2

(1)

 

 

 

 

 

K K K K K

 

 

&

= fn (t, x1

, x2

, K, xn ).

 

xn

 

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

Пример 1. Решить систему

x& =y =&z& =

y,

z,

xy + z.

--------------------------------------------------

© Тимофеев В.А., Тимофеев А.А., 2011

66

Решение. Продифференцируем первое уравнение по t и заменим производную y& ее выражением из второго уравнения

&x&= y& = z .

Продифференцируем еще раз и заменим z& ее выражением из третьего уравнения

&x&&= z& = x y + z .

Т.к. y = x&, z = &x&, то окончательно получим

&x&&− &x&+ x&− x = 0 – линейное дифференциальное уравнение 2-го порядка. Характеристическое уравнение r3 r 2 + r 1 = 0 , или (r 1)(r 2 +1) = 0 , имеет корни r1 =1, r2,3 = ±i . Следовательно,

x = C1 et +C2 cost +C3 sin t .

Т.к. y = x& и z = &x&, то

y= C1 et C2 sin t +C3 cost ,

z= C1 et C2 cos t C3 sin t .

Общее решение нормальной системы (1) имеет вид

x

=ϕ (t, C

, C

, K, C

),

1

1

1

 

2

 

n

 

 

x2 = ϕ2 (t, C1 , C2 , K, Cn ),

K K K K K K

x

=ϕ

(t, C

 

, C

 

, K, C

 

),

n

n

1

2

n

 

где C1 , C2 , …, Cn – произвольные постоянные.

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

x1 t =t0 = x10 , x2 t =t0 = x20 , …, xn t =t0 = xn0 .

Подставляя начальные условия в общее решение, получаем систему уравнений для определения произвольных постоянных

 

ϕ

(t, C

, C

, K, C

) = x

 

,

 

1

1

2

 

n

10

 

 

 

ϕ2 (t, C1 , C2 , K, Cn ) = x20 ,

 

 

 

 

 

 

 

 

 

 

K K K K K K

 

 

 

ϕ

(t, C

, C

 

, K, C

) = x

n0

.

 

n

1

2

n

 

 

Теорема. Если правые части нормальной системы непрерывны вместе со своими частными производными в окрестности значений t0 , x10 , x20 , …,

xn0 , то существует единственная система функций x1 (t) , x2 (t) , …, xn (t) , яв-

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

67

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

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

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

x1 = a11 x1 + a12 x2 +K+ a1n xn ,

 

 

 

&

 

 

 

&

= a21x1 + a22 x2 +K+ a2n xn

,

 

x2

(2)

 

 

 

 

K K K K K K K

 

 

 

&

= an1 x1 + an2 x2 +K+ ann xn

,

 

xn

 

или в векторной форме

x&1 x& = Ax , где x& = M ,

x&n

x

 

a

a

K a

 

 

11

12

1n

 

1

 

a21

a22

K a2n

x = M

 

, A =

K

K

K K

.

 

 

 

 

xn

 

 

 

an2

 

 

 

 

an1

K ann

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

 

a11 λ

a12

K

a1n

 

 

 

 

 

 

 

a21

a22 λ K a2n

 

 

= 0 .

 

K

K

K

K

 

 

 

 

an1

an2

K ann λ

 

 

 

Пусть характеристическое уравнение имеет n различных корней λ1 ,

λ2 , …, λn , Пусть каждому λk

соответствует

собственный вектор

( p1k , p2k , K, pnk ) , где k =1, 2, K, n .

Тогда система дифференциальных уравнений имеет n решений: 1-е решение, соответствующее корню λ = λ1

x11 = p11eλ1t , x21 = p21eλ1t , …, xn1 = pn1eλ1t ;

2-е решение, соответствующее корню λ = λ2

x12 = p12eλ2t , x22 = p22eλ2t , …, xn2 = pn2eλ2t ;

… … … … … … … … … … …

n -е решение, соответствующее корню λ = λn

x1n = p1neλnt , x2n = p2neλnt , …, xnn = pnneλnt .

Мы получили фундаментальную систему решений. Общее решение системы таково

 

x

= C

x

 

+C

x

 

+K+C

x

 

,

 

1

1

11

2

12

n

1n

 

 

x2

= C1 x21 +C2 x22 +K+Cn x2n ,

 

K K K K K K K K

 

 

 

x

= C

x

n1

+C

x

n2

+K+C

x

nn

.

 

n

1

 

2

 

n

 

 

68

Если среди простых корней характеристического уравнения есть комплексные, то и решение получится в виде комплексных функций. Если при этом коэффициенты системы (2) вещественны, то решение можно выразить только через вещественные функции. Для этого надо воспользоваться тем, что вещественная и мнимая части комплексного решения, соответствующего корню λ = α + βi ( β 0 ), является линейно независимыми решениями.

 

 

 

Если характеристическое уравнение имеет корень λ кратности m , то

этому

корню

соответствует

решение

x

= P (t) eλt ,

x

2

= P (t) eλt , …,

 

 

 

 

 

 

 

1

1

 

2

x

n

= P (t) eλt ,

где P (t) ,

P (t) , …, P (t)

многочлены степени не выше

 

 

n

1

2

n

 

 

 

 

 

m 1.

Пример 2. Решить систему

x&

y&

=7x + y,

=2x 5y.

Решение.

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

 

 

 

 

 

 

 

 

 

 

7 λ

1

 

 

 

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

5 λ

 

 

 

 

 

 

 

 

 

(7 λ)(5 λ) + 2 = 0 , (7 + λ)(5 + λ) + 2 = 0 ,

 

λ2 +12λ +37 = 0 .

 

Решим полученное уравнение:

 

 

 

 

12 ± 2i

 

 

 

 

 

 

D =144 4 37 =144 148 = −4 ,

λ

=

= −6 ±i .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,2

 

 

 

2

 

 

 

 

 

 

 

 

Возьмем корень λ = −6 +i . Найдем собственный вектор μ(a;b) :

 

 

( A λE)μ = 0

(1i)a +b = 0,

 

b = (1+i)a,

 

a

=1,

 

2a + (1i)b

= 0,

 

 

 

 

 

 

 

=1+i.

 

 

 

 

 

 

2a = 2a,

 

 

 

b

Решение, соответствующее корню λ = −6 + i

имеет вид

 

 

 

 

 

x = e(6+i)t = e6t (cost +i sin t) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

y = e(6+i)t (1+i) = e6t (cost sin t +i(sin t + cos t)) .

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

= Re e(6+i)t

= e6t

cos t,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

= Re(1+i)e(6+i)t = e6t (cost sin t),

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

= Im e(6+i)t

= e6t

sin t,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= Im(1+i)e(6+i)t = e6t (sin t

+ cos t).

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

Общее решение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x = C1

x1 +C2 x2

,

 

 

x =

C e6t

cos t +C

 

e6t sin t,

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

y = C1

x1 +C2 y2 ,

 

 

 

y =

C

e6t (cos t sin t) +C

2

e6t (sin t + cos t).

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

69

Пример 3.

Решить систему

 

 

 

+3x ,

 

 

 

 

 

 

 

 

 

 

x = 7x

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(*)

 

 

 

 

 

&1

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

x2 = 6x1 + 4x2 .

 

 

 

 

 

Решение. Составим характеристическое уравнение:

 

 

 

 

 

 

 

7 λ

 

3

 

 

 

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

4 λ

 

 

 

 

 

 

 

 

(7 λ)(4 λ) 18 = 0 , λ2 11λ +10 = 0 .

 

 

 

 

 

 

Решим полученное уравнение:

 

11±9

 

 

 

 

 

 

 

 

 

 

 

 

 

D =1214 10 = 81,

λ

=

,

 

λ =1,

λ

2

=10 .

 

 

 

 

 

 

 

 

 

1,2

 

 

2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Возьмем корень

λ =1. Найдем собственный вектор μ(a;b) :

 

 

( A λE)μ = 0

 

6a +3b = 0,

 

 

2a +b = 0

a

=1,

 

 

 

 

 

 

 

 

 

 

= −2.

 

 

 

6a +3b = 0,

 

 

 

 

 

 

 

 

 

 

b

Решение, соответствующее корню λ =1 имеет вид

 

 

 

 

 

 

x = et

,

 

x

21

 

= −2et .

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Возьмем корень

λ =10 . Найдем собственный вектор μ(a;b) :

 

 

( A λE)μ = 0

3a +

3b = 0,

 

 

 

a b = 0

a =1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6a 3b = 0,

 

 

 

 

 

 

 

 

 

 

b =1.

Решение, соответствующее корню λ =10 имеет вид

 

 

 

 

 

x = e10t ,

x

22

= e10t .

 

 

 

 

Общее решение

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

= C

et +C

2

e10t

,

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2-й способ.

 

x2 = −2C1 et +C2 e10t .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, корни характеристического уравнения

λ1 =1, λ2

=10 . Исходя из

структуры общего решения, будем искать общее решение в виде

 

 

 

 

x = C

et +C

2

e10t ,

 

 

 

(**)

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 = a et +b e10t .

 

 

 

 

 

Выразим a и b через C1 и C2 . Для этого подставим (**) в (*):

 

 

C

et +10C

2

e10t

= 7C et + 7C

2

et +3a et +3b e10t ,

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

a et +10b e10t = 6C1et +6C2 e10t + 4a et + 4b e10t .

 

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

C

= 7C

+3a,

 

a = −2C

,

 

1

1

 

 

1

 

 

10C2 = 7C2 +3b,

 

b = C2 .

 

 

Таким образом, общее решение имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

70

 

 

 

x = C et +C

2

e10t ,

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

x2 = −2C1 et +C2 e10t .

 

 

Пример 4. Решить систему

 

 

 

 

 

 

 

 

 

 

 

 

x = y,

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

&

 

t

+ e

t

.

 

 

 

 

 

y = x + e

 

 

 

 

Решение. Решим однородную систему

 

 

 

 

 

 

 

 

 

 

 

 

x = y,

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

(3*)

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

y = x.

 

 

 

 

 

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

 

 

 

 

 

 

 

 

λ 1

 

= 0 , λ2 1 = 0 , λ =1, λ

2

= −1.

 

 

 

1 λ

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение будем искать в виде

 

 

 

 

 

 

 

 

 

 

 

 

x

= C

et +C

2

et ,

 

 

 

 

 

 

1

 

 

 

 

 

 

(4*)

 

 

 

= a et +b

et .

 

 

 

 

y

 

 

Выразим a и b через C1

и C2 . Для этого подставим (4*) в (3*) и приравняем

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

 

 

 

 

 

 

C

et C

2

et = a et +b et ,

 

a

= C1 ,

 

 

 

 

 

 

 

 

1

 

 

 

 

C

et +C

 

 

et ,

 

 

= −C2 .

 

 

 

 

 

 

 

 

a et b et =

2

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x =

C

et

+C

2

et

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y =

1

et

C

et

– общее решение однородной системы.

 

 

 

C

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть

C1 = C1 (t) , C2 = C2 (t) . Тогда

 

 

 

 

 

 

 

 

 

 

 

x = C (t) et +C

2

(t) et

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

(t) et

 

– подставим в исходную систему уравнений.

 

 

y = C (t) et C

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C(t)

et

+C (t) et

+C(t) et C

2

(t) et

= C (t) et

C

2

(t) et

,

 

 

1

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

1

 

 

 

C1(t) et +C1 (t) et C2(t) et +C2 (t) et = C1 (t) et +C2 (t) et

+ et + et ,

 

 

C

(t) et +C(t)

et = 0,

 

 

 

 

 

C(t) = −C(t) e2t ,

 

 

 

1

 

 

 

 

2

 

 

 

 

 

 

 

 

+ et

,

 

2

1

 

 

 

 

 

C1(t) et C2(t) et = et

 

 

2C1(t) et = et + et ,

 

 

C

(t) = 0,5(1+ e2t ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

+1).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

C2(t) = −0,5(e2t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1+ e

2t

) dt =

 

(t 0,5e

2t

) +C1 ,

 

 

 

 

 

 

 

C1 (t) =

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

71

C2

1

(e

2t

+1) dt = −

1

(t + 0,5e

2t

) +C2 .

(t) = −

2

 

2

 

 

 

 

 

 

 

 

Таким образом, общее решение неоднородной системы уравнений

x = C et

+ 0,5(t 0,5e2t )et +C

2

et

0,5(t + 0,5e2t )et =

1

 

 

 

 

 

 

 

 

 

= C

et +C

2

et

+t et et

1 et + et

,

 

 

1

 

 

2

 

2

2

 

 

 

 

 

 

 

 

 

 

 

y = C et +0,5(t 0,5e2t )et +0,5(t + 0,5e2t )et C

2

et =

1

 

 

 

 

 

 

 

 

 

= C

et C

2

et

+t et + et

+ 1 et et .

 

 

1

 

 

2

 

2

2