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

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

Линейное ОДУ второго порядка a·y″ + b·y′ + c·y = g(x)

Быстрые примеры:

Без демпфирования С недодемпфированием Критическое демпфирование С передемпфированием Вынужденные (const) Затухающая сила Неустойчивое

Коэффициенты

a·y″

b·y′

c·y

Уравнение имеет вид a·y″ + b·y′ + c·y = g(x), где a, b, c — вещественные константы. Коэффициент a не должен быть равен нулю.

Вынуждающий член (правая часть)

g(x) =

Аналитическое решение y_p выводится для констант, многочленов (до 2-й степени) и A·exp(k·x). Для любых других выражений строится высокоточное численное решение.

Начальные условия и диапазон графика

x₀

y₀

y′₀

x min

x max

Требуется два начальных условия: y(x₀) = y₀ и y′(x₀) = y′₀. Решатель выполняет интегрирование от x₀ в обоих направлениях x.

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

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

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

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

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

a · y″(x) + b · y′(x) + c · y(x) = g(x)

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

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

y(x) = y_h(x) + y_p(x)

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

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

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

a · r² + b · r + c = 0

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

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

Дискриминант Δ = b² − 4ac	Корни a·r² + b·r + c = 0	Однородное решение y_h(x)	Физический режим
Δ > 0	Два различных вещественных корня r₁, r₂	C₁·e^(r₁·x) + C₂·e^(r₂·x)	Передемпфирование — медленный монотонный возврат к равновесию (если оба корня отрицательны).
Δ = 0	Один кратный вещественный корень r = −b/(2a)	(C₁ + C₂·x)·e^(r·x)	Критическое демпфирование — максимально быстрый возврат без колебаний.
Δ < 0	Комплексно-сопряженные корни α ± β·i, где α = −b/(2a), β = √(−Δ)/(2a)	e^(α·x)·(C₁·cos(β·x) + C₂·sin(β·x))	Недодемпфирование (α < 0) — затухающие колебания; без демпфирования (α = 0) — чистые гармонические колебания; неустойчивое (α > 0) — нарастающие колебания.

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

Когда 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 — вынужденная система с недодемпфированием.

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

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

Введите коэффициенты a, b, c в верхней строке. a не должно быть равно нулю (иначе уравнение станет первым порядком).
Введите вынуждающий член g(x) или оставьте 0 для однородной задачи. Аналитические частные решения выводятся для констант, многочленов до 2-й степени и одиночных экспонент A·e^(k·x), включая случай резонанса.
Укажите начальные условия (x₀, y₀, y′₀). Необходимо задать и y, и y′ в точке x₀, так как это уравнение второго порядка.
Выберите диапазон x для графиков. Решатель выполняет интегрирование от x₀ в обоих направлениях x с помощью RK4.
Нажмите "Решить и визуализировать". Вы получите характеристическое уравнение с его корнями на комплексной плоскости, классификацию режима демпфирования, аналитические однородное и частное решения, временной график 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′ = v v′ = ( g(x) − b·v − c·u ) / a

с начальными условиями 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 подшагов на каждом интервале обеспечивают точность около шести знаков для нежестких уравнений в выбранном окне.

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

Ссылайтесь на этот контент, страницу или инструмент так:

"Решатель ОДУ второго порядка" на сайте https://ru.miniWebtool.com/решатель-оду-второго-порядка/ от MiniWebtool, https://MiniWebtool.com/

от команды miniwebtool. Обновлено: 22 апреля 2026 г.

Вы также можете попробовать наш AI Решатель Математических Задач GPT, чтобы решить ваши математические проблемы с помощью вопросов и ответов на естественном языке.