
- •Теоремы Ферма, Ролля, Лагранжа, Коши
- •Сравнение функций по скорости роста
- •Формулы Маклорена и Тейлора
- •Разложение в ряд Маклорена элементарных функций
- •Вопросы для повторения
- •Условия возрастания и убывания функции
- •Понятие экстремума
- •Необходимое условие экстремума
- •Выпуклость функции. Точки перегиба
- •Схема исследования функции на выпуклость
- •Асимптоты графика функции
- •Исследование функций и построение их графиков
- •Геометрическая интерпретация
- •Свойства эластичности функции
- •Эластичность элементарных функций
- •Вопросы для повторения
- •3 Функции многих переменных
- •Дифференцируемость функции многих переменных
- •Экстремум функции многих переменных
- •Достаточные условия экстремума
- •Метод наименьших квадратов
- •4 Неопределенный интеграл
- •Неопределенный интеграл
- •Основные методы интегрирования
- •5-6 Определенный интеграл
- •Площадь криволинейной трапеции
- •Определение определенного интеграла
- •Площадь плоской фигуры
- •Длина дуги плоской кривой
- •Объем тела вращения
- •Несобственные интегралы
- •Двойные интегралы
- •7-8 Дифференциальные уравнения
- •Задача Коши. Теорема Коши. Понятие общего решения.
- •8(доп) Комплексные числа
- •Краткие сведения теории ЛДУ 2-го порядка
- •Краткие сведения теории ЛНДУ 2-го порядка
- •Модель гонки вооружений Ричардсона
- •Модель ведения боевых действий Ланчестера
- •Числовые ряды
- •Понятие числового ряда и его сходимости
- •Некоторые примеры
- •Знакопеременные ряды
- •Знакочередующиеся ряды. Признак Лейбница
- •Свойства степенных рядов
- •Примеры
- •Приложения рядов

© БГЭУ Лекция № 9 Дифференциальные уравнения 2-го порядка проф. Дымков М.П. 1
Краткие сведения теории ЛНДУ 2-го порядка
Определение Дифференциальное уравнение вида y′′+ py′+ qy = f (x), (1)
где p, q—константы, f (x)—заданная функция называется линейным неоднородным дифференциаль-
ным уравнением второго порядка с постоянными коэффициентами.
Как устроено общее решение ЛНДУ ?
Теорема (структура общего решения ЛНДУ).
Общее решение линейного неоднородного дифференциального уравнения ЛНДУ (1) есть сумма
|
(*) |
y(x) = y *(x) + yобщ(x) |
некоторого частного решения y *(x) неоднородного ДУ
(1) и общего решения |
yобщ(x) = C1 y1(x) +C2 y2 (x) |
||||
однородной части ДУ (1), т.е. решения для ДУ вида |
|||||
|
|
y′′+ py′+ qy = 0. |
(**) |
||
|
|
|
y2 (x) линейно независимые решения для |
||
|
Здесь y1(x), |
||||
|
ДУ (**), а |
C1,C2 —произвольные постоянные. |
Как искать частное решение неоднородного ДУ?
Укажем для специального вида правой части f (x)!!

© БГЭУ Лекция № 9 Дифференциальные уравнения 2-го порядка проф. Дымков М.П. 2
1.Пусть
f (x) =Pn(x) ≡a0xn +a1xn−1 +...+an−1x+an
а) Предположим, что число λ = 0 не является корнем
характеристического уравнения |
2 + p λ+ q = λ0 |
Тогда частное решение для неоднородного уравнения ищем тоже в виде многочлена той же степени, (но с неизвестными пока коэффициентами):
y*(x) =Rn(x) ≡b0xn +b1xn−1 +...+bn−1x+bn (а1)
Подставим (а1) в ДУ (1):
|
для x |
(в1) |
Rn′′(x) + pRn′ (x) + qRn (x) ≡ Pn (x) |
Сравнивая в тождестве (в1) коэффициенты при одинаковых степенях x , получим алгебраическую систему из (n +1) уравнений относительно неизвестных коэффициентов bi
б) Пусть число λ = 0 является корнем (однократным !)
характеристического уравнения |
2 + p λ+ q = λ0 |
(в этом случае оно имеет вид λ |
(λ + p) = 0)! |
Тогда частное решение для неоднородного уравнения
ищем в виде |
|
y*(x) =x Rn(x) |
(а2) |
Далее используется процедура пункта а).

© БГЭУ Лекция № 9 Дифференциальные уравнения 2-го порядка проф. Дымков М.П. 3
в) Пусть число λ = 0 является корнем (двухкратным !)
характеристического уравнения |
λ2 + pλ + q = 0 |
(в этом случае оно имеет вид λ2 = 0)!
Тогда частное решение для неоднородного уравнения
ищем в виде |
|
y*(x) =x2 Rn(x) |
(а3) |
Далее используется процедура пункта а).
Пример. Найти общее решение дифференциального
уравнения y′′− 2 y′+ 2 y = x2.
Решение. 1)Ищем общее решение для y′′− 2 y′+ 2 y = 0.
Характеристическое уравнение: λ2 − 2λ + 2 = 0, |
|
|||
D = 4 −8 = −4 = 4i2 , λ = |
2 ± 2i |
=1 ± i, |
α =1, |
β =1, |
|
||||
1,2 |
2 |
|
|
|
|
|
|
|
Значит, общее решение однородного уравнения есть yобщ = ex (C1 sin x +C2 cos x).
2) Для нахождения частного решения неоднородного уравнения используем специальный вид правой части
Pn (x) = x2 , n = 2.
Так как среди корней характеристического уравнения нет λ = 0, то множители x и x2 отсутствуют, значит, ищем частное решение y в виде многочлена второй степени с неопределенными коэффициентами
y = ax2 + bx + c, y′ = 2ax + b, |
y′′ = 2a. |

© БГЭУ Лекция № 9 Дифференциальные уравнения 2-го порядка проф. Дымков М.П. 4
Подставляем y , y′ , y′′ в исходное уравнение y′′ − 2 y′ + 2 y = x2
Имеем
2a − 4ax − 2b + 2ax2 + 2bx + 2c = x2;
x2 (2a) + x(−4a + 2b) + 2a − 2b + 2c = x2.
Составляем систему для нахождения a, b и c, приравнивая коэффициенты при соответствующих степенях x в левой и правой части уравнения:
x2 : 2a =1, |
|
|
|
|
|
a =1 2, |
|||||||
|
1 |
|
|
|
|
|
|
b =1, |
|
|
|||
x |
: − 4a + 2b = 0, |
|
|
|
|
|
|
||||||
x |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
: 2a − 2b + 2c = 0. |
|
|
|
|
|
c =1 2. |
||||||
Итак, частное решение |
y |
|
= |
1 |
x |
2 |
+ x + |
1 |
|
. |
|||
|
2 |
|
2 |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Следовательно, общее решение неоднородного ДУ есть
Y= yобщ + y = ex (C1 sin x +C2 cos x) + 12 x2 + x + 12 .
2.Пусть правая часть имеет вид
f (x) =ekx Pn(x) ≡ekx (a0xn +a1xn−1 +...+an−1x+an)
Важную роль играет взаимодействие корней характеристического уравнения с правой частью!!!

© БГЭУ Лекция № 9 Дифференциальные уравнения 2-го порядка проф. Дымков М.П. 5
Пусть λi ,i =1,2 |
— корни характеристического |
уравнения |
λ2 + pλ + q = 0 |
Возможны случаи:
а) k ≠ λi
Тогда частное решение ищем в виде y*(x) =ekx Rn(x)
где Rn (x)— многочлен n-ой степени с неизвестными
коэффициентами.
б) число k является простым корнем (кратности один) характеристического уравнения,
т.е. k = λ1 или k = λ2
Частное решение ищем в виде y*(x) =x ekx Rn(x) где Rn (x)— многочлен n-ой степени с неизвестными коэффициентами.
в) число k является двукратным корнем характеристического уравнения, т.е. k = λ1 = λ2 = λ
Частное решение ищем в виде y*(x) =x2 ekx Rn(x)
где Rn (x)— многочлен n -ой степени с неизвестными коэффициентами.
Пример. Найти общее решение дифференциального
уравнения y′′− 2 y′+ y = xex.
Решение. Характеристическое уравнение имеет вид
λ2 − 2λ +1 = 0, (λ −1)2 = 0, |
λ1 = λ2 =1. |
|||
Значит, общее решение есть |
y = (C + C |
2 |
x)ex . |
|
|
1 |
|
|
© БГЭУ Лекция № 9 Дифференциальные уравнения 2-го порядка проф. Дымков М.П. 6
Так как правая часть имеет вид xe1 x и k =1 совпадает с |
|
||||||||||||||||||||||||
λ1 = λ2 =1, |
|
|
то |
|
в |
частном решении |
|
появляется |
|
||||||||||||||||
множитель x2: |
|
|
|
|
|
|
y = x2 (ax + b)ex. |
|
|
|
|||||||||||||||
Тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y = ex (ax3 +bx2 ); |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
′ |
x |
(ax |
3 |
+bx |
2 |
) +e |
x |
(3ax |
2 |
+ 2bx) |
= e |
x |
(ax |
3 |
+(b +3a)x |
2 |
+ 2bx); |
|||||||
|
y * = e |
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
′′ |
|
x |
(ax |
3 |
+(b +3a)x |
2 |
+ 2bx) +e |
x |
(3ax |
2 |
+ 2(b +3a)x |
+ 2b = |
||||||||||||
|
y * = e |
|
|
|
|
|
|
||||||||||||||||||
|
= ex (ax3 +(b +6a)x2 +(4b +6a)x + 2b) |
|
|
|
|
|
|||||||||||||||||||
|
Подставляем y , y ″, y ″ в исходное уравнение, |
|
|||||||||||||||||||||||
группируя слагаемые по степеням |
x |
|
и вынося ex |
за |
|
||||||||||||||||||||
скобки: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ex (ax3 −2ax3 + ax3 +bx2 −2(b +3a)x2 +(b +6a)x2
−4bx +(4b +6a)x + 2b) = xex .
Приравнивая коэффициенты при x3, x2 , x и x0в левой и правой частях последнего уравнения, получаем систему для нахождения неопределенных коэффициентов a и b:
x |
3 |
: |
0a = 0, |
|
|
|
|
0a = 0, |
|
|
|
|
|
||||||
|
|
|
|
|
+6a = |
0, |
|
|
|
|
|
|
1 |
|
|
||||
x2 : |
b −2b −6a +b |
0a +0b = 0, |
|
a = |
, |
b = 0. |
|||||||||||||
x : |
|
−4b +b +16a |
=1, |
|
|
|
6a =1, |
6 |
|||||||||||
|
|
|
|
|
|
|
|
||||||||||||
x |
0 |
: |
|
2b = 0 |
|
|
|
|
|
2b = 0 |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
Значит, y = |
x3 |
ex. |
|
|
Общее решение есть |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
Y = (C + C |
2 |
x)ex + |
x3 |
ex . |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
1 |
|
|
6 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|

© БГЭУ Лекция № 9 Дифференциальные уравнения 2-го порядка проф. Дымков М.П. 7
3. Трудный случай представляет правая часть вида
f (x) =ekx(Pn(x) cospx+Qm(x) sinpx)
где Pn (x) — многочлен n-ой степени, Qm (x) —
многочлен m-ой степени.
а) Пусть комплексное число µ = k ± p i не является корнем характеристического уравнения
λ2 + pλ + q = 0
Тогда частное решение ищем в виде
y*(x) =ekx(Rq(x) cospx+Tq(x) sinpx)
где степень q = max{m, n}
б) Пусть комплексное число µ = k ± p i является корнем характеристического уравнения
λ2 + pλ + q = 0
Тогда частное решение ищем в виде
y*(x) =x ekx(Rq(x) cospx+Tq(x) sinpx)
где степень q = max{m, n}
Пример. Найти решение ДУ
y′′+ y′− 2 y = cos x −3sin x,
удовлетворяющее начальным условиям y(0) =1, y′(0) = 2.
Решение. Прежде находим общее решение
однородного уравнения
y′′+ y′− 2 y = 0.

© БГЭУ Лекция № 9 Дифференциальные уравнения 2-го порядка проф. Дымков М.П. 8
Характеристическое уравнение :
λ 2+ λ − 2 = 0, D =1+ 4 2 = 9 = 32 , λ 1= −2, λ 2= +1. (*)
Значит, общее решение однородного уравнения есть y = C1e−2x +C2ex.
Корней вида k + pi = i в уравнении (*) нет.
Значит, частное решение y неоднородного уравнения будем искать в виде y = a cos x + b sin x.
Тогда y ′ = −a sin x + b cos x, y ″ = −acos x − bsin x.
Подставляем y , y ″, y ″ в исходное уравнение и, приводя подобные члены при cos x и sin x, получаем:
−a cos x −b sin x −a sin x +b cos x + a cos x +b sin x
=cos x −3sin x,
(b −3a) cos x + (−3b − a) sin x ≡ cos x −3sin x .
Приравнивая коэффициенты при cos x и при sin x в правой и левой частях последнего уравнения, получим:
b −3a =1, |
3b − 9a = 3, |
−10a = 0, a = 0,b =1 |
|
|
|
−3b −a = −3 |
−3b − a = −3, |
|
Значит, частное решение неоднородного уравнения имеет вид y = sin x, а общее решение неоднородного
уравнения есть |
Y = C e−2x + C |
2 |
ex + sin x; |
|
|
|
1 |
|
|
|
|
|
|
|
′ |
= 2. |
|
Учтем начальные условия y(0) =1, y (0) |
Найдем Y ′ = −2C1e−2x + C2ex + cos x.

© БГЭУ Лекция № 9 Дифференциальные уравнения 2-го порядка проф. Дымков М.П. 9
Отсюда |
|
1 = C1 +C2 и |
2 = −2C1 + C2 +1. |
|
|
|
||||||||||||
Решаем полученную систему уравнений: |
|
|
|
|
||||||||||||||
C |
+ C |
2 |
=1, |
|
|
C |
= 0, |
|
|
|
|
|
|
|
|
|||
|
|
Ответ: |
Y = ex + sin x. |
|
|
|||||||||||||
1 |
|
|
|
1 |
|
|
|
|||||||||||
− 2C1 + C2 =1, |
|
C2 =1. |
|
|
|
|
|
|
|
|
||||||||
|
|
Сводная таблица: |
|
|
|
|
|
|
|
|
|
|
||||||
Пусть функция |
f (x) имеет вид |
Pn (x), |
или |
Pn (x)eαx , |
||||||||||||||
или |
|
|
(Pn |
(x)sin βx+ Pn |
(x)cos βx)eαx , |
|
|
|
||||||||||
где Pn (x), Pn |
1 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|||||
(x), Pn |
(x) − многочлены степени n или не |
|||||||||||||||||
|
|
|
|
|
1 |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
меньше n . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Тогда частное решение можно найти в виде |
|
|
|
|||||||||||||||
|
Qn (x) , |
Qn (x)eαx , eαx (Qn |
(x)sin βx+Qn |
(x) sin βx), |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
2 |
|
|
|
где |
Qn (x), |
Qn (x), Qn |
(x) − многочлены |
такой |
же |
|||||||||||||
|
|
|
|
|
|
|
1 |
|
|
2 |
|
|
|
|
|
|
|
|
степени, |
что |
|
и |
Pn (x), Pn |
(x), Pn |
2 |
(x), |
но |
с |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
неопределенными коэффициентами, или в таком же виде, но с множителем x или x2 в зависимости от
соотношения корней |
λ1 и λ2 характеристического |
|||||
уравнения и числа α (α ± iβ) : |
|
x2 |
||||
1. если |
λ ≠α, λ |
2 |
≠ α, |
то множители x |
и |
|
|
1 |
|
|
|
|
отсутствуют;
2.если λ1 =α, (или λ2 =α), но λ1 ≠ λ2 , то появляется множитель x ;
3.если λ1 = λ2 = λ, то появляется множитель x2;
4. если ни один их корней λ1,2 =α1 ± β1 i характеристического уравнения не равен α ± β i , то множитель x отсутствует; если λ1 =α, β1 = β , то появится множитель x .

© БГЭУ Лекция № 10 |
Приложения ДУ в экономике |
проф. Дымков М.П. 1 |
||||
|
Найти общее решение дифференциального |
|||||
Пример. |
||||||
уравнения |
y′′−2y′+ y = |
3e2x |
+ |
xex |
(Super) |
Замечание (принцип суперпозиции решений ДУ)
|
Пусть правая часть ДУ |
y′′+ py′+ qy = |
|
f (x) |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
(*) |
|
|
|||||||||||||||||||||||||||||
|
равна сумме двух функций |
|
|
|
|
|
|
f (x) = f1(x) + f2 ( |
x) |
, |
||||||||||||||||||||||
|
а y1(x), y2 (x) |
есть решения с одной |
левой частью, |
но |
с |
|||||||||||||||||||||||||||
|
разными правыми частями |
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
y′′+ py′+ qy = f1(x) |
|
y′′+ py′+ qy = f2 (x) |
|
|
|
|
|
|
||||||||||||||||||||||
|
УПР |
*(доказать) |
|
|
|
|
|
|
|
есть решение |
|
|
|
|
|
|
|
|
||||||||||||||
|
|
Тогда |
y |
(x) = y1(x) + y2 (x) |
(*) |
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
Решение. В нашем случае |
|
f |
(x) = 3e2x , f |
2 |
(x) = xex |
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
1) Решение |
|
y′′−2y′+ y = |
|
xex |
( для |
f2 (x) = xex ) |
|||||||||||||||||||||||||
|
|
|
|
имеет вид |
y |
|
(x) = (C |
+C |
|
|
x)ex + |
x3 |
ex . |
|||||||||||||||||||
|
(см. Лекция 9) |
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
1 |
|
|
2 |
|
|
|
|
6 |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
2) Для f1(x) = 3e2x имеем ДУ |
|
|
|
y′′−2y′+ y = 3e2x |
|||||||||||||||||||||||||||
|
Характеристическое уравнение имеет вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
λ2 − 2λ +1 = 0, (λ −1)2 = 0, |
|
λ1 = λ2 =1. |
||||||||||||||||||||||||||||
|
Значит, общее решение есть |
|
|
y = (C + C |
2 |
x)ex . |
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Так как правая часть имеет вид |
|
|
3e2 x |
|
|
и |
|
|
k = 2 |
|
|
не |
|||||||||||||||||||
совпадает с λ1 = λ2 =1, то |
|
частное решение следует |
||||||||||||||||||||||||||||||
искать в виде: |
|
|
|
|
|
y = ae2x . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Тогда

© БГЭУ Лекция № 10 |
Приложения ДУ в экономике |
|
проф. Дымков М.П. |
2 |
|||||||
y |
|
= e |
x |
(ax +b); |
′ |
x |
(ax +b) + e |
x |
a; |
′′ |
2x |
|
|
y * = e |
|
|
y * = 4ae |
|
Подставляем y , y ″, y ″ в исходное уравнение,
группируя слагаемые по степеням x и вынося ex за скобки:
e2x (4ax − 4a + a) = 3e2x . a = 3 y *(x) = 3e2x
Общее решение есть y1(x) = (C1 +C2 x)ex +3e2x .
Итак, по принципу суперпозиции общее решение для
ДУ(Super) |
равно |
|
|
x3 |
|
|||
y(x) = y (x) + y |
2 |
(x) = (C +C |
2 |
x)ex +3e2x + |
ex . |
|||
|
||||||||
1 |
|
1 |
6 |
|
||||
|
|
|
|
|
|
|
|
|
y′′+ py′+ qy = |
|
f (x), |
|
(1) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Как быть, если f (x) −произвольная функция? |
|
|
||||||||
|
|
|
( p, q ≠ const ?) |
|
|
|
|
|
|
|
||
|
Для нахождения частного |
решения |
y (x) |
|
можно |
|||||||
|
использовать метод вариации |
произвольных |
||||||||||
|
постоянных (метод Лагранжа) |
|
|
|
|
|
|
|
||||
|
|
однородное уравнение |
|
|
|
|
|
|||||
|
Пусть |
(2) |
||||||||||
|
|
y′′+ py′+ qy = 0 |
|
|||||||||
|
имеет |
общее решение |
|
|
|
|
|
|
|
|
|
|
|
|
|
y(x) = C1 y1(x) +C2 y2 (x) |
(3) |
|
Будем искать решение |
неоднородного |
|
уравнения |
(1) |
|||
|
в виде (метод Лагранжа |
! ) |
|
|
|
|
|
|
|
|
y(x) = C1(x) y1(x) +C2 (x) y2 (x) |
|
|
(*) |
|
||
|
где C1(x |
),C2 (x) −пока неизвестные функции. |
|
|
|
|||
|
|
|
|
|
|
|
|
|

© БГЭУ Лекция № 10 Приложения ДУ в экономике |
проф. Дымков М.П. 3 |
Дифференцируем |
|
y′(x) = C1′(x) y1(x) +C1 (x) y1′(x) +C2′ (x) y2 (x) +C2 (x) y2′ (x)
|
(это можем !), чтобы выражение для y′ |
|
|
|
|
|
|||||||||||||||||||
Потребуем |
|
|
|
|
|
||||||||||||||||||||
имело такой |
же вид как и для y . |
Для этого положим |
|
|
|||||||||||||||||||||
Тогда имеем |
|
|
C1′(x) y1(x) +C2′ (x) y2 (x) = 0 |
. |
|
|
|
|
(А) |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
y′(x) = C1 (x) y1′(x) +C2 (x) y2′ (x) |
|
|
|
|
|
|
|
|
|||||||||||||
Дифференцируем полученное выражение |
|
|
|
|
|
|
|
|
|||||||||||||||||
′′ |
′ |
|
|
|
′ |
(x) |
|
|
|
|
′′ |
|
′ |
(x) y |
′ |
(x) +C |
|
(x) y |
′′ |
(x) |
|||||
y (x) |
= C |
(x) y |
+C (x) y (x) +C |
|
|
|
2 |
|
|||||||||||||||||
|
1 |
|
|
1 |
|
|
|
1 |
1 |
|
2 |
|
2 |
|
|
2 |
|
|
|
|
|
||||
Найденные |
|
|
′ |
|
|
′′ |
|
подставим в исходное ДУ |
(1): |
|
|
|
|||||||||||||
|
y, y , y |
|
|
|
|
|
|||||||||||||||||||
y′′+ py′+ qy =C1′y1′ +C1 y1′′+C2′ y2′ |
+C2 y2′′+ |
|
|
|
|
|
|
|
|
+ pC1 y1′ + pC2 y2′ +qC1 y1 + qC2 y2 =(перегруппируем)= =C1′y1′ +C2′ y2′ +C1( y1′′+ py1′ + qy1)+C2 ( y2′′ + py2′ + qy2 )= ( y1, y2 − решения !) = C1′y1′ +C2′ y2′ +C1 0 +C2 0 = f (x)
Итак, получили систему уравнений
C1′y1 +C2′ y2 = 0C1′y1′ +C2′ y2′ = f (x)
Определитель системы(вронскиан) W (x) =
Следовательно, можно найти C1′(x),C2′ (x)
(**)
y1 |
y2 |
|
≠ 0 |
|
|||
y′1 |
y2′ |
|
|
,

© БГЭУ Лекция № 10 Приложения ДУ в экономике |
проф. Дымков М.П. 4 |
|
Проинтегрировав их найдем |
искомые |
функции |
C1(x),C2 (x) и тем самым найдем требуемое решение неоднородного ДУ :
y(x) = C1(x) y1(x) +C2 (x) y2 (x)
|
|
|
|
|
|
|
|
Пример |
y′′+ y = tgx |
|
|
|
|
|
|
Найдем общее решение однородного ДУ |
y′′+ y = 0 |
||||||
Характеристическое уравнение λ2 +1 = 0,λ |
|
= ±i . |
|||||
|
1,2 |
|
|
|
|||
Общее решение ОДУ |
y(x) = e0 x (C sin x +C |
2 |
cos x) |
|
|||
|
|
1 |
|
|
|
По методу Лагранжа будем искать частное решение для неоднородного ДУ в виде
y(x) = C1(x) sin x +C2 (x) cos x)
где C1(x),C2 (x) -- подлежащие определению функции.
Система уравнений (**) в этом случае имеет вид
C1′ cos x +C2′ sin x = 0
−C1′ sin x1 +C2′ cos x = tgx
Отсюда найдем C1′ |
= − |
sin2 x |
, |
C2′ = sin x |
|
|
|
|
|||||
cos x |
|
|
|
|
|||||||||
Интегрируем: |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
||||
C1 = −∫ |
sin2 x |
dx = ∫ |
(cos x − |
1 |
|
)dx |
= sin x −ln tg( |
π |
+ |
x |
) + K1 |
||
cos x |
cos x |
4 |
2 |
||||||||||
|
|
|
|
|
|
|
|
C2 = ∫sin xdx = −cos x + K2
Окончательно, общее решение неоднородного ДУ есть
y(x) = C1(x) sin x +C2 (x) cos x=
=(sin x −ln tg(π4 + 2x) + K1) sin x + (−cos x +K 2) cos x

© БГЭУ Лекция № 10 Приложения ДУ в экономике |
проф. Дымков М.П. 5 |
Если одно дифференциальное уравнение описывает изменение во времени состояния одного объекта, то система подобных уравнений описывает динамику состояния системы взаимосвязанных объектов.
Однородные системы
Однородной называется система следующего вида:
|
dy1 |
= a |
y |
....+ + a |
y |
n |
|
||
|
|
||||||||
dt |
11 |
1 |
1n |
|
|
||||
|
|
|
|
|
|
|
|
(0.1) |
|
.......... |
|
|
.......... |
|
.................... |
|
|
||
dy |
n |
|
|
|
|
|
|
||
|
|
= an1 y1 |
....+ + ann yn |
|
|||||
|
|
|
|||||||
dt |
|
|
|
|
|
|
В матричной форме систему можно переписать следующим образом:
dy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
a ....a |
y |
|
|
|
|
|
|
|||
dt |
|
|
|
|
|
|||||||||
|
|
|
|
11 |
1n |
|
1 |
|
или |
d y |
= Ay |
|||
.... |
|
|
= |
................. |
.... |
|||||||||
|
dt |
|||||||||||||
|
dyn |
|
|
|
|
|
|
|
|
|
|
|||
|
|
an1....ann |
yn |
|
|
|
|
|
||||||
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где br - |
||||||||||
Будем искать решение |
в виде (Эйлер): |
|
y =b eλt |
вектор постоянных множителей. Подставляя это решение в (0.1), получаем:
bλ eλt = Abeλt ( A −λE)beλt = 0 (eλt ≠ 0) ( A −λE)b = 0

© БГЭУ Лекция № 10 Приложения ДУ в экономике |
проф. Дымков М.П. 6 |
То есть мы получаем так называемую задачу о
собственных числах и собственных векторах матрицы A .
Цель—определить такие λ и b, при которых задача имеет нетривиальное решение (b≠0).
Решая матричное характеристическое уравнение: det(A −λE)= 0 получаем собственные числа λi .
Далее можно обычным способом находить для каждого найденного собственного числа соответствующий ему
собственный вектор bi , (i =1,2,..., n).
Тогда общее решение системы ДУ имеет вид
|
|
|
y(t) =C b eλ1t |
+C b eλ2t |
+....+C b eλnt |
|
|||||
|
|
|
|
|
1 1 |
2 |
2 |
|
n n |
|
|
|
|
|
|
|
|
|
|
dx |
= x + 2 y |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||
|
Пример 1 |
. Решить систему |
dt |
|
|
||||||
|
|
|
|
|
|
|
|
|
dy |
= 2x + y |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
Решим характеристическое уравнение: |
|
||||||||||
|
|
1−λ |
2 |
|
= 0 (1−λ)2 − 4 = 0 λ |
2 − 2λ −3 = 0 λ = 3, λ = −1 |
|||||
|
|
|
|||||||||
|
|
2 |
1−λ |
|
|
|
|
|
1 |
2 |
|
|
|
|
|
|
|
|
|