How to use Eigenvector and Eigenvalues of a matrix to formulate Entropy equation?

3 views (last 30 days)
Dear Matlab experts,
I have a matrix T = [T11, T12 ; T21, T22] of size , where all elements in T are 126*126.
After using this function [Val, Vect] = eig(T); I obtained matrices of Val() , and Vect (digonal).
Now I have eigenvactors and eigenvalues. I need to implement following experssions.
I have attached T matrix and crossponding eigenvalues and eigenvectors, I need to estimates both (1) and (2)
Thank you so much.

Accepted Answer

Bjorn Gustavsson
Bjorn Gustavsson on 6 Sep 2022
First you should extract the eigenvalues from the diagonal matrix (mainly for convenience):
vLambda = diag(Vect);
Then you want the sum of the "first two" for your P_i. Presumably "first two" means the two largest, though that's not made explicitly clear. Let's check where those are:
plot(vLambda,'.-')
So the last eigenvalues are the biggest. This gives us for P:
P = vLambda(end-1:end)/vLambda(end) + vLambda(end-1)
Then you're asked for the sum of P multiplied with acos(|u_i|). You should be able to figure that one out. Read the help and documentation of eig and think about what more you know about the eigenvectors (write these facts down in a list) and one fact of those can be used to some insight about acos.
Also since these questions make it seem as you are very new to matlab, I recommend you browse through the on-ramp material. It is designed to get you up and running as fast as possible.
HTH
  3 Comments
Bjorn Gustavsson
Bjorn Gustavsson on 6 Sep 2022
You can extract an eigenvector in your case like this:
u_1 = Val(:,end);
You can calculate it's norm like this:
n_u_1 = norm(u_1);
and I assume you already know about the acos-function.
Amjad Iqbal
Amjad Iqbal on 7 Sep 2022
Thank you so much for your detailed inputs. It helps a lot to understand to use eigen values and vecotors.

Sign in to comment.

More Answers (0)

Categories

Find more on Creating and Concatenating Matrices in Help Center and File Exchange

Products


Release

R2021a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!