где y1, ..., yn – искомые функции независимой переменной x. Будем предполагать, что функции f1, ..., fn, стоящие в правых частях уравнений системы (2.1), заданы в некоторой области D (n+1)-мерного пространства переменных x, y1, ..., yn. Упрощая изложение и несколько ограничивая общность, будем считать, что область D задается неравенствами
a0 < x < b0, a1 < y1 < b1,
. . . . . . . . . . . . .
an < yn < bn
(D = (a0, b0) × (a1, b1) × ... × (an, bn), где −∞ ≤ ai < bi ≤ +∞ для всех i = 0, 1, ..., n).
Решением системы (2.1) на промежутке (a, b) (a0, b0) называется набор n функций y1(x), ..., yn(x), определенных и дифференцируемых на (a, b), обладающих свойством
(x, y1(x), ..., yn(x)) D при x (a, b),
и обращающих все уравнения (2.1) в тождества для x (a, b). Геометрически решение системы (2.1) представляет собой кривую в
(n + 1)-мерном пространстве переменных (x, y1, ..., yn). Эта кривая называется интегральной кривой системы (2.1).
Рассмотрим, например, систему
dy |
|
|
|
|
|
||
|
1 |
= 3√3 |
|
y2, |
|||
y1 |
|||||||
dx |
|||||||
dy2 |
= y1 . |
||||||
|
|
|
|
|
|||
dx |
x3 |
||||||
В правых частях уравнений системы стоят функции
√ y1 f1(x, y1, y2) = 3 3 y1 y2, f2(x, y1, y2) = x3 ,
определенные в области D, задаваемой неравенствами
0 < x < +∞, −∞ < y1 < +∞, −∞ < y2 < +∞
˜
(можно рассматривать систему и в области D, заданной неравенствами −∞ < x < 0, −∞ < y1 < +∞, −∞ < y2 < +∞). Нетрудно убедиться, что решением этой системы является, например, пара функций y1(x) = x3, y2(x) = x при 0 < x < +∞.
Для удобства будем записывать систему (2.1) в векторном виде. Рассмотрим вектор
y1
Y = ... Rn. yn
25
Напомним, что его производной служит вектор
dY
dx
dy1
dx
=. . . .
dyn
|
|
dx |
|
Рассмотрим, кроме того, вектор-функцию |
|
||
F(x, Y) = f.1.(.x,. . .Y. .) |
, |
||
|
|
fn(x, Y) |
|
определенную в D. Ясно, что система равенств (2.1) равносильна одному |
|||
векторному равенству |
|
||
|
dY |
= F(x, Y), |
(2.2) |
|
|
||
|
dx |
|
|
которое и будем называть векторной формой записи нормальной системы дифференциальных уравнений. Решением системы (2.2) на (a, b) является вектор-функция
y1(x)
Y(x) = . . . . . , yn(x)
компоненты которой – функции набора y1(x), ..., yn(x), служащие решением системы (2.1).
Как и в случае дифференциального уравнения первого порядка, система (2.2) имеет обычно бесконечное множество решений.
Для выделения из этого множества одного решения ставится задача Коши: фиксируется точка (x0, Y0) D и ищется решение Y(x) на проме-
жутке (a, b), содержащем точку x0, и такое, что |
|
Y(x0) = Y0. |
(2.3) |
Геометрически это означает, что ищем интегральную кривую, проходящую через точку (x0, Y0).
Если
y10
Y0 = ... , yn0
то нахождению решения Y(x) задачи Коши (2.2)–(2.3) соответствует нахождение решения y1(x), ..., yn(x) системы (2.1), обладающего свойством
y1(x0) = y10, ..., yn(x0) = yn0.
26
Определим матрицу Якоби (матрицу частных производных) векторфункции F по переменной Y:
|
|
|
|
∂f1 |
∂f1 |
|
||
∂F |
|
|
... |
|
||||
|
∂y1 |
∂yn |
||||||
|
= |
|
∂f. . .n. . . .∂f. .n. |
. |
||||
∂Y |
||||||||
|
|
|
|
|
... |
|
|
|
|
|
|
|
∂yn |
||||
|
|
|
|
∂y1 |
|
|||
|
|
|
|
|
|
|
|
|
Будем говорить, что вектор-функция или матрица непрерывна, если непрерывны все ее компоненты.
Сформулируем теорему, дающую достаточные условия существования и единственности решения задачи (2.2)–(2.3).
Теорема 2.1. Если вектор-функция F и матрица ∂F непрерывны в
∂Y
области D, то для любой точки (x0, Y0) D существует единственное решение Y = Y(x) задачи Коши (2.2)–(2.3), заданное на некотором интервале (a, b), содержащем точку x0.
Эта теорема означает, что при указанных условиях через каждую точку области D проходит интегральная кривая и притом единственная. Отметим, что в общем случае интервалы определения решений различных задач Коши для системы (2.2) различны.
Рассмотрим теперь дифференциальное уравнение порядка n, разре-
шенное относительно старшей производной |
|
y(n) = f(x, y, y0, ..., y(n−1)) |
(2.4) |
(о таких уравнениях уже упоминалось во введении (см. (B.2)).
Будем предполагать, что функция f определена и непрерывна в области D (n + 1)-мерного пространства своих агрументов.
В качестве примера рассмотрим уравнение третьего порядка
y000 = x2 + y00(1 + ln(y0)).
Функция, стоящая в правой части, определена и непрерывна в области D
пространства переменных x, y, y0, y00, задаваемой неравенствами |
|
|
|
||||||||||||||
−∞ |
< x < |
+∞ |
, |
−∞ |
< y < |
+∞ |
, |
0 |
< y0 |
< |
+∞ |
, |
−∞ |
< y00 |
< |
+∞ |
. |
|
|
|
|
|
|
|
|
|
|
||||||||
Решением уравнения (2.4) на промежутке (a, b) будем называть функцию y(x), имеющую все производные до порядка n на (a, b), такую, что
(x, y(x), y0(x), ..., y(n−1)(x)) D при x (a, b),
и обращающую уравнение (2.4) в верное тождество на (a, b). 27
Вместе с уравнением (2.4) будем рассматривать нормальную систему
с неизвестными функциями z1, ..., zn следующего вида: |
|
|||||
|
dz |
|
|
|
|
|
|
1 |
|
= z2, |
|
||
dx |
|
|||||
dz2 |
= z3, |
|
||||
|
|
|
|
|
|
|
dx |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.5) |
. . . . . . . . . . . . . . . . . . . . . |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dzn |
|
1 |
|
|
||
|
|
− |
|
= zn, |
|
|
|
dx |
|
|
|
||
dzn
= f(x, z1, ..., zn),
dx
где функция f та же, что в правой части уравнения (2.4).
Рассмотрим некоторое решение y(x) уравнения (2.4) на (a, b) и постро-
им n функций: |
|
z1(x) = y(x), z2(x) = y0(x), ..., zn(x) = y(n−1)(x). |
(2.6) |
Очевидно, что набор z1(x), ..., zn(x) является решением системы (2.5) на
(a, b).
Пусть теперь z1(x), ..., zn(x) – решение системы (2.5) на (a, b). Тогда функция y(x) = z1(x) – решение уравнения (2.4) на (a, b), при этом выполнены соотношения (2.6).
Говорят, что уравнение (2.4) сводится к нормальной системе (2.5). При этом неправильно утверждать, что уравнение (2.4) и система (2.5) эквивалентны, так как обычно в математике эквивалентность уравнений или систем означает, что множества их решений совпадают. В данном случае решения уравнения (2.4) и системы (2.5) имеют различную природу – решением (2.4) является одна функция, а решением (2.5) - набор n функций.
Приведем пример сведения уравнения вида (2.4) к системе. Уравнение второго порядка
сводится к системе |
|
y00 = x sin(y) − (y0)3 |
|
|
||||||
dz |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
||
|
1 |
= z2, |
|
|
|
|
|
|
||
|
dx |
|
|
|
|
|
|
|||
|
|
dz2 |
= x sin(z |
) |
− |
z3. |
|
|
||
|
|
|
|
|||||||
|
dx |
|
1 |
|
2 |
, y |
, y0 |
, ..., y(n−1)) области D. |
||
Рассмотрим для |
уравнения |
(2.4) точку (x |
||||||||
|
|
|
|
|
0 |
0 |
0 |
0 |
||
Пусть
z1
z2
Z = ... ,
zn
|
|
|
y0 |
|
|
|
Z0 |
= |
y.00 |
R |
n. |
||
|
|
.. |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y0(n−1)
28
Из ранее изложенного следует, что задаче Коши с начальным условием Z(x0) = Z0 для системы (2.5) соответствует задача нахождения решения y(x) уравнения (2.4), для которого
y(x0) = y0 |
, y0(x0) = y0 |
, |
..., y(n−1)(x0) = y(n−1) |
(2.7) |
|
0 |
|
0 |
|
(эта задача также называется задачей Коши). Таким образом, в задаче Коши для дифференциального уравнения порядка n задаются значения самого решения y и его производных y0, ..., y(n−1) в некоторой точке x0.
Сформулированная ранее теорема существования и единственности для нормальной системы позволяет доказать аналогичный результат для уравнения порядка n.
Теорема 2.2. Если функция f в правой части уравнения (2.4) непрерывна в области D вместе со своими частными производными по переменным y, y0, ..., y(n−1), то для любой точки (x0, y0, y00 , ..., y0(n−1)) D существует единственное решение задачи Коши (2.4)–(2.7), заданное на некотором интервале, содержащем точку x0.
В некоторых случаях удается понизить порядок уравнения (2.4). Отметим два из них.
Случай 1. Уравнение не содержит y, y0, ..., y(k), т. е. имеет вид
y(n) = f(x, y(k+1), ..., y(n−1)). |
(2.8) |
|
В этом случае перейдем от переменных x, |
y к переменным x, |
z, где |
z = y(k+1). Так как |
|
|
z0 = y(k+2), ..., z(n−k−2) = y(n−1), |
z(n−k−1) = y(n), |
|
уравнение (2.8) сводится к уравнению |
|
|
z(n−k−1) = f(x, z, ..., z(n−k−2) ), |
(2.9) |
|
порядок которого равен n − k − 1. Если решение z(x) уравнения (2.9) найдено, y(x) получится из него последовательным интегрированием.
Случай 2. Уравнение не содержит x, т. е. имеет вид
y(n) = f(y, y0, ..., y(n−1) ). |
(2.10) |
Перейдем от переменных x, y к переменным y, z = y0 (y становится независимой переменной, т. е. ищется y0 = z(y) ). Пересчитаем производные:
y0 = z, y00 = |
dy0 |
= |
d(z(y)) |
= |
d(z(y)) |
|
dy |
|
= |
dz |
z |
(2.11) |
dx |
dx |
dy dx |
|
|||||||||
|
|
|
|
dy |
|
|||||||
29