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

2006-matematicheskiy-analiz-differencial-nye-uravneniya-2mb

.pdf
Скачиваний:
326
Добавлен:
23.03.2016
Размер:
2.95 Mб
Скачать

Тогда

x = u + 1,

dx = du,

 

 

y = v + 1

dy = dv;

y′ = dy = dv . dx du

Исходное уравнение примет вид

dv = u + 1+ 3(v + 1) 4 ; du 3(u + 1) (v + 1) 2

dv = u + 3v . du 3u v

Получили однородное уравнение. Сделаем еще одну замену

t = v v = tu . u

Тогда

v′ = tu + t.

Значит, уравнение примет вид

tu + t = u + 3tu ; 3u tu

tu = u + 3tu t; 3u tu

tu = 1+ 3t t; 3 t

tu =

1+ 3t 3t + t2

;

 

 

 

 

 

 

 

 

3 t

 

 

 

u

dt

=

1+ t2

;

 

du

3 t

 

 

 

 

 

 

udt =

1+ t2

 

du.

3 t

 

 

 

 

 

 

 

31

Получили уравнение с разделяющимися переменными. Разделим обе части уравнения на u 0:

 

 

 

 

 

 

 

 

 

 

 

dt =

t2 + 1

 

 

 

 

du

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 t

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

+ 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее разделим на

 

 

 

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 t

dt =

 

 

du

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2 + 1

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрируем обе части полученного уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

du

= ln

 

 

u

 

+ с;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 t

dt = −

 

t 3

dt = −

1

 

 

d (t2 + 1)

+ 3

 

 

 

 

dt

= −

1

ln(t2 + 1) + 3arctgt + с.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2 + 1

t2 + 1

2

 

 

t2 + 1

 

 

 

 

 

 

 

 

 

 

t2 + 1

2

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3arctgt

2

ln(t

 

+ 1) = ln

u

 

+ с;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6arctgt = ln(t

2

 

+

1) + 2ln

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 2с;

 

 

 

 

 

 

6arctgt = ln (t2 + 1) + ln

 

u2

+ с;

 

 

 

 

 

 

 

6arctgt = ln (u2 (t2 + 1)) + c.

 

Подставим t =

v

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6arctg

 

= ln u

2

v

 

 

+

1

+ c;

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6arctg v = ln (v2 + u2 ) + c. u

Подставив u = x – 1, v = y – 1, получим общий интеграл данного уравнения

6arctg y 1 = ln (( y 1)2 + (x 1)2 ) + c . x 1

32

3.5. Линейные дифференциальные уравнения

Линейным дифференциальным уравнением первого порядка назы-

вается уравнение вида

y+ P(x)y = Q(x),

где P(x), Q(x) – непрерывные на некотором интервале (α, β) извест- ные функции.

Заметим, что для линейного дифференциального уравнения пер- вого порядка выполнены условия теоремы существования и единст- венности решения задачи Коши, т.е. данное уравнение всегда имеет единственное решение, удовлетворяющее заданному начальному ус-

ловию y(х0) = у0, где х0 (α, β), у0 R.

Если Q(x) = 0, то уравнение называется однородным. Если Q(x) 0, то уравнение называется неоднородным.

Общее решение однородного уравнения получается разделением переменных:

у ′ + P(x)y = 0;

dy + P(x) y = 0 . dx

Умножим обе части уравнения на dx:

dy + P(x)ydx = 0.

Разделим обе части уравнения на y (при делении на у можем поте- рять решение у = 0):

dy + P(x)dx = 0; y

dy = − P(x)dx . y

Интегрируя левую часть уравнения по y, а правую часть по x по- лучаем

 

 

 

 

 

 

 

 

 

ln

 

y

 

 

= −P(x)dx + с.

 

 

Выразим y:

 

 

 

 

 

 

 

 

 

P( x)dx+ с

 

 

 

 

 

 

 

= e

 

 

 

 

 

y

 

;

33

y = eс eP( x)dx .

Пусть eс = с , тогда общее решение однородного уравнения

y = сeP( x)dx .

Заметим, что решение у = 0 мы не теряем, так как оно получается при с = 0.

Для решения неоднородного уравнения применяем метод вариа- ции произвольной постоянной, который заключается в следующем. Сначала находим общее решение соответствующего однородного линейного уравнения. Затем, заменив постоянную с на функцию c(x), ищем частное решение неоднородного уравнения в виде

yч.н = с(x)eP( x)dx .

Подставив в исходное уравнение yч.н и yч.н, находим функцию с(x). Общее решение неоднородного уравнения будет равно сумме общего решения однородного уравнения уо.о и частного решения не- однородного уравнения уч.н:

у = уо.о + уч.н.

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

y(x) = u(x)v(x) ,

где u(x) и v(x) – неизвестные функции переменной x. Тогда имеем

y'(x) = u'(x) v(x) + u(x) v'(x).

Подставляя значения y и y' в уравнение y+ P(x)y = Q(x), получаем u'(x) v(x) + u(x) v'(x)+ P(x) u(x) v(x)= Q(x).

Сгруппируем слагаемые:

u'(x) v(x) + u(x)(v'(x)+ P(x) v(x))= Q(x).

Подберем функцию v(x) так, чтобы v'(x)+ P(x) v(x)= 0. Относи- тельно v(x) имеем линейное однородное уравнение, следовательно, как было показано ранее

34

v = сeP( x)dx .

Затем найдем функцию u(x) из условия u'(x) v(x) = Q(x).

Замечание

Иногда удобнее рассматривать уравнение x+ P(y)x = Q(y) – ли- нейное относительно переменной y.

Пример 3.13. Решить задачу Коши: y '

y

= x ln x , y(e) =

e2

.

 

x ln x

2

 

Решение

Данное уравнение линейное неоднородное. Рассмотрим два спо- соба решения линейного уравнения: метод вариации произвольной постоянной и метод Бернулли.

Метод вариации произвольной постоянной

Рассмотрим однородное уравнение y '

y

= 0,

 

x ln x

соответствующее данному неоднородному уравнению. Это урав- нение с разделяющимися переменными.

Подставим y′ = dy : dx

dy = y . dx x ln x

Умножим обе части уравнения на dx и разделим на y:

dy

=

dx

.

 

 

y

x ln x

Интегрируем обе части полученного уравнения:

 

 

dy

= ln

 

 

y

 

+ c;

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

=

d ln x

= ln

 

ln x

 

+ c.

 

 

 

 

 

Тогда

x ln x

 

 

 

ln x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

y

= ln

ln x

+ c;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

35

ln y = ln ln x + ln c ; y = c ln x .

Итак, общее решение однородного уравнения

yо.о = c ln x .

Найдем частное решение неоднородного уравнения методом ва- риации произвольной постоянной. Решение будем искать в виде

yч.н = c(x)ln x .

Тогда

y′ = c(x)ln x + c(x) 1 .

ч.н

x

Подставим уч.н и y

в исходное уравнение:

ч.н

 

 

 

 

 

с′(x)ln x + с(x)

1

с(x)ln x

= x ln x;

x

 

 

 

 

x ln x

 

с′(x)ln x = x ln x;

с′(x) = x.

Тогда

с(x) = xdx = x2 + с. 2

При с = 0 получим частное решение неоднородного уравнения

 

2

 

уч.н

= c(x)ln x =

х

ln x .

 

 

2

 

Так как у = уо.о + уч.н, то общее решение неоднородного уравнения

у = с ln x + x2 ln x . 2

Найдем частное решение, удовлетворяющее заданному начальному

2

условию у(е) = е . Для этого найдем константу с, подставив x = e и

2

2

у = е в найденное общее решение неоднородного уравнения:

2

36

 

 

 

 

e2

= с ln e +

e2

ln e;

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

e2

 

= с +

e2

с = 0.

 

 

2

 

 

 

 

 

2

 

 

 

Тогда y =

x2

ln x решение задачи Коши.

 

2

 

 

 

 

 

 

 

 

 

Метод Бернулли

Решение ищем в виде произведения двух функций y(x) = u(x)v(x). Тогда

y'(x) = u'(x) v(x) + u(x) v'(x).

Подставляя значения y и y' в уравнение y '

y

= x ln x , получаем

x ln x

 

 

 

 

 

 

 

 

 

uv + uv′ −

uv

 

= x ln x.

 

 

 

 

 

 

 

 

x ln x

 

 

 

Сгруппируем слагаемые:

 

 

 

 

 

 

 

 

 

 

v

 

 

 

uv + u v′ −

 

 

 

 

 

= x ln x.

 

 

 

 

 

 

 

 

 

x ln x

 

 

 

Подберем функцию v(x) так, чтобы

 

 

 

v′ −

 

v

 

= 0 .

 

 

 

 

 

 

 

 

 

 

 

 

x ln x

Это уравнение с разделяющимися переменными. Разделим пере- менные:

dv

=

 

v

;

 

 

 

dx

 

x ln x

dv

=

 

dx

.

 

 

v

 

x ln x

Интегрируем обе части полученного уравнения:

 

 

dv

= ln

 

 

v

 

+ c;

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

dx

=

d ln x

= ln

 

ln x

 

+ c.

 

 

 

 

 

x ln x

 

 

ln x

 

 

 

 

37

Тогда

ln v = ln ln x + c; ln v = ln ln x + ln c ; v = c ln x .

Так как нам достаточно найти какое-нибудь одно ненулевое ре- шение уравнения, то при с = 1 получим v = ln x .

Найдем функцию u(x) из условия uv = x ln x . Подставим v = ln x :

uln x = x ln x; u′ = x;

u = xdx =

x2

+ c.

2

 

 

 

 

 

Тогда

 

 

 

 

 

 

2

 

 

y = uv = ln x

x

 

 

+ c .

 

 

 

 

2

 

 

Получили общее решение исходного неоднородного уравнения

y = x2 ln x + c ln x .

2

Подставив начальные условия, получим с = 0. Итак, решение задачи Коши:

y =

x2

ln x .

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

Пример 3.14. Решить задачу Коши: y′ − 3y

 

1

 

= −

5arctgx

e3arctgx ,

x2

+ 1

 

 

 

 

 

x2 + 1

y(0) = 2.

Решение

Данное уравнение линейное неоднородное. Рассмотрим одно- родное уравнение

1

 

y′ − 3y

 

= 0,

x2 + 1

соответствующее данному неоднородному уравнению. Это уравне-

ние с разделяющимися переменными. Подставим y′ = dy : dx

38

dy

=

3y

.

 

 

dx

x2 + 1

Умножим обе части уравнения на dx и разделим на y:

dy

=

3dx

.

y

x2 + 1

 

Интегрируем обе части полученного уравнения:

dyy = ln y + c;

3dx = +

x2 + 1 3arctg x c.

Тогда

ln y = 3arctgx + c;

y = ec+3arctgx = ec e3arctgx .

Получили общее решение однородного уравнения

yo.o = ce3arctgx .

Найдем частное решение неоднородного уравнения методом ва- риации произвольной постоянной. Решение будем искать в виде

 

 

 

yч.н

= с(x)e3arctgx .

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

y

= с′(x)e3arctgx + c(x)e3arctgx

3

 

.

 

 

 

 

 

 

 

1+ x2

ч.н

 

 

 

 

 

 

 

 

Подставим уч.н и y

 

в исходное уравнение:

 

 

 

 

 

ч.н

 

 

 

 

 

 

 

 

 

с′(x)e3arctgx + c(x)e3arctgx

 

3

3с(x)e3arctgx

1

 

= −

5arctgx

e3arctgx ;

1+ x2

 

 

 

 

1+ x2

 

 

x2 + 1

(x)e3arctgx = − 5arctgx e3arctgx ;

с

x2 + 1

с′(x) = − 5arctgx . x2 + 1

39

Тогда

 

5arctgx

 

arctgx = t;

 

 

 

t2

5arctg2 x

 

 

 

 

 

 

 

с(x) = −

 

 

 

dx =

dx

 

= −5tdt = −5

 

+ c = −

 

+ c.

 

2

+ 1

= dt

 

 

 

x

 

 

 

2

 

2

 

 

 

 

1+ x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При с = 0 получим частное решение

 

 

 

 

 

 

yч.н = с(x)e3arctgx

= −

5arctg2 x

e3arctgx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

Так как у = уо.о + уч.н, то общее решение неоднородного уравнения

y = ce3arctgx 5 e3arctgxarctg2 x .

2

Найдем частное решение, удовлетворяющее начальному условию у(0) = 2. Для этого найдем константу с, подставив x = 0 и у = 2 в най- денное общее решение неоднородного уравнения:

2 = ce0 5 arctg2 0 e0 c = 2. 2

Итак, решение задачи Коши:

y = 2e3arctgx 5 e3arctgxarctg2 x .

2

Пример 3.15. Найти общее решение уравнения

y′ +

2y

=

4ln(2x 5)

.

 

2x 3

 

2x 3

 

Решение

Данное уравнение линейное неоднородное. Рассмотрим одно- родное уравнение

y′ +

2y

= 0,

 

 

2x 3

соответствующее данному неоднородному уравнению. Это уравне-

ние с разделяющимися переменными. Подставим y′ = dy : dx

40