Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Математика_Семестр3_МетодПособие.pdf
Скачиваний:
157
Добавлен:
05.06.2015
Размер:
1.75 Mб
Скачать

 

 

 

 

 

 

 

 

 

 

 

 

14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u = a1x + b1y.

 

 

 

 

 

 

 

Пример. Решить уравнение y

=

 

x + y +

2

.

 

 

 

 

 

 

 

 

x +1

 

 

 

 

 

 

В этом уравнении a1 = b1 =1,

 

c1 = 2,

a2

 

=1, b2

= 0,

c2

=1. Поэтому

δ =

 

1

1

 

0 . Полагая

x = u +α,

 

 

y = v + β,

находим α

и

β из системы

 

 

 

 

 

 

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

системы уравнений:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

+ β + 2 = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α +1= 0.

 

 

 

 

 

 

 

 

Следовательно, α = −1,

 

β = −1,

 

 

и формулы перехода от старых перемен-

ных к новым и обратно примут вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x = u 1,

y = v 1,

 

 

 

 

 

 

 

 

 

 

 

u = x +1,

v = y +1.

 

 

 

 

 

 

 

В результате уравнение приводится к однородному

 

 

 

 

 

 

 

dv

= u 1+ v 1+ 2

= u + v =1+

v

.

 

 

 

 

 

 

 

du

 

u 1+1

 

 

v

 

u

u

 

 

 

 

 

 

Полагая далее t = v / u,

 

v = tu,

=t u + t,

приходим к уравнению

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с разделяющимися переменными относительно функции t:

 

 

 

 

 

 

 

 

 

 

du

,

 

t = lnu + lnC = lnCu.

 

 

 

 

 

 

t u +t =1+t, d t =

u

 

 

Возвращаясь к старой переменной, получим

 

 

 

 

 

 

 

 

 

 

v

= lnCu,

 

 

y +1

= lnC(x +1),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x +1

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

y = (x +1) lnC(x +1) 1.

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

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

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

(1.12)

15

содержащее искомую функцию и ее производную в первой степени. Функции P(x) и Q(x) предполагаются непрерывными.

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

Рассмотрим интегрирование этого уравнения методом вариации про-

извольной постоянной (методом Лагранжа). В соответствии с этим мет о-

дом сначала ищется решение соответствующего линейного однородного

уравнения: y′+ P(x) y = 0 .

Разделяя в нем переменные, получим его общее

решение в виде

 

 

 

 

dy

= −P(x) y,

dy

= −P(x)dx, ln y = −P(x)dx +lnC,

 

dx

 

y

 

 

 

 

 

y =C e−∫P(x)dx .

(1.13)

Далее ищется общее решение исходного уравнения с ненулевой правой частью в виде (1.13), но произвольная постоянная в (1.13) заменяется неизвестной функцией:

y = v(x) e−∫P(x)dx .

(1.14)

Подставляя (1.14) в (1.12), получаем дифференциальное уравнение с разделяющимися переменными относительно варьируемой функции v(x), интегрируя которое, находим эту функцию.

В результате, общее решение уравнения (1.12) может быть представлено в виде

y = e−∫P(x)dx [ Q(x)eP(x)dx dx + C].

(1.15)

 

Формула (1.15) дает общее решение линейного дифференциального

уравнения первого порядка в форме Коши.

 

Пример. Решить задачу Коши для уравнения y′− 4y = x

при началь-

ном условии y(0)=1.

 

Находим сначала общее решение линейного однородного уравнения y′− 4y = 0 .

Оно имеет вид

16

y = C e−∫P(x)dx = C e−∫−4dx = C e4x .

Заменяем произвольную постоянную C в этом решении неизвестной функцией v(x):

y = v(x)e4x .

Вычисляем производную y

 

 

 

 

 

 

y

 

 

 

 

 

 

4x

+ 4v(x)e

4x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= v (x)e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и подставляя

y и

y

в исходное уравнение, получим

 

 

 

 

 

 

 

 

4x

+

4v(x)e

4x

4v(x)e

4x

 

 

 

 

 

 

 

 

4x

= x,

 

 

v (x)e

 

 

 

 

 

= x, v (x)e

 

 

 

 

4x

, v(x)

 

= xe

4x

dx + C

= −

1

e

4x

 

+

1

 

+ C.

 

v

(x) = xe

 

 

 

 

4

 

 

x

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение уравнения примет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4x

 

 

 

 

 

1

 

4x

 

 

 

1

 

 

4x

 

 

4x

 

 

 

x +1/ 4

 

 

y =ν(x)e

 

=

 

 

e

 

x

+

 

+ C e

 

 

=Ce

 

 

 

 

.

 

 

 

4

 

4

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Находим произвольную постоянную C из начального условия:

при x = 0

 

 

 

 

C

 

 

1

 

=1,

 

C =

17.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, решение задачи Коши будет

y =161 [17e4x (4x +1)].

Решение линейного дифференциального уравнения (1.12) может быть также получено, если искомую функцию представить в виде произведения двух произвольных функций (метод Бернулли [4]):

 

 

y =u (x) v(x).

(1.16)

Тогда

 

y

=u v + v u.

(1.17)

 

 

 

 

Подставляя (1.16) и (1.17) в (1.12), получим

(1.18)

[u

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

 

 

 

 

 

17

Функцию u(x) подбираем так, чтобы она была одним из решений урав-

нения

u′+ P(x)u = 0 .

 

После разделения переменных получим

 

u = e−∫ P(x)dx .

(1.19)

Тогда уравнение (1.18) примет вид

ve−∫P(x)dx = Q(x).

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

dv =Q(x)eP(x)dx dx.

Интегрируя это уравнение с разделенными переменными, находим функцию v:

v = Q(x)eP(x)dx dx + C .

(1.20)

Подставляя (1.19) и (1.20) в (1.16), получим общее решение уравнения

(1.12) в виде (1.15).

1.7. Уравнение Бернулли

Уравнением Бернулли называется уравнение вида

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

(1.21)

где n 0;1.

Уравнение Бернулли является нелинейным, но оно приводится к линейному следующим преобразованием:

1) Обе части уравнения умножаются на yn , тогда

yn y′+ P(x) y1n = Q(x).

2) Далее применяется подстановка

z = y1n .

Тогда по правилу дифференцирования сложной функции получим

dz

= (1n) y

n

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

y

n

y′ =

1 dz

dx

 

 

 

dx .

 

 

1n

В результате уравнение становится линейным относительно функции z:

18

 

1

 

1n z′ + P(x)z =Q(x).

(1.22)

Уравнение (1.22) может быть решено методом вариации произвольной постоянной или методом Бернулли.

Пример. Решить уравнение

 

 

 

 

 

 

 

 

 

 

 

 

y′+

y

lnx

y

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

y2 :

 

 

 

y2 y′+

y1

=

lnx

.

 

 

 

 

 

 

 

 

 

x

x

 

 

Положим

z = y

1

,

 

 

 

z

= −y

2

 

и

уравнение преобразуется в линей-

 

 

 

 

 

 

 

 

 

 

y ,

ное:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z′−

z

 

 

 

 

 

lnx

.

 

 

 

 

 

 

 

 

(1.23)

 

 

 

 

 

 

 

 

 

 

 

 

 

= − x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

Находим сначала решение соответствующего линейного однородного

уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

z

= 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

=

z

,

 

 

 

ln

 

z

 

= ln

 

x

 

+ ln C,

 

 

z = Cx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение неоднородного уравнения (1.23) отыскиваем в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z = v(x) x,

 

 

 

 

 

 

 

 

 

 

 

 

тогда

 

 

 

 

 

 

 

 

 

 

 

z

 

= v (x) x + v(x).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dv(x) x + v(x) v(x) x = −lnx ,

 

 

 

 

 

dv(x) x = −lnx

,

 

dv(x) = −

lnx

dx +C.

 

 

 

 

 

 

2

dx

 

 

x

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

x

 

 

 

 

 

 

 

x

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v(x) =

lnx

+

1

+ C,

 

 

 

 

 

 

 

 

z = v(x) x

lnx

+

1

 

 

 

 

x

x

 

 

 

 

 

 

 

 

=

x

 

x

+ C x,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

поэтому общее решение исходного уравнения будет иметь вид