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

Дифференциальные уравнения. Учебное пособие-3

.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
1 Мб
Скачать

3. ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ N-ГО ПОРЯДКА

Дифференциальные уравнения n-го порядка имеют вид

F (x, y, y , , y(n) ) 0 ,

(1)

если они не разрешены относительно старшей производной и

y(n) f (x, y, y ,

, y(n 1) ) ,

(2)

еслиониразрешеныотносительностаршейпроизводной. Сделавзамену

переменных

в уравнении (2): y y , y

y

2

,

, y

y

n 1

, y

 

f (x, y, y ,

1 1

 

 

n 2

 

n 1

 

, yn 1) , получаем систему ДУ.

 

 

 

 

 

 

 

Замечание. Так как ДУ (2) эквивалентно системе линейных дифференциальных уравнений, то все теоремы для СЛНДУ верны для ДУ n-го порядка.

ПРОСТЕЙШИЕ СЛУЧАИ ПОНИЖЕНИЯ ПОРЯДКА ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ

1.  Уравнение не содержит искомой функции и ее производных до порядка (k 1) включительно: F(x, y(k) , y(k 1) , , y(n) ) 0 . Порядок

уравнения можно понизить до (n k) заменой y(k) (x) p(x) . 2.  Уравнение не содержит независимого переменного x:

F ( y, y , y , , y(n) ) 0 . Порядок уравнения можно понизить на единицу заменой y (x) p( y) .

41

3.  Левая часть уравнения (1) есть производная некоторого дифференциального выражения (n – 1)-го порядка.

F(x, y, y , , y(n) ) dxd (x, y, y , , y(n 1) ) .

Тогда первый интеграл этого уравнения имеет вид

(x, y, y , , y(n 1) ) C.

Иногда левую часть уравнения (1) привести к такому виду можно, домножив ее на некоторую функцию (x, y, y , , y(n 1) ) . Введение та-

кого множителя приводит к появлению лишних решений, обращающих множитель в нуль [1]. Если множитель разрывен, то возможна потеря решения.

4.  Уравнение (1) однородно относительно аргументов

x, y, y , ,

y(n) , т. е.

для него выполняется тождество F(x, ky, ky ,

, ky(n) )

k p F (x,

y, y , , y(n) ) [2]. Тогда порядок уравнения (1) может быть

понижен на единицу заменой y(x) exp z(x)dx , где z(x)

– новая не-

известная функция.

 

ЛИНЕЙНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ N-ГО ПОРЯДКА

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

n 1

 

y(n) Pi (x)y(n i) .

(3)

i 1

 

Если коэффициенты Pi непрерывны на [a,

b], то в окрестности лю-

бых начальных значений вида y(x0) y0, y (x0) y1,0, ..., y(n 1) (x0)

yn 1,0 , где x0 (a, b) , выполненыусловиятеоремыосуществовании

иединственности для СНДУ, что гарантирует существование единственного решения ДУ (3). Линейность и однородность уравнения (3)

42

сохраняется при любом преобразовании x (t) , где – произвольная n раз дифференцируемая функция. Причем (t) 0 t. Линейность

и однородность ДУ (3) сохраняется при линейном и однородном преобразовании неизвестной функции y(x) y(x) (x)z(x) .

Введем

линейный дифференциальный

оператор:

L[ y] y(n)

P (x) y(n 1)

P (x) y , тогда (3) можно

записать

так: L[ y] 0 .

1

n

 

 

Определитель Вронского для L[ y] будет иметь вид

 

y1

y2

 

y

y

 

1

2

W (x)

 

 

 

 

y(n 1)

y(n 1)

 

1

2

 

yn

 

 

 

y

 

 

 

n

 

,

 

 

 

y(n 1)

 

 

 

 

 

n

 

 

где yi – линейно независимые решения ДУ (3).

Теорема1. Еслилинейнонезависимыефункции y1, y2 , , yn – это

решение

линейного однородного уравнения (3) с непрерывными

на [a, b]

коэффициентами Pi (x) , то определитель Вронского W (x)

не обращается в нуль ни в одной точке отрезка [a, b] (доказывается аналогично случаю системы линейных дифференциальных уравнений).

Теорема2. Общимрешениемлинейногооднородногоуравнения(3)

с непрерывными на [a, b]

коэффициентами Pi (x) будет линейная ком-

 

n

бинация решений yi , т. е.

y(x) Ci yi (x) , где yi (x) – линейно неза-

 

i 1

висимые на отрезке [a, b] частные решения (3) (доказывается анало-

гично случаю системы линейных дифференциальных уравнений).

Следствие. Максимальное число линейно независимых решений ДУ (3) равно его порядку.

Зная одно нетривиальное частное решение уравнения (3) y1(x) 0 , можносделатьподстановку y y1(x)U (x)dx ипонизитьпорядокурав-

нения, сохранив его линейность и неоднородность. Обычно эту подстановку разбивают на две: y y1(x)z(x), z (x) U (x) . Поскольку такая

43

подстановка – линейна и однородна, она сохраняет линейность и однородность (3), а значит, (3) может быть приведена к виду

a (x)z(n) a (x)z(n 1) a (x)z(x) 0 . Решение y y

в силу под-

0

1

n

1

становки

y y1z

соответствует

решению z 1, следовательно, коэф-

фициент

an 0 .

Сделав замену

z U, получим уравнение с поряд-

ком n 1 .

 

 

 

 

ЛИНЕЙНЫЕ ОДНОРОДНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ N-ГО ПОРЯДКА С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ

Определение. ДУ вида

a

y(n) a y(n 1)

a

n

y 0 ,

(4)

0

1

 

 

 

где ai – постоянныекоэффициенты, называютсялинейнымидифферен-

циальными уравнениями с постоянными коэффициентами. Решение ДУ (4) будем искать в виде y(x) exp(kx) . Подставляя решение в урав-

нение (4), получаем характеристическое уравнение

a kn a kn 1 a

0 .

(5)

0

1

n

 

 

1.  Если корни уравнения (5) различны и действительны, то общее

решение ДУ (4) можно

представить

в

виде

y(x) C1 exp(k1x)

C2 exp(k2 x) Cn exp(kn x) , где Ci – произвольные постоянные; ki – корни уравнения (5); yi (x) exp(ki x) – частные решения ДУ (4).

2.  Если коэффициенты уравнения

(4) – действительные числа,

акорни ki

i и k j i – комплексно-сопряженные, тосогласно

теореме

о

действительном линейном

операторе

[1] решения вида

y1(x) cos x exp x,

y2 (x) sin x exp x тоже будут частными реше-

 

 

 

 

 

 

 

 

 

ниями ДУ (4).

 

 

 

 

 

 

3.  Еслинекоторыйкорень k ki характеристическогоуравнения(5)

имеет кратность

i , то частные решения, соответствующие этому

корню,

определяются

следующим

образом:

y1(x) exp (ki x),

y (x) xexp(k x),

, y

i

(x) x i 1 exp(k x) .

 

2

 

i

 

 

 

i

 

44

Уравнение Эйлера: a xn y(n) a xn 1y(n 1)

a

n

y 0 . Это урав-

0

1

 

 

нение сводится к ДУ (4) заменой x expt, частные решения можно искать в виде y(x) xk .

ЛИНЕЙНЫЕ НЕОДНОРОДНЫЕ ДУ N-ГО ПОРЯДКА. МЕТОД ВАРИАЦИИ ПОСТОЯННЫХ

Пусть дано линейное неоднородное уравнение n-го порядка

 

y(n) p (x) y(n 1) p

n

(x) y f (x) .

(6)

 

1

 

 

Будем искать общее решение уравнения (6) в виде

 

 

n

 

 

 

 

y(x) Ci (x) yi (x) ,

(7)

 

i 1

 

 

 

где Ci (x)

– неизвестные функции; yi (x)

– частные решения однород-

ного ДУ,

соответствующего уравнению (6). Функции Ci (x)

должны

удовлетворять уравнению (6), но так как всего неизвестных функций n, то дополнительные (n – 1) условия, которым будут удовлетворять функции Ci (x) , можно выбрать произвольно [3]. Выберем эти допол-

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

n

n

Например, y (x) Ci (x) yi (x) Ci (x) yi (x) , но при постоянных Ci

i 1

i 1

n

n

y (x) Ci (x) yi (x) и, значит,

Ci (x) yi (x) 0 . Аналогично получаем

i 1

i 1

остальные ограничения:

n

Ci (x) yi (x) 0 ,

i 1 n

Ci (x) yi (x) 0 ,

i 1

45

...

n C (x) y(n 2) (x) 0

i i , i 1

n C (x) y(n 1) (x) f (x)

i i . i 1

Система линейных алгебраических уравнений относительно неизвестных функций Ci (x) имеет единственное решение.

ПРИМЕР РЕШЕНИЯ НЕОДНОРОДНОГО УРАВНЕНИЯ N-ГО ПОРЯДКА

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

1)   y 3y 2 y 2x .

y 3y 2 y 0 ищем

 

 

Решение однородного уравнения

в

виде

y(x) exp kx.

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

k2

3k 2 0

имеет

корни k1 2,

 

 

k2 1 . Тогда общее

решение

однородного

ДУ

y(x) C e2x C

2

ex . Решение неоднородного

ДУ:

y(x) C (x)e2x

1

 

 

 

 

1

 

 

C2 (x)ex . В соответствии с идеей метода вариации постоянной [1] получаем условия для определения неизвестных функций Ci (x), i 1, 2 :

C

(x)e2x C

(x)ex 0,2C

(x)e2x C

(x)ex 2x . Решая эту систему,

1

2

1

2

 

определяем C

(x) e 2x 2x , C

(x)

 

1

 

e 2x 2x

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

ln 2 2

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

C

(x)

 

 

e x 2x C

.

 

 

 

 

 

 

 

 

 

 

 

1

ln 2

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение неоднородного ДУ:

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

x

 

2x

 

 

 

 

 

y(x)

 

 

 

 

 

 

 

 

2

 

C1e

 

 

 

 

 

 

 

 

1

ln 2

 

 

 

 

 

 

 

ln 2 2

 

 

 

 

 

 

 

C1 , C (x) e x 2x ,

2

C2ex .

46

2)  x2 y 3xy 5y 3x2 .

Это уравнение Эйлера [3], сделаем замену x et , тогда ДУ примет вид y 4 y 5y 3e2t . Будем искать решение соответствующего одно-

родного уравнения в виде

y(t) ekt .

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

k2 4k 5 0, его корни k

2 i, k

2

2 i. Согласно теореме о дей-

1

 

 

ствительном операторе[1] частныерешенияоднородногоДУопределя-

ются соотношениями y (t) e2t cost,

y

2

(t) e2t sin t. Тогда общее ре-

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

шение неоднородного ДУ

y(t) C (t)e2t cost C

2

(t)e2t sin t. Получаем

 

 

 

 

1

 

 

 

 

 

 

 

 

Ci (t), i 1, 2 :

условия

для

определения

неизвестных функций

C (t)e2t

cost C (t)e2t sin t 0,

C (t)e2t (2cost sin t) C (t)e2t (2sin t

1

 

2

 

1

 

 

 

 

2

 

 

 

 

 

cost) 3e2t .

Решение

 

этой

системы:

 

 

 

C (t) 3cost C ,

 

 

 

 

 

 

 

 

 

1

1

C2 (t) 3sin t C2 . Общее

решение неоднородного ДУ:

y(t) 3e2t

e2t C1 cost C2 sin t или

y(x) 3x2

e2t C1 cosln

 

x

 

C2 sin ln

 

x

 

.

 

 

 

 

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

1)   y 2 y 3y exp 4x ;

2)   y y sin1 x ;

3)   y y 2exp x x2 ; 4)   y 2y y expx x ; 5)   y 3y 2 y sin x ; 6)   y y 2sec3 x ;

7)   y 5y 4 y 4x2 exp 2x ; 8)   y 4 y 4 y x exp 2x ;

9)   y 3y 4 y exp( 4x) x exp( x) ;

47

10)   y 2 y y 6x exp x ; 11)   y 2 y 3y x2 exp x ;

12)   y 4 y 8 y exp 2x sin 2x ; 13)   y 4 y 4 y x exp 2x ;

14)   y 9 y exp3x cos x ; 15)   y y xsin x ;

16)   y 5y 3x2 sin 5x ; 17)   y y 4xexp x ;

18)   y y 2 y 3x exp x ;

 

 

 

 

 

1

 

 

19)   y

3y

2y exp x 1

;

 

 

20)   y 3y 2 y x cos x ;

 

21)   y y 4shx ;

 

22)   y 2 y y 2xex ex sin 2x ; 23)   y y 4sin x ;

24)  x2 y xy y 8x3 ; 25)  x2 y 3xy 5y 3x2 .

48

ИНТЕГРИРОВАНИЕ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ С ПОМОЩЬЮ РЯДОВ

Теоремаобаналитичностирешения. Пустьданолинейноедиффе-

ренциальное уравнение второго порядка

p0 (x) y p1(x) y p2 (x) y 0 .

(1)

Если pi (x) – аналитические функции аргумента x в окрестности точки x0 , то решения уравнения (1) тоже будут аналитическими функциями в окрестности точки x0 , и, следовательно, их можно представить в виде ряда:

 

 

y(x) an (x x0 )n .

(2)

n 0

Теоремаоразложимостирешенияв обобщенныйстепенной ряд.

Если (1) удовлетворяет условиям теоремы об аналитичности, но при этом x x0 являетсянулемпорядкаs дляфункции p0 (x) , гдеs – конеч-

ное число, нулем порядка s – 1 или выше для функции p1(x) и нулем порядка не ниже s – 2 для функции p2 (x) (s > 2), то существует хотя бы

одно нетривиальное решение (1), которое может быть представлено в виде обобщенного степенного ряда:

 

 

y(x) an (x x0 )n ,

(3)

n 0

где – некоторое действительное число.

49

Замечание. Второе линейно независимое решение (1) тоже имеет вид(3) илиможетбытьпредставленоввидепроизведенияобобщенного степенного ряда и ln(x x0 ) .

Замечание. В конкретных задачах подбирают степенной или обобщенный степенной ряд, формально удовлетворяющий уравнению (1). При подстановке решения в уравнение (1) мы должны получить тождество, котороепозволитопределитькоэффициенты an или an и . Далее

полученный ряд исследуется на сходимость и вычисляется его сумма, которая и будет частным решением (1). Линейная комбинация частных решений в виде бесконечных рядов – общее решение (1).

Перепишем (1) в виде

y y f1(x) yf2 (x) 0 .

(4)

Замечание. С помощью замены переменных можно свести процесс поиска решения в окрестности точки x x0 к поиску решения в окрест-

ности точки x = 0. В дальнейшем будем, не нарушая общности, считать, что решение ищется в точке x = 0. Это решение имеет вид

 

 

y(x) an xn .

(5)

n 0

 

Решение (5) подставим в уравнение (4), приведем подобные и приравняем коэффициенты приодинаковыхстепеняхx. Такимобразом, получим рекуррентные соотношения, позволяющие определить коэффициенты an . Первые одно или два уравнения позволяют вычислить .

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

ТеоремаФукса. Пустьданодифференциальноеуравнение(4) такое, чтофункции f1, f2 имеютполюсывточке x x0 . Тогдарешенияэтого уравнения можно искать в виде обобщенного степенного ряда (3) при

условии, что f (x)(x x )

и f

2

(x)(x x

)2 остаются конечными вели-

1

0

 

 

 

 

0

 

чинами в точке x0 .

 

 

 

 

 

 

 

 

Пример. Пусть дано ДУ (4) и f1,

f2

удовлетворяют условиям тео-

ремы Фукса, т. е. могут быть представлены в виде

 

 

 

1

 

 

 

 

 

f1(x)

 

bn xn bn xn 1,

 

 

 

 

 

x n 0

n 0

50

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]