где K - положительное число, M(x, y1, , yn) D. Тогда система имеет единственное решение
y1 1(x), , yn n(x),
удовлетворяющее начальным условиям y1 y1(0), , yn yn(0) при x x0. Это решение определено и непрерывно дифференцируемо (все k (x) имеют непрерывные производные) в интервале x x0 h, где h min a,b
M .
3. Понятие общего, частного и особого решений.
Вупрощенной формулировке общее решение системы (5.1)
определяется как семейство решений этой системы, зависящее от n
произвольных параметров C1,C2, ,Cn |
|
|
y1 1(x,C1, ,Cn), |
|
|
y2 2(x,C1, ,Cn), |
(5.8) |
|
|
||
|
yn n(x,C1, ,Cn).
Геометрически общее решение есть семейство интегральных кривых в (n 1)-мерном пространстве x, y1, , yn, зависящее от n параметров
C1,C2, ,Cn, причем уравнения этого семейства линий разрешены
относительно y1, ,yn.
Частное решение системы – это такое решение, которое состоит из точек единственности решения задачи Коши.
Решение, получаемое из общего решения при определенных числовых значениях параметров C1,C2, ,Cn, включая , является частным решением.
Особым решением называется такое, в каждой точке которого нарушается единственность решения задачи Коши.
Пример. Найти решение системы
dy |
|
2 |
|
|
|
||
x |
y |
z, |
|||||
|
|
|
|||||
|
|
||||||
dx |
|
x |
|
|
|||
|
|
|
|
(5.9) |
||
|
dz |
|
|
|
|
|
2 |
z, |
где x 0. |
||||
|
|
|||||
dx |
|
|
|
|
||
Решение. Заданная система, представляет собой совокупность зацепляющихся уравнений – одно из них можно решить независимо от другого. Интегрируем второе уравнение, разделяя переменные
|
|
|
dz |
|
|
x C , |
z 0. |
|
|
|
|
|
|
||||
|
|
2 z |
1 |
|
|
|||
|
|
|
|
|
|
|
||
Откуда |
|
|
|
|
|
z (x C )2. |
||
|
|
x C |
(x C ), или |
|||||
|
z |
|||||||
1 |
|
|
1 |
|
1 |
|||
|
|
|
|
|
120 |
|
|
|
Подставим найденную функцию в первое уравнение
dy 2 y C1. dx x
Полученное уравнение является линейным дифференциальным уравнением 1-го порядка. Его решение найдем методом Бернулли: y uv, y u v uv ,
dv |
|
|
2 |
v 0, |
v x2, |
|
|
|||
|
|
|
|
|
|
|||||
|
|
|
|
|||||||
dx |
|
|
x |
|
|
|
|
|
||
du |
|
|
|
|
|
C |
|
|
||
|
|
|
v C , |
u |
1 |
C |
2 |
, |
||
|
|
|
||||||||
|
|
|
1 |
|
x |
|
||||
dx |
|
|
|
|
|
|
|
|||
откуда: |
y C x C |
2 |
x2 |
, |
z (x C )2 |
(x C ) |
- общее решение системы. |
|
1 |
|
|
1 |
1 |
|
Второе уравнение системы (5.9) имеет особое решение z 0. Подставим его в первое уравнение, получим
dy x 2 y. dx x
Решая это уравнение методом Бернулли, получим y x2 C ln x . Итак,
система (5.9) помимо общего решения имеет еще семейство решений
z 0,
y x2 C ln x .
5.3.Интегрирование системы дифференциальных уравнений сведением
кодному уравнению n–го порядка
Впредыдущих разделах 5–й главы рассматривались примеры нахождения решений систем уравнений. Однако используемые там приемы интегрирования не были систематическими, а носили скорее характер ”ad hoc”.
Вэтом разделе будет рассмотрен один из основных методов интегрирования системы дифференциальных уравнений.
Суть используемого здесь метода заключается в следующем: из уравнений системы и из уравнений, получающихся дифференцированием уравнений системы, исключаются все искомые функции, кроме одной. Для определения этой функции получают одно дифференциальное уравнение n–го порядка (в случае системы из n уравнений). Из полученного уравнения находится одна неизвестная функция, а остальные искомые функции находятся из исходных уравнений и уравнений, получившихся в результате их дифференцирования.
1. Общий случай нормальной системы.
Рассмотрим схему приведения системы к одному дифференциальному уравнению n – го порядка в случае системы вида (5.1):
121
|
(k 1,2, ,n). |
(5.10) |
yk fk (x, y1, y2, , yn) |
||
Метод исключения неизвестных реализуется в предположении, что все |
||
функции fk (k 1,2, ,n) имеют непрерывные частные |
производные до |
|
(n 1)- го порядка включительно по всем аргументам x, y1, y2, , yn . |
||
Дифференцируя, например, первое из уравнений (5.10) последовательно |
||
(n 1) раз и подставляя на каждом шаге вместо производных |
|
|
yk (k 1,2, ,n) |
||
их значения, взятые из правых частей уравнений (5.10), мы получим следующие уравнения:
|
|
|
|
|
f1 |
|
|
|
|
|
f1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
f1 |
|
|
|
|
|
f1 |
|
|
|
|
|
f1 |
|
|
|
f1 |
|
|
fn, |
|
|
||||||||||||||
|
|
|
|
x |
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
y |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
y1 |
|
|
|
|
|
y1 |
|
y |
n |
|
yn |
|
|
f1 y |
n |
|
|
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|||||||
или |
y |
F (x, y , y |
2 |
, , y |
n |
) . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
1 |
2 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
y1 |
F2 |
|
|
F2 |
|
y1 |
|
F2 |
|
yn |
|
F2 |
|
|
|
|
F2 |
f1 |
Fn |
fn |
, |
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
x |
|
|
y |
|
|
|
|
y |
n |
|
|
x |
|
|
|
y |
|
y |
n |
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|||||||
или |
y |
|
F (x, y , y |
2 |
, , y |
n |
). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
1 |
3 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
(m) |
|
|
|
|
F |
|
|
|
|
|
|
|
n |
F |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
y1 |
|
|
|
|
|
|
|
|
|
|
|
m 1 |
|
|
|
|
|
m 1 |
fi, |
(m 2,3,...,n 1), |
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
(m) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
i 1 |
|
|
yi |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
или |
y |
F (x, y , y |
2 |
, , y |
n |
) |
|
(m 2, ,n 1) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
1 |
m |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F |
|
|
|
|
|
|
n |
|
F |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y1(n) |
|
|
|
|
n 1 |
|
|
|
n 1 |
fl , |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
yl |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
или |
y(n) F (x, y , , y |
n |
) . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
1 |
n |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
При |
получении |
|
|
|
|
|
|
этих |
|
|
|
|
уравнений |
|
|
|
использовали |
правило |
|||||||||||||||||||||||||||||||||||||||||
дифференцирования |
|
|
|
|
|
сложной |
|
|
|
|
функции |
|
|
многих |
переменных: |
если |
||||||||||||||||||||||||||||||||||||||||||||
u f (x, y1, , yn), |
|
где |
|
|
|
|
y1 g1(x), |
|
|
y2 g2(x), |
|
, |
yn gn(x), функция |
|||||||||||||||||||||||||||||||||||||||||||||||
u f (x, y1, , yn ) дифференцируема в точке |
|
(x, y1, , yn ), а функции |
gk (x) |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(k 1,2, ,n) имеют производную в точке |
|
x, то производная по x |
от сложной |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||
функции u F(x) f (x,g1(x), ,gn(x)) вычисляется по формуле |
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
dF |
|
|
|
|
|
f |
|
|
f |
|
g1(x) |
|
f |
|
gn (x). |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
|
|
|
x |
|
y |
|
y |
n |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Таким образом, получена система (n 1) |
уравнений |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||
122
|
|
dy1 |
f (x, y , , y ), |
|
|
|
||||||||||
|
|
|
|
|
||||||||||||
|
dx |
|
1 |
|
1 |
|
n |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
d2 y1 |
F (x, y , , y |
n |
), |
|
|
|||||||||
|
|
|
|
|||||||||||||
|
|
|
|
2 |
|
1 |
|
|
|
|
|
|
||||
|
dx2 |
|
|
|
|
|
|
|
|
|
|
|
(5.11) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
dn 1y |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
1 |
|
F |
|
(x, y , , y |
|
) |
|||||
|
|
|
|
n |
||||||||||||
|
|
|
n 1 |
1 |
|
|
|
|
|
|
||||||
|
dxn 1 |
|
|
|
|
|
|
|
|
|
|
|
||||
и еще одно уравнение |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
d n y |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
1 |
|
F (x, |
y , , y |
n |
). |
|
(5.12) |
|||||||
|
|
|
|
|||||||||||||
|
dxn |
n |
|
1 |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Дальнейшее продвижение по намеченной схеме заключается в |
||||||||||||||||
исключении из уравнения (5.12) |
функций |
|
|
y1,y2, , yn . Для этого следует |
||||||||||||
использовать систему (5.11), рассматривая ее как систему алгебраических уравнений относительно неизвестных y2, , yn . Требуем, чтобы эта система была разрешима относительно указанных неизвестных. Из математического анализа известно, что для разрешимости такой системы уравнений требуется выполнение условия
|
|
|
|
|
|
|
|
|
f1 |
|
|
|
f1 |
|
|
|
|
f1 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
y2 |
|
|
y3 |
|
|
|
yn |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
F2 |
|
|
|
|
F2 |
|
|
|
|
F2 |
|
|
0. |
|
|
|
|||
|
|
|
|
|
|
|
|
|
y2 |
|
|
y3 |
|
|
|
|
yn |
|
|
|
(5.13) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
Fn 1 |
|
|
Fn 1 |
|
Fn 1 |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
y2 |
|
|
y3 |
|
|
|
|
|
y3 |
|
|
|
|
|
|
||||
|
Предполагая, что условие (5.13) выполняется, выразим из уравнений |
|||||||||||||||||||||||||||
(5.11) |
y |
2 |
, , y |
n |
через |
x, y y |
, , y(n 1) |
и |
|
|
подставим эти |
выражения |
в |
|||||||||||||||
|
|
|
|
|
|
|
|
|
1, |
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|||||
уравнение (5.12). В результате для функции |
|
y1 получим дифференциальное |
||||||||||||||||||||||||||
уравнение n – го порядка: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
y(n) f (x, y |
y |
, , y(n 1) ). |
|
(5.14) |
|||||||||||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
1, |
1 |
|
|
|
1 |
|
|
|
|
|
|||||
|
Решив |
уравнение (5.14) |
получим |
функцию y1(x), |
а |
затем |
ее |
|||||||||||||||||||||
производные |
y |
, , y(n 1) . |
Это |
|
позволит |
найти остальные |
неизвестные |
|||||||||||||||||||||
|
|
|
|
1 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x, y1, y1, , y1(n 1) |
|
|
|||
y2, , yn , |
поскольку они выражаются |
через |
|
|
из системы |
|||||||||||||||||||||||
(5.11).
Можно показать (мы на этом не останавливаемся), что при сделанных предположениях, решение y1(x), найденное из (5.14) а также функции y2, , yn , найденные из (5.11), в совокупности будут являться решением
123
системы уравнений (5.10). И наоборот – решение системы (5.10) будет являться решением уравнения (5.14) и системы (5.11).
Пример. Решить систему уравнений
dy |
|
1 |
|
|
|||
|
|
|
1 |
|
, |
||
|
|
|
|||||
dx |
|
|
z |
||||
|
|
1 |
|
||||
|
dz |
|
. |
||||
|
|
||||||
|
|
y x |
|||||
dx |
|
||||||
Решение. Дифференцируя обе части первого уравнения, получим:
|
|
|
|
|
d |
2 y |
1 |
|
|
dz |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
z2 |
|
|
|
|||||
Из уравнений системы имеем: |
|
dx2 |
|
|
|
dx |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
||||
|
dz |
|
|
1 |
|
|
|
1 |
|
|
dy |
|
||||
|
|
|
|
|
|
, |
|
|
|
|
|
|
|
1 . |
||
|
dx |
|
y x |
|
|
z2 |
|
|
|
|||||||
|
|
|
|
|
|
|
|
dx |
|
|
||||||
Подставим эти выражения в предыдущее уравнение и присоединим выражение для z, найденное из 1-го уравнения:
d |
2 y |
1 dy |
|
2 |
z |
1 |
|
|
||||
|
|
|
|
|
|
1 |
, |
|
|
. |
||
|
|
|
|
|
|
|||||||
dx2 |
y x dx |
|
|
|
1 |
dy |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
dx
Для решения уравнения 2-го порядка представим его в виде:
|
|
|
|
d dy |
|
|
|
|
dy |
|
|
|
|
|
|
|
|
dy |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
dx dx |
|
|
dx |
|
|
, |
где |
|
|
|
dx |
|
|
|
d |
ln |
|
y x |
|
. |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
dy |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
y x |
|
|
|
|
|
y x |
|
|
dx |
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Получим уравнение |
|
|
d |
|
|
dy |
|
|
d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
ln |
1 |
|
ln |
|
y x |
|
. |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
dy |
|
|
|
|
|
|
|
|
dx |
|
dx |
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
Отсюда |
1 C1(y x), |
C1 0. Преобразуем его к виду: |
||||||||||||||||||||||||||||||||||
|
||||||||||||||||||||||||||||||||||||
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
d (y x) C1(y x). dx
Интегрируя, имеем: y x C2eC1x, C1 0, C2 , . Находим вторую неизвестную функцию:
z |
1 |
|
|
1 |
|
1 |
|
|
. |
|
|
dy |
|
C (y x) |
|
|
|
||||
|
|
|
|
|
C C |
e |
C1x |
|||
|
1 |
|
|
|
1 |
|
|
|
||
|
|
|
|
|
|
1 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
dx
Таким образом, решение системы имеет вид:
C x |
|
1 |
|
C x |
|
y x C2e 1 , |
z |
|
e |
1 . |
|
C1C2 |
|||||
|
|
|
|
||
|
124 |
|
|
|