CODEWITHSUNDEEP
matlabeigenvectoreigenvalue
user1234440
user1234440

Reputation: 23567

Eigen Values and Eigen Vectors Matlab

I have a matrix A

A = [ 124.6,95.3,42.7 ; 95.3,55.33,2.74 ; 42.7,2.74,33.33 ]

The eigenvalues and vectors:

[V,D] = eig(A)

How do I show the eigenvalues are mutually perpendicular? I've tried that if the dot product of the eigenvalues are zero, this demonstrates they are mutually perpendicular, but how would you compute this in MATLAB? I tried the following code

transpose(diag(D)) * diag(D)  %gives 4.1523e+04

Also, how can I verify the definition of eigenvalues and vector holds:

A e_i - L_i e_i = 0

The above equation: for i equal 1 to 3. For a real, symmetric matrix all eigenvales are positive numbers and eigenvectors for a basis for the matrix

I tried the following code, but it doesn't seem to give me 0. Any ideas?

A*V(1) - D(1)*V(1)

Upvotes: 0

Views: 637

Answers (2)

Luis Mendo
Luis Mendo

Reputation: 112659

To compute all dot products between columns of V:

M = squeeze(sum(bsxfun(@times, conj(V), permute(V, [1 3 2]))));

The columns of V (eigenvectors) will be orthogonal if the above M is diagonal.

Upvotes: 0

Matt J
Matt J

Reputation: 1137

To show orthogonality

>> V'*V-eye(size(V))

ans =

   1.0e-15 *

    0.2220    0.1110    0.2498
    0.1110   -0.4441    0.1388
    0.2498    0.1388    0.4441

To show that the definition of the eigendecomposition is satisfied,

>> A*V - V*D

ans =

   1.0e-13 *

    0.4086    0.0400    0.8527
    0.3908    0.0355    0.5684
    0.1954    0.0355         0

The results won't be exactly zero, because digital computers don't do exact math, but you can see that they're pretty close.

Upvotes: 2

Related Questions