Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
BORZENKOV_MATLAB ОДУ.pdf
Скачиваний:
169
Добавлен:
11.05.2015
Размер:
1.91 Mб
Скачать
A(t).

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

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

Системой линейных однородных дифференциальных уравнений называется система следующего вида:

n

 

 

 

 

 

 

 

y'i aik (t) yk ,

i

 

 

 

(1.1)

1,n,

 

k 1

 

 

 

 

 

 

 

где функции aik (t) непрерывны на некотором интервале T .

 

 

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

 

Введем

 

T

 

 

a11 (t) a1n (t)

 

следующие обозначения: y(t) (y1(t),..., yn(t))

A(t)

 

 

. Тогда

,

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

 

 

 

 

 

 

 

 

систему (1.1) можно записать в виде матричной формы

an1 (t) ann (t)

 

 

 

 

 

y' A(t)y.

 

 

 

 

 

 

(1.2)

Решением системы (1.2) называется непрерывно дифференцируемая векторфункция y(t), которая при подстановке в систему обращает все уравнения в

тождества. Если векторы y1(t),...,ym(t)

являются решениями системы (1.2), то

m

(t) , где C1,..., Cm – действительные

любая линейная комбинация Ci yi

i 1

 

произвольные постоянные, вновь является решением системы (1.2).

Если требуется найти решение системы (1.2), удовлетворяющее условию y(t0 ) y0 , то говорят, что для системы (1.2) поставлена задача Коши и

записывают ее в виде

 

 

y' A(t)y ,

y(t0 ) y0 .

(1.3)

Замечание. Из общей теоремы Коши о существовании и единственности решения нормальной системы дифференциальных уравнений сразу следует, что, если элементы матрицы A(t) – функции aik (t) непрерывны на отрезке

T [t0,t*], то на нем существует единственное решение задачи Коши (1.3). Теорема Коши не только дает условия существования и единственности

решения задачи

Коши.

При выполнении этих

условий решение

y(t) (y (t), , y

n

(t))

определено

на отрезке t

,t* , т.е. в концевых точках t t

0

и

1

 

 

0

 

 

 

t t* функции y (x), i 1,2, ,n, имеют односторонние производные.

 

 

 

 

i

 

 

 

 

 

 

Теорема

 

Коши дает

достаточные условия

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

 

и

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

73

Напомним,

что векторы y1(t),..., yn(t)

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

интервале

T , если тождество

C1y1(t) Cnyn(t) 0 выполняется тогда и

только тогда, когда все коэффициенты Ci равны нулю одновременно.

 

Пусть

векторы 1(t),..., n(t)

непрерывны на интервале

T .

 

 

1, , n

 

11(t) 1n(t)

 

 

 

 

 

 

 

 

Функциональный определитель

 

 

 

, где ik (t) –

 

 

 

 

 

n1(t) nn(t)

 

 

 

координаты вектора k (t) , называется определителем Вронского, либо

вронскианом векторов 1(t),..., n(t) .

Теорема 1. Если определитель Вронского решений y1(x), , yn(x) системы

(1.2) тождественно равен нулю на интервале T , тогда эти решения линейно зависимы на T . Если определитель Вронского решений y1(x), , yn(x)

уравнения (1.2) не равен нулю ни в одной точке интервала T , тогда эти решения линейно независимы на T .

Замечание. Если 1(t),..., n(t) – произвольные векторы, то из равенства нулю их определителя Вронского, вообще говоря, не следует их линейная

зависимость. Действительно, рассмотрим два вектора

1(t) ( 0 1)T ,

2(t) (t 0 )T . Векторы линейно-независимы, например, на

отрезке

1,1 ,

так как условие

C (t) C

2

(t) 0 выполняется только при

C 0 и

C

2

0.

 

1

1

2

 

1

 

 

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

Часто бывает полезен другой критерий линейной независимости

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема 2. Для того чтобы произвольные

векторы i

(t), i 1,

2,…,n были

линейно независимы на a,b , необходимо и достаточно,

чтобы определитель

 

 

 

 

1, , n

 

( 1, 1)( 1, 2) ( 1, n)

 

 

 

 

 

 

 

 

 

 

Грама этих векторов

 

 

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

 

 

был отличен

 

 

 

 

 

 

 

( n, 1)( n, 2) ( n, n)

 

 

 

 

, ,

 

0,

 

 

 

 

 

 

 

 

b

 

(t),

 

(t))dt

 

от нуля:

n

где

(

i

,

j

) (

i

j

– скалярное

1

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

произведение функций i(t) , j(t) i,

j 1,n

на отрезке.

 

 

 

 

 

74

1.1. Фундаментальная матрица

 

 

 

Определение. Фундаментальной

системой решений y1(t),..., yn(t) системы

(1.2) называются любые n линейно независимых решений системы (1.2).

Если коэффициенты aij(t),

i, j 1,2, ,n, в

системе

уравнений (1.2)

непрерывны на t0,t* ,

то она

имеет систему

решений

y1(t), , yn(t) ,

фундаментальную на t0,t*

.

 

 

 

Определение. Матрица W(t), столбцами которой являются координаты векторов, образующих фундаментальную систему решений, называется фундаментальной матрицей системы (1.2).

Определитель матрицы W(t) – это определитель Вронского системы n линейно независимых решений системы (1.2). Определитель не равен нулю,

поэтому матрица W(t) имеет обратную W 1 t .

 

 

Пусть каждый из линейно независимых векторов

y1(t),...,yn(t) является

решением

матричного уравнения (1.2). Эту совокупность n

векторных

уравнений можно кратко записать в следующем виде:

 

 

 

W (t)

A(t) W (t) .

 

(1.4)

Столбцами матрицы W являются координаты векторов

y1(t),..., yn(t). Каждый

столбец

матрицы W (t) равняется

произведению

матрицы

A t на

одноименный столбец матрицы W(t) . Уравнение (1.4) называется матричным уравнением, сопоставленным векторному уравнению (1.2).

Очевидно, что фундаментальная матрица W(t) есть решение матричного уравнения (1.4).

Теорема 3. Линейная однородная система всегда имеет фундаментальную систему решений, а значит, и фундаментальную матрицу. По заданной системе n линейно независимых векторов y1(t),...,yn (t) можно найти единственную систему (1.2), для которой эти векторы образуют фундаментальную систему решений.

Пусть матрица W t , столбцами

которой

являются

координаты

этих

векторов, удовлетворяет системе (1.4). Требуется найти матрицу

A t .

Умножая тождество (1.4) справа на

обратную

матрицу

W 1 t ,

получаем

W'(t)W 1(t) A(t)W(t)W 1(t) A(t), где

A t

матрица

искомой

системы

уравнений.

 

 

 

 

 

Замечание. Множество решений линейной однородной дифференциальной системы образует линейное пространство функций. Любая фундаментальная система решений является базисом этого пространства. Существует бесконечно

75

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

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

y W(t)C ,

(1.5)

где W (t) – фундаментальная матрица, а C

– произвольный постоянный

вектор. Правую часть (1.5) удобно представлять как скалярное произведение строки (y1(t),...,yn(t)) векторов, столбцы координат которых образуют матрицу

W t , на столбец C : y(t) C1y1(t) Cn yn (t). Получается линейная комбинация векторов фундаментальной системы, коэффициентами которой служат координаты вектора C . По аналогии строилось общее решение линейного однородного уравнения n-го порядка.

1.2. Матрицант. Матрица Коши. Пусть фундаментальная матрица W(t)

удовлетворяет условию W (t0 ) E , т.е. векторы фундаментальной системы имеют единичную матрицу начальных значений. Тогда фундаментальная матрица W (t) называется матрицантом системы (1.2). Матрицант – аналог нормальной фундаментальной системы решений линейного однородного уравнения n-го порядка. Матрицант, очевидно, является решением следующей однородной дифференциальной системы:

 

d

 

 

 

 

W (t) A(t) W (t), W (t0 )

E .

(1.6)

 

 

 

dt

 

 

Пусть W (t) – произвольная фундаментальная матрица. Тогда матрица

 

 

 

K(t,t0 ) W(t)W 1(t0 ),

 

(1.7)

зависящая от двух аргументов t и t0 , также является матрицантом системы

(1.2). Действительно, матрица K(t,t0 ) удовлетворяет следующей системе:

d

 

 

K(t,t0 ) A(t) K(t,t0 ) , K(t0 ,t0 ) E.

(1.8)

 

dt

 

Поскольку K'(t,t0) (W(t) W 1(t0))' W'(t)W 1(t0) A(t)W(t)W 1(t0) A(t)K(t,t0) и

при t t0 , выполнено начальное условие K (t0 ,t0 ) W(t0 )W 1(t0 ) E.

Принято также другое название матрицы, определяемой формулой (1.7), – матрица Коши системы (1.2).

Замечание. Матрица Коши K(t,t0 ) определяется по формуле (1.7)

единственным образом, несмотря на то, что матрица W (t) – произвольная фундаментальная матрица. Это следует из того, что матричное уравнение (1.2) с единичной матрицей начальных значений имеет единственное решение.

76

Если известна матрица Коши, то решение однородной задачи Коши

y' A(t)y , y(t0 ) y0 имеет следующий вид: y K(t,t0 )y0 . Справедливо и обратное. Если для произвольного вектора начальных значений y0 решение однородной задачи Коши y' A(t)y , y(t0 ) y0 может быть записано в форме

y K(t,t0 )y0 , то входящая в нее матрица является матрицей Коши. Замечание. Полученное представление является аналогом формулы, выражающей решение задачи Коши для линейного однородного уравнения n-го порядка через нормальную фундаментальную систему решений.

Поясним физический смысл матрицы Коши. Для этого рассмотрим следующую однородную дифференциальную систему:

 

 

 

 

d

y A(t) y,t [t0,t

*

],

 

 

 

 

 

 

 

 

 

(1.9)

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

y(t

0

) y .

 

 

 

 

 

 

 

 

 

 

0

 

 

 

Пусть

вектор

Ki (t, )

является

решением

 

следующего

векторного

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

уравнения:

 

 

d

Ki(t, ) A(t) Ki(t, ) с начальным условием

 

 

 

 

 

 

ei (0,...,

 

dt

 

 

 

 

K i(t0,t0) ei. Здесь

0,1.., 0,..., 0 )T – нулевой вектор с единицей

на i

 

 

 

 

 

 

 

 

 

i

K(t, )

 

 

 

месте.

Обозначим

через

n n-матричную

функцию,

составленную из векторов Ki(t, ) как из столбцов, где i принимает значения от 1 до n.

Функция K(t, ) однозначно определена для любого значения аргумента t. Эта функция абсолютно непрерывна и удовлетворяет следующей системе дифференциальных уравнений:

d

 

 

 

K(t, ) A(t) K(t, ),

(1.10)

 

dt

 

 

 

 

K(t0,t0 ) E.

 

Отсюда следует физический смысл матрицы Коши K(t, )

(функции Коши).

Каждый i-й столбец функции

Ki(t, ) представляет собой отклик в момент

времени t однородной системы (1.9), которую возмутили в момент времени t0

единичным вектором y0 ei . Поэтому функции Коши часто называют еще функциями точечного источника (функциями Грина).

1.3. Линейные однородные системы с постоянными коэффициентами.

Построение ФСР линейной однородной системы с постоянными коэффициентами сводится к решению характеристического уравнения для матрицы системы и построению соответствующих собственных векторов.

77

 

 

d y1

 

7y 3y

2

,

 

 

 

 

dt

1

 

 

 

 

 

 

 

Пример 1. Найти общее решение и ФСР системы уравнений

 

 

 

 

 

 

d y2

 

6y1 4y2 ,

 

 

 

 

 

dt

 

 

 

 

 

d y1(t)

7 3

y1(t)

т. е.

 

 

y

 

 

6 4

 

 

. Выпишем общее решение системы и

 

 

dt

2 (t)

 

 

y2

(t)

найдем векторы ФСР. Общее решение считается через собственные числа и соответствующие им векторы матрицы A. Получаем: A x x ; Ax Ex 0;

(A E)x 0;

det(A E) 0. Отсюда

7 3

0,

или 2

11 10 0,

6

4

 

 

 

 

 

 

1 1, 2 10.

Собственные

числа являются

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

и простыми.

Найдем собственный вектор

x1 , соответствующий собственному значению 1.

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

 

(A 1E)

x1 0.

Получаем x

 

1

 

 

1. Для второго вектора аналогично x

2

 

1

2 10.

1

 

 

 

1

 

 

 

 

 

 

 

 

2

 

 

 

 

1

 

 

Выписываем общее решение однородной дифференциальной системы в виде

y(t) C1

 

1

 

e

t

C2

 

1

10t

. В матричной форме относительно

C

и C

 

 

 

 

 

e

2

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

1

 

 

 

 

 

 

et

 

получаем y(t)

 

t

2e

 

e10t

C1

 

10t

 

. Векторы ФСР

e

C2

 

et

 

e10t

 

 

 

t

и 10t

.

2e

 

e

 

dy

1

 

4y1 3y2 ,

 

 

 

 

 

dt

 

 

Пример 2. Найти общее решение и ФСР системы уравнений

 

 

 

dy2

 

3y1 4y2 ,

 

 

 

 

dt

 

 

d y1(t)

т.е.

dt y2 (t)

 

4

3

y (t)

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

 

 

4

 

1

 

 

3

y2 (t)

 

относительно матрицы A: det(A E) 0. Отсюда следует

4

3

0;

3

4

 

 

 

 

(4 )2

9; 4 3i;

4 3i . Найдем собственный вектор x ,

 

1,2

 

 

1

 

78

соответствующий собственному значению 1. Для этого решаем систему

 

 

 

линейных алгебраических уравнений (A 1E)

 

x1 0. Получаем

 

 

 

 

1

 

 

 

4 3i. Для второго вектора аналогично x2

1

2 4 3i.

 

 

 

x1 1

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

Общее

решение

нашей

системы

имеет

вид

 

 

 

 

 

 

 

 

 

 

 

 

 

1

(4 3i)t

1

(4 3i)t

. Переписывая в матричной форме

 

 

 

y(t) c1

e

 

 

c2 i e

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

e(4 3i)t

e(4 3i)t

 

C1

 

 

 

 

и c2 , получаем y(t)

 

(4 3i)t

 

 

(4 3i)t

 

. Векторы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ie

 

 

 

ie

 

 

 

C2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e(4 3i)t

 

e(4 3i)t

 

 

 

 

 

 

 

 

 

 

 

 

 

ФСР соответственно будут

(4 3i)t

 

и

 

(4 3i)t . Выделим в общем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ie

 

 

 

 

ie

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решении действительную и мнимую части. Имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e(4 3i)t

 

e4t

cos3t ie4t sin3t, ie(4 3i)t

e4t sin3t ie4t cos3t.

 

 

 

Отсюда

 

получаем

y

(t)

e4t cos3t

e4t sin 3t

c1

 

Выписываем

векторы

 

 

4t

 

e

4t

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

e

 

sin 3t

 

cos3t

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

4t

 

 

 

 

 

e4t sin3t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ФСР e

 

cos3t

 

и

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e4t sin3t

 

 

e4t cos3t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy1

 

 

5y

y

2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

Найти общее решение и ФСР системы уравнений dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy2

y1 3y2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d y1(t)

5

1

y1(t)

 

 

 

 

 

 

 

 

 

Перепишем

систему

в виде

 

 

 

 

 

 

 

 

 

 

 

.

Решаем уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

y2(t)

1

3 y2(t)

 

 

 

 

 

 

 

 

 

det(A E) 0:

 

5

1

 

 

0;

2 8 16 0;

 

1 2

4.

Если 1

– корень

 

 

 

 

 

 

 

 

 

1

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

характеристического уравнения кратностью m, то этому корню соответствует

решение

x

p

(t)e 1t , x

2

p

2

(t)e 1t ,

…,

x

n

p

n

(t)e 1t ,

где p (t),

p

2

(t), …,

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

pn (t) – многочисленные

степени

не

 

выше

 

m 1. Таким

образом,

действительному

 

корню

 

4 кратностью

два

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

решения

y

e4t (a t a ),

y

2

e4t (b t b ). Дифференцируя y

и y

2

,

получаем тождества

1

 

1

2

 

 

1

 

2

 

 

 

 

 

 

1

 

 

 

 

 

 

79

dy1

ae4t 4(at a

)e4t ,

dy2

be4t 4(bt b )e4t . Значения

x ,

x

,

dx1

,

dx2

dt

dt

1

1

2

 

1

1

2

1

2

 

dt

dt

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

a1 4(a1t a2 ) 5(a1t a2 ) (b1t b2 ); b1 4(b1t b2 ) a1t a2 3(b1t b2 ).

Приравнивая коэффициенты при t

и свободные члены, получаем следующие

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

4a 5a b,

a 4a 5a b ,

 

1

 

1

1

 

1

2

2

2 . Тогда

 

 

 

 

 

4b1 a1 3b1,

b1 4b2 a2 3b2.

следует, что

a1 b1,

a2 b2

a1

b1.

Полагая

 

a1 C1,

a2

C2

(C1,C2

произвольные постоянные), находим b1 C1,

b2 C2

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

 

 

x e4t (Сt С

),

x

2

e4t (Сt С

2

С ).

 

 

 

 

1

1

2

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

4 t

 

4 t

 

C

 

 

 

 

В матричном

виде имеем y(t)

e t

 

e

 

1

.

Соответствующими

4 t

(t 1)

4 t

 

 

 

 

 

e

e

C2

 

 

 

 

e4tt

 

e

векторами ФСР будут

4t

 

 

и

e

 

(t 1)

e

4t

4t

.

Замечание. Систему можно решить методом исключения. Действительно,

выразив из первого уравнения

x2 и продифференцировав, подставим затем

значения

x

 

и

 

x

во второе

 

уравнение.

В

 

результате получим

линейное

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

однородное уравнение второго порядка относительно x1.

 

Выражаем

x

 

из первого уравнения:

x

2

5x

x .

Дифференцируя,

получаем

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

x

5x x . Подставляем полученные значения во второе уравнение системы

 

2

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и находим общее решение для x1:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5x

x

x 15x 3x ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

1

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1 8x1 16x1 0.

 

 

 

 

 

Решая характеристическое уравнение, получаем 1 2

 

4. Отсюда общее

решение для

 

y:

y e4t(C C t)

. Теперь подставляем

x

 

в выражение для x :

 

 

 

 

 

1

 

1

 

1

2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

x

2

5x x

,

x

2

5e4t (C C

t) 4e4t (C C

t) e

4tC

2

. В итоге имеем общее

 

1

1

 

 

 

 

1

2

 

 

 

 

 

 

1

2

 

 

 

 

решение для

x :

x

2

e4t (C

С

2

C

2

t). Окончательно

 

 

 

 

 

 

 

 

 

2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e4tt

e4t

C1

 

e4tt

e

y(t)

e

4t

(t 1)

e

4t

 

. Векторы ФСР

 

4t

(t

 

и

 

 

 

C

 

e

 

1)

e

 

 

 

 

2

 

 

 

 

 

4t

4t

.

80

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