function dplxy(k,r,A,B)
%The function to draw the frequency response of discrete system
p=roots(A);                     %求極點
q=roots(B);						%求零點
figure(1)
ljdt(A,B)						%畫零極點圖
w=0:l*pi/k:r*pi;				
y=exp(i*w);						%定義單位圓上的k個頻率等分點
N=length(p);					%求極點個數
M=length(q);					%求零點個數
yp=ones(N,1)*y;					%定義行數為極點個數的單位圓向量
yq=ones(M,1)*y;					%定義行數為零點個數的單位圓向量
vp=yp-p*ones(1,k+1);			%定義極點到單位圓上各點的向量
vq=yq-q*ones(1,k+1);			%定義零點到單位圓上各點的向量
Ai=abs(vp);						%求出極點到單位圓上各點的向量的模
Bj=abs(vq);						%求出零點到單位圓上各點的向量的模
Ci=angle(vp);					%求出極點到單位圓上各點的向量的相角
Dj=angle(vq);					%求出零點到單位圓上各點的向量的相角
fai=sum(Dj,1)-sum(Ci,1);		%求系統相頻響應
H=prod(Bj,1)./prod(Ai,1);		%求系統幅頻響應
figure(2)
plot(w,H);						%繪制幅頻特性曲線
title('離散系統幅頻特性曲線')
xlabel('角頻率')
ylabel('幅度')
figure(3)
plot(w,fai)						
title('離散系統的相頻特性曲線')
xlabel('角頻率')
ylabel('相位')