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

даишев лекция по дифурам

.pdf
Скачиваний:
96
Добавлен:
10.02.2015
Размер:
959 Кб
Скачать
☆

Можно доказать (мы этого делать не будем), что если x1 = ϕ1(t) удовлетворяет уравнению (12.8), то, определив из (12.7) x2 = ϕ2(t), x3 = ϕ3(t), ... , xn = ϕn(t), получим, что набор функций { ϕ1(t), ϕ2(t), ... , ϕn(t) } есть решение нашей системы уравнений (12.4).

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

Замечание. Если указанный выше процесс применить к

линейной однородной системе

dxi

=

n

aik(t)xk, i = 1, 2, ..., n,

 

 

 

 

 

 

dt

 

k=1

 

то (12.8) тоже будет линейным однородным уравнением n-го по-

рядка. Если при этом aik = const,

то и (12.8) будет линейным

однородным уравнением с постоянными коэффициентами. Аналогичное замечание справедливо и для неоднородной системы

dxi

 

n

 

 

=

 

aik(t) xk + ψi(t).

dt

k=1

 

 

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

Определение. Интегрируемой комбинацией называется дифференциальное уравнение, являющееся следствием комбинирования уравнений системы (12.4), но уже легко интегриру-

91

ющееся: d Φ(t, x1, x2, ..., xn) = 0. Отсюда

 

Φ( t, x1, x2, ... , xn ) = C.

(12.9)

Данное выражение называется первым интегралом системы (12.4). Другими словами, первым интегралом Φ(t, x1, x2, ... , xn) =

C системы (12.4) называется соотношение, обращающееся в тождество при некотором C, если x1, x2, ... , xn заменить решением x1 = ϕ1(t), x2 = ϕ2(t), ... , xn = ϕn(t) системы (12.4).

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

Если найдено

s интегрируемых комбинаций, то может

быть получено s

первых интегралов Φ1(t, x1, x2, ..., xn) = C1,

Φ2(t, x1, x2, ..., xn)

= C2, ..., Φs(t, x1, x2, ..., xn) = Cs, и ес-

ли эти первые интегралы функционально независимы, то есть

D(Φ1, Φ2, ..., Φs) = 0, то s независимых функций из набора

D(xj1, xj2, ..., xjs)

x1, x2, ... , xn можно выразить через остальные и, подставляя их в систему (12.4), придем к системе уравнений с меньшим числом неизвестных. При s = n, а также когда все интегралы независимы, все неизвестные функции могут быть определены из системы Φ1(t, x1, x2, ..., xn) = C1, Φ2(t, x1, x2, ..., xn) = = C2, ..., Φn(t, x1, x2, ..., xn) = Cn.

Систему (12.4) иногда удобно представить в виде

dx1

 

dx2

 

dxn

 

dt

 

=

 

=

... = f (t, x)

=

 

.

f (t, x)

f (t, x)

1

1

2

 

n

 

 

 

92

Обозначая t = x0 и переобозначая удобным образом переменные, можно получить систему, записанную в симметричной форме

dx0

 

dx1

 

 

dxn

 

=

 

=

... =

An(x0, x1, ..., xn).

A0(x0, x1, ..., xn)

A1(x0, x1, ..., xn)

Преимуществом такой формы записи системы является то, что в эту систему все переменные x0, x1, x2, ..., xn входят равноправно, тогда как в нормальной системе (12.4) такого равноправия нет: x1, x2, ..., xn рассматриваются как функции, а t — как независимая переменная.

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

 

a1

=

a2

= ... =

an

= t,

b1

 

 

 

 

b2

bn

то при любых k1, k2, ... , kn

имеем

 

 

k1a1 + k2a2 + ... + knan = t. k1b1 + k2b2 + ... + knbn

ЛЕКЦИЯ 13

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

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

ний

dxi

 

n

 

 

 

=

 

aik(t)xk + fi(t) (i = 1, 2, ..., n).

(13.1)

dt

k=1

 

 

 

93

Эта система может быть записана в матричной форме

 

 

 

 

 

 

 

 

dX

= AX + F,

 

 

 

 

(13.1 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a11(t) a12(t) ... a1n(t)

 

 

x1(t)

 

 

f1(t)

 

 

a

21

(t) a

22

(t) ... a

2n

(t)

 

 

x

(t)

 

 

f (t)

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

A =

 

 

 

 

 

 

 

 

 

 

, X =

 

 

 

 

, F =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

........

 

 

........

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an1(t) an2(t) ... ann(t)

 

 

xn(t)

 

 

fn(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если все fi(t) ≡ 0, то есть матрица-столбец

F (t) ≡ 0, система

называется однородной. При

F (t) = 0, то есть когда хотя бы

одна fi(t) = 0, система называется неоднородной.

 

Введём в рассмотрение линейный оператор

 

 

 

 

 

 

 

 

def

d

 

 

 

 

 

 

 

 

L ( )

=

dt − A.

 

 

 

(13.2)

Тогда систему уравнений (13.1 ) можно записать в виде

L (X) = 0

для однородной системы и L (X) = F (t)

для неод-

нородной.

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим некоторые свойства оператора L(

).

Свойство 1. Оператор

L(

) линеен: L (C1X1 + C2X2) =

= C1L(X1) + C2L(X2).

 

 

 

 

 

 

 

 

 

 

Действительно,

 

d

 

 

 

 

 

 

 

 

 

(C1X1 + C2X2) −A (C1X1

+ C2X2) =

 

dt

= C1

dX1

−

AX1 + C2

 

dX2

−

AX2 = C1 L(X1) + C2 L(X2).

dt

dt

 

 

 

 

 

m

 

 

Свойство 2. Линейная комбинация

CsXs

с посто-

s=1

 

 

 

 

 

 

 

 

 

 

 

 

линейной

янными коэффициентами решений X1, X2,

... , Xs

однородной системы также является решением этой системы дифференциальных уравнений.

94

Доказательство. L m

CsXs =

m

Cs L (Xs) = 0.

 

 

s=1

 

s=1

 

 

 

 

 

 

 

 

 

=0

 

 

 

 

 

 

$

 

 

 

 

 

 

 

 

%& '

Свойство 3. Если линейная однородная система имеет

комплексное решение

X = U (t) + iV (t), то

 

 

 

 

u1(t)

 

 

v1(t)

 

 

u

(t)

 

 

v

(t)

 

 

2

 

 

 

 

2

 

 

U (t) =

 

 

 

, V (t) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.......

 

 

.......

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

un(t)

 

 

vn(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— также решения этой системы.

Доказательство. L(X) = L(U + iV ) = L(U ) + iL(V ) = = 0 L(U ) = 0, L(V ) = 0, поскольку комплексная величина равна нулю тогда и только тогда, когда действительная и мнимая части этой комплексной величины равны нулю. Таким образом, U и V — решения системы.

Определение. Векторы X1, X2, ... , Xn называются линейно зависимыми на сегменте t [ a, b ], если существуют такие постоянные λ1, λ2, ... , λn, среди которых хотя бы одно число λk = 0, что линейная комбинация λ1X1 + λ2X2 + ...+ +λnXn = 0. Если же это соотношение выполнено тогда и только

тогда, когда все

λk = 0, то векторы

X1, X2,

... , Xn линейно

независимы.

 

 

 

 

 

 

 

 

 

 

 

 

Из матриц-столбцов

 

 

 

 

 

 

 

 

 

x11

 

 

 

x12

 

 

x1n

 

 

x

21

 

 

 

x

22

 

 

x

2n

 

 

 

 

 

 

 

 

 

 

 

X1 =

 

 

 

, X2

=

 

 

 

, ... , Xn =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

....

 

 

 

....

 

 

....

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn1

 

 

 

xn2

 

 

xnn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

95

составим квадратную матрицу

 

 

 

 

 

x11

x12

... x1n

 

 

 

 

 

 

 

 

 

 

 

 

x

21

x

22

... x

2n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W (t) =

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn2

... xnn

 

 

 

xn1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

которая называется матрицей Вронского. Здесь, как всегда, первый индекс обозначает строку, а второй — столбец. Определитель этой матрицы называется определителем Вронского или

вронскианом.

Теорема 1 (о линейной зависимости системы решений).

Если определитель Вронского решений X1, X2, ... , Xn однородной системы линейных дифференциальных уравнений с непрерывными на сегменте [ a, b ] коэффициентами aik(t) ра-

вен нулю хотя бы в одной точке t0

[ a, b

], то решения

X1, X2, ... , Xn линейно зависимы на

[ a, b ]

и, следователь-

но, det W (t) ≡ 0 на [ a, b ].

 

 

Доказательство. Вследствие непрерывности коэффициентов aik(t), выполняется теорема существования и единственности решения для системы дифференциальных уравнений. Сле-

довательно, начальное значение

X(t0) = 0 определяет един-

ственное решение X(t) ≡ 0 однородной системы.

 

Поскольку

det W (t0) = 0,

то существует такая система

постоянных C1,

C2, ... , Cn, среди которых хотя бы одна Ci

отлична от нуля, что

 

 

C1X1(t0) + C2X2(t0) + ... + CnXn(t0) = 0,

(13.3)

96

что эквивалентно системе алгебраических уравнений

 

C1 x11(t0) + C2 x12(t0) + ... + Cn x1n(t0) = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

x

21

(t ) + C

x

22

(t ) + ... + C

x

2n

(t ) = 0,

 

 

1

0

2

0

n

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

C1 xn1(t0) + C2 xn2(t0) + ... + Cn xnn(t0) = 0

с определителем равным нулю. Соответствующее этой нетриви-

альной системе постоянных

 

C1,

C2, ... , Cn решение X(t) =

n

CiXi(t)

в силу (13.3) удовлетворяет начальному условию

i=1

 

 

 

 

 

 

 

 

 

 

 

X(t ) = 0. Но тогда X(t)

≡

0,

и мы на всём сегменте

[a, b]

 

0

 

n

 

 

 

 

 

имеем X(t) =

i=1

CiXi(t)

 

0, при этом хотя бы одна из

Ci

 

 

 

≡

 

 

 

 

 

 

 

 

 

 

 

 

 

 

отлична от нуля. Это означает линейную зависимость на

[a,

b]

системы решений

X1, X2, ... , Xn. Теорема доказана.

 

 

 

Теорема 2(об общем решении линейной однородной систе-

мы.) Линейная комбинация

 

n

CiXi(t) с произвольными по-

 

i=1

 

 

 

 

 

 

 

 

 

 

 

стоянными

n

 

линейно независимых решений Xi(t)

(i

=

1, 2, ..., n) однородной системы линейных уравнений с непрерывными на сегменте [a, b] коэффициентами aik(t) является общим решением этой системы дифференциальных уравнений.

Доказательство. Покажем, что какое бы частное реше-

ние, отвечающее начальным значениям

 

 

X(t0) = X0 (то есть

 

 

 

 

 

 

 

 

 

 

x1(t0) = x10,

x2(t0) = x20,

... , xn(t0) = xn0, ), ни взять, всегда

 

 

 

 

что

 

 

n

 

(Иначе гово-

найдутся такие Ci,

X(t)

=

 

Ci Xi(t).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

Ci в решении

ря, надо показать, что подбором постоянных

X(t) =

n

Ci Xi(t) можно удовлетворить любым произвольно

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

X(t0) = X0. )

 

выбранным начальным условиям

 

В силу теоремы существования и единственности решения,

97

CiXi(t)

это требование эквивалентно тому, что должно быть удовлетво-

n

рено векторное уравнение Ci Xi(t0) = X0, где Ci считаются

i=1

неизвестными величинами, или подробно

C1 x11(t0) + C2 x12(t0) + ... + Cn

x1n(t0) = x1,

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C x

21

(t ) + C

x

22

(t ) + ... + C

 

x

2n

(t ) = x

,

 

1

0

2

0

n

0

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1 xn1(t0) + C2 xn2(t0) + ... + Cn xnn(t0) = xn .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определитель данной системы — это определитель Вронского нашей линейно независимой системы решений. Поэтому он отличен от нуля, det W (t) = 0. Следовательно, система всегда имеет единственное решение, которое может быть найдено, например, по формулам Крамера.

Теорема 3(об общем решении линейной неоднородной системы). Если Y (t) – частное решение неоднородной системы

L (X(t)) = F (t), то X(t) = Y (t) +

n

Ci Xi(t),

где Xi(t)

 

 

 

 

i=1

— линейно независимые решения соответствующей однородной системы L (X(t)) = 0, а Ci — произвольные постоянные, есть общее решение линейной неоднородной системы дифференциальных уравнений.

Доказательство. Вместо неизвестного вектора

X(t) вве-

дём новый вектор Z(t) = X(t) − Y (t). Тогда

L (Z(t))

=

L (X(t)) −L (Y (t)) = 0, поскольку L (X(t)) = L (Z(t) + Y (t))

=

L (Z(t)) + L (Y (t)) = F (t) L (Z(t)) = 0.

 

 

$%& '

≡F (t)

Если X1(t), X2(t), ..., Xn(t) — система линейно независимых

n

решений однородной системы L (Z(t)) = 0, то Z(t) =

i=1

98

— общее решение этой однородной системы. Отсюда

X(t) =

Y (t) +

n

Ci Xi(t)

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

˜

 

Оно действительно общее, поскольку какое бы решение

 

Y (t),

отвечающее произвольным начальным данным

 

˜

 

=

˜

Y (t0)

Y0,

 

 

 

 

 

 

˜

˜

 

 

˜

что

мы не взяли, найдутся такие постоянные C1, C2, ... , Cn,

˜

 

 

 

n

˜

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y (t) = Y (t) +

 

CiXi(t), так как система алгебраических урав-

 

 

 

 

 

 

 

 

 

 

 

 

 

нений

 

 

i=1

 

 

 

 

 

 

 

 

 

x11(t0) + C˜2

x12(t0) + ... + C˜n x1n(t0) =

 

y˜1

 

 

y1(t0),

C˜1

 

 

 

 

 

 

 

 

 

 

 

0

 

−

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

˜

 

 

 

˜

˜

 

y˜2

 

 

y2(t0),

 

C1

x21(t0) + C2

x22(t0) + ... + Cn x2n(t0) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

˜

 

 

 

˜

˜

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1 xn1(t0) + C2 xn2(t0) + ... + Cn xnn(t0) =

 

y˜n

 

yn(t0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Теорема 4(о суперпозиции решений.)Решением системы

L (X(t)) = F1(t)+F2(t)+...+Fs(t) является сумма всех решений

X1(t), X2(t), ... , Xs(t) уравнений L (X(t)) = F1(t), L (X(t)) = F2(t), ... , L (X(t)) = Fs(t).

Доказательство. Дано:

L (X1(t)) ≡ F1(t), L (X2(t)) ≡ F2(t), ... , L (Xs(t)) ≡ Fs(t).

Тогда L (X1(t) + X2(t) + ... + Xs(t)) =

= L (X1(t))+L (X2(t))+...+L (Xs(t)) ≡ F1(t)+F2(t)+ ... +Fs(t),

то есть сумма решений X1(t)+X2(t)+ ... +Xs(t) действительно

99

является решением системы L (X(t)) = F1(t) + F2(t) + ...+ Fs(t). Как следствие этой теоремы имеем, что если система

L (X(t)) = U (t) + i V (t)

имеет решение ˜ ˜ то ˜ ≡ ˜ ≡

Z = X(t) + i Y (t), L X(t) U (t), L Y (t)

≡ V (t). Иначе говоря, функции

˜

и

˜

являются реше-

X(t)

Y (t)

ниями уравнений L (X(t)) = U (t) и

L (Y (t)) = V (t) .

13.2. Метод вариации произвольных постоянных

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

 

dX

− A(t)X = F (t)

(13.4)

 

 

 

dt

ищем в виде

 

n

 

 

 

 

 

 

 

X(t) =

 

Ci(t) Xi(t),

(13.5)

 

i=1

где X1(t), X2(t), ... , Xn(t) — система линейно независимых

решений однородной системы dXdt − A(t)X = 0, а Ci(t) — неизвестные пока функции. Для нахождения этих функций под-

ставим решение (13.5) в исходную систему (13.4):

n

Ci(t)Xi(t) +

n

 

 

dXi(t)

 

A(t)

n

Ci(t)Xi(t) = F (t).

 

 

Ci(t)

−

 

 

 

 

 

 

dt

 

 

i=1

 

 

i=1

 

 

 

 

 

i=1

 

Но Xi(t)

— решения однородной системы, то есть

 

 

 

 

 

dX

 

 

 

 

 

 

 

 

 

i

≡ A(t)Xi(t),

 

 

 

 

 

 

dt

 

и, следовательно, второе и третье слагаемые взаимно уничтожа-

100

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