Линейные ОДУ высших порядков: структура общего решения

Физическая мотивация: колебания и волны

Закон Ньютона для пружинного маятника с демпфированием: mẍ + cẋ + kx = F(t). Это линейное ОДУ второго порядка с постоянными коэффициентами. Параметры m (масса), c (коэффициент демпфирования) и k (жёсткость пружины) полностью определяют поведение системы. Три качественно разных режима — колебательный, переходный (критический) и апериодический — напрямую связаны с типом корней характеристического уравнения. Понимание структуры решений этого уравнения — ключ ко всей теории колебаний в физике и инженерии.

Общая теория линейных ОДУ

Линейное ОДУ n-го порядка: L[y] = y^(n) + p_{n-1}(x) y^{(n-1)} + ... + p₁(x) y' + p₀(x) y = f(x).

Оператор L линеен: L[c₁y₁ + c₂y₂] = c₁L[y₁] + c₂L[y₂].

Ключевая теорема об общем решении: Общее решение неоднородного уравнения имеет вид y = y_ч + y_о, где y_ч — любое частное решение неоднородного, а y_о — общее решение соответствующего однородного уравнения L[y] = 0.

Из этого следует, что нужно решить две задачи: (1) найти пространство решений однородного уравнения; (2) найти хотя бы одно частное решение неоднородного.

Определитель Вронского

Если y₁, ..., yₙ — n решений однородного уравнения, то определитель Вронского:

W(y₁, ..., yₙ)(x) = det[yₖ^(j-1)(x)] для k = 1, ..., n и j = 1, ..., n.

Теорема Абеля: W(x) = W(x₀) · exp(−∫{x₀}^x p{n-1}(t) dt). Это означает: W либо нигде не равен нулю, либо тождественно равен нулю. Первый случай — решения линейно независимы (образуют фундаментальную систему решений, ФСС).

Практическое значение: Если вронскиан не равен нулю, все n решений «вносят независимый вклад» в общее решение — ни одно из них нельзя выразить через остальные. Именно ФСС формирует «базис» пространства решений.

Уравнения с постоянными коэффициентами

L[y] = y^(n) + a_{n-1} y^{(n-1)} + ... + a₁ y' + a₀ y = 0.

Пробное решение y = e^{λx}: каждая производная y^(k) = λᵏ e^{λx}. Подставляем:

(λⁿ + a_{n-1} λ^{n-1} + ... + a₁ λ + a₀) e^{λx} = 0.

Поскольку e^{λx} ≠ 0, получаем характеристическое уравнение: P(λ) = λⁿ + a_{n-1} λ^{n-1} + ... + a₀ = 0.

Случай 1. Различные вещественные корни λ₁, ..., λₙ. ФСС = {e^{λ₁x}, ..., e^{λₙx}}. Общее решение: y = c₁ e^{λ₁x} + ... + cₙ e^{λₙx}.

Случай 2. Кратный вещественный корень λ кратности k. Отвечающие ему k линейно независимых решений: e^{λx}, xe^{λx}, x²e^{λx}, ..., x^{k-1}e^{λx}.

Смысл: кратный корень «порождает» полиномиальный множитель. Это математически объясняет резонансные явления в физике.

Случай 3. Комплексные корни α ± βi. Вещественные решения: e^{αx} cos(βx) и e^{αx} sin(βx).

Пружинный маятник: полный разбор

Уравнение: ẍ + 2bẋ + ω₀²x = 0, где b = c/(2m) — коэффициент затухания, ω₀² = k/m — квадрат собственной частоты.

Характеристическое уравнение: λ² + 2bλ + ω₀² = 0. Корни: λ = −b ± √(b² − ω₀²).

Случай b < ω₀ (слабое демпфирование): Корни комплексные: λ = −b ± iω, где ω = √(ω₀² − b²). Решение: x(t) = e^{−bt}(C₁ cos ωt + C₂ sin ωt) = Ae^{−bt} cos(ωt + φ). Затухающие колебания с амплитудой Ae^{−bt} и частотой ω < ω₀.

Случай b = ω₀ (критическое демпфирование): Кратный корень λ = −b. Решение: x(t) = (C₁ + C₂t)e^{−bt}. Апериодический возврат в равновесие — самый быстрый возврат без перехода нуля.

Случай b > ω₀ (сильное демпфирование): Два различных вещественных корня λ₁, λ₂ < 0. Решение: x(t) = C₁ e^{λ₁t} + C₂ e^{λ₂t}. Апериодическое затухание.

Пример с числами: Пружина с k = 4 Н/м, масса m = 1 кг, демпфирование c = 2 Нс/м. Уравнение: ẍ + 2ẋ + 4x = 0. Характеристическое: λ² + 2λ + 4 = 0. Дискриминант: 4 − 16 = −12 < 0. Корни: λ = −1 ± i√3. Решение: x(t) = e^{−t}(C₁ cos(√3 t) + C₂ sin(√3 t)). При начальных условиях x(0) = 0.1 м, ẋ(0) = 0: C₁ = 0.1, C₂ = 0.1/√3 ≈ 0.058. Амплитуда через 1 с: 0.1e^{−1} ≈ 0.037 м — уменьшилась в 2.7 раза.

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

Метод вариации параметров для неоднородных уравнений

Для уравнения y'' + py' + qy = f(x) с известным фундаментальным набором {y₁, y₂} метод вариации параметров ищет решение y = C₁(x)y₁ + C₂(x)y₂. Условие C₁'y₁ + C₂'y₂ = 0 и C₁'y₁' + C₂'y₂' = f дают систему для C₁' и C₂' через определитель Вронского W = y₁y₂' − y₂y₁'. Итог: C₁'= −y₂f/W, C₂' = y₁f/W. Преимущество метода: работает для любой правой части f(x), даже если метод неопределённых коэффициентов неприменим (например, для f = tan x или f = e^x/x, где антипроизводная не выражается элементарными функциями).