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

3447

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

(x)

A(x)Ф(x) ,

(5.37)

 

dx

 

 

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

Отметим: правило дифференцирования суммы и произведения сохраняются и для матриц; однако при

дифференцировании

произведения матриц

необходимо

сохранять

 

порядок

сомножителей

-

(A(x)B(x))

A (x)B(x) A(x)B (x) .

 

 

Общее решение однородной линейной системы

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

 

 

 

Как

и

в

случае

однородного

линейного

дифференциального

уравнения

n -го порядка,

знание

фундаментальной системы решений позволяет построить

общее

решение

однородной

линейной

системы

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

 

 

 

Теорема 6 (об общем решении). Если система функций

 

 

y11 ,

y12 ,

y1n,

 

 

 

y21,

y22,

y2n ,

(5.38)

 

 

 

 

 

 

 

 

 

 

 

yn1,

yn2,

 

ynn

 

является фундаментальной системой решений однородной линейной системы дифференциальных уравнений (5.23) в интервале (a,b) , то функции

y1 (x)

C1 y11

C2 y21 Cn yn1,

 

y2 (x)

C1 y21

C2 y22

Cn yn2 ,

(5.39)

 

 

yn (x)

C1 yn1

C2 yn2

Cn ynn ,

 

или в сокращенной форме

 

 

 

n

 

 

 

yk (x) Ci yik (x),

(k 1,2,..., n),

(5.40)

 

i 1

 

 

 

180

(где C1, C2 ,..., Cn

произвольные

постоянные) дают общее

решение системы (5.23) в области

 

 

 

 

 

 

 

 

 

 

 

a

x

b,

yk

(k

1,2,..., n) .

(5.41)

*Замечание. В векторной форме теорема формулируется

так:

если

система

вектор

-

функций

решений

Yk (x)

colon(yk1, yk 2 ,..., yk n)

(k

1,2, n.)

. является. ,

фундаментальной системой решений в интервале (a,b)

однородного

уравнения

Y

A(x)Y ,

то вектор-функция

 

 

 

 

 

 

n

 

Y (x)

colon(y1, y2 ,..., yn ) ,

где

ym (x)

Ci yim

дает общее

 

 

 

 

 

 

i 1

 

решение

этого

уравнения

в

области

 

 

 

 

 

 

a x

b,

ym

(m 1,2,..., n) .

 

 

 

 

Доказательство. Совокупность функций (5.39) является

решением системы (5.23)

при

всех значениях

постоянных

C1, C2 ,..., Cn , поскольку представляет собой линейную

комбинацию решений системы.

Система (5.39) представляет собой линейную алгебраическую систему уравнений относительно C1 , C2 ,..., Cn . Ее определителем является вронскиан W (x) (5.32), который не равен нулю при любых x (a,b) , так как совокупность

функций (5.38) является по условию фундаментальной системой решений.

В соответствии с определением общего решения нормальной системой дифференциальных уравнений, учитывая вышеприведенные рассуждения, можно сделать вывод: совокупность функций (5.39) является общим решением системы (5.23) в области изменения переменных

(5.41).

*Покажем, что формула (5.39) содержит все решения системы (5.23). Для этого следует доказать, что при любых начальных условиях, наложенных на функции y1 ,..., yn

y1 y1(0) , y2 y2(0) ,..., yn yn(0) при x x0 ,

(5.42)

181

где (x0 , y1(0) , y2(0) ,..., yn(0) ) - произвольная фиксированная точка области (5.41), можно единственным образом найти значения постоянных C1(0) ,C2(0) ,..., Cn(0) таких, что полученное частное решение системы (5.23) с начальными условиями (5.42) будет

включаться

 

 

 

 

в

 

 

 

 

 

решение

 

 

 

 

 

 

 

 

 

(5.39)

при

C

C(0) ,C

2

C

(0)

,..., C

n

C(0) .

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

2

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Действительно, полагая в (5.39) x

 

x0 , учитывая (5.42)

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(0)

C ( y

 

)

0

C

2

( y

21

)

0

C

n

( y

n1

)

0

,

 

 

 

1

 

1 11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(0)

C ( y

 

)

 

C

 

( y

 

 

)

 

C

n

( y

2n

)

 

 

,

(5.39)

 

2

 

1

21

 

0

 

2

 

 

22

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

y(0)

C ( y

n1

)

0

C

2

( y

n2

)

0

C

( y

 

 

)

0

,

 

 

n

 

1

 

 

 

 

 

 

 

 

n

 

 

nn

 

 

 

где

 

( yij )0

 

yij (x0 ).

Полученная

 

 

система

линейных

неоднородных

уравнений

относительно

 

неизвестных

C1 , C2 ,..., Cn

 

 

 

имеет единственное решение, поскольку ее

определитель

равен W (x0 ) -

вронскиану

фундаментальной

системы функций (5.38) при x

x0 , которой отличен от нуля.

Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

C (0) (x

0

, y(0)

,..., y

(0) ), C

2

C (0) (x

0

, y(0)

,..., y (0) ),...,

1

1

 

1

n

2

1

 

n

(5.44)

 

 

C (0) (x

 

 

, y(0) ,..., y(0) ).

 

 

 

 

 

 

C

n

0

 

 

 

 

 

 

 

 

n

 

1

 

n

 

 

 

 

 

 

 

 

 

Подставляя эти значения постоянных в общее решение

(5.39) найдем решение задачи Коши (5.23), (5.42).

 

 

 

Таким образом, чтобы линейная комбинация n решений

однородной

 

 

 

линейной

системы (5.23)

с

произвольными

постоянными C1 , C2 ,..., Cn

давала

общее

решение этой

системы, необходимо и достаточно, чтобы эти решения были линейно независимы, т.е. чтобы они составляли фундаментальную систему решений.

182

*З а м е ч а н и е. Общее решение линейной однородной системы дифференциальных уравнений можно выразить через фундаментальную матрицу (5.36). Формулу (5.39) можно записать в виде

 

 

Y(x) Ф(x)C,

 

(5.45)

где

Y (x)

colon(y1 (x), y2 (x),..., yn (x)) ,

yk (x) - функции из

(5.39),

Ф(x) -фундаментальная

матрица

(5.36),

C

colon(C1, C2 ,..., Cn ) - постоянный вектор с произвольными

координатами.

 

 

 

Если положить в (5.45) x x0 , получим Y (x0 )

Ф(x0 )C .

Отсюда,

поскольку матрица Ф(x0 )

не вырождена, имеем

C

Ф 1(x0 )Y (x0 ) . Подставляя C в формулу (5.45), получаем

 

 

Y (x) Ф(x)Ф(x0 )Y (x0 ) .

(5.46)

Полученная формула дает решение начальной задачи (5.23), (5.42), т.е. решение задачи Коши для линейной однородной системы дифференциальных уравнений.

Матрица

K(x, x

)

Ф(x)Ф 1(x )

(5.47)

0

 

0

 

называется матрицей Коши. Используются также другие названия этой матрицы – импульсная матрица, а также

матрицант.

Таким образом решение задачи Коши (5.23), (5.42) можно записать с помощью этой матрицы в виде

 

 

Y (x) K (x, x0 )Y (x0 ) .

(5.48)

Пример. Найти общее решение линейной однородной

системы

 

 

 

 

dy

(x x 2 ) y (x3

x 2

x 1)z,

 

 

 

dx

 

 

 

 

 

dz

(1 x) y (x 2

x 1)z.

 

 

 

 

 

dx

 

 

 

 

 

Решение. Из вида системы ясно, что одно из его решений представляет собой многочлен. Поскольку в правой

183

части уравнений системы коэффициент при функции z(x) имеет степень на единицу большую, чем коэффициент при y(x) , то частное решение ищем в виде

 

y (x)

a b x c x2 ,

z (x)

a

2

 

b x .

 

 

 

1

 

 

1

1

1

 

 

1

 

 

 

2

 

 

Подставляя

y1(x) и z1 (x) в уравнения заданной системы

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

x ,

получим уравнения для определения коэффициентов:

 

c1

b2

 

0, c1

b1

 

a2

b2

0, b1

a1

a2

b2 0,

 

a1

c1

 

a2 2c1

 

0,

b1

a2 , a1

 

 

a2

b2

0.

 

Уравнения удовлетворяются, если b1

a2

0, a1

b2

c1.

Полагая a1

b2

c1

1,

находим одно решение системы:

 

 

 

 

 

 

y

 

x2

1,

z

 

x.

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

Второе решение системы будем искать в виде

 

 

 

 

 

y

2

( A B x)e x ,

z

2

( A B

2

x)e x .

 

 

 

 

 

1

 

1

 

 

 

2

 

 

 

 

 

 

Подставим эти функции в уравнения системы, разделим

обе части уравнений на e x , затем приравняем коэффициенты в правой и левой частях полученных уравнений при одинаковых степенях x . Получим уравнения:

A1 B1

A2 ,

A1

B2 A1

B2 ;

B1

A1

A2

B2 ,

B2

A1

B1

A2

B2 ;

0

A1

B1

A2

B2 ,

0

B1

A2

B2 ; 0

B1

A2 ,

0 B2 .

 

Уравнения удовлетворяются, если B2

A1

0 ,

A2

B1 1.

Поэтому, второе решение системы имеет вид

 

 

 

 

y2

xex ,

z2

e x .

 

 

 

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

W (x)

x 2

1 xex

e x

0 .

x

e x

 

 

 

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

y C y

1

C

2

y

2

C (x2

1) C

2

xex ,

1

 

 

1

 

 

184

z C z

1

C

2

z

2

C x C

e x .

1

 

 

1

2

 

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

Рассмотрим далее неоднородную линейную систему

(5.22)

dyk

n

 

ak i (x) yi

fk (x), (k 1,2,..., n) ,

dx

i 1

 

 

 

которую в векторной и операторной формах мы записывали следующим образом

 

 

 

dY

A(x)Y F,

 

L[Y ]F.

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

Установим наиболее важные свойства решений

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

 

 

 

 

 

 

Теорема

7

 

(принцип

 

суперпозиции).

 

Пусть

 

m

 

 

 

 

 

 

 

 

 

 

F (x)

i Fi (x),

где

i -

постоянные

 

числа,

 

i 1

 

 

 

 

 

 

 

 

 

 

F (x)

colon( f

(i) (x),..., f (i) (x)) ,

 

и

пусть

Y (x)

-

решение

i

1

 

 

n

 

 

 

i

 

 

уравнения Yi

A(x)Yi

Fi (x) . Тогда вектор-функция

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

Y (x)

 

iYi (x)

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

является решением уравнения (5.22).

Доказательство вытекает из линейности оператора L (5.23). Действительно,

 

d m

 

m

m

dY

 

m

L[Y ]

 

 

iYi A(x)

iYi

i

i

A(x)Yi

i Fi .

dx i 1

dx

 

 

i 1

i 1

 

i 1

Таким образом, Y (x) является решением уравнения (5.22).

Физический смысл принципа суперпозиции заключается в следующем. Пусть имеется физическая или техническая система, внутренние свойства которой описываются линейным оператором L . Пусть на вход такой системы подается сигнал Fk (x) (k 1,2,..., m) , а на выходе фиксируется

185

сигнал Yk (x)

(k 1,2,..., m) . Если

на вход такой системы

 

 

m

 

подать сигнал

F(x)

 

i Fi (x) , то на выходе из устройства

 

 

i 1

 

 

 

 

m

 

будет наблюдаться сигнал Y (x)

iYi (x) (рис.11).

 

 

 

i 1

 

 

 

 

 

 

Fk (x)

 

L

 

Yk (x)

1F1

 

 

 

 

 

 

 

 

 

 

2 F2

 

 

L

Y (x)

mFm

Рис. 11

Принцип суперпозиции имеет эмпирическое происхождение. Например, в теории тяготения согласно этому принципу гравитационное поле, возбуждаемое какой-либо массой, не зависит от наличия других масс. И кроме того, гравитационное поле, создаваемое несколькими телами, равно геометрической сумме гравитационных полей, возбуждаемых этими телами в отдельности.

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

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

Теорема 8. Если система линейных уравнений

dY

dx A(x)Y F1 (x) iF2 (x) ,

186

где i

-

 

мнимая

единица,

 

F1 (x)

colon(u1 ,..., un ) ,

F2 (x)

colon(v1 ,..., vn )

с

действительными

функциями aij (x) ,

ui (x) ,

vi (x) (i, j

1,2,..., n) имеет решение

Y (x)

Y1(x) Y2 (x),

Y (x)

colon(y(1)

,..., y(n) ) ,

 

 

Y (x) colon(y(1)

,..., y(n) ),

то

1

 

 

 

1

 

1

 

 

 

2

 

 

2

 

2

 

действительная

часть

решения

Y1 (x)

и

его мнимая часть

Y2 (x) соответственно являются решениями уравнений

 

 

 

dY1

A(x)Y1

F1 (x) ,

 

dY2

A(x)Y2

F2 (x) .

 

 

 

dx

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

По условию L[Y ]

F1

iF2 , где

L -

линейный

оператор

(5.23).

 

 

Докажем,

что

L[Y1 ]

F1 ,

L[Y2 ]

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

в силу линейности оператора

L ,

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L[Y1

iY2 ]

L[Y1] iL[Y2 ]

F1

iF2 .

 

 

Учитывая определение равенства комплексных чисел,

имеем L[Y1 ]

 

F1 ,

L[Y2 ]

F2 .

 

 

 

 

 

 

 

 

 

Теорема 9.

Если Y (n) (x)

является решением линейной

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L[Y]

 

F ,

 

 

 

 

(5.41)

а Y (0) (x) - решением соответствующей однородной системы

L[Y]

, то сумма Y (n) (x) Y (0) (x)

также будет решением

неоднородной системы (5.49).

 

 

 

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

По условию теоремы

L[Y (n) ]

F и

L[Y (0) ]

,

где

- нулевой

вектор.

Покажем,

что

L[Y (n)

Y (0) ]

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

используя линейность

оператора L , будем иметь:

 

 

 

 

L[Y (n) Y (0) ]

L[Y (n) ] L[Y (0) ] F

F.

 

Таким образом, теорема доказана.

Теорема 10 (об общем решении неоднородной системы).

Общее решение на отрезке x [a,b] неоднородной системы

187

(5.49) с непрерывными на [a,b]

коэффициентами aij (x) и

непрерывными правыми частями

fi (x) равно сумме общего

 

 

n

 

решения

Yoo (x)

CiYi , соответствующей неоднородной

 

 

i 1

 

системы ( Yi (x)

- линейно независимые решения однородной

системы)

и частного решения

Yчн (x) рассматриваемой

неоднородной системы.

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

 

 

 

dyk

fk (x, y1,..., yn )

(k

1,2,..., n),

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с

начальными

 

условиями

yk (x0 )

yk 0 ,

являются:

непрерывность всех

функций f k

в окрестности

начальных

данных

и

ограниченность

 

fk

 

.

Очевидно,

что для

 

yk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

нормальной

линейной системы

с непрерывными aij (x)

и

fi (x) эти условия выполняются при

 

x [a,b] .

Поэтому для

доказательства

сформулированной

 

теоремы

 

достаточно

показать,

что

подбором

постоянных

Ci

в

решении

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

CiYi

Yчн

можно удовлетворить произвольно заданным

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

начальным условиям Y (x0 )

Y0

colon(y10 , y20 ,..., yn0 ) .

 

Покажем, что векторное уравнение

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CiYi (x0 ) Yчн (x0 ) Y0

 

 

 

 

 

 

 

 

i

1

 

 

 

 

 

 

 

 

 

 

или

эквивалентная

система

алгебраических

линейных

по

Ci (i

1,..., n) уравнений

 

 

 

 

 

 

 

 

 

 

188

 

 

n

 

 

 

 

 

Ci y1i (x0 )

yчн1 (x0 )

y10 ,

 

 

i

1

 

 

 

 

 

n

 

 

 

 

 

Ci y2i (x0 )

yчн2 (x0 ) y20 ,

(5.50)

 

i

1

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

Ci yni (x0 )

yчнn (x0 )

yn0

 

 

i

1

 

 

 

имеет

решение C(0) ,..., C(0) , при

любых

правых частях

 

 

1

n

 

 

yk 0

yчнк (x0 )

системы (5.50). Но определитель этой системы

равен определителю Вронского в точке x

x0 . Для линейно

независимых

решений

Y1,Y2 ,..., Yn

соответствующей

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

отличен

от нуля. Следовательно,

система (5.50) имеет решение C1(0) ,C2(0) ,..., Cn(0) при любых

правых частях. Теорема доказана.

Метод вариации произвольных постоянных (метод Лагранжа).

В том случае, когда удалось найти общее решение однородной системы линейных уравнений, соответствующей заданной неоднородной системы, для нахождения частного решения последней (а вместе с тем и общего ее решения) применяется метод Лагранжа.

Теорема 11. Пусть на сегменте [a,b] матрица A(x) и вектор F(x) непрерывны (т.е. матричные элементы aij (x) и координатные функции fi (x) непрерывны при x [a,b] ).

Пусть известна

фундаментальная система решений для

соответствующей

однородной

системы

уравнений

Y (x) A(x)Y (x) .

Тогда общее решение неоднородной

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

Y (x) A(x)Y(x)

F(x) может быть

найдено с помощью квадратур.

 

 

189

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