Лекция 9. Простейшие методы интегрирования: метод исключения, метод интегрируемых комбинаций.

                  Системы линейных дифференциальных уравнений с постоянными коэффициентами.

 

Пример 1.

Решить систему дифференциальных уравнений Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image024.gif 

 Решение:  Решим систему методом исключения.

Алгоритм решения:

1) Берем второе уравнение системы Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image030.gif и выражаем из него Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image032.gif:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image034.gif   

2) Дифференцируем по Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image010_0000.gif обе части полученного уравнения Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image037.gif:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image039.gif

Со «штрихами» процесс выглядит так:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image041.gif
3) Подставим Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image037_0000.gif и Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image039_0000.gif в первое уравнение системы Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image043.gif:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image045.gif

И проведём максимальные упрощения:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image047.gif

Получено однородное уравнение второго порядка с постоянными коэффициентами. Со «штрихами» оно записывается так: Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image049.gif.

Составим и решим характеристическое уравнение:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image051.gif
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image053.gif – получены различные действительные корни, поэтому:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image055.gif.

4) Идём за функцией Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image006_0002.gif. Для этого берём уже найденную функцию Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image058.gif и находим её производную. Дифференцируем по Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image010_0001.gif:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image060.gif

Подставим Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image058_0000.gif и Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image062.gif в уравнение (*):

Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image064.gif

Или короче: Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image066.gif

5) Обе функции найдены, запишем общее решение системы:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image068.gif

Неоднородная система дифференциальных уравнений имеет следующий вид:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image121.gif

По сравнению с однородной системой в каждом уравнении дополнительно добавляется некоторая функция, зависящая от «тэ». Функции Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image123.gif могут быть константами (причем, по крайне мере одна из них не равна нулю), экспонентами, синусами, косинусами и т.д.

Пример 2. Найти частное решение системы линейных ДУ, соответствующее заданным начальным условиям
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image125.gif

Решение: Дана линейная неоднородная система дифференциальных уравнений, в качестве «добавок» Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image123_0000.gif выступают константы. Используем метод исключения.

1) Из первого уравнения системы выражаем:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image127.gif  

2) Дифференцируем по Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image010_0002.gif обе части:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image129.gif

3) Подставим Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image131.gif и Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image129_0000.gif во второе уравнение системы Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image134.gif:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image136.gif

Сразу после подстановки целесообразно избавиться от дробей, для этого каждую часть уравнения умножаем на 5:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image138.gif

Теперь проводим упрощения:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image140.gif

В результате получено линейное неоднородное уравнение второго порядка с постоянными коэффициентами

Найдем общее решение соответствующего однородного уравнения:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image142.gif

Составим и решим характеристическое уравнение:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image144.gif
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image146.gif – получены сопряженные комплексные корни, поэтому:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image148.gif.

Частное решение неоднородного уравнения ищем в виде Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image150.gif.
Найдем первую и вторую производную:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image152.gif

Подставим Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image154.gif в левую часть неоднородного уравнения:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image156.gif

Таким образом: Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image158.gif

В результате: Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image161.gif

4) Ищем функцию Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image008_0003.gif. Сначала находим производную от уже найденной функции Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image006_0004.gif:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image165.gif

Подставим Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image167.gif 
и Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image169.gif в уравнение (*):
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image171.gif

5) Общее решение системы:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image173.gif

 

 Рассмотрим метод характеристического уравнения (метод Эйлера) на примере:

Пример 3. Дана линейная однородная система дифференциальных уравнений

Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image187.gif

Найти общее решение системы уравнений с помощью характеристического уравнения

Решение: Смотрим на систему уравнений и составляем определитель второго порядка:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image189.gif
Составим характеристическое уравнение, для этого из каждого числа, которое располагается на главной диагонали, вычитаем некоторый параметр Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image191.gif:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image193.gif

Раскрываем определитель:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image195.gif

Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image197.gif

Если характеристическое уравнение имеет два различных действительных корня, то общее решение системы дифференциальных уравнений имеет вид:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image199.gif

Коэффициенты в показателях экспонент Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image201.gif нам уже известны, осталось найти коэффициенты Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image203.gif

1) Рассмотрим корень Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image205.gif и подставим его в характеристическое уравнение:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image207.gif
Из чисел  определителя составим систему двух линейных уравнений с двумя неизвестными:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image209.gif

Из обоих уравнений следует одно и то же равенство:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image211.gif

Теперь нужно подобрать наименьшее значение Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image213.gif, такое, чтобы значение Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image215.gif было целым. Очевидно, что следует задать Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image217.gif. А если Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image217_0000.gif, то Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image220.gif

2) Аналогично рассмотрим корень Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image222.gif и устно подставим его в характеристическое уравнение:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image224.gif

Из чисел  определителя составим систему:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image226.gif

Из обоих уравнений следует равенство:
Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image228.gif

Подбираем наименьшее значение Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image230.gif, таким образом, чтобы значение Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image232.gif было целым. Очевидно, что Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image234.gif.

Все четыре коэффициента Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image203_0000.gif найдены, осталось их подставить в общую формулу Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image236.gif

Ответ: общее решение: Описание: http://mathprofi.ru/h/sistemy_differencialnyh_uravnenij_clip_image238.gif