Eigs passes the wrong StartVector

Question

Sukhi Singh 2021-12-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1602075-eigs-passes-the-wrong-startvector

回答： Christine Tobler 2021-12-3

Hi all,

I'm using MATLAB 2019b 64-bits. I'm calling eigs() to find the largest eigenvector of a matrix by passing it a function handle and specifying a 'StartVector' as follows:

A = reshape(1:9,3,3); % create some matrix for testing
StartVector = [0.4 0.9 1.7]';  % specify a start vector for eigs()
[vec,ev] = eigs(@(vec)test_eigs(vec,A),3,1,'largestabs','Display',0,'IsFunctionSymmetric',0,'MaxIterations',300,'StartVector',StartVector);

The definition of test_eigs() is as follows:

function vec = test_eigs(vec,A)
vec = A*vec;
end

However, eigs does not pass the specified start vector to the first call of test_eigs(). Instead, it passes the vector [1 0 0]' on all runs of the program. What am I missing? I'd appreciate any pointers!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Christine Tobler 2021-12-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1602075-eigs-passes-the-wrong-startvector#answer_846810

在 MATLAB Online 中打开

The standard method used in EIGS is only efficient for matrices that are quite large. When eigs detects that the input matrix is so small that it will be more accurate and faster to just call EIG instead and extract the requested values from that, it does this instead and most options that affect the standard method's behavior are ignored as that method isn't used.

When a function handle is passed, this means that the function handle is just called on each column of the identity matrix one at a time, to get a matrix representation that can be passed to EIG.

If you increase the size of the input matrix, you'll see the 'StartVector' be passed in on the first call to the function handle:

n = 21;
A = ones(n);
StartVector = flip(1:n)';
[vec,ev] = eigs(@(vec)test_eigs(vec,A),n,1,'largestabs','MaxIterations',300,'StartVector',StartVector);