Решатель ОДУ второго порядка

О Решатель ОДУ второго порядка

Решатель ОДУ второго порядка принимает линейное обыкновенное дифференциальное уравнение вида a·y″ + b·y′ + c·y = g(x) с постоянными вещественными коэффициентами, автоматически выводит его характеристическое уравнение, классифицирует режим демпфирования (с передемпфированием, критическое, с недодемпфированием, без демпфирования или неустойчивое) и выдает как аналитическое решение в замкнутой форме, так и высокоточное численное решение. Интерактивный вывод сочетает временной график с двумя кривыми y(x) и y′(x) и траекторию на фазовой плоскости (y, y′) — вид, который позволяет мгновенно определить режим: сходящаяся спираль при недодемпфировании, узел при передемпфировании, замкнутая петля при отсутствии демпфирования, расходящаяся спираль при неустойчивости.

Что такое линейное ОДУ второго порядка с постоянными коэффициентами?

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

где a ≠ 0, b, c — вещественные константы, а g(x) — вынуждающий член. Два начальных условия y(x₀) = y₀ и y′(x₀) = y′₀ превращают его в задачу Коши (начальную задачу) с единственным решением в окрестности x₀ — это следует из теоремы Пикара-Линделёфа, примененной к эквивалентной системе первого порядка.

Если g(x) = 0, уравнение называется однородным. В противном случае оно является неоднородным, и полное решение представляется в виде суммы

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

Характеристическое уравнение

Поиск решения в виде y = e^(r·x) для однородного уравнения приводит к характеристическому (или вспомогательному) уравнению

Это квадратное уравнение, чей дискриминант Δ = b² − 4ac определяет всё качественное поведение системы:

Три случая корней и режимы демпфирования

Метод неопределенных коэффициентов (неоднородный случай)

Когда g(x) имеет одну из следующих простых форм, метод неопределенных коэффициентов позволяет найти частное решение, предполагая вид решения аналогичным g с неизвестными коэффициентами:

• Константа g(x) = k. Пробная функция: y_p = K. Если c = 0, умножьте на x; если также b = 0, умножьте на x еще раз.
• Многочлен степени n. Пробная функция: общий многочлен степени n. Умножьте на x или x², если константа или линейный член резонируют с корнями.
• Экспонента g(x) = A·e^(k·x). Пробная функция: y_p = K·e^(k·x). Если k совпадает с характеристическим корнем, умножьте на x (простой корень) или x² (кратный корень) — это резонанс.
• Синусоида g(x) = A·cos(ω·x) + B·sin(ω·x). Пробная функция: y_p = K₁·cos(ω·x) + K₂·sin(ω·x). Умножьте на x, если iω является корнем (резонанс на чистой частоте).
• Произведения и суммы обрабатываются согласно принципу линейности и правилу произведения.

Анализ фазовой плоскости

Эквивалентная система первого порядка имеет вид u = y, v = y′, где u′ = v и v′ = (g(x) − b·v − c·u)/a. Параметрический график v от u по переменной x дает траекторию на фазовой плоскости. Для однородных автономных систем (где g не зависит от x) орбиты однозначно определяются начальной точкой (y₀, y′₀) и наглядно показывают режим:

• Недодемпфирование: траектория закручивается в спираль к началу координат.
• Передемпфирование: траектория приближается к началу координат вдоль инвариантной линии (медленного собственного вектора).
• Критическое демпфирование: вырожденный узел, траектория касается единственного собственного вектора.
• Без демпфирования: замкнутый эллипс вокруг начала координат — вечные колебания.
• Неустойчивость: траектория закручивается наружу или уходит в бесконечность.

Пример решения: Вынужденный затухающий гармонический осциллятор

Рассмотрим уравнение y″ + 2·y′ + 5·y = 10 с начальными условиями y(0) = 0, y′(0) = 0 — вынужденная система с недодемпфированием.

1. Характеристическое уравнение: r² + 2r + 5 = 0 → Δ = 4 − 20 = −16 → r = −1 ± 2i.
2. Однородное решение: y_h = e^(−x)·(C₁·cos 2x + C₂·sin 2x).
3. Частное решение для постоянного воздействия g = 10: ищем y_p = K, тогда 5K = 10, откуда y_p = 2.
4. Применяем НУ: y(0) = 0 → C₁ + 2 = 0 → C₁ = −2. y′(0) = 0 → −C₁ + 2C₂ = 0 → C₂ = −1.
5. Итоговый ответ: y(x) = 2 − e^(−x)·(2·cos 2x + sin 2x) — колебания с затухающей огибающей и пределом y → 2.

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

1. Введите коэффициенты a, b, c в верхней строке. a не должно быть равно нулю (иначе уравнение станет первым порядком).
2. Введите вынуждающий член g(x) или оставьте 0 для однородной задачи. Аналитические частные решения выводятся для констант, многочленов до 2-й степени и одиночных экспонент A·e^(k·x), включая случай резонанса.
3. Укажите начальные условия (x₀, y₀, y′₀). Необходимо задать и y, и y′ в точке x₀, так как это уравнение второго порядка.
4. Выберите диапазон x для графиков. Решатель выполняет интегрирование от x₀ в обоих направлениях x с помощью RK4.
5. Нажмите "Решить и визуализировать". Вы получите характеристическое уравнение с его корнями на комплексной плоскости, классификацию режима демпфирования, аналитические однородное и частное решения, временной график y и y′, а также траекторию на фазовой плоскости.

Типичные области применения

• Механические системы "пружина-масса-демпфер": m·x″ + c·x′ + k·x = F(t). Передемпфированный, критический и недодемпфированный режимы соответствуют различным коэффициентам затухания ζ = c/(2·√(m·k)).
• Электрические RLC-цепи: последовательные RLC-цепи описываются уравнением L·Q″ + R·Q′ + Q/C = V(t) — идентичная структура, другие символы.
• Маятник (малые углы): θ″ + (g/L)·θ = 0 дает гармонические колебания; учет сопротивления воздуха дает затухающие колебания.
• Реакция зданий на землетрясения: конструкция с одной степенью свободы, где ускорение основания выступает в качестве вынуждающей силы.
• Сервосистемы с ПИД-регулированием: динамика ошибки в замкнутом контуре сводится к ОДУ второго порядка, где коэффициент демпфирования определяет перерегулирование.
• Популяционные модели с инерцией: экономический рост с лагом накопления капитала или экологические модели с задержкой реакции.

Численный метод — классический метод Рунге-Кутты (RK4) для 2D-системы

Инструмент сводит уравнение a·y″ + b·y′ + c·y = g(x) к системе первого порядка

с начальными условиями u(x₀) = y₀, v(x₀) = y′₀. Четырехстадийный метод Рунге-Кутты применяется к вектору состояния (u, v). Метод RK4 имеет локальную погрешность O(h⁵) и глобальную погрешность O(h⁴); стандартные 400 подшагов в каждом направлении обеспечивают точность около шести значащих цифр для нежестких задач.

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

Что такое линейное ОДУ второго порядка с постоянными коэффициентами?

Линейное ОДУ второго порядка с постоянными коэффициентами имеет вид a·y″ + b·y′ + c·y = g(x), где a, b, c — вещественные константы, а g(x) — вынуждающий (неоднородный) член. При двух начальных условиях y(x₀) = y₀ и y′(x₀) = y′₀ решение единственно. Однородный случай g(x) = 0 всегда имеет аналитическое решение через характеристическое уравнение a·r² + b·r + c = 0; неоднородный случай решается как y(x) = y_h(x) + y_p(x).

Что такое характеристическое уравнение?

Для уравнения a·y″ + b·y′ + c·y = 0 подстановка вида y = e^(r·x) дает a·r² + b·r + c = 0 — это характеристическое (вспомогательное) уравнение. Его корни определяют вид однородного решения: два различных вещественных корня дают y_h = C₁·e^(r₁·x) + C₂·e^(r₂·x); один кратный корень r дает y_h = (C₁ + C₂·x)·e^(r·x); комплексно-сопряженные корни α ± β·i дают y_h = e^(α·x)·(C₁·cos(β·x) + C₂·sin(β·x)).

Что означают термины недодемпфированный, критический и передемпфированный?

Эти термины взяты из модели пружина-масса-демпфер m·x″ + c·x′ + k·x = 0. Передемпфирование (дискриминант > 0, два вещественных корня) означает, что система медленно возвращается к равновесию без колебаний. Критическое демпфирование (дискриминант = 0, кратный корень) — это самый быстрый возврат без "перелета" через ноль. Недодемпфирование (дискриминант < 0, комплексные корни) дает затухающие колебания. Без демпфирования (b = 0, c/a > 0) система совершает чистые синусоидальные колебания бесконечно долго.

Что такое метод неопределенных коэффициентов?

Для простых функций g(x) — констант, многочленов, экспонент, синусов, косинусов и их произведений — частное решение y_p ищется в том же виде, что и g, с неизвестными коэффициентами, которые находятся подстановкой в ОДУ. Если g(x) резонирует с характеристическим корнем, пробную функцию нужно умножить на x (или x² для кратных корней).

Что такое фазовая плоскость?

Для уравнения второго порядка, сведенного к 2D-системе (y, y′), фазовая плоскость отображает зависимость y′ от y по мере изменения x. Кривые на этой плоскости наглядно показывают режим: сходящиеся спирали для недодемпфирования, входящие узлы для передемпфирования, замкнутые эллипсы для гармонических колебаний и расходящиеся спирали для неустойчивости. Это геометрический аналог диаграммы корней характеристического уравнения.

Какой численный метод использует этот инструмент?

Применяется классический метод Рунге-Кутты четвертого порядка (RK4) для эквивалентной системы первого порядка u = y, v = y′, где u′ = v и v′ = (g(x) − b·v − c·u)/a. RK4 обладает локальной погрешностью аппроксимации O(h⁵), а 400 подшагов на каждом интервале обеспечивают точность около шести знаков для нежестких уравнений в выбранном окне.

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

• Линейное дифференциальное уравнение — Википедия
• Характеристическое уравнение — Википедия
• Метод неопределённых коэффициентов — Википедия
• Гармонический осциллятор — Википедия
• Фазовая плоскость — Википедия
• Метод Рунге — Кутты — Википедия