How do i keep track of eigenvalue on pzmap for a pariticular system matrix A containing a variable D?

5 vues (au cours des 30 derniers jours)
I would appreciate if you can help me with this problem. I have a system matrix A containing a damping variable D of which I want to see how it affects the position of eigenvalues on pzmap.
  1. Please how do I link the eigen values of the same value of D together by a line?
  2. How do i link only corresponding eigenvalue Lambda1 of different value oF D together, and the same thing applicable to other corresponding values of lampdas on the pz maz
  3. with arrow showing possible direction of motion of eigenvalue as a result of change in D
Ta1 = 24; Ta2 = 27; Ta3= 20;
H11 = -0.0641; H12 = 0.0359;
H21 = 0.1176; H22 = -0.2057;
H31 = 0.2077; H32 = 0.1961;
for D = [0 7 10 15 20 25 40 70]
A = [0 0 1 0 -1;0 0 0 1 -1;(-H11/Ta1) (-H12/Ta1) (-D/Ta1) 0 0;...
(-H21/Ta2) (-H22/Ta2) 0 (-D/Ta2) 0;(-H31/Ta3) (-H32/Ta3) 0 0 (-D/Ta3)];
Eig = eig(A);
a = Eig(1,1);
b = Eig(2,1);
c = Eig(3,1);
d = Eig(4,1);
e = Eig(5,1);
s = tf('s');
T = (1)/((s-a)*(s-b)*(s-c)*(s-d)*(s-e));
P = pole(T);
if D == 7
pzmap(T)
end
hold on
if D == 15
pzmap(T)
end
hold on
if D == 25
pzmap(T)
end
end

Réponses (1)

John D'Errico
John D'Errico le 1 Mai 2023
Modifié(e) : John D'Errico le 1 Mai 2023
You can use the tool I wrote, eigenshuffle. It sorts the eigenvalues in as consistent an order as possible, by assuming the corresponding eigenvectors are similar. You can find eigenshuffle on the file exchange.
  1 commentaire
Kamilu Sanusi
Kamilu Sanusi le 1 Mai 2023
Modifié(e) : Kamilu Sanusi le 1 Mai 2023
@John D'Errico, thank you for the link. I am just trying to learn MATLAB, i have checked your code over and over but i really do not understand how it be useful in my work. Please kindly guide in a way I can understand
Thank you and kind regards

Connectez-vous pour commenter.

Catégories

En savoir plus sur Linear Algebra dans Help Center et File Exchange

Produits


Version

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by