Метод исключения
Нормальная система обыкновенных дифференциальных уравнений (СОДУ), являющаяся линейной однородной с постоянными коэффициентами, имеет следующий вид: {y′1=a11⋅y1+a12⋅y2+…+a1n⋅yny′2=a21⋅y1+a22⋅y2+…+a2n⋅yn…y′n=an1⋅y1+an2⋅y2+…+ann⋅yn.
Здесь y1(x),y2(x),…,yn(x) -- искомые функции независимой переменной x, коэффициенты ajk,1≤j,k≤n -- заданные действительные числа.
Для решения СОДУ такого вида применим метод исключения, состоящий в преобразовании её в одно дифференциальное уравнение (ДУ) n-го порядка, которое затем решим каким-либо из известных методов.
Решить систему ДУ {dy1dx=2⋅y1+y2dy2dx=3⋅y1+4⋅y2.
Систему решаем исключением неизвестной функции y2.
Шаг 1. Из первого уравнения находим y2: y2=dy1dx−2⋅y1.
Шаг 2. Подставляем y2 во второе уравнение:
dy2dx=3⋅y1+4⋅(dy1dx−2⋅y1);dy2dx=4⋅dy1dx−5⋅y1.Шаг 3. Дифференцируем первое уравнение по x: d2y1dx2=2⋅dy1dx+dy2dx.
Шаг 4. Подставляем выражение, полученное на шаге 2, в выражение, полученное на шаге 3:
d2y1dx2=2⋅dy1dx+4⋅dy1dx−5⋅y1;d2y1dx2−6⋅dy1dx+5⋅y1=0.Шаг 5. Решаем линейное однородное ДУ второго порядка с постоянными коэффициентами:
- характеристическое уравнение k2−6⋅k+5=0;
- корни характеристического уравнения k1=1, k2=5 -- действительные, разные;
- искомая функция y1=C1⋅ex+C2⋅e5⋅x.
Шаг 6. Находим функцию y2:
- производная dy1dx=C1⋅ex+5⋅C2⋅e5⋅x;
- результат подстановки в выражение, полученное на шаге 1:
Общее решение данной системы:
y1=C1⋅ex+C2⋅e5⋅x;y2=−C1⋅ex+3⋅C2⋅e5⋅x.Решить систему ДУ
{dy1dx=3⋅y1−y2dy2dx=4⋅y1−y2.
Систему решаем исключением неизвестной функции y2.
Шаг 1. Из первого уравнения находим y2: y2=−dy1dx+3⋅y1.
Шаг 2. Подставляем y2 во второе уравнение:
dy2dx=4⋅y1+dy1dx−3⋅y1;dy2dx=dy1dx+y1.Шаг 3. Дифференцируем первое уравнение по x: d2y1dx2=3⋅dy1dx−dy2dx.
Шаг 4. Подставляем выражение, полученное на шаге 2, в выражение, полученное на шаге 3:
d2y1dx2=3⋅dy1dx−dy1dx−y1;d2y1dx2−2⋅dy1dx+y1=0.Шаг 5. Решаем линейное однородное ДУ второго порядка с постоянными коэффициентами:
- характеристическое уравнение k2−2⋅k+1=0;
- корни характеристического уравнения k1=1, k2=1 -- действительные, равные;
- искомая функция y1=C1⋅ex+C2⋅x⋅ex.
Шаг 6. Находим функцию y2:
- производная dy1dx=C1⋅ex+C2⋅(ex+x⋅ex);
- результат подстановки в выражение, полученное на шаге 1:
Общее решение данной системы:
y1=C1⋅ex+C2⋅x⋅ex;y2=2⋅C1⋅ex−C2⋅ex+2⋅C2⋅x⋅ex.Решить систему ДУ {dy1dx=y1−3⋅y2dy2dx=3⋅y1+y2.
Систему решаем исключением неизвестной функции y2.
Шаг 1. Из первого уравнения находим y2: y2=13⋅(−dy1dx+y1).
Шаг 2. Подставляем y2 во второе уравнение:
dy2dx=3⋅y1+13⋅(−dy1dx+y1);dy2dx=−13⋅dy1dx+103⋅y1.Шаг 3. Дифференцируем первое уравнение по x: d2y1dx2=dy1dx−3⋅dy2dx.
Шаг 4. Подставляем выражение, полученное на шаге 2, в выражение, полученное на шаге 3:
d2y1dx2=dy1dx−3⋅(−13⋅dy1dx+103⋅y1);d2y1dx2−2⋅dy1dx+10⋅y1=0.Шаг 5. Решаем линейное однородное ДУ второго порядка с постоянными коэффициентами:
- характеристическое уравнение k2−2⋅k+10=0;
- корни характеристического уравнения k1=1+3⋅i, k2=1−3⋅i -- комплексные;
- искомая функция y1=ex⋅(C1⋅cos(3⋅x)+C2⋅sin(3⋅x)).
Шаг 6. Находим функцию y2:
- производная
- результат подстановки в выражение, полученное на шаге 1: y2=13⋅ex⋅(−C1⋅cos(3⋅x)−C2⋅sin(3⋅x))+
dy1dx=ex⋅(C1⋅cos(3⋅x)+C2⋅sin(3⋅x))+ +ex⋅(−3⋅C1⋅sin(3⋅x)+3⋅C2⋅cos(3⋅x));
Общее решение данной системы:
y1=ex⋅(C1⋅cos(3⋅x)+C2⋅sin(3⋅x));