Математика / Матрицы, определители

Степень диагонализируемой матрицы

Если A=PDP^{-1}, то степень A^k вычисляется как PD^kP^{-1}. Диагональную матрицу D возводят в степень по диагональным элементам.

Опубликовано: Обновлено:

Формула

$$A^k=PD^kP^{-1},\quad D^k=\operatorname{diag}(\lambda_1^k,\ldots,\lambda_n^k)$$
matrix-power-diagonal Степени действуют по собственным координатам

Схема показывает, как компоненты в собственном базисе умножаются на lambda_i^k.

Многократное применение оператора становится независимым ростом или затуханием собственных компонент.

Обозначения

$k$
целая неотрицательная степень, число шагов
$D^k$
степень диагональной матрицы, диагональная матрица
$\lambda_i^k$
k-я степень i-го собственного значения, число
$P$
матрица собственных векторов, матрица перехода

Условия применения

  • Матрица A должна быть диагонализируемой: A=PDP^{-1}.
  • Для отрицательных степеней дополнительно нужно, чтобы все собственные значения были ненулевыми и A была обратима.
  • Порядок собственных значений в D должен соответствовать порядку столбцов P.

Ограничения

  • Если A не диагонализируема, формула с диагональной D неприменима.
  • При больших k численные ошибки могут сильно расти, если P плохо обусловлена.
  • Нулевые собственные значения запрещают отрицательные степени.

Подробное объяснение

Формула следует из сокращения множителей. Если A=PDP^{-1}, то A^2=(PDP^{-1})(PDP^{-1})=PD(P^{-1}P)DP^{-1}=PD^2P^{-1}. По индукции получается A^k=PD^kP^{-1}. Главная причина простоты - P^{-1}P=I между соседними множителями.

Диагональная матрица возводится в степень покомпонентно. Если D=diag(lambda1,...,lambdan), то D^k=diag(lambda1^k,...,lambdan^k). Поэтому сложная операция с матрицей превращается в степени чисел. Это одно из главных практических преимуществ диагонализации.

Геометрически A^k означает k-кратное применение оператора. В собственном базисе это k-кратное масштабирование каждой координаты своим lambda_i. Если |lambda_i|<1, соответствующая компонента затухает. Если |lambda_i|>1, растет. Если lambda_i отрицательно, знак компоненты чередуется. Так спектр управляет долгосрочным поведением системы.

Формула не только ускоряет вычисления, но и дает качественный анализ. В рекурсиях, графах, марковских цепях и линейных моделях часто важнее понять доминирующее собственное значение, чем просто получить одну степень матрицы.

Как пользоваться формулой

  1. Проверьте, что A диагонализируема.
  2. Найдите P и D такие, что A=PDP^{-1}.
  3. Возведите каждое диагональное значение D в степень k.
  4. Соберите D^k.
  5. Вычислите A^k=PD^kP^{-1} или примените эту формулу к вектору.

Историческая справка

Использование диагонализации для степеней матриц связано с развитием спектральных методов в линейных системах и динамике. Как только матрица стала рассматриваться как оператор, возник естественный вопрос о повторном применении: что происходит после k шагов. Собственные значения дали ответ через независимые режимы роста и затухания. Эта идея стала важной в рекуррентных соотношениях, механике, вероятности, численных методах и анализе графов. Исторически она опирается на ту же линию, что и диагонализация: характеристические корни, матричная алгебра и приведение оператора к простому виду.

В XIX веке диагонализация выросла из задач о приведении линейных преобразований и квадратичных форм к более простому виду. Смысл был не в красивой записи ради записи, а в том, чтобы заменить связанную систему координат независимыми направлениями. Позже тот же язык стал базовым для дифференциальных уравнений, механики, статистики и численных методов. Поэтому страница "Степень диагонализируемой матрицы" находится на границе между абстрактной алгеброй и вычислительной практикой: она объясняет, когда сложный оператор можно читать как набор независимых масштабирований.

Историческая линия формулы

Формула A^k=PD^kP^{-1} является прямым следствием диагонализации и композиции линейных отображений. Ее исторический контекст связан с матричной алгеброй Кэли, операторной линией Гамильтона и Фробениуса и развитием спектральных методов.

Пример

Пусть A=PDP^{-1}, где D=diag(3,1). Тогда A^5=PD^5P^{-1}=P diag(3^5,1^5) P^{-1}=P diag(243,1) P^{-1}. Не нужно пять раз перемножать A на себя. Если вектор x разложен в собственном базисе как x=c1v1+c2v2, то A^5x=3^5 c1v1+1^5 c2v2. То есть каждая собственная компонента развивается независимо. Это особенно удобно в рекурсиях x_{k+1}=Ax_k: после диагонализации x_k=A^k x_0=PD^kP^{-1}x_0. Дополнительная проверка для этой страницы: после преобразования нужно перемножить матрицы обратно и убедиться, что получается исходный оператор. Для темы "Степень диагонализируемой матрицы" это особенно полезно, потому что диагональная форма часто выглядит убедительно сама по себе, но ошибка в порядке базисных векторов или в одном собственном значении сразу ломает равенство. В учебной задаче удобно отдельно выписать матрицу перехода, диагональную матрицу и обратную матрицу перехода, а затем проверить хотя бы один столбец произведения. Такой контроль показывает не только численный ответ, но и то, какие направления пространства стали главными.

Частая ошибка

Частая ошибка - возводить P и P^{-1} в степень отдельно и писать P^kD^k(P^{-1})^k. Это неверно: промежуточные P^{-1}P сокращаются при перемножении, поэтому остается PD^kP^{-1}. Вторая ошибка - применять формулу без проверки диагонализируемости. Третья ошибка - забыть, что D^k имеет степени собственных значений на диагонали, а не kD. Еще одна ошибка - использовать отрицательные степени при нулевом собственном значении.

Практика

Задачи с решением

Степень через D

Условие. Если D=diag(2,-1,3), чему равно D^4?

Решение. Возводим диагональные элементы в четвертую степень: 2^4=16, (-1)^4=1, 3^4=81.

Ответ. D^4=diag(16,1,81).

Использовать диагонализацию

Условие. A=P diag(5,2) P^{-1}. Запишите A^3.

Решение. A^3=P diag(5^3,2^3) P^{-1}=P diag(125,8) P^{-1}.

Ответ. A^3=P diag(125,8) P^{-1}.

Дополнительные источники

  • MIT OpenCourseWare 18.06SC, Matrix powers and diagonalization
  • Jim Hefferon, Linear Algebra, powers of diagonalizable matrices
  • TUDelft Interactive Linear Algebra, diagonalization applications

Связанные формулы

Математика

Диагонализация матрицы

$A=PDP^{-1}$

Диагонализация представляет квадратную матрицу A в виде A=PDP^{-1}, где D диагональна, а столбцы P являются собственными векторами A. Это переводит действие оператора в собственный базис.

Математика

Функция от диагонализируемой матрицы

$f(A)=P f(D) P^{-1},\quad f(D)=\operatorname{diag}(f(\lambda_1),\ldots,f(\lambda_n))$

Если A диагонализируема, функцию от матрицы можно вычислить через функцию от ее собственных значений. Для диагональной D функция применяется к каждому диагональному элементу.

Математика

Спектр матрицы

$\sigma(A)=\{\lambda:\det(A-\lambda I)=0\}$

Спектр матрицы - это множество ее собственных значений. Для конечной квадратной матрицы он состоит из корней характеристического уравнения.