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

3 visualizaciones (últimos 30 días)
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

Respuestas (1)

John D'Errico
John D'Errico el 1 de Mayo de 2023
Editada: John D'Errico el 1 de Mayo de 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 comentario
Kamilu Sanusi
Kamilu Sanusi el 1 de Mayo de 2023
Editada: Kamilu Sanusi el 1 de Mayo de 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

Iniciar sesión para comentar.

Categorías

Más información sobre Linear Algebra en Help Center y File Exchange.

Productos


Versión

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by