gden=[1 16 51 36];                        %開環傳遞函數分母的表達式
a=2:0.1:10;n=length(a);                   %n為參數a向量的長度
for i=1:n                                    
    gnum=5*[1 a(i)];                      %gnum為a(i)時開環傳遞函數分子的表達式
    [bnum,bden]=feedback(gnum,gden,1,1);
    gen(i,1:3)=roots(bden)';              %求出a(i)時的根向量，共有三個根 
    clpole(i,1:4)=[gen(i,:),a(i)];     %clpole存放根向量與a參變量對應值組成矩陣
end
plot(gen,'k.');grid                    %繪制具有參變量a的根軌跡圖 
axis ([-4 0 -2 2])                      
%設置X軸，得出虛軸附近區域的根軌跡圖