adoty = [Na 0 0 -1 ];
adotp = [0 Na 1 0 ];
qdot = [Mpa*wr Ma Mq -I*wr ];
rdot = [-Ma -Mpa*wr I*wr -Mq ];
A = [adoty; adotp; qdot; rdot];
EigA = eig(A)
So, for this matrix A, I can find the eigenvalues using this function. However, I would like to be able to get the eigenvalues as a function of wr. I can make a for loop to post a new eigenvalue per unit wr, but I don't have the ability to find the correlation between all those outputs since they're outputted as 4 different values. I tried to use eig(sym(A)), but the results end up being about 5 pages per answer, so that is useless to me. How do I find the correlation between my eigenvalues as a function of that wr value?
Attached is the program with all the numerical values if that is of any help.
