Дифференциальные уравнения. Учебное пособие-3
.pdf3. ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ 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
