решений однородной линейной системы n дифференциальных уравнений,
приведенного выше, следует: система n решений будет фундаментальной системой решений в интервале (a,b) тогда и только тогда, когда определитель Вронского этих решений не равен нулю хотя бы в одной точке интервала (a,b).
Матрица вида
y11 |
y12 |
|
y1n |
|
|
|
|
|
|
|
|
Ф(x) y21 |
y22 |
|
y2n |
(5.37) |
|
|
|
|
|
||
|
|||||
|
yn2 |
|
|
|
|
yn1 |
ynn |
|
столбцы которой сформированы из координатных функций линейно
независимых векторов Yk (x) |
(k 1,2,...,n) - |
решений однородной системы, |
|||
называется фундаментальной |
матрицей |
этой системы. |
Очевидно, что |
||
detФ(x) W(x). |
|
|
|
|
|
Можно показать, что фундаментальная матрица (5.37) является |
|||||
решением матричного дифференциального уравнения |
|
||||
|
dФ(x) |
A(x)Ф(x), |
(5.38) |
||
|
|
||||
|
|
dx |
|
|
|
где производная матрицы равна по определению матрице, элементы которой являются производными от соответствующих элементов исходной матрицы.
Отметим: правило дифференцирования суммы и произведения сохраняются и для матриц; однако при дифференцировании произведения
матриц |
необходимо |
сохранять |
порядок |
сомножителей |
- |
(A(x)B(x)) A (x)B(x) A(x)B (x) . |
|
|
|
||
Общее решение однородной линейной системы дифференциальных уравнений.
Как и в случае однородного линейного дифференциального уравнения n-го порядка, знание фундаментальной системы решений позволяет построить
общее решение однородной линейной системы дифференциальных уравнений.
Теорема 6 (об общем решении). |
Если система функций |
|
||
y11, |
y12, |
y1n, |
|
|
|
y22, |
y2n, |
|
|
y21, |
|
|||
|
|
|
|
(5.39) |
|
||||
|
yn2, |
ynn |
|
|
yn1, |
|
|||
является фундаментальной системой решений однородной линейной системы дифференциальных уравнений (5.24) в интервале (a,b), то функции
135
y1(x) C1y11 C2 y21 Cn yn1,
y2(x) C1y21 C2 y22 Cn yn2,
|
(5.40) |
|
|
||
yn(x) C1yn1 C2 yn2 Cn ynn, |
|
|
или в сокращенной форме |
|
|
n |
|
|
yk (x) Ci yik (x), |
(k 1,2,...,n), |
(5.41) |
i 1 |
|
|
(где C1,C2,...,Cn произвольные постоянные) дают общее решение системы (5.24) в области
a x b, |
yk |
|
(k 1,2,...,n). |
(5.42) |
Замечание. В векторной форме теорема формулируется |
так: если |
|||
система вектор - функций решений |
Yk (x) colon(yk1, yk2,..., ykn) |
(k 1,2,...,n) |
||
является фундаментальной системой решений в интервале (a,b) однородного
уравнения |
Y A(x)Y , |
то |
вектор-функция Y(x) colon(y , y |
2 |
,...,y |
n |
), где |
||||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
||
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
ym(x) Ci yim |
дает |
общее решение этого |
уравнения |
|
в |
|
области |
||||||||
i 1 |
|
|
|
|
|
|
|
|
|
|
|||||
a x b, |
|
|
ym |
|
(m 1,2,...,n). |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||||||
Доказательство. |
Совокупность |
функций (5.39) является |
решением |
||||||||||||
системы (5.24) |
при всех |
значениях |
постоянных |
C1,C2,...,Cn, |
поскольку |
||||||||||
представляет собой линейную комбинацию решений системы.
Система (5.40) представляет собой линейную алгебраическую систему уравнений относительно C1,C2,...,Cn. Ее определителем является вронскиан W(x) (5.33), который не равен нулю при любых x (a,b), так как совокупность функций (5.39) является по условию фундаментальной системой решений.
В соответствии с определением общего решения нормальной системой дифференциальных уравнений, учитывая вышеприведенные рассуждения, можно сделать вывод: совокупность функций (5.40) является общим решением системы (5.24) в области изменения переменных (5.42).
Покажем, что формула (5.40) содержит все решения системы (5.24). Для этого следует доказать, что при любых начальных условиях, наложенных на
функции y1,..., yn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
y |
y(0) |
, |
y |
2 |
y(0),..., y |
n |
y |
(0) при |
x x |
0 |
, |
|
|
|
(5.43) |
|
|
|
|
1 |
1 |
|
|
2 |
|
n |
|
|
|
|
|
|
|||
где (x |
0 |
, y(0) |
, y |
(0),...,y(0)) |
|
- произвольная фиксированная точка области (5.42), |
|||||||||||||
|
1 |
|
2 |
n |
|
|
|
|
|
|
|
|
|
|
C(0) |
|
(0) |
,...,C(0) |
|
можно единственным образом найти значения постоянных |
,C |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
2 |
n |
таких, что полученное частное решение системы (5.24) с начальными
условиями |
(5.43) |
будет |
включаться |
в |
решение |
(5.40) |
при |
|
|
|
136 |
|
|
|
|
C C(0) |
,C |
2 |
C(0),...,C |
n |
C(0) . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
1 |
1 |
|
2 |
|
|
|
n |
|
|
|
|
x x0, учитывая (5.43) получим |
|
||||||||||||||
|
Действительно, полагая в (5.40) |
|
|||||||||||||||||||||||||
|
|
|
|
y(0) |
C (y |
|
) |
0 |
C |
(y |
21 |
) |
0 |
C |
(y |
n1 |
) |
, |
|
|
|||||||
|
|
|
|
1 |
1 |
11 |
|
2 |
|
|
|
n |
|
0 |
|
|
|
||||||||||
|
|
|
|
y(0) |
C (y |
21 |
) |
0 |
C |
(y |
22 |
) |
0 |
C |
|
(y |
|
) |
|
, |
|
||||||
|
|
|
|
2 |
1 |
|
|
2 |
|
|
|
|
n |
|
|
2n 0 |
|
(5.44) |
|||||||||
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
y(0) |
C (y |
n1 |
) |
0 |
C |
(y |
n2 |
) |
0 |
C |
|
(y |
|
) |
|
, |
|
||||||
|
|
|
|
n |
1 |
|
|
2 |
|
|
|
|
n |
|
|
nn 0 |
|
|
|||||||||
где (yij )0 yij (x0). Полученная система линейных неоднородных уравнений относительно неизвестных C1,C2,...,Cn имеет единственное решение, поскольку ее определитель равен W(x0) - вронскиану фундаментальной системы функций
(5.39) при x x0, которой отличен от нуля. Поэтому
C1 C1(0)(x0, y1(0),...,yn(0)),C2 C2(0)(x0, y1(0),...,yn(0)),...,
Cn Cn(0)(x0, y1(0),...,yn(0)).
(5.45)
Подставляя эти значения постоянных в общее решение (5.44) найдем решение задачи Коши (5.24), (5.43).
Таким образом, чтобы линейная комбинация n решений однородной линейной системы (5.24) с произвольными постоянными C1,C2,...,Cn давала общее решение этой системы, необходимо и достаточно, чтобы эти решения были линейно независимы, т.е. чтобы они составляли фундаментальную систему решений.
З а м е ч а н и е. Общее решение линейной однородной системы дифференциальных уравнений можно выразить через фундаментальную матрицу (5.37). Формулу (5.44) можно записать в виде
|
Y(x) Ф(x)C, |
|
(5.46) |
|
где Y(x) colon(y1(x), y2(x),..., yn(x)) , yk (x) - функции из (5.44), |
Ф(x)- |
|||
фундаментальная матрица (5.37), |
C colon(C1,C2,...,Cn) - постоянный вектор с |
|||
произвольными координатами. |
|
|
|
|
Если положить в |
(5.46) |
x x0, получим |
Y(x0) Ф(x0)C . |
Отсюда, |
поскольку матрица Ф(x0) |
не вырождена, имеем C Ф 1(x0)Y(x0). Подставляя C |
|||
в формулу (5.46), получаем |
|
|
|
|
|
Y(x) Ф(x)Ф(x0)Y(x0). |
|
(5.47) |
|
Полученная формула дает решение начальной задачи (5.24), (5.43), т.е. решение задачи Коши для линейной однородной системы дифференциальных уравнений.
Матрица
K(x,x0) Ф(x)Ф 1(x0) |
(5.48) |
называется матрицей Коши. Используются также другие названия этой
137
матрицы – импульсная матрица, а также матрицант.
Таким образом решение задачи Коши (5.24), (5.43) можно записать с помощью этой матрицы в виде
Y(x) K(x,x0)Y(x0). |
(5.49) |
Пример. Найти общее решение линейной однородной системы
dy |
(x x2)y (x3 x2 x 1)z, |
||||
|
|
|
|||
|
|
||||
|
|
|
|
||
dx |
|
2 |
|
||
dz |
(1 x)y (x |
x 1)z. |
|||
|
|
|
|
||
dx |
|
|
|
||
Решение. Из вида системы ясно, что одно из его решений представляет собой многочлен. Поскольку в правой части уравнений системы коэффициент при функции имеет степень на единицу большую, чем коэффициент при 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)ex, |
z |
2 |
(A B |
2 |
x)ex . |
|
|||||||||
|
|
1 |
|
1 |
|
|
|
|
|
2 |
|
|
|
|
||||
Подставим эти функции в |
уравнения системы, разделим обе части |
|||||||||||||||||
уравнений на |
ex , затем приравняем коэффициенты в правой и левой частях |
|||||||||||||||||
полученных уравнений при одинаковых степенях 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 ex. |
|
|
|
|
|
|||||||||
Проверим линейную независимость решений. Вычислим определитель |
||||||||||||||||||
Вронского |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
W |
(x) |
|
x2 |
1 |
|
xex |
|
e |
x |
|
0. |
|
|||
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
x |
|
ex |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Таким образом, общее решение системы имеет вид
138
y C1y1 C2 y2 C1(x2 1) C2xex,
zC1z1 C2z2 C1x C2ex.
3.Свойства решений неоднородной линейной системы
дифференциальных уравнений.
Рассмотрим далее неоднородную линейную систему (5.23)
dyk |
n |
|
|
aki(x)yi fk (x), |
(k 1,2,...,n) , |
||
dx |
|||
i 1 |
|
которую в векторной и операторной формах мы записывали следующим образом
|
|
|
dY |
A(x)Y F, |
|
L[Y]F. |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
dx |
|
|
|
|
|
|
|
Установим наиболее важные свойства решений линейной неоднородной |
||||||||||
системы. |
|
|
|
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
где i - |
|
Теорема 7 (принцип суперпозиции). |
Пусть F(x) iFi (x), |
|||||||||
|
|
|
|
|
(i)(x),..., f |
|
|
i 1 |
|
|
постоянные |
числа, |
F (x) colon(f |
(i) |
(x)), и пусть |
Y (x) |
- |
решение |
|||
|
|
i |
1 |
n |
|
i |
|
|
||
уравнения Y A(x)Y F (x). Тогда вектор-функция |
|
|
|
|||||||
i |
i |
i |
|
|
|
|
|
|
||
m
Y(x) iYi(x)
i 1
является решением уравнения (5.23). |
|
|
|
|
|
|
|
|
|||
Доказательство |
вытекает |
из |
линейности |
оператора L (5.23). |
|||||||
Действительно, |
|
m |
m |
|
m |
|
dY |
|
m |
||
|
d |
|
|
|
|||||||
L[Y] |
|
iYi |
A(x) iYi |
i |
i |
A(x)Yi |
iFi. |
||||
|
|
||||||||||
|
dx i 1 |
i 1 |
|
i 1 |
|
dx |
|
i 1 |
|||
Таким образом, Y(x) является решением уравнения (5.23).
Физический смысл принципа суперпозиции заключается в следующем. Пусть имеется физическая или техническая система, внутренние свойства которой описываются линейным оператором L. Пусть на вход такой системы
подается сигнал |
Fk (x) (k 1,2,...,m), |
а на |
выходе фиксируется |
сигнал |
Yk (x) (k 1,2,...,m). |
Если на вход |
такой |
системы подать |
сигнал |
m
F(x) iFi (x) , то на выходе из устройства будет наблюдаться сигнал
i 1 m
Y(x) iYi (x) (рис. 12).
i 1
139