Сандаков Решение ОДУ
.pdfyoo = c1e−x +c2e3x .
Остается найти частное решение неоднородного уравнения (2.6). Сравнивая его правую часть f = 9х2 с (2.3), видим, что п = 2, α = 0, β = 0. Число λ0 = α + iβ = 0 не является корнем характеристического уравнения, следовательно, s = 0. Поэтому уч.н согласно (2.4)
следует искать в виде
учн = ах2 + bx + с.
Подставляя эту функцию в (2.6), получим
2a – 2(2ax + b) – 3(ax2 + bx + c) = 9x2
или
–3aх2 – (4a + 3b)х + 2а –2b – 3c = 9x2.
Приравнивая коэффициенты при одинаковых степенях х, имеем для нахождения коэффициентов а, b, с систему
−3a =9;4a +3b = 0;
2a −2b −3c = 0.
Из этой системы находим: а = –3, b = 4, c = −143 . Следовательно,
yчн = −3x2 + 4x −143 . Тогда общее решение неоднородного уравнения (2.6) имеет вид
y |
он |
= c e−x +c e3x −3x2 |
+ 4x −14 . |
|
|
|
1 |
2 |
3 |
|
|
Пример 2.2. Решить уравнение |
|
||||
|
|
||||
|
|
y′′+ y′−2y = 4e−x . |
(2.7) |
Решение. Сначала найдем общее решение соответствующего однородного уравнения
y′′+ y′−2y = 0.
Его характеристическое уравнение
λ2 + λ – 2 = 0
имеет корни: λ1 = –2 и λ2 = 1. Следовательно, общее решение однородного уравнения получаем в виде yoo = c1e−2 x +c2ex .
11
Остается найти частное решение неоднородного уравнения (2.7). Сравнивая его правую часть f = 4е–х с (2.3), видим, что п = 0, α = –1, β = 0. Число λ0 = α + iβ = –1 не является корнем характеристического уравнения, следовательно, s = 0. Поэтому учн согласно (2.4) следует искать в виде учн = ае–х. Подставляя эту функцию в
(2.7), получим
aе–х – aе–х – 2aе–х = 4е–х,
откуда следует, что а = –2. Следовательно, учн = –2е–х.
Тогда общее решение неоднородного уравнения (2.7) имеет вид
y = c e−2 x +c ex −2e−x . |
|
||
он |
1 |
2 |
|
Пример 2.3. Решить уравнение |
|
|
|
|
y′′+9y = 6cos3x. |
(2.8) |
Решение. Сначала найдем общее решение соответствующего однородного уравнения
y′′+9y = 0.
Его характеристическое уравнение
λ2 + 9 = 0
имеет два комплексно-сопряженных корня: λ1 = 3i и λ2 = –3i. Следовательно, общее решение однородного уравнения получаем в виде
yoo = c1 cos3x +c2 sin3x.
Остается найти частное решение неоднородного уравнения (2.8). Сравнивая его правую часть f = 6cos3х с (2.3), видим, что п = = 0, α = 0, β = 3. Число λ0 = α + iβ = 3i является корнем характеристического уравнения, следовательно, s = 1. Поэтому уч.н согласно (2.4) следует искать в виде
учн = х(acos 3x + bsin3x).
По формуле Лейбница находим
yчн′′ = 2(−3asin3x +3bcos3x)+ x(−9a cos3x −9bsin3x)
и вместе с уч.н , подставляя в (2.8), получим
x(−9a cos3x −9bsin3x)+2(−3asin3x +3bcos3x)+
+9x(a cos3x +bsin3x) = 6cos3x.
После упрощения получаем
−6asin3x +6bcos3x = 6cos3x ,
12
откуда следует, что а = 0, b = 1. Поличили уч.н = хsin3х. Тогда общее решение неоднородного уравнения (2.8) имеет вид
yон = c1 cos3x +c2 sin3x + xsin3x .
Пример 2.4. Решить уравнение |
|
y′′−6y′+10y = 2e3x sinx. |
(2.9) |
Решение. Сначала найдем общее решение соответствующего однородного уравнения
y′′−6y′+10y = 0.
Его характеристическое уравнение
λ2 – 6λ + 10 = 0
имеет два комплексно-сопряженных корня: λ1 = 3 – i и λ2 = 3 + i. Следовательно, общее решение однородного уравнения имеет вид
yoo = c1e3x cosx +c2e3x sinx.
Остается найти частное решение неоднородного уравнения (2.9). Сравнивая его правую часть f = 2е3хcosх с (2.3), видим, что п = = 0, α = 3, β = 1. Число λ0 = α + iβ = 3 + i является корнем характеристического уравнения кратности один, следовательно, s = 1. По-
этому учн согласно (2.4) следует искать в виде
учн = х(A cos x + B sin x) е3х.
Найдем yчн и yчн : |
|
|
|
|
|
|
||
|
′ |
′′ |
|
|
|
|
|
|
′ |
=[(Bx + |
A)cosx +(B − Ax)sinx]e |
3x |
+ |
3(Axcosx + Bxsinx)e |
3x |
||
yчн |
|
|
||||||
′′ |
|
|
|
|
3x |
+6[(Bx + A)cosx + |
|
|
yчн =[(2B − Ax)cosx +(Bx −2A)sinx]e |
|
|
|
+(B − Ax)sinx]e3x +9(Axcosx + Bxsinx)e3x ,
и подставим их вместе с учн в (2.9):
(2B − Ax)cosx −(Bx +2A)sinx +6[(Bx + A)cosx +(B − Ax)sinx]
+9(Ax cosx + Bxsinx)−6[(Bx + A)cosx +(B − Ax)sinx]− −18(Ax cosx + Bxsinx)+10(axcosx + Bxsinx) = sin x
или
2B cosx −2Asin x = 2sin x ,
откуда получим А = –1, В = 0. Следовательно, учн = –хе3хcosх.
Таким образом, окончательно имеем
yон = c1e3x cosx +c2e3x sinx − xe3x cosx.
,
+
13
Пример 2.5. Решить уравнение |
|
y′′+ y′−2y =18xex +3e−2 x . |
(2.10) |
Решение. Сначала найдем общее решение соответствующего однородного уравнения
y′′+ y′−2y = 0.
Его характеристическое уравнение
λ2 + λ – 2 = 0
имеет два корня: λ1 = –2 и λ2 = 1. Следовательно, общее решение однородного уравнения имеет вид
yoo = c1e−2 x +c2ex .
Остается найти частное решение неоднородного уравнения
(2.10). Правая часть уравнения (2.10) f = 18хех + 3е–2х может быть представлена в виде f = f1 + f2, где f1 = 18хех, f2 = 3е–2х.
Тогда для нахождения учн уравнения (2.10) воспользуемся принципом суперпозиции (утверждение 2.2), согласно которому учн = = у1 + у2, где у1 – частное решение уравнения
y′′+ y′−2y =18xex , |
(2.11) |
а у2 – частное решение уравнения |
|
y′′+ y′−2y =3e−2 x . |
(2.12) |
Будем искать у1 в виде у1 = х(ах + b)ех. Подставляя у1 в (2.11), по-
лучим
2a + 2(2ax + b) + (ax2 + bx) + (2ax + b) + (ax2 + bx) – 2(ax2 + bx) = =18x
или |
2а + 3b + 6ах = 18х. |
Следовательно, 2a +3b = 0; откуда а = 3, b = –2.
6a =18,
Итак, у1 = х(3х – 2)ех.
Будем искать у2 в виде у2 = ахе–2х. Подставляя у2 в (2.12), получим
–4a + 4ax + а – 2ax – 2ax = 3,
откуда а = –1. Следовательно, у2 = –хе–2х. Тогда
учн = у1 + у2 = х(3х – 2)ех – хе–2х.
14
Таким образом, общее решение уравнения (2.10) будет иметь
вид
yон = c1e−2 x +c2ex +(3x2 −2x)ex − xe−2 x .
Задачи для самостоятельного решения
1.у″ – 4у′ + 4у = хе2х + 2cos2x.
2.у″ – 4у′ = хsin2x + x2е–4х.
3.у″ + 2у′ + у = 2хе–х + 3sinx.
4.у″ – 9у′ = хе3х + cos3x.
5.у″ + 2у′ + 2у = е–х + cosx + xex.
6.у″ + 6у′ + 9у = 2е–3х + xsinx.
7.у″ – 4у′ + 8у = е2хsin2x + 2xcosx.
8.у″ – 2у′ + у = 3хех + 2cosx.
9.у″ – 4у′ + 3у = 2хех + e3x.
10.у″′ – 3у′′ + 3у′– y = 2хех + 3e–2x.
11.у″ – 6у′ + 9у = 4хе3х + sin2x.
12.у″ – 2у′ + 2у = хехcosx + 2e–x.
13.у″′ + 2у′′ – у′– 2y = хе–2x + chx.
14.у″ – 10у′ + 25у = хе5х + xsin2x.
15.у′′ – 8у′+ 16y = хе4x + ch2x.
16.у′′ + 25y = хsin5x + x2е2x.
17.у′′ – 4у′ + 8y = е2xcos2x + 2xsh2x.
18.у(IV) + 4у″ + 4у = xcos2x + x2chx.
19.у′′ – 9y = xе3x + cos3x.
20.у″ + 8у′ + 16у = 2хе–4х + xsin2x.
21.у″ + 16у = хsin4x + 2x2е2х.
22.у″ + 4у′ + 8у = е–2хsin2x + 3xex.
23.у(IV) + 2у″′ + у″ = 2xe–x + x2.
24.у(IV) + 2у″ + у = xsinx + 2x2e2x.
25.у″′ + 6у′′ + 9у′ = хе–3x + 3ch2x.
15
3. Неоднородные дифференциальные уравнения n-го порядка с постоянными коэффициентами, решаемые методом вариации постоянных
В теме 2 был рассмотрен метод нахождения частного решения неоднородного уравнения (2.1) со специальной правой частью. Метод вариации постоянных, который будет рассмотрен здесь, применим для любого вида линейных неоднородных дифференциальных уравнений независимо от вида правой части и позволяет найти общее решение неоднородного уравнения во всех случаях, если известно общее решение соответствующего однородного уравнения.
Рассмотрим дифференциальное уравнение
a0 (x) y(n) +a1(x) y(n−1) +... +an−1(x) y′+an (x) y = f (x) , |
(3.1) |
где ak (x) и f (x) – заданные непрерывные функции (k = 0, 1, 2, ..., n), определенные на <a; b> , а0(х) ≠ 0. (Всюду ниже через <a, b> будем обозначать любое из множеств (a, b), (a, b], [a, b), [a, b], −∞ ≤ a < < b ≤ +∞.)
Пусть у1(х), у2(х), ..., уп(х) – ФСР соответствующего однородного
уравнения |
|
a0 (x) y(n) +a1(x) y(n−1) +... +an−1(x) y′+an (x) y = 0 , |
(3.2) |
Тогда общее решение однородного уравнения (3.2) имеет вид |
|
уоо = с1у1(х) + с2у2(х) + ... + спуп(х),
где с1, с2, ..., сп – произвольные постоянные.
Утверждение 3.1. Решение неоднородного уравнения (3.1) мо-
жет быть найдено в виде |
|
уон = с1(х) у1 + с2(х) у2 + ... + сп (х) уп, |
(3.3) |
где с1(х), с2(х), ..., сп(х) – неизвестные функции, которые определяются из системы:
c′y |
|
+c′y |
2 |
+... +c′y |
n |
= 0; |
|
||
1 1 |
|
2 |
|
n |
|
|
|||
c1′y1′ +c2′y2′ +... +cn′yn′ = 0; |
|
||||||||
|
|
|
|
|
|
|
|
|
(3.4) |
...; |
|
|
+c2′y2 |
+... +cn′yn |
|||||
c1′y1 |
|
= 0; |
|||||||
|
(n−2) |
|
|
(n−2) |
|
|
(n−2) |
|
|
c1′y1 |
|
+c2′y2 |
+... +cn′yn |
= f (x). |
|||||
|
(n−1) |
|
|
(n−1) |
|
|
(n−1) |
|
16
Замечание. Определитель системы (3.4) всегда отличен от нуля на всём <a, b> так как он совпадает с определителем Вронского от ФСР у1, у2, ..., уп. Тогда из теоремы Крамера следует, что система (3.4) имеет единственное решение ck′(x) = ϕk (x), k =1,2,..., n. От-
куда находим
ck (x) = ∫ϕk (x)dx +ck , k =1,2,..., n,
где ck – произвольные постоянные.
Подставляя найденные сk (х) в (3.3), получим общее решение исходного неоднородного уравнения (3.1).
Рассмотрим примеры применения этого метода. Пример 3.1. Решить уравнение
|
′′ |
|
′ |
|
1 |
|
(3.5) |
|
y |
+3y |
+ 2y =1+ex . |
||||||
|
|
Решение. Сначала найдем общее решение соответствующего однородного уравнения
y′′+3y′+ 2y = 0.
Его характеристическое уравнение
λ2 + 3λ + 2 = 0
имеет два корня: λ1 = –2 и λ2 = –1. Следовательно, общее решение однородного уравнения имеет вид yoo = c1e−2 x +c2e−x , где с1, с2 –
произвольные постоянные.
Общее решение неоднородного уравнения (3.5) будем искать в виде
y |
= c (x)e−2 x +c (x)e−x , |
(3.6) |
|
он |
1 |
2 |
|
где с1(х), с2(х) – неизвестные функции, которые согласно утвержде-
c1′e−2 x +c2′e−x = 0;
нию 3.1 могут быть найдены из системы
−2c1′e−2 x −c2′e−x =1+1ex .
Решая эту систему, находим c1′ = − |
e2 x |
и c2′ |
|
ex |
|
|
= |
|
. Откуда полу- |
||
1+ex |
1+ex |
||||
чаем: |
|
|
|
|
|
17
|
c |
|
(x) = − |
∫ |
e2 x dx |
= −e |
x |
+ln(1+e |
x |
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
)+c |
, |
|
|
||||||||||||||||||||||
|
1 |
|
|
|
|
|
1+ex |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
∫ |
|
ex dx |
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
c (x) = |
|
|
|
|
|
= ln(1 |
+e |
|
)+c . |
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
2 |
|
|
|
|
1+ex |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Подставляя найденные с1(х), с2(х) |
|
в (3.6), получим общее реше- |
||||||||||||||||||||||||||||||||
ние неоднородного уравнения (3.5): |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
y |
|
−2 x |
−e |
−x |
+e |
−2 x |
ln(1+e |
x |
|
|
|
−x |
+e |
−x |
|
ln(1+e |
x |
) |
||||||||||||||||
= c e |
|
|
|
|
|
|
|
|
|
)+c e |
|
|
|
|
||||||||||||||||||||
он |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
|
|
|
−2 x |
|
−x |
+ |
(e |
−2 x |
+e |
x |
)ln(1+e |
x |
), |
|
|
||||||||||||||||
|
|
|
= c e |
|
|
|
+c e |
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
он |
|
1 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где c1 и c2 – произвольные постоянные. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
Пример 3.2. Решить уравнение |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
y′′+ 4y = 4ctg 2x. |
|
|
|
|
|
|
|
|
(3.7) |
||||||||||||||
Решение. Сначала найдем общее решение соответствующего |
||||||||||||||||||||||||||||||||||
однородного уравнения |
|
y′′+ 4y = 0. |
|
|
Его характеристическое урав- |
нение λ2 + 4 = 0 имеет два комплексно-сопряженных корня: λ1 = –2i и λ2 = 2i. Следовательно, общее решение однородного уравнения имеет вид yoo = c1 cos2x +c2 sin2x , где с1, с2 – произвольные посто-
янные.
Общее решение неоднородного уравнения (3.7) будем искать в
виде |
|
yон = c1 (x)cos2x +c2 (x)sin2x , |
(3.8) |
где с1(х), с2(х) – неизвестные функции, которые согласно утверждению 3.1 могут быть найдены из системы
|
|
c′cos2x +c′ sin2x = 0; |
|
||||||||
|
|
|
1 |
2 |
|
|
|
|
|
||
|
|
−2c1′sin2x + 2c2′ cos2x = 4ctg2x. |
|
||||||||
Решая эту систему, находим c1′ = −2cos 2x и c2′ = 2 |
cos2 2x |
. Откуда |
|||||||||
sin2x |
|||||||||||
получаем: |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
c1 (x) = −2∫cos2xdx = −sin2x +c1 , |
|
||||||||
c (x) = 2 |
|
cos2 2x |
dx = 2 |
|
dx |
−2 |
|
sin2xdx =ln|tgx |+cos2x +c . |
|||
∫ |
|
∫sin2x |
∫ |
||||||||
2 |
sin2x |
|
|
|
|
2 |
|||||
|
|
|
|
|
|
18
Подставляя найденные с1(х), с2(х) в (3.8), получим общее решение неоднородного уравнения (3.7):
yон = (−sin2x +c1)cos2x +(ln|tgx |+cos2x +c2)sin2x или yон = c1 cos2x +c2 sin2x +sin2xln|tgx |.
Задачи для самостоятельного решения
1.y′′+ y′ =1+1ex .
2.y′′+ y′−2y =1−1ex .
3.y′′+ y′−2y =sech x .
e−2 x
4. y′′+ 4y′+ 4y = ch 2x . 5. y′′−4y = 4 th 2x .
6. y′′+ 2y′ = sh22 x .
7. y′′−2y′−3y = (1+e2exx )2 .
14.y′′+ 4у′ = 2ctg2x.
15.y′′+ 4y′+ 4y = e−2 x xx −2 .
|
y′′ |
−2y′ = |
|
ex |
|||||||
16. |
|
|
|
. |
|
||||||
1−e2 x |
|||||||||||
|
|
′′ |
|
′ |
|
|
|
|
ex |
|
|
17. |
y |
−3y |
+ |
2y = (1−ex )2 . |
|||||||
|
|
18.y′′+9y = sin32 3x .
19.y′′−2y′+ y =1xe+ xx 2 .
20.y′′+3y′+ 2y = 2 th x .
8. |
y′′− y′−6y |
= |
|
|
ex |
|
y′′ |
+3y′ |
|
|
ex |
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
. |
21. |
= |
|
|
. |
|
|
|
|
|
|
|
||||||||||||||||
|
1+e2 x |
sh2 x |
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
y |
′′ |
−9y =3e |
3x |
cth 3x . |
|
|
′′ |
|
′ |
|
|
|
|
|
|
|
e2 x |
|
|||||||||||||||
9. |
|
22. |
y |
−4y |
+3y = (1+ex )2 . |
|||||||||||||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||||||||||
10. |
|
y′′ |
+ 4y′ = |
e−4 x |
|
y′′ |
−6y′ |
|
|
|
|
|
|
xe2 x |
||||||||||||||||||||
|
|
|
|
. |
|
|
23. |
+9y = |
|
|
|
|
|
. |
||||||||||||||||||||
|
ch2x |
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
x2 +1 |
||||||||||||||||||||||||||||||||
|
|
|
′′ |
|
′ |
|
|
|
|
|
e−x |
|
|
|
′′ |
|
|
|
1 |
|
|
|
|
|
|
|
|
|||||||
11. |
|
y |
+ 2y |
+ y =1+ x2 . |
24. |
y |
+16y = sin2 4x . |
|||||||||||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||||||||||
12. |
|
y′′ |
−3y′ |
= |
|
|
1 |
|
|
|
. |
|
25. |
y′′ |
+ 4y′+3y = e |
−x |
cth x . |
|||||||||||||||||
|
1+e3x |
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
′′ |
|
′ |
|
|
|
|
|
|
|
2e2 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
13. |
|
y |
−4y |
+ 4y = x + 2 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
19
4. Уравнения Эйлера
Уравнения
a0 xn y(n) + a1xn−1 y(n−1) +... + an−2 x2 y′′+ an−1xy′+ an y = f (x) , (4.1)
где а0, а1, а2, ..., ап – постоянные (а0 ≠ 0) называется уравнением Эйлера. Оно сводится к линейному уравнению с постоянными ко-
эффициентами заменой независимого переменного х = еt при х > 0 (или х = –еt при х < 0).
Вычисляя производные
dy |
= |
dy |
|
dt |
= |
1 dy |
, |
d 2 y |
= |
1 d 2 y |
− |
1 |
dy |
, |
||||||||||
dx |
dt |
dx |
x dt |
dx2 |
x2 |
|
dt2 |
|
x2 |
dt |
||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
d |
3 |
|
|
|
|
1 |
|
|
3 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
y3 |
= |
|
|
d |
3y |
−3 d |
|
2y |
+2 dy |
|
и т.д. |
|
|
||||||||
|
|
dx |
|
3 |
|
|
|
|
||||||||||||||||
|
|
|
|
|
x |
dt |
|
dt |
|
|
dt |
|
|
|
|
|
|
и подставляя их в уравнение Эйлера, получим линейное уравнение с постоянными коэффициентами, характеристическое уравнение которого имеет вид
a0λ(λ−1)(λ−2)...(λ−n +1) +... +an−2λ(λ−1) +an−1λ+an = 0 . (4.2)
При составлении этого уравнения каждое произведение xk y(k )
в (4.1) заменяется на произведение к убывающих на единицу чисел
λ(λ – 1) ... (λ – k + 1).
Пример 4.1. Решить уравнение
x3 y′′′+2xy′−2y = 0. |
(4.3) |
Решение. Это однородное уравнение Эйлера. После замены независимого переменного х = еt оно сводится к линейному уравнению с постоянными коэффициентами, характеристическое уравнение которого имеет вид
λ(λ – 1)(λ – 2) + 2λ – 2 = 0 . |
(4.4) |
У него корень λ1 = 1 и два комплексно-сопряженных корня λ2 = = –1 –i, λ3 = –1 + i. При таких λ общее решение однородного уравнения с постоянными коэффициентами (с характеристическим уравнением (4.4)) имеет вид
20