Plotting graph with Alpha and spectral radius
显示 更早的评论
采纳的回答
Alan Stevens
2021-1-10
0 个投票
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
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
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 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Agriculture 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
