№179.11.Ястребов
.pdf31
Наложим на функции C1, C2 условие: |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
C ′ y + C |
′ y |
2 |
= 0; |
|
|
|
(23) |
|||||
|
|
|
|
1 |
1 |
2 |
|
|
|
|
|
|
|
|
|
тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
′ |
|
|
|
′ . |
|
|
|
|
|
|
|
|
yчн |
= C1 y1 |
+ C2 y2 |
|
|
|
|
|
|
|||||
Дифференцируем повторно: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
′′ |
′ ′ |
|
′′ |
|
|
′ ′ |
|
′′ . |
|
|
|
|||
|
yчн = C1 y1 |
+ C1 y1 + C2 y2 |
+ C2 y2 |
|
|
|
|
||||||||
Подставим |
выражения |
для |
yчн, yчн, yчн |
в исходное |
уравнение и |
||||||||||
|
|
′ |
|
|
′′ |
||||||||||
сгруппируем по отдельности слагаемые с C1 и с C2 : |
|
|
|
|
|||||||||||
′′ |
′ |
|
|
′′ |
′ |
|
|
|
|
|
′ ′ |
′ |
′ |
f ( x) |
. |
C1( y1 + py1 + qy1) + C2 ( y2 + py2 |
+ qy2 ) + C1 y1 |
+ C2 y2 = |
|
||||||||||||
Множители при C1 и C2 равны тождественно нулю, поскольку функции |
|||||||||||||||
y1 и y2 являются решениями однородного уравнения, так что |
|
|
|||||||||||||
|
|
|
C1 y1 |
+ C2 y2 = f (x) |
. |
|
|
(24) |
|||||||
|
|
|
|
′ ′ |
′ |
|
′ |
|
|
|
|
|
Соотношения (23) и (24) дают систему двух уравнений относительно неиз-
вестных C′, C′ :
1 2
|
|
|
|
|
+ C2′ y2 |
= 0, |
||||
C1′ y1 |
||||||||||
|
|
′ |
|
′ |
|
|
′ |
|
′ |
(25) |
C |
y |
+ C |
|
y |
= f (x). |
|||||
|
|
2 |
|
2 |
||||||
|
1 |
1 |
|
|
|
|
Решая эту систему, получим выражения для C′, C′ через уже известные
1 2
функции |
y1, y2 |
, y1 |
, y2 |
C1(x), C2 |
(x) |
на- |
|
′ |
′ , после чего сами функции |
|
|
ходятся интегрированием.
В качестве итога сформулируем алгоритм решения неоднородного уравнения методом вариации произвольных постоянных:
1. Решение соответствующего однородного уравнения, получение функций y1, y2 .
32
2. Запись общего решения соответствующего однородного уравнения в виде yоо = C1 y1(x) + C2 y2 ( x) .
3.Вычисление производных y′, y′ .
12
4.Запись системы (25) для отыскания C′, C′ .
12
5.Решение системы, получение функций C′, C′ .
12
6.Нахождение каких-либо первообразных C1(x), C2 (x) интегри-
рованием функций |
C1 |
, C2 |
|
|
′ |
′ . |
|
7. Запись частного |
|
решения неоднородного уравнения в виде |
yчн = C1(x) y1(x) + C2 (x) y2 (x) .
8. Запись общего решения неоднородного уравнения: yон = yоо + yчн.
Пример. Решим методом вариации произвольных постоянных урав-
нение y′′ + 2 y′ = x. Следуя алгоритму, последовательно получаем:
1. y′′ + 2 y′ = 0; z2 + 2z = 0; z1 = −2, z2 = 0 ;
y1 = e−2 x ; y2 = 1.
2. |
y = C e−2 x |
+ C |
2 |
. |
|
|
|
||||||||
|
оо |
|
1 |
|
|
|
|
|
|
|
|
|
|
||
3. |
y1 = −2e |
−2 x |
; y2 = 0 |
. |
|
|
|||||||||
′ |
|
|
|
|
|
|
|
|
′ |
|
|
|
|
||
|
|
′ |
−2 x |
|
|
|
′ |
|
1 = 0 |
|
|
|
|||
4. |
C1e |
|
|
+ C2 |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
+ 0 = x. |
|
|
|||||
|
C′(−2e−2 x ) |
|
|
||||||||||||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
5. |
C1 |
= − |
1 |
|
xe |
2 x |
; |
C2 |
= |
1 |
x |
. |
|||
′ |
|
2 |
|
|
|
|
|
|
′ |
|
2 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6. Интегрируя по частям, имеем:
|
|
|
|
33 |
|
|
|
|
C (x) = − 1 |
∫ |
xe2 xdx = − |
1 1 xe2 x − |
1 e2 x |
= |
|||
1 |
2 |
|
|
2 |
4 |
|
|
|
|
|
2 |
|
|
= − |
1 |
xe |
2 x |
+ |
1 |
e |
2 x |
= e |
2 x |
− |
1 |
x + |
1 |
|
|
|
|
|
|
|
|
|
|||||||
4 |
|
8 |
|
|
|
|
|
4 |
8 |
. |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
1 |
∫ |
|
|
|
|
|
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
C2 (x) = |
2 |
xdx = |
4 |
x . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
7. yчн = e2 x |
|
− |
1 |
x + |
1 |
|
|
|
|
1 |
x2 |
1 |
= |
1 |
x2 |
− |
1 |
x + |
1 |
|
|||||||||
|
4 |
8 |
e−2 x + |
4 |
4 |
4 |
8 |
. |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
8. yон = yоо + yчн = C1e−2 x + C2 + 14 x2 − 14 x + 18 .
10. Метод неопределенных коэффициентов
Это метод отыскания частного решения линейного неоднородного уравнения с постоянными коэффициентами
y′′ + py′ + qy = f (x)
в случае, когда правая часть f (x) имеет специальный вид, самая общая запись которого:
f (x) = eα x (P (x) cos β x + Q (x)sin β x), |
(26) |
||||||
|
m |
|
n |
|
|
|
|
где Pm (x), Qn (x) — многочлены степени m и n соответственно: |
|
||||||
P (x) = p xm |
+ p xm−1 |
+ ... + p |
|
|
x + p |
|
|
m |
0 |
1 |
m−1 |
m |
|
||
Q (x) = q xn + q xn−1 |
+ ... + q |
|
x + q . |
|
|||
n |
0 |
1 |
n−1 |
|
n |
|
|
Укажем характерные частные случаи функции (26). |
|
|
|
|
|||
1. α = 0; β = 0 . В этом случае правая часть |
f (x) = Pm (x) является |
многочленом (поскольку e0 = 1; cos0 = 1, sin 0 = 0). Например:
34
1) y′′ + py′ + qy = 3x2 − 4
y′′ + py′ + qy = e0 x [(3x2 − 4) cos(0x) + 0sin(0x)]; 2) y′′ + py′ + qy = −3
y′′ + py′ + qy = e0 x [−3cos(0x) + 0sin(0x)].
2. α = 0. В этом случае f (x) = Pm (x) cos β x + Qn ( x)sin β x .
Например:
y′′ + py′ + qy = 4 cos x;
y′′ + py′ + qy = −3cos x + 7x3 sin x.
3. β = 0; Pm (x) = 1. В этом случае f (x) = eα x . Например:
y′′ + py′ + qy = e− x ; y′′ + py′ + qy = eπ x .
4. Qn (x) = 0. В этом случае f (x) = eα x Pm (x) cos β x . В част-
ности, если степень многочлена m = 0, то многочлен является постоянным числом: P0 (x) = p0 ≠ 0
Например:
y′′ + py′ + qy = −8e−9 x cos x;
y′′ + py′ + qy = e2 x (5x + 1) cos12x.
5. Pm (x) = 0 . В этом случае f (x) = eα xQn (x)sin β x . Напри-
мер:
y′′ + py′ + qy = sin 2x;
y′′ + py′ + qy = e2 x (5x + 1)sin12x.
35
Рассмотрим комплексное число r = α + β i , где α и β берутся из записи (26) правой части f (x) неоднородного уравнения. Корни характе-
ристического уравнения z2 + pz + q = 0 записываются в виде
|
z |
|
= − |
p |
± |
|
|
D |
|
|
i , |
|
|
|
|
|
|||||||
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
||||
|
1,2 |
2 |
|
|
2 |
|
|
|
|||
|
|
|
|
|
|
|
|
||||
где D — дискриминант уравнения. |
|
|
|
|
|
||||||
Кратность k числа r |
|
как корня |
характеристического уравнения |
||||||||
может иметь одно из трех значений: |
|
|
|
|
|
||||||
1) k = 0; это означает, |
что |
|
r не является корнем характеристического |
||||||||
уравнения; |
|
|
|
|
|
|
|
|
|
|
|
2) k = 1; это означает, |
что |
r является простым корнем характеристиче- |
|||||||||
ского уравнения; при этом дискриминант уравнения D отличен от нуля; |
3) k = 2 , то есть r — кратный корень; последнее возможно в том и толь-
ко том случае, когда D = 0; корень r при этом вещественный, так как его
мнимая часть β = |
|
|
D |
|
|
= 0. |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|||||||
|
2 |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
На практике для определения кратности k нужно решить (квадрат- |
||||||||||||||
ное) характеристическое уравнение и сравнить его корни с числом r . |
|
|||||||||||||
Пусть l |
— максимальная из степеней m и n многочленов Pm и Qn |
|||||||||||||
в записи правой части (26). |
|
|
|
|
|
|
||||||||
Можно убедиться непосредственной подстановкой в дифференци- |
||||||||||||||
альное уравнение, что функция |
|
|
|
|
|
|
||||||||
y |
чн |
= xk eα x ( A (x) cos β x + B (x)sin β x) , |
(27) |
|||||||||||
|
|
|
|
|
|
|
l |
|
|
|
l |
|
|
|
где Al и Bl — многочлены степени l : |
|
|
|
|||||||||||
A = a xl + a xl −1 |
+ a |
2 |
xl−2 |
+ ... + a |
x + a , |
|
||||||||
|
l |
0 |
1 |
|
|
l−1 |
l |
|
36
Bl = b0 xl + b1xl−1 + b2 xl−2 + ... + bl −1x + bl ,
при некоторых значениях коэффициентов ai , bi этих многочленов явля-
ется частным решением неоднородного уравнения. Для отыскания этих (неопределенных вначале) коэффициентов у функции (27) вычисляются
первая и вторая производные; затем |
yчн, yчн, |
yчн |
|
′ |
′′ |
подставляются в ис- |
ходное уравнение, после чего приравниваются коэффициенты при одинаковых функциях-слагаемых обеих частей равенства. Это дает систему уравнений для отыскания коэффициентов ai , bi .
В итоге отыскание частного решения проводится по следующему ал-
горитму: |
|
|
|
|
1. |
Нахождение корней характеристического уравнения. |
|||
2. |
Определение величин α , β , |
m, n, l . |
||
3. |
Выявление кратности k числа α + β i как корня характеристиче- |
|||
ского уравнения. |
|
|
|
|
4. |
Запись частного решения |
yчн в виде (27) с неопределенными ко- |
||
эффициентами многочленов Al , Bl . |
|
|||
5. |
Вычисление производных |
yчн, yчн |
||
′ |
′′ . |
|||
6. |
Подстановка |
yчн, yчн, yчн |
в дифференциальное уравнение. |
|
′ |
′′ |
7.Приравнивание коэффициентов при одинаковых функциях в обеих частях равенства.
8.Решение получившейся системы уравнений, получение коэффици-
ентов многочленов Al , Bl .
9. Запись частного решения с найденными коэффициентами.
После этого можно записать общее решение в виде
yон = yоо + yчн
37
Примеры. 1. Рассмотрим уравнение y′′ + 2 y′ = 20x + 2 .
Выполняем последовательно инструкции алгоритма:
1. z2 + 2z = 0; z1 = −2; z2 = 0 . Общее решение соответствующего однородного уравнения: yоо = C1e−2 x + C2e0 x = C1e−2 x + C2 .
2. Находим параметры правой части:
20x + 2 = e0 x ((20x + 2) cos(0x) + 0sin(0x)); α = 0; β = 0; m = 1; n = 0; l = 1.
3.α + β i = 0 + 0i = 0; кратность k = 1.
4.yчн = xA1(x) = x(a0 x + a1) = a0 x2 + a1x .
5. |
′ |
′′ |
. |
|
yчн = 2a0 x + a1; yчн = 2a0 |
|
6.Подставляем в уравнение: 2a0 + 2(2a0 x + a1) = 20x + 2 .
7.4a0 x + 2(a0 + a1) = 20x + 2;
4a0 |
= 20 |
; a0 |
= 5; a1 = −4 . |
8. |
|
||
2(a0 + a1) = 2 |
|
|
9. yчн = x(5x − 4) = 5x2 − 4x.
Общее решение неоднородного уравнения:
yоо = C1e−2 x + C2 + 5x2 − 4x .
2. Рассмотрим уравнение
y′′ − 10 y′ + 21y = e6 x (−3x2 + 4x + 2).
Выполняем последовательно инструкции алгоритма:
38
1. z2 − 10z + 21 = 0; z1 = 3; z2 = 7 . Общее решение соответствую-
щего однородного уравнения: yоо = C1e3x + C2e7 x . 2. Находим параметры правой части:
e6x (−3x2 + 4x + 2) = e6x ((−3x2 + 4x + 2)cos(0x) + 0sin(0x));
α = 6; β = 0; m = 2; n = 0; l = 2.
3.α + β i = 6 + 0i = 6; кратность k = 0.
4.yчн = x0e6 x A2 (x) = e6 x (a0 x2 + a1x + a2 );
5.yчн′ = 6e6 x (a0 x2 + a1x + a2 )+ e6 x (2a0 x + a1 ) =
=e6 x (6a0 x2 + (2a0 + 6a1)x + (a1 + 6a2 ));
′′ |
|
6 x |
2 |
. |
|
|
(36a0 x + (24a0 + 36a1)x + (2a0 + 12a1 + 36a2 )) |
||
yчн = e |
|
|
||
6. Подставляем в уравнение: |
|
|||
|
e6 x ( |
36a0 x2 + (24a0 + 36a1)x + (2a0 + 12a1 + 36a2 ))− |
−10e6 x (6a0 x2 + (2a0 + 6a1)x + (a1 + 6a2 ))+ +21e6 x (a0 x2 + a1x + a2 ) = e6 x (−3x2 + 4x + 2).
7. Приводим подобные члены и приравниваем коэффициенты при x2e6 x , при xe6 x и при e6 x :
−3a0 = −3
4a0 − 3a1 = 4
2a0 + 2a1 − 3a2 = 2.
39
8. |
Решаем систему: a0 = 1; |
|
a1 = 0; |
a2 = 0. |
|
|||
9. |
yчн = e6 x (x2 + 0x + 0) = x2e6 x . |
|
||||||
Общее решение неоднородного уравнения: |
|
|||||||
y |
оо |
= C e3x + C e7 x + x2e6 x . |
|
|
||||
|
1 |
2 |
|
|
|
|
|
|
|
|
3. Рассмотрим уравнение y′′ + 4 y = 2cos x − 5sin x . |
||||||
1. |
z2 + 4 = 0; |
z = ± |
−4 = ± |
4i = 0 ± 2i . |
Общее решение соот- |
|||
|
|
|
1,2 |
|
|
|
|
|
ветствующего однородного уравнения: |
|
|||||||
y |
оо |
= e0 x (C cos 2x + C |
2 |
sin 2x) = C cos 2x + C sin 2x . |
||||
|
1 |
|
|
|
1 |
2 |
2. Находим параметры правой части:
2 cos x − 5sin x = e0 x (P0 (x) cos x + Q0 (x)sin x); α = 0; β = 1; m = 0; n = 0; l = 0 .
3.α + β i = 0 + 1i ; кратность k = 0.
4.yчн = x0e0 x (a0 cos x + b0 cos x) = a0 cos x + b0 cos x ;
5.yчн′ = −a0 sin x + b0 cos x ;
′′ |
. |
yчн = −a0 cos x − b0 sin x |
|
6. Подставляем в уравнение: |
|
−a0 cos x − b0 sin x + 4(a0 cos x + b0 sin x) = 2 cos x − 5sin x ; 3a0 cos x + 3b0 sin x = 2 cos x − 5sin x .
7. Приравниваем коэффициенты при cos x и при sin x :
3a0 = 23b0 = −5.
8. a0 = 23 ; b0 = − 53 .
40
9. yчн = 23 cos x − 53 sin x .
Общее решение неоднородного уравнения:
yоо = C1 cos 2x + C2 sin 2x + 23 cos x − 53 sin x .