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

Математический_анализ_2_семестр

.pdf
Скачиваний:
23
Добавлен:
24.05.2015
Размер:
2.05 Mб
Скачать

Теорема. Пусть задана система n дифференциальных уравнений первого порядка

y

 

 

 

f

1

x, y ,... y

n

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

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

 

 

 

 

y

 

 

 

f

n

x,

y ,... y

n

 

 

 

 

 

n

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

Обозначим

 

 

 

 

 

 

 

f k

 

 

 

F2 f1 f1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

k 1 yk

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

n

F

 

 

 

 

 

 

 

 

F3

 

 

 

 

 

2

 

 

 

2

 

 

f k

 

 

 

 

 

 

 

 

 

yk

 

 

 

 

 

 

 

 

 

 

x

 

 

 

k 1

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

F

 

 

n

F

2

 

Fn 1

 

 

 

 

 

 

n 2

 

 

 

 

n

f k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

k 1

 

 

yk

Потребуем, чтобы

 

функция Fn 1 была бы дифференцируемой по

совокупности переменных. Потребуем, чтобы определитель

 

 

 

 

f1

 

...

 

 

f1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

yn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F2

...

 

F2

 

0

 

 

y2

 

yn

 

 

 

 

 

 

 

 

 

 

 

 

... ... ...

 

 

 

 

 

 

Fn 1

...

Fn 1

 

 

 

 

 

 

 

y2

 

yn

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда система n дифференциальных уравнений эквивалентна одному дифференциальному уравнению n-ого порядка.

Доказательство. Метод доказательства называется методом исключения

переменных и применяется на практике при сведении системы к одному уравнению. Продифференцируем Fn 1 :

Fn

 

Fn 1 Fn 1

yk

 

 

 

n

 

 

 

 

x

k 1

yk

 

1) Построим алгоритм метода исключения.

Пусть y1 , ... yn - решения системы ( y1 f1 ,...yn fn ), тогда уравнения системы F2 ..., Fn 1 ... представляют собой тождества

 

 

f1

n

f1

 

 

 

F2

 

 

 

 

 

yk

y1

 

x

yk

 

 

 

k 1

 

 

 

 

 

 

n

 

 

yk y1

 

F3

 

F2

F2

'''

 

 

x

k 1 yk

 

 

 

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

Fn 1

Fn 2 Fn 2

yk y1 n 1

 

 

n

 

 

Fn

x

k 1 yk

yk y1 n

Fn 1 Fn 1

 

 

n

 

 

 

x

k 1 yk

 

 

71

Получены выражения производных

y1 f1 x, y1 , y2 ,... yn , y1'' F2 x, y1 , y2 ,... yn , y1''' F3 x, y1 , y2 ,... yn ,

...

 

y n 1

F

1

x, y ,

y

2

,...

y

n

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

n

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из этих уравнений можно выразить

y

2

 

,... y

n

через

y ,... y n 1

,

так как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

определитель системы этих уравнений 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

y

2

,... y

n

через

y

 

 

,... y n 1

 

в последнее уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

n

 

 

y1 ,

y2

 

 

 

 

 

n 1

 

 

 

 

 

 

n 1

 

 

 

Так

как

y1 , ... yn -

 

решения

y1

Fn x,

y1

... y1

 

,... yn y1

...y1

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

системы

y

 

f

1

,...y

 

f

n

,

 

то

они

являются

и

решениями

полученного

 

 

1

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уравнения.

 

Следовательно,

 

система

y

 

f

1

,...y

 

f

n

сведена к

 

одному

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

n

 

 

 

 

 

 

уравнению n-ого порядка.

2)Покажем эквивалентность решений. Предположим, что y1 , ... yn - решения полученного уравнения, покажем, что y1 , ... yn - решения системы.

 

 

 

 

 

 

 

 

 

 

 

 

f1

 

f1

 

n

f1

 

 

 

 

 

 

 

y1

f1

 

,

 

y1

 

 

 

 

 

 

 

 

 

 

yk .

Обозначим

y1

 

F2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

y

y1

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 2

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

'''

 

F2

 

F2

 

 

n

F2

 

 

 

 

 

 

 

 

 

 

 

 

'

 

 

 

y1

 

 

 

 

 

 

 

 

 

yk .

 

 

 

Обозначим

 

y1 F3 ,

и

 

т.д.

x

y

y1

y

 

 

 

 

 

 

 

 

 

 

 

k 2

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

F

 

 

F

 

 

n

 

F

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

n 1

 

 

n 1

 

 

n 1

 

 

 

 

 

 

 

 

Fn .

 

 

 

y1

x

y

y1

y

 

yk . Обозначим y1

 

 

 

 

 

 

 

 

k 2

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Приравниваем полученные здесь функции F2 , F3 ,...Fn введенным ранее, сокращая первые и вторые слагаемые, получаем систему уравнений

n

f1

 

 

 

 

 

 

 

 

 

yk

 

 

f k

 

0

yk

 

k 2

 

 

 

 

 

 

n

F2

 

 

 

 

 

 

 

 

 

 

yk

 

 

f k

 

0

 

yk

 

k 2

 

 

 

 

 

 

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

n

Fn 1

 

 

 

 

 

 

 

 

 

 

yk

 

f k

0 .

 

yk

 

 

k 2

 

 

 

 

 

 

 

Определитель этой системы равен единственного решения системы имеем y2 эквивалентны. Теорема доказана.

x 2x y

Пример. y x 2 y

0 , следовательно, в качествеf2 ,...yn fn . Поэтому решения

72

x 2x y 2x x 2 y 2x x 2 x 2x 4x 3x ,

x 4x 3x 0, k

2 4k 3 0,

k 1, k

2

3,

 

 

 

 

 

x t C et

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

C

e3t ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y t x t 2x t C et 3C

e3t 2C et

2C

e3t C et

C

e3t

 

 

 

 

 

 

 

 

1

 

2

 

 

1

 

 

2

 

 

 

 

1

2

 

 

 

Функция

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

x, C называется общим решением системы, если

1)

 

 

 

 

 

 

 

 

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

 

 

 

 

 

для любого C

y

x, C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)

для

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

начальных

условий

 

x0 , y0

найдется

C0 , что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y0 x0 , C0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если зафиксировать

C

в общем решении, получим частное решение

системы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача Коши.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найти

решение

системы

 

 

 

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

заданным

y

f x, y ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

y x0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема Коши о существовании и единственности решения задачи Коши

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть функция f x,

y непрерывна по совокупности переменных. Пусть

существуют и непрерывны частные производные

 

 

f k

,

k 1,...n,

s 1,...n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ys

 

 

 

 

 

Тогда существует и единственно решение задачи Коши.

Первые интегралы.

Пусть выполнены условия теоремы Коши. Рассмотрим решение задачи

 

 

 

 

Коши y x, x0 при заданных начальных условиях

y x0

y0 . По теореме Коши

 

 

 

 

оно существует и единственно. Это решение y x, x0 , y0 можно представить

 

 

 

 

себе как некоторую интегральную кривую, соединяющие точки x0 , y0

, x, y .

 

 

 

 

Если в качестве начальных условий выбрать y x

y , то по теореме Коши

через эту точку проходит та же единственная интегральная кривая, ее уравнение

 

 

 

 

 

 

 

 

 

 

можно записать в виде

y0

 

x0 , x, y .

Зафиксируем

x0 , обозначим C y0

,

получим соотношение

 

 

 

 

 

общий

интеграл системы

 

 

x, y

C

дифференциальных уравнений (векторное соотношение). Первый интеграл системы дифференциальных уравнений – скалярная составляющая общего интеграла. Общий интеграл системы дифференциальных уравнений – векторная функция, сохраняющая свое значение на решениях системы. Первый интеграл системы дифференциальных уравнений – скалярная функция, сохраняющая свое значение на решениях системы.

Знание одного первого интеграла позволяет понизить порядок системы на единицу. Знание общего интеграла дает общее решение системы, если только

 

 

 

 

можно разрешить уравнение

x, y

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

y .

73

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

 

f k .

 

n

 

 

x

k 1 yk

 

Скалярная функция x, y1,... yn является первым интегралом, если

 

f k 0 .

 

n

 

 

x

k 1 yk

 

Симметричная форма записи системы.

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

 

y

 

 

dy1

 

f

 

 

x, y ,... y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

n

 

 

 

 

 

 

 

 

1

 

 

 

 

dx

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

dyn

f

 

 

x, y ,... y

 

 

 

 

 

 

 

 

 

 

n

 

n

n

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и запишем эти уравнения в симметричном виде

 

 

 

dx

 

 

 

 

 

 

 

 

dy1

 

 

 

 

...

 

 

dyn

 

 

.

 

 

 

1

 

f

1

x, y ,... y

n

 

f

n

x, y ,... y

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

Или,

 

 

 

заменяя

 

 

 

 

переменные

и

правые

части

xxn 1 , y1 x1 ,... yn xn , f1 X1 ,... f n X n , 1 X n 1 ,

получим симметричную форму записи системы

 

 

dx1

 

 

...

 

 

dxn 1

 

 

.

X

1

x ,... x

n 1

 

X

n 1

x ,... x

n 1

 

 

1

 

 

 

1

 

 

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

Пример. x

 

 

x

 

 

 

y

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 2 y 2

 

 

x 2

y 2

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

x 2

y 2

 

 

 

dy

 

 

 

y

 

 

 

 

 

 

dx

1

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

y Cx

 

 

 

 

 

 

,

xdx

 

 

dy

 

 

 

 

y

 

 

dx

x

 

dt

x 1 C 2

 

1 C 2

 

dt

 

 

x 2

y 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 2

 

 

 

 

1

t C ,

 

 

y 2

 

 

C 2

t C C 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

1 C 2

 

 

 

1

 

2

 

 

1 C 2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Автономные системы и свойства их решений.

Система называется автономной, если в ее правую часть не входит явно

 

 

 

 

.

независимая переменная: y

f

y

Решение автономной системы можно рассматривать в пространстве координат y1 ,... yn , которое принято называть фазовым пространством. Проекция интегральной кривой на это пространство называется фазовой

74

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

независимую переменную

yn 1

x и дополнительное уравнение

 

x

 

1.

yn 1

 

Фазовое пространство такой системы принято называть расширенным фазовым пространством.

Свойства решений автономных систем.

1)

Если

 

 

 

 

 

 

 

 

 

 

 

x C тоже решение.

y

x - решение системы,

то и y

 

 

x C

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

d x C

 

 

 

 

 

 

 

 

 

 

 

 

 

f

x C .

 

 

 

 

 

dx

 

 

d x C

 

 

 

Следствие.

 

Фазовая

траектория

 

 

- это та же фазовая

 

y

x C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

траектория, что и y x .

 

 

 

 

В самом деле, любая точка x C,

первой фазовой траектории является

y

 

 

 

 

 

 

второй фазовой траектории и наоборот.

точкой x, y

2) Две

фазовых траектории либо не имеют общих точек, либо

 

 

совпадают.

 

 

 

 

 

 

 

 

Пусть две различных фазовых траектории x , x имеют общую точку

x1 x2 . Рассмотрим решение v x x x2 x1 .

v x1 x1 x2

x1 x2 x1 .

Следовательно, по теореме Коши

v x

x .

 

Но

v x

- это траектория x ,

 

сдвинутая на x2 x1 по

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

Следствие. Множество фазовых траекторий автономной системы в фазовом пространстве представляет собой совокупность непересекающихся кривых.

Точка a называется точкой покоя (точкой равновесия) автономной

системы, если

 

 

 

 

 

 

 

 

 

 

f

 

a 0 .

 

 

 

 

 

 

 

 

 

 

- точка покоя, то

 

 

- решение системы.

3) Если точка a

y

a

В самом деле,

 

 

 

 

 

 

.

 

 

 

y

a 0

f

a

 

 

 

4)Любая фазовая траектория автономной системы есть траектория одного из трех типов:

-гладкая, не самопересекающаяся кривая,

-замкнутая гладкая кривая,

-точка покоя.

Фазовый поток.

75

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим

решение

задачи

 

Коши

 

автономной

системы

y x, y0 .

Определим фазовый поток

как оператор

g x сдвига (по

аргументу

x ) по

 

 

 

 

 

 

 

.

 

 

 

фазовым траекториям системы g x y

0

= y x, y

0

 

 

 

Рассмотрим

некоторую

область

D фазового пространства

(фазовым)

объемом V0 . Фазовый поток переводит эту область в область Dx объемом Vx .

 

 

 

 

 

 

 

 

 

 

dV

x

 

 

 

Справедлива теорема Лиувилля

 

divf

y

d y .

dx

 

 

 

 

 

 

 

 

 

 

Dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь мерой y в фазовом пространстве может служить фазовый объем

 

 

f

 

 

f

 

 

 

 

 

 

 

 

 

 

1

 

n

 

f

 

 

 

 

 

 

V , divf

y

 

...

 

Tr

 

(дивергенция векторного поля правых частей

 

 

 

 

 

 

y1

 

yn

 

y

 

 

 

 

 

 

системы или след матрицы Якоби). Левая часть этой формулы представляет собой изменение фазового объема в единицу «времени» – аргумента, т.е. известный из теории поля поток векторного поля правых частей системы – фазовых скоростей. Приведенная формула аналогична формуле Остроградского

– Гаусса в теории поля.

 

 

 

 

 

 

 

 

 

 

 

 

Если divf

y 0

, то Vx

const .

 

 

 

 

 

 

 

 

 

 

 

 

x ... ann x ,

 

 

 

 

Если

f y A x y , то divf

y TrA x a11

что

 

дает

 

 

 

 

фазового объема V x C exp

 

 

 

 

 

формулу

для

определения

 

 

 

,

что

TrA x dx

 

 

 

 

 

 

 

 

Dx

 

 

 

 

совпадает с формулой Остроградского – Лиувилля определителя Вронского для линейных автономных систем. Поэтому определитель Вронского имеет смысл фазового объема (определитель всегда имеет смысл некоторого объема, вспомним хотя бы смысл смешанного произведения векторов).

Лекция 21. Системы линейных дифференциальных уравнений.

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

 

 

 

 

x .

y

A x y

f

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

записать в виде

 

 

 

 

 

 

 

 

y

A x y .

 

 

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

Теоремы о свойствах решений однородной и неоднородной системы.

 

 

 

 

 

 

 

 

Если yo1

, yo 2

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

yo1

y02

, yo1

, y02

-

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

76

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если

y0,

yн

- решения однородной и неоднородной систем,

то

y0

yн -

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если

yн1

, yн 2 -

 

решения неоднородной

системы,

то

y

н1 yн 2 -

решение

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

 

 

 

 

 

 

 

 

 

 

 

 

 

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

y02

A x yo1

A x y02

A x y01 y02

,

 

 

 

yo1 yo2

 

y01

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y01

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y01 A x y01 A x y01

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

y0

yн

 

yo

yн

A x y0 A x yн

f

x A x y0 yн f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yн1

yн2

 

yн1

yн2

A x yн1

f x A x yн2 f x

A x yн1

yн2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема. Множество решений линейной однородной системы есть линейное пространство.

Из теорем о свойствах решений видно, что операции сложения и умножения на число на решениях однородной системы определены корректно.

Легко проверяется ассоциативность по сложению, существования «нуля»

 

существование

«противоположного

– тривиального решения y x 0 ,

 

по сложению.

Отсюда следует, что

элемента» y x , коммутативность

решения однородной системы образуют коммутативную группу по сложению (абелев модуль) (4 аксиомы линейного пространства). Существует единица – число, справедлива ассоциативность по умножению на число (еще 2 аксиомы).

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

 

 

 

x называются линейно независимыми, если

Функции y1 x ,...

yn

1 y1 x ... n yn x 0

1 0,... n 0 .

 

 

 

 

x называются линейно зависимыми, если

Функции y1 x ,...

yn

1 ,.. k

 

 

 

 

 

0,... n : 1 y1

... n yn 0 .

 

 

 

 

 

y11 x ...

yn1 x

Введем

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

Вронского W ... ...

... , по столбцам

 

 

 

 

y1n x ...

ynn x

 

y

x

 

y

 

x

 

11

 

 

 

n1

 

 

 

которого расположены векторы y1

x ...

,... yn

x

...

 

, введем также

 

y

x

 

y

nn

x

 

 

1n

 

 

 

 

 

 

 

 

y11 x ...

yn1

x

 

матрицу Y x

...

...

 

...

,

W x detY x .

 

 

 

 

 

 

 

 

 

 

y

x ...

y

nn

x

 

 

1n

 

 

 

 

 

 

77

 

 

 

x линейно зависимы, то W x 0 .

Теорема. Если функции

y1

x ,... yn

Доказательство. Так как функции линейно зависимы, то одна из них линейно выражается (тождественно) через остальные, поэтому соответствующий столбец определителя Вронского линейно выражается через остальные. Тогда по свойству определителя W x 0 .

 

 

 

 

 

 

 

 

Теорема. Пусть

y1 x ,... yn x

-

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

x0

 

 

 

 

 

x линейно зависимы.

: W x0 0 , тогда решения y1

x ,... yn

 

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

W x0 0 ,

то его столбцы в x0 линейно зависимы,

 

 

x0 ...

 

x0 0 .

т.е.

1 ,.. k 0,... n : 1 y1

n yn

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим решение

y x

1 y1

... n yn (с теми же коэффициентами).

y x - решение однородной системы как линейная комбинация решений

однородной системы (теоремы о свойствах решений). Начальные условия для

этого решения в точке

x0 ,

как показано выше, нулевые.

Но есть решение

однородной

системы

(тривиальное

решение

 

 

имеющее

те же

y x 0 ),

 

 

 

 

 

 

 

 

 

 

 

 

начальные условия. Следовательно,

 

 

 

 

 

 

и есть

по теореме Коши решение y x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,

k 0 ,

 

 

тривиальное решение. Тогда

y x 1 y1 ... n yn

следовательно,

 

x ,...

 

 

 

 

 

 

 

 

 

 

решения y1

yn x линейно зависимы.

 

 

 

 

 

 

Следствие.

Равенство

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

Вронского нулю

для

решений

однородной системы хотя бы в одной точке – критерий линейной зависимости решений, отличие определителя Вронского от нуля для решений однородной системы хотя бы в одной точке – критерий линейной независимости решений.

Доказательство. Пусть x0 : W x0

 

 

 

x

0

, тогда решения y1

x ,... yn

 

 

x линейно зависимы, то W x 0

линейно зависимы. Если решения y1

x ,... yn

по теореме о равенстве определителя Вронского нулю для системы линейно зависимых функций. Заметим, что тогда x0 : W x0 0 W x 0 .

 

x0 : W x0 0 ,

 

 

 

Пусть

если решения y1 x ,... yn x линейно зависимы, то

W x 0

(противоречие).

Пусть

решения

линейно независимы. Если

 

 

 

 

 

x0 : W x0 0 , тогда решения y1 x ,... yn x линейно зависимы (противоречие).

Теорема. Размерность пространства решений однородной системы равна

n.

 

 

 

 

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

доказать 1) существуют n линейно независимых

решений однородной системы, 2) любое решение однородной системы линейно выражается через эти линейно независимые решения.

1) В любой точке x0, y10, ... yn0 для однородной системы выполнены условия

теоремы Коши, следовательно, через любую такую точку пройдет единственная интегральная кривая – график решения однородной системы. Зададим такие точки – начальные условия, которые по теореме Коши

определят

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решения

 

 

 

 

 

 

1

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x ,

 

x0

 

 

0

 

 

x ,

 

x0

 

 

1

 

 

x ,

yn x0

 

0

 

 

y1

y1

 

 

,

y2

y2

 

 

...yn

 

...

.

 

 

 

 

 

...

 

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

78

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0 ...

0

 

 

Эти решения линейно независимы, так как

W

x

0

 

0

1 ...

0

1 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0 ...

1

 

 

Существование n линейно независимых решений однородной системы

доказано.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) Рассмотрим произвольное решение однородной системы

y x . В точке x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вектор

y x0 разлагается по естественному базису

 

 

 

 

 

 

 

1

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

1

 

 

 

0

 

 

 

 

 

 

 

 

 

 

y1 x0

 

 

...

, y2 x0

 

 

 

...yn

x0

 

 

...

.Поэтому

y x0

C1 y1

x0 ... Cn yn x0

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим решение y x C1 y1

... Cn yn x - линейную комбинацию этих

линейно независимых решений.

Оно имеет те же начальные условия,

что и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

выбранное

 

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

 

 

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

по теореме

Коши

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

выбранное

 

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

решение

 

 

 

и

есть

(тождественно

равно)

 

 

y x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y x C1 y1

x ... Cn yn

Поэтому

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

решение линейно

выражается через выбранные линейно независимые решения. Теорема доказана.

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

Матрица Y x , составленная из этих решений detY x W x 0 , называется

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

Теорема о структуре общего решения однородной системы.

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

 

 

 

 

 

 

 

 

 

yoo

x C1 y1 x ... Cn yn x .

 

 

 

x является

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

Проверим,

что

yoo x C1 y1 x ... Cn yn

общим решением, исходя из определения общего решения.

 

 

 

 

 

 

x -

 

 

 

 

1)

yoo x

C1 y1

x ... Cn yn

решение

однородной

системы как

линейная комбинация ее решений (теорема о свойствах решений).

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

01

 

 

 

 

2) Зададим произвольные начальные условия

y0

...

 

и покажем, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y0n

 

 

 

 

можно единственным образом выбрать набор констант

C1 ,...Cn ,

при

 

 

 

 

 

 

 

 

 

 

 

котором

yoo x0

C1 y1 x0 ... Cn yn x0

y0 .

 

Запишем

это

соотношение

покоординатно как

систему

уравнений

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

C1 ,...Cn .

C1 y11 x0 ...Cn yn1 x0 y01

79

C1 y1n x0 ...Cn ynn x0 y0n
..........................................

C1 y12 x0 ...Cn yn2 x0 y02

Определитель этой системы равен W x0 0 , так как решения линейно независимы. Поэтому набор констант C1 ,...Cn определяется из системы уравнений единственным образом. Теорема доказана.

Следствие. Общее решение однородной системы можно записать в виде

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yoo x Y x C,

 

C

, , , .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Cn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Матрица Коши (матрициант).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

условиям y x0

y0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y x

 

 

 

 

 

 

,

 

 

 

 

K x, x

 

y x

y

0

 

0

C C Y 1 x

0

y

 

y x Y x C Y x Y 1 x y

0

y

0

0

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

0

 

Матрица K x, x

 

Y x Y 1 x

 

называется матрицей Коши.

 

K x, x

 

 

 

.

 

 

0

0

y x

0

y

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема.

 

Фундаментальная

 

матрица

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

 

однородной

 

 

системе,

 

A x Y x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство. Столбцы фундаментальной матрицы являются решениями

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x A x yk x ,

k 1,...n в матрицу,

 

 

однородной системы. Объединяя запись yk

 

 

получим утверждение теоремы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формула Остроградского – Лиувилля.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

TrA x a

x a

 

x ... a

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W x C exp TrA x dx W x

0

exp

 

TrA x dx ,

22

nn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выведем формулу Остроградского – Лиувилля.

Фундаментальная матрица системы является решением однородной системы. Запишем уравнение для k –го столбца фундаментальной матрицы – координат решения yk :

y 'k

...

yk' n

 

a

 

11

 

...

 

 

 

an1

...

a

 

y

 

 

 

 

 

 

 

 

1n

 

 

k1

 

 

 

 

 

 

 

...

...

 

...

 

 

 

 

 

 

.

...

 

 

yk

 

 

ann

n

 

 

 

 

 

 

Отсюда yks' as1 yk1 as 2 yk 2 ... asn ykn .

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

80