EIGENVECTORS FROM EIGENVALUES论文matlab代码

本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

论文原文见预印本论坛：题目Eigenvectors from Eigenvalues
PDF: https://arxiv.org/pdf/1908.03795.pdf

原来例子的python版本见链接：https://blog.csdn.net/appleyuchi/article/details/103092344

本文给出一个简单的matlab实现，不具有通用性，慎用，学习讨论之用。

% ocnzhao
% 20191117
clc
clear
format rat
%A=[-1 4 1; 4 3 0; 1 0 2];                    % A=[a, M*;X,M\in R^{2\times 2}
randc=rand(1,2);
A=[rand, randc;randc',diag(rand(1,2))]
[m,n]=size(A);
[AV,AD,~] = eig(A); %full matrix W whose columns are the corresponding left eigenvectors,
% returns diagonal matrix D of eigenvalues and matrix V whose columns are the
% corresponding right eigenvectors, so that A*V = V*D.

%% verify whether the vector is a unit vector
for ii=1:m
  disp(sprintf("The norm of %d th vector is: %d",ii, norm(AV(:,ii),2) ))
end
% The formula
% \left|v_{i, 1}\right|^{2}=\frac{\prod_{k=1}^{n-1}\left(\lambda_{i}(A)-\lambda_{k}(M)\right)}
% {\prod_{k=1 ; k \neq i}^{n}\left(\lambda_{i}(A)-\lambda_{k}(A)\right)}

% \prod_{k=1}^{n-1}\left(\lambda_{i}(A)-\lambda_{k}(M)\right)
indexi=1;     % The ith eigenvector
indexj=1;     % The jth entry of the ith eigenvector v_{i,j}
indexdifrow=setdiff(1:m,indexi);
indexdifcol=setdiff(1:n,indexj);

M=A(indexdifrow,indexdifcol)
[eigmV,eigmD]=eig(M);
num=1;
for knum=1:n-1
   num=num*(AD(indexi,indexj)-eigmD(knum,knum));
end
% \prod_{k=1 ; k \neq i}^{n}\left(\lambda_{i}(A)-\lambda_{k}(A)\right)
den=1;
for knum=1:n
    if indexi~=knum
        den=den*(AD(indexi,indexj)-AD(knum,knum));
    end
end

% test the formula
if norm(num/den-AV(indexi,indexj)^2,'fro')<10^(-9)
    disp(sprintf("The %d th vector v_i%d is true", indexi, indexj))
end