k=1
num=[k];
den=[1 5 8 6 0];
sys=tf(num,den);
%求零极点
[p z]=pzmap(num,den)
%根轨迹图
rlocus(sys);
p =
0.0000 + 0.0000i
-3.0000 + 0.0000i
-1.0000 + 1.0000i
-1.0000 - 1.0000i