pz_plxy.m

the code of the book come form the book of marlix laboratory with the book of singal

function pz_plxy(f1,f2,k,p,z)
%根据系统零极点分布绘制系统频率响应曲线程序
%f1、f2：绘制频率响应曲线的频率范围（即频率起始和终止点，单位为赫兹）
%p、z：系统函数极点和零点位置行向量
%k：绘制频率响应曲线的频率取样间隔
p=p';
z=z';
f=f1:k:f2;									%定义绘制系统频率响应曲线的频率范围
w=f*(2*pi);
y=i*w;
n=length(p);
m=length(z);
if n==0											%如果系统无极点
   yq=ones(m,1)*y; 
   vq=yq-z*ones(1,length(w));
   bj=abs(vq);
   ai=1;
elseif m==0									%如果系统无零点
   yp=ones(n,1)*y;
   vp=yp-p*ones(1,length(w));
   ai=abs(vp);
   bj=1;
else 
   yp=ones(n,1)*y;
   yq=ones(m,1)*y;
   vp=yp-p*ones(1,length(w));
   vq=yq-z*ones(1,length(w));
   ai=abs(vp);
   bj=abs(vq);
end
Hw=prod(bj,1)./prod(ai,1);
plot(f,Hw);
title('   系统幅频响应曲线')
xlabel(' 频率w（单位：赫兹）')
ylabel('H(jw)')