Я предполагаю, что вы имели в виду x(A-I)=0
, поскольку то, что вы написали, не имеет для меня никакого смысла.Уравнение, которое я написал, подразумевает, что x*A^n=x
, что обычно подразумевается под устойчивым состоянием.Решением уравнения является левый собственный вектор A
с собственным значением 1
.
Вы можете получить собственные векторы и собственные значения A
, используя функцию eig
.
A = [0.4, 0.1, 0.2; 0.3, 0.7, 0.7; 0.3, 0.2, 0.1];
% Get the eigenvalues (D) and left eigenvectors (W)
[~,D,W] = eig(A);
% Get the index of the eigenvalue closest to 1
[~,idx] = min(abs(diag(D)-1));
% Get associated eigenvector
x = W(:,idx).';
Проверка решения
>> all(abs(x*(A-eye(size(A)))) < 1e-10)
ans =
logical
1