Материал: Обыкновенные дифференциальные уравнения с приложениями к задачам механики, физики, термодинамики и экологии. Ряжских В.И., Бырдин А.П

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

где 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,bM .

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

 

 

 

Источник: https://studfile.net/preview/16568742/