Я проверяю теорему, что A = Q * Лямбда * Q_inverse, где Q Матрица с собственными векторами, а лямбда - диагональная матрица, имеющая собственные значения в диагонали.
Мой код следующий:
import numpy as np
from numpy import linalg as lg
Eigenvalues, Eigenvectors = lg.eigh(np.array([
[1, 3],
[2, 5]
]))
Lambda = np.diag(Eigenvalues)
Eigenvectors @ Lambda @ lg.inv(Eigenvectors)
Что возвращает:
array([[ 1., 2.],
[ 2., 5.]])
не должна ли возвращенная матрица совпадать с исходной матрицей, которая была разложена?
Ваш совет будет оценен.