Матрица Коши¶
Матрица Коши — это фундаментальная матрица системы линейных дифференциальных уравнений, решение которой в каждый момент времени равно единичной матрице. Она является решением матричного дифференциального уравнения и играет ключевую роль в теории обыкновенных дифференциальных уравнений, теории управления и математической физике. Матрица Коши позволяет выразить решение системы с начальными условиями в произвольный момент времени через начальное состояние, не требуя повторного решения задачи Коши.
¶Определение и основные свойства
Пусть дана система линейных дифференциальных уравнений первого порядка:
\[ \frac{d\mathbf{x}}{dt} = A(t) \mathbf{x}(t), \]
где \(\mathbf{x}(t) \in \mathbb{R}^n\) — вектор состояния, \(A(t)\) — квадратная матрица размером \(n \times n\), элементы которой являются непрерывными функциями времени на интервале \(I \subseteq \mathbb{R}\).
Матрицей Коши (или фундаментальной матрицей, нормированной в начальной точке) называется матричная функция \(C(t, s)\) размером \(n \times n\), удовлетворяющая следующим условиям:
- Для каждого фиксированного \(s \in I\) матрица \(C(t, s)\) как функция от \(t\) удовлетворяет матричному дифференциальному уравнению:
\[ \frac{\partial}{\partial t} C(t, s) = A(t) C(t, s). \]
- Начальное условие: \(C(s, s) = I\), где \(I\) — единичная матрица.
Таким образом, матрица Коши является решением задачи Коши для матричного уравнения. Она определена для всех \(t, s \in I\) и обладает следующими ключевыми свойствами:
- Групповое свойство: Для любых \(t, s, r \in I\) выполняется равенство:
\[ C(t, s) = C(t, r) C(r, s). \] В частности, \(C(t, s) = C(t, 0) C(0, s)\) и \(C(t, s) = [C(s, t)]^{-1}\).
- Обратимость: Матрица \(C(t, s)\) невырождена для всех \(t, s \in I\), то есть \(\det C(t, s) \neq 0\). Более того, \(\det C(t, s) = \exp\left( \int_s^t \operatorname{tr} A(\tau) d\tau \right)\), где \(\operatorname{tr} A\) — след матрицы \(A\).
- Связь с решением задачи Коши: Решение задачи Коши для неоднородной системы
\[ \frac{d\mathbf{x}}{dt} = A(t) \mathbf{x}(t) + \mathbf{f}(t), \quad \mathbf{x}(s) = \mathbf{x}_0 \] выражается через матрицу Коши по формуле Коши: \[ \mathbf{x}(t) = C(t, s) \mathbf{x}_0 + \int_s^t C(t, \tau) \mathbf{f}(\tau) d\tau. \]
¶Связь с фундаментальной матрицей
В теории дифференциальных уравнений часто используется понятие фундаментальной матрицы \(\Phi(t)\) — матрицы, столбцы которой образуют линейно независимые решения однородной системы. Фундаментальная матрица не нормирована в начальной точке: \(\Phi(s)\) — произвольная невырожденная матрица.
Матрица Коши \(C(t, s)\) является частным случаем фундаментальной матрицы, нормированной так, что в момент \(t = s\) она равна единичной матрице. Связь между ними: \[ C(t, s) = \Phi(t) \Phi^{-1}(s). \]
Таким образом, если известна любая фундаментальная матрица \(\Phi(t)\), матрица Коши строится по этой формуле. Это свойство широко используется для аналитического и численного построения матрицы Коши.
¶Вычисление матрицы Коши
¶Постоянная матрица \(A\)
Если матрица \(A\) не зависит от времени (\(A(t) \equiv A\)), то матрица Коши выражается через матричную экспоненту: \[ C(t, s) = e^{A(t-s)}. \] В этом случае она зависит только от разности аргументов, что отражает свойство инвариантности относительно сдвига времени.
¶Переменная матрица \(A(t)\)
В общем случае, когда \(A(t)\) зависит от времени, матрица Коши не выражается в виде простой экспоненты. Её можно найти следующими способами:
- Через фундаментальную матрицу: Если удаётся найти \(n\) линейно независимых решений однородной системы (например, методом вариации постоянных или методом подбора), то строится \(\Phi(t)\) и затем \(C(t, s) = \Phi(t) \Phi^{-1}(s)\).
- Через матричный ряд: Решение матричного уравнения можно представить в виде ряда, аналогичного экспоненциальному, но с учётом некоммутативности \(A(t)\) в разные моменты времени. Этот ряд называется магнусовским разложением или экспонентой, упорядоченной по времени:
\[ C(t, s) = I + \int_s^t A(\tau_1) d\tau_1 + \int_s^t A(\tau_1) \int_s^{\tau_1} A(\tau_2) d\tau_2 d\tau_1 + \dots \] Этот ряд сходится, если \(A(t)\) ограничена.
- Численное интегрирование: Матрицу Коши можно вычислить численно, решая матричное дифференциальное уравнение с начальным условием \(C(s, s) = I\) на интервале \([s, t]\). Для этого используются стандартные методы (Рунге-Кутты, Адамса и т.д.).
¶Применение
¶Теория управления
В теории линейных систем управления матрица Коши (часто называемая переходной матрицей состояния) является центральным объектом. Она позволяет:
- Выразить реакцию системы на произвольное входное воздействие (формула Коши).
- Анализировать устойчивость системы (по поведению \(C(t, s)\) при \(t \to \infty\)).
- Решать задачи наблюдаемости и управляемости.
- Строить дискретные модели непрерывных систем.
¶Математическая физика
В уравнениях в частных производных, сводящихся к системам ОДУ (например, метод разделения переменных, метод характеристик), матрица Коши используется для распространения начальных данных вдоль характеристик.
¶Численные методы
Матрица Коши применяется в методах стрельбы (shooting methods) для решения краевых задач. В этих методах краевая задача сводится к последовательности задач Коши, и матрица Коши позволяет линеаризовать отображение начальных условий в конечные.
¶Теория возмущений
При анализе систем с малыми параметрами матрица Коши разлагается в ряд по малому параметру, что даёт приближённые решения.
¶Пример
Рассмотрим систему второго порядка, сводящуюся к системе первого порядка: \[ \ddot{x} + \omega^2 x = 0, \quad \omega = \text{const}. \] Введём переменные: \(x_1 = x\), \(x_2 = \dot{x}\). Тогда система: \[ \frac{d}{dt} \begin{pmatrix} x_1 \\ x_2 \end{pmatrix} = \begin{pmatrix} 0 & 1 \\ -\omega^2 & 0 \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \end{pmatrix}. \] Матрица \(A\) постоянна. Матрица Коши: \[ C(t, s) = e^{A(t-s)} = \begin{pmatrix} \cos \omega (t-s) & \frac{1}{\omega} \sin \omega (t-s) \\ -\omega \sin \omega (t-s) & \cos \omega (t-s) \end{pmatrix}. \] Нетрудно проверить, что \(C(s, s) = I\) и что она удовлетворяет групповому свойству.
¶Интересные факты
- Матрица Коши названа в честь французского математика Огюстена Луи Коши, который внёс фундаментальный вклад в теорию дифференциальных уравнений.
- В русскоязычной литературе термин «матрица Коши» часто используется как синоним «фундаментальной матрицы, нормированной в точке», в то время как в англоязычной литературе чаще говорят о «state transition matrix» (матрица перехода состояния).
- Для систем с периодической матрицей \(A(t)\) (теория Флоке) матрица Коши имеет специальную структуру: \(C(t, s) = P(t) e^{R(t-s)} P^{-1}(s)\), где \(P(t)\) — периодическая матрица, а \(R\) — постоянная.
¶Критика и ограничения
Основное ограничение матрицы Коши связано с тем, что её аналитическое построение возможно лишь для узкого класса систем (постоянные, некоторые переменные). Для большинства реальных систем её приходится вычислять численно, что может быть ресурсоёмко для больших размерностей. Кроме того, при больших интервалах времени численное интегрирование может накапливать ошибки, особенно для жёстких систем.
¶Источники
- Понтрягин Л. С. «Обыкновенные дифференциальные уравнения». — М.: Наука, 1974.
- Арнольд В. И. «Обыкновенные дифференциальные уравнения». — М.: МЦНМО, 2012.
- Калман Р., Фалб П., Арбиб М. «Очерки по математической теории систем». — М.: Мир, 1971.
- Хартман Ф. «Обыкновенные дифференциальные уравнения». — М.: Мир, 1970.
- Гантмахер Ф. Р. «Теория матриц». — М.: Наука, 1967.