2

I'm trying to convert a generalized eigenvalue problem into a normal eigenvalue calculation.

I have this code:

[V,D,flag] = eigs(A, T);

Now I convert it into:

A1 = inv(T)*A;
[V1,D1,flag1] = eigs(A1);

Shouldn't I get the same result? From what I understand in the Matlab documentation, the first equation solves:

A*V = B*V*D

and the second one solves:

A*V = V*D

am I missing something?

Thanks!!

4

2 に答える 2

5

簡単な例：

A = rand(4); B = randn(4); B = B'*B;         %'# some matrices
[VV,DD] = eig(B\A);
[V,D] = eigs(A,B);
V = bsxfun(@rdivide, V, sqrt(sum(V.*V)));    %# make: norm(V(:,i))==1

結果：

V =
     -0.64581       0.8378      0.77771      0.50851
      0.70571     -0.51601     -0.32503     -0.70623
      0.27278     0.076874     -0.51777      0.25359
      0.10245      0.16095     -0.14641     -0.42232
VV =
     -0.64581       0.8378     -0.77771     -0.50851
      0.70571     -0.51601      0.32503      0.70623
      0.27278     0.076874      0.51777     -0.25359
      0.10245      0.16095      0.14641      0.42232
D =
       17.088            0            0            0
            0      0.27955            0            0
            0            0     -0.16734            0
            0            0            0     0.027889
DD =
       17.088            0            0            0
            0      0.27955            0            0
            0            0     -0.16734            0
            0            0            0     0.027889

注：固有値は常に同じようにソートされるとは限りません。また、符号の規則が異なる場合があります...

于 2011-07-19T14:00:40.023 に答える
0

まずTが可逆かどうかを確認します。第二に、私は確信D = D1しており、それV = V1はスケールファクターまでです。の各列がの対応する列と倍率までV1同じかどうかを確認します (つまり、 を見てください)。VV./V1

于 2011-07-19T13:48:37.617 に答える