Решатель систем ОДУ

О Решатель систем ОДУ

Решатель систем ОДУ — это универсальный набор инструментов для работы со связанными линейными и нелинейными дифференциальными уравнениями. Вставьте матрицу коэффициентов 2×2 или 3×3, и инструмент выполнит полный анализ собственных значений и векторов, запишет общее и частное решение в LaTeX, классифицирует равновесие в начале координат (седло, узел, спираль или центр) и построит интерактивный фазовый портрет с анимированной траекторией. Для нелинейных планарных систем вы можете ввести произвольные правые части \(f(x,y)\) и \(g(x,y)\), и инструмент создаст высокоточный фазовый портрет методом RK4.

Что такое система ОДУ?

Система обыкновенных дифференциальных уравнений связывает несколько неизвестных функций одной переменной — обычно времени \(t\) — через их производные. В наиболее компактном виде:

Когда \(\mathbf{F}(t, \mathbf{x}) = A\mathbf{x}\) для постоянной матрицы \(A\), система является линейной и автономной. Именно здесь теория наиболее изящна: всё долгосрочное поведение определяется собственными значениями матрицы \(A\).

Алгоритм нахождения собственных значений для линейных систем

Для \(\mathbf{x}' = A\mathbf{x}\) стандартный метод включает:

1. Вычисление характеристического многочлена \(\det(\lambda I - A) = 0\).
2. Нахождение корней — собственных значений \(\lambda_1, \lambda_2, \dots\).
3. Для каждого значения — поиск собственного вектора \(v\) из уравнения \((A - \lambda I) v = 0\).
4. Составление общего решения как линейной комбинации: \(\mathbf{x}(t) = c_1 v_1 e^{\lambda_1 t} + c_2 v_2 e^{\lambda_2 t} + \cdots\).
5. Определение констант \(c_i\) путем подстановки начального условия \(\mathbf{x}(0)\) в общее решение.

Три случая для систем 2×2

Плоскость след-детерминант

Для матрицы 2×2 со следом \(T = a_{11} + a_{22}\) и детерминантом \(D = a_{11} a_{22} - a_{12} a_{21}\) вся классификация умещается в одну диаграмму:

Вот почему в панели результатов на видном месте отображаются \(T\), \(D\) и \(\Delta = T^2 - 4D\) — этих трех чисел достаточно, чтобы назвать тип равновесия.

Нелинейные системы и фазовый портрет

Большинство реальных ОДУ нелинейны и не имеют решения в замкнутой форме. Инструмент обрабатывает их путем численного интегрирования уравнений методом Рунге–Кутты 4-го порядка (RK4), который имеет локальную погрешность \(O(h^5)\) и является стандартным «рабочим инструментом» для гладких векторных полей.

Фазовый портрет включает:

• Векторное поле, построенное на сетке 13×13, показывающее направление потока в каждой точке.
• Траекторию от вашего начального условия, выделенную красным цветом с анимированным оранжевым бегунком, указывающим направление времени.
• Несколько вспомогательных линий потока, дающих глобальную картину динамики.
• Для линейных систем 2×2 — оси собственных векторов (пунктирные голубые линии). Это инвариантные направления, вдоль которых решения движутся экспоненциально.

Как пользоваться этим решателем

1. Выберите режим — Линейная 2×2, Линейная 3×3 или Нелинейная 2D — с помощью вкладок в верхней части формы.
2. Заполните коэффициенты или уравнения. Нажмите на любой «Быстрый пример», чтобы предварительно заполнить каноническую систему (устойчивый узел, центр, седло, маятник, Ван дер Поль и т. д.).
3. Введите начальное условие \((x_0, y_0)\) и интервал времени \(T\). Типичные значения \(T\) составляют 6–20 для осцилляторов и 3–6 для быстро затухающих устойчивых систем.
4. Нажмите «Решить». Появится полная страница результатов с классификацией, собственными значениями, векторами, решением в замкнутой форме (для линейных режимов), анимированным фазовым портретом и графиком временных рядов.
5. Повторите траекторию с помощью кнопки под фазовым портретом, если хотите снова увидеть, как бегунок проходит кривую решения.

Пример — затухающий гармонический осциллятор

Затухающий осциллятор \(\ddot{x} + 2\zeta \omega \dot{x} + \omega^2 x = 0\) можно переписать в виде 2D-системы, приняв \(y = \dot{x}\):

Для \(\omega = 1\) и \(\zeta = 0.2\) (слабое затухание) матрица примет вид \(A = \begin{pmatrix} 0 & 1 \\ -1 & -0.4 \end{pmatrix}\). След \(T = -0.4\), детерминант \(D = 1\), дискриминант \(\Delta = 0.16 - 4 = -3.84 < 0\), таким образом, мы получаем устойчивую спираль с собственными значениями \(-0.2 \pm 0.9798\,i\). Траектория закручивается в начало координат, а временные ряды показывают экспоненциально затухающие синусоиды.

Применение

• Механические системы — связанные системы груз-пружина, маятники, гироскопы.
• Электрические цепи — RLC-контуры, фильтры на ОУ, управление в пространстве состояний.
• Популяционная динамика — модель хищник-жертва Лотки–Вольтерры, конкурирующие виды, эпидемиология (SIR, SIS).
• Химическая кинетика — реакционные сети, осцилляторы Белоусова–Жаботинского.
• Нейробиология — модель нейрона ФитцХью–Нагумо, редукции Ходжкина–Хаксли.
• Теория управления — линеаризованные модели объектов, проектирование наблюдателей, запасы устойчивости.

Советы и нюансы

• Если траектория слишком быстро уходит в бесконечность, уменьшите интервал времени T — неустойчивая система может выйти за пределы видимости за несколько единиц времени.
• Для кратного собственного значения решатель автоматически находит обобщенный собственный вектор \(w\), решая \((A - \lambda I)w = v\), поэтому вы получаете член \(tv\) без ручных вычислений.
• Для нелинейных систем стрелки векторного поля также выявляют положения равновесия, отличные от начала координат, в виде голубых точек — следите за областями с нулевой амплитудой на портрете.
• Для систем 3×3 фазовый портрет отсутствует (3D-графику сложно отобразить на 2D-странице), но расчет решения и вердикт об устойчивости по-прежнему актуальны.
• Начальные условия и интервалы времени не связаны с классификацией: их изменение перемещает только красную траекторию, но не меняет вердикт о собственных значениях.

Часто задаваемые вопросы

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

Система обыкновенных дифференциальных уравнений (ОДУ) — это набор связанных уравнений, описывающих производные нескольких неизвестных функций от одной независимой переменной (обычно времени). Классическая форма: \( \mathbf{x}'(t) = F(t, \mathbf{x}(t)) \), где \( \mathbf{x} \) — вектор состояний, а \(F\) — векторное поле. Линейные системы записываются как \( \mathbf{x}' = A\mathbf{x} + \mathbf{b} \), а их динамика определяется собственными значениями матрицы \(A\).

Как собственные значения классифицируют равновесие системы 2×2?

Для системы \( \mathbf{x}' = A\mathbf{x} \) тип равновесия в нуле определяется следом \(T\) и детерминантом \(D\): \(D < 0\) дает седло; \(D > 0\) при \(T^2 > 4D\) дает узел; \(D > 0\) при \(T^2 < 4D\) дает спираль (или центр при \(T = 0\)). Случай \(T^2 = 4D\) соответствует вырожденному узлу.

Как выглядит решение, если собственные значения комплексные?

При комплексно-сопряженных значениях \( \alpha \pm i\beta \) решение имеет вид \( \mathbf{x}(t) = e^{\alpha t} \left[ c_1 (p \cos\beta t - q \sin\beta t) + c_2 (p \sin\beta t + q \cos\beta t) \right] \). Экспонента определяет рост или затухание, а тригонометрические функции — вращение.

Что происходит при кратном собственном значении?

Если у матрицы кратное значение \(\lambda\) и только один собственный вектор \(v\), используется обобщенный вектор \(w\). Решение принимает вид \( \mathbf{x}(t) = c_1 v e^{\lambda t} + c_2 (tv + w) e^{\lambda t} \). Если же векторов два, решение упрощается до суммы обычных экспонент.

Может ли этот инструмент решать нелинейные системы символьно?

Нелинейный режим использует численный метод Рунге–Кутты (RK4). Большинство таких систем не имеют аналитического решения. Для понимания локального поведения можно линеаризовать систему вблизи точки равновесия (вычислить якобиан) и ввести полученную матрицу в линейный режим 2×2.

Что такое фазовый портрет?

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

Дополнительная литература

• Система дифференциальных уравнений — Википедия
• Фазовый портрет — Википедия
• Собственные числа и векторы — Википедия
• Метод Рунге — Кутты — Википедия
• Осциллятор Ван дер Поля — Википедия