Приложение 3
Вычисление матричной экспоненты
Матричную экспоненту eAt, где А — квадратная матрица порядка n, можно вычислить из степенного разложения
,
откуда следует, что и eAt представляет собой квадратную матрицу порядка n.
Численное определение элементов eAt осуществляется с помощью следующей процедуры: 1) с помощью приведенного отрезка ряда Тейлора находят достаточно точное значение j0 = eAh для малого h; 2) последовательно находят величины jk+1 = jk2 (k = 0, 1, 2, …, n), определяющие при H = 2h, 4h, 16h, …
Аналитические выражения элементов матрицы eAt можно получить с помощью формулы Сильвестра
,
где λ —корни характеристического уравнения — собственные числа матрицы А.
Для n = 2, когда
,
элементы матрицы eAt выражаются следующим образом:
,
где λ1,2 — корни характеристического уравнения
.
Для комплексно сопряженной пары корней λ1,2 = – δ ± jω приведенное выражение преобразуется к виду
.