%===========================================================
%  程序名：zero_pole 
%  本程序繪制系統函數的零點、極點在s平面上的分布圖。
%  A:系統函數分母多項式系數向量，實驗者從鍵盤上輸入
%  B:系統函數分子多項式系數向量，實驗者從鍵盤上輸入
%  p:函數返回的系統函數極點位置行向量
%  q:函數返回的系統函數零點位置行向量
%===========================================================
b=input('請輸入系統函數分子多項式系數：');
a=input('請輸入系統函數分母多項式系數：');
p=roots(a);                       %求系統極點
q=roots(b);                       %求系統零點
p=p';					          %將極點列向量轉置為行向量
q=q';					          %將零點列向量轉置為行向量
x=max(abs([p q]));				  %確定縱坐標范圍
x=x+1;
y=x;						      %確定橫坐標范圍
clf
%subplot(2,1,1);
axis([-x x -y y]);		          %確定坐標軸顯示范圍
axis('square')

%===========================================================
plot([-x x],[0 0]);hold on; 	       %畫橫坐標軸
plot([0 0],[-y y]);hold on;		       %畫縱坐標軸
plot(real(p),imag(p),'x');hold on;     %畫極點
plot(real(q),imag(q),'o');hold on;     %畫零點
grid on;
title('系統函數零極點分布');
ylabel('虛軸');
xlabel('實軸');