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

Дифференциальные уравнения. Учебное пособие-3

.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
1 Мб
Скачать

творяющая однородным начальным условиям и, значит, равная нулю в некоторой окрестности точки x0 . Следовательно, решение СЛДУ –

тривиальное решение, т. е. Yi линейно зависимы и определитель Вронского равен нулю в любой точке отрезка [a, b].

Определение. Совокупность решений Y1, Y2 , , Yn системы (9)

называется фундаментальнойсистемойрешенийна [a, b], еслиопре-

делитель Вронского не обращается в нуль ни в одной точке [a, b].

Определение. Если для однородной системы ДУ (9) начальные условия определены следующим образом: Y (x0 ) E, то система реше-

ний Y1, Y2 , , Yn называется нормальной фундаментальной систе-

мой решений.

Замечание. Если Y1, Y2 , , Yn – фундаментальная система или нормальная фундаментальная система, то линейная комбинация

n

CiYi (x) – общее решение (9).

i 1

n

Теорема 2. Линейная комбинация CiYi (x) линейно независимых

i 1

решений Yi , i 1, ..., n, однороднойсистемыДУ(9) снепрерывнымина отрезке [a, b] коэффициентами aij (x) будет общим решением (9) на

этом же отрезке.

Доказательство. Так как коэффициенты aij (x) непрерывны на

[a, b], то система удовлетворяет условиям теоремы о существовании и

единственности. Следовательно, для доказательства теоремы достаточно показать, что подбором постоянных Сi можно удовлетворить не-

которому произвольно выбранному начальному условию (7), т. е.

n

можно удовлетворить векторному уравнениюCiYi (x0 ) Y0 . Так как

i 1

n

CiYi (x) – общее решение (9), то система разрешима относительно Сi ,

i 1

посколькувсе Yi линейнонезависимыи W (x0 ) 0. Однозначноопределяем Сi ; так как Сi линейно независимы, то W (x) 0 x [a, b] .

21

Теорема 3. Если Y – решение неоднородной системы ДУ (8), а Y

решение однородной системы ДУ (9), то Y + Y будет решением неоднородной системы (8).

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

L[Y Y ] L[Y ] L[Y ] F.

Теорема 4. Общее решение неоднородной системы ДУ (8) на отрезке [a, b] с непрерывными на этом отрезке коэффициентами aij (x)

и правыми частями fi (x) равно сумме общего решения соответствую-

щей однородной системы (9) и частного решения Y неоднородной системы (8).

Доказательство. Так как условия теоремы о существовании и един-

n

ственности выполнены, то остается доказать, что Y (x) CiYi Y бу-

i 1

дет удовлетворять произвольно заданным значениям (7), т. е.

n

 

Y (x0 ) Y0 CiYi (x0 ) Y (x0 ) .

(11)

i 1

Для системы (11) всегда можно определитьзначения Сi . Это можно сделать, так как Yi , i = 1, …, n, – фундаментальная система решений.

m

Теорема 5 (принцип суперпозиции). Решение СЛДУ L[Y ] Fi

 

 

 

m

 

i 1

 

 

 

 

 

может быть представлено в виде Yi , где каждое Yi

удовлетворяет

 

 

 

i 1

 

 

 

 

 

m

m

m

уравнению L[Yi ] Fi

i

1,m

, тогда L[Yi ] L[Yi ]

Fi .

 

 

 

i 1

i 1

i 1

Замечание. Принцип суперпозиций можно распространить на случай m при условии, что ряд, составленный из Yi , равномерно схо-

дится.

22

ОДНОРОДНЫЕ СИСТЕМЫ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ

Рассмотрим однородную систему дифференциальных уравнений (ОСДУ)

dY

AY

(12)

dx

 

 

где A – матрица с постоянными коэффициентами (n n); Y – неизвест-

ная вектор-функция (n). Общее решение системы (12) имеет вид Y (x) exp(Ax)C, где C – вектор произвольных постоянных [3].

Согласно теореме о приведении матрицы к квазидиагональному

виду [3] матрица A может быть представлена в виде A SJS 1 , где J – жордановаформаматрицыA; S – матрица, составленнаяизсобственных и присоединенных векторов матрицы A. Тогда систему (12) можно за-

писать в виде

dY

SJS 1Y, или

d(S 1Y )

JS 1Y. Сделаем замену пе-

dx

dx

ременных Z S 1Y и запишем ОЛСДУ покомпонентно для случая, когда J – одна жорданова клетка размерностью (n n) :

dzdx1 z1 z2 ; dzdx2 z2 z3 ;

.... ....

dzn 1 zn 1 zn ; dx

 

dzn

zn .

 

 

 

dx

 

 

 

 

 

 

Решая

эту систему,

получаем: zn (x) Cn exp( x) ,

zn 1(x)

(Cn x Cn 1)exp ( x) , ...,

z1(x) Cn xn 1 Cn 1xn 2 C1

exp ( x).

Тогда общее решение ОЛСДУ: Y (x) SZ(x) .

23

Если все собственные значения матрицы A действительны и раз-

n

 

личны, то общее решение ОЛСДУ Y (x) Ciei exp( i x) , где Ci

– про-

i 1

 

извольные постоянные; ei – собственные векторы матрицы A; i

– соб-

ственные значения матрицы A.

 

ЕслиматрицаA вещественнаяисреди еесобственных значенийесть комплексно-сопряженные, то частными решениями ОЛСДУ будут не только решения Yj (x) ej exp(( p iq)x) , но и решения вида

Y1 j (x) Re{ej exp(( p iq)x)} и Y2 j (x) Im{ej exp(( p iq)x)} [1].

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

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

Рассмотрим неоднородную линейную систему ДУ (НЛСДУ):

 

 

 

dY

AY F(x) .

(13)

 

 

 

dx

 

 

 

 

 

Пусть Y (x) [Y

Y

Y ]

(Y – частные решения ОЛСДУ) – фун-

1

2

 

n

i

 

даментальнаясистемарешенийоднороднойсистемыДУ(12). Будемискать решение неоднородной системы (13) в виде Y (x) Y (x)C(x) , где C(x) – неизвестная вектор-функция. Подставляя решение в си-

стему (13), получаем соотношения для C(x)

dC(x)

Y

1F(x) , тогда

dx

C(x) Y 1(x)F(x)dx C;

C – вектор произвольных

постоянных.

Тогда общее решение неоднородной системы ДУ будет иметь вид

Y (x) Y (x) Y 1(x)F (x)dx C .

24

МЕТОД КОШИ

Рассмотрим однородную систему ДУ (12) и сопряженную к ней систему:

 

 

dZ

 

 

 

 

A Z.

(14)

 

 

dx

Пусть 1(x), , n (x)

– нормальная фундаментальная система ре-

шений СЛДУ (12), 1(x),

, n (x) , – нормальная фундаментальная

системарешенийСЛДУ(14). Дляэтихсистемрешенийвыполненыусловия взаимной ортогональности: ( i (x), j (x)) 0, i j [1]. Это условие

выполняется для всех x [a, b] . Будем искать решение неоднородной

 

 

 

n

 

 

 

 

 

 

системы ДУ (13) в виде Y (x) ui (x) i (x) , где ui (x)

– неизвестные

 

 

i 1

 

 

 

 

 

 

функции.

Подставляя решение

 

в систему

(13)

и

учитывая, что

i (x), i 1,

 

 

 

 

 

 

n

 

 

 

, n, – решения(12), получаем ui (x) i (x) F(x) . Умно-

 

 

 

 

 

 

 

i 1

 

 

 

жая последнее соотношение скалярно на

j (x) , получаем уравнения

для определения функций ui (x) :

 

 

 

 

 

 

 

 

 

u j (x)

(F, j )

dx C j ,

j 1, , n .

 

 

( j , j )

 

 

 

 

 

 

 

 

 

 

Тогда общее решение НЛСДУ (13) будет иметь вид

 

 

 

 

n

 

(F, j )

 

 

 

 

 

Y (x) j (x)

dx

C j

.

 

 

 

 

 

j 1

 

( j , j )

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание. Метод Коши применим для СЛНДУ с матрицей А, приводящейся к диагональному виду.

25

МАТРИЧНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ

Пусть дано матричное дифференциальное уравнение

dY (x)

dx

где

y11(x)

Y (x) y21(x)

yn1(x)

A(x)Y (x) ,

y12 (x) y1n (x)

y22 (x)

 

 

y2n (x)

 

 

 

yn2 (x)

 

 

ynn (x)

(15)

(16)

с начальными значениями

 

Y (x0 ) Y0 .

(17)

Теорема. Если матрица A(x) непрерывна на [a, b], а определитель матрицы Y0 0 , то на отрезке [a, b] существует единственное решение Y (x) уравнения (15) и определитель Вронского этого решения не обращается в нуль ни в одной точке [a, b].

Будем рассматривать одновременно с системой (15) систему вида

dZ(x)

Z(x)A(x) .

(18)

dx

 

 

Определение. СЛДУ (18) называется сопряженной к СЛНДУ (15).

Теорема. Пусть Y (x) – решение (15), а матрица A(x) непрерывна

на [a, b], тогда Y 1(x) существует на [a, b] и является решением си-

стемы (16).

Доказательство: по предположению теоремы решение Y (x) на отрезке [a, b] существует и определитель этого решения не равен нулю в любой точке [a, b], следовательно, существует обратная матри-

26

ца Y 1(x) . А так как Y (x)Y 1(x) E, то, продифференцировав это соотношение, получаем:

dxd Y (x)Y 1(x) 0,

Y (x) dY 1(x) dY (x) Y 1(x) 0, dx dx

где Y (x) – решение уравнения (15), следовательно,

dY

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

dY 1(x)

A(x)Y (x)Y 1(x) 0,

dx

dx

 

 

dY 1(x) Y 1(x)A(x). dx

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

dY (x)

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

(19)

dx

 

 

Будем искать решение этой системы методом вариации постоянной. Обозначим Y (x) – решение (15). Будем искатьрешение системы (19)

в виде

 

 

 

 

Y (x) Y (x)C(x) ,

(20)

где C(x) – неизвестная вектор-функция.

 

 

Подставим (20) в (19):

 

 

 

 

 

 

dY

 

 

dC(x)

 

 

 

 

 

dx C(x) Y (x)

dx

A(x)Y (x)C(x) F(x);

 

dC(x)

 

dC

 

 

 

 

x

Y

F(x),

Y 1(x)F(x),

C(x)

Y 1( )F( )d C.

dx

 

 

 

dx

 

 

 

x0

 

 

 

 

 

 

 

 

27

Тогда общее решение неоднородной системы будет иметь вид

x

Y (x) Y (x) Y 1

x0

( )F( )d C . (21)

Теорема. Общее решение неоднородной системы (19) можно представить в виде

x

Y (x)Y 1

 

 

 

Y (x) (x) Y (x)

( )F( )d

,

(22)

 

 

 

 

 

x0

 

 

 

 

где (x) – общее решение соответствующей однородной системы.

Замечание. ПустьматрицаАпостояннаиначальныеусловияимеют

вид

dY

AY (x),

 

(23)

dx

Y (0) E.

Покажем, что решение CЛНДУ (23) – Y (x) удовлетворяет функци-

ональному уравнению

 

 

 

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

 

(24)

При любом фиксированном

матрицы Y (x )

и Y (x)Y ( )

будут

решением (23). При x 0 эти матрицы будут совпадать в силу начальных условий и, следовательно, по теореме о существовании и единственности они будут совпадать для любых x . Если воспользоваться (24), то можно записать(сделавзамену x x иумноживсоот-

ношение (24) справа на Y 1( )) :

Y (x)Y 1( ) Y (x )Y ( )Y 1( ) Y (x ) .

(25)

Подставляем (25) в(21) иполучаемвидобщегорешения НСДУспостоянной матрицей А:

x

 

Y (x) (x) Y (x )F( )d .

(26)

x0

 

28

ПРИМЕР РЕШЕНИЯ СНЛДУ

1. Решить неоднородную систему дифференциальных уравнений методом вариации постоянной:

 

 

dY

 

 

 

0

3

3

 

0

 

 

 

AY F(x) , где

A

1

8

6

,

F (x) 0

.

 

 

 

 

 

dx

 

 

 

2

14

 

 

 

 

 

 

 

 

 

10

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Собственные значения матрицы A: 1

0, 2,3

1. Собственный век-

тор для

e (2, 1,

1)T

, собственный и присоединенный векторы

 

 

1

 

 

1

4)T ,

e (3, 0, 2)T . Тогда

 

 

для

2

 

e

(3, 3,

матрица перехода

 

 

2

 

3

 

 

 

 

 

 

S

J

dY

dx

2

3

3

 

 

 

 

1

3

 

и жорданова

форма

матрицы A имеет вид

 

0

 

1

4

2

 

 

 

 

 

 

 

 

 

 

0

0

0

 

 

1

 

 

1

 

Сделаем замену

Z S

0

1 .

Y в однородной СЛДУ

0

0

1

 

 

 

 

 

 

 

 

 

 

AY и получим систему линейных дифференциальных уравнений:

dz1

0,

dz2

z

2

z ,

dz3

z .

 

 

 

dx

dx

3

dx

3

 

 

 

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

z1(x) C1,

z2 (x) (C3x C2 )exp( x),

z3 (x) C3 exp( x) .

 

2

3

3

 

 

 

C1

 

 

Y (x) SZ (x)

1

3

0

(C

x C

2

)exp ( x)

 

 

 

 

1

4

2

 

3

 

 

 

 

 

 

 

 

 

C exp ( x)

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

29

 

2

 

3e x

 

3(x 1)e x

C1

 

 

 

 

3e x

 

 

3xe x

 

 

 

1

 

 

 

 

 

C

2

.

 

 

 

 

x

 

 

x

 

 

 

 

1

4e

( 4x 2)e

 

 

 

 

 

 

 

C3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

2

 

3e x

3(x 1)e x

 

 

 

 

 

3e x

3xe x

 

 

 

Y (x)

1

 

 

 

.

 

 

 

 

 

4e x

( 4x 2)e x

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

системы

 

ДУ: Y (x) Y (x)C(x) ,

C(x) Y 1(x)F(x)dx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 2

 

 

5x

7

 

 

1

(x)F(x)

 

 

x

 

x

Y

 

3

 

 

e

 

3

 

 

e

 

 

 

 

1

 

x

 

5

 

x

 

 

 

 

 

e

 

e

 

 

 

 

 

 

3

 

 

 

 

3

 

 

 

Тогда

3

 

 

 

 

 

0

 

 

 

 

 

(1 x)ex

0

 

 

 

 

 

 

 

 

 

 

 

 

ex

x

 

 

 

 

 

 

 

 

 

 

 

 

3x

 

 

 

 

(1

x)ex .

 

xex

 

 

 

C1(x) 3xdx 32x2 C1 ,

C2 (x) (x x2 )exdx (x2 2x 2)ex C2 , C3(x) xexdx (x 1)ex C3 .

30

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