- •Лекция 2. Дифференциальные уравнения высшего порядка. Задача Коши. Теорема существования и единственности решения задачи Коши. Общее решение и общий интеграл. Методы понижения порядка дифференциального уравнения
- •Лекция 4. Неоднородные линейные дифференциальные уравнения. Структура общего решения. Принцип суперпозиции решения. Метод Лагранжа вариации произвольных постоянных
- •Лекция 5. Линейные дифференциальные уравнения с постоянными коэффициентами. Метод Эйлера решения однородных уравнений. Метод подбора неоднородных уравнений
- •Лекция 6. Комплексные числа, последовательности комплексных чисел
- •Лекция 7. Функции комплексной переменной. Предел. Непрерывность. Дифференцирование функции комплексной переменной. Понятие аналитической функции комплексной переменной
- •Лекция 8. Ряды с комплексными числами. Степенные ряды. Элементарные функции комплексной переменной
- •Лекция 9. Интеграл от функции комплексной переменной по кривой на комплексной плоскости. Интеграл от аналитической функции. Формула Ньютона-Лейбница
- •Лекция 10. Интегральная формула Коши. Интегральная теорема Коши. Теоремы Морера и Лиувилля. Разложимость аналитической функции в степенной ряд. Ряд Тейлора
- •Лекция 11. Ряд Лорана. Теорема Лорана. Изолированные особые точки однозначной аналитической функции. Вычет в изолированной особой точке
- •Лекция 12. Вычисление контурных интегралов с помощью вычетов
- •Лекция 13. Основные понятия операционного исчисления. Преобразование Лапласа, его свойства
- •Лекция 14. Обращение преобразования Лапласа (формула Меллина). Восстановление оригиналов по известным изображениям. Решение дифференциальных уравнений операционным методом
- •Лекция 16. Автономная нормальная система двух линейных дифференциальных уравнений. Фазовые траектории, фазовый портрет. Понятие об устойчивости точки покоя системы
Лекция 4. Неоднородные линейные дифференциальные уравнения. Структура общего решения. Принцип суперпозиции решения. Метод Лагранжа вариации произвольных постоянных
Неоднородные линейные дифференциальные уравнения
Линейным дифференциальным уравнением n–го порядка называется уравнение
pn ( x) y(n ) + pn-1 (x) y( n-1) +... + p1 (x) y¢ + p0 ( x) y = f (x), |
(4.1) |
в котором неизвестная функция y = y(x) и все ее производные входят линейным образом
(т.е. с целой |
неотрицательной |
степенью |
не |
выше |
).первойПри |
этом |
функции |
|
pn (x),..., p0 (x) |
называются коэффициентами |
уравнения (4.1), |
а правая |
часть f (x) – |
||||
неоднородностью этого уравнения. Если в (4.1) |
f (x) º/ 0, то уравнение (4.1) называется |
|||||||
неоднородным дифференциальным уравнением. |
|
|
|
|
|
|||
Структура общего решения. |
|
|
|
|
|
|
|
|
Теорема 4.1 (о структуре |
общего |
решения |
неоднородного уравнения). Если в |
|||||
уравнении (4.11) все коэффициенты ai (x) и правая часть h(x) |
непрерывны на отрезке |
|||||||
[a;b] , то общее решение уравнения (4.11) (на этом отрезке) имеет вид |
|
|
||||||
y = y(x;=C1,..., Cn ) C1 y1 ( x) +... + Cn yn (x) + y* ( x) Û yо.н. = yо.о. + yч.н. , |
(4.2) |
где y1 (x),..., yn (x) – фундаментальная система решений соответствующего однородного уравнения Ly = 0, а y = y* (x) – частное решение неоднородного уравнения (4.11), C1,..., Cn
– произвольные постоянные. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Доказательство. Применяя оператор L к функции (4.2), будем иметь |
|
|||||||||||||
|
Ly = C1Ly1 (x) +... + Cn Lyn (x) + Ly* (x) º 0 +... + 0 + h(x) º h(x). |
|
||||||||||||
Это означает, |
что функция (4.2) |
является решением уравнения(4.1) при произвольных |
||||||||||||
значениях постоянных C ,..., C |
n |
. |
Пусть |
теперь (x , y |
, y0 ,..., y0 |
) |
– произвольная точка в |
|||||||
|
1 |
|
|
|
0 0 |
|
1 |
n -1 |
|
|
|
|
|
|
Rn +1 ( x Î[a;b]). Покажем, что решение y = y(x) |
задачи Коши |
|
|
|
|
|
|
|||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ly = h( x), y(x0 ) |
|
|
0 |
(n |
-1) |
(x0 ) |
0 |
|
= |
|
|
(4.3) |
|
|
|
y0 ,=y¢(x0 ) y1 |
,..., y = |
yn-1 |
|
|
||||||||
можно получить из (4.2) выбором определенных значений C = C 0 |
,..., C |
n |
= C 0 |
постоянных. |
||||||||||
Подчиняя (4.2) условиям (4.3), будем иметь |
|
|
|
|
1 |
|
n |
|
||||||
|
|
|
|
|
|
|
|
|
||||||
C1 y1 (x0 ) +¼+ Cn y=n (x0 ) |
|
y0 - y* (x0 ), |
|
|
|
|
|
|
|
|
|
|||
¢ |
¢ |
|
|
0 |
¢ |
|
|
|
|
|
|
|
|
|
C1 y1 |
(x0 ) +¼+ Cn y=n (x0 ) |
y1 - y* (x0 ), |
|
|
|
|
|
|
|
|
|
. . .
C y(n-1) |
(x |
) +¼+ C |
y(n-1) |
(x=) |
y0 |
- y(n-1) |
(x ) Û |
1 1 |
0 |
n |
n |
0 |
n-1 |
* |
0 |
æ |
y1 (x0 ) |
ç |
y1¢(x0 ) |
Û ç |
¼ |
ç |
çè y1(n-1) (x0 )
¼ yn (x0 ) |
|
öæ C1 |
ö |
æ y0 - y* (x0 ) |
|||||||
|
yn¢ (x0 ) |
|
֍ |
|
÷ |
ç |
|
0 |
- y*¢(x0 ) |
||
¼ |
|
÷çC2 ÷ |
= ç |
|
y1 |
||||||
¼ |
¼ |
|
֍ M |
÷ |
ç |
|
|
M |
|
||
¼ |
y(n-1) |
(x |
) |
֍C |
n |
÷ |
ç y0 |
|
- y( n-1) |
(x |
|
|
n |
0 |
|
øè |
ø |
è |
n-1 |
* |
0 |
ö
÷
÷÷.
÷
) ø
Принцип суперпозиции решения. |
|
Если правая частьh(x) |
уравнения (4.1) является линейной |
непрерывных на отрезке [a;b] функций h1 (x),..., hm (x), т.е.
(4.4)
комбинацией
24
то линейная комбинация |
h(x) º a1h1(x) +... +amhm (x), |
|
|
|
|
|
|
|||||||||
y º a1 y1 (x) +... +am ym (x), |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||
решений y1(x),..., ym (x) |
уравнений Ly1 = h1 (x),..., Lym = hm (x) |
|
является |
|
решением |
|||||||||||
уравнения (4.1) (здесь числа a1 ,...,am могут быть и комплексными). |
|
|
|
|
|
|||||||||||
Действительно, в силу линейности оператора L из тождеств |
|
|
|
|
|
|||||||||||
вытекает тождество |
Ly1 º h1 (x), ... , Lym º hm (x) |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
L[a1 y1(x) +... +am ym (x)] º a1Ly1 (x) +... +am Lym (x) º a1h1 (x) +... +am hm (x). |
|
|
|||||||||||||
Пример 4.1. Найти общее решение уравнения |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
Ly º y¢¢+ y |
|
=x2 +1. |
|
|
|
|
|
(4.5) |
|||
Так |
как в |
правой |
части(4.5) нет ea x , |
то a = 0, и |
поскольку |
в правой |
части |
нет |
||||||||
синусов |
и |
косинусов, то |
b = 0. |
Значит, |
спектральное |
значение |
|
правой |
части |
равно |
||||||
a + ib = 0. |
Характеристическое уравнение l2 +1 = 0 |
имеет два различных чисто мнимых |
||||||||||||||
корня l1,2 = ±i и |
спектральное значениеa + ib = 0 |
отличается от |
них, поэтому |
частное |
||||||||||||
решение неоднородного уравнения (4.5) следует искать в виде |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
y = y |
|
= |
Ax2 + Bx + C |
|
|
|
|
(4.6) |
||
|
|
|
|
|
|
|
ч.н. |
|
|
|
|
|
|
|
|
|
(т.е. в таком же виде, что и правая часть уравнения(4.5)). |
Для |
вычисления |
||||||||||||||
неопределенных |
коэффициентов A,B,C |
|
надо подставить(4.6) в |
(4.5) и |
произвести |
|||||||||||
приравнивание |
коэффициентов |
при |
|
|
одинаковых |
степеняхx. Вычислим |
сначала |
|||||||||
производные от многочлена (4.6): |
y¢ = 2 Ax + B, y¢¢ = 2 A. |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Значит, |
Lyч.н. º 2 A + Ax2 + Bx + C =x2 +1. Приравнивание коэффициентов дает: |
|
|
|
||||||||||||
|
|
|
x2 |
|
A =1, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
x1 |
|
B = 0, Û =A= =1, B 0, C 1- 2=A -1. |
|
|
|
|
|
||||||
|
|
|
x0 |
|
2 A + C =1 |
|
|
|
|
|
|
|
|
|
|
Подставляя найденные коэффициенты в(4.6), найдем окончательно частное решение
неоднородного |
уравнения (4.5): |
y = y |
= |
x2 -1. |
Поскольку |
|
корниl |
= ±i |
|||||||||||||
|
|
|
|
|
|
|
ч.н. |
|
|
|
|
|
|
|
|
|
|
|
|
1,2 |
|
характеристического |
уравнения (*) |
различны, |
то |
общее |
|
решение |
соответствующего |
||||||||||||||
однородного |
уравнения |
имеет видyo.o. = C1 cos x + C2 sin x, |
|
а |
|
значит общее |
|
решение |
|||||||||||||
исходного уравнения (4.5) запишется в форме y = y |
= |
C cos x + C |
2 |
sin x + x2 -1. |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
o.н. |
1 |
|
|
|
|
|
|
|
|
|
||
Пример 4.2. Вычислить общее решение уравнения |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
Ly º y¢¢ + y¢= |
|
x + 2. |
|
|
|
|
|
|
|
|
|
(4.7) |
|||
Найти |
общее решение |
соответствующего |
однородного |
|
|
|
¢¢ |
+ y |
¢ |
= 0 не |
|||||||||||
|
уравненияy |
|
|||||||||||||||||||
составляет |
труда, так |
как |
здесь корни |
характеристического |
|
|
2 |
|
|
||||||||||||
уравненияl + l = 0 |
|||||||||||||||||||||
различны: |
l |
= 0,= l |
-1. |
Оно |
имеет |
видy |
|
|
= C e0 x + C |
e=-1x |
|
C + C |
e- x . |
Займемся |
|||||||
|
1 |
2 |
|
|
|
|
o.o. |
|
|
1 |
2 |
|
|
1 |
2 |
|
|
|
|
|
вычислением частного решения неоднородного уравнения (4.7). Поскольку в правой части нет экспоненты ea x и синусов и косинусов, то спектральное значение l* = a + ib = 0. Оно
является корнем характеристического уравнения кратностиr=1. Частное решение неоднородного уравнения (4.7) следует искать в виде
y = yч.н. = x( Ax + B) º Ax2 + Bx. |
(4.8) |
25
Вычисляя производные y¢ = 2Ax + B, y¢¢ = 2 A функции (4.8) и подставляя ее в (4.7), будем иметь
2A + 2Ax + B= x + 2.
|
|
Приравнивая здесь коэффициенты при одинаковых степенях x, получаем |
|
||||||||
|
|
|
x1 |
|
2A =1, |
Û A =1/ 2,=B 2 - 2 A =1. |
|
||||
|
|
|
|||||||||
|
|
|
x0 |
|
2A + B = |
2 |
|
|
|
|
|
|
|
Подставляя A и B |
|
в (4.8), |
найдем частное |
решение уравнения(4.7) в |
виде |
||||
yч.н. |
= |
1 |
x2 + x, а значит общее решение этого уравнения запишется в форме |
|
|||||||
|
|
||||||||||
|
2 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
y = C + C e-x + |
x2 + x. |
|
|
|||
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
1 |
2 |
2 |
|
|
|
|
|
Пример 4.3. Решить уравнение |
|
|
|
|
|||||
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
Ly º y¢¢ -3y¢ + 2 y =(x2 + x)e3x . |
(4.9) |
||||
|
|
Спектральное значение правой |
частиl* = a + ib |
равно l* = 3, так как a = 3, а |
синусы и косинусы отсутствуют. Характеристическое уравнение l2 - 3l + 2 = 0 имеет два различных действительных корня: l1 =1 и l2 = 2. Так как спектральное значение l* = 3 не
является корнем характеристического уравнения, то частное решение уравнения(4.9) следует искать в виде:
|
|
|
|
|
|
|
|
|
|
y = y |
|
= e3x ( Ax2 + Bx + C ). |
|
(4.10) |
|
||||
|
|
|
|
|
|
|
|
|
|
|
ч.н. |
|
|
|
|
|
|
|
|
Для составления уравнений относительно неопределенных коэффициентовA, B и |
|
||||||||||||||||||
C, найдем производные функции (4.10): |
|
|
|
|
|
|
|
|
|
||||||||||
|
y¢ = e3x (3Ax2 + 3Bx + 3C) + e3 x (2 Ax + B) e3 x (3=Ax2 + 3Bx + 2 Ax + 3C + B), |
|
|
||||||||||||||||
|
y¢¢ = e3 x (9 Ax2 + 9Bx + 6 Ax + 9C + 3B) + e3 x (6 Ax + 3B + 2 A) = |
|
|
|
|||||||||||||||
|
= e3 x (9 Ax2 + 9Bx +12 Ax + 9C + 6B + 2A). |
|
|
|
|
|
|
|
|||||||||||
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Lyч.н. º e3x (9Ax2 + 9Bx +12Ax + 9C + 6B + 2 A) - 3e3 x (3Ax2 + 3Bx + 2Ax + 3C + B) + |
|
|
|||||||||||||||||
+2e3x ( Ax2 + Bx + C) º (x2 + x)e3x . |
|
|
|
|
|
|
|
|
|
||||||||||
Сокращая здесь экспоненту, получаем тождество |
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
2Ax2 + 2Bx + 6Ax + 2C + 3B + 2 A º x2 + x. |
|
|
|
|||||||||||
Приравнивание коэффициентов дает: |
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
x2 |
|
|
2 A =1, |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
x1 |
|
|
2B + 6 A =1, |
Û A =1/ 2, B= -1, C =1. |
|
|
|
||||||||
|
|
|
|
x0 |
|
2C + 3B + 2 A = 0 |
|
|
|
|
|
|
|
|
|
||||
Подставляя |
A, B |
и |
C |
в (4.10), |
получаем |
частное |
решение |
в |
виде |
||||||||||
y = y= |
e3 x (x2 / 2 - x +1). |
При |
|
этом |
общее |
решение |
соответствующего |
однородного |
|||||||||||
ч.н. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
yравнения строится по корням l1 =1 и l2 |
= 2 характеристического уравнения следующим |
|
|||||||||||||||||
образом: |
y |
o.o. |
= C ex |
+ C |
e2 x , |
а |
значит |
общее |
решение |
неоднoродного уравнения(4.9) |
|
||||||||
|
|
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
запишется в виде |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|||
|
|
|
|
|
|
y = y= |
|
C ex + C |
e2 x + e3 x æ |
x2 - x |
+1ö. |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
o.н. |
|
1 |
2 |
|
ç |
|
|
÷ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
è 2 |
|
ø |
|
|
|
Метод Лагранжа вариации произвольных постоянных.
Поиск общего решения неоднородного дифференциального уравнения(4.1) сводится к двум процедурам:
26
1) |
построение |
фундаментальной |
системы |
решенийy (x),..., y |
n |
(x) |
соответствующего однородного уравнения; |
|
1 |
|
|||
|
|
|
|
|||
2) вычисление частного решения y = y* (x) неоднородного уравнения (4.1). |
|
|
||||
Самым трудным |
является осуществление |
первой |
процедуры. Однако |
|
для |
уравнений с постоянными коэффициентами(см. следующий раздел) ее можно всегда
реализовать. Если же найдена фундаментальная система решенийy (x),..., y |
n |
(x) |
1 |
|
однородного уравнения Ly = 0, то реализовать вторую процедуру не составляет особого труда.
Теорема 4.2. Пусть y1 (x),..., yn (x) – фундаментальная система решений однородного уравнения (4.5) с непрерывными на отрезке[a;b] коэффициентами ai (x). Если правая часть h(x) соответствующего неоднородного уравнения (4.1) непрерывна на отрезке [a;b] , то его частное решение можно вычислить в виде
y = y *=(x) C1 (x) y1 (x) +... + Cn (x) yn (x), |
(4.11) |
где функции C1(x),..., Cn (x) (представляющие собой варьированные постоянные общего решения однородного уравнения (4.5)) находятся из системы
æ y1 (x) |
¼ yn (x) |
||
ç |
y1¢(x) |
¼ |
yn¢ (x) |
ç |
|||
ç |
¼ |
¼ |
¼ |
ç |
|
¼ yn(n-1) (x) |
|
è y1(n-1) (x) |
öæ C¢ ö |
æ |
0 |
ö |
|
||
֍ |
1 |
÷ |
ç |
|
÷ |
|
÷çC2¢ |
÷ |
= ç |
0 |
÷. |
(4.12) |
|
֍ |
M |
÷ |
ç |
M |
÷ |
|
֍ |
|
÷ |
ç |
|
÷ |
|
øèCn¢ |
ø |
è h(x) ø |
|
Пример |
4.4. |
Проверить, |
|
что |
функции |
|
y1 = cos x, y2 |
= sin x |
|
образуют |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
¢¢ |
+ y = 0 и |
найти |
общее |
решение |
|||||||
фундаментальную систему решения уравненияy |
|||||||||||||||||||||||
неоднородного уравнения y¢¢ + y = 1 / cos x. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Поскольку |
y¢(x) = -sin x и |
y ¢¢ |
= -cos x, |
то |
функцияy = cos x |
удовлетворяет |
|||||||||||||||||
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
уравнению |
y¢¢ + y = 0 . |
Точно |
так |
же |
убеждаемся, |
что |
функция |
y2 (x) = sin x |
также |
||||||||||||||
удовлетворяет уравнению y¢¢ + y = 0 . Вычисляем вронскиан |
|
|
|
|
|
|
|||||||||||||||||
|
|
W[ y1 |
, y2 ] = |
|
y1 |
y2 |
|
|
|
cos x |
sin x |
|
|
cos |
2 |
x + sin |
2 |
x =1. |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
y1¢ |
= |
|
|
|
= |
cos x |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
y2¢ |
|
|
-sin x |
|
|
|
|
|
|
|
|
|
|
|
||||
Видим, что |
он |
не |
обращается |
|
в нуль |
на |
промежутке(-¥; +¥) , значит |
функции |
|||||||||||||||
y1 = cos x, y2 |
= sin x образуют фундаментальную систему решений уравнения y¢¢ + y = 0 . |
||||||||||||||||||||||
Найдем |
теперь |
частное |
|
решениеy = y* (x) |
неоднородного |
|
уравнения |
||||||||||||||||
y¢¢ + y¢ = 1 / cos x |
в форме y = C1 (x) cos x + C2 (x) sin x. |
При |
этом |
|
функцииC1 (x) |
и C2 (x) |
|||||||||||||||||
должны удовлетворять системе |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
æ C¢ ö |
æ |
cos x |
|
ç 1 |
÷ |
= ç |
-sin x |
èC2¢ |
ø |
è |
C1¢ = - sin x cos x
C2¢ = ò1=dx
sin x ö-1 æ |
0 |
|
÷ |
ç |
|
cos x ø |
è1 / cos |
,C2¢ =1 Û C1
x+ C2 .
ö |
Û |
æ C¢ |
ö |
÷ |
ç 1 |
÷ |
|
x ø |
|
èC2¢ |
ø |
-ò |
sin xdx |
|
|
= |
|
|
|
|
cos x |
|
æcos x |
-sin x öæ |
= |
0 ö |
Û |
|
ç |
cos x |
֍ |
÷ |
||
è sin x |
øè1 / cos x ø |
|
= |
d cos x |
|
|
|
|
=ln | cos x | +C |
, |
||||
|
|||||
ò |
cos x |
1 |
|||
|
|
|
Поскольку нас интересует частное |
решение неоднородного уравнения |
|||||
C (x) |
и C |
2 |
(x) можно |
взять |
в видеC (x) = ln | cos x |,C (x) = x. |
|
1 |
|
|
|
1 |
2 |
|
y = C1 (x) cos x + C2 (x) sin x, |
получаем частное решение в виде |
|
y¢¢ + y¢ = 1 / cos x , то Подставляя их в
27
y = y* (x) (ln | cos x=|) cos x + x sin x,
а значит общее решение неоднородного уравнения запишется в форме
y = C1 cos x + C2 sin x + (ln | cos x |) cos x + x sin=x (C1 + ln | cos x |) cos x + (C2 + x) sin x.
28