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

5 次查看（过去 30 天）
Kamilu Sanusi
Kamilu Sanusi 2023-5-1
编辑： Kamilu Sanusi 2023-5-1
Ran in:
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

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

回答（1 个）

John D'Errico
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 个评论
Kamilu Sanusi
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

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

产品


版本

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by