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

5 Ansichten (letzte 30 Tage)
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

Antworten (1)

John D'Errico
John D'Errico am 1 Mai 2023
Bearbeitet: John D'Errico am 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 Kommentar
Kamilu Sanusi
Kamilu Sanusi am 1 Mai 2023
Bearbeitet: Kamilu Sanusi am 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

Melden Sie sich an, um zu kommentieren.

Kategorien

Mehr zu Mathematics finden Sie in Help Center und File Exchange

Produkte


Version

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by