Plotting graph with Alpha and spectral radius

9 次查看（过去 30 天）
HAHA
HAHA 2021-1-10
编辑： HAHA 2021-1-21
hello

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

采纳的回答

Alan Stevens
Alan Stevens 2021-1-10
It works if you set i = 0.1:0.1:10. With i = 0; the inverse of S is full of Inf's.
  2 个评论
Alan Stevens
Alan Stevens 2021-1-10
You have to leave out i = 0. The following
A = [5 -1 0 0 0; -2 5 -1 0 0; 0 -2 5 -1 0; 0 0 -2 5 -1; 0 0 0 -2 5];
H = 0.5*(A+A');
S = 0.5*(A-A');
P_M=[];
alpha=[];
for i = [-0.5:0.1:-0.1 0.1:0.1:10] % leave out i = 0
alpha=[alpha i];
s_radius=norm(eig(inv(i.*eye(5)+S)*(i.*eye(5)-H)*inv(i.*eye(5)+H)*(i.*eye(5)-S)),inf);
P_M=[P_M s_radius];
end
plot(alpha,P_M), grid
xlabel('alpha'),ylabel('P M')
produces this for me:
Alan Stevens
Alan Stevens 2021-1-10
In a previous comment you wrote "Unfortunately there is the same error but anyway i think i should set it as i = -0.5:0.1:10 as alpha is in the open set of 0-10.", which is why I set alpha from -0.5. However, if you just start i from 0.1 you should be able to get an appropriate graph.

请先登录，再进行评论。

更多回答（0 个）

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by