Решение систем из трех дифференциальных уравнений

Системы дифференциальных уравнений с примерами решения и образцами выполнения

Также как и обыкновенные дифференциальные уравнения, системы дифференциальных уравнений применяются для описания многих процессов реальной действительности. В частности, к ним относятся различного рода физические и химические процессы, процессы нефте- и газодобычи, геологии, экономики и т.д. Действительно, если некоторые физические величины (перемещение тела, пластовое давление жидкости в фиксированной точке с тремя координатами, концентрация веществ, объемы продаж продуктов) оказываются меняющимися со временем под воздействием тех или иных факторов, то, как правило, закон их изменения по времени описывается именно системой дифференциальных уравнений, т.е. системой, связывающей исходные переменные как функции времени и производные этих функций. Независимой переменной в системе дифференциальных уравнений может выступать не только время, но и другие физические величины: координата, цена продукта и т.д.

Решение систем из трех дифференциальных уравнений

Видео:Решение системы дифференциальных уравнений методом ЭйлераСкачать

Решение системы дифференциальных уравнений методом Эйлера

Решение систем дифференциальных уравнений

К системе дифференциальных уравнений приводит уже простейшая задача динамики точки: даны силы, действующие на материальную точку; найти закон движения, т. е. найти функции Решение систем из трех дифференциальных уравненийвыражающие зависимость координат движущейся точки от времени. Система, которая при этом получается, в общем случае имеет вид

Решение систем из трех дифференциальных уравнений

Здесь x, у, z — координаты движущейся точки, t — время, f, g, h — известные функции своих аргументов.

Система вида (1) называется канонической. Обращаясь к общему случаю системы т дифференциальных уравнений с т неизвестными функциями Решение систем из трех дифференциальных уравненийаргумента t, назовем канонической систему вида

Решение систем из трех дифференциальных уравнений

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

Решение систем из трех дифференциальных уравнений

Если Решение систем из трех дифференциальных уравненийв (2) принять за новые вспомогательные функции, то общую каноническую систему (2) можно заменить эквивалентной ей нормальной системой, состоящей из Решение систем из трех дифференциальных уравненийуравнений. Поэтому достаточно рассматривать лишь нормальные системы.

Например, одно уравнение

Решение систем из трех дифференциальных уравнений

является мастным случаем канонической системы. Положив Решение систем из трех дифференциальных уравненийв силу исходного уравнения будем иметь

Решение систем из трех дифференциальных уравнений

В результате получаем нормальную систему уравнений

Решение систем из трех дифференциальных уравнений

эквивалентную исходному уравнению.

Определение:

Решением нормальной системы (3) на интервале (а, Ь) изменения аргумента t называется всякая система n функций

Решение систем из трех дифференциальных уравнений

дифференцируемых на интервале а Решение систем из трех дифференциальных уравнений

Теорема:

Существования и единственности решения задачи Коши. Пусть имеем нормальную систему дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

и пусть функции Решение систем из трех дифференциальных уравненийопределены в некоторой (n + 1) — мерной области D изменения переменных Решение систем из трех дифференциальных уравненийЕсли существует окрестность Решение систем из трех дифференциальных уравненийточки Решение систем из трех дифференциальных уравненийв которой функции fi непрерывны по совокупности аргументов и имеют ограниченные частные производные по переменным Решение систем из трех дифференциальных уравненийто найдется интервал Решение систем из трех дифференциальных уравненийизменения t, на котором существует единственное решение нормальной системы (3), удовлетворяющее начальным условиям

Решение систем из трех дифференциальных уравнений

Определение:

Система n функций

Решение систем из трех дифференциальных уравнений

зависящих от t и n произвольных постоянных Решение систем из трех дифференциальных уравненийназывается общим решением нормальной системы (3) в некоторой области Решение систем из трех дифференциальных уравненийсуществования и единственности решения задачи Коши, если

1) при любых допустимых значениях Решение систем из трех дифференциальных уравненийсистема функций (6) обращает уравнения (3) в тождества,

2) в области Решение систем из трех дифференциальных уравненийфункции (6) решают любую задачу Коши.

Решения, получающиеся из общего при конкретных значениях постоянных Решение систем из трех дифференциальных уравненийназываются частными решениями.

Обратимся для наглядности к нормальной системе двух уравнений,

Решение систем из трех дифференциальных уравнений

Будем рассматривать систему значений t, x1, х2 как прямоугольные декартовы координаты точки трехмерного пространства, отнесенного к системе координат Решение систем из трех дифференциальных уравненийРешение

Решение систем из трех дифференциальных уравнений

системы (7), принимающее при Решение систем из трех дифференциальных уравненийзначения Решение систем из трех дифференциальных уравненийопределяет в пространстве некоторую линию, проходящую через точку Решение систем из трех дифференциальных уравненийЭта линия называется интегральной кривой нормальной системы (7). Задача Коши для системы (7) получает следующую геометрическую формулировку: в пространстве переменных t, x1, х2 найти интегральную кривую, проходящую через данную точку Решение систем из трех дифференциальных уравнений(рис. 1). Теорема 1 устанавливает существование и единственность такой кривой.

Решение систем из трех дифференциальных уравнений

Нормальной системе (7) и ее решению можно придать еще такое истолкование: будем независимую переменную t рассматривать как параметр, а решение

Решение систем из трех дифференциальных уравнений

системы — как параметрические уравнения кривой на плоскости Решение систем из трех дифференциальных уравненийЭту плоскость переменных х1х2 называют фазовой плоскостью. В фазовой плоскости решение Решение систем из трех дифференциальных уравненийсистемы (7), принимающее при t = to начальные значения Решение систем из трех дифференциальных уравненийизображается кривой АВ, проходящей через точку Решение систем из трех дифференциальных уравнений(рис. 2). Эту кривую называют траекторией системы (фазовой траекторией). Траектория системы (7) есть проекция интегральной кривой на фазовую плоскость. По интегральной кривой фазовая траектория определяется однозначно, но не наоборот.

Методы интегрирования систем дифференциальных уравнений

Метод исключения

Один из методов интегрирования — метод исключения. Частным случаем канонической системы является одно уравнение n-го порядка, разрешенное относительно старшей производной

Решение систем из трех дифференциальных уравнений

Введя новые функции Решение систем из трех дифференциальных уравненийзаменим это уравнение следующей нормальной системой n уравнений:

Решение систем из трех дифференциальных уравнений

т. е. одно уравнение n-го порядка эквивалентно нормальной системе (1)

Можно утверждать и обратное, что, вообще говоря, нормальная система п уравнений первого порядка эквивалентна одному уравнению порядка n. На этом и основан метод исключения для интегрирования систем дифференциальных уравнений.

Делается это так. Пусть имеем нормальную систему

Решение систем из трех дифференциальных уравнений

Продифференцируем первое из уравнений (2) по t. Имеем

Решение систем из трех дифференциальных уравнений

Заменяя в правой части производные Решение систем из трех дифференциальных уравненийих выражениями Решение систем из трех дифференциальных уравненийполучим

Решение систем из трех дифференциальных уравнений

Уравнение (3) снова дифференцируем по t. Принимая во внимание систему (2), получим

Решение систем из трех дифференциальных уравнений

Продолжая этот процесс, найдем

Решение систем из трех дифференциальных уравнений

Предположим, что определитель

Решение систем из трех дифференциальных уравнений

(якобиан системы функций Решение систем из трех дифференциальных уравненийотличен от нуля при рассматриваемых значениях Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

Тогда система уравнений, составленная из первого уравнения системы (2) и уравнений

Решение систем из трех дифференциальных уравнений

будет разрешима относительно неизвестных Решение систем из трех дифференциальных уравненийПри этом Решение систем из трех дифференциальных уравненийвыразятся через Решение систем из трех дифференциальных уравнений

Внося найденные выражения в уравнение

Решение систем из трех дифференциальных уравнений

получим одно уравнение n-го порядка

Решение систем из трех дифференциальных уравнений

Из самого способа его построения следует, что если Решение систем из трех дифференциальных уравненийесть решения системы (2), то функция х1(t) будет решением уравнения (5).

Обратно, пусть Х1(t) — решение уравнения (5). Дифференцируя это решение по t, вычислим Решение систем из трех дифференциальных уравненийи подставим найденные значения как известные функции

Решение систем из трех дифференциальных уравнений

от t в систему уравнений

Решение систем из трех дифференциальных уравнений

По предположению эту систему можно разрешить относительно Решение систем из трех дифференциальных уравненийт. е найти Решение систем из трех дифференциальных уравненийкак функции от t.

Можно показать, что так построенная система функций

Решение систем из трех дифференциальных уравнений

составляет решение системы дифференциальных уравнений (2). Пример:

Требуется проинтегрировать систему

Решение систем из трех дифференциальных уравнений

Дифференцируя первое уравнение системы, имеем

Решение систем из трех дифференциальных уравнений

откуда, используя второе уравнение, получаем

Решение систем из трех дифференциальных уравнений

— линейное дифференциальное уравнение второго порядка с постоянными коэффициентами с одной неизвестной функцией. Его общее решение имеет вид

Решение систем из трех дифференциальных уравнений

В силу первого уравнения системы находим функцию

Решение систем из трех дифференциальных уравнений

Найденные функции x(t), y(t), как легко проверить, при любых значениях С1 и С2 удовлетворяют заданной системе.

Функции x(t), y(t) можно представить в виде

Решение систем из трех дифференциальных уравнений

откуда видно, что интегральные кривые системы (6) — винтовые линии с шагом Решение систем из трех дифференциальных уравненийи с общей осью х = у = 0, которая также является интегральной кривой (рис. 3).

Решение систем из трех дифференциальных уравнений

Исключая в формулах (7) параметр t, получаем уравнение

Решение систем из трех дифференциальных уравнений

так что фазовые траектории данной системы суть окружности с центром в начале координат — проекции винтовых линий на плоскость хОу.

При А = 0 фазовая траектория состоит из одной точки х = 0, у = 0, называемой точкой покоя системы.

Замечание:

Может оказаться, что функции Решение систем из трех дифференциальных уравненийнельзя выразить через Решение систем из трех дифференциальных уравненийТогда уравнения n-го порядка, эквивалентного исходной системе, мы не получим. Вот простой пример. Систему уравнений

Решение систем из трех дифференциальных уравнений

нельзя заменить эквивалентным уравнением второго порядка относительно х1 или x2. Эта система составлена из пары уравнений 1-го порядка, каждое из которых интегрируется независимо, что дает

Решение систем из трех дифференциальных уравнений

Метод интегрируемых комбинаций

Интегрирование нормальных систем дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

иногда осуществляется методом интегрируемых комбинаций.

Интегрируемой комбинацией называется дифференциальное уравнение, являющееся следствием уравнений (8), но уже легко интегрирующееся.

Пример:

Решение систем из трех дифференциальных уравнений

Складывая почленно данные уравнения, находим одну интегрируемую комбинацию:

Решение систем из трех дифференциальных уравнений

Вычитая почленно из первого уравнения системы второе, получаем вторую интегрируемую комбинацию:

Решение систем из трех дифференциальных уравнений

Мы нашли два конечных уравнения

Решение систем из трех дифференциальных уравнений

из которых легко определяется общее решение системы:

Решение систем из трех дифференциальных уравнений

Одна интегрируемая комбинация дает возможность получить одно уравнение

Решение систем из трех дифференциальных уравнений

связывающее независимую переменную t и неизвестные функции Решение систем из трех дифференциальных уравненийТакое конечное уравнение называется первым интегралом системы (8). Иначе: первым интегралом системы дифференциальных уравнений (8) называется дифференцируемая функция Решение систем из трех дифференциальных уравненийне равная тождественно постоянной, но сохраняющая постоянное значение на любой интегральной кривой этой системы.

Если найдено п первых интегралов системы (8) и все они независимы, т. е. якобиан системы функций Решение систем из трех дифференциальных уравненийотличен от нуля:

Решение систем из трех дифференциальных уравнений

то задача интефирования системы (8) решена (так как из системы

Решение систем из трех дифференциальных уравнений

определяются все неизвестные функции Решение систем из трех дифференциальных уравнений

Системы линейных дифференциальных уравнений

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

Решение систем из трех дифференциальных уравнений

или, в матричной форме,

Решение систем из трех дифференциальных уравнений

Теорема:

Если все функции Решение систем из трех дифференциальных уравненийнепрерывны на отрезке Решение систем из трех дифференциальных уравненийто в достаточно малой окрестности каждой точки Решение систем из трех дифференциальных уравненийгде Решение систем из трех дифференциальных уравненийвыполнены условия теоремы существования и единственности решения задачи Коши, следовательно, через каждую такую точку проходит единственная интегральная кривая системы (1).

Действительно, в таком случае правые части системы (1) непрерывны по совокупности аргументов t, Решение систем из трех дифференциальных уравненийи их частные производные по Решение систем из трех дифференциальных уравненийограничены, так как эти производные равны непрерывным на отрезке [а,b] коэффициентам Решение систем из трех дифференциальных уравнений

Введем линейный оператор

Решение систем из трех дифференциальных уравнений

Тогда система (2) запишется в виде

Решение систем из трех дифференциальных уравнений

Если матрица F — нулевая, т. е. Решение систем из трех дифференциальных уравненийна интервале (а,b), то система (2) называется линейной однородной и имеет вид

Решение систем из трех дифференциальных уравнений

Приведем некоторые теоремы, устанавливающие свойства решений линейных систем.

Теорема:

Если X(t) является решением линейной однородной системы

Решение систем из трех дифференциальных уравнений

то cX(t), где с — произвольная постоянная, является решением той же системы.

Теорема:

Решение систем из трех дифференциальных уравнений

двух решений Решение систем из трех дифференциальных уравненийоднородной линейной системы уравнений является решением той же системы.

Следствие:

Решение систем из трех дифференциальных уравнений

с произвольными постоянными коэффициентами сi решений Решение систем из трех дифференциальных уравненийлинейной однородной системы дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

является решением той же системы.

Теорема:

Если Решение систем из трех дифференциальных уравненийесть решение линейной неоднородной системы

Решение систем из трех дифференциальных уравнений

a Xo(t) — решение соответствующей однородной системы

Решение систем из трех дифференциальных уравнений

будет решением неоднородной системы Решение систем из трех дифференциальных уравнений

Действительно, по условию,

Решение систем из трех дифференциальных уравнений

Пользуясь свойством аддитивности оператора Решение систем из трех дифференциальных уравненийполучаем

Решение систем из трех дифференциальных уравнений

Это означает, что сумма Решение систем из трех дифференциальных уравненийесть решение неоднородной системы уравнений Решение систем из трех дифференциальных уравнений

Определение:

Решение систем из трех дифференциальных уравнений

называются линейно зависимыми на интервале a Решение систем из трех дифференциальных уравнений

при Решение систем из трех дифференциальных уравненийпричем по крайней мере одно из чисел аi, не равно нулю. Если тождество (5) справедливо только при Решение систем из трех дифференциальных уравненийто векторы Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравненийназываются линейно независимыми на (а, b).

Заметим, что одно векторное тождество (5) эквивалентно n тождествам:

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

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

Определение:

Пусть имеем линейную однородную систему

Решение систем из трех дифференциальных уравнений

где Решение систем из трех дифференциальных уравненийматрица с элементами Решение систем из трех дифференциальных уравненийСистема n решений

Решение систем из трех дифференциальных уравнений

линейной однородной системы (6), линейно независимых на интервале а Решение систем из трех дифференциальных уравнений

с непрерывными на отрезке Решение систем из трех дифференциальных уравненийкоэффициентами Решение систем из трех дифференциальных уравненийявляется линейная комбинация п линейно независимых на интервале а Решение систем из трех дифференциальных уравнений

(Решение систем из трех дифференциальных уравнений) — произвольные постоянные числа).

Пример:

Решение систем из трех дифференциальных уравнений

имеет, как нетрудно проверить, решения

Решение систем из трех дифференциальных уравнений

Эти решения линейно независимы, так как определитель Вронского отличен от нуля:

Решение систем из трех дифференциальных уравнений

Общее решение системы имеет вид

Решение систем из трех дифференциальных уравнений

(с1, с2 — произвольные постоянные).

Фундаментальная матрица

Квадратная матрица

Решение систем из трех дифференциальных уравнений

столбцами которой являются линейно независимые решения Решение систем из трех дифференциальных уравненийсистемы (6), называется фундаментальной матрицей этой системы. Нетрудно проверить, что фундаментальная матрица удовлетворяет матричному уравнению

Решение систем из трех дифференциальных уравнений

Если Х(t) — фундаментальная матрица системы (6), то общее решение системы можно представить в виде

Решение систем из трех дифференциальных уравнений

— постоянная матрица-столбец с произвольными элементами. Полагая в (7) t = t0, имеем

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

Матрица Решение систем из трех дифференциальных уравненийназывается матрицей Коши. С ее помощью решение системы (6) можно представить так:

Решение систем из трех дифференциальных уравнений

Теорема:

О структуре общего решения линейной неоднородной системы дифференциальных уравнений. Общее решение в области Решение систем из трех дифференциальных уравненийлинейной неоднородной системы дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

с непрерывными на отрезке Решение систем из трех дифференциальных уравненийкоэффициентами aij(t) и правыми частями fi(t) равно сумме общего решения

Решение систем из трех дифференциальных уравнений

соответствующей однородной системы и какого-нибудь частного решения Решение систем из трех дифференциальных уравненийнеоднородной системы (2):

Решение систем из трех дифференциальных уравнений

Метод вариации постоянных

Если известно общее решение линейной однородной системы (6), то частное решение неоднородной системы можно находить методом вариации постоянных (метод Лагранжа).

Решение систем из трех дифференциальных уравнений

есть общее решение однородной системы (6), тогда

Решение систем из трех дифференциальных уравнений

причем решения Xk(t) линейно независимы.

Будем искать частное решение неоднородной системы

Решение систем из трех дифференциальных уравнений

где Решение систем из трех дифференциальных уравненийнеизвестные функции от t. Дифференцируя Решение систем из трех дифференциальных уравненийпо t, имеем

Решение систем из трех дифференциальных уравнений

Подставляя Решение систем из трех дифференциальных уравненийв (2), получаем

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

то для определения Решение систем из трех дифференциальных уравненийполучаем систему

Решение систем из трех дифференциальных уравнений

или, в развернутом виде,

Решение систем из трех дифференциальных уравнений

Система (10) есть линейная алгебраическая система относительно Решение систем из трех дифференциальных уравненийопределителем которой является определитель Вронского W(t) фундаментальной системы решений Решение систем из трех дифференциальных уравнений. Этот определитель отличен от нуля всюду на интервале a Решение систем из трех дифференциальных уравнений

где Решение систем из трех дифференциальных уравнений— известные непрерывные функции. Интегрируя последние соотношения, находим

Решение систем из трех дифференциальных уравнений

Подставляя эти значения Решение систем из трех дифференциальных уравненийв (9), находим частное решение системы (2)

Решение систем из трех дифференциальных уравнений

(здесь под символом Решение систем из трех дифференциальных уравненийпонимается одна из первообразных для функции Решение систем из трех дифференциальных уравнений

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

Рассмотрим линейную систему дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

в которой все коэффициенты Решение систем из трех дифференциальных уравнений— постоянные. Чаще всего такая система интегрируется сведением ее к одному уравнению более высокого порядка, причем это уравнение будет также линейным с постоянными коэффициентами. Другой эффективный метод интегрирования систем с постоянными коэффициентами — метод преобразования Лапласа.

Мы рассмотрим еще метод Эйлера интегрирования линейных однородных систем дифференциальных уравнений с постоянными коэффициентами. Он состоит в следующем.

Метод Эйлера

Будем искать решение системы

Решение систем из трех дифференциальных уравнений

где Решение систем из трех дифференциальных уравнений— постоянные. Подставляя Xk в форме (2) в систему (1), сокращая на Решение систем из трех дифференциальных уравненийи перенося все члены в одну часть равенства, получаем систему

Решение систем из трех дифференциальных уравнений

Для того, чтобы эта система (3) линейных однородных алгебраических уравнений с n неизвестными Решение систем из трех дифференциальных уравненийимела нетривиальное решение, необходимо и достаточно, чтобы ее определитель был равен нулю:

Решение систем из трех дифференциальных уравнений

Уравнение (4) называется характеристическим. В его левой части стоит многочлен относительно Решение систем из трех дифференциальных уравненийстепени n. Из этого уравнения определяются те значения Решение систем из трех дифференциальных уравнений, при которых система (3) имеет нетривиальные решения Решение систем из трех дифференциальных уравнений. Если все корни Решение систем из трех дифференциальных уравненийхарактеристического уравнения (4) различны, то, подставляя их по очереди в систему (3), находим соответствующие им нетривиальные решения Решение систем из трех дифференциальных уравненийэтой системы n, следовательно, находим п решений исходной системы дифференциальных уравнений (1) в виде

Решение систем из трех дифференциальных уравнений

где второй индекс указывает номер решения, а первый — номер неизвестной функции. Построенные таким образом п частных решений линейной однородной системы (1)

Решение систем из трех дифференциальных уравнений

образуют, как можно проверить, фундаментальную систему решений этой системы.

Следовательно, общее решение однородной системы дифференциальных уравнений (1) имеет вид

Решение систем из трех дифференциальных уравнений

где Решение систем из трех дифференциальных уравненийпроизвольные постоянные.

Случай, когда характеристическое уравнение имеет кратные корни, мы рассматривать не будем.

Пример:

Решение систем из трех дифференциальных уравнений

Ищем решение в виде

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

имеет корни Решение систем из трех дифференциальных уравнений

Система (3) для определения a1, а2 выглядит так:

Решение систем из трех дифференциальных уравнений

Подставляя в (*) Решение систем из трех дифференциальных уравненийполучаем

Решение систем из трех дифференциальных уравнений

откуда а21 = а11. Следовательно,

Решение систем из трех дифференциальных уравнений

Полагая в Решение систем из трех дифференциальных уравненийнаходим a22 = — a12, поэтому

Решение систем из трех дифференциальных уравнений

Общее решение данной системы:

Решение систем из трех дифференциальных уравнений

Матричный метод

Изложим еще матричный метод интегрирования однородной системы (1). Запишем систему (1) в виде

Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравненийматрица с постоянными действительными элементами Решение систем из трех дифференциальных уравнений

Напомним некоторые понятия из линейной алгебры. Вектор Решение систем из трех дифференциальных уравненийназывается собственным вектором матрицы А, если

Решение систем из трех дифференциальных уравнений

Число Решение систем из трех дифференциальных уравненийназывается собственным значением матрицы А, отвечающим собственному вектору g, и является корнем характеристического уравнения

Решение систем из трех дифференциальных уравнений

где I — единичная матрица.

Будем предполагать, что все собственные значения Решение систем из трех дифференциальных уравненийматрицы А различны. В этом случае собственные векторы g1, g2, …gn линейно независимы и существует Решение систем из трех дифференциальных уравненийматрица Т, приводящая матрицу А к диагональному виду, т. е. такая, что

Решение систем из трех дифференциальных уравнений

Столбцами матрицы Т являются координаты собственных векторов g1, g2 …, gn матрицы А.

Введем еще следующие понятия. Пусть В(t) — Решение систем из трех дифференциальных уравненийматрица, элементы Решение систем из трех дифференциальных уравненийкоторой суть функции аргумента t, определенные на множестве Решение систем из трех дифференциальных уравнений. Матрица В(t) называется непрерывной на Решение систем из трех дифференциальных уравнений, если непрерывны на Решение систем из трех дифференциальных уравненийвсе ее элементы Решение систем из трех дифференциальных уравнений. Матрица В(t) называется дифференцируемой на Решение систем из трех дифференциальных уравнений, если дифференцируемы на Решение систем из трех дифференциальных уравненийвсе элементы Решение систем из трех дифференциальных уравненийэтой матрицы. При этом производной матрицы Решение систем из трех дифференциальных уравненийназывается матрица, элементами которой являются производные Решение систем из трех дифференциальных уравненийу соответствующих элементов матрицы В(t).

Пусть B(t) — n х n-матрица,

Решение систем из трех дифференциальных уравнений

— вектор-столбец. Учитывая правила алгебры матриц, непосредственной проверкой убеждаемся в справедливости формулы

Решение систем из трех дифференциальных уравнений

В частности, если В — постоянная матрица, то

Решение систем из трех дифференциальных уравнений

так как Решение систем из трех дифференциальных уравненийесть нуль-матрица.

Теорема:

Если собственные значения Решение систем из трех дифференциальных уравненийматрицы А различны, то общее решение системы (7) имеет вид

Решение систем из трех дифференциальных уравнений

где g1, g2,…, gn — собственные векторы-столбцы матрицы А, Решение систем из трех дифференциальных уравненийпроизвольные постоянные числа.

Введем новый неизвестный вектор-столбец Y(t) по формуле

Решение систем из трех дифференциальных уравнений

где Т — матрица, приводящая матрицу А к диагональному виду. Подставляя X(t) из (11) в (7), получим систему

Решение систем из трех дифференциальных уравнений

Умножая обе части последнего соотношения слева на Решение систем из трех дифференциальных уравненийи учитывая, что Решение систем из трех дифференциальных уравненийпридем к системе

Решение систем из трех дифференциальных уравнений

Мы получили систему из n независимых уравнений, которая без труда интегрируется:

Решение систем из трех дифференциальных уравнений

Здесь Решение систем из трех дифференциальных уравнений— произвольные постоянные числа.

Вводя единичные n-мерные векторы-столбцы

Решение систем из трех дифференциальных уравнений

решение Y(t) можно представить в виде

Решение систем из трех дифференциальных уравнений

В силу (11) Х(t) = TY(t). Так как столбцы матрицы Т есть собственные векторы матрицы Решение систем из трех дифференциальных уравненийсобственный вектор матрицы А. Поэтому, подставляя (13) в (11), получим формулу (10):

Решение систем из трех дифференциальных уравнений

Таким образом, если матрица А системы дифференциальных уравнений (7) имеет различные собственные значения, для получения общего решения этой системы:

1) находим собственные значения Решение систем из трех дифференциальных уравненийматрицы как корни алгебраического уравнения

Решение систем из трех дифференциальных уравнений

2) находим все собственные векторы g1, g2,…, gn;

3) выписываем общее решение системы дифференциальных уравнений (7) по формуле (10).

Пример:

Решение систем из трех дифференциальных уравнений

Матрица А системы имеет вид

Решение систем из трех дифференциальных уравнений

1) Составляем характеристическое уравнение

Решение систем из трех дифференциальных уравнений

Корни характеристического уравнения Решение систем из трех дифференциальных уравнений

2) Находим собственные векторы

Решение систем из трех дифференциальных уравнений

Для Решение систем из трех дифференциальных уравнений= 4 получаем систему

Решение систем из трех дифференциальных уравнений

откуда g11 = g12, так что

Решение систем из трех дифференциальных уравнений

Аналогично для Решение систем из трех дифференциальных уравнений= 1 находим

Решение систем из трех дифференциальных уравнений

3) Пользуясь формулой (10), получаем общее решение системы дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

Корни характеристического уравнения могут быть действительными и комплексными. Так как по предположению коэффициенты Решение систем из трех дифференциальных уравненийсистемы (7) действительные, то характеристическое уравнение

Решение систем из трех дифференциальных уравнений

будет иметь действительные коэффициенты. Поэтому наряду с комплексным корнем Решение систем из трех дифференциальных уравненийоно будет иметь и корень Решение систем из трех дифференциальных уравнений*, комплексно сопряженный с Решение систем из трех дифференциальных уравнений. Нетрудно показать, что если g — собственный вектор, отвечающий собственному значению Решение систем из трех дифференциальных уравнений, то Решение систем из трех дифференциальных уравнений* — тоже собственное значение, которому отвечает собственный вектор g*, комплексно сопряженный с g.

При комплексном Решение систем из трех дифференциальных уравненийрешение

Решение систем из трех дифференциальных уравнений

системы (7) также будет комплексным. Действительная часть

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

этого решения являются решениями системы (7). Собственному значению Решение систем из трех дифференциальных уравнений* будет отвечать пара действительных решений X1 и -Х2, т. е. та же пара, что и для собственного значения Решение систем из трех дифференциальных уравнений. Таким образом, паре Решение систем из трех дифференциальных уравнений, Решение систем из трех дифференциальных уравнений* комплексно сопряженных собственных значений отвечает пара действительных решений системы (7) дифференциальных уравнений.

Пусть Решение систем из трех дифференциальных уравнений— действительные собственные значения, Решение систем из трех дифференциальных уравненийРешение систем из трех дифференциальных уравнений— комплексные собственные значения. Тогда всякое действительное решение системы (7) имеет вид

Решение систем из трех дифференциальных уравнений

где сi — произвольные постоянные.

Пример:

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

1) Характеристическое уравнение системы

Решение систем из трех дифференциальных уравнений

Его корни Решение систем из трех дифференциальных уравнений

2) Собственные векторы матриц

Решение систем из трех дифференциальных уравнений

3) Решение системы

Решение систем из трех дифференциальных уравнений

где а1, а2 — произвольные комплексные постоянные.

Найдем действительные решения системы. Пользуясь формулой Эйлера

Решение систем из трех дифференциальных уравнений

Следовательно, всякое действительное решение системы имеет

Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений

где с1, с2 — произвольные действительные числа.

Видео:15. Однородная система линейных уравнений / фундаментальная система решенийСкачать

15. Однородная система линейных уравнений / фундаментальная система решений

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

Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений

Решение заданий и задач по предметам:

Дополнительные лекции по высшей математике:

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений Решение систем из трех дифференциальных уравнений

Образовательный сайт для студентов и школьников

Копирование материалов сайта возможно только с указанием активной ссылки «www.lfirmal.com» в качестве источника.

© Фирмаль Людмила Анатольевна — официальный сайт преподавателя математического факультета Дальневосточного государственного физико-технического института

Видео:Видеоурок "Системы дифференциальных уравнений"Скачать

Видеоурок "Системы дифференциальных уравнений"

Системы дифференциальных уравнений

Этот раздел мы решили посвятить решению систем дифференциальных уравнений простейшего вида d x d t = a 1 x + b 1 y + c 1 d y d t = a 2 x + b 2 y + c 2 , в которых a 1 , b 1 , c 1 , a 2 , b 2 , c 2 — некоторые действительные числа. Наиболее эффективным для решения таких систем уравнений является метод интегрирования. Также рассмотрим решение примера по теме.

Решением системы дифференциальных уравнений будет являться пара функций x ( t ) и y ( t ) , которая способна обратить в тождество оба уравнения системы.

Рассмотрим метод интегрирования системы ДУ d x d t = a 1 x + b 1 y + c 1 d y d t = a 2 x + b 2 y + c 2 . Выразим х из 2 -го уравнения системы для того, чтобы исключить неизвестную функцию x ( t ) из 1 -го уравнения:

d y d t = a 2 x + b 2 y + c 2 ⇒ x = 1 a 2 d y d t — b 2 y — c 2

Выполним дифференцирование 2 -го уравнения по t и разрешим его уравнение относительно d x d t :

d 2 y d t 2 = a 2 d x d t + b 2 d y d t ⇒ d x d t = 1 a 2 d 2 y d t 2 — b 2 d y d t

Теперь подставим результат предыдущих вычислений в 1 -е уравнение системы:

d x d t = a 1 x + b 1 y + c 1 ⇒ 1 a 2 d 2 y d t 2 — b 2 d y d t = a 1 a 2 d y d t — b 2 y — c 2 + b 1 y + c 1 ⇔ d 2 y d t 2 — ( a 1 + b 2 ) · d y d t + ( a 1 · b 2 — a 2 · b 1 ) · y = a 2 · c 1 — a 1 · c 2

Так мы исключили неизвестную функцию x ( t ) и получили линейное неоднородное ДУ 2 -го порядка с постоянными коэффициентами. Найдем решение этого уравнения y ( t ) и подставим его во 2 -е уравнение системы. Найдем x ( t ) . Будем считать, что на этом решение системы уравнений будет закончено.

Найдите решение системы дифференциальных уравнений d x d t = x — 1 d y d t = x + 2 y — 3

Начнем с первого уравнения системы. Разрешим его относительно x :

x = d y d t — 2 y + 3

Теперь выполним дифференцирование 2 -го уравнения системы, после чего разрешим его относительно d x d t : d 2 y d t 2 = d x d t + 2 d y d t ⇒ d x d t = d 2 y d t 2 — 2 d y d t

Полученный в ходе вычислений результат мы можем подставить в 1 -е уравнение системы ДУ:

d x d t = x — 1 d 2 y d t 2 — 2 d y d t = d y d t — 2 y + 3 — 1 d 2 y d t 2 — 3 d y d t + 2 y = 2

В результате преобразований мы получили линейное неоднородное дифференциальное уравнение 2 -го порядка с постоянными коэффициентами d 2 y d t 2 — 3 d y d t + 2 y = 2 . Если мы найдем его общее решение, то получим функцию y ( t ) .

Общее решение соответствующего ЛОДУ y 0 мы можем найти путем вычислений корней характеристического уравнения k 2 — 3 k + 2 = 0 :

D = 3 2 — 4 · 2 = 1 k 1 = 3 — 1 2 = 1 k 2 = 3 + 1 2 = 2

Корни, которые мы получили, являются действительными и различными. В связи с этим общее решение ЛОДУ будет иметь вид y 0 = C 1 · e t + C 2 · e 2 t .

Теперь найдем частное решение линейного неоднородного ДУ y

d 2 y d t 2 — 3 d y d t + 2 y = 2

Правая часть записи уравнения представляет собой многочлен нулевой степени. Это значит, что частное решение будем искать в виде y

= A , где А – это неопределенный коэффициент.

Определить неопределенный коэффициент мы можем из равенства d 2 y

= 2 :
d 2 ( A ) d t 2 — 3 d ( A ) d t + 2 A = 2 ⇒ 2 A = 2 ⇒ A = 1

Таким образом, y

= 1 и y ( t ) = y 0 + y

= C 1 · e t + C 2 · e 2 t + 1 . Одну неизвестную функцию мы нашли.

Теперь подставим найденную функцию во 2 -е уравнение системы ДУ и разрешим новое уравнение относительно x ( t ) :
d ( C 1 · e t + C 2 · e 2 t + 1 ) d t = x + 2 · ( C 1 · e t + C 2 · e 2 t + 1 ) — 3 C 1 · e t + 2 C 2 · e 2 t = x + 2 C 1 · e t + 2 C 2 · e 2 t — 1 x = — C 1 · e t + 1

Так мы вычислили вторую неизвестную функцию x ( t ) = — C 1 · e t + 1 .

Ответ: x ( t ) = — C 1 · e t + 1 y ( t ) = C 1 · e t + C 2 · e 2 t + 1

Видео:Решение систем Д/У: 1. Знакомство с функциями odeXYСкачать

Решение систем Д/У: 1. Знакомство с функциями odeXY

Системы дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

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

Для этого необходимо располагать, вообще говоря, таким же числом уравнений. Если каждое из этих уравнений является дифференциальным, то есть имеет вид соотношения, связывающего неизвестные функции и их производные, то говорят о системе дифференциальных уравнений.

1. Нормальная система дифференциальных уравнений первого порядка. Задача Коши.

Определение. Системой дифференциальных уравнений называется совокупность уравнений, содержащих несколько неизвестных функций и их производные, причём в каждое из уравнений входит хотя бы одна производная.

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

Линейная система называется нормальной, если она разрешена относительно всех производных

Решение систем из трех дифференциальных уравнений(1)

В нормальной системе правые части уравнений не содержат производных искомых функций.

Решением системы дифференциальных уравнений называется совокупность функций Решение систем из трех дифференциальных уравненийудовлетворяющих каждому из уравнений этой системы.

Равенства при Решение систем из трех дифференциальных уравненийназываются начальными условиями системы дифференциальных уравнений.

Часто начальные условия записывают в виде

Решение систем из трех дифференциальных уравнений

Общим решением (интегралом) системы дифференциальных уравнений называется совокупность «n» функций от независимой переменной x и «n» произвольных постоянных C1 , C2 , …,Cn:

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений(2)

Решение систем из трех дифференциальных уравнений

которые удовлетворяют всем уравнениям этой системы.

Чтобы получить частное решение системы, удовлетворяющее заданным начальным условиям Решение систем из трех дифференциальных уравнений, надо из уравнений (2) определить соответствующие начальным условиям значения постоянных C10 , C20 , …,Cn0 .

Задача Коши для системы дифференциальных уравнений состоит в том, чтобы найти такое решение, которое при Решение систем из трех дифференциальных уравненийпринимало бы заданные значения Решение систем из трех дифференциальных уравнений.

Записывается задача Коши для нормальной системы дифференциальных уравнений следующим образом

Решение систем из трех дифференциальных уравнений

Теорема существования и единственности решения задачи Коши.

Для нормальной системы дифференциальных уравнений (1) теорема Коши существования и единственности решения формулируется следующим образом:

Теорема. Пусть правые части уравнений системы (1), т. е. функции Решение систем из трех дифференциальных уравнений, (i=1,2,…,n) непрерывны по всем переменным в некоторой области D и имеет в ней непрерывные частные производные Решение систем из трех дифференциальных уравнений.

Тогда каковы бы ни были значения Решение систем из трех дифференциальных уравнений, принадлежащие области D, существует единственное решение системы (1) Решение систем из трех дифференциальных уравнений, удовлетворяющее начальным условиям

Решение систем из трех дифференциальных уравнений.

2. Решение нормальной системы методом исключения.

Для решения нормальной системы дифференциальных уравнений используется метод исключения неизвестных или метод Коши.

Пусть дана нормальная система

Решение систем из трех дифференциальных уравнений

Дифференцируем по х первое уравнение системы

Решение систем из трех дифференциальных уравнений

Заменяя производные Решение систем из трех дифференциальных уравненийих выражениями Решение систем из трех дифференциальных уравненийиз системы уравнений (1), будем иметь

Решение систем из трех дифференциальных уравнений

Дифференцируем полученное уравнение и поступая аналогично предыдущему, найдём

Решение систем из трех дифференциальных уравнений

Продолжая далее таким же образом, получим уравнение

Решение систем из трех дифференциальных уравнений

Итак, получили систему

Решение систем из трех дифференциальных уравнений(2)

Из первых п-1 уравнений определим y2 , y3 , … , yn , выразив их через

Решение систем из трех дифференциальных уравненийи Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений(3)

Подставляя эти выражения в последнее из уравнений (2), получим уравнения п-го порядка для определения y1 :

Решение систем из трех дифференциальных уравнений(4)

Решив это уравнение, найдём y1

Решение систем из трех дифференциальных уравнений(5)

Дифференцируя последнее выражение п-1 раз, найдём производные

Решение систем из трех дифференциальных уравнений

как функции от Решение систем из трех дифференциальных уравнений. Подставляя эти функции в уравнения (4), определим y2 , y3 , … , yn .

Решение систем из трех дифференциальных уравнений

Итак, получили общее решение системы (1)

Решение систем из трех дифференциальных уравнений(6)

Чтобы найти частное решение системы (1) удовлетворяющее начальным условиям при Решение систем из трех дифференциальных уравнений

надо найти из уравнения (6) соответствующие значения произвольных постоянных С1 , С2 , … , Сn .

Найти общее решение системы уравнений:

Решение систем из трех дифференциальных уравнений

Продифференцируем первое уравнение: Решение систем из трех дифференциальных уравнений

Подставим в это выражение производную у¢ =2x + 2y из второго уравнения.

Решение систем из трех дифференциальных уравнений

Подставим сюда у, выраженное из первого уравнения:

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

Решение систем из трех дифференциальных уравнений

получаем решение системы: Решение систем из трех дифференциальных уравнений

3. Преобразование дифференциального уравнения порядка п к нормальной системе Коши.

Всякое уравнение п-го порядка

Решение систем из трех дифференциальных уравнений

можно привести к системе уравнений первого порядка, если принять Решение систем из трех дифференциальных уравнений

за новые неизвестные функции.

С системами дифференциальных уравнений встречаются при изучении процессов, для описания которых одной функции недостаточно. Например, отыскание векторных линий поля требует реше­ния системы дифференциальных уравнений. Решение задач динамики криволинейного движения при­водит к системе трех дифференциальных уравнений, в которых неиз­вестными функциями являются проекции движущейся точки на оси координат, а независимой переменной — время. Позже вы узнаете, что решение задач электротехники для двух электрических цепей, нахо­дящихся в электромагнитной связи, потребует решения системы двух дифференциальных уравнений. Количество подобных примеров легко можно увеличить.

🔍 Видео

Системы дифференциальных уравнений. Часть 2Скачать

Системы дифференциальных уравнений. Часть 2

Системы дифференциальных уравненийСкачать

Системы дифференциальных уравнений

Метод Крамера за 3 минуты. Решение системы линейных уравнений - bezbotvyСкачать

Метод Крамера за 3 минуты. Решение системы линейных уравнений - bezbotvy

Система с тремя переменнымиСкачать

Система с тремя переменными

Матричный метод решения систем уравненийСкачать

Матричный метод решения систем уравнений

Решение системы уравнений методом ГауссаСкачать

Решение системы уравнений методом Гаусса

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

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

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

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

Лабораторная работа 1. Решение систем обыкновенных дифференциальных уравненийСкачать

Лабораторная работа 1. Решение систем обыкновенных дифференциальных уравнений

ДУ Линейные системыСкачать

ДУ Линейные системы

18+ Математика без Ху!ни. Дифференциальные уравнения.Скачать

18+ Математика без Ху!ни. Дифференциальные уравнения.

Численное решение системы дифференциальных уравнений(задачи Коши)Скачать

Численное решение системы дифференциальных уравнений(задачи Коши)

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

Система неоднородных дифференциальных уравнений
Поделиться или сохранить к себе: