% ex4_1
% 繪制奈魁斯特圖
figure('pos',[30 50 260 200],'color','w');
axes('pos',[0.15 0.2 0.7 0.7]);
G1=tf([1],[10 1]);
w=[0 logspace(-2,3,200)];
G1v=polyval(G1.num{1},w*j)./polyval(G1.den{1},w*j);
G2=zpk([],[-1+0.2j -1-0.2j],1);
G2t=tf(G2);
G2v=polyval(G2t.num{1},w*j)./polyval(G2t.den{1},w*j);
plot(real(G1v),imag(G1v),'k',real(G2v),imag(G2v),'k');