Материал: Обыкновенные дифференциальные уравнения

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

1) y dp

− p − 4y = 0;

dp

= 4 +

p

;

dy

y

dy

 

 

 

Для решения полученного дифференциального уравнения произведем замену переменной: u = py .

 

 

 

u + du y = 4 +u;

du = 4 dy

;

 

 

 

 

 

 

 

dy;

dy

 

y

 

 

 

 

 

∫

du = 4

∫

u = 4ln

 

y

 

+ 4ln C ;

u = 4ln

 

C y

 

;

 

 

 

 

 

y

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p= 4yln C1 y;

Сучетом того, что p = dydx , получаем:

 

dy

= 4yln

 

 

 

C1 y

 

;

 

 

 

 

 

∫

 

 

dy

 

 

 

= ∫dx;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

4yln

 

C y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

d(ln

 

C1 y

 

)

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

= 4 ∫

 

 

 

 

 

 

 

 

 

=

4 ln

 

ln

 

C1 y

 

+C2 ;

 

 

ln

 

C y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общий интеграл имеет вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

ln

 

C1 y

 

 

 

= 4x +C;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) p = 0;

y′ = 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = C;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, получили два общих решения.

Линейные дифференциальные уравнения высших порядков.

Определение. Линейным дифференциальным уравнением n – го порядка

называется любое уравнение первой степени относительно функции у и ее производных y′, y′′,..., y(n) вида:

p0 y(n) + p1 y(n−1) + p2 y(n−2) +... + pn−1 y′+ pn y = f (x);

где p0, p1, …,pn – функции от х или постоянные величины, причем p0 ≠ 0.

Левую часть этого уравнения обозначим L(y).

p0 y(n) + p1 y(n−1) + p2 y(n−2) +... + pn−1 y′+ pn y = L(y);

Определение. Если f(x) = 0, то уравнение L(y) = 0 называется линейным

однородным уравнением, если f(x) ≠ 0, то уравнение L(y) = f(x) называется линейным неоднородным уравнением, если все коэффициенты p0, p1, p2, … pn – постоянные числа, то уравнение L(y) = f(x) называется линейным дифференциальным уравнением высшего порядка с постоянными коэффициентами.

Отметим одно важное свойство линейных уравнений высших порядков, которое отличает их от нелинейных. Для нелинейных уравнений частный интеграл находится из общего, а для линейных – наоборот, общий интеграл составляется из частных. Линейные

уравнения представляют собой наиболее изученный класс дифференциальных уравнений высших порядков. Это объясняется сравнительной простотой нахождения решения. Если при решении каких – либо практических задач требуется решить нелинейное дифференциальноеуравнение,точастоприменяютсяприближенныеметоды,позволяющие заменить такое уравнение “близким” к нему линейным.

Рассмотрим способы интегрирования некоторых типов линейных дифференциальных уравнений высших порядков.

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

Рассмотрим уравнение вида p0 y(n) + p1 y(n−1) + p2 y(n−2) +... + pn−1 y′+ pn y = 0

Определение. Выражение p0 y(n) + p1 y(n−1) + p2 y(n−2) +... + pn−1 y′+ pn y = L(y)

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

Линейный дифференциальный оператор обладает следующими свойствами:

1)L(Cy) = CL(y);

2)L(y1 + y2 ) = L(y1 ) + L(y2 );

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

1)Если функция у1 является решением уравнения, то функция Су1, где С – постоянное число, также является его решением.

2)Если функции у1 и у2 являются решениями уравнения, то у1 +у2 также является его решением.

Структура общего решения.

Определение. Фундаментальной системой решений линейного однородного дифференциального уравнения n –го порядка на интервале (a, b) называется всякая система n линейно независимых на этом интервале решений уравнения.

Определение. Если из функций yi составить определитель n – го порядка

 

y1

y2

...

yn

 

 

 

W =

y1′

y2′

...

yn′

,

 

...

...

...

...

 

 

y(n−1)

y(n−1)

... y(n−1)

 

 

1

2

 

n

 

то этот определитель называется определителем Вронского.

( Юзеф Вроньский (1776 – 1853) – польский математик и философ - мистик)

Теорема. Если функции y1 , y2 ,..., yn линейно зависимы, то составленный для них определитель Вронского равен нулю.

Теорема. Если функции y1 , y2 ,..., yn линейно независимы, то составленный для них определитель Вронского не равен нулю ни в одной точке рассматриваемого интервала.

Теорема. Для того, чтобы система решений линейного однородного дифференциального уравнения y1 , y2 ,..., yn была фундаментальной необходимо и

достаточно, чтобы составленный для них определитель Вронского был не равен нулю.

Теорема. Если y1 , y2 ,..., yn - фундаментальная система решений на интервале (a,

b), то общее решение линейного однородного дифференциального уравнения является линейной комбинацией этих решений.

y = C1 y1 +C2 y2 +... +Cn yn ,

где Ci –постоянные коэффициенты.

Применение приведенных выше свойств и теорем рассмотрим на примере линейных однородных дифференциальных уравнений второго порядка.

Общее решение линейного однородного дифференциального уравнения второго порядка.

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

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

Тем не менее, если известно одно ненулевое частное решение, то задача может быть решена.

Теорема. Если задано уравнение вида

y′′+ p1 (x)y′+ p2 (x)y = 0 и известно одно

ненулевое решение у = у1, то общее решение может быть найдено по формуле:

1

 

−∫p1

( x)dx

 

y = C2 y1 ∫

 

e

 

 

dx +C1 y1.

y2

 

 

1

 

 

 

 

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

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

Решение дифференциального уравнения вида y(n) + a1 y(n−1) +... + an y = 0 или, короче, L(y) = 0 будем искать в виде y = ekx , где k = const.

Т.к. y′ = kekx ;

y′′ = k2ekx ; ... y(n) = kn ekx , то

 

L(ekx ) = ekx (kn + a1kn−1 +... + an ).

При этом многочлен F (k) = kn + a1kn−1 +... + an называется характеристическим многочленом дифференциального уравнения.

Для того, чтобы функция y = ekx являлась решением исходного дифференциального уравнения, необходимо и достаточно, чтобы

L(ekx ) = 0; т.е. ekx F (k) = 0.

Т.к. ekx ≠ 0, то F (k) = 0 - это уравнение называется характеристическим уравнением.

Как и любое алгебраическое уравнение степени n, характеристическое уравнение kn + a1kn−1 +... + an = 0 имеет n корней. Каждому корню характеристического уравнения ki соответствует решение дифференциального уравнения.

Взависимостиоткоэффициентовk характеристическоеуравнениеможетиметьлибо n различных действительных корней, либо среди действительных корней могут быть кратные корни, могут быть комплексно – сопряженные корни, как различные, так и кратные.

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

1)Составляем характеристическое уравнение и находим его корни.

2)Находим частные решения дифференциального уравнения, причем:

a)каждому действительному корню соответствует решение ekx;

б) каждому действительному корню кратности m ставится в соответствие m решений:

ekx ; xekx ; ...

xm−1ekx .

в) каждой паре комплексно –

сопряженных корней α ±iβ

характеристического уравнение ставится в соответствие два решения: eαx cosβx и eαx sin βx.

г) каждой паре m – кратных комплексно – сопряженных корней α ±iβ характеристического уравнения ставится в соответствие 2m решений:

eαx cosβx,

xeαx cosβx,

... xm−1eαx cosβx,

eαx sin βx,

xeαx sin βx,

...xm−1eαx sin βx.

3) Составляем линейную комбинацию найденных решений.

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

Пример. Решить уравнение y′′′− y = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Составим характеристическое уравнение: k3 −1 = 0;

 

 

 

 

 

 

 

 

 

 

 

 

(k −1)(k2 + k +1) = 0;

k

=1;

 

k2

+ k +1 = 0;

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= − 1

 

 

 

 

 

 

 

 

 

 

= − 1

 

−

 

 

 

i;

D =1− 4 = −3;

k

2

+

 

 

 

3

 

i;

k

3

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

3

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

C

 

 

 

 

 

 

 

 

sin

 

 

 

x .

 

Общее решение имеет вид:

y = C ex + e

−

cos

 

 

 

3

 

x +C

 

3

 

2

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример. Решить уравнение (1− x2 )y′′− 2xy′+ 2y = 0.

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

Таким частным решением будет являться функция y1 = x.

 

 

 

 

 

 

 

y1′ =1; y1′′ = 0;

 

 

 

 

 

 

0 − 2x + 2x = 0;

 

 

Исходное дифференциальное уравнение можно преобразовать:

 

 

 

 

 

 

 

 

 

 

 

y′′−

 

2x

y′

+

 

2y

 

 

 

= 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1− x2

 

1− x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 ∫

2x

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение имеет вид:

 

 

 

 

 

 

∫x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = C x

 

 

 

 

 

e 1−x

 

 

 

dx +C

2

x;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∫

e−ln(1−x2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = C x

 

 

x2

 

 

dx +C

2

x;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = C x

 

dx

 

+C

 

x;

 

 

 

 

 

 

y = C

 

 

x

+C x

 

 

 

 

 

1

+

1

+

1

dx;

∫x2 (1− x2 )

2

 

 

 

 

 

 

2

∫

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 2(1− x)

2(1+ x)

 

 

 

 

 

y = C

2

x +C x −

 

1 + 1 ln

 

1+ x

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

x 2

 

 

 

 

1

− x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончательно:

y = C2 x +C3 xln

 

1+ x

 

 

+C4 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1− x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример. Решить уравнение yIV − y = 0.

 

 

 

 

 

 

 

 

 

 

 

Составим характеристическое уравнение: k4

−1 = 0.

 

 

 

 

 

 

 

 

 

 

(k2 −1)(k2

+1) = 0;

 

 

 

k =1; k

2

= −1;

k

3

= i;

k

4

= −i.

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

y = C ex +C

2

e−x +C

3

cos x +C

4

sin x.

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример. Решить уравнение y′′− 4y′+ 4y = 0.

 

 

 

 

 

 

 

 

 

Характеристическое уравнение: k2

− 4k + 4 = 0;

 

k

 

= k

2

= 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

y = C e2x +C

2

xe2x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

′′

 

 

′

 

5y =

0.

 

 

 

 

 

 

 

 

 

y

 

 

+ 2y +

 

 

 

 

 

 

 

 

 

Характеристическое уравнение: k2

+ 2k +5 = 0;

 

D

= −16;

 

 

k

= −1+ 2i;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k2 = −1− 2i.

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

y = e−x (C

cos 2x +C

2

sin 2x).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример. Решить уравнение y′′′−7 y′′+ 6y′ = 0.

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