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

Question

Kamilu Sanusi 2023-5-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1955924-how-do-i-keep-track-of-eigenvalue-on-pzmap-for-a-pariticular-system-matrix-a-containing-a-variable-d

编辑： Kamilu Sanusi 2023-5-1

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.

Please how do I link the eigen values of the same value of D together by a line?
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
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

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

John D'Errico 2023-5-1

编辑：John D'Errico 2023-5-1

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.

https://www.mathworks.com/matlabcentral/fileexchange/22885-eigenshuffle?s_tid=srchtitle

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Kamilu Sanusi 2023-5-1

编辑：Kamilu Sanusi 2023-5-1

@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

请先登录，再进行评论。