Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Математика. Дифференциальные уравнения

.pdf
Скачиваний:
5
Добавлен:
12.11.2023
Размер:
833.93 Кб
Скачать

б) Найдем общее решение соответствующего однородного уравнения. Составляем характеристическое уравнение:

k4 + 2k2 +1 = 0 ,

 

 

 

 

(

)

 

 

 

 

 

 

 

 

k2 +1 2 = 0 .

 

 

 

Тогда, k1,2 = ±i, k3,4

= ±i .

 

 

 

Согласно п. 6 §2

 

 

 

 

 

 

y = e0 x cos x = cos x, y

2

= e0 x sin x =sin x, y = x cos x,

y

4

= xsin x .

1

 

 

 

3

 

 

 

y =C1 cos x +C2 sin x +C3 x cos x +C4 xsin x

 

 

 

 

 

 

 

 

или

 

 

 

 

 

=C1 cos x +C2 sin x + x(C3 cos x +C4 sin x) .

 

 

 

 

y

 

 

 

f (x) =5sin x xcos x

или f (x) = e0 x (5sin x x cos x) .

Так как α ±iβ = 0 ±i

являются корнем характеристического

уравнения, кратности λ = 2 (см. замечание 1, гл. II, §3, п. 3), то

y* = ( Ax + B)sin x +(Cx + D)cos x x2 .

 

 

71

Глава IV

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

§1. Общие понятия

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

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

Обычно приходится иметь дело не с произвольными системами дифференциальных уравнений, а с так называемыми нормальными системами.

 

 

Основные формулы

 

 

Определения и замечания

 

dy1

 

 

= f1 (x, y1 , y2 ,..., yn )

(4.1)

 

Система дифференциальных уравне-

 

 

 

 

 

 

 

 

ний первого порядка, относительно

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

искомых функций

y1 (x), y2 (x), ...,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy2

 

 

= f

2 (x, y1 , y2 ,..., yn )

 

 

 

 

 

1.

 

 

 

 

 

 

 

yn (x) , причем правые части уравне-

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

ний не содержат производных иско-

 

 

dyn

 

 

 

 

 

 

 

 

 

 

мых функций, называется системой

 

 

= f

 

(x, y

, y

 

,..., y

 

)

 

 

 

n

2

n

 

 

в нормальной

форме или нор-

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

мальной системой.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1 = ϕ1 (x,C1,C2 ,...,Cn )

 

(4.2)

 

(4.2) – общее решение нормальной

 

 

 

 

 

 

= ϕ2

(x,C1,C2 ,...,Cn )

 

 

системы (4.1),

где

C1 ,C2 ,...,Cn

2.

y2

 

 

произвольные постоянные.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ϕn

(x,C1 ,C2 ,...,Cn )

 

 

 

 

 

 

yn

 

 

 

 

 

72

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Основные формулы

 

Определения и замечания

y1 (x0 ) = y10

(4.3)

(4.3) – начальные условия, при по-

 

(x0 ) = y20

 

мощи которых из общего решения

y2

 

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

 

 

 

 

...

 

 

 

yn (x0 ) = yn0

 

 

ϕ1 ( x0 ,C1,C2 ,...,Cn ) = y10

(4.4)

Подставляя начальные условия в об-

 

( x0 ,C1,C2 ,...,Cn ) = y20

 

щее решение (4.2), получаем систему

ϕ2

 

алгебраических уравнений для опре-

 

 

 

...

 

 

деления произвольных постоянных.

ϕn (x0 ,C1 ,C2 ,...,Cn ) = yn0

 

 

§2. Метод исключения неизвестных

Для интегрирования системы (4.1) можно применить метод, с помощью которого данная система, содержащая n уравнений относительно n искомых функций, сводится к одному уравнению n-го порядка относительно одной неизвестной функции. Этот метод называется методом исключения неизвестных.

Покажем его применение на примере. Пусть дана система уравнений:

dx

=3x

+ y,

 

dt

 

 

 

 

 

dy

= −x

+ y.

 

 

 

dt

 

 

 

Из всех решений выделить то решение, которое удовлетворяет

условию: x(0) =1, y (0) = 0.

 

 

 

Решение

 

 

 

Дифференцируя первое из уравнений системы по t,

находим:

d 2 x =3 dx

+ dy

(4.5)

dt2

dt

dt

 

 

 

 

73

Подставляя в это равенство выражение dydt из второго уравне-

ния системы, получим:

d 2 x

=3 dx

x + y .

(4.6)

dt2

dt

 

 

Заменяя, наконец, функцию y ее выражением из первого уравнения системы

y = dx

3x ,

(4.7)

dt

 

 

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

d 2 x

=3 dx

x + dx 3x

dt2

 

dt

dt

или

 

 

 

 

d 2 x

4 dx + 4x = 0.

 

 

 

dt2

dt

Составляем характеристическое уравнение

 

 

k2 4k + 4 = 0,

k = 2, x = e2t , x =te2t .

1,2

 

1

2

Общее решение:

x =C1e2t +C2te2t

или

x = e2t (C1 +C2t ) .

Дифференцируя равенство (4.9), находим:

dxdt = 2e2t (C1 +C2t) +C2e2t

(4.8)

(4.9)

74

 

или

 

dx

= e2t (2C1 + 2C2t +C2 ) .

(4.10)

dt

 

 

Подставляя выражения (4.9) и (4.10) в равенство (4.7), получим:

y = e2t (2C1 + 2C2t +C2 ) 3e2t (C1 +C2t )

или

y = −C1e2t +C2e2t (1t) .

Общее решение системы:

 

 

2t

 

 

2t

 

 

x =C1e

+C2te

,

 

 

 

 

y = −C e2t +C

e2t

(1t ).

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

Решим задачу Коши. Пусть x(0) =1 . Тогда

C1e0 +C2 0 e0 =1,

следовательно,

C1 =1 .

Пусть y (0) = 0. Тогда:

C1e0 +C2e0 1 = 0,

C1 +C2 = 0,

C2 =1 .

Искомое решение системы

 

2t

(1

+t )

 

x = e

 

 

.

 

 

 

 

y = −te2t

 

 

 

 

 

 

 

(4.11)

75

ПРИЛОЖЕНИЯ

Задача 1. Найти общий интеграл дифференциального уравнения.

1.01. ( x +1)3 dy ( y 2)2 dx = 0

1.02. sec2 x sec ydx +ctg x sin ydy = 0 1.03. x 1+ y2 + yy2 + x2 = 0

1.04. 3x2 ydy 2xy2dx = 6xdx 6 ydy

1.05. ( xy + x ) y′− y = 0 1.06. ye2 xdx +(e2 x +5)dy = 0

1.07.

2x+y 3x2 y y′= 0

1.08.

x 4 + y2 dx + y 3 + x2 dy = 0

1.09.2(1+ex ) yy′= ex

1.10.4xdx 3ydy =3x2 ydy 2xy2 dx

1.11. yy

1x2

+1 = 0

 

1y2

 

1.12.3ex sin ydx =(ex 1)sec ydy

1.13.(1+ x2 ) y3dx ( y2 1)x3dy = 0

1.14.5 + y2 dx ydy = x2 ydy

1.15. 2x + 2xy2 + 2 x2 y′= 0

1.16.y ln y + xy′= 0

1.17.ysin2 x ln y + y = 0

1.18.( xy x)2 dy + y(1x)dx = 0

76

1.19.y = 2(1+ x2 y)+ xy

1.20.( y 1) xdy y2 ln xdx = 0

1.21.(4 +ex ) y′= 2 yex

1.22.x + xy2 + 1x2 y′= 0

1.23.2 ydy 6xdx =3xy2dx 2 yx2 dy

1.24.y (1+ln y) + xy′= 0

1.25.y3 ln2 xdx ( y 1) xdy = 0

1.26.2x2 yy′+ y2 = 2

1.27.5 + y2 dx + 4(x2 y + y)dy = 0

1.28.ycos x ( y +1)sin x = 0

1.29. 2(x2 y + y)dy + 3 + y2 dx = 0

1.30.4 + y2 dx ydy = x2 ydy

Задача 2. Найти общий интеграл дифференциального уравнения.

2.01. (x2 + y2 )dx 2xydy = 0 2.02. y xy′= y ln xy

2.03. xdy ydx = x tg xy dx

2.04.3y′= y2 +8 y + 4

x2 x

2.05. xy′ =3 2x2 + y2 + y

2.06. y xy′= x + yy

77

2.07. y = x(y′− x ey )

2.08. ( y2 3x2 )dy + 2xydx = 0

2.09.xy′+ y ln xy 1 = 0

2.10.ydx +(2 xy x)dy = 0

2.11.y(2x2 2xy) = x2 + 2xy y2

2.12.4 y′= y2 +10 y +5

x2 x

2.13.x3 y′= xy2 x2 y

2.14.xy′= y + xsin xy

2.15.xy′− y(ln y ln x) = 0

2.16. ycos

y

y

 

cos

y

 

+1 = 0

x

 

x

 

 

 

 

x

 

 

 

 

y = xy′− xe

y

 

 

 

 

 

 

2.17.

x

 

 

 

 

 

 

2.18.

y (2

x

y )dx + xdy = 0

2.19. y xy′= 2(x+ yy)

2.20.

y′=

x2 + xy 3y2

 

 

 

 

 

 

 

 

x2 4xy

 

2.21.

y xy′= xsec

y

 

 

x

 

 

 

 

 

 

 

 

 

 

 

2.22.

xy′= 3x3 +10 yx2

 

 

 

 

2 y2 +5x2

 

2.23. ( x y) y′= x + y

78

2.24.xy′− y =(x + y)ln x +x y

2.25.dydx xy (1+ln y ln x) = 0

2.26. xy′− y =3 x2 + y2

2.27.x(2 y2 + 4x2 ) y′=3y3 +8yx2

2.28.y(2x y) = x + 2 y

2.29.2 y′ = y2 +8 y +8

x2 x

2.30.xy′− y = x tg xy

Задача 3. Найти решение задачи Коши.

3.01.

y′+

2xy

 

=

2x2

 

1+ x2

 

 

 

 

 

 

 

 

1+ x2

3.02.

ysin x y cos x =1

3.03.

xy′+ y + xex2

= 0

 

3.04.

xy′− y x3 = 0

 

 

 

3.05.

y′−

y

 

= tg

 

x

 

sin x

 

 

 

 

 

 

2

 

 

3.06.

x( y′− y) = ex

 

 

 

3.07.

y′+

2xy

 

= 2

 

 

 

 

1 x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.08.

xy′+( x +1) y =3x2ex

3.09.

y′+ 2xy = xex2 sin x

y(0) = 23

 

π

=

1

y

4

 

 

2

 

 

 

y(1) = 12 е1 y(2) = 4

y π = 02

y(1) = 0

y(0) = 2

y(1) = 0 y(0) =1

79

3.10.y′− x2+y1 = ex (x +1)2

3.11.y′− y ctgx = 2xsin x

3.12.y′+ y tgx = cos2 x

3.13.xy′+ y = ln x +1

3.14.y′+ 2 yx1 = x2

3.15.(x2 1) y′− xy = x3 x

3.16.(xy′−1)ln x = 2 y

3.17.y = x( y′− xcos x)

3.18.( x + 2) y′− y =(x2 + 2x)(x + 2)

3.19.y′− 2xx2 5 y =5

3.20.xy′+ y = xsin x

3.21.( x +1) y′+ y = x3 + x2

3.22.y′+ 2yx = x2

3.23.y′− y cos x = sin 2x

3.24.y′− x2+y1 =(x +1)3

3.25.y′− 2xy = −lnxx

3.26.y′− y ctgx = sin1 x

80

y(0) =1

 

π

= 0

y

2

 

 

 

 

 

 

π

=

1

y

4

 

2

 

 

 

y(1) = 0

 

y(2) =1

 

y(

2 ) =1

y(e) = 0

y π = 02

y(1) = 32 y(2) = 4

y(π) = π1

y(0) = 0 y(1) =1

y(0) = −1 y(0) = 12 y(1) =1

y π =12