Однородные и неоднородные системы линейных дифференциальных уравнений (стр. 1 из 4)
Федеральное Агентство по образованию
государственное Образовательное Учреждение высшего профессионального образования
«Тюменский Государственный Нефтегазовый Университет»
Институт Нефти и Газа
Кафедра «математические методы в экономике»
Курсовая работа
по математическому анализу
Однородные и неоднородные системы линейных дифференциальных уравнений
Проверил: старший преподаватель
Тюмень 2007
Содержание
ВВедение
1 Системы линейных дифференциальных уравнений.
1.1 Общие сведения о линейных системах.
1.2 Метод сведения линейной системы к одному уравнению более высокого порядка.
1.3 Методы решения однородных линейных систем дифференциальных уравнений.
1.4 Методы решения неоднородных линейных систем дифференциальных уравнений.
2. Решение линейных систем дифференциальных уравнений.
2.1.Решение методом сведения линейной системы к одному уравнению более высокого порядка.
2.2. Решение однородных линейных систем дифференциальных уравнений.
2.2.1. Решение видоизмененным методом Эйлера
2.3.2. Решение методом неопределенных коэффициентов
ВВЕдение
1. Системы линейных дифференциальных уравнений.
1.1 Общие сведения о линейных системах.
Линейные системы – это системы дифференциальных уравнений вида
(1)Где коэффициенты aij и fi – некоторые функции независимой переменной x. Будем считать их непрерывными; тогда для данной системы заведомо выполняются условия теоремы о существование и единственности решения задачи Коши. Если все fi=0, то система называется однородной, в противном случае она называется неоднородной. Система
(2)Называется однородной системой, соответствующей неоднородной системе (1).
При изучении линейных систем удобно использовать матричные обозначения
Позволяющие записать систему (1) в виде одного матричного уравнения
(3)Так же, как и в случае линейных уравнений, общее решение неоднородной системы представляет собой сумму частного решения этой системы и общего решения соответствующей ей однородной системы. В свою очередь, общее решение однородной системы имеет вид
(4)Где С1,…,Сn- произвольные постоянные, а
-произвольные линейно независимые решения, называемые фундаментальным набором решений этой системы. Критерием линейной независимости этих решений является неравенство нулю определителя Вронского
(5)1.2 Метод сведения линейной системы к одному уравнению более высокого порядка.
(Этот метод применим как для однородной, так и для неоднородной системы линейных дифференциальных уравнений.)
Один из методов интегрирования линейной системы заключается в сведении системы к одному уравнению n-ого порядка с одной неизвестной функцией. Продемонстрируем это на примере системы двух уравнений.
(6)Дифференцируя (по x) обе части первого уравнения системы (6), находим
откуда, заменяя производные y1', y2' их выражениями из самой системы, имеем
.Группируя в правой части, получим уравнение вида
(7)Где коэффициенты b1, b2 и d1 определенным образом выражаются через коэффициенты aij и q1 и их производные. Комбинируя уравнение (7) с первым уравнением системы (6), получим
(8)Предположим, что в рассматриваемой области изменения x определитель
отличен от нуля. Тогда систему (8) можно решить относительно y1 и y2, т.е. выразить y1и y2 через y’1 и y”2.В результате приходим к уравнениям вида
(9) 
. (10)Первое из них представляет собой линейное дифференциальное уравнение второго порядка с одной неизвестной функцией y1(t). Заметим, что если в исходной системе (6) все коэффициенты aij постоянны, то уравнение (9) также является уравнением с постоянными коэффициентами. [ 3 стр 509-510]
1.3 Методы решения однородных линейных систем дифференциальных уравнений.
1) Сведение к одному уравнению n-ого порядка. (Этот метод мы разбирали выше)
2) Решение ЛОСДУ с постоянными коэффициентами с помощью матриц (видоизмененный метод Эйлера).
Пусть дана система n линейных дифференциальных уравнений с n неизвестными функциями, коэффициенты которой постоянные:
(11)Эту систему можно записать в виде одного матричного дифференциального уравнения
.Здесь
[2 стр 169]
Ищем решение системы в виде:
X1=α1ekt , X2=α2ekt,……. Xn=αnekt . (12)
Требуется определить постоянные α1, α2,…, αn и k так, чтобы функции α1ekt, α2ekt,…,αnekt удовлетворяли системе уравнений (11). Подставим их в систему(1), получим:
Сократим на ekt. Перенося все члены в одну сторону и собирая коэффициенты при α1, α2,…., αn , получим систему уравнений
Выберем α1, α2,…., αn и k такими, чтобы удовлетворялась система (13).Эта система есть система линейных однородных алгебраических уравнений относительно α1, α2,…., αn. Составим определитель системы (13):
(14)Если k таково, что определитель ∆ отличен от нуля, то система (13) имеет только нулевые решения α1=α2=…=αn=0,а следовательно, формулы (12) дают только тривиальные решения:
X1(t)=X2(t)=…=X(t)=0.
Таким образом, нетривиальные решения (12) мы получим только при таких k,при которых определитель (14) обращается в нуль. Мы приходим к уравнению n-ого порядка для определения k:
(15)Это уравнение называется характеристическим уравнением для системы (1),его корни называются корнями характеристического уравнения.
Рассмотрим несколько случаев.
Случай 1. Корни характеристического уравнения действительны и различны. Обозначим через k1, k2,….kn корни характеристического уравнения. Для каждого корня kj напишем систему (13) и определим коэффициенты
α1(i),α2(i),…, αn(i).
Можно показать, что один из них произвольный, его можно считать равным единице. Таким образом, получаем:
для корня k1 решение системы (11)
Для корня k2 решение системы (1)
……………………………………………………….
для корней kn решение системы (1)
Путем непосредственной подстановки в уравнения можно убедиться, что система функций
(16)где С1, С2,….,Сn -произвольные постоянные, тоже является решением системы дифференциальных уравнений (11). Это есть общее решение системы (11). Легко показать, что можно найти такие значения постоянных, при которых решение будет удовлетворять заданным начальным условиям.
Случай 2. Корни характеристического уравнения различные, но среди них есть комплексные. Пусть среди корней характеристического уравнения имеется два комплексных сопряженных корня:
k1 = α+iβ, k2 = α-iβ.
Этим корням будут соответствовать решения
(j = 1, 2, …,n), (17)