Системы линейных ОДУ: матричный подход

От уравнения к системе

Одиночное ОДУ высокого порядка y^(n) = f(x, y, y', ..., y^(n-1)) эквивалентно системе первого порядка. Вводим: x₁ = y, x₂ = y', ..., xₙ = y^(n-1). Тогда x₁' = x₂, x₂' = x₃, ..., xₙ' = f(x, x₁, ..., xₙ). Это стандартная форма системы. Переход к матричной записи — не просто удобная нотация, но мощный инструмент: алгебра матриц позволяет применить весь арсенал линейной алгебры.

В механике системы ОДУ описывают движение нескольких взаимодействующих тел. В электронике — ток и напряжение в сложных схемах. В экологии — совместную динамику нескольких популяций (модель хищник-жертва).

Матричная запись системы

Система: x'₁ = a₁₁x₁ + ... + a₁ₙxₙ + f₁, ..., x'ₙ = aₙ₁x₁ + ... + aₙₙxₙ + fₙ.

В векторной форме: x' = A(t) x + f(t), где x = (x₁, ..., xₙ)ᵀ — вектор-столбец состояния, A — матрица n×n коэффициентов, f — вектор неоднородности.

Задача Коши: x(t₀) = x₀.

Фундаментальная матрица и принцип суперпозиции

Фундаментальная матрица Φ(t) — матрица n×n, столбцы которой образуют ФСС (n линейно независимых решений однородной системы x' = Ax). Определитель Вронского: det Φ(t) ≠ 0.

Общее решение однородной системы: x(t) = Φ(t) c, где c — произвольный вектор констант.

Метод вариации постоянных (Дюамель): Для неоднородной системы:

x(t) = Φ(t)c + Φ(t) ∫_{t₀}^t Φ⁻¹(s) f(s) ds.

Второй член — принцип Дюамеля: отклик системы на непрерывное воздействие f складывается из откликов на мгновенные импульсы.

Системы с постоянной матрицей: матричная экспонента

Для x' = Ax (A — постоянная матрица) решение: x(t) = e^{At} x₀, где матричная экспонента:

e^{At} = I + At + (At)²/2! + (At)³/3! + ... = Σₖ₌₀^∞ (At)ᵏ/k!.

Вычисление e^{At} через спектральное разложение. Если A = PDP⁻¹, где D = diag(λ₁, ..., λₙ) — диагональная матрица из собственных значений, P — матрица из собственных векторов, то:

e^{At} = P e^{Dt} P⁻¹ = P · diag(e^{λ₁t}, ..., e^{λₙt}) · P⁻¹.

Смысл: в базисе собственных векторов компоненты развиваются независимо, каждая как e^{λᵢt}.

Полный пример: система хищник-жертва (линеаризованная)

Модель Лотки–Вольтерры вблизи точки равновесия линеаризуется. Пусть A = [[0, 1], [−1, 0]] (матрица малых колебаний).

Собственные значения: det(A − λI) = λ² + 1 = 0 → λ = ±i.

Собственные векторы: для λ = i: (A − iI)v = 0. [[−i, 1], [−1, −i]]v = 0 → v₁ = (1, i)ᵀ.

Матричная экспонента: e^{At} = [[cos t, sin t], [−sin t, cos t]] (матрица поворота!).

Решение: (x(t), y(t))ᵀ = [[cos t, sin t], [−sin t, cos t]] · (x₀, y₀)ᵀ. Это эллиптические орбиты вокруг равновесия — популяции совершают замкнутые колебания вокруг точки равновесия, как предсказывает модель Лотки–Вольтерры для случая без потерь.

Жорданов случай: кратные собственные значения

Если A неприводима к диагональной форме (кратные собственные значения), используем жорданову форму A = PJP⁻¹.

Для жордановой клетки J(λ, k) размера k:

e^{J(λ,k)t} = e^{λt} · [[1, t, t²/2!, ..., t^{k-1}/(k-1)!], [0, 1, t, ..., t^{k-2}/(k-2)!], ..., [0, 0, ..., 1]].

Физически: кратные корни порождают множители вида tᵏe^{λt} — именно это и есть условие резонанса в матричном контексте.

Численный пример: A = [[1, 1], [0, 1]], λ = 1 (кратность 2).

e^{At} = e^t · [[1, t], [0, 1]]. Если x(0) = (1, 0)ᵀ: x(t) = (eᵗ, 0)ᵀ. Если x(0) = (0, 1)ᵀ: x(t) = (teᵗ, eᵗ)ᵀ. Первая компонента нарастает линейно — признак резонанса.

Вопрос для размышления: Почему собственные значения с отрицательной вещественной частью соответствуют затуханию, а с нулевой — нейтральной устойчивости? Что происходит при λ = 0?

Практические выводы для инженера

Матричный подход позволяет сразу по знаку вещественных частей собственных значений оценить поведение всей системы, не решая систему явно. Это главный инструмент инженера-управленца: перед настройкой регулятора нужно вычислить собственные значения матрицы замкнутой системы и убедиться, что все они лежат в левой полуплоскости (Re λ < 0). Именно поэтому знание линейной алгебры — обязательная составляющая теории управления.

Метод вариации параметров для систем

Для неоднородной системы x' = Ax + g(t) метод вариации параметров даёт: x(t) = Φ(t)c + Φ(t)∫Φ⁻¹(s)g(s) ds, где Φ(t) — фундаментальная матрица. Это матричная формула Дюамеля. Передаточная матрица: G(t,s) = Φ(t)Φ⁻¹(s) описывает реакцию системы на импульсный вход в момент s. В теории управления это матрица импульсных переходных функций. Устойчивость и G: система асимптотически устойчива тогда и только тогда, когда ∫₀^∞ ‖G(t,0)‖ dt < ∞ — интегральная характеристика затухания переходных процессов, важная для оценки качества регулятора.